ReceivedHere we present a quantum-statistical formulation of relaxation dynamics of superconductors related to pump-probe (PP) spectroscopy. The method is based on the perturbation expansion of the non-equilibrium density matrix for calculation of multi-time correlation functions, and the corresponding response function. As a model for our study, the high-temperature superconductor MgB 2 has been selected. Knowledge of the electronic structure and of the corresponding Eliashberg function of electron-phonon (EP) interactions enabled us to distinguish nonequilibrium processes in a sequence of interactions with laser pulses on a microscopic level. Temperature-dependent relaxation dynamics as a function of delay time between the pump and probe pulses have been derived. We have shown that an abrupt increase of the relaxation time at T c is the direct consequence of sudden changes in the character of EP coupling in transition from an adiabatic to a stabilized superconducting anti-adiabatic state, as it predicts the anti-adiabatic theory of electron-vibration interactions. The BCS model, which is based on adiabatic character of EP coupling, is basically unable to reflect the enormous sudden increase of the relaxation time at T c . Differences in the optical pump-optical probe and the optical pump-terahertz probe settings of PP relaxation dynamics are discussed within diagrammatic perturbation theory.Keywords: pump-probe spectroscopy; relaxation dynamics of superconductors; electron-phonon coupling; theory of superconductivity; anti-adiabatic theory of superconductivity
IntroductionPump-probe (PP) spectroscopy, a method of nonlinear optics, is an important experimental tool in the study of ultrafast relaxation dynamics of excited states in condensed matter systems [1,2]. The interpretation of experimental results relies on various types of phenomenological equations of motion for the relaxation of "hot electrons" and/or high-energy phonons in coupling to bosonic degrees of freedom [3][4][5][6][7]. The assumption of quasi-equilibrium, which is often behind these methods, is discussed in [8], where the authors present an analytic solution of excited state relaxation dynamics within the Boltzmann kinetic equation, the method of which deals with electrons and phonons out of equilibrium.In general, however, phenomenological treatments not only ignore nonlinear aspects of processes running in PP experiments, but ignore also the microscopic character of nonequilibrium states -coherences and/or populations, which can evolve in the studied system over a sequence of different time periods via the interaction of matter with the incident fields. In these circumstances, mainly in the case of superconductors, it is difficult to disclose the relationship between the derived relaxation parameters, i.e. relaxation times, or electron-phonon (EP) coupling constants (), and the microscopic mechanism of the superconducting state transition. The theory of high-temperature superconductors (e.g. cuprates, pnictides) is far from being settled,...