We investigate the nonequilibrium behavior of dilute, attractively interacting Fermi gases subjected to finite-duration ramps of their internal interaction strength. We identify three dynamical regimes as a function of ramp duration using a time-dependent version of the Bardeen-Cooper-Schrieffer theory of superconductivity to model these systems. For short ramp durations, these systems become nonsuperconducting; however, fermions with opposite momenta remain paired albeit with reduced amplitudes, and the associated pair amplitude distribution is nonthermal. In this first regime, the disappearance of superconductivity is due to the loss of phase coherence between pairs. By contrast, for intermediate ramp durations, superconductivity survives but the magnitude of the order parameter is reduced and presents long-lived oscillations. Finally, for long ramp durations, phase coherence is almost fully retained during the finite-duration interaction quench, and the steady state is characterized by a thermal-like pair amplitude distribution. Using this approach, one can therefore dynamically tune the coherence between pairs to control the magnitude of the superconducting order parameter and even engineer a nonequilibrium state made of preformed pairs.
We investigate the nonequilibrium behavior of dilute, attractively interacting Fermi gases subjected to finite-duration ramps of their internal interaction strength. We identify three dynamical regimes as a function of ramp duration using a time-dependent version of the Bardeen-Cooper-Schrieffer theory of superconductivity to model these systems. For short ramp durations, these systems become nonsuperconducting; however, fermions with opposite momenta remain paired albeit with reduced amplitudes, and the associated pair amplitude distribution is nonthermal. In this first regime, the disappearance of superconductivity is due to the loss of phase coherence between pairs. By contrast, for intermediate ramp durations, superconductivity survives but the magnitude of the order parameter is reduced and presents long-lived oscillations. Finally, for long ramp durations, phase coherence is almost fully retained during the finite-duration interaction quench, and the steady state is characterized by a thermal-like pair amplitude distribution. Using this approach, one can therefore dynamically tune the coherence between pairs to control the magnitude of the superconducting order parameter and even engineer a nonequilibrium state made of preformed pairs.
Using the time-dependent density matrix renormalization-group method and exact diagonalization, we study the nonequilibrium dynamics of the one-dimensional Fermi-Hubbard model following a quantum quench or a ramp of the on-site interaction strength. We are particularly interested in the nonequilibrium evolution of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) correlations, which, at finite spin polarizations and for attractive interactions, are the dominant two-body correlations in the ground state. For quenches from the noninteracting to the attractive regime, we investigate the dynamical emergence of FFLO correlations and their signatures in the pair quasi--momentum distribution function. We observe that the postquench double occupancy exhibits a maximum as the interaction strength becomes of the order of the bandwidth. Finally, we study quenches and ramps from attractive to repulsive interactions, which imprint FFLO correlations onto repulsively bound pairs. We show that a quite short ramp time is sufficient to wipe out the characteristic FFLO features in the postquench pair momentum distribution functions.
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.
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.
We study the crossover between the sudden quench limit and the adiabatic dynamics of superconducting states in the attractive Hubbard model. We focus on the dynamics induced by the change of the attractive interaction during a finite ramp time. The ramp time is varied in order to track the evolution of the dynamical phase diagram from the sudden quench to the equilibrium limit. Two different dynamical regimes are realized for quenches towards weak and strong coupling interactions. At weak coupling the dynamics depends only on the energy injected into the system, whereas a dynamics retaining memory of the initial state takes place at strong coupling. We show that this is related to a sharp transition between a weak and a strong coupling quench dynamical regime, which defines the boundaries beyond which a dynamics independent from the initial state is recovered. Comparing the dynamics in the superconducting and nonsuperconducting phases, we argue that this is due to the lack of an adiabatic connection to the equilibrium ground state for nonequilibrium superconducting states in the strong coupling quench regime.
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.