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

Energy-level statistics of interacting trapped bosons

Barnali Chakrabarti, Anindya Biswas, V. K. B. Kota, Kamalika Roy, Sudip Kumar Haldar

Published Jul 23, 2012• 10 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreva.86.013637

Related papers

Spectral fluctuation and 1/fα noise in the energy level statistics of interacting trapped bosons

Kamalika Roy, Barnali Chakrabarti, Anindya Biswas +2 more

It has been recently shown numerically that the transition from integrability to chaos in quantum systems and the corresponding spectral fluctuations are characterized by 1/f^{α} noise with 1≤α≤2. The system of interacting trapped bosons is inhomogeneous and complex. The presence of an external harmonic trap makes it more interesting as, in the atomic trap, the bosons occupy partly degenerate single-particle states. Earlier theoretical and experimental results show that at zero temperature the low-lying levels are of a collective nature and high-lying excitations are of a single-particle nature. We observe that for few bosons, the P(s) distribution shows the Shnirelman peak, which exhibits a large number of quasidegenerate states. For a large number of bosons the low-lying levels are strongly affected by the interatomic interaction, and the corresponding level fluctuation shows a transition to a Wigner distribution with an increase in particle number. It does not follow Gaussian orthogonal ensemble random matrix predictions. For high-lying levels we observe the uncorrelated Poisson distribution. Thus it may be a very realistic system to prove that 1/f^{α} noise is ubiquitous in nature.

arXiv:1112.3718•Physical Review E•Jun 20, 2012•8 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

Spectral fluctuation and correlation structure of δn statistics in the spectra of interacting trapped bosons

Kamalika Roy, Barnali Chakrabarti, V. K. B. Kota

It is a well-known fact that the statistical behaviors of level fluctuation and level correlation in the energy-level spectra are the most efficient tool to characterize quantum chaos in nonintegrable quantum systems. The system of interacting trapped bosons is a complex system where the low-lying energy levels are highly influenced by the level repulsion. In this case, interatomic interaction is a dominating fact with strong level correlation between distant levels. Here we numerically calculate the correlation function, number variance, and Dyson-Mehta least-square deviation for the low-lying levels for a few thousand interacting trapped bosons, and our data show good analogy with the Gaussian orthogonal ensemble (GOE) results with a signature of chaos. In the next part of our study, the energy spectrum of these low-lying levels is considered as a discrete signal and the fluctuation of the excitation energy is considered as discrete time series. Then we calculate numerically the height-height correlation function for different order of momentum. In our study logarithmic correlation structure is found instead of multiscaling structure, and we observe that spectral statistics are compatible with those of GOE.

Physical Review E•Jun 3, 2013•2 citations

Statistical properties of spectral fluctuations of N interacting bosons in a harmonic trap

Kamalika Roy, Barnali Chakrabarti, V. K. B. Kota

Spectral fluctuations of a system of N weakly interacting bosons in an isotropic harmonic trap are studied, with the focus on the deviations from Poisson spectral statistics, typical of a quantum integrable systems. We have utilized the ideas formulated by Makino et al. [Phys. Rev. E 67, 066205 (2003)1063-651X10.1103/PhysRevE.67.066205] who have extended the approach of Berry and Robnik [J. Phys. A 17, 2413 (1984)JPHAC50305-447010.1088/0305-4470/17/12/013]. Earlier investigations of the Berry conjecture [Proc. R. Soc. London, Ser. A 356, 375 (1977)1364-502110.1098/rspa.1977.0140] of Poisson spectral statistics mainly considered quantum systems whose classical counterparts are integrable. However, the system of N weakly interacting bosons in the external trap has no classical counterpart. Also, it is a realistic and experimentally achievable system with close relation to Bose-Einstein condensation. Thus, a stringent analysis of the applicability of the Berry conjecture to this kind of systems is indeed required. We observe that for small boson number, the system is close to integrability and the nearest-neighbor level spacing distribution and the level number variance exhibit deviations from Poisson statistics similar to those of rational rectangular billiards.

Physical Review E•Nov 18, 2014•1 citation

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

Energy fluctuation of a finite number of interacting bosons: A correlated many-body approach

Satadal Bhattacharyya, M. L. Lekala, Barnali Chakrabarti +1 more

We calculate the energy fluctuation of a truly finite number of interacting bosons and study the role of interaction. Although the ideal Bose gas in thermodynamic limit is an exactly solvable problem and analytic expression of various fluctuation measures exists, the experimental Bose-Einstein condensation (BEC) is a nontrivial many-body problem. We employ a two-body correlated basis function and utilize the realistic van der Waals interaction. We calculate the energy fluctuation $(△{E}^{2})$ of the interacting trapped bosons and plot $\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}$ as a function of $\frac{T}{{T}_{c}}$. In the classical limit $△{E}^{2}$ is related to the specific heat per particle ${c}_{v}$ through the relation $△{E}^{2}={k}_{B}{T}^{2}{c}_{v}$. We have obtained a distinct hump in $\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}$ around the condensation point for three-dimesional harmonically trapped Bose gas when the particle number $N\ensuremath{\simeq}5000$ and above which corresponds to the second-order phase transition. However for finite-size interacting bosons ($N\ensuremath{\simeq}$ a few hundred) the hump is not sharp, and the maximum in $\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}$ can be interpreted as a smooth increase in the scaled fluctuation below ${T}_{c}$ and then a decrease above ${T}_{c}$. To illustrate the justification we also calculate ${c}_{v}$, which exhibits the same feature, which leads to the conjecture that for finite-sized interacting bosons phase transition is ruled out.

