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

Analysis of self-thermalization dynamics in the Bose-Hubbard model by using the pseudoclassical approach

Andrey R. Kolovsky

quant-phSubmitted Jan 30, 2026• Announced Feb 2, 2026• 0 citations
v1 · Jan 30, 2026v2 · Sep 19, 2026
arXivHTMLPDF

Abstract

We analyze the self-thermalization dynamics of the $M$-site Bose-Hubbard model in terms of the single-particle density matrix that is calculated by using the pseudoclassical approach. It is shown that a weak inter-particle interaction, which suffices to convert the integrable system of non-interacting bosons into a chaotic system, has a negligible effect on the thermal density matrix given by the Bose-Einstein distribution. This opens the door for equilibration where the two coupled Bose-Hubbard systems, which are initially in different thermal states, relax to the same thermal state. When we couple these two subsystems by using a lattice of the length $L\ll M$, we numerically calculate the quasi-stationary current of Bose particles across the lattice and show that its magnitude is consistent with the solution of the master equation for the boundary driven $L$-site Bose-Hubbard model.

Comments
7 pages, 6 figures

Related papers

Open Bose-Hubbard chain: Pseudoclassical approach

Anna A. Bychek, Pavel S. Muraev, Дмитрий Н. Максимов +1 more

We analyze the stationary current of bosonic carriers in the Bose-Hubbard chain of length L where the first and the last sites of the chain are attached to reservoirs of Bose particles acting as a particle source and sink, respectively. The analysis is carried out by using the pseudoclassical approach which reduces the original quantum problem to the classical problem for L coupled nonlinear oscillators. It is shown that an increase of oscillator nonlinearity (which is determined by the strength of interparticle interactions) results in a transition from the ballistic transport regime, where the stationary current is independent of the chain length, to the diffusive regime, where the current is inversely proportional to L.

arXiv:1909.06213•Physical review. E•Jan 13, 2020•19 citations

Parametric Bose-Hubbard Hamiltonians: Quantum Dissipation, Irreversibility, and Pumping

Moritz Hiller

In this thesis, we study interacting bosons on lattices consisting of a few sites. The main focus is on their response to external driving fields, their transport and decay properties. We address fundamental questions of quantum classical correspondence and investigate the implications of interactions in the transport properties of (ultra-)cold atom devices. Quantum mechanically, we describe bosons on lattices using a parametric Bose-Hubbard Hamiltonian while the corresponding classical limit is described by the discrete nonlinear Schrödinger equation.In the first part of this thesis we focus on the chaotic regime of the Bose-Hubbard Hamiltonian. We investigate the parametric change of eigenstates as the coupling strength between neighboring sites is varied, using perturbative and semiclassical methods. We then turn to the quantum evolution of driven interacting bosonic systems and study the resulting quantum dynamics and the irreversibility (fidelity) of the quantum motion.In the second part we initiate the study of quantum pumping/stirring in Bose-Einstein condensates with the aim to identify optimal pumping cycles and propose such a device in order to probe the interatomic interactions. Finally, in the last part we focus on the interplay of intrinsic dynamics with coupling to the continuum and investigate the structure of the resonance widths of a Bose-Hubbard dimer which is coupled to the continuum at one of the sites using an effective non-Hermitian Hamiltonian formalism.

•Jan 1, 2008•3 citations

Entanglement and thermalization in the extended Bose–Hubbard model after a quantum quench: A correlation analysis

Xiaoqiang Su, Zong-Ju Xu, You-Quan Zhao

Exploring the role of entanglement in quantum nonequilibrium dynamics is important to understand the mechanism of thermalization in an isolated system. We study the relaxation dynamics in a one-dimensional extended Bose–Hubbard model after a global interaction quench by considering several observables: the local Boson numbers, the nonlocal entanglement entropy, and the momentum distribution functions. We calculate the thermalization fidelity for different quench parameters and different sizes of subsystems, and the results show that the degree of thermalization is affected by the distance from the integrable point and the size of the subsystem. We employ the Pearson coefficient as the measurement of the correlation between the entanglement entropy and thermalization fidelity, and a strong correlation is demonstrated for the quenched system.

Chinese Physics B•Jun 8, 2022•2 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

Expansion dynamics in two-dimensional Bose-Hubbard lattices: Bose-Einstein condensate and thermal cloud

Mauricio Trujillo-Martinez, Anna Posazhennikova, Johann Kroha

