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:0805.3388 · Physical Review Letters

Threshold for Chaos and Thermalization in the One-Dimensional Mean-Field Bose-Hubbard Model

Amy Cassidy, Douglas Mason, Vanja Dunjko, Maxim Olshanii

Published Jan 16, 2009• 49 citations
arXivHTMLPDF

Abstract

We study the threshold for chaos and its relation to thermalization in the 1D mean-field Bose-Hubbard model, which, in particular, describes atoms in optical lattices. We identify the threshold for chaos, which is finite in the thermodynamic limit, and show that it is indeed a precursor of thermalization. Far above the threshold, the state of the system after relaxation is governed by the usual laws of statistical mechanics.

DOI
10.1103/physrevlett.102.025302

Related papers

Chaos and Thermalization in the one-dimensional Bose-Hubbard model in the classical-field approximation

Amy Cassidy

In this thesis, we present a comprehensive study of chaos and thermalization of the one-dimensional Bose-Hubbard Model (BHM) within the classical field approximation. Two quantitative measures are compared: the ensemble-averaged Finite-time Maximal Lyapunov exponent, a measures of chaos and the normalized spectral entropy, a measure of the distance between the numerical time-averaged momentum distribution and the one predicted by thermodynamics. A threshold for chaos is found, which depends on two parameters, the nonlinearity and the total energy-per-particle. Below the threshold, the dynamics are regular, while far above the threshold, complete thermalization is observed, as measured by the normalized spectral entropy. We study individual resonances in the Bose-Hubbard model to determine the criterion for chaos. The criterion based on Chirikov's method of overlapping resonances diverges in the thermodynamic limit, in contrast to the criterion parameters inferred from numerical calculations, signifying the failure of the standard Chirikov's approach. The Ablowitz-Ladik lattice is one of several integrable models that are close to the BHM. We outline the method of Inverse Scattering Transform and generate the integrals of motion of the Ablowitz-Ladik lattice. Furthermore, we discuss the possible role of these quantities in the relaxation dynamics of the BHM.

arXiv:1003.5206•Mar 23, 2010•1 citation

Threshold coupling strength for equilibration between small systems

Ralf Bürkle, J. R. Anglin

In this paper we study the thermal equilibration of small bipartite Bose-Hubbard systems, both quantum mechanically and in mean-field approximation. In particular we consider small systems composed of a single-mode ``thermometer'' coupled to a three-mode ``bath,'' with no additional environment acting on the four-mode system, and test the hypothesis that the thermometer will thermalize if and only if the bath is chaotic. We find that chaos in the bath alone is neither necessary nor sufficient for equilibration in these isolated four-mode systems. The two subsystems can thermalize if the combined system is chaotic even when neither subsystem is chaotic in isolation, and under full quantum dynamics there is a minimum coupling strength between the thermometer and the bath below which the system does not thermalize even if the bath itself is chaotic. We show that the quantum coupling threshold scales like $1/N$ (where $N$ is the total particle number), so that the classical results are obtained in the limit $N\ensuremath{\rightarrow}\ensuremath{\infty}$.

arXiv:1903.12083•Physical Review A•Jun 19, 2019•4 citations

Superdiffusion, normal diffusion, and chaos in semiclassical Bose-Hubbard chains

Dragan Marković, Mihailo Čubrović

We study the evolution of two-point correlation functions of the one-dimensional Bose-Hubbard model in the semiclassical regime in the framework of truncated Wigner approximation with quantum jumps as first-order corrections. At early times, the correlation functions show strong superdiffusion with universal integer exponents determined solely by the initial conditions and completely insensitive to system parameters and chaos. Only after a long time does this regime crosse over to the normal diffusion regime which is most robust when nonintegrability is strong. For strong nonintegrability, the system ends up in a homogeneous state, while for weak nonintegrability, the oscillations and inhomogeneities persist, despite the fact that chaos is nearly always strong and only weakly depends on the nonintegrability parameter. We conclude that the superdiffusive regime is neither prethermalized nor a precursor to thermalization but an early-time phenomenon related to a special scaling symmetry of the Bose-Hubbard Hamiltonian.

arXiv:2502.19584•Physical review. E•Jul 15, 2025•1 citation

Stochastic mean-field theory for the disordered Bose-Hubbard model

Ulf Bissbort, Walter Hofstetter

We investigate the effect of diagonal disorder on bosons in an optical lattice described by an Anderson-Hubbard model at zero temperature. It is known that within Gutzwiller mean-field theory spatially resolved calculations suffer particularly from finite system sizes in the disordered case, while arithmetic averaging of the order parameter cannot describe the Bose glass phase for finite hopping J >0. Here we present and apply a new stochastic mean-field theory which captures localization due to disorder, includes non-trivial dimensional effects beyond the mean-field scaling level and is applicable in the thermodynamic limit. In contrast to fermionic systems, we find the existence of a critical hopping strength, above which the system remains superfluid for arbitrarily strong disorder.

arXiv:0804.0007•Europhysics Letters (EPL)•Jun 1, 2009•52 citations

Multiparticle quantum chaos in tilted optical lattices

Andrey R. Kolovsky, Andreas Buchleitner

