Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0701333 · Physical Review A

Rapidly rotating boson molecules with long- or short-range repulsion: An exact diagonalization study

Leslie O. Baksmaty, Constantine Yannouleas, Uzi Landman

Published Feb 26, 2007• 31 citations
arXivHTMLPDF

Abstract

The Hamiltonian for a small number $N\ensuremath{\leqslant}11$ of bosons in a rapidly rotating harmonic trap, interacting via a short range (contact potential) or a long range (Coulomb) interaction, is studied via an exact diagonalization in the lowest Landau level. Our analysis shows that, for both low and high fractional fillings, the bosons localize and form rotating boson molecules (RBMs) consisting of concentric polygonal rings. Focusing on systems with the number of trapped atoms sufficiently large to form multiring bosonic molecules, we find that, as a function of the rotational frequency and regardless of the type of repulsive interaction, the ground-state angular momenta grow in specific steps that coincide with the number of localized bosons on each concentric ring. Comparison of the conditional probability distributions (CPDs) for both interactions suggests that the degree of crystalline correlations appears to depend more on the fractional filling $\ensuremath{\nu}$ than on the range of the interaction. The RBMs behave as nonrigid rotors, i.e., the concentric rings rotate independently of each other. At filling fractions $\ensuremath{\nu}<1∕2$, we observe well developed crystallinity in the CPDs (two-point correlation functions). For larger filling fractions $\ensuremath{\nu}>1∕2$, observation of similar molecular patterns requires consideration of even higher-order correlation functions.

DOI
10.1103/physreva.75.023620

Related papers

Bosonic Molecules in Rotating Traps

Igor Romanovsky, Constantine Yannouleas, Leslie O. Baksmaty +1 more

We present a variational many-body wave function for repelling bosons in rotating traps, focusing on rotational frequencies that do not lead to restriction to the lowest Landau level. This wave function incorporates correlations beyond the Gross-Pitaevskii (GP) mean-field approximation, and it describes rotating boson molecules (RBMs) made of localized bosons that form polygonal-ring-like crystalline patterns in their intrinsic frame of reference. The RBMs exhibit characteristic periodic dependencies of the ground-state angular momenta on the number of bosons in the polygonal rings. For small numbers of neutral bosons, the RBM ground-state energies are found to be always lower than those of the corresponding GP solutions, in particular, in the regime of GP vortex formation.

arXiv:cond-mat/0608571•Physical Review Letters•Aug 30, 2006•28 citations

Exact diagonalization study of energy-level statistics in harmonically confined interacting bosons

Mohd Talib, M. A. H. Ahsan

We present an exact diagonalization study of the spectral properties of bosons harmonically confined in a quasi-2D plane and interacting via repulsive Gaussian potential. We consider the lowest 100 energy levels for systems of N=12,16, and 20 bosons in two distinct regimes: (a) when the interaction energy is small compared to the trap energy (moderate interaction) and (b) when the interaction energy is comparable to the trap energy (strong interaction), for the nonrotating (L_{z}=0) as well as the rotating single-vortex state (L_{z}=N). For higher angular momenta, L_{z}=2N and L_{z}=3N, only the strong interaction regime is considered. While the nearest-neighbor spacing distribution (NNSD) P(s) and the ratios of consecutive level spacings distribution P(r) are used to study the short-range correlations, the Dyson-Mehta Δ_{3} statistic and the level number variance Σ^{2}(L) are used to examine the long-range correlations. In the moderate interaction regime, the nonrotating system exhibits Poisson distribution, a characteristic of the regular energy spectra. In the strong interaction regime, the nonrotating system exhibits chaotic behavior signified by GOE distribution. Furthermore, in the rotating case for the single-vortex state (L_{z}=N) in the moderate interaction regime, the system exhibits signatures of weak chaos with some degree of regularity in the energy-level spectra. However, in the strong interaction regime for the rotating case with L_{z}=N, 2N and 3N, the system exhibits strong chaotic behavior. The rotation is found to contribute to enhancement of chaotic behavior in the system for both the moderate and the strong interaction regimes. Our results of NNSD analysis are supported by the analysis of the ratios of consecutive level spacings distribution P(r), which does not involve unfolding.

arXiv:2508.01729•Physical review. E•Jan 8, 2026•2 citations

Effect of Finite-Range Interactions on Rapidly Rotating Ultracold Bosonic Atoms

Nobukuni Hamamoto

