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:2211.10868 · Physical review. B./Physical review. B

Quench dynamics in the one-dimensional mass-imbalanced ionic Hubbard model

Zhuotao Xie, Ming Lei Zhao, Hantao Lu, Zhongbing Huang, Gregory A. Fiete, Xiang Hu, Liang Du

Published May 24, 2023• 2 citations
arXivHTMLPDF

Abstract

Using the time-dependent Lanczos method, we study the nonequilibrium dynamics of the one-dimensional ionic-mass imbalanced Hubbard chain driven by a quantum quench of the on-site Coulomb interaction, where the system is prepared in the ground state of the Hamiltonian with a different Hubbard interaction. The exact diagonalization method (Lanczos algorithm) is adopted to study the zero temperature phase diagram in equilibrium, which is shown to be in good agreement with previous studies using density matrix renormalization group (DMRG). We then study the nonequilibrium quench dynamics of the spin and charge order parameters by fixing the initial and final Coulomb interactions while changing the quenching time protocols. Our study shows that the time evolution of the charge and spin order parameters strongly depend on the quenching time protocols. In particular, the effective temperature of the system will decrease monotonically as the quenching time is increased. By taking the final Coulomb interaction strength to be in the strong coupling regime, we find that the oscillation frequency of the charge order parameter increases monotonically with the Coulomb interaction. By contrast, the frequency of the spin order parameter decreases monotonically with increasing Coulomb interaction. We explain this result using an effective spin model and a two-site Hubbard model in the strong coupling limit. Finally, we take the final Coulomb interaction strength to be in the weak coupling regime and find that the oscillation frequency of both the charge and spin order parameters increases monotonically with decreasing Coulomb interaction. Our study suggests strategies to engineer the relaxation behavior of interacting quantum many-particle systems.

DOI
10.1103/physrevb.107.195147

Related papers

Dynamical recovery of SU(2) symmetry in the mass-quenched Hubbard model

Liang Du, Gregory A. Fiete

We use nonequilibrium dynamical mean-field theory with iterative perturbation theory as an impurity solver to study the recovery of SU(2) symmetry in real time following a hopping integral parameter quench from a mass-imbalanced to a mass-balanced single-band Hubbard model at half filling. A dynamical order parameter $\ensuremath{\gamma}(t)$ is defined to characterize the evolution of the system towards SU(2) symmetry. By comparing the momentum-dependent occupation from an equilibrium calculation [with the SU(2) symmetric Hamiltonian after the quench at an effective temperature] with the data from our nonequilibrium calculation, we conclude that the SU(2) symmetry recovered state is a thermalized state. Further evidence from the evolution of the density of states supports this conclusion. We find the order parameter in the weak Coulomb interaction regime undergoes an approximate exponential decay. We numerically investigate the interplay of the relevant parameters (initial temperature, Coulomb interaction strength, initial mass-imbalance ratio) and their combined effect on the thermalization behavior. Finally, we study evolution of the order parameter as the hopping parameter is changed with either a linear ramp or a pulse. Our results can be useful in strategies to engineer the relaxation behavior of interacting quantum many-particle systems.

arXiv:1711.05390•Physical review. B./Physical review. B•Feb 27, 2018•3 citations

Spin-selective thermalization plateau in the mass-imbalanced Hubbard model

Liang Du, Li Huang, Gregory A. Fiete

We use time-dependent nonequilibrium dynamical mean-field theory with weak-coupling auxiliary-field continuous time quantum Monte Carlo as an impurity solver to study the thermalization behavior of the mass-imbalanced single-band Hubbard model after a quench of the Coulomb interaction from the noninteracting limit to a finite positive value. When the Coulomb interaction in our model is increased under equilibrium conditions, the quasiparticle weights for spin-up and spin-down (the mass imbalance) electrons approach zero simultaneously, indicating the absence of a spin-selective Mott transition. By contrast, our out-of-equilibrium study of the mass-imbalanced Hubbard model suggests that there exists the spin-selective dynamical phase transition (one spin orientation undergoes a fast thermalization at its critical Coulomb interaction strength while the other spin orientation shows prethermalization behavior). The spin-selective dynamical phase transition is characterized by the relaxation behavior of the spin-resolved kinetic energy and the spin-resolved momentum-dependent occupation. To connect with possible experiments, we calculate the spin-resolved two-time optical conductivity, which confirms the spin-selective thermalization plateau. We find that the critical Coulomb interaction of each spin orientation for the spin-selective thermalization grows as the mass imbalance decreases.

arXiv:1707.09692•Physical review. B./Physical review. B•Oct 30, 2017•6 citations

Quantum critical behavior of the one-dimensional ionic Hubbard model

Salvatore R. Manmana, V. Meden, R. M. Noack +1 more