Physical Review A•Mar 14, 2016•4 citations
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.
It has been recently shown numerically that the transition from integrability to chaos in quantum systems and the corresponding spectral fluctuations are characterized by 1/f^{α} noise with 1≤α≤2. The system of interacting trapped bosons is inhomogeneous and complex. The presence of an external harmonic trap makes it more interesting as, in the atomic trap, the bosons occupy partly degenerate single-particle states. Earlier theoretical and experimental results show that at zero temperature the low-lying levels are of a collective nature and high-lying excitations are of a single-particle nature. We observe that for few bosons, the P(s) distribution shows the Shnirelman peak, which exhibits a large number of quasidegenerate states. For a large number of bosons the low-lying levels are strongly affected by the interatomic interaction, and the corresponding level fluctuation shows a transition to a Wigner distribution with an increase in particle number. It does not follow Gaussian orthogonal ensemble random matrix predictions. For high-lying levels we observe the uncorrelated Poisson distribution. Thus it may be a very realistic system to prove that 1/f^{α} noise is ubiquitous in nature.
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.
It is a well-known fact that the statistical behaviors of level fluctuation and level correlation in the energy-level spectra are the most efficient tool to characterize quantum chaos in nonintegrable quantum systems. The system of interacting trapped bosons is a complex system where the low-lying energy levels are highly influenced by the level repulsion. In this case, interatomic interaction is a dominating fact with strong level correlation between distant levels. Here we numerically calculate the correlation function, number variance, and Dyson-Mehta least-square deviation for the low-lying levels for a few thousand interacting trapped bosons, and our data show good analogy with the Gaussian orthogonal ensemble (GOE) results with a signature of chaos. In the next part of our study, the energy spectrum of these low-lying levels is considered as a discrete signal and the fluctuation of the excitation energy is considered as discrete time series. Then we calculate numerically the height-height correlation function for different order of momentum. In our study logarithmic correlation structure is found instead of multiscaling structure, and we observe that spectral statistics are compatible with those of GOE.
Spectral fluctuations of a system of N weakly interacting bosons in an isotropic harmonic trap are studied, with the focus on the deviations from Poisson spectral statistics, typical of a quantum integrable systems. We have utilized the ideas formulated by Makino et al. [Phys. Rev. E 67, 066205 (2003)1063-651X10.1103/PhysRevE.67.066205] who have extended the approach of Berry and Robnik [J. Phys. A 17, 2413 (1984)JPHAC50305-447010.1088/0305-4470/17/12/013]. Earlier investigations of the Berry conjecture [Proc. R. Soc. London, Ser. A 356, 375 (1977)1364-502110.1098/rspa.1977.0140] of Poisson spectral statistics mainly considered quantum systems whose classical counterparts are integrable. However, the system of N weakly interacting bosons in the external trap has no classical counterpart. Also, it is a realistic and experimentally achievable system with close relation to Bose-Einstein condensation. Thus, a stringent analysis of the applicability of the Berry conjecture to this kind of systems is indeed required. We observe that for small boson number, the system is close to integrability and the nearest-neighbor level spacing distribution and the level number variance exhibit deviations from Poisson statistics similar to those of rational rectangular billiards.
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 calculate the energy fluctuation of a truly finite number of interacting bosons and study the role of interaction. Although the ideal Bose gas in thermodynamic limit is an exactly solvable problem and analytic expression of various fluctuation measures exists, the experimental Bose-Einstein condensation (BEC) is a nontrivial many-body problem. We employ a two-body correlated basis function and utilize the realistic van der Waals interaction. We calculate the energy fluctuation (△E2)(△{E}^{2})(△E2) of the interacting trapped bosons and plot △E2kB2T2\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}kB2​T2△E2​ as a function of TTc\frac{T}{{T}_{c}}Tc​T​. In the classical limit △E2△{E}^{2}△E2 is related to the specific heat per particle cv{c}_{v}cv​ through the relation △E2=kBT2cv△{E}^{2}={k}_{B}{T}^{2}{c}_{v}△E2=kB​T2cv​. We have obtained a distinct hump in △E2kB2T2\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}kB2​T2△E2​ around the condensation point for three-dimesional harmonically trapped Bose gas when the particle number $N\ensuremath{\simeq}5000$ and above which corresponds to the second-order phase transition. However for finite-size interacting bosons ($N\ensuremath{\simeq}$ a few hundred) the hump is not sharp, and the maximum in △E2kB2T2\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}kB2​T2△E2​ can be interpreted as a smooth increase in the scaled fluctuation below Tc{T}_{c}Tc​ and then a decrease above Tc{T}_{c}Tc​. To illustrate the justification we also calculate cv{c}_{v}cv​, which exhibits the same feature, which leads to the conjecture that for finite-sized interacting bosons phase transition is ruled out.