By combining the Baeriswyl wave function with equilibrium and time-dependent variational principles, we develop a nonequilibrium formalism to study quantum quenches for two-dimensional spinless fermions with nearest-neighbor hopping and repulsion. The variational ground-state energy, the charge-density wave (CDW) order parameter, and the short-time dynamics agree convincingly with the results of numerically exact simulations. We find that, depending on the initial and final interaction strength, the quenched system either exhibits oscillatory behavior or relaxes to a time-independent steady state. The time-averaged expectation value of the CDW order parameter rises sharply when crossing from the steady-state regime to the oscillating regime, indicating that the system, being nonintegrable, shows signs of thermalization with an effective temperature above or below the equilibrium critical temperature, respectively.
By combining the Baeriswyl wave function with equilibrium and time-dependent variational principles, we develop a nonequilibrium formalism to study quantum quenches for two-dimensional spinless fermions with nearest-neighbor hopping and repulsion. The variational ground-state energy, the charge-density wave (CDW) order parameter, and the short-time dynamics agree convincingly with the results of numerically exact simulations. We find that, depending on the initial and final interaction strength, the quenched system either exhibits oscillatory behavior or relaxes to a time-independent steady state. The time-averaged expectation value of the CDW order parameter rises sharply when crossing from the steady-state regime to the oscillating regime, indicating that the system, being nonintegrable, shows signs of thermalization with an effective temperature above or below the equilibrium critical temperature, respectively.
We investigate the formation of charge and spin ordering by starting from a noninteracting state and studying how it evolves in time under a Hamiltonian with finite electronic interactions. We consider the one-dimensional, half-filled extended Hubbard model, which we solve within time-dependent density matrix renormalization group. By employing linear finite-time quenches in on-site and nearest-neighbor interactions, we find the existence of impulse, intermediate, and adiabatic regimes of time evolution. For the quenches we analyze, we observe that the adiabatic regime is reached with distinct ramping timescales depending on whether the charge density wave (CDW) or the spin density wave (SDW) is formed. The former needs to be slower than the latter to prevent entangled excited states from being accessed during the quench. More interestingly, in the intermediate regime, we observe an enhancement of the entanglement entropy with respect to its initial value, which precedes the formation of the CDW ordering; a similar enhancement is not seen in the quench toward SDW. Our findings also show that the breaking of the system integrability, by turning on the nearest-neighbor interactions, does not give rise to significant changes in the nonequilibrium behavior within the adiabatic approximation.
In this Thesis we discussed several issues related to the dynamics of isolated quantum many-body systems. The brief overview of Chapter 1 highlighted a series of open fundamental theoretical questions which are crucial for the investigation of nonequilibrium phenomena. Among them, we mainly focused on the emergence of thermal behaviors and on the characterization of dynamical phase transitions and of the associated universal dynamics. To this end, we analyzed some paradigmatic models, such as the Fermi-Hubbard model, the quantum Ising chain, and the O(N) model. -Chapter 2: Absence of thermalization in a Fermi liquid. In this Chapter we report two results that question the common expectation of thermalization in a Landau-Fermi liquid after an interaction quench. We first consider the perturbative expansion in the interaction strength of the long-wavelength structure factor S(q) and show that it does not satisfy the hypothesis that steady-state averages correspond to thermal ones. In particular, S(q) has an analytical component ∼ const. + O(q^2) for q → 0, compatible with thermalization, but it retains also a nonanalytic term ∼ |q| characteristic of a Fermi liquid at zero-temperature. In real space, this nonanalyticity turns in a power-law decay of the density-density correlations, in contrast with the exponential decay associated to thermalization. We next analyze the low-density case, where one can obtain results valid at any order in interaction but at leading in the density, and find that in the steady state the momentum distribution jump at the Fermi surface is strictly finite, though smaller than at equilibrium. -Chapter 3: Nonadiabatic stationary behavior in a driven Ising chain. In this Chapter we discuss the emergence of a nonadiabatic behavior in the dynamics of the order parameter in a low-dimensional system subject to a linear ramp of one of its parameters within a gapped phase, which should be the most favorable situation for an adiabatic evolution. We study in details the dynamics of the spontaneous magnetization m_x(t) in a quantum Ising chain after a linear variation in time of the transverse field within the ordered phase. In particular, focusing on the value of the order parameter at the end of the ramp m_x(τ), we show that the smaller the switching rate of the transverse field is the closer m_x(τ) gets to its ground state value. Nonetheless, this apparently small correction to adiabaticity eventually leads to the disruption of the order exponentially fast in the subsequent time evolution, no matter how slowly the ramp is performed. -Chapter 4: Aging and coarsening in isolated quantum systems after a quench. This Chapter analyzes the nonequilibrium dynamics of an isolated quantum system after a sudden quench to the dynamical critical point, where the emergence of scaling and universal exponents is expected, due to the absence of time and length scales. We explore these features for a bosonic interacting scalar field theory with O(N) symmetry in the large-N limit, where the model is exactly solvable and the exponents and scaling forms of the relevant two-time correlation functions can be analytically calculated. Moreover, we provide evidence of the emergence of a dynamical scaling behavior also for quenches below the dynamical critical point, associated with coarsening. We find that the latter case is characterized by the same scaling functions as those describing the critical case, yet with different exponents. -Chapter 5: Dynamical transitions and statistics of excitations. In this Chapter we study the dynamics of the O(N) model (introduced in the previous Chapter) resulting from a different protocol: a linear ramp of one of its parameters. We find that the presence of a dynamical phase transition, as well as its critical properties, are robust against the change of the protocol. We show that a characterization based on the critical dimensions and exponents would suggest that the dynamical phase transition is analogous to the equilibrium thermal one. However, its nonequilibrium nature becomes evident analyzing the statistics of excitations produced in the ramp process. In particular, the critical properties are encoded in the fluctuations in the number of excitations, which display qualitatively different behaviors depending on the ramp being performed above, at, or below the dynamical critical point. These behaviors bear no dependence on the duration of the protocol.
We study spectral properties and the dynamics after a quench of one-dimensional spinless fermions with short-range interactions and long-range random hopping. We show that a sufficiently fast decay of the hopping term promotes localization effects at finite temperature, which prevents thermalization even if the classical motion is chaotic. For slower decays, we find that thermalization does occur. However, within this model, the latter regime falls in an unexpected universality class, namely, observables exhibit a power-law (as opposed to an exponential) approach to their thermal expectation values.
Motivated by the nonequilibrium dynamics experiment, we study the quantum quench dynamics of the JCH model with weak nearest-neighbor hopping strength by exact diagonalization. One of the crucial insights is to focus on the comparison between the time-averaged density matrix and the canonical ensemble in a wide range of values for the coupling strength between the cavity and the atom, and significant differences are found. Thus, this phenomenon can be used to infer that the system with weak nearest-neighbor hopping strength may be integrable. In order to verify our prediction, we perform the level spacing distribution of the system and find that they conform with the characteristics of the integrable system. Further, a nearly integrable system is presented by removing the degenerate energy level and fitting Brody distribution. Besides, the evolution of the momentum distribution functions of the photon and the equilibrium value predicted by the canonical ensemble are followed, which demonstrates that the system cannot be thermalized. Finally, the way how the system avoids thermalization is also discussed by the finite-size scaling of the fluctuations of the diagonal and the off-diagonal matrix elements.
We present a pseudospin formulation for the post-quench dynamics of charge-density-wave (CDW) order in the half-filled spinless Holstein model on a square lattice, assuming spatially homogeneous evolution. This Anderson pseudospin description captures the coherent nonequilibrium dynamics of the coupled electron-lattice system. Numerical simulations reveal three distinct dynamical regimes of the CDW order parameter following a quench-locked oscillations, Landau-damped dynamics, and overdamped relaxation-closely paralleling quench dynamics in BCS superconductors and other electronically driven symmetry-breaking phases. Crucially, however, the presence of dynamical lattice degrees of freedom leads to qualitatively different long-time behavior. In particular, while the oscillation amplitude is reduced in the damped regimes, CDW oscillations do not fully decay but instead persist indefinitely due to feedback from the lattice field. We further show that these persistent oscillations are characterized by a nonequilibrium electronic distribution, which provides an intuitive understanding of both their amplitude and the renormalization of the oscillation frequency relative to the bare Holstein phonon frequency. Our results highlight the essential role of lattice dynamics in nonequilibrium ordered phases and establish a clear distinction between electron-lattice-driven CDW dynamics and their purely electronic counterparts.
The real-time dynamics of the Fermi-Hubbard model, driven out of equilibrium by quenching or ramping the interaction parameter, is studied within the framework of the nonequilibrium self-energy functional theory. A dynamical impurity approximation with a single auxiliary bath site is considered as a reference system, and the time-dependent hybridization is optimized as prescribed by the variational principle. The dynamical two-site approximation turns out to be useful to study the real-time dynamics on short and intermediate time scales. Depending on the strength of the interaction in the final state, two qualitatively different response regimes are observed. For both weak and strong couplings, qualitative agreement with previous results of nonequilibrium dynamical mean-field theory is found. The two regimes are sharply separated by a critical point at which the low-energy bath degree of freedom decouples in the course of time. We trace the dependence of the critical interaction of the dynamical Mott transition on the duration of the interaction ramp from sudden quenches to adiabatic dynamics and therewith link the dynamical to the equilibrium Mott transition.