Abstract The ground and low-lying collective states of a rotating system of N = 3 bosons harmonically confined in quasi-two-dimension and interacting via repulsive finite-range Gaussian potential is studied in weakly to moderately interacting regime. The N-body Hamiltonian matrix is diagonalized in subspaces of quantized total angular momenta 0 ≥ L ≥ 4N to obtain the ground and low-lying eigenstates. Our numerical results show that breathing modes with N-body eigenenergy spacing of 2ħω⊥, known to exist in strictly 2D system with zero-range (δ-function) interaction potential, may as well exist in quasi-2D system with finite-range Gaussian interaction potential. To gain an insight into the many-body states, the von Neumann entropy is calculated as a measure of quantum correlation and the conditional probability distribution is analyzed for the internal structure of the eigenstates. In the rapidly rotating regime the ground state in angular momentum subspaces L = (q/2)N (N − 1) with q = 2, 4 is found to exhibit the anticorrelation structure suggesting that it may variationally be described by a Bose-Laughlin like state. We further observe that the first breathing mode exhibits features similar to the Bose-Laughlin state in having eigenenergy, von Neumann entropy and internal structure independent of interaction for the three-boson system considered here. On the contrary, for eigenstates lying between the Bose-Laughlin like ground state and the first breathing mode, values of eigenenergy, von Neumann entropy and internal structure are found to vary with interaction.
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 investigate a system of
N spinless bosons confined in quasi-two-dimensional harmonic trap with repulsive two-body finite-range Gaussian interaction potential of large
s-wave scattering length. Exact diagonalization of the Hamiltonian matrix is carried out to obtain the
N-body ground state as well as low-lying excited states, using Davidson algorithm in beyond lowest-Landau-level approximation. We examine the finite-range effects of the interaction potential on the many-body ground state energy as also the degree of condensation of the Bose-condensed gas. The results obtained indicate that the finite-range Gaussian interaction potential enhances the degree of condensation compared to the zero-range interaction potential. We further analyze the effect of finite-range interaction potential on the breathing mode collective excitation. Our theoretical results may be relevant for experiments currently conducted on quasi-two-dimensional Bose gas with more realistic interaction potential.
We investigate the ground and low excited states of a rotating, weakly interacting Bose-Einstein condensed gas in a harmonic trap for a given angular momentum. Analytical results in various limits as well as numerical results are presented, and these are compared with those of previous studies. Within the mean-field approximation and for repulsive interaction between the atoms, we find that for very low values of the total angular momentum per particle, $L/\stackrel{\ensuremath{\rightarrow}}{N}0,$ where $L\ensuremath{\Elzxh}$ is the angular momentum and N is the total number of particles, the angular momentum is carried by quadrupole
(∣m∣=2) surface modes. For
L/N=1, a vortexlike state is formed and all the atoms occupy the
m=1 state. For small negative values of $L/N\ensuremath{-}1,$ the states with
m=0 and
m=2 become populated, and for small positive values of $L/N\ensuremath{-}1,$ atoms in the states with
m=5 and
m=6 carry the additional angular momentum. In the whole region
0< L/N< 1, we have strong analytic and numerical evidence that the interaction energy drops linearly as a function of
L/N. We have also found that an array of singly quantized vortices is formed as
L/N increases. Finally, we have gone beyond the mean-field approximation and have calculated the energy of the lowest state up to order N for small negative values of $L/N\ensuremath{-}1,$ as well as the energy of the low-lying excited states.
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.
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.