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

Quantum quench in the infinitely repulsive Hubbard model: the stationary state

Bruno Bertini, Elena Tartaglia, Pasquale Calabrese

Published Oct 1, 2017• 44 citations
arXivHTMLPDF

Abstract

Abstract We use the quench action approach to study the non-equilibrium dynamics after a quantum quench in the Hubbard model in the limit of infinite interaction. We identify a variety of low-entangled initial states for which we can directly compute the overlaps with the Hamiltonian’s eigenstates. For these initial states, we analytically find the rapidity distributions of the stationary state characterising the expectation values of all local observables. Some of the initial states considered are not reflection symmetric and lead to non-symmetric rapidity distributions. To study such cases, we have to introduce a generalised form for the reduced entropy which measures the entropy restricted to states with non-zero overlap. The initial states considered are of direct experimental realisability and also represent ideal candidates for studying non-equilibrium dynamics in the Hubbard model for finite interactions.

DOI
10.1088/1742-5468/aa8c2c

Related papers

Integrable quenches in the Hubbard model

Colin Rylands, Bruno Bertini, Pasquale Calabrese

Abstract We study the quench dynamics of the one-dimensional Hubbard model through the quench action formalism. We introduce a class of integrable initial states—expressed as product states over two sites—for which we can provide an exact characterisation of the late-time regime. This is achieved by finding a closed-form expression for the overlaps between our states and the Bethe ansatz eigenstates, which we check explicitly in the limits of low densities and infinite repulsion. Our solution gives access to the stationary values attained by local observables (we show the explicit example of the density of doubly occupied sites) and the asymptotic entanglement dynamics directly in the thermodynamic limit. Interestingly, we find that for intermediate interaction strength Rényi entropies display a double-slope structure.

arXiv:2206.07985•Journal of Statistical Mechanics Theory and Experiment•Oct 1, 2022•27 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

On computing non-equilibrium dynamics following a quench

Neil J. Robinson, Albertus J. J. M. de Klerk, Jean-Sébastien Caux

Computing the non-equilibrium dynamics that follows a quantum quench is difficult, even in exactly solvable models. Results are often predicated on the ability to compute overlaps between the initial state and eigenstates of the Hamiltonian that governs time evolution. Except for a handful of known cases, it is generically not possible to find these overlaps analytically. Here we develop a numerical approach to preferentially generate the states with high overlaps for a quantum quench starting from the ground state or an excited state of an initial Hamiltonian. We use these preferentially generated states, in combination with a "high overlap states truncation scheme" and a modification of the numerical renormalization group, to compute non-equilibrium dynamics following a quench in the Lieb-Liniger model. The method is non-perturbative, works for reasonable numbers of particles, and applies to both continuum and lattice systems. It can also be easily extended to more complicated scenarios, including those with integrability breaking.

arXiv:1911.11101•SciPost Physics•Dec 9, 2021•1 citation

On computing non-equilibrium dynamics following a quench

Albertus J. J. M. de Klerk

Computing the non-equilibrium dynamics that follows a quantum quench is difficult, even in exactly solvable models. Results are often predicated on the ability to compute overlaps between the initial state and eigenstates of the Hamiltonian that governs time evolution. Except for a handful of known cases, it is generically not possible to find these overlaps analytically. Here we develop a numerical approach to preferentially generate the states with high overlaps for a quantum quench starting from the ground state or an excited state of an initial Hamiltonian. We use these preferentially generated states, in combination with a "high overlap states truncation scheme" and a modification of the numerical renormalization group, to compute non-equilibrium dynamics following a quench in the Lieb-Liniger model. The method is non-perturbative, works for reasonable numbers of particles, and applies to both continuum and lattice systems. It can also be easily extended to more complicated scenarios, including those with integrability breaking.

•Dec 1, 2021•20 citations

The quench action approach in finite integrable spin chains

Vincenzo Alba, Pasquale Calabrese

Abstract We consider the problem of constructing the stationary state following a quantum quench, using the exact overlaps for finite size integrable models. We focus on the isotropic Heisenberg spin chain with initial state Néel or Majumdar–Ghosh (dimer), although the proposed approach is valid for an arbitrary integrable model. We consider only eigenstates which do not contain zero-momentum strings because the latter are affected by fictitious singularities that are very difficult to take into account. We show that the fraction of eigenstates that do not contain zero-momentum strings is vanishing in the thermodynamic limit. Consequently, restricting to this part of the Hilbert space leads to vanishing expectation values of local observables. However, it is possible to reconstruct the asymptotic values by properly reweighting the expectations in the considered subspace, at the price of introducing finite-size corrections. We also develop a Monte Carlo sampling of the Hilbert space which allows us to study larger systems. We accurately reconstruct the expectation values of the conserved charges and the root distributions in the stationary state, which turn out to match the exact thermodynamic results. The proposed method can be implemented even in cases in which an analytic thermodynamic solution is not obtainable.

arXiv:1512.02213•Journal of Statistical Mechanics Theory and Experiment•Apr 1, 2016•51 citations

Quantum quenches in the sinh-Gordon model: steady state and one-point correlation functions

Bruno Bertini, Lorenzo Piroli, Pasquale Calabrese