It is shown that, in the parameter regime of state of the art experiments on Bose-Einstein condensates loaded into optical lattices, the energy spectrum of the 1D Bose-Hubbard model amended by a static field exhibits unambiguous signatures of quantum chaos. In the dynamics, this leads to the irreversible decay of Bloch oscillations.

arXiv:cond-mat/0312041•Journal of Modern Optics•Apr 1, 2004•2 citations

Bose–Hubbard Hamiltonian: Quantum chaos approach

Andrey R. Kolovsky

We discuss applications of the theory of quantum chaos to one of the paradigm models of many-body quantum physics — the Bose–Hubbard (BH) model, which describes, in particular, interacting ultracold Bose atoms in an optical lattice. After preliminary, pure quantum analysis of the system we introduce the classical counterpart of the BH model and the governing semiclassical equations of motion. We analyze these equations for the problem of Bloch oscillations (BOs) of cold atoms where a number of experimental results are available. The paper is written for nonexperts and can be viewed as an introduction to the field.

arXiv:1507.03413•International Journal of Modern Physics B•Mar 28, 2016•40 citations
We study the threshold for chaos and its relation to thermalization in the 1D mean-field Bose-Hubbard model, which, in particular, describes atoms in optical lattices. We identify the threshold for chaos, which is finite in the thermodynamic limit, and show that it is indeed a precursor of thermalization. Far above the threshold, the state of the system after relaxation is governed by the usual laws of statistical mechanics.
In this thesis, we present a comprehensive study of chaos and thermalization of the one-dimensional Bose-Hubbard Model (BHM) within the classical field approximation. Two quantitative measures are compared: the ensemble-averaged Finite-time Maximal Lyapunov exponent, a measures of chaos and the normalized spectral entropy, a measure of the distance between the numerical time-averaged momentum distribution and the one predicted by thermodynamics. A threshold for chaos is found, which depends on two parameters, the nonlinearity and the total energy-per-particle. Below the threshold, the dynamics are regular, while far above the threshold, complete thermalization is observed, as measured by the normalized spectral entropy. We study individual resonances in the Bose-Hubbard model to determine the criterion for chaos. The criterion based on Chirikov's method of overlapping resonances diverges in the thermodynamic limit, in contrast to the criterion parameters inferred from numerical calculations, signifying the failure of the standard Chirikov's approach. The Ablowitz-Ladik lattice is one of several integrable models that are close to the BHM. We outline the method of Inverse Scattering Transform and generate the integrals of motion of the Ablowitz-Ladik lattice. Furthermore, we discuss the possible role of these quantities in the relaxation dynamics of the BHM.
In this paper we study the thermal equilibration of small bipartite Bose-Hubbard systems, both quantum mechanically and in mean-field approximation. In particular we consider small systems composed of a single-mode ``thermometer'' coupled to a three-mode ``bath,'' with no additional environment acting on the four-mode system, and test the hypothesis that the thermometer will thermalize if and only if the bath is chaotic. We find that chaos in the bath alone is neither necessary nor sufficient for equilibration in these isolated four-mode systems. The two subsystems can thermalize if the combined system is chaotic even when neither subsystem is chaotic in isolation, and under full quantum dynamics there is a minimum coupling strength between the thermometer and the bath below which the system does not thermalize even if the bath itself is chaotic. We show that the quantum coupling threshold scales like 1/N1/N1/N (where NNN is the total particle number), so that the classical results are obtained in the limit $N\ensuremath{\rightarrow}\ensuremath{\infty}$.
We study the evolution of two-point correlation functions of the one-dimensional Bose-Hubbard model in the semiclassical regime in the framework of truncated Wigner approximation with quantum jumps as first-order corrections. At early times, the correlation functions show strong superdiffusion with universal integer exponents determined solely by the initial conditions and completely insensitive to system parameters and chaos. Only after a long time does this regime crosse over to the normal diffusion regime which is most robust when nonintegrability is strong. For strong nonintegrability, the system ends up in a homogeneous state, while for weak nonintegrability, the oscillations and inhomogeneities persist, despite the fact that chaos is nearly always strong and only weakly depends on the nonintegrability parameter. We conclude that the superdiffusive regime is neither prethermalized nor a precursor to thermalization but an early-time phenomenon related to a special scaling symmetry of the Bose-Hubbard Hamiltonian.
We investigate the effect of diagonal disorder on bosons in an optical lattice described by an Anderson-Hubbard model at zero temperature. It is known that within Gutzwiller mean-field theory spatially resolved calculations suffer particularly from finite system sizes in the disordered case, while arithmetic averaging of the order parameter cannot describe the Bose glass phase for finite hopping J >0. Here we present and apply a new stochastic mean-field theory which captures localization due to disorder, includes non-trivial dimensional effects beyond the mean-field scaling level and is applicable in the thermodynamic limit. In contrast to fermionic systems, we find the existence of a critical hopping strength, above which the system remains superfluid for arbitrarily strong disorder.
We discuss applications of the theory of quantum chaos to one of the paradigm models of many-body quantum physics — the Bose–Hubbard (BH) model, which describes, in particular, interacting ultracold Bose atoms in an optical lattice. After preliminary, pure quantum analysis of the system we introduce the classical counterpart of the BH model and the governing semiclassical equations of motion. We analyze these equations for the problem of Bloch oscillations (BOs) of cold atoms where a number of experimental results are available. The paper is written for nonexperts and can be viewed as an introduction to the field.