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:1405.1312 · EPJ Quantum Technology

Propagation of quantum correlations after a quench in the Mott-insulator regime of the Bose-Hubbard model

Konstantin V Krutitsky, Patrick Navez, Friedemann Queisser, Ralf Schützhold

Published Sep 3, 2014• 37 citations
arXivHTMLPDF

Abstract

We study a quantum quench in the Bose-Hubbard model where the tunneling rate J is suddenly switched from zero to a finite value in the Mott regime. In order to solve the many-body quantum dynamics far from equilibrium, we consider the reduced density matrices for a finite number of lattice sites and split them up into on-site density operators, i.e., the mean field, plus two-point and three-point correlations etc. Neglecting three-point and higher correlations, we are able to numerically simulate the time-evolution of the on-site density matrices and the two-point quantum correlations (e.g., their effective light-cone structure) for a comparably large number of lattice sites.

DOI
10.1140/epjqt12

Related papers

Equilibration and prethermalization in the Bose-Hubbard and Fermi-Hubbard models

Friedemann Queißer, Konstantin V. Krutitsky, Patrick Navez +1 more

We study the Bose-Hubbard and Fermi-Hubbard models in the (formal) limit of large coordination numbers $Z\ensuremath{\gg}1$. Via an expansion into powers of $1/Z$, we establish a hierarchy of correlations which facilitates an approximate analytical derivation of the time evolution of the reduced density matrices for one and two sites, etc. With this method, we study the quantum dynamics (starting in the ground state) after a quantum quench, i.e., after suddenly switching the tunneling rate $J$ from zero to a finite value, which is still in the Mott regime. We find that the reduced density matrices approach a (quasi)equilibrium state after some time. For one lattice site, this state can be described by a thermal state (within the accuracy of our approximation). However, the (quasi)equilibrium state of the reduced density matrices for two sites including the correlations can not be described by a thermal state. Thus, real thermalization (if it occurs) should take a much longer time. This behavior has already been observed in other scenarios and is sometimes called ``prethermalization''. Finally, we compare our results to numerical simulations for finite lattices in one and two dimensions and find qualitative agreement.

arXiv:1311.2212•Physical Review A•Mar 12, 2014•51 citations

Light-Cone and Diffusive Propagation of Correlations in a Many-Body Dissipative System

Jean-Sébastien Bernier, Ryan Tan, Lars Bonnes +3 more

We analyze the propagation of correlations after a sudden interaction change in a strongly interacting quantum system in contact with an environment. In particular, we consider an interaction quench in the Bose-Hubbard model, deep within the Mott-insulating phase, under the effect of dephasing. We observe that dissipation effectively speeds up the propagation of single-particle correlations while reducing their coherence. In contrast, for two-point density correlations, the initial ballistic propagation regime gives way to diffusion at intermediate times. Numerical simulations, based on a time-dependent matrix product state algorithm, are supplemented by a quantitatively accurate fermionic quasiparticle approach providing an intuitive description of the initial dynamics in terms of holon and doublon excitations.

arXiv:1702.04136•Physical Review Letters•Jan 12, 2018•40 citations

Semiclassical quench dynamics of Bose gases in optical lattices

Kazuma Nagao, Masaya Kunimi, Yosuke Takasu +2 more

We analyze the time evolution of the Bose-Hubbard model after a sudden quantum quench to a weakly interacting regime. Specifically, motivated by a recent experiment at Kyoto University, we numerically simulate redistribution of the kinetic and on-site interaction energies at an early time, which was observed in nonequilibrium dynamics of ultracold Bose gases in a cubic optical lattice starting with a singly occupied Mott-insulator state. In order to compute the short-time dynamics corresponding to the experimental situation, we apply the truncated-Wigner approximation to the Bose-Hubbard model on a cubic lattice. We show that our semiclassical approach quantitatively reproduces the fast redistribution dynamics. We further analyze spatial spreading of density-density correlations at equal time in the Bose-Hubbard model on a square lattice with a large filling factor. When the system is initially prepared in a coherent state, we find that a propagation velocity of the correlation wave packet in the correlation function strongly depends on the final interaction strength, and it is bounded by twice the maximum group velocity of the elementary excitations. In contrast, when the system is initially in a Mott-insulator state, the propagation velocity of the wave packet is approximately independent of the final interaction strength.

