A rotating bosonic many-body system in a harmonic trap is studied with the three-dimensional cranked Hartree-Fock-Bogoliubov method at zero temperature, which has been applied to nuclear many-body systems at high spin. This method is a variational method extended from Hartree-Fock theory, which can treat the pairing correlations in a self-consistent manner. An advantage of this method is that a finite-range interaction between constituent particles can be used in the calculation, unlike the original Gross-Pitaevskii approach. To demonstrate the validity of our method, we present a calculation for a toy model---that is, a rotating system of ten bosonic particles interacting through the repulsive quadrupole-quadrupole interaction in a harmonic trap. It is found that the yrast states, the lowest-energy states for the given total angular momentum, do not correspond to the Bose-Einstein condensate, except for a few special cases. One such case is a vortex state, which appears when the total angular momentum
L is twice the particle number
N (i.e.,
L=2N).
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.
We derive exact analytical results for the many-body wave functions and energies of harmonically trapped two-component Bose-Einstein condensates with weak repulsive interactions under rotation. The wave functions are universal as they do not depend on the matrix elements of the two-body interaction. A comparison with the results from numerical diagonalization shows that the analytical solutions describe the ground state, as well as a set of low-lying excited states for repulsive contact interactions, Coulomb interactions, and the repulsive interactions between aligned dipoles.
Using Perron-Frobenius theorem, we prove that the results by Wilkin, Gunn and Smith [1] for the ground states of N Bose atoms rotating at the angular momentum L in a harmonic atomic trap with frequency omega interacting via attractive delta^2(r) forces, are valid for a broad class of predominantly attractive interactions V(r), not necessarily attractive for any r. The sufficient condition for the interaction is that all the two-body matrix elements allowed by the conservation of angular momentum k+l = m+n, are negative. This class includes, in particular, the Gaussian attraction of arbitrary radius, -1/r - Coulomb and log(r)-Coulomb forces, as well as all the short-range R << omega^{-1/2} interactions satisfying inequality int d^2r V(r) < 0. There is no condensation at L>> 1, and the angular momentum is concentrated in the collective ``center-of-mass'' mode.
In a k-dimensional system of weakly interacting Bose atoms trapped by a spherically symmetric and harmonic external potential, an exact expression is obtained for the rotating ground states at a fixed angular momentum. The result is valid for arbitrary interactions obeying minimal physical requirements. Depending on the sign of a modified scattering length, it reduces to either a collective rotation or a condensed vortex state, with no alternative. The ground state can undergo a kind of quantum phase transition when the shape of the interaction potential is smoothly varied.