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:1903.04974 · Physical Review A

Energy-level splitting for weakly interacting bosons in a harmonic trap

Ben Craps, Marine De Clerck, Oleg Evnin, Surbhi Khetrapal

Published Aug 6, 2019• 14 citations
arXivHTMLPDF

Abstract

We consider identical quantum bosons with weak contact interactions in a two-dimensional isotropic harmonic trap. When the interactions are turned off, the energy levels are equidistant and highly degenerate. At linear order in the coupling parameter, these degenerate levels split, and we study the patterns of this splitting. It turns out that the problem is mathematically identical to diagonalizing the quantum resonant system of the two-dimensional Gross-Pitaevskii equation, whose classical counterpart has been studied previously in the mathematical literature on turbulence. Our purpose is to explore the implications of the symmetries and energy bounds of this resonant system, previously studied for the classical case, for the quantum level splitting. Simplifications in computing the splitting spectrum numerically result from exploiting the symmetries. The highest energy state emanating from each unperturbed level is explicitly described by our analytics. We furthermore discuss the energy-level spacing distributions in the spirit of quantum chaos theory. After separating the eigenvalues into blocks with respect to the known conservation laws, we observe the Wigner-Dyson statistics within specific large blocks, which leaves little room for further integrable structures in the problem beyond the symmetries that are already explicitly known.

DOI
10.1103/physreva.100.023605

Related papers

Phase-space localization at the lowest Landau level

Ben Craps, Marine De Clerck, Oleg Evnin +1 more

We consider bosons with weak contact interactions in a harmonic trap and focus on states at the lowest Landau level. Motivated by the known nontrivial phase-space topography of the energy functional of the corresponding Gross-Pitaevskii equation, we explore Husimi distributions of quantum energy eigenstates in the classical phase space of the Schr\"odinger field. With interactions turned off, the energy levels are highly degenerate and the Husimi distributions do not manifest any particular localization properties. With interactions turned on, the degeneracy is lifted, and a selection of energy levels emerges whose Husimi distributions are localized around low-dimensional surfaces in the phase space.

arXiv:2407.07675•Physical Review A•Dec 27, 2024•2 citations

Energy-level statistics of interacting trapped bosons

Barnali Chakrabarti, Anindya Biswas, V. K. B. Kota +2 more

It is a well-established fact that statistical properties of energy-level spectra are the most efficient tool to characterize nonintegrable quantum systems. The statistical behavior of different systems such as complex atoms, atomic nuclei, two-dimensional Hamiltonians, quantum billiards, and noninteracting many bosons has been studied. The study of statistical properties and spectral fluctuations in interacting many-boson systems has developed interest in this direction. We are especially interested in weakly interacting trapped bosons in the context of Bose-Einstein condensation (BEC) as the energy spectrum shows a transition from a collective nature to a single-particle nature with an increase in the number of levels. However this has received less attention as it is believed that the system may exhibit Poisson-like fluctuations due to the existence of an external harmonic trap. Here we compute numerically the energy levels of the zero-temperature many-boson systems which are weakly interacting through the van der Waals potential and are confined in the three-dimensional harmonic potential. We study the nearest-neighbor spacing distribution and the spectral rigidity by unfolding the spectrum. It is found that an increase in the number of energy levels for repulsive BEC induces a transition from a Wigner-like form displaying level repulsion to the Poisson distribution for $P(s)$. It does not follow the Gaussian orthogonal ensemble prediction. For repulsive interaction, the lower levels are correlated and manifest level-repulsion. For intermediate levels $P(s)$ shows mixed statistics, which clearly signifies the existence of two energy scales: external trap and interatomic interaction, whereas for very high levels the trapping potential dominates, generating a Poisson distribution. Comparison with mean-field results for lower levels are also presented. For attractive BEC near the critical point we observe the Shnirelman-like peak near $s=0$, which signifies the presence of a large number of quasidegenerate states.

arXiv:1101.5469•Physical Review A•Jul 23, 2012•10 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

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

Spin-orbit-coupled bosons interacting in a two-dimensional harmonic trap

Pere Mujal, A. Polls, Bruno Juliá-Díaz

A system of bosons in a two-dimensional harmonic trap in the presence of Rashba-type spin-orbit coupling is investigated. An analytic treatment of the ground state of a single atom in the weak-coupling regime is presented and used as a basis for a perturbation theory in the interacting two-boson system. The numerical diagonalization of both the single-particle and the two-boson Hamiltonian matrices allows us to go beyond those approximations and obtain not only the ground state, but also the low-energy spectra and the different energy contributions separately. We show that the expectation value of the spin-orbit term is related to the expectation value of ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\sigma}}}_{z}{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{L}}_{z}$ for the eigenstates of the system, regardless of the trapping potential. The low-energy states of the repulsively interacting two-boson system are characterized. In the presence of a sufficiently strong interaction and spin-orbit coupling strength, there is a direct energy-level crossing in the ground state of the system between states of different ${J}_{z}$, the third component of the total angular momentum, that changes its structure. This is reflected in a discontinuity in the different energy terms and it is signaled in the spatial density of the system.

