Quench dynamics of fermionic superfluids are an active topic both experimentally and theoretically. Using the BCS theory, such non-equilibrium problems can be reduced to nearly independent spin dynamics, only with a time-dependent mean-field pairing term. This results in persisting oscillations of the pairing strength in certain parameter regimes. However, experiments have observed that the oscillations decay rapidly when the interaction becomes strong, such as in the unitary Fermi gas [Phys. Rev. Res. 3, 023205 (2021)]. Theoretical analysis on this matter is still absent. In this work, we construct an SYK-like model to analyze the effect of strong interactions in a one-dimensional BCS system. We employ the large- N N approximation and a Green’s function-based technique to solve the equilibrium problem and quench dynamics. Our findings reveal that a strong SYK interaction suppresses the pairing order. Additionally, we verify that the system quickly thermalizes with SYK interactions, whether it involves intrinsic pairing order or proximity effect, resulting in a rapid decay of the oscillation strength. The decay rates exhibit different scaling laws against SYK interaction, which can be understood in terms of the Boltzmann equation. This work represents a first step towards understanding the attenuating dynamics of strongly interacting fermionic superfluids.
Quench dynamics of fermionic superfluids are an active topic both experimentally and theoretically. Using the BCS theory, such non-equilibrium problems can be reduced to nearly independent spin dynamics, only with a time-dependent mean-field pairing term. This results in persisting oscillations of the pairing strength in certain parameter regimes. However, experiments have observed that the oscillations decay rapidly when the interaction becomes strong, such as in the unitary Fermi gas [Phys. Rev. Res. 3, 023205 (2021)]. Theoretical analysis on this matter is still absent. In this work, we construct an SYK-like model to analyze the effect of strong interactions in a one-dimensional BCS system. We employ the large- N N approximation and a Green’s function-based technique to solve the equilibrium problem and quench dynamics. Our findings reveal that a strong SYK interaction suppresses the pairing order. Additionally, we verify that the system quickly thermalizes with SYK interactions, whether it involves intrinsic pairing order or proximity effect, resulting in a rapid decay of the oscillation strength. The decay rates exhibit different scaling laws against SYK interaction, which can be understood in terms of the Boltzmann equation. This work represents a first step towards understanding the attenuating dynamics of strongly interacting fermionic superfluids.
Non-equilibrium phenomena in many-body quantum systems involves the interplay of correlations and dynamics spanning different spatio-temporal scales. We present an experimental study of these dynamics in a two-component Fermi gas following a quench from weak to strong interactions. Starting with a weakly attractive gas in the normal phase, interactions are ramped to the unitarity limit at a range of rates, which involves crossing the normal to superfluid phase transition. We measure the subsequent dynamics of the pair momentum distribution which reveals the timescales for both the formation and condensation of pairs. Furthermore, we show that the contact parameter responds very rapidly to changes in the interaction strength, demonstrating that short-range correlations, based on the occupation of high-momentum modes, evolve far more rapidly than correlations in low-momentum modes necessary for pair condensation.
A new field of research, which gained interest in the past years, is the field of non- equilibrium physics of strongly interacting fermionic systems. In this thesis we present a novel apparatus to study an ultracold strongly interacting superfluid Fermi gas driven out of equilibrium. In more detail, we study the excitation of a collective mode, the Higgs mode, and the dynamics occurring after a rapid quench of the interaction strength. Ultracold gases are ideal candidates to explore the physics of strongly interacting Fermi gases due to their purity and the possibility to continuously change many different system parameters, such as the interaction strength. The novel apparatus enabling the experiments outlined below is built from scratch. The setup is described and characterized in detail. We produce a fermionic superfluid of 4 × 10^6 6Li atoms at a temperature of T/T_F = 0.07 ± 0.02, where T_F is the Fermi temperature. This is achieved within 23 s by a combination of laser cooling, radio frequency evaporation of 23Na and simultaneous sympathetic cooling of 6Li in an optical plugged magnetic trap and subsequent evaporative cooling of 6Li in a dipole trap. In this thesis we provide the first experimental evidence of the Higgs mode in a strongly interacting Fermi gas in the crossover from a BCS superfluid to a molecular Bose-Einstein condensate. We develop a novel excitation method, which directly couples to the amplitude of the order parameter. This is achieved by continuously changing the population of one spin component by driving a radio frequency transition to a previously unoccupied third hyperfine state. This effectively modulates the interaction strength and the amplitude of the order parameter. We spectroscopically observe a resonance behavior at twice the gap frequency. For strong coupling, the peak width broadens and eventually the Higgs mode disappears when the Cooper pairs turn into tightly bound dimers signaling the instability of the Higgs mode. It has been suggested that the Higgs mode frequency is a precise measure of the superconducting gap in the BEC-BCS crossover, where the exact value of the gap is yet unknown and numerical calculations are challenging. Hence, our novel method provides a unique technique to determine the superconducting gap. Moreover, we perform rapid variations of the interaction parameter and vary both the initial interaction strength and the amplitude of the quench. The rapid quenches are performed by rapidly transferring one spin component of a two-component spin mixture into a third state with a different interaction strength by using a radio frequency transition. Using this novel method, we can perform the quenches in half the Fermi time, which is experimentally very close to a sudden change and to date faster than any other research group performing interaction quenches with ultra cold gases. In the experiment, we observe a fast relaxation to a zero order parameter for large quenches, whereas for small quenches we observe, after a sudden drop, a revival of the order parameter and equilibration to a long-term superfluid steady state. Our measurement provides the first evidence of the collapse and subsequent revival of order in a strongly interacting fermionic system.
We revisit the study of amplitude oscillations in a pair condensate of fermions after an interaction quench, and generalize it to nonzero temperature. For small variations of the order parameter, we show that the energy transfer during the quench determines both the asymptotic pseudo-equilibrated value of the order parameter and the magnitude of the oscillations, after multiplication by, respectively, the static response of the order parameter and spectral weight of the pair-breaking threshold. Since the energy transferred to the condensed pairs decreases with temperature as the superfluid contact, the oscillations eventually disappear at the critical temperature. For deeper quenches, we generalize the regimes of persistent oscillations and monotonic decay to nonzero temperatures, and explain how they become more abrupt and are more easily entered at high temperatures when the ratio of the initial to final gap either diverges, when quenching towards the normal phase, or tends to zero, when quenching towards the superfluid phase. Our results are directly relevant for existing and future experiments on the non-equilibrium evolution of Fermi superfluids near the phase transition.
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 theoretically investigate a quartet superfluid state in fermionic matter by using the quartet Bardeen-Cooper-Schrieffer (BCS) variational theory and the Green's function method. We demonstrate that the quartet BCS theory with the multiple-infinite-product ansatz successfully reproduces an exact four-body result in a one-dimensional four-component Fermi gas at the dilute limit, in contrast to the single-infinite-product ansatz. To see the validity of the quartet BCS state, we derive the self-consistent equation for the quartet superfluid order parameter within the generalized imaginary-time Nambu-Gor'kov formalism, which is found to be consistent with the quartet BCS variational equation. Moreover, by numerically computing the momentum-resolved single-particle spectral function in a one-dimensional system, we discuss how the single-particle spectra evolve with increasing the strength of the four-body cluster formation. We show that a coherent BCS-like quasiparticle branch on the weak-coupling side evolves into a strongly damped, continuum-dominated spectrum in the strong-coupling side, while nonzero quartet superfluid order parameter persists throughout the crossover regime. Our results would be useful for understanding beyond-BCS pairing effects and four-body cluster formations in fermionic systems in an interdisciplinary way.
We study the response of a Fermi superfluid to a rapid change of the interaction strength. For a broad range of quench parameters, the order parameter exhibits damped oscillations, however, we also identify quench regimes in which these oscillations are absent. By comparing our data to a numerical model we find that the damping time constants are explainable by a dephasing in a local density approximation, however, the oscillation frequencies, while being comparable to the superconducting order parameter, display in detail significant differences to the integrable BCS (Bardeen-Cooper-Schrieffer) model.