We study the zero-temperature phase diagram of the half-filled one-dimensional ionic Hubbard model. This model is governed by the interplay of the on-site Coulomb repulsion and an alternating one-particle potential. Various many-body energy gaps, the charge-density-wave and bond-order parameters, the electric as well as the bond-order susceptibilities, and the density-density correlation function are calculated using the density-matrix renormalization group method. In order to obtain a comprehensive picture, we investigate systems with open as well as periodic boundary conditions and study the physical properties in different sectors of the phase diagram. A careful finite-size scaling analysis leads to results which give evidence in favor of a scenario with two quantum critical points and an intermediate spontaneously dimerized phase. Our results indicate that the phase transitions are continuous. Using a scaling ansatz we are able to read off critical exponents at the first critical point. In contrast to a bosonization approach, we do not find Ising critical exponents. We show that the low-energy physics of the strong-coupling phase can only partly be understood in terms of the strong-coupling behavior of the ordinary Hubbard model.

arXiv:cond-mat/0307741•Physical Review B•Oct 29, 2004•117 citations

Time-resolved single-particle spectrum of the one-dimensional extended Hubbard model after interaction quenches

Yong-Guang Su, Ruifeng Lu, Hantao Lu +1 more

Abstract We investigate the non-equilibrium dynamics of the one-dimensional extended Hubbard model after interaction quenches. In strong-coupling regime with large on-site interaction, the ground states of this model with small and large nearest-neighbor interactions are in spin-density-wave and charge-density-wave phases, respectively. Combining twisted boundary conditions with the time-dependent Lanczos method, we obtain snapshots of the time-dependent single-particle spectrum after quenches. We find that for quench within the same phase, the single-particle spectrum becomes close to that of the quenched Hamiltonian immediately after the quench. While for quench across the critical point, the afterward evolution process depends mainly on the distribution of the initial state among the eigenstates of the quenched Hamiltonian. Our finding may serve as a way to detect the phase transition in ultracold atom systems with interactions.

arXiv:2303.12322•Journal of Physics B Atomic Molecular and Optical Physics•Mar 15, 2023•1 citation

Quench dynamics of the spin-imbalanced Fermi-Hubbard model in one dimension

Xiao Ling Yin, Leo Radzihovsky

We study a nonequilibrium dynamics of a one-dimensional spin-imbalanced Fermi-Hubbard model following a quantum quench of on-site interaction, realizable, for example, in Feshbach-resonant atomic Fermi gases. We focus on the post-quench evolution starting from the initial BCS and Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) ground states and analyze the corresponding spin-singlet, spin-triplet, density-density, and magnetization-magnetization correlation functions. We find that beyond a light-cone crossover time, rich post-quench dynamics leads to thermalized and pre-thermalized stationary states that display strong dependence on the initial ground state. For initially gapped BCS state, the long-time stationary state resembles thermalization with the effective temperature set by the initial value of the Hubbard interaction. In contrast, while the initial gapless FFLO state reaches a stationary pre-thermalized form, it remains far from equilibrium. We suggest that such post-quench dynamics can be used as a fingerprint for identification and study of the FFLO phase.

