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:1004.2203 · Journal of Statistical Mechanics Theory and Experiment

Statistical properties of the spectrum of the extended Bose–Hubbard model

Corinna Kollath, Guillaume Roux, Giulio Biroli, Andreas Martin Läuchli

Published Aug 1, 2010• 98 citations
arXivHTMLPDF

Abstract

Motivated by the role that spectral properties play for the dynamical evolution of a quantum many-body system, we investigate the level spacing statistics of the extended Bose–Hubbard model. In particular, we focus on the distribution of the ratio of adjacent level spacings, useful at large interactions, to distinguish between chaotic and non-chaotic regimes. After revisiting the bare Bose–Hubbard model, we study the effect of two different perturbations: next-nearest-neighbor hopping and nearest-neighbor interaction. The system size dependence is investigated together with the effect of the proximity to integrable points or lines. Lastly, we discuss the consequences of a cutoff in the number of onsite bosons on the level statistics.

DOI
10.1088/1742-5468/2010/08/p08011

Related papers

Spectral analysis of two-dimensional Bose-Hubbard models

David Fischer, Darius Hoffmann, Sandro Wimberger

One-dimensional Bose-Hubbard models are well known to obey a transition from regular to quantum-chaotic spectral statistics. We are extending this concept to relatively simple two-dimensional many-body models. Also in two dimensions a transition from regular to chaotic spectral statistics is found and discussed. In particular, we analyze the dependence of the spectral properties on the bond number of the two-dimensional lattices and the applied boundary conditions. For maximal connectivity, the systems behave most regularly in agreement with the applicability of mean-field approaches in the limit of many nearest-neighbor couplings at each site.

arXiv:1604.05664•Physical Review A•Apr 22, 2016•16 citations

Nonergodic behavior of the clean Bose-Hubbard chain

Angelo Russomanno, Michele Fava, Rosario Fazio

We study ergodicity breaking in the clean Bose-Hubbard chain for small hopping strength. We see the existence of a nonergodic regime by means of indicators as the half-chain entanglement entropy of the eigenstates, the average level spacing ratio, the properties of the eigenstate-expectation distribution of the correlation and the scaling of the inverse participation ratio averages. We find that this ergodicity breaking is different from many-body localization because the average half-chain entanglement entropy of the eigenstates obeys volume law. This ergodicity breaking appears unrelated to the spectrum being organized in quasidegenerate multiplets at small hopping and finite system sizes, so in principle, it can survive also for larger system sizes. We find that some imbalance oscillations in time which could mark the existence of glassy behavior in space are well described by the dynamics of a single symmetry-breaking doublet and quantitatively captured by a perturbative effective XXZ model. We show that the amplitude of these oscillations vanishes in the large-size limit. Our findings are numerically obtained for systems for $L<12$. Extrapolations of our scalings to larger system sizes should be taken with care, as discussed in the paper.

arXiv:2005.09892•Physical review. B./Physical review. B•Oct 9, 2020•24 citations

Statistical and dynamical aspects of quantum chaos in a kicked Bose-Hubbard dimer

Chenguang Liang, Yu Zhang, Shu Chen

Systems of interacting bosons in double-well potentials, modeled by two-site Bose-Hubbard models, are of significant theoretical and experimental interest and attracted intensive studies in contexts ranging from many-body physics and quantum dynamics to the onset of quantum chaos. In this work we systematically study a kicked two-site Bose-Hubbard model (Bose-Hubbard dimer) with the on-site potential difference being periodically modulated. Our model can be equivalently represented as a kicked Lipkin-Meshkov-Glick model and thus displays different dynamical behaviors from the kicked-top model. By analyzing spectral statistics of Floquet operator, we unveil that the system undergoes a transition from regularity to chaos with increasing the interaction strength. Then based on semiclassical approximation and the analysis of R\'enyi entropy of coherent states in the basis of Floquet operator eigenstates, we reveal the local chaotic features of our model, which indicate the existence of integrable islands even in the deep chaotic regime. The semiclassical analysis also suggests that the system in chaotic regime may display different dynamical behavior depending on the choice of initial states. Finally, we demonstrate that dynamical signatures of chaos can be manifested by studying dynamical evolution of local operators and out-of-time-order correlation function as well as the entanglement entropy. Our numerical results exhibit the richness of dynamics of the kicked Bose-Hubbard dimer in both regular and chaotic regimes as the initial states are chosen as coherent spin states located in different locations of phase space.