Abstract. We consider quantum quenches to the sinh-Gordon integrable quantum field theory from a particular class of initial states. Our analysis includes the case of mass and interaction quenches starting from a non-interacting theory. By means of the recently developed quench action method, we fully characterize the stationary state reached at long times after the quench in terms of the corresponding rapidity distribution. We also provide exact results for the expectation values of arbitrary vertex operators in the post-quench stationary state by proposing a formula based on the analogy with the standard thermodynamic Bethe ansatz. Finally, we comment on the behavior of the post-quench stationary state under the mapping between the sinh-Gordon field theory and the one-dimensional Lieb–Liniger model.

arXiv:1602.08269•Journal of Statistical Mechanics Theory and Experiment•Jun 1, 2016•98 citations
Abstract We use the quench action approach to study the non-equilibrium dynamics after a quantum quench in the Hubbard model in the limit of infinite interaction. We identify a variety of low-entangled initial states for which we can directly compute the overlaps with the Hamiltonian’s eigenstates. For these initial states, we analytically find the rapidity distributions of the stationary state characterising the expectation values of all local observables. Some of the initial states considered are not reflection symmetric and lead to non-symmetric rapidity distributions. To study such cases, we have to introduce a generalised form for the reduced entropy which measures the entropy restricted to states with non-zero overlap. The initial states considered are of direct experimental realisability and also represent ideal candidates for studying non-equilibrium dynamics in the Hubbard model for finite interactions.
Abstract We study the quench dynamics of the one-dimensional Hubbard model through the quench action formalism. We introduce a class of integrable initial states—expressed as product states over two sites—for which we can provide an exact characterisation of the late-time regime. This is achieved by finding a closed-form expression for the overlaps between our states and the Bethe ansatz eigenstates, which we check explicitly in the limits of low densities and infinite repulsion. Our solution gives access to the stationary values attained by local observables (we show the explicit example of the density of doubly occupied sites) and the asymptotic entanglement dynamics directly in the thermodynamic limit. Interestingly, we find that for intermediate interaction strength Rényi entropies display a double-slope structure.
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.
Computing the non-equilibrium dynamics that follows a quantum quench is difficult, even in exactly solvable models. Results are often predicated on the ability to compute overlaps between the initial state and eigenstates of the Hamiltonian that governs time evolution. Except for a handful of known cases, it is generically not possible to find these overlaps analytically. Here we develop a numerical approach to preferentially generate the states with high overlaps for a quantum quench starting from the ground state or an excited state of an initial Hamiltonian. We use these preferentially generated states, in combination with a "high overlap states truncation scheme" and a modification of the numerical renormalization group, to compute non-equilibrium dynamics following a quench in the Lieb-Liniger model. The method is non-perturbative, works for reasonable numbers of particles, and applies to both continuum and lattice systems. It can also be easily extended to more complicated scenarios, including those with integrability breaking.
Computing the non-equilibrium dynamics that follows a quantum quench is difficult, even in exactly solvable models. Results are often predicated on the ability to compute overlaps between the initial state and eigenstates of the Hamiltonian that governs time evolution. Except for a handful of known cases, it is generically not possible to find these overlaps analytically. Here we develop a numerical approach to preferentially generate the states with high overlaps for a quantum quench starting from the ground state or an excited state of an initial Hamiltonian. We use these preferentially generated states, in combination with a "high overlap states truncation scheme" and a modification of the numerical renormalization group, to compute non-equilibrium dynamics following a quench in the Lieb-Liniger model. The method is non-perturbative, works for reasonable numbers of particles, and applies to both continuum and lattice systems. It can also be easily extended to more complicated scenarios, including those with integrability breaking.
Abstract We consider the problem of constructing the stationary state following a quantum quench, using the exact overlaps for finite size integrable models. We focus on the isotropic Heisenberg spin chain with initial state Néel or Majumdar–Ghosh (dimer), although the proposed approach is valid for an arbitrary integrable model. We consider only eigenstates which do not contain zero-momentum strings because the latter are affected by fictitious singularities that are very difficult to take into account. We show that the fraction of eigenstates that do not contain zero-momentum strings is vanishing in the thermodynamic limit. Consequently, restricting to this part of the Hilbert space leads to vanishing expectation values of local observables. However, it is possible to reconstruct the asymptotic values by properly reweighting the expectations in the considered subspace, at the price of introducing finite-size corrections. We also develop a Monte Carlo sampling of the Hilbert space which allows us to study larger systems. We accurately reconstruct the expectation values of the conserved charges and the root distributions in the stationary state, which turn out to match the exact thermodynamic results. The proposed method can be implemented even in cases in which an analytic thermodynamic solution is not obtainable.
Abstract. We consider quantum quenches to the sinh-Gordon integrable quantum field theory from a particular class of initial states. Our analysis includes the case of mass and interaction quenches starting from a non-interacting theory. By means of the recently developed quench action method, we fully characterize the stationary state reached at long times after the quench in terms of the corresponding rapidity distribution. We also provide exact results for the expectation values of arbitrary vertex operators in the post-quench stationary state by proposing a formula based on the analogy with the standard thermodynamic Bethe ansatz. Finally, we comment on the behavior of the post-quench stationary state under the mapping between the sinh-Gordon field theory and the one-dimensional Lieb–Liniger model.