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.
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.
Motivated by recent experiments in ultracold atomic gases that explore the nonequilibrium dynamics of interacting quantum many-body systems, we investigate the opposite limit of Landau's Fermi-liquid paradigm: We study a Hubbard model with a sudden interaction quench, that is, the interaction is switched on at time t=0. Using the flow equation method, we are able to study the real time dynamics for weak interaction U in a systematic expansion and find three clearly separated time regimes: (i) An initial buildup of correlations where the quasiparticles are formed. (ii) An intermediate quasi-steady regime resembling a zero temperature Fermi liquid with a nonequilibrium quasiparticle distribution function. (iii) The long-time limit described by a quantum Boltzmann equation leading to thermalization of the momentum distribution function with a temperature T proportional, variantU.
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.
In the last years the impressive progress on the experimental side led to a variety of new experiments allowing to address systems out of equilibrium.In this way the behavior of such systems far from equilibrium is no longer a purely theoretical issue but indeed observable.New experimental techniques, like particles trapped in optical lattices, render a realization of quantum systems with nearly arbitrary system parameters possible and provide a possibility to study their time evolution.Systems out of equilibrium are characterized by the fact, that these systems are in highly excited states giving rise to totally new fascinating properties.In the present thesis one-and two-dimensional fermionic Hubbard models out of equilibrium are discussed.The system is taken out of equilibrium by a so-called interaction quench.At the beginning the system is prepared in the groundstate of the non-interacting Hamiltonian.At a time t the interaction between the fermions is suddenly turned on so that the time evolution is governed by the whole, interacting Hamiltonian.Hence the system is prepared in the groundstate of one Hamiltonian but evolves according to a di_erent Hamiltonian.Consequently the system ends up in a highly excited state.To describe such a system a method based on an expansion of the Heisenberg equations of motion to highest order possible is developed in this thesis.This method provides an exact description of the time evolution on short and intermediate time scales after the quench.As the method reveal exact results and does not rely on any perturbative assumption, a study of arbitrarily large interaction strengths is possible.Besides, the method is one of the few methods capable of two-dimensional systems.In the following the method used in this thesis is explained and advantages and disadvantages of the approach are thematized.For this purpose the results of the developed iterated equation of motion approach are compared to results obtained in calculation up to second order in the interaction.Then the dynamics of the onedimensional model is discussed with a focus on the relation of the results derived by the iterated equations of motion approach to results obtained by bosonization theory, the behavior for strong interactions and the dynamical transition from the weak to the strong quench regime.Furthermore, the Hubbard model is studied on a two-dimensional square lattice.This model is fundamentally di_erent from the one-dimensional model: In contrast to the one-dimensional model the two-dimensional model is not integrable allowing a true relaxation of the system.For this system a calculation up to second order in the interaction is performed and compared to the results of the iterated equation of motion approach.Besides, the time evolution of the momentum distribution and the influence of doping on the dynamics is studied.Moreover, a first estimate for relaxation times is provided without relying on the assumption of a mixed state.
Strongly correlated quantum systems display a myriad of complex, fascinating behaviour. The competition between charge, spin, orbital, or lattice degrees of freedom and their respective energy scales gives rise to rich and complex emergent phases. More recently, the non-equilibrium properties of correlated quantum many-body systems have come into the focus of active research, with remarkable advances. While the advent of tensor network based simulations on the theoretical side has allowed access to the full time-evolution of the quantum state as it traverses the Hilbert space, the experimental progress in the ultracold atom community has allowed physicists to enter the paradigm of analog quantum simulation. Their tunability and isolation from the environment makes them ideal platforms to emulate and study open problems in condensed matter physics and go beyond the simulational capabilities of state of the art classical simulations. In this thesis, we explore the dynamics of correlated Fermi gases. We develop and present different facets of applicability of radiofrequency modulation techniques to drive the system out of equilibrium. We begin by considering a homogeneous, three-dimensional Fermi gas in the BCS-BEC crossover and develop an off-resonant radiofrequency transfer scheme to excite the Higgs mode of the superfluid. In a subsequent study we investigate the tunability of the coherence between Cooper pairs by tuning the duration of interaction ramps of the internal interaction strength. These works highlight the possibility to externally tune and stabilise complex quantum many-body states away from equilibrium. Using quasi-exact matrix product state simulations, we study the response of an attractively interacting, one-dimensional Fermi-Hubbard model to weak radiofrequency perturbations. We reveal the emergence of two distinct dynamical regimes in the time evolution and are able to relate it back to the underlying excitation spectrum of the system. Finally, we explore the possibility to perform quantum quenches through radiofrequency π-pulses, achieving nearly complete population transfer between different interacting fermionic states. We reveal non-trivial dynamical effects in the pair correlation of the final state, as well as the excitation of a collective trap mode of the system.
In this thesis, strongly correlated quantum many-body systems in equilibrium and in out-of-equilibrium situations are investigated. This is done by applying and developing well established numerical methods. The focus of the thesis lies in the development and application of the density matrix renormalization group method (DMRG) to quantum many-body systems out of equilbrium. In this thesis, in addition to the DMRG, methods for the exact diagonalization of the Hamiltonian of the system, like the Lanczos- or the Jacobi-Davidson method, are treated. An extension of the Lanczos method makes it possible to treat the time evolution of strongly correlated quantum systems with an accuracy comparable to machine precision. This method is the basis for a possible extension of the DMRG for the treatment of systems out of equilibrium, the so-called adaptive time-dependent DMRG ("adaptive t-DMRG"). A second variant of the adaptive t-DMRG uses the Suzuki-Trotter decomposition of the time-evolution operator. An error analysis demonstrates that both methods, for a suitable choice of control parameters, have errors < 1% at the end of the time evolution in complicated quantities like, e.g., the momentum distribution function. We apply these numerical methods to investigate the quantum critical behavior of a variant of the Hubbard model and to treat two non-equilibrium situations of current interest. Extensive use of the DMRG makes it possible to clarify the quantum-critical behavior of the ionic Hubbard model; in particular, the numerical results demonstrate quite clearly that the correct scenario has two critical points; at the first critical point, only the charge degrees of freedom are critical, while at the second one only the spin degrees of freedom are critical. Quite surprisingly, the Mott-insulator phase in the strong-coupling regime shows a divergent electrical susceptibility. Next, the dynamics of a system of so-called soft-core bosons is treated. The particles initially are trapped in a deep box-potential and the dynamics is investigated after releasing them from this trapping potential. Similar to exact results obtained for so-called hard-core bosons, (quasi-)coherent matter waves emerge, demonstrating the possibility of realizing an atom laser in experiments on optical lattices. We find that the wave vector of the emerging matter wave can be tuned by changing the strength of the interaction between the particles. The second non-equilibrium situation is a so-called "quantum-quench" of a system of strongly correlated fermions, i.e., the evolution of the system after suddenly changing an intrinsic parameter like, e.g., the interaction strength between the particles, is investigated. In this case, we are mainly interested in the long-time behavior; in particular, the question arises whether, relying on Boltzmann's ergodic hypothesis, relaxation to a thermal state is obtained. We find that, in general, non-thermal final states are reached, which can be described by generalized Gibbs-Boltzmann ensembles.
The recent experimental implementation of condensed matter models in optical lattices has motivated research on their non-equilibrium behavior. Predictions about the dynamics of superconductors following a sudden quench of the pairing interaction have been made based on the effective Bardeen–Cooper–Schrieffer (BCS) Hamiltonian; however, their experimental verification requires the preparation of a suitable excited state of the Hubbard model along a twofold constraint: (i) a sufficiently non-adiabatic ramping scheme is essential to excite the non-equilibrium dynamics and (ii) overheating beyond the critical temperature of superconductivity must be avoided. For commonly discussed interaction ramps, there is no clear separation of the corresponding energy scales. Here we show that the matching of both conditions is simplified by the intrinsic relaxation behavior of ultracold fermionic systems: for the particular example of a linear ramp we examine the transient regime of prethermalization (Moeckel and Kehrein 2008 Phys. Rev. Lett. 100 175702) under the crossover from sudden to adiabatic switching using Keldysh perturbation theory. A real-time analysis of the momentum distribution exhibits a temporal separation of an early energy relaxation and its later thermalization by scattering events. For long but finite ramping times this separation can be large. In the prethermalization regime, the momentum distribution resembles a zero-temperature Fermi liquid as the energy inserted by the ramp remains located in high-energy modes. Thus ultracold fermions prove to be robust to heating, which simplifies the observation of non-equilibrium BCS dynamics in optical lattices.