Attenuating dynamics of strongly interacting fermionic superfluids in SYK solvable models

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.