We investigate the effects of the finite-range interactions of six rotating ultracold bosonic atoms using a Gaussian-type interatomic interaction model. The model is analyzed numerically by exact diagonalization within the Lowest Landau Level (LLL) approximation and semiclassical approximation. The result of exact diagonalization shows that the ground-state angular momentum changes discretely with increasing angular velocity. For the short-range limit, the ground-state angular momentum and wavefunctions agree with those of the delta interaction evaluated by Bertsch and Papenbrock [Phys. Rev. A 63, 023616 (2001)]. Different from the delta interaction, the ground-state angular momenta higher than 30, i.e., N(N − 1), are observed at a high angular frequency as a result of the finite-range two-body interactions. For the intermediate-range interaction, the sequence of ground-state angular momenta increases in steps of five, which was not found in previous works on the Gaussian interaction. For the long-range limit of Gaussian interaction, we find that the ground-state angular momenta increase in steps of six. These steps of the ground-state angular momentum according to the width of the Gaussian interactions are explained by semiclassical and classical analysis based on the rovibrating molecule picture. The increments of the ground-state angular momentum of five and six are explained by the semiclassical quantization condition of the rotational and vibrational modes of fivefold and sixfold molecules, respectively. Our analysis based on the classical model also confirms that the fivefold molecule picture is more stable than the sixfold molecule picture in the intermediate range of the Gaussian interaction. These results suggest that the Gaussian interaction model can be used to emulate and characterize interactions by their width as the model can reproduce various rotational states including the ground states of rotating systems interacting with delta-type and Coulomb type interactions.

Journal of the Physical Society of Japan•Nov 15, 2017•0 citations

Ground state properties of trapped boson system with finite-range Gaussian repulsion: exact diagonalization study

Mohd. Imran, M. A. H. Ahsan

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.

arXiv:1602.03811•Journal of Physics B Atomic Molecular and Optical Physics•Mar 13, 2020•3 citations

Exact solutions for interacting boson systems under rotation

T. Papenbrock, G. F. Bertsch

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

•Jan 1, 2000•6 citations

Crystalline Boson Phases in Harmonic Traps: Beyond the Gross-Pitaevskii Mean Field

Igor Romanovsky, Constantine Yannouleas, Uzi Landman

Strongly-interacting bosons in two-dimensional harmonic traps are described through breaking of rotational symmetry at the Hartree-Fock level and subsequent symmetry restoration via projection techniques, thus incorporating correlations beyond the Gross-Pitaevskii (GP) solution. The bosons localize and form polygonal-ringlike crystalline patterns, both for a repulsive contact potential and a Coulomb interaction, as revealed via conditional-probability-distribution analysis. For neutral bosons, the total energy of the crystalline phase saturates in contrast to the GP solution, and its spatial extent becomes smaller than that of the GP condensate. For charged bosons, the total energy and dimensions approach the values of classical pointlike charges in their equilibrium configuration.