arXiv:1810.04347•Physical Review A•Feb 20, 2019•21 citations

Equilibration versus (pre) thermalization in the Bose and Fermi Hubbard models

Friedemann Queißer, Konstantin V. Krutitsky, Patrick Navez +1 more

We study the Bose and Fermi Hubbard model in the (formal) limit of large coordination numbers $Z\gg1$. Via an expansion into powers of $1/Z$, we establish a hierarchy of correlations which facilitates an approximate analytical derivation of the time-evolution of the reduced density matrices for one and two sites etc. With this method, we study the quantum dynamics (starting in the ground state) after a quantum quench, i.e., after suddenly switching the tunneling rate $J$ from zero to a finite value, which is still in the Mott regime. We find that the reduced density matrices approach a (quasi) equilibrium state after some time. For one lattice site, this state can be described by a thermal state (within the accuracy of our approximation). However, the (quasi) equilibrium state of the reduced density matrices for two sites including the correlations cannot be described by a thermal state. Thus, real thermalization (if it occurs) should take much longer time. This behavior has already been observed in other scenarios and is sometimes called ``pre-thermalization.'' Finally, we compare our results to numerical simulations for finite lattices in one and two dimensions and find qualitative agreement.

•Nov 9, 2013•0 citations

Space-time correlations in the Bose Hubbard model after a quantum quench

Matthew Fitzpatrick, Malcolm P. Kennett

We study correlations in space and time after a quantum quench in the Bose Hubbard model. We obtain a two particle irreducible effective action within the Schwinger-Keldysh formalism that allows us to study dynamics in both the superfluid and Mott insulating phases. We derive equations of motion for both the superfluid order parameter and two-point correlation functions. We discuss the numerical solutions of these equations and calculate the quasimomentum distribution for quenches in the Mott phase. We demonstrate light-cone like spreading of correlations in the Mott phase in both one and two dimensions and calculate the propagation velocity in each dimension.

•Jun 13, 2016•1 citation

Quantum quench dynamics of the Bose-Hubbard model at finite temperatures

J. M. Zhang, Chao Shen, W. M. Liu

We study quench dynamics of the Bose-Hubbard model by exact diagonalization. Initially, the system is at thermal equilibrium and of a finite temperature. The system is then quenched by changing the on-site interaction strength $U$ suddenly. Both the single-quench and double-quench scenarios are considered. In the former case, the time-averaged density matrix and the real-time evolution are investigated. It is found that though the system thermalizes only in a very narrow range of the quenched value of $U$, it does equilibrate or relax well into a much larger range. Most importantly, it is proven that this is guaranteed for some typical observables in the thermodynamic limit. In order to test whether it is possible to distinguish the unitarily evolving density matrix from the time-averaged (thus time-independent), fully decohered density matrix, a second quench is considered. It turns out that the answer is affirmative or negative depending on whether the intermediate value of $U$ is zero or not.