We study the temporal expansion of an ultracold Bose gas in two-dimensional square optical lattices. The gas is described by the Bose-Hubbard model deep in the superfluid regime, with initially all bosons condensed in the central site of the lattice. We use the previously developed nonequilibrium propagator method for capturing the time evolution of an interacting bosonic system, where the many-body Hamiltonian is represented in an appropriate local basis and the corresponding field operators are separated into the classical [Bose-Einstein condensate (BEC)] part and quantum mechanical fluctuations. After a quench, i.e., after a sudden switch of the lattice nearest-neighbor hopping, the expanding bosonic cloud separates spatially into a fast ballistic forerunner and a slowly expanding central part controlled by self-trapping. We show that the forerunner expansion is driven by the coherent dynamics of the BEC and that its velocity is consistent with the Lieb-Robinson bound. For smaller lattices we analyze how quasiparticle collisions lead to enhanced condensate depletion and oscillation damping.

arXiv:2007.06331•Physical Review A•Mar 15, 2021•8 citations

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
We analyze the self-thermalization dynamics of the MMM-site Bose-Hubbard model in terms of the single-particle density matrix that is calculated by using the pseudoclassical approach. It is shown that a weak inter-particle interaction, which suffices to convert the integrable system of non-interacting bosons into a chaotic system, has a negligible effect on the thermal density matrix given by the Bose-Einstein distribution. This opens the door for equilibration where the two coupled Bose-Hubbard systems, which are initially in different thermal states, relax to the same thermal state. When we couple these two subsystems by using a lattice of the length L≪ML\ll ML≪M, we numerically calculate the quasi-stationary current of Bose particles across the lattice and show that its magnitude is consistent with the solution of the master equation for the boundary driven LLL-site Bose-Hubbard model.
We analyze the stationary current of bosonic carriers in the Bose-Hubbard chain of length L where the first and the last sites of the chain are attached to reservoirs of Bose particles acting as a particle source and sink, respectively. The analysis is carried out by using the pseudoclassical approach which reduces the original quantum problem to the classical problem for L coupled nonlinear oscillators. It is shown that an increase of oscillator nonlinearity (which is determined by the strength of interparticle interactions) results in a transition from the ballistic transport regime, where the stationary current is independent of the chain length, to the diffusive regime, where the current is inversely proportional to L.
In this thesis, we study interacting bosons on lattices consisting of a few sites. The main focus is on their response to external driving fields, their transport and decay properties. We address fundamental questions of quantum classical correspondence and investigate the implications of interactions in the transport properties of (ultra-)cold atom devices. Quantum mechanically, we describe bosons on lattices using a parametric Bose-Hubbard Hamiltonian while the corresponding classical limit is described by the discrete nonlinear Schrödinger equation.In the first part of this thesis we focus on the chaotic regime of the Bose-Hubbard Hamiltonian. We investigate the parametric change of eigenstates as the coupling strength between neighboring sites is varied, using perturbative and semiclassical methods. We then turn to the quantum evolution of driven interacting bosonic systems and study the resulting quantum dynamics and the irreversibility (fidelity) of the quantum motion.In the second part we initiate the study of quantum pumping/stirring in Bose-Einstein condensates with the aim to identify optimal pumping cycles and propose such a device in order to probe the interatomic interactions. Finally, in the last part we focus on the interplay of intrinsic dynamics with coupling to the continuum and investigate the structure of the resonance widths of a Bose-Hubbard dimer which is coupled to the continuum at one of the sites using an effective non-Hermitian Hamiltonian formalism.
Exploring the role of entanglement in quantum nonequilibrium dynamics is important to understand the mechanism of thermalization in an isolated system. We study the relaxation dynamics in a one-dimensional extended Bose–Hubbard model after a global interaction quench by considering several observables: the local Boson numbers, the nonlocal entanglement entropy, and the momentum distribution functions. We calculate the thermalization fidelity for different quench parameters and different sizes of subsystems, and the results show that the degree of thermalization is affected by the distance from the integrable point and the size of the subsystem. We employ the Pearson coefficient as the measurement of the correlation between the entanglement entropy and thermalization fidelity, and a strong correlation is demonstrated for the quenched system.
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 study the temporal expansion of an ultracold Bose gas in two-dimensional square optical lattices. The gas is described by the Bose-Hubbard model deep in the superfluid regime, with initially all bosons condensed in the central site of the lattice. We use the previously developed nonequilibrium propagator method for capturing the time evolution of an interacting bosonic system, where the many-body Hamiltonian is represented in an appropriate local basis and the corresponding field operators are separated into the classical [Bose-Einstein condensate (BEC)] part and quantum mechanical fluctuations. After a quench, i.e., after a sudden switch of the lattice nearest-neighbor hopping, the expanding bosonic cloud separates spatially into a fast ballistic forerunner and a slowly expanding central part controlled by self-trapping. We show that the forerunner expansion is driven by the coherent dynamics of the BEC and that its velocity is consistent with the Lieb-Robinson bound. For smaller lattices we analyze how quasiparticle collisions lead to enhanced condensate depletion and oscillation damping.
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.