Abstract Using exact diagonalization, we investigate the many-body ground state for regular vortex patterns in a rotating Bose-condensed gas of N spinless particles, confined in a quasi-two-dimensional harmonic trap and interacting repulsively via finite-range Gaussian potential. The N -body Hamiltonian matrix is diagonalized in given subspaces of quantized total angular momentum L z , to obtain the lowest-energy eigenstate. Further, the internal structure of these eigenstates is analyzed by calculating the corresponding conditional probability distribution. Specifically, the quantum mechanically stable as well as unstable states in a co-rotating frame are examined in the moderately rotating regime corresponding to angular momenta 4 N ≤ L z < 5 N for N = 16 bosons. In response to externally impressed rotation, the patterns of singly quantized vortices are formed, shaping into canonical polygons with a central vortex at the trap center. The internal structure of unstable states reveals the mechanism of entry, nucleation and pattern formation of vortices with structural phase transition, as the condensate goes from one stable vortical state to the other. The stable polygonal vortex patterns having discrete p -fold rotational symmetry with p = 5 and p = 6 are observed. The hexagonal vortex pattern with p = 6 symmetry is a precursor to the triangular vortex lattice of singly quantized vortices in the thermodynamic limit. For unstable states, quantum melting of vortex patterns due to uncertainty in positions of individual vortices, is also briefly discussed.
Abstract Using exact diagonalization, we investigate the many-body ground state for regular vortex patterns in a rotating Bose-condensed gas of N spinless particles, confined in a quasi-two-dimensional harmonic trap and interacting repulsively via finite-range Gaussian potential. The N -body Hamiltonian matrix is diagonalized in given subspaces of quantized total angular momentum L z , to obtain the lowest-energy eigenstate. Further, the internal structure of these eigenstates is analyzed by calculating the corresponding conditional probability distribution. Specifically, the quantum mechanically stable as well as unstable states in a co-rotating frame are examined in the moderately rotating regime corresponding to angular momenta 4 N ≤ L z < 5 N for N = 16 bosons. In response to externally impressed rotation, the patterns of singly quantized vortices are formed, shaping into canonical polygons with a central vortex at the trap center. The internal structure of unstable states reveals the mechanism of entry, nucleation and pattern formation of vortices with structural phase transition, as the condensate goes from one stable vortical state to the other. The stable polygonal vortex patterns having discrete p -fold rotational symmetry with p = 5 and p = 6 are observed. The hexagonal vortex pattern with p = 6 symmetry is a precursor to the triangular vortex lattice of singly quantized vortices in the thermodynamic limit. For unstable states, quantum melting of vortex patterns due to uncertainty in positions of individual vortices, is also briefly discussed.
We present an exact diagonalization study of the breathing mode collective excitations for a rotating Bose-Einstein condensate of N=10 spinless bosons interacting via repulsive finiterange Gaussian potential and harmonically confined in quasi-two-dimension. The yrast state and the low-lying excited states are variationally obtained in given subspaces of the quantized total angular momentum L employing the beyond lowest Landau level approximation in slowly rotating regime with 0≤L<2N. For a given L, the low-energy eigenspectra (bands) are obtained in weakly to moderately interacting regime. Further, for a given interaction, the split in low-lying eigenenergies with increasing L is the precursor to spontaneous symmetry breaking of the axisymmetry associated with the entry of the first vortex. With increase in repulsive interaction, the value of the first breathing mode increases for stable total angular momentum states L=0 and N, but decreases for intermediate 0<L<N metastable states. The position of the observed first breathing modes in the eigenspectrum remains unchanged as the interaction is varied over several orders of magnitude.
We consider a gas of N(<15) Bose particles with hard-core repulsion, contained in a quasi-two-dimensional harmonic trap and subjected to an overall angular velocity $\ensuremath{\Omega}$ about the z axis. Exact diagonalization of the $n\ifmmode\times\else\texttimes\fi{}n$ many-body Hamiltonian matrix in given subspaces of the total (quantized) angular momentum Lz, with $n\ensuremath{\sim}{10}^{5}$ (e.g., for Lz=N=15,n=240782) was carried out using Davidson's algorithm. The many-body variational ground-state wave function, as also the corresponding energy and the reduced one-particle density-matrix $\ensuremath{\rho}(\mathbf{r},{\mathbf{r}}^{\ensuremath{'}})={\ensuremath{\sum}}_{\ensuremath{\mu}}{\ensuremath{\lambda}}_{\ensuremath{\mu}}{\ensuremath{\chi}}_{\ensuremath{\mu}}^{*}(\mathbf{r}){\ensuremath{\chi}}_{\ensuremath{\mu}}({\mathbf{r}}^{\ensuremath{'}})$ were determined. With the usual identification of $\ensuremath{\Omega}$ as the Lagrange multiplier associated with Lz for a rotating system, the ${L}_{z}\ensuremath{-}\ensuremath{\Omega}$ phase diagram (or the stability line) was determined that gave a number of critical angular velocities ${\ensuremath{\Omega}}_{\mathrm{ci}}, i=1,2,3,\dots{},$ at which the ground-state angular momentum and the associated condensate fraction, given by the largest eigenvalue of the reduced one-particle density matrix, undergo abrupt jumps. For a given N, a number of (total) angular momentum states were found to be stable at successively higher critical angular velocities ${\ensuremath{\Omega}}_{\mathrm{ci}}, i=1,2,3,\dots{}.$ All the states in the regime N>Lz>0 are metastable. For Lz>N, the Lz values for the stable ground states generally increased with increasing critical angular velocities ${\ensuremath{\Omega}}_{\mathrm{ci}},$ and the condensate was strongly depleted. The critical ${\ensuremath{\Omega}}_{\mathrm{ci}}$ values, however, decreased with increasing interaction strength as well as the particle number, and were systematically greater than the nonvariational yrast-state values for the Lz=N single vortex state. We have also observed that the condensate fraction for the single vortex state (as also for the higher vortex states) did not change significantly even as the two-body interaction strength was varied over several $(\ensuremath{\sim}4)$ orders of magnitude in the moderately to the weakly interacting regime.
Bose gases in rotating optical lattices combine two important topics in quantum physics: superfluid rotation and strong correlations. In this paper, we examine square two-dimensional systems at zero temperature comprised of strongly repulsive bosons with filling factors of up to one atom per lattice site. The entry of vortices into the system is characterized by jumps of $2\ensuremath{\pi}$ in the phase winding of the condensate wave function. A lattice of size $L\ifmmode\times\else\texttimes\fi{}L$ can have at most $L\ensuremath{-}1$ quantized vortices in the lowest Bloch band. In contrast to homogeneous systems, angular momentum is not a good quantum number since the continuous rotational symmetry is broken by the lattice. Instead, a quasiangular momentum captures the discrete rotational symmetry of the system. Energy level crossings indicative of quantum phase transitions are observed when the quasiangular momentum of the ground state changes.
In this thesis, I explore the behavior of rotating ultra-cold Bose gases, by diagonalizing the Hamiltonian. This method has the advantage of being exact in the limit of weak interactions, and thus is a useful complement to the more commonly used Gross-Pitaevskii (GP) equation, which relies on a mean-field approximation. It is generally known that an ultra-cold Bose gas shows some remarkable properties under rotation, such as quantized vortices and persistent currents, and these are one main reason why we choose to study rotation. Abstract A repulsively interacting Bose-Einstein condensate (BEC) in a harmonic trapping potential will form singly quantized vortices when rotated. Papers I-V are examples of different kinds of generalizations to this basic system, each showing its particular non-trivial new features. Abstract In Paper I we make the trapping potential weakly anharmonic, and confirm the existence of multiply quantized vortices in both repulsively and attractively interacting condensates, which were previously predicted in Gross-Pitaevskii calculations. Abstract In Papers II and III, we generalize to a two-component system in a harmonic trap. The vortices that form under rotation will in general be coreless, i.e. one component will form a vortex and rotate around the other component. In Paper II we give some analytical expressions for the dispersion relation and occupation numbers at low angular momenta, and in Paper III we show how multiply quantized vortices can be stable in these systems. Abstract Finally we investigate bosons in potentials which can support persistent currents, i.e. rotating states that do not decay in a finite amount of time. These are possible if there is an energy barrier separating the rotating state from the non-rotating ground state. In Paper IV, we have a two-component condensate in a one-dimensional potential and we calculate among other things the minimum interaction strength where persistent currents can occur. In Paper V, we study an experimentally relevant two-dimensional annulus. Again we find that persistent currents are possible above a certain interaction strength, and we find this strength as a function of the width of the annulus, as well as how it depends on the relative population of the two species.
In this dissertation, we study rotational properties of ultracold bosonic quantum gases in two trapped configurations, a quasi-two-dimensional gas in a harmonic trap and a quasi-one-dimensional gas in a toroidal trap. First, we investigate the effects of correlations on the properties of the ground state of the rotating harmonically-trapped Bose gas by adding Bogoliubov fluctuations to the mean-field ground state of an N-particle single-vortex system. We demonstrate that the fluctuation-induced correlations lower the energy compared to that of the mean-field ground state, that the vortex core is pushed slightly away from the center of the trap, and that an unstable mode with negative energy (for rotations slower than a critical frequency) emerges in the energy spectrum, thus, pointing to a better state for slow rotation. We construct mean-field ground states of 0-, 1-, and 2-vortex states as a function of rotation rate and determine the critical frequencies for transitions between these states, as well as the critical frequency for appearance of a metastable state with an off-center vortex and its image vortex in the evanescent tail of the cloud. Then, we show how the configuration-space form of the above-mentioned Bogoliubov ground-state wave function of a bosonic condensate with a single vortex in a harmonic trap can be described in terms of bosonic Jastrow correlations. We then generalize this result to study the effects of such correlations on a mean-field vortex lattice state and show that the included correlations lower the energy below that of the mean-field state. Although the reduction is relatively small, it is a precursor of the more general expected effect of correlations in describing the melting of the vortex lattice at a high angular momentum per particle. Finally, we study the stability and dynamics of an ultracold bosonic gas trapped in a toroidal geometry and driven by rotation in the absence of dissipation. We first delineate, via the Bogoliubov mode expansion, the regions of stability and the nature of instabilities of the system for both repulsive and attractive interaction strengths. To study the response of the system to variations in the rotation rate, we introduce a “disorder” potential, breaking the rotational symmetry. We demonstrate the breakdown of adiabaticity as the rotation rate is slowly varied and find forced tunneling between the eigenstates of the system. The nonadiabaticity is signaled by the appearance of a swallow-tail loop in the lowest-energy level, a general sign of hysteresis. Then, we show that this system is in one-to-one correspondence with a trapped gas in a double-well potential and thus exhibits macroscopic quantum self-trapping. Finally, we show that self-trapping is a direct manifestation of the behavior of the lowest-energy level.
We present exact diagonalization study on a system of 10≤N≤24 spinless bosons interacting via repulsive Gaussian potential, harmonically confined in xy-plane with an externally impressed rotation about the z-axis. The two-body interaction strength in the Gaussian potential is taken in the strongly interacting regime with values of interaction range in the regime 0≤σ≤1. The diagonalization of the N-body Hamiltonian matrix, in subspaces of total angular momentum in the regime 0≤Lz≤4N corresponds to the filling fraction ν≲3.2 is carried out to obtain the variationally exact ground-state wavefunction and the corresponding eigenvalue. It is found that an increase in interaction range σ leads to (a) a systematic decrease in energy, (b) an increase in the critical angular velocity Ωci of the i{th} vortex state and (c) an increase in the largest eigenvalue λ1, (condensate fraction), of one-particle reduced density matrix (OPRDM). The von Neumann entropy S1(Lz,σ), quantifying the quantum entanglement between the particles in the many-body ground state, is largely found to decrease with increase in σ. Crossings in von Neumann entropy for several of the angular momentum states are observed with variation in σ. The response of the Bose-Einstein condensate to rotation is examined through Lz(σ)−Ω(σ) stability graph for several values of σ. A vortex state with larger plateau length on the Lz−Ω graph is considered to be more stable. The internal structure of the condensate, as depicted by the conditional probability distribution (CPD) in the body-fixed frame, exhibits characteristic features with interaction range σ. One such feature is the merging of the cores of the two-vortex state.