arXiv:1610.00670•Physical Review A•Dec 27, 2016•17 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
Using the time-dependent Lanczos method, we study the nonequilibrium dynamics of the one-dimensional ionic-mass imbalanced Hubbard chain driven by a quantum quench of the on-site Coulomb interaction, where the system is prepared in the ground state of the Hamiltonian with a different Hubbard interaction. The exact diagonalization method (Lanczos algorithm) is adopted to study the zero temperature phase diagram in equilibrium, which is shown to be in good agreement with previous studies using density matrix renormalization group (DMRG). We then study the nonequilibrium quench dynamics of the spin and charge order parameters by fixing the initial and final Coulomb interactions while changing the quenching time protocols. Our study shows that the time evolution of the charge and spin order parameters strongly depend on the quenching time protocols. In particular, the effective temperature of the system will decrease monotonically as the quenching time is increased. By taking the final Coulomb interaction strength to be in the strong coupling regime, we find that the oscillation frequency of the charge order parameter increases monotonically with the Coulomb interaction. By contrast, the frequency of the spin order parameter decreases monotonically with increasing Coulomb interaction. We explain this result using an effective spin model and a two-site Hubbard model in the strong coupling limit. Finally, we take the final Coulomb interaction strength to be in the weak coupling regime and find that the oscillation frequency of both the charge and spin order parameters increases monotonically with decreasing Coulomb interaction. Our study suggests strategies to engineer the relaxation behavior of interacting quantum many-particle systems.
We use nonequilibrium dynamical mean-field theory with iterative perturbation theory as an impurity solver to study the recovery of SU(2) symmetry in real time following a hopping integral parameter quench from a mass-imbalanced to a mass-balanced single-band Hubbard model at half filling. A dynamical order parameter $\ensuremath{\gamma}(t)$ is defined to characterize the evolution of the system towards SU(2) symmetry. By comparing the momentum-dependent occupation from an equilibrium calculation [with the SU(2) symmetric Hamiltonian after the quench at an effective temperature] with the data from our nonequilibrium calculation, we conclude that the SU(2) symmetry recovered state is a thermalized state. Further evidence from the evolution of the density of states supports this conclusion. We find the order parameter in the weak Coulomb interaction regime undergoes an approximate exponential decay. We numerically investigate the interplay of the relevant parameters (initial temperature, Coulomb interaction strength, initial mass-imbalance ratio) and their combined effect on the thermalization behavior. Finally, we study evolution of the order parameter as the hopping parameter is changed with either a linear ramp or a pulse. Our results can be useful in strategies to engineer the relaxation behavior of interacting quantum many-particle systems.
We use time-dependent nonequilibrium dynamical mean-field theory with weak-coupling auxiliary-field continuous time quantum Monte Carlo as an impurity solver to study the thermalization behavior of the mass-imbalanced single-band Hubbard model after a quench of the Coulomb interaction from the noninteracting limit to a finite positive value. When the Coulomb interaction in our model is increased under equilibrium conditions, the quasiparticle weights for spin-up and spin-down (the mass imbalance) electrons approach zero simultaneously, indicating the absence of a spin-selective Mott transition. By contrast, our out-of-equilibrium study of the mass-imbalanced Hubbard model suggests that there exists the spin-selective dynamical phase transition (one spin orientation undergoes a fast thermalization at its critical Coulomb interaction strength while the other spin orientation shows prethermalization behavior). The spin-selective dynamical phase transition is characterized by the relaxation behavior of the spin-resolved kinetic energy and the spin-resolved momentum-dependent occupation. To connect with possible experiments, we calculate the spin-resolved two-time optical conductivity, which confirms the spin-selective thermalization plateau. We find that the critical Coulomb interaction of each spin orientation for the spin-selective thermalization grows as the mass imbalance decreases.
We study the zero-temperature phase diagram of the half-filled one-dimensional ionic Hubbard model. This model is governed by the interplay of the on-site Coulomb repulsion and an alternating one-particle potential. Various many-body energy gaps, the charge-density-wave and bond-order parameters, the electric as well as the bond-order susceptibilities, and the density-density correlation function are calculated using the density-matrix renormalization group method. In order to obtain a comprehensive picture, we investigate systems with open as well as periodic boundary conditions and study the physical properties in different sectors of the phase diagram. A careful finite-size scaling analysis leads to results which give evidence in favor of a scenario with two quantum critical points and an intermediate spontaneously dimerized phase. Our results indicate that the phase transitions are continuous. Using a scaling ansatz we are able to read off critical exponents at the first critical point. In contrast to a bosonization approach, we do not find Ising critical exponents. We show that the low-energy physics of the strong-coupling phase can only partly be understood in terms of the strong-coupling behavior of the ordinary Hubbard model.
Abstract We investigate the non-equilibrium dynamics of the one-dimensional extended Hubbard model after interaction quenches. In strong-coupling regime with large on-site interaction, the ground states of this model with small and large nearest-neighbor interactions are in spin-density-wave and charge-density-wave phases, respectively. Combining twisted boundary conditions with the time-dependent Lanczos method, we obtain snapshots of the time-dependent single-particle spectrum after quenches. We find that for quench within the same phase, the single-particle spectrum becomes close to that of the quenched Hamiltonian immediately after the quench. While for quench across the critical point, the afterward evolution process depends mainly on the distribution of the initial state among the eigenstates of the quenched Hamiltonian. Our finding may serve as a way to detect the phase transition in ultracold atom systems with interactions.
We study a nonequilibrium dynamics of a one-dimensional spin-imbalanced Fermi-Hubbard model following a quantum quench of on-site interaction, realizable, for example, in Feshbach-resonant atomic Fermi gases. We focus on the post-quench evolution starting from the initial BCS and Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) ground states and analyze the corresponding spin-singlet, spin-triplet, density-density, and magnetization-magnetization correlation functions. We find that beyond a light-cone crossover time, rich post-quench dynamics leads to thermalized and pre-thermalized stationary states that display strong dependence on the initial ground state. For initially gapped BCS state, the long-time stationary state resembles thermalization with the effective temperature set by the initial value of the Hubbard interaction. In contrast, while the initial gapless FFLO state reaches a stationary pre-thermalized form, it remains far from equilibrium. We suggest that such post-quench dynamics can be used as a fingerprint for identification and study of the FFLO phase.
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.