Nonequilibrium dynamical mean field theory and nonequilibrium self-consistent strong coupling expansion[1] are used to study the relaxation of correlated quantum systems driven out of equilibrium by DC electric fields[2]. Both the FalicovKimball and the Hubbard model are found to exhibit regimes of monotonic or oscillatory themalization as well as regimes where they evolve in a monotonic or oscillatory manner towards a non-thermal state. This suggests that driven quantum systems have a richer behavior than their quenched counterparts and that integrability does not play as critical a role. In the monotonic thermalization scenario, the system evolves through successive quasi-thermal states and it is possible to extrapolate its long time properties from its transient; bridging the gap between the transient and the steady state with very little computational cost. Furthermore, regardless of the relaxation scenario, it is interesting to ask how the particles are distributed as the system evolves in time. We will show that non-trivial parameter-dependent patterns are formed when the system is visualized in momentum space [3]. These features should be observable in current cold atom experiments.
Nonequilibrium dynamical mean field theory and nonequilibrium self-consistent strong coupling expansion[1] are used to study the relaxation of correlated quantum systems driven out of equilibrium by DC electric fields[2]. Both the FalicovKimball and the Hubbard model are found to exhibit regimes of monotonic or oscillatory themalization as well as regimes where they evolve in a monotonic or oscillatory manner towards a non-thermal state. This suggests that driven quantum systems have a richer behavior than their quenched counterparts and that integrability does not play as critical a role. In the monotonic thermalization scenario, the system evolves through successive quasi-thermal states and it is possible to extrapolate its long time properties from its transient; bridging the gap between the transient and the steady state with very little computational cost. Furthermore, regardless of the relaxation scenario, it is interesting to ask how the particles are distributed as the system evolves in time. We will show that non-trivial parameter-dependent patterns are formed when the system is visualized in momentum space [3]. These features should be observable in current cold atom experiments.
There is much interest in how quantum systems thermalize after a sudden change, because unitary evolution should preclude thermalization. The eigenstate thermalization hypothesis resolves this because all observables for quantum states in a small energy window have essentially the same value; it is violated for integrable systems due to the infinite number of conserved quantities. Here, we show that when a system is driven by a DC electric field there are five generic behaviors: (i) monotonic or (ii) oscillatory approach to an infinite-temperature steady state; (iii) monotonic or (iv) oscillatory approach to a nonthermal steady state; or (v) evolution to an oscillatory state. Examining the Hubbard model (which thermalizes under a quench) and the Falicov-Kimball model (which does not), we find both exhibit scenarios (i-iv), while only Hubbard shows scenario (v). This shows richer behavior than in interaction quenches and integrability in the absence of a field plays no role.
th and O Sts. NW, Washington, DC 20057 USA. There is much interest in how quantum systems thermalize after a sudden change, because unitary evolution should preclude thermalization. The eigenstate thermalization hypothesis resolves this because all observables for quantum states in a small energy window have essentially the same value; it is violated for integrable systems due to the infinite number of conserved quantities. Here, we show that when a system is driven by a DC electric field there are five generic behaviors: (i) monotonic or (ii) oscillatory approach to an infinite-temperature steady state; (iii) monotonic or (iv) oscillatory approach to a nonthermal steady state; or (v) evolution to an oscillatory state. Examining the Hubbard model (which thermalizes under a quench) and the Falicov-Kimball model (which does not), we find both exhibit scenarios (i–iv), while only Hubbard shows scenario (v). This shows richer behavior than in interaction quenches and integrability in the absence of a field plays no role.
Recent experimental advances in ultrafast phenomena have triggered renewed interest in the dynamics of correlated quantum systems away from equilibrium. We review nonequilibrium dynamical mean-field theory studies of both the transient and the steady states of a DC field-driven correlated quantum system. In particular, we focus on the nonequilibrium behavior and how it relates to the fluctuation-dissipation theorem. The fluctuation-dissipation theorem emerges as an indicator for how the system thermalizes and for how it reaches a steady state. When the system thermalizes in an infinite temperature steady state it can pass through a succession of quasi-thermal states that approximately obey the fluctuation-dissipation theorem. We also discuss the Wigner distribution and what its evolution tells us about the nonequilibrium many-body problem.
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.
Many-body quantum systems in nonequilibrium remain one of the frontiers of many-body physics. While there has been significant advances in describing the short-time evolution of these systems using a variety of different numerical algorithms, it has been quite difficult to evolve a system from an equilibrium state prior to the application of a driving field, to the long-time steady (or periodically oscillating) state. These dynamics are complex: the retarded quantities tend to approach their long-time limit much faster than the lesser (or greater) quantities. Recent work on strongly correlated electrons in DC electric fields illustrated that the system may evolve through successive quasi-thermal states obeying an effective fluctuation-dissipation theorem in time. We demonstrate an extrapolation scheme that uses the short-time transient calculation to obtain the retarded quantities and to extract how the lesser/greater quantities vary with time and then extend the numerical solutions all the way to the steady state, with minimal additional computational cost. Our approach focuses on extrapolating the electronic self-energy and then employing that to determine the Green's function and various experimentally relevant expectation values.
The study of the interplay between interaction-induced correlations and nonequilibrium initial conditions in many-body systems has recently attracted a lot of attention. New experimental techniques provide high control over many-body systems and systematic access to their nonequilibrium regime. In many systems the nonequilibrium properties of a Fermi liquid can be relevant. A first approach to their understanding is the main content of this thesis. At the beginning, nonequilibrium phenomena, Landau's theory of a Fermi liquid, the Hubbard model, time evolution in quantum mechanics, related experiments, and the flow equation method are reviewed. The key observation of this thesis, namely a characteristic mismatch of expectation values in equilibrium and nonequilibrium, is first illustrated for the squeezed oscillator. For this one-particle model Hamiltonian the perturbative approach can be compared with the exact solution and it is shown that the mismatch holds even beyond perturbation theory. Afterwards these observations are generalized to a larger class of one-particle models. Then the nonequilibrium behavior of a Fermi liquid is examined by analyzing the Fermi liquid phase of the Hubbard model in more than one dimension. After a sudden switch-on of a weak two-particle interaction to the noninteracting Fermi gas the relaxation of the many-body system is observed. For this purpose, the flow equation transformation is implemented for the Hubbard Hamiltonian. Then the discussion of the momentum distribution function and of the kinetic energy displays a three-step relaxation behavior of the Fermi liquid from the initial perturbation until thermalization is reached. Firstly, the sudden switch inserts excitation energy into the system which drives the following dynamics. Then a rapid initial build-up of correlations is caused by initial dephasing and leads to the establishment of a quasiparticle picture. By that the system enters into a quasi-stationary state which can be long-lasting for weak interaction. This state shows prethermalization of the kinetic energy which already has relaxed to its final values. However, the momentum distribution still resembles a zero temperature Fermi liquid. Its later relaxation on a second time scale is caused by a residual two-particle interaction which allows for scattering processes and can be described by a quantum Boltzmann equation. The physical origin of this delayed relaxation can be traced back to an interplay of translational invariance of the Hubbard model and the Pauli principle for fermions; together they restrict the phase space for scattering events. Comparing with similar work I conjecture on the generic nature of the findings made for the quenched Fermi liquid in other many-body systems. Finally I point out to the potential relevance of the delayed relaxation for the observation of further nonequilibrium phenomena, for instance in BCS systems. In order to extend the study of sudden switching to arbitrary switching processes the calculation is repeated using the Keldysh perturbation theory. First evaluations for a linear increase of the interaction strength are given. Extensions of this work are suggested and may motivate further research.