Response of superfluid fermions at finite temperature

A consistent finite-temperature microscopic theory for the response of strongly coupled superfluid fermionic systems is formulated. We start from the general many-body Hamiltonian with the vacuum (bare) two-fermion interaction and derive the equation of motion (EOM) for the thermally averaged two-time two-fermion correlation function, which determines the spectrum of the system under study. The superfluidity is introduced via the Bogoliubov transformation of the fermionic field operators, and the entire formalism is carried out in the basis of Bogoliubov's quasiparticles, keeping the complete $4\times 4$ block matrix structure of the two-fermion EOM. Static and dynamical interaction kernels of the resulting EOM are discussed. The least correlated static kernel reproduces the known finite-temperature quasiparticle random phase approximation, while the dynamical kernel was previously largely unknown in the finite-temperature superfluid theory. The major focus is then placed on the latter kernel, which is completed and advanced to a factorized form, enabling a minimal truncation of the many-body problem with keeping important effects of emergent collectivity and mapping to the quasiparticle-vibration coupling (qPVC). As in the zero-temperature and non-superfluid cases, the qPVC can be associated with a new order parameter, qPVC vertex. In the thermal superfluid theory, the latter vertex as well as the components of the dynamical kernel acquire an extended form including thermally unblocked transition amplitudes.