arXiv:cond-mat/0410598•Physical Review Letters•Dec 1, 2004•51 citations
The Hamiltonian for a small number N\ensuremath⩽11N\ensuremath{\leqslant}11N\ensuremath⩽11 of bosons in a rapidly rotating harmonic trap, interacting via a short range (contact potential) or a long range (Coulomb) interaction, is studied via an exact diagonalization in the lowest Landau level. Our analysis shows that, for both low and high fractional fillings, the bosons localize and form rotating boson molecules (RBMs) consisting of concentric polygonal rings. Focusing on systems with the number of trapped atoms sufficiently large to form multiring bosonic molecules, we find that, as a function of the rotational frequency and regardless of the type of repulsive interaction, the ground-state angular momenta grow in specific steps that coincide with the number of localized bosons on each concentric ring. Comparison of the conditional probability distributions (CPDs) for both interactions suggests that the degree of crystalline correlations appears to depend more on the fractional filling \ensuremathν\ensuremath{\nu}\ensuremathν than on the range of the interaction. The RBMs behave as nonrigid rotors, i.e., the concentric rings rotate independently of each other. At filling fractions \ensuremathν<1∕2\ensuremath{\nu}<1∕2\ensuremathν<1∕2, we observe well developed crystallinity in the CPDs (two-point correlation functions). For larger filling fractions \ensuremathν>1∕2\ensuremath{\nu}>1∕2\ensuremathν>1∕2, observation of similar molecular patterns requires consideration of even higher-order correlation functions.
We present a variational many-body wave function for repelling bosons in rotating traps, focusing on rotational frequencies that do not lead to restriction to the lowest Landau level. This wave function incorporates correlations beyond the Gross-Pitaevskii (GP) mean-field approximation, and it describes rotating boson molecules (RBMs) made of localized bosons that form polygonal-ring-like crystalline patterns in their intrinsic frame of reference. The RBMs exhibit characteristic periodic dependencies of the ground-state angular momenta on the number of bosons in the polygonal rings. For small numbers of neutral bosons, the RBM ground-state energies are found to be always lower than those of the corresponding GP solutions, in particular, in the regime of GP vortex formation.
We present an exact diagonalization study of the spectral properties of bosons harmonically confined in a quasi-2D plane and interacting via repulsive Gaussian potential. We consider the lowest 100 energy levels for systems of N=12,16, and 20 bosons in two distinct regimes: (a) when the interaction energy is small compared to the trap energy (moderate interaction) and (b) when the interaction energy is comparable to the trap energy (strong interaction), for the nonrotating (L_{z}=0) as well as the rotating single-vortex state (L_{z}=N). For higher angular momenta, L_{z}=2N and L_{z}=3N, only the strong interaction regime is considered. While the nearest-neighbor spacing distribution (NNSD) P(s) and the ratios of consecutive level spacings distribution P(r) are used to study the short-range correlations, the Dyson-Mehta Δ_{3} statistic and the level number variance Σ^{2}(L) are used to examine the long-range correlations. In the moderate interaction regime, the nonrotating system exhibits Poisson distribution, a characteristic of the regular energy spectra. In the strong interaction regime, the nonrotating system exhibits chaotic behavior signified by GOE distribution. Furthermore, in the rotating case for the single-vortex state (L_{z}=N) in the moderate interaction regime, the system exhibits signatures of weak chaos with some degree of regularity in the energy-level spectra. However, in the strong interaction regime for the rotating case with L_{z}=N, 2N and 3N, the system exhibits strong chaotic behavior. The rotation is found to contribute to enhancement of chaotic behavior in the system for both the moderate and the strong interaction regimes. Our results of NNSD analysis are supported by the analysis of the ratios of consecutive level spacings distribution P(r), which does not involve unfolding.
We investigate the effects of the finite-range interactions of six rotating ultracold bosonic atoms using a Gaussian-type interatomic interaction model. The model is analyzed numerically by exact diagonalization within the Lowest Landau Level (LLL) approximation and semiclassical approximation. The result of exact diagonalization shows that the ground-state angular momentum changes discretely with increasing angular velocity. For the short-range limit, the ground-state angular momentum and wavefunctions agree with those of the delta interaction evaluated by Bertsch and Papenbrock [Phys. Rev. A 63, 023616 (2001)]. Different from the delta interaction, the ground-state angular momenta higher than 30, i.e., N(N − 1), are observed at a high angular frequency as a result of the finite-range two-body interactions. For the intermediate-range interaction, the sequence of ground-state angular momenta increases in steps of five, which was not found in previous works on the Gaussian interaction. For the long-range limit of Gaussian interaction, we find that the ground-state angular momenta increase in steps of six. These steps of the ground-state angular momentum according to the width of the Gaussian interactions are explained by semiclassical and classical analysis based on the rovibrating molecule picture. The increments of the ground-state angular momentum of five and six are explained by the semiclassical quantization condition of the rotational and vibrational modes of fivefold and sixfold molecules, respectively. Our analysis based on the classical model also confirms that the fivefold molecule picture is more stable than the sixfold molecule picture in the intermediate range of the Gaussian interaction. These results suggest that the Gaussian interaction model can be used to emulate and characterize interactions by their width as the model can reproduce various rotational states including the ground states of rotating systems interacting with delta-type and Coulomb type interactions.
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 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
Strongly-interacting bosons in two-dimensional harmonic traps are described through breaking of rotational symmetry at the Hartree-Fock level and subsequent symmetry restoration via projection techniques, thus incorporating correlations beyond the Gross-Pitaevskii (GP) solution. The bosons localize and form polygonal-ringlike crystalline patterns, both for a repulsive contact potential and a Coulomb interaction, as revealed via conditional-probability-distribution analysis. For neutral bosons, the total energy of the crystalline phase saturates in contrast to the GP solution, and its spatial extent becomes smaller than that of the GP condensate. For charged bosons, the total energy and dimensions approach the values of classical pointlike charges in their equilibrium configuration.