arXiv:1907.13355•Physical Review A•Apr 24, 2020•9 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
We consider identical quantum bosons with weak contact interactions in a two-dimensional isotropic harmonic trap. When the interactions are turned off, the energy levels are equidistant and highly degenerate. At linear order in the coupling parameter, these degenerate levels split, and we study the patterns of this splitting. It turns out that the problem is mathematically identical to diagonalizing the quantum resonant system of the two-dimensional Gross-Pitaevskii equation, whose classical counterpart has been studied previously in the mathematical literature on turbulence. Our purpose is to explore the implications of the symmetries and energy bounds of this resonant system, previously studied for the classical case, for the quantum level splitting. Simplifications in computing the splitting spectrum numerically result from exploiting the symmetries. The highest energy state emanating from each unperturbed level is explicitly described by our analytics. We furthermore discuss the energy-level spacing distributions in the spirit of quantum chaos theory. After separating the eigenvalues into blocks with respect to the known conservation laws, we observe the Wigner-Dyson statistics within specific large blocks, which leaves little room for further integrable structures in the problem beyond the symmetries that are already explicitly known.
We consider bosons with weak contact interactions in a harmonic trap and focus on states at the lowest Landau level. Motivated by the known nontrivial phase-space topography of the energy functional of the corresponding Gross-Pitaevskii equation, we explore Husimi distributions of quantum energy eigenstates in the classical phase space of the Schr\"odinger field. With interactions turned off, the energy levels are highly degenerate and the Husimi distributions do not manifest any particular localization properties. With interactions turned on, the degeneracy is lifted, and a selection of energy levels emerges whose Husimi distributions are localized around low-dimensional surfaces in the phase space.
It is a well-established fact that statistical properties of energy-level spectra are the most efficient tool to characterize nonintegrable quantum systems. The statistical behavior of different systems such as complex atoms, atomic nuclei, two-dimensional Hamiltonians, quantum billiards, and noninteracting many bosons has been studied. The study of statistical properties and spectral fluctuations in interacting many-boson systems has developed interest in this direction. We are especially interested in weakly interacting trapped bosons in the context of Bose-Einstein condensation (BEC) as the energy spectrum shows a transition from a collective nature to a single-particle nature with an increase in the number of levels. However this has received less attention as it is believed that the system may exhibit Poisson-like fluctuations due to the existence of an external harmonic trap. Here we compute numerically the energy levels of the zero-temperature many-boson systems which are weakly interacting through the van der Waals potential and are confined in the three-dimensional harmonic potential. We study the nearest-neighbor spacing distribution and the spectral rigidity by unfolding the spectrum. It is found that an increase in the number of energy levels for repulsive BEC induces a transition from a Wigner-like form displaying level repulsion to the Poisson distribution for P(s)P(s)P(s). It does not follow the Gaussian orthogonal ensemble prediction. For repulsive interaction, the lower levels are correlated and manifest level-repulsion. For intermediate levels P(s)P(s)P(s) shows mixed statistics, which clearly signifies the existence of two energy scales: external trap and interatomic interaction, whereas for very high levels the trapping potential dominates, generating a Poisson distribution. Comparison with mean-field results for lower levels are also presented. For attractive BEC near the critical point we observe the Shnirelman-like peak near s=0s=0s=0, which signifies the presence of a large number of quasidegenerate states.
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.
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.
A system of bosons in a two-dimensional harmonic trap in the presence of Rashba-type spin-orbit coupling is investigated. An analytic treatment of the ground state of a single atom in the weak-coupling regime is presented and used as a basis for a perturbation theory in the interacting two-boson system. The numerical diagonalization of both the single-particle and the two-boson Hamiltonian matrices allows us to go beyond those approximations and obtain not only the ground state, but also the low-energy spectra and the different energy contributions separately. We show that the expectation value of the spin-orbit term is related to the expectation value of ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\sigma}}}_{z}{\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{L}}_{z}$ for the eigenstates of the system, regardless of the trapping potential. The low-energy states of the repulsively interacting two-boson system are characterized. In the presence of a sufficiently strong interaction and spin-orbit coupling strength, there is a direct energy-level crossing in the ground state of the system between states of different Jz{J}_{z}Jz​, the third component of the total angular momentum, that changes its structure. This is reflected in a discontinuity in the different energy terms and it is signaled in the spatial density of the system.
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