arXiv:2312.08159•Physical Review A•Mar 22, 2024•7 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

Spectrum statistics in the integrable Lieb-Liniger model

Samy Mailoud, F. Borgonovi, Felix M. Izrailev

We address the old and widely debated question of the spectrum statistics of integrable quantum systems, through the analysis of the paradigmatic Lieb-Liniger model. This quantum many-body model of one-dimensional interacting bosons allows for the rigorous determination of energy spectra via the Bethe ansatz approach and our interest is to reveal the characteristic properties of energy levels in dependence of the model parameters. Using both analytical and numerical studies we show that the properties of spectra strongly depend on whether the analysis is done for a full energy spectrum or for a single subset with fixed total momentum. We show that the Poisson distribution of spacing between nearest-neighbor energies can occur only for a set of energy levels with fixed total momentum, for neither too large nor too weak interaction strength, and for sufficiently high energy. By studying long-range correlations between energy levels, we found strong deviations from the predictions based on the assumption of pseudorandom character of the distribution of energy levels.

arXiv:2105.02967•Physical review. E•Sep 15, 2021•6 citations

Krylov complexity of a tilted extended Bose-Hubbard chain with Rydberg-dressed interactions

Yifan Chen, Tianyi Yan, Lu Qin +1 more

We investigate Krylov state complexity in a tilted extended Bose-Hubbard chain in which both the on-site interaction $U$ and nearest-neighbor interaction $V$ are present. The tilted extended Bose-Hubbard model can be realized with Rydberg-dressed interactions in optical lattices. Using exact diagonalization and Lanczos recursion from three physically motivated Fock states, we compute the time-dependent complexity, its long-time saturation value, and the fluctuations of the Lanczos coefficients for an open chain at unit filling. Crucially, we identify a pronounced quasi-chaotic diagonal regime near $U\simeq V$ in which the Krylov complexity saturation exhibits a sharp ridge while the spectral level-spacing statistics remain partially chaotic. This dissociation arises from the competition between the on-site and nearest-neighbor interaction, which generates a dense manifold of near-resonant many-body configurations that is efficiently coupled by the state-dependent dynamics but retains residual spectral structure. Our results show that Krylov state complexity provides complementary, state-sensitive information about many-body chaos that is not captured by the spectral diagnostics alone, establishing Krylov complexity as a useful diagnostic for studying thermalization and information scrambling in experimentally accessible quantum simulators.