arXiv:1103.1539•Physical Review A•Jun 17, 2011•17 citations
We study a quantum quench in the Bose-Hubbard model where the tunneling rate J is suddenly switched from zero to a finite value in the Mott regime. In order to solve the many-body quantum dynamics far from equilibrium, we consider the reduced density matrices for a finite number of lattice sites and split them up into on-site density operators, i.e., the mean field, plus two-point and three-point correlations etc. Neglecting three-point and higher correlations, we are able to numerically simulate the time-evolution of the on-site density matrices and the two-point quantum correlations (e.g., their effective light-cone structure) for a comparably large number of lattice sites.
We study the Bose-Hubbard and Fermi-Hubbard models in the (formal) limit of large coordination numbers $Z\ensuremath{\gg}1$. Via an expansion into powers of 1/Z1/Z1/Z, we establish a hierarchy of correlations which facilitates an approximate analytical derivation of the time evolution of the reduced density matrices for one and two sites, etc. With this method, we study the quantum dynamics (starting in the ground state) after a quantum quench, i.e., after suddenly switching the tunneling rate JJJ from zero to a finite value, which is still in the Mott regime. We find that the reduced density matrices approach a (quasi)equilibrium state after some time. For one lattice site, this state can be described by a thermal state (within the accuracy of our approximation). However, the (quasi)equilibrium state of the reduced density matrices for two sites including the correlations can not be described by a thermal state. Thus, real thermalization (if it occurs) should take a much longer time. This behavior has already been observed in other scenarios and is sometimes called ``prethermalization''. Finally, we compare our results to numerical simulations for finite lattices in one and two dimensions and find qualitative agreement.
We analyze the propagation of correlations after a sudden interaction change in a strongly interacting quantum system in contact with an environment. In particular, we consider an interaction quench in the Bose-Hubbard model, deep within the Mott-insulating phase, under the effect of dephasing. We observe that dissipation effectively speeds up the propagation of single-particle correlations while reducing their coherence. In contrast, for two-point density correlations, the initial ballistic propagation regime gives way to diffusion at intermediate times. Numerical simulations, based on a time-dependent matrix product state algorithm, are supplemented by a quantitatively accurate fermionic quasiparticle approach providing an intuitive description of the initial dynamics in terms of holon and doublon excitations.
We analyze the time evolution of the Bose-Hubbard model after a sudden quantum quench to a weakly interacting regime. Specifically, motivated by a recent experiment at Kyoto University, we numerically simulate redistribution of the kinetic and on-site interaction energies at an early time, which was observed in nonequilibrium dynamics of ultracold Bose gases in a cubic optical lattice starting with a singly occupied Mott-insulator state. In order to compute the short-time dynamics corresponding to the experimental situation, we apply the truncated-Wigner approximation to the Bose-Hubbard model on a cubic lattice. We show that our semiclassical approach quantitatively reproduces the fast redistribution dynamics. We further analyze spatial spreading of density-density correlations at equal time in the Bose-Hubbard model on a square lattice with a large filling factor. When the system is initially prepared in a coherent state, we find that a propagation velocity of the correlation wave packet in the correlation function strongly depends on the final interaction strength, and it is bounded by twice the maximum group velocity of the elementary excitations. In contrast, when the system is initially in a Mott-insulator state, the propagation velocity of the wave packet is approximately independent of the final interaction strength.
We study the Bose and Fermi Hubbard model in the (formal) limit of large coordination numbers Z≫1Z\gg1Z≫1. Via an expansion into powers of 1/Z1/Z1/Z, we establish a hierarchy of correlations which facilitates an approximate analytical derivation of the time-evolution of the reduced density matrices for one and two sites etc. With this method, we study the quantum dynamics (starting in the ground state) after a quantum quench, i.e., after suddenly switching the tunneling rate JJJ from zero to a finite value, which is still in the Mott regime. We find that the reduced density matrices approach a (quasi) equilibrium state after some time. For one lattice site, this state can be described by a thermal state (within the accuracy of our approximation). However, the (quasi) equilibrium state of the reduced density matrices for two sites including the correlations cannot be described by a thermal state. Thus, real thermalization (if it occurs) should take much longer time. This behavior has already been observed in other scenarios and is sometimes called ``pre-thermalization.'' Finally, we compare our results to numerical simulations for finite lattices in one and two dimensions and find qualitative agreement.
We study correlations in space and time after a quantum quench in the Bose Hubbard model. We obtain a two particle irreducible effective action within the Schwinger-Keldysh formalism that allows us to study dynamics in both the superfluid and Mott insulating phases. We derive equations of motion for both the superfluid order parameter and two-point correlation functions. We discuss the numerical solutions of these equations and calculate the quasimomentum distribution for quenches in the Mott phase. We demonstrate light-cone like spreading of correlations in the Mott phase in both one and two dimensions and calculate the propagation velocity in each dimension.
We study quench dynamics of the Bose-Hubbard model by exact diagonalization. Initially, the system is at thermal equilibrium and of a finite temperature. The system is then quenched by changing the on-site interaction strength UUU suddenly. Both the single-quench and double-quench scenarios are considered. In the former case, the time-averaged density matrix and the real-time evolution are investigated. It is found that though the system thermalizes only in a very narrow range of the quenched value of UUU, it does equilibrate or relax well into a much larger range. Most importantly, it is proven that this is guaranteed for some typical observables in the thermodynamic limit. In order to test whether it is possible to distinguish the unitarily evolving density matrix from the time-averaged (thus time-independent), fully decohered density matrix, a second quench is considered. It turns out that the answer is affirmative or negative depending on whether the intermediate value of UUU is zero or not.