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:1003.5206

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

Amy Cassidy

Published Mar 23, 2010• 1 citation
arXivHTMLPDF

Abstract

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.

DOI
10.25549/usctheses-m2217

Related papers

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

Amy Cassidy, Douglas Mason, Vanja Dunjko +1 more

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.

arXiv:0805.3388•Physical Review Letters•Jan 16, 2009•49 citations

Relaxation and thermalization in the one-dimensional Bose-Hubbard model: A case study for the interaction quantum quench from the atomic limit

Stefan Sorg, Lev Vidmar, Lode Pollet +1 more

Motivated by recent experiments, we study the relaxation dynamics and thermalization in the one-dimensional Bose-Hubbard model induced by a global interaction quench. Specifically, we start from an initial state that has exactly one boson per site and is the ground state of a system with infinitely strong repulsive interactions at unit filling. Using exact diagonalization and the density-matrix renormalization-group method, we compute the time dependence of such observables as the multiple occupancy and the momentum distribution function. Typically, the relaxation to stationary values occurs over just a few tunneling times. The stationary values are identical to the so-called diagonal ensemble on the system sizes accessible to our numerical methods and we further observe that the microcanonical ensemble describes the time averages of many observables reasonably well for small and intermediate interaction strength. The expectation values of observables in the canonical ensemble agree quantitatively with the time averages obtained from the quench at small interaction strengths, and qualitatively provide a good description even in parameter regimes where the microcanonical ensemble is not applicable due to finite-size effects. We discuss our numerical results in the framework of the eigenstate thermalization hypothesis. Moreover, we also observe that the diagonal and the canonical ensembles are practically identical for our initial conditions already on the level of their respective energy distributions for small interaction strengths. Finally, we discuss implications of our results for the interpretation of a recent sudden expansion experiment [Phys. Rev. Lett. 110, 205301 (2013)], in which the same interaction quench was realized.

arXiv:1405.5404•Physical Review A•Sep 8, 2014•129 citations

Chaos in the three-site Bose-Hubbard model: Classical versus quantum

Goran Nakerst, Masudul Haque

We consider a quantum many-body system-the Bose-Hubbard system on three sites-which has a classical limit, and which is neither strongly chaotic nor integrable but rather shows a mixture of the two types of behavior. We compare quantum measures of chaos (eigenvalue statistics and eigenvector structure) in the quantum system, with classical measures of chaos (Lyapunov exponents) in the corresponding classical system. As a function of energy and interaction strength, we demonstrate a strong overall correspondence between the two cases. In contrast to both strongly chaotic and integrable systems, the largest Lyapunov exponent is shown to be a multivalued function of energy.

arXiv:2203.09953•Physical review. E•Feb 21, 2023•25 citations

Out-of-equilibrium dynamics of the Bose-Hubbard model in the strong coupling regime

Matthew Fitzpatrick

Experimental advances have made ultracold atoms in optical lattices a favourable setting to study out-of-equilibrium phenomena and attracted considerable attention in recent years.These systems are highly versatile in that experimental parameters can be tuned over a wide range of values in real time.This facilitates the study of quantum quenches, in which parameters in the corresponding Hamiltonian are varied in time faster than the system can respond adiabatically.Such protocols open the door to a rich range of many-body physics and have been studied intensely both theoretically and experimentally.The Bose-Hubbard model (BHM) has been shown to describe interacting ultracold bosons in an optical lattice, allowing the opportunity for experiments to probe the out-of-equilibrium dynamics of the model.The BHM is a particularly convenient context for studying quantum quenches as it displays a quantum phase transition between superfluid and Mott-insulator phases.In this thesis, we develop a strong-coupling approach that allows the study of correlations in space and time in both the superfluid and Mott-insulating phases of the BHM.Specifically, we obtain a two-particle irreducible (2PI) effective action within the contourtime formalism that allows for a description of both equilibrium and out-of-equilibrium phenomena.We derive equations of motion for both the superfluid order parameter and the single-particle many-body Green's functions.First, we assess the accuracy of our formalism by studying the equilibrium solution for the homogeneous BHM and comparing our results to existing strong-coupling methods as well as to exact methods where possible.We then consider homogeneous systems that are initially thermalized in the Mott phase, and which are then subjected to quenches.We solve numerically the equations of motion for this scenario and calculate the single-particle density matrix.We demonstrate a Lieb-Robinson-like maximal propagation velocity for the spreading of single-particle correlations in one, two, and three dimensions.We study the dependence of the maximal propagation velocity on the quench protocol, chemical potential, temperature, and dimensionality.We compare our results to exact methods, existing strong-coupling approaches, and experiments where possible.Lastly, we extend our strong-coupling approach to the disordered BHM and derive equations of motion for the disorder-averaged single-particle Green's function.We discuss how these equations of motion can be used to study the phase stability of many-body localization in the disordered BHM for dimensions higher than one.