arXiv:2610.00712•Sep 30, 2026•0 citations
Motivated by the role that spectral properties play for the dynamical evolution of a quantum many-body system, we investigate the level spacing statistics of the extended Bose–Hubbard model. In particular, we focus on the distribution of the ratio of adjacent level spacings, useful at large interactions, to distinguish between chaotic and non-chaotic regimes. After revisiting the bare Bose–Hubbard model, we study the effect of two different perturbations: next-nearest-neighbor hopping and nearest-neighbor interaction. The system size dependence is investigated together with the effect of the proximity to integrable points or lines. Lastly, we discuss the consequences of a cutoff in the number of onsite bosons on the level statistics.
One-dimensional Bose-Hubbard models are well known to obey a transition from regular to quantum-chaotic spectral statistics. We are extending this concept to relatively simple two-dimensional many-body models. Also in two dimensions a transition from regular to chaotic spectral statistics is found and discussed. In particular, we analyze the dependence of the spectral properties on the bond number of the two-dimensional lattices and the applied boundary conditions. For maximal connectivity, the systems behave most regularly in agreement with the applicability of mean-field approaches in the limit of many nearest-neighbor couplings at each site.
We study ergodicity breaking in the clean Bose-Hubbard chain for small hopping strength. We see the existence of a nonergodic regime by means of indicators as the half-chain entanglement entropy of the eigenstates, the average level spacing ratio, the properties of the eigenstate-expectation distribution of the correlation and the scaling of the inverse participation ratio averages. We find that this ergodicity breaking is different from many-body localization because the average half-chain entanglement entropy of the eigenstates obeys volume law. This ergodicity breaking appears unrelated to the spectrum being organized in quasidegenerate multiplets at small hopping and finite system sizes, so in principle, it can survive also for larger system sizes. We find that some imbalance oscillations in time which could mark the existence of glassy behavior in space are well described by the dynamics of a single symmetry-breaking doublet and quantitatively captured by a perturbative effective XXZ model. We show that the amplitude of these oscillations vanishes in the large-size limit. Our findings are numerically obtained for systems for L<12L<12L<12. Extrapolations of our scalings to larger system sizes should be taken with care, as discussed in the paper.
Systems of interacting bosons in double-well potentials, modeled by two-site Bose-Hubbard models, are of significant theoretical and experimental interest and attracted intensive studies in contexts ranging from many-body physics and quantum dynamics to the onset of quantum chaos. In this work we systematically study a kicked two-site Bose-Hubbard model (Bose-Hubbard dimer) with the on-site potential difference being periodically modulated. Our model can be equivalently represented as a kicked Lipkin-Meshkov-Glick model and thus displays different dynamical behaviors from the kicked-top model. By analyzing spectral statistics of Floquet operator, we unveil that the system undergoes a transition from regularity to chaos with increasing the interaction strength. Then based on semiclassical approximation and the analysis of R\'enyi entropy of coherent states in the basis of Floquet operator eigenstates, we reveal the local chaotic features of our model, which indicate the existence of integrable islands even in the deep chaotic regime. The semiclassical analysis also suggests that the system in chaotic regime may display different dynamical behavior depending on the choice of initial states. Finally, we demonstrate that dynamical signatures of chaos can be manifested by studying dynamical evolution of local operators and out-of-time-order correlation function as well as the entanglement entropy. Our numerical results exhibit the richness of dynamics of the kicked Bose-Hubbard dimer in both regular and chaotic regimes as the initial states are chosen as coherent spin states located in different locations of phase space.
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 address the old and widely debated question of the spectrum statistics of integrable quantum systems, through the analysis of the paradigmatic Lieb-Liniger model. This quantum many-body model of one-dimensional interacting bosons allows for the rigorous determination of energy spectra via the Bethe ansatz approach and our interest is to reveal the characteristic properties of energy levels in dependence of the model parameters. Using both analytical and numerical studies we show that the properties of spectra strongly depend on whether the analysis is done for a full energy spectrum or for a single subset with fixed total momentum. We show that the Poisson distribution of spacing between nearest-neighbor energies can occur only for a set of energy levels with fixed total momentum, for neither too large nor too weak interaction strength, and for sufficiently high energy. By studying long-range correlations between energy levels, we found strong deviations from the predictions based on the assumption of pseudorandom character of the distribution of energy levels.
We investigate Krylov state complexity in a tilted extended Bose-Hubbard chain in which both the on-site interaction UUU and nearest-neighbor interaction VVV are present. The tilted extended Bose-Hubbard model can be realized with Rydberg-dressed interactions in optical lattices. Using exact diagonalization and Lanczos recursion from three physically motivated Fock states, we compute the time-dependent complexity, its long-time saturation value, and the fluctuations of the Lanczos coefficients for an open chain at unit filling. Crucially, we identify a pronounced quasi-chaotic diagonal regime near U≃VU\simeq VU≃V in which the Krylov complexity saturation exhibits a sharp ridge while the spectral level-spacing statistics remain partially chaotic. This dissociation arises from the competition between the on-site and nearest-neighbor interaction, which generates a dense manifold of near-resonant many-body configurations that is efficiently coupled by the state-dependent dynamics but retains residual spectral structure. Our results show that Krylov state complexity provides complementary, state-sensitive information about many-body chaos that is not captured by the spectral diagnostics alone, establishing Krylov complexity as a useful diagnostic for studying thermalization and information scrambling in experimentally accessible quantum simulators.