A system of N two-dimensional weakly interacting bosons in a harmonic trap is considered. When the two-particle potential is a \ensuremath{\delta} function, Smith and Wilkin have proved analytically that the elementary symmetric polynomials of particle coordinates measured from the center of mass are exact eigenstates. In this study, we point out that their proof works equally well for an arbitrary two-particle potential that possesses translational and rotational symmetries. We find that the interaction energy associated with the eigenstate with angular momentum L is equal to $aN(N\ensuremath{-}1)/2+(b\ensuremath{-}a)NL/2,$ where a and b are the interaction energies of two bosons in the lowest-energy one-particle state with zero and one unit of angular momentum, respectively. Additionally, we study briefly the case of attractive quartic interactions. We prove rigorously that the lowest-energy state is the one in which all angular momentum is carried by the center of mass motion.
We consider a model of N two-dimensional bosons in a harmonic trap with translational and rotational invariant, weak two-particle interaction. We present in configuration space a systematical recursive method for constructing all wave functions with angular momentum L and corresponding energies and apply it to
L< 6 for all N. The lower and the upper bounds for interaction energy are estimated. We analytically confirm the conjecture of Smith et al. that the elementary symmetric polynomial is the ground state for repulsive \ensuremath{\delta} interaction for all
N> L up to
L< 6. Additionally, we find that there exist vanishing-energy solutions for $L>~N(N\ensuremath{-}1),$ signaling the exclusive statistics. Finally, we consider briefly the case of attractive powerlike potential
rk, $k>\ensuremath{-}2,$ and prove that the lowest-energy state is still the one in which the whole angular momentum is absorbed by the center-of-mass motion.
We study a class of interacting, harmonically trapped boson systems at angular momentum L. The Hamiltonian leaves a L-dimensional subspace invariant, and this permits an explicit solution of several eigenstates and energies for a wide class of two-body interactions. PACS numbers: 05.30.Jp, 03.75.Fi, 03.65.Fd, 67.40.Db Typeset using REVTEX 1 The study of vortices in Bose-Einstein condensates of trapped atomic vapors is of considerable experimental [1,2] and theoretical [3–13] interest. Present experiments and much theoretical work [3–7] focus on the Thomas-Fermi regime of short coherence length [14]. The opposite limit of perturbatively weak contact interactions is not yet realized experimentally but is of theoretical interest as well [8–13]. The case of attractive interactions was studied by Wilkin et al. [9]. Mottelson developed a theory for repulsive interactions [10], and Kavoulakis et al. [12] and Jackson et al. [13] compared mean field and exact numerical results. Exact diagonalization techniques by us showed that the ground state energy depends linearly on angular momentum and led to an analytical expression for the ground state wave function [11]. Recently, the form of this wave function and its eigenenergy were confirmed in analytic calculations [15–17]. However, we still lack insight about why certain eigenstates turn out to be simple analytical functions. We show in this paper that there is a subspace structure that explains these findings. Interestingly, this structure is not limited to the two-body contact interaction considered in previous investigations but present in a rather large class of Hamiltonians. Let us consider a system of N harmonically trapped bosons in two spatial dimensions 1 at angular momentum
For a system of N bosons in a two-dimensional harmonic trap with frequency $\ensuremath{\omega},$ interacting via forces $V\ensuremath{\ll}\ensuremath{\Elzxh}\ensuremath{\omega},$ we develop a method to find the ground states of rotating Bose condensate as a function of two quantum numbers, the total angular momentum and the angular momentum of internal excitations (generalized yrast states). The energies of these condensed vortex states are expressed through the single two-body matrix element of interaction,
V. A broad universality class of predominantly repulsive interactions for which these results hold is described by a simple integral condition on
V. It includes Gaussian, $\ensuremath{\delta}$ function, and log-Coulomb forces.
Abstract We use exact diagonalization to study an interacting system of N spinless bosons with finite-range Gaussian repulsion, confined in a quasi-two-dimensional harmonic trap with and without an introduced rotation. The diagonalization of the Hamiltonian matrix using Davidson algorithm in subspaces of quantized total angular momentum L z is carried out to obtain the N -body lowest eigenenergy and eigenstate. To bring out the effect of quantum (Bose) statistics and consequent phase stiffness (rigidity) of the variationally obtained many-body wavefunction on various physical quantities, our study spans from few-body ( N = 2) to many-body ( N = 16) systems. Further, to examine the finite-range effect of the repulsive Gaussian potential on many-body ground state properties of the Bose-condensate, we obtain the lowest eigenstate, the critical angular velocity of single vortex state and the quantum correlation (measured) in terms of von Neumann entanglement entropy and degree of condensation. It is found that for small values of the range (measured by the parameter σ ) of Gaussian potential, the ground state energy increases for few-boson (2 ⩽ N ⩽ 8) systems but decreases for many-boson ( N > 8) systems. On the other hand for relatively large values of the range of Gaussian potential, the ground state energy exhibits a monotonic decrease, regardless of the number of bosons N . For a given N , there is found an optimal value of the range of Gaussian potential for which the first vortex (with L z = N ) nucleates at a lower value of the rotational angular velocity Ω c1 compared to the zero-range ( δ -function) potential. Further, we observe that the inter-particle interaction and the introduced rotation are competing effects with latter being dominant over the former. With increase in the range of Gaussian potential, the value of von Neumann entropy decreases and the degree of condensation increases implying an enhanced quantum correlation and phase rigidity.
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.