•Jul 5, 2019•1 citation

Dynamical thermalization in Bose-Hubbard systems

Peter Schlagheck, Dima L. Shepelyansky

We numerically study a Bose-Hubbard ring of finite size with disorder containing a finite number of bosons that are subject to an on-site two-body interaction. Our results show that moderate interactions induce dynamical thermalization in this isolated system. In this regime the individual many-body eigenstates are well described by the standard thermal Bose-Einstein distribution for well-defined values of the temperature and the chemical potential, which depend on the eigenstate under consideration. We show that the dynamical thermalization conjecture works well at both positive and negative temperatures. The relations to quantum chaos, quantum ergodicity, and the Åberg criterion are also discussed.

arXiv:1510.01864•Physical review. E•Jan 19, 2016•11 citations

Dynamical witnesses and universal behavior across chaos and non-ergodicity in the tilted Bose-Hubbard model

Carlos Diaz-Mejia, Sergio Lerma-Hernandez, Jorge G. Hirsch

Quantum chaos in isolated quantum systems is intimately linked to thermalization and the rapid relaxation of observables. Although the spectral properties of the chaotic phase in the tilted Bose-Hubbard model have been well characterized, the corresponding dynamical signatures across the transition to regularity remain less explored . In this work, we investigate this transition by analyzing the time evolution of the survival probability, the single-site entanglement entropy, and the half-chain imbalance. Our results reveal a clear hierarchy in the sensitivity of these observables: the relaxation value of the entanglement entropy varies smoothly as a function of the Hamiltonian parameters across the chaos-regular transition, while the imbalance exhibits a more pronounced distinction. Most notably, the survival probability emerges as the most robust indicator of the transition between chaos and regularity. When appropriately scaled, all three observables converge onto a common behavior as a function of the Hamiltonian parameters for different numbers of sites and bosons,enabling a universal characterization of the transition between chaotic and regular dynamics.

quant-pharXiv:2602.00369•Feb 3, 2026•0 citations
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.
Motivated by recent experiments, we study the relaxation dynamics and thermalization in the one-dimensional Bose-Hubbard model induced by a global interaction quench. Specifically, we start from an initial state that has exactly one boson per site and is the ground state of a system with infinitely strong repulsive interactions at unit filling. Using exact diagonalization and the density-matrix renormalization-group method, we compute the time dependence of such observables as the multiple occupancy and the momentum distribution function. Typically, the relaxation to stationary values occurs over just a few tunneling times. The stationary values are identical to the so-called diagonal ensemble on the system sizes accessible to our numerical methods and we further observe that the microcanonical ensemble describes the time averages of many observables reasonably well for small and intermediate interaction strength. The expectation values of observables in the canonical ensemble agree quantitatively with the time averages obtained from the quench at small interaction strengths, and qualitatively provide a good description even in parameter regimes where the microcanonical ensemble is not applicable due to finite-size effects. We discuss our numerical results in the framework of the eigenstate thermalization hypothesis. Moreover, we also observe that the diagonal and the canonical ensembles are practically identical for our initial conditions already on the level of their respective energy distributions for small interaction strengths. Finally, we discuss implications of our results for the interpretation of a recent sudden expansion experiment [Phys. Rev. Lett. 110, 205301 (2013)], in which the same interaction quench was realized.
We consider a quantum many-body system-the Bose-Hubbard system on three sites-which has a classical limit, and which is neither strongly chaotic nor integrable but rather shows a mixture of the two types of behavior. We compare quantum measures of chaos (eigenvalue statistics and eigenvector structure) in the quantum system, with classical measures of chaos (Lyapunov exponents) in the corresponding classical system. As a function of energy and interaction strength, we demonstrate a strong overall correspondence between the two cases. In contrast to both strongly chaotic and integrable systems, the largest Lyapunov exponent is shown to be a multivalued function of energy.
Experimental advances have made ultracold atoms in optical lattices a favourable setting to study out-of-equilibrium phenomena and attracted considerable attention in recent years.These systems are highly versatile in that experimental parameters can be tuned over a wide range of values in real time.This facilitates the study of quantum quenches, in which parameters in the corresponding Hamiltonian are varied in time faster than the system can respond adiabatically.Such protocols open the door to a rich range of many-body physics and have been studied intensely both theoretically and experimentally.The Bose-Hubbard model (BHM) has been shown to describe interacting ultracold bosons in an optical lattice, allowing the opportunity for experiments to probe the out-of-equilibrium dynamics of the model.The BHM is a particularly convenient context for studying quantum quenches as it displays a quantum phase transition between superfluid and Mott-insulator phases.In this thesis, we develop a strong-coupling approach that allows the study of correlations in space and time in both the superfluid and Mott-insulating phases of the BHM.Specifically, we obtain a two-particle irreducible (2PI) effective action within the contourtime formalism that allows for a description of both equilibrium and out-of-equilibrium phenomena.We derive equations of motion for both the superfluid order parameter and the single-particle many-body Green's functions.First, we assess the accuracy of our formalism by studying the equilibrium solution for the homogeneous BHM and comparing our results to existing strong-coupling methods as well as to exact methods where possible.We then consider homogeneous systems that are initially thermalized in the Mott phase, and which are then subjected to quenches.We solve numerically the equations of motion for this scenario and calculate the single-particle density matrix.We demonstrate a Lieb-Robinson-like maximal propagation velocity for the spreading of single-particle correlations in one, two, and three dimensions.We study the dependence of the maximal propagation velocity on the quench protocol, chemical potential, temperature, and dimensionality.We compare our results to exact methods, existing strong-coupling approaches, and experiments where possible.Lastly, we extend our strong-coupling approach to the disordered BHM and derive equations of motion for the disorder-averaged single-particle Green's function.We discuss how these equations of motion can be used to study the phase stability of many-body localization in the disordered BHM for dimensions higher than one.
We numerically study a Bose-Hubbard ring of finite size with disorder containing a finite number of bosons that are subject to an on-site two-body interaction. Our results show that moderate interactions induce dynamical thermalization in this isolated system. In this regime the individual many-body eigenstates are well described by the standard thermal Bose-Einstein distribution for well-defined values of the temperature and the chemical potential, which depend on the eigenstate under consideration. We show that the dynamical thermalization conjecture works well at both positive and negative temperatures. The relations to quantum chaos, quantum ergodicity, and the Åberg criterion are also discussed.
Quantum chaos in isolated quantum systems is intimately linked to thermalization and the rapid relaxation of observables. Although the spectral properties of the chaotic phase in the tilted Bose-Hubbard model have been well characterized, the corresponding dynamical signatures across the transition to regularity remain less explored . In this work, we investigate this transition by analyzing the time evolution of the survival probability, the single-site entanglement entropy, and the half-chain imbalance. Our results reveal a clear hierarchy in the sensitivity of these observables: the relaxation value of the entanglement entropy varies smoothly as a function of the Hamiltonian parameters across the chaos-regular transition, while the imbalance exhibits a more pronounced distinction. Most notably, the survival probability emerges as the most robust indicator of the transition between chaos and regularity. When appropriately scaled, all three observables converge onto a common behavior as a function of the Hamiltonian parameters for different numbers of sites and bosons,enabling a universal characterization of the transition between chaotic and regular dynamics.