Non-equilibrium Phenomena in Superconductors Probed by Femtosecond Time-Domain Spectroscopy

Development of ultrafast lasers and nonlinear optical techniques over the last two decades provides tools to access real-time dynamics of low energy excitations in superconductors. For example, time-resolved THz spectroscopy and time- and angular-resolved photoemission spectroscopy provide access to the real-time dynamics of the superconducting gap amplitude. Such studies enable determination of microscopic parameters like quasi-particle recombination rates, pair-breaking rates and electron–boson coupling constants. Recently, intense THz pulses have been used to probe the nonlinear dynamics, including observation of collective modes. Moreover, using low-frequency electromagnetic pulses, there are several reports of amplification of superconductivity in both conventional and unconventional superconductors. Starting with a brief historical overview of the pioneering work, where non-equilibrium phenomena in superconductors were investigated using quasi-continuous excitation, we review some of the insights that are provided by using real-time approaches. We focus on conventional BCS superconductors, whose ground state is reasonably well understood, and address similarities and open questions related to the corresponding studies in high-Tc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${T}_{c}$$\end{document} superconductors.


Introduction
The field of non-equilibrium superconductivity [1,2] started soon after the first observation of the superconducting gap using optical spectroscopy [3] and the development of the BCS theory. Given that small energy gap Δ on the order of a meV (depending on the material), superconductors were first considered as potential detectors for far-infrared radiation [4]. For such applications, understanding the Journal of Low Temperature Physics (2020) 201:  dynamical processes in a superconductor as a function of base temperature, photon energy and excitation density is crucial. On the other hand, even more fascinating observations followed in early 1960s, where superconductivity in Al and Sn was enhanced (increasing the critical current, the superconducting gap and even the critical temperature) under illumination with electromagnetic radiation at sub-gap frequencies [5][6][7][8]. Most of the early studies were performed in a quasi-continuous excitation regime, with the characteristic timescales of the processes only estimated. In the last two decades or so, however, the time-resolved techniques, with their ability to probe the evolution of optical constants, electronic band structure, magnetization and crystal structure with femtosecond time resolution contributed to further our understanding of out-of-equilibrium superconductors.
As noted, non-equilibrium superconductivity has been an intriguing field of the modern solid state physics since the 1960s [1]. Here, we first briefly review the results of the fascinating studies performed in 1960s and 1970s-for more details we refer the reader to Refs. [1,2]. Then, we turn to studies of real-time dynamics, driven by the rapid development of femtosecond lasers and the related nonlinear optical methods. While the first real-time experiments with picosecond time resolution were performed on conventional superconductor Pb [9], the field started to blossom following the early studies on cuprate superconductors [10][11][12][13][14][15][16][17][18]. The first systematic real-time studies of dynamics in conventional superconductors were performed substantially later [19][20][21][22]. Such studies still represent a small fraction of the body of work performed on superconductors, still largely dominated by studies of cuprate and pnictide high-temperature superconductors. Nevertheless, we focus on the results obtained on standard s-wave superconductors and address the differences and peculiarities of the results obtained on unconventional high-temperature superconductors after that. For a review, focusing on ultrafast spectroscopy in correlated high-temperature superconductors, we refer the interested reader to a recent review [23].

Dynamics Following the Above Gap Excitation of a Superconductor
We should first note that, in general, there are two types of non-equilibrium that can be realized in superconductors [1], often referred to as the energy-mode and chargemode imbalance. The two types are related to the nature of excitations in superconductors, which continuously changes from a hole-like (far below the Fermi level) to an electron-like (far above the Fermi level) within the energy range of ≈ 2Δ around the Fermi level.
The first type of nonequilibrium is realized by perturation with charged particles, e.g. by quasiparticle injection through a superconductor-metal junction, resulting in net quasiparticle charge density (note the fractional charge of quasiparticles close to Δ ). In this case, the populations of the electron-like and the hole-like quasiparticle branches do not coincide, and such a state is referred to as charge imbalance (or branch-imbalance). Correspondingly, there is a difference in the chemical potentials of quasiparticles and the condensate (see Figure 11.1 of Ref. [1]).

3
The energy-mode imbalance is realized by exciting a superconductor with charge-neutral drive, like electromagnetic radiation or ultrasound, or with symmetrical tunnel injection in a superconductor-insulator-superconductor junction. Here, the nonequilibrium state is characterized by a (non-equilibrium) distribution function, which is particle-hole symmetric; i.e. the populations of the electron-like and the hole-like branches of the quasiparticle spectrum are the same (see Figure 11.1 of Ref. [1]). Often, this state can be approximated by an effective quasiparticle temperature T * . It is the energy-mode imbalance that is considered in the case of photoexcitation discussed in this work.
The first non-equilibrium studies/effects were observed using injection of quasiparticles through tunnel barriers [24], or by (quasi-)continuous laser irradiation [25]. Here, photo-excitation with light, whose frequency is larger than 2Δ , breaks Cooper pairs and gives rise to excess quasiparticle (QP) density.
Experiments in multi-tunnel-junction configuration provided first estimates of lifetimes of excited QPs. While radiative recombination lifetime was theoretically estimated for Pb to ≈ 0.4 s at 2 K [4], tunnel-junction experiments provided the upper limit for the average QP lifetime in Pb to be of the order of 0.2 μs at 1.44 K [24]. It was argued that the relaxation proceeds via QP recombination, where QP pairs recombine to form Cooper pairs by the emission of phonons with � > 2Δ [24,26]. We note again that these experiments were performed under the constant drive, with the relaxation time estimates based on a series of underlying assumptions [24,26]. One of the first systematic temperature-dependent studies using QP injection was performed in Al [27]. The authors noted that the errors of the extracted relaxation times are relative only (due to the fact that the extracted absolute values of relaxation time depend on the choice/uncertainty of the model parameters used) [27,28]. Despite the uncertainty in the absolute values, the results revealed that at temperatures far below T c the recombination rate followed an exponential temperature dependence with −1 R ∝ exp(−0.3Δ∕k B T) [27]. This behavior (apart from the factor 0.3 in the exponent) could be accounted for by considering that the relaxation of the QP density, n QP , reflects the bi-particle recombination process (see Fig. 1a).
In this case, the dynamics of excess QP's is given by the bi-molecular rate equation where the recombination rate R is a microscopic parameter. These studies were performed in the limit of weak perturbations, where the density of injected (excess) QPs, n ′ QP , is small compared to the density of thermally excited QPs, n T QP . With n QP = n T QP + n � QP and n ′ QP ≪ n T QP , the above rate equation can be approximated by dn � QP dt = −2Rn T QP n � QP . It is thus clear that the exponential temperature dependence of relaxation time can be attributed to the temperature dependence of the thermally excited QPs, n T QP ∝ . We should note that most of these initial studies assumed that the re-absorption of � > 2Δ phonons (Fig. 1b) emitted during the QP recombination may be neglected. If these high energy phonons rapidly decay into the substrate, or via anharmonic phonon decay (Fig. 1c), it will be the bi-particle recombination process that will be governing the recovery of the SC state. It was Rothwarf and Taylor [29] that pointed out that in general, phonon re-absorption cannot be neglected (including the cases of high excitation densities). The slow decay of the high frequency phonon populations leads in this case to the so-called phonon bottleneck [29].
One of the first real-time studies of suppression of superconductivity with light pulses was performed by Testardi [30]. He studied changes in resistivity in optically thin lead films upon excitation with s pulses. He demonstrated that the absorbed energy density required to suppress superconductivity is substantially lower than the thermal energy, required to heat up the superconductor up to T c . To account for the observation, Owen and Scalapino proposed the so-called * model [31], where under continuous excitation the QPs, the condensate and the lattice are in thermal equilibrium, yet the chemical potentials of QPs and the condensate are different. This assumption was motivated by considering the thermalization times being faster than the QP recombination times. According to the * model, upon increasing the excitation density, the superconductor should undergo a first order phase transition into the normal state. Considering the idea that it is the escape of high energy phonons that governs the recovery of the superconducting state [29], and following the experimental work of Sai-Halasz at al., who failed to observe the predicted first-order phase transition [32], Parker suggested an alternative model [33]. He argued that in out-of-equilibrium the QPs and the condensate share the same chemical potential. However, their temperature T * , which is the same as the temperature of � > 2Δ phonons, differs from the temperature of low energy ( � < 2Δ ) phonons (base temperature). While it was argued that the above mentioned T * and * models may not adequately describe the behavior in the vicinity of the phase transition [34], the solutions of kinetic equations clearly show non-equilibrium effects in phonon population, consistent with the idea put forward by Parker (see, e.g. Figure 8 of Ref. [34]).

Enhancement of Superconductivity with Sub-gap Excitation
As mentioned above, studies of superconducting properties under continuous irradiation with sub-gap radiation (e.g. microwaves) showed surprising results [5][6][7][8]. These Fig. 2, summarizing observations of microwave enhancement of superconductivity in superconducting thin-film strips of Al [7]. It was demonstrated that at temperatures close to the equilibrium critical temperature and for frequencies of microwave radiation above some characteristic value, one can observe a measurable increase in the superconducting critical current (Fig. 2a), as well as an increase in the superconducting critical temperature and the magnitude of the superconducting gap. In these experiments, the order parameter jumped discontinuously to zero at the superconducting transition (see Fig. 2b) [7,35]. With microwave irradiation, superconductivity could be observed above the equilibrium critical temperature, a phenomenon that is commonly referred to as photo-induced superconductivity.
These observations were explained by Eliashberg [36,37], who realized that in the vicinity of the critical temperature, the effect is mainly related to a non-equilibrium distribution of quasiparticles induced by the microwave field. In particular, at finite temperatures, sub-gap radiation couples to thermally excited QPs at the gap edge, giving rise to excitation of quasiparticles to energies further from the gap edge. Within the BCS theory, the superconducting gap is determined by the self-consistent temperaturedependent gap equation: Fig. 2 a Temperature dependence of the critical current in Al for various levels of microwave power at 3 GHz (attenuations are shown on the plot, with ∞ corresponding to the absence of microwave radiation). The results show not only a microwave enhanced critical current, but also an increase in the superconducting transition temperature. In equilibrium T c = 1.170 K, while under microwave radiation the sample superconducts already at 1.184 K. b The calculated temperature dependence of the superconducting gap in Al under 3 GHz radiation using the Eliashberg theory. Inset presents the resistivity curves with and without microwave radiation (for the latter, the hysteresis demonstrates the first order nature of the phase transition under microwave excitation). The presented data were adapted from Ref. [7] (Color  figure online) where is the kinetic energy relative to the Fermi level. The maximum gap enhancement would result by exciting thermally excited QPs as far as possible from the gap-edge. Using Eq. (2) and expanding for small changes of the order parameter |Δ| 2 = |Δ 0 | 2 (1 + Δ∕|Δ 0 |) as a result of changes in the Fermi-Dirac distribution, where f FD is the Fermi-Dirac distribution and f ( ) is the excitation induced change in the distribution function, one obtains (in the limit of small relative gap changes) [38] Here, C is a positive dimensionless constant that generally depends on the base temperature. From Eq. (3), it follows that pushing the gap-edge quasiparticles to higher energies, f (0) < 0 , results in Δ > 0.
Since these experiments were performed using continuous microwave radiation, while inducing changes in the distribution functions require slow thermalization of QPs, it is no surprise that enhancement was observed only for pump frequencies exceeding the characteristic QP scattering rate. Moreover, the observed enhancements of the gap and critical temperature were quite small (e.g. an increase of T c from 1.170 to 1.184 K in Al-see Fig. 2).
Detailed experimental studies as a function of frequency, power and base temperature showed some departures from the Eliashberg theoretical predictions [35,39]. Chang and Scalapino [40][41][42] proposed an extension of the Eliashberg theory by considering the energy dependence of the QP recombination rate (in Fig. 1a the recombination rate R is considered to be a constant, but in general it depends on QP energy). The high energy quasiparticles recombine into Cooper pairs faster, since there are more available phonon states to relax the energy to. By analyzing the coupled Boltzmann equations for phonon and quasiparticle distributions under sub-gap excitation [42], Chang and Scalapino demonstrated cooling of a superconductor upon photo-excitation.
Finally, not only microwave radiation, but also excitation with sub-gap microwave phonons was demonstrated to give rise to enhancement of superconductivity [43][44][45], which could be explained within the generalized Eliashberg model [46].
For further reading on experimental and theoretical approaches to realize enhancement of superconductivity using quasi-equilibrium sub-gap excitation, we refer to several reviews to [47,48].

Real-Time Spectroscopy of the Superconducting Gap Dynamics
To study the dynamics of light-induced suppression and recovery of the superconducting state in real-time the common approach is to use the so-called pump-probe spectroscopy [49]. In the simplest version, femtosecond optical pulses of the same photon energy are used for excitation (pump) and for probing the induced changes in optical constants (probe) as a function of time delay between the two. Most of the early studies were performed in this configuration [10-15, 18, 50-57]. While the photo-induced changes at optical (probe) frequencies do reflect, via Kramers-Kronig relations, the photo-induced changes and dynamics of the low-energy excitation spectrum (e.g. the superconducting gap), probing the gap dynamics using far-infrared (THz) pulses [9,16,17,[19][20][21][22][58][59][60][61][62] or with photoemission [63][64][65][66][67][68][69], is much more straightforward to interpret. In the following, we first focus on phenomena observed on a conventional BCS superconductor NbN obtained using time-resolved THz spectroscopy [22]. We note that similar results were obtained in another conventional BCS superconductor MgB 2 [20], as well as in an electron-doped cuprate superconductor Pr 2−x Ce x CuO 4 [61]. The latter is especially interesting, given the d-wave symmetry of the order parameter in cuprates, as elaborated in Sect. 2.2.
Before discussing the real-time dynamics of superconductors driven by pulsed laser excitation, we should briefly address the characteristic excitation densities. To do so, let us assume a conventional BCS s-wave superconductor, where most of the absorbed energy of optical pulse is used to break Cooper pairs (as elaborately discussed in the following sections this condition seems to be fulfilled in conventional superconductors, while in the high-T c cuprate superconductors the situation may be more complicated). Under this assumption, the absorbed energy density, A, should be lower than the superconducting condensation energy, E c , if we are to discuss the dynamics in the perturbative regime (for A ≈ E c excitation results in melting of superconductivity). Here, , where F is the fluence of the incoming optical pulse, typically in the μJ/cm 2 range, R is the reflectivity of the sample and l opt is the optical penetration depth. The condensation energy, E c , on the other hand, is-for a conventional BCS s-wave superconductor-given by E c = B 2 c ∕2 0 = N(0)Δ 2 ∕2 , where B c is the thermodynamical critical field and N(0) is the normal state density of states. The above condition, however, does not yet imply the dynamics to be linear, i.e. excitation density independent. Since the relaxation processes contain bimolecular terms (see Eq. 1), the recovery dynamics will be linear only for excitation densities, where the density of photo-excited quasiparticles is small compared to the density of thermally excited quasiparticles, given by n T ≃ N(0) √ 2 Δk B Texp(−Δ∕T) .

Photo-Induced Dynamics of the Superconducting Gap in NbN
Thin NbN thin films on MgO substrates, with T c of 15.4 K were investigated using near-infrared pump THz-probe spectroscopy [70]. Following linear THz spectroscopy studies (Fig. 3a), where the temperature dependence of the superconducting gap, Δ(T) , is extracted using the fit to the complex optical conductivity [71], ( ) , the dynamics of the gap following excitation with a 50 fs near-infrared (800 nm) pulse can be studied. The time evolution of the SC state as a function of time delay, t, after optical excitation, was studied by either directly measuring spectrally resolved ( , t) or by measuring the spectrally integrated response by tracking the time evolution of the induced changes in the transmitted electric field, ΔE tr (t � , t) , at a fixed t ′ (see Fig. 3c). The latter approach is particularly useful for studying dynamics at low excitations. Given the fact that films are optically thin, they are homogeneously excited. The absorbed energy density, A (in units meV per unit cell volume, meV/ucv), can be determined from the reflectivity and transmission at 800 nm (pump).
There are several interesting observations in Fig. 3 that should be addressed: (a) The time-evolution of ( ) shows that even at quite high excitation densities as presented here, ( , t) can be well described by the time-evolution of the superconducting gap, Δ(t) . In fact, the data suggest that ( , t) can be approximated by ( , T * (t)) , where at any given time delay (apart from perhaps on the sub-ps timescale) the superconducting state can be approximated by an effective temperature T * , similar to the early description by Parker [33].
(b) Despite the fact that the sample is excited by a 50 fs optical pulse, the pairbreaking process (initial decrease of Δ ) proceeds on a long 10 ps timescale. The observation is surprising, since the e-e thermalization time in metals is assumed to proceed on a sub-ps timescale [72][73][74]. Thus, one would naively expect the Cooper-pair breaking to be finished within a few 100 fs. The observation of slow pair-breaking suggests an alternative scenario, where bosonic excitations (e.g. phonons), generated during the initial decay of hot carriers towards the gap-edge, are responsible for pair-breaking. (c) Finally, in NbN the superconducting state is almost completely recovered on a timescale of a few hundred ps. While the timescale is much shorter than the SC recovery timescales deduced from continuous excitation experiments on Al or Sn [27,28], this observation does not prove the early estimates wrong. Based on available data from time-resolved studies on different superconductors, ranging from conventional BCS superconductors to high-T c cuprates, one can conclude that the characteristic SC recovery timescale is roughly inversely proportional to the size of the superconducting gap.
The most interesting and insightful is however following changes in dynamics as a function of base temperature and excitation density. Here, the data reveal that both pairbreaking rate and the superconducting state recovery timescale depend on both, the excitation density and temperature [20,22,61,62]. These dependencies provide access to the nature of the underlying processes and material specific microscopic parameters illustrated in Fig. 1.
From Fig. 3 it follows that in NbN [22,62] (the same is true for MgB 2 [20] and Pr 2−x Ce x CuO 4 [61]), the Cooper pair-breaking process ( ≈ 10 ps) and SC state recovery ( ≈ 100 ps) are well separated in time, which allows for a detailed view of the underlying microscopic processes based on their temperature and excitation density dependencies. Below we address the two processes in more detail.

Cooper Pair Breaking Dynamics: Determination of the Electron-Phonon Coupling
The rise-time dynamics is clearly reflecting the pair-breaking processes (i.e. the initial reduction of condensate density and an increase in the density of QPs). The pairbreaking was found to be temperature dependent (the process is faster as T is increased) [20], and photo-excitation density dependent (see Fig. 4) [20,22]. Especially the latter provides strong indications that the delayed pair-breaking is caused by high-frequency phonons, released during the sub-ps relaxation of high energy electrons [74][75][76]. As shown in Fig. 4, in both NbN and MgB 2 , at low excitation densities the pairbreaking can take tens of ps, speeding up upon increasing the excitation density. To account for these observations, in particular the excitation density dependence of the pair-breaking process, we consider the interplay between the condensate/QPs and highfrequency ( � > 2Δ ) phonons as proposed by Rothwarf and Taylor [29]. Within the model, the dynamics of the QP and high frequency phonon (HFP) densities, n and N, are described by a set of two coupled differential equations [29]: Here is the probability for pair-breaking by HFP absorption, R is the bare QP recombination rate with the creation of a HFP, N T is the concentration of HFP in thermal equilibrium at temperature T, and their decay rate. The rate could be governed either by anharmonic decay of HFP [15], since < 2Δ phonons do not have sufficient energy to break Cooper-pairs, or by the diffusion of HFP into the substrate [77]. Finally, I 0 and J 0 represent the initial densities of excess QPs and HFPs, which are generated during the avalanche relaxation of hot QPs towards the gap edge.
Now, let us imagine that the decay of HFP is the limiting step, governing the recovery of SC, i.e. the so-called phonon-bottleneck scenario by Rothwarf and Taylor [29]. Since the SC recovery dynamics proceed on a much longer timescale (100s of ps), one can-in the first approximation-neglect the decay of HFP given by the last term in Eq. (4), to reduce the equations to In this case, the early stage of relaxation will result in thermalization of QPs and HFP, which will reach a quasi-equilibrium state, described by a detailed balance equation, Rn 2 T * = N T * , where T * will be some new effective temperature (as in the model proposed by Parker [33]). For a SC gap much smaller than the Debye energy and the phonon density of states D( ) = 9 2 ∕ 3 D the HFP density is given by [15]. On the other hand, the QP density is given by . This means that at temperatures, where k B T ≪ Δ , in equilibrium most of the energy is actually in the QP channel (note the 2Δ∕k B T vs Δ∕T arguments in the exponents). Moreover, from the detailed balance equation it follows that ∕R = n 2 It is instructive to investigate the solutions of Eq. (5) for different excitation conditions. For simplicity, we can consider the low temperature limit, where the density of thermally excited QP's and HFP can be neglected. Secondly, we note that the ratio  [22]. The dashed lines are fits to the data with Eq. (7). b Superconducting depletion dynamics 7 K taken at various fluences (in μJ/cm 2 ) on MgB 2 [20]. Solid circles represent the data obtained by energy integrated approach (see Fig. 3b), while the open circles correspond to Δ ( ) measured directly. Again, the dashed lines are fits to the data using Eq. (7) (Color figure online) of microscopic constants ∕R has the dimensionality of concentration. Hence, one can introduce dimensionless QP and HFP concentrations, q ≡ Rn∕ and p ≡ RN∕ , while ≡ t is the dimensionless time [78]. With this, Eq. (5) read: with the initial conditions p(0) = p 0 , q(0) = q 0 . These coupled equations have analytic solutions q( ) and p( ) [20,78]. Here, we limit the discussion to the dynamics of QP density, q( ) , which can be directly compared to the experimental data Δ(t) in Fig. 4. The analytic solution for q( ) is . It is further instructive to take a closer look at the meaning of the dimensionless parameter K. Obviously, K = 0 corresponds to the case when immediately after photo-excitation the stationary solution is realized, i.e. when . In this case, the time evolution of the QP density, or of the reduction of Δ , would be a step function. The regime 0 < K ≤ 1 corresponds to the situation when q 0 > q(∞) , while −1 ≤ K < 0 represents the q 0 < q(∞) . For the latter case, p 0 > p(∞) , which results in an additional pair-breaking until q 2 (∞) = p(∞) is reached. Figure 5 presents a few simulations of q( ) which illustrate different limiting cases. Note that the function (7) was convoluted with the Heaviside step-function with the rise time of 100 fs, which accounts for photo-excitation and avalanche initial relaxation/multiplication process of hot quasiparticles. We consider a superconductor with Δ = 4 meV and the absorbed energy density A [x eV∕ucv] distributed between QPs and HFPs such that 2Δ and q 0 = R Ax Δ , where x is the fraction of the absorbed energy density in the QP channel ( 1 ≥ x ≥ 0). Figure 5a displays the evolution of the QP density for a constant excitation density and varying x, demonstrating the slowly rising signal due to the excess density of HFP (compared to p(∞) ). It is noteworthy that only for x → 1 the situation is realized that the QP density is decaying with time in this so-called pre-bottleneck regime. This actually reflects the ratio of n T ∝ e −Δ∕k B T and N T ∝ e −2Δ∕k B T discussed above. In other words, for k B T << Δ in quasi-equilibrium n T ≫ N T , and nearly all of the absorbed energy ends up in the QP channel. Therefore it is no surprise that even if only 20% of the absorbed energy is initially in the HFP channel, the delayed buildup of QP density (or suppression of Δ ) can be observed.
Panel (b) presents the dependence of pair-breaking on excitation density, with all other parameters being constant. With signal normalized to A the dependence of the apparent pair-breaking time on A reveals the behavior observed in experiments (see Fig. 4). Sub-linearity in the value of q(∞) , i.e. the photo-induced reduction of the gap, at high excitation densities is also revealed. Finally, panel (c) presents the (6) dq∕d = p − q 2 and dp∕d = −p∕2 + q 2 ∕2, evolution of the dynamics when varying the recombination rate (R) while keeping the pair-breaking rate constant. Using this analysis on NbN (the fits to the data are shown by dashed lines in Fig. 4a, where Δ is known, and it is demonstrated (see Sect. 2.1.2) that the absorbed energy can indeed be simply distributed between the QPs and HFP, one can extract and K and the fraction of the energy density initially stored in the QP channel x. The best fit is obtained when x = 0.09 , giving the values of the microscopic constants −1 = 6 ± 1 ps, R = 160 ± 20 ps −1 unit cell.
Since the pair-recombination rate R can be expressed in terms of microscopic [50,79], where Ω c is the phonon cut-off frequency and N(0) is the density of states, the value of can be determined. Taking the known values for Δ , N(0) and Ω c = 16 THz [80], we obtain = 1.1 ± 0.12 , which is in excellent agreement with the theoretical estimates, = 1 − 1.12 [81,82]. Importantly, such an approach does not suffer from the questionable assumptions [74][75][76] of the so-called two-temperature model, the model commonly used to determine from time-resolved experiments [72,73].

Energetics of the Gap Suppression in Conventional Superconductors
As addressed above, when at k B T << Δ the absorbed energy is distributed between the QPs and HFP and the quasi-equilibrium is reached, most of the energy is in the QP system. This fact provides the means to estimate the condensation energy of a superconductor. Figure 6a presents the dependence of the (maximal) relative change of the SC gap ( Δ∕Δ ) in NbN recorded at 4.3 K as a function of absorbed energy density (in mJ/cm 3 ). At low excitation densities Δ∕Δ increases linearly with A, followed by Simulation of the pair-breaking dynamics using Eq. (7) in the T → 0 K limit (also valid when photo-excited carrier density much higher than the thermal one). For all simulations Δ = 4 meV and the pair breaking rate is = 0.1 ps −1 . a presents the dynamics for constant excitation density, R = 100 ps −1 ucv and different fractions of the absorbed energy in the QP channel ( 0 < x < 1 ). b Excitation dependence of dynamics for x = 0.1 and R = 100 ps −1 ucv , note that the signal is normalized to excitation density. c The dependence of dynamics on the recombination rate and at constant A (Color figure online) a saturation (plateau at Δ∕Δ = 1 is reached when superconductivity is completely suppressed). Following the discussion of the coupled QP and HFP system above, it becomes clear that for the complete collapse of the gap, where QPs would be in a quasi-thermal equilibrium with HFP, the required energy would equal the thermal energy to heat up the sample from the base temperature to above T c . This follows from the fact that -in the absence of the gap-all phonons are the HFP (the energy stored in the HFP system in this case will be dominant). On the other hand, the rate at which the gap (or condensate density) gets suppressed in the limit of weak perturbations should provide an estimate of the condensation energy. Alternatively, this may be particularly relevant for bulk samples where the absorbed energy density varies with the depth, an estimate of the characteristic A required to suppress SC, one commonly uses a phenomenological saturation model, Here, A dep is the characteristic absorbed energy density required to deplete SC. It follows from Fig. 6a Here, B c is the thermodynamic critical field and N(0) is the electronic density of states at the Fermi level. For NbN B c = 0.234 T [80], giving E c = 22 mJ/cm 3 , nearly identical to A dep (4.3 K) . Moreover, as shown in Fig. 6b, the T-dependence of A dep is found to follow the T-dependence of Δ 2 (dashed line). We note that the observation in NbN is in strong contrast to similar studies in cuprates [52,[60][61], where A dep is found to be about one order of magnitude higher than E c and is nearly independent on temperature [60,61].

Superconducting State Recovery Dynamics at Temperatures Far Below T c
Now, let us turn our attention to the dependence of the SC recovery time on temperature and excitation density, which is for NbN shown in Fig. 7.  The T-dependence of the energy density required to deplete superconductivity, A dep , compared to the T-dependence of Δ 2 (dashed black line) (Color figure online) Figure 7a presents the relaxation dynamics recorded at 4K over 2 orders of magnitude in excitation density. While at the lowest excitation densities the recovery time, rec , substantially exceeds the data acquisition window, one can clearly see speeding up of relaxation dynamic upon increasing excitation density. As shown in Fig. 7b the dynamics is below ≈ 10 K both temperature and excitation density dependent.
As shown in Fig. 7c, at 4 K −1 rec first increases linearly with A, mimicking the intrinsic bi-molecular kinetics of the QP recombination discussed in Sect. 1.1. However, given the fact that the early time dynamics in NbN reveals a delayed pair-breaking with build-up of quasi-equilibrium between QPs and HFP (see Sect. 2.1.1), it is clear that the recovery dynamics is determined by the decay of the HFP population. That is, the SC recovery dynamics should be described within the phonon bottleneck scenario of Rothwarf and Taylor [29,78].
Before addressing this, we should note that such excitation density and temperature dependence of SC state recovery was first observed in cuprate superconductors [58,83,84], where it was argued to be an evidence for the absence of phonon (boson) bottleneck. It was stated that in the case of a phonon bottleneck the relaxation dynamics should be, at low temperatures, roughly constant and excitation density independent [83,85]. Surprisingly, and perhaps counterintuitive, considering the full solution of the coupled Rothwarf-Taylor equations, Eq. (4), it has been demonstrated that even in the case of phonon-bottleneck the relaxation should display excitation density and temperature dependence [78] as observed in NbN (Fig. 7).
While for the complete derivation of analytic solutions for SC state recovery we refer to Ref. [78], we briefly sketch the main arguments here. We remain with the description using dimensionless quantities as above. First of all, including the phonon decay term, the Rothwarf-Taylor equations read where ̃≡ ∕ is the dimensionless HFP decay rate and p T is the thermal density of HFP. At large time delays, ≥ 1∕̃ , both p( ) and q( ) are obviously slowly decaying, with the decay being determined by the HFP decay. Thus, at large times the difference q 2 − p ≈̃(p − p T ) . Considering that the difference q( ) 2 − p( ) ≪ p( ), q( ) and that it is decaying slowly, an analytic approximate solution was derived [78], describing the recovery dynamics: Here q T corresponds to the thermal QP density at base temperature and q(∞) is the value after the initial pre-bottleneck dynamics. Combining solutions for pre-bottleneck case and recovery dynamics an approximate solution for q( ) , valid over the entire time delay has been derived [78], where with q 1 ( ) given by Eq. (7) and q 2 ( ) by Eq. (9). The comparison of the analytic solution, given by Eq. (10), to the numerical solution of Eq. (8), shown in Fig. 8, manifests an excellent agreement for all different limiting cases. Figure 8 shows that, irrespective of the regime (strong or weak bottleneck), the SC state recovery dynamics displays excitation density dependence at low temperatures. Indeed, at T → 0 K, the time dependence of the inverse of the photo-induced QP density is a linear function of , mimicking the pure bimolecular kinetics. At higher temperatures, the relaxation becomes nearly exponential (dashed curve in Fig. 8d presents the exponential decay fit). We note that similar behavior was reported in cuprates [58].  For the case of a strong bottleneck ( ∕ ≡̃≪ 1 ), the initial relaxation rate −1 rec is given by [78] (the expression here is valid for k B T ≪ Δ) It is clear that this expression qualitatively accounts for the observations in NbN (Fig. 7). In the limit of very weak excitations, where q(∞) ≈ q T , the relaxation time should display exponential temperature dependence with −1 rec ∝ q T . Indeed, this accounts also for the early studies of temperature dependence of the relaxation time using biased tunnel junctions [27,28].
Moreover, the intensity dependence is observed if q(∞) ≫ q T . In this case −1 rec is proportional to q(∞) , which is in turn proportional to the excitation intensity (see Fig. 5). It follows that the excitation density dependence of the relaxation time should be observed at low-temperatures, as demonstrated in cuprates [58,61,83] as well as in conventional SCs like NbN [22], and MgB 2 [86]. Further argument for the validity of the Rothwarf-Taylor description of the light-induced gap dynamics in superconductors follows from the comparison of the experimental results on NbN (Fig. 7) with the model predictions. The data on NbN show that the excitation and temperature dependence of rec is observed for A ≤ 5 mJ/cm 3 and for T ≤ 10 K. For higher A and T ≥ 10 K −1 rec is constant at −1 rec (0) ≈ 0.01 ps −1 , likely governed by the escape of the HFP into the substrate (since it was found to be inversely proportional to the film thickness [87]). In a phonon bottleneck case −1 rec is proportional to the density of photo-excited QPs, n PI , if n PI falls into the range n T < n PI ≤ ∕R . Here, ∕R is the material-dependent characteristic QP density, below which the Rothwarf-Taylor perturbative description is applicable [78]. Indeed, A = 5 mJ/cm 3 corresponds to n PI ≈ 0.001 unit cell −1 , identical to ∕R ≈ 0.001 unit cell −1 determined from the analysis of the pair-breaking dynamics [22]. Similarly, the fact that −1 rec is constant above ≈ 10 K follows from the fact that n T (10K) ≈ 0.0003 ucv −1 is comparable to ∕R . Therefore, as experimentally demonstrated for a standard BCS superconductor, the A-dependent SC state recovery at low temperatures seems to be an expected behavior for a superconductor driven out of equilibrium.
We note that, while the phenomenological Rothwarf-Taylor has been derived to describe the low-temperature dynamics in superconductors [29,78], it was demonstrated to be able to qualitatively describe also the dynamics in other systems with small energy gaps near the Fermi level , i.e. when the gap is compared to the cut-off frequency of the coupled bosonic excitation [88,89]. For example, in heavy fermion metals a dramatic (exponential) slowing down of carrier relaxation was observed at low temperatures [88,90], originally attributed to slowing down of electron-phonon thermalization [90,91]. Systematic studies on both heavy fermion metals and Kondo insulators revealed that the dynamics can be accounted for by the presence of the hybridization gap near the Fermi level, with the dynamics of QPs and coupled HFP following the Rothwarf-Taylor bottleneck scenario [89,[91][92][93][94]. (1 + 2̃) .

Superconducting State Recovery Dynamics in the Limit T → T c
In the early subsections, we were focusing on the dynamics at low temperatures, i.e. for k B T ≪ Δ . Now we address the behavior near T c in the limit when Δ → 0 . As shown in Fig. 9a, already early experiments using real-time response of the current-voltage characteristic of a tunnel junction to fast laser pulses [96] revealed the anomalous slowing down of relaxation upon approaching T c . This can be interpreted by considering the relaxation in terms of the Ginzburg-Landau theory of the secondorder phase transition. Here, the restoring force is given by the derivative of the free energy with respect to the order parameter, giving rise to a divergence of rec ∝ 1∕Δ . Considering the phonons being in thermal equilibrium, Schmid and Schoen derived the expression to describe the dependence of the superconducting state recovery time near T c [97]: where E is the inelastic electron scattering time. Similarly to early measurements using tunnel junctions, the divergence of superconducting state recovery time was observed in superconductors using optical pump-probe spectroscopy. Typical results obtained on YBa 2 Cu 3 O 7 are shown in Fig. 9b. Assuming the SC state recovery being governed by the decay of HFPs, an alternative scenario, base on the quasi-equilibrium between QPs and HFP was put forward [15]. Here, Kabanov et al. [15] argued that the relaxation rate of the photoinduced QPs is dominated by the energy transfer from HFP ( � > 2Δ ) to phonons with � < 2Δ , i.e. by an anharmonic decay of HFP. Considering the kinetic equation for phonons, taking into account phonon-phonon scattering (anharmonicity), the expression describing the relaxation of the equilibrated QP-HFP temperature was derived: Here Γ is anharmonic phonon decay rate given, e.g. by the phonon linewidth in Raman experiments, N(0) is the normal state electronic density of states, and A the absorbed energy density. Indeed, for YBa 2 Cu 3 O 7 a good agreement of this formula with data (Fig. 9b) is obtained. Finally, we should mention another aspect of Kabanov's theory [15] , addressing the temperature dependence of the photo-induced QP density in the vicinity of the critical temperature, when QPs and HFP are in quasi-thermal equilibrium. Assuming that the energy gap is more or less isotropic, the temperature dependence of the photo-excited QP density n pe (which corresponds to R [q(∞) − q T ] discussed in Sect.

2.1.1) is given by:
where ℏΩ c is the phonon cut-off energy and the effective number of phonon modes per unit cell (or phonon branches) participating in the relaxation. Using this expression, and considering the BCS dependence of the superconducting gap Δ(T) , this expression is often used to determine the gap magnitude even for systems, where Δ(T → 0) may not be spectroscopically determined. Examples of such analysis are studies of the superconducting gap in cuprates [14,50,[98][99][100][101][102] (see also Sect. 2.2), pnictides [54,55,103] and organic superconductors [104]. Furthermore, a similar expression [15] has also been derived for the case where the energy gap is temperature independent, e.g. when discussing pseudo-gaps in cuprate superconductors [14,50,[98][99][100][101][102].

Superconducting Gap Dynamics in High-T c Cuprate Superconductors
In the above sections we focused on the real-time dynamics in conventional BCS superconductors with s-wave (isotropic) gap. Cuprate high-T c superconductors, however, are characterized by an anisotropic d-wave gap with nodes (directions in the momentum space with vanishing gap). Thus, one may expect their behavior under non-equilibrium conditions may be substantially different. Surprisingly, however, there seem to be more similarities than differences between the two classes of superconductors in terms of dynamics. Instead of a detailed review of the field of dynamics in cuprate superconductors, which can be found in Ref. [23], we briefly review the similarities and anomalies with respect to the observations in conventional BCS s-wave superconductors (discussed in sections above).
, Figure 10a presents the typical data obtained on cuprate superconductors in the low excitation density limit. It shows the temperature dependence of photo-induced reflectivity ( ΔR∕R ) traces on La 1.9 Sr 0.1 CuO 4 single crystal ( T c = 30 K) [99]. The data are often described in terms of a two-component response, with component A (see Fig. 10a) present only below T c . Component B, however, is present to temperatures far above T c . Probe polarization [98,102] and chemical doping [14,99] dependent studies revealed the coexistence of the two components below T c . On the other hand, the correlation of the (doping-dependent) onset temperature of the component B [14,99] has been made with the so called "pseudo-gap temperatures" determined by other methods, sensitive to the changes in the low energy density of states [105]. These suggest the origin of the component B being in the response of the competing order, like, e.g. (fluctuating) charge stripes or density waves.
From here on, we focus on the component (A) that can be clearly linked to the dynamics of the superconducting condensate, as demonstrated in Fig. 11. Here, a comparison of the dynamics, recorded by (a) probing the changes in reflectivity at 800 nm and (b) the THz conductivity as a function of excitation density (T = 4 K) is presented. At low excitation densities, where the fast relaxing component (B) associated to the dynamics of the competing order is small compared to the response of the condensate, both responses (reflectivity changes and changes in the THz conductivity) are comparable. Only for high excitation densities, when superconductivity becomes strongly suppressed, i.e. when the photoinduced change is saturating, the difference between the two responses becomes obvious. We note that the fast relaxing component B is absent (or dramatically reduced) in the THz response. The reason for that lies in the nature of the THz response of a superconductor, which is governed by the change of the phase of the transmitted THz electric field pulse The temperature dependence of relaxation time of the reflectivity signal. c The temperature dependence of amplitudes of the two components, including fits with the model by Kabanov et al. [15]. For details, see the original publication [99], from which the data are reproduced (Color figure online) (caused by the inductive character of the complex optical conductivity in the superconducting state)-see inset to Fig. 11b. Now, let us address the similarities (and some surprising inconsistencies) of the condensate dynamics in high-T c cuprates with conventional superconductors.

Pair-Breaking Dynamics in Cuprates
Looking at Figs. 10a and 11a, b, one realizes a similar delayed pair-breaking dynamics as in NbN and MgB 2 shown in Fig. 4. While the timescale is substantially shorter (depends also on the compound), one can argue in favor of pair-breaking by absorption of bosonic excitations created during the avalanche relaxation of hot QPs. This result suggests the same type of boson-bottleneck governing the relaxation as in conventional superconductors. Here, considering the possible scenarios of superconductivity in cuprates being mediated by phonons or magnetic excitations, we cannot specify a priori which bosonic excitation may be responsible for pair-breaking (or pairing). Note that the THz electric field is mainly sensitive to superconductivity, as exemplified in inset to b. At long time delays (100 ps) the system is already thermalized, and the corresponding response is a bolometric one. Indeed, the values for characteristic fluence require to heat up the film to above T c , F th , perfectly match the thermodynamic calculation [60]. In d, f the presented data are taken from Ref. [60] (Color figure online) We should further mention the study on YBa 2 Cu 3 O 7 crystals, where the dynamics of both, the condensate and specific IR-active c axis phonons were studied using time-resolved broadband THz conductivity [59]. The dynamics in the superconducting state is presented in Fig. 12. Panel (a) presents the in-plane response, showing a rapid filling-in of the superconducting gap (equilibrium 1 ( ) in the superconducting and normal state are shows at the top) and its recovery on the timescale of a few ps. The dynamics is similar to the results of earlier studies using near-infrared probe [14] of 0.2-1 THz pulses [17]. Panel (b) presents the out-of-plane (c axis) response, involving the dynamics of two c axis polarized phonons. Especially the apex oxygen mode at ≈ 17 THz is found to be strongly affected by superconductivity, giving rise to a highly asymmetric line-shape in the superconducting state. The analysis of the time-evolution of the phonon line-shape in panel (c) provides access to the dynamics of the condensate, related to the phonon asymmetry parameter shown in (d), and the phonon frequency. The later, presented in panel (e), displays a rapid red-shift of the mode and its sub-picosecond recovery. Since the red-shift is characteristic of an anharmonic lattice potential, the shift of the apex phonon frequency can be used as a probe of the vibrational occupation. Indeed, the presented frequency shift in Fig. 12e suggests a rapid heating of the mode, reaching occupations compared to those achieved by heating up to 200 K [59]. The results indicate highly excited lattice modes on the timescale comparable to pair-breaking.

Superconducting State Recovery Dynamics in Cuprates
The relaxation time in many of the cuprate superconductors shows similar temperature [50,83,99,106] (see also Fig. 10b) and excitation density dependence [58,61,83] as observed in NbN [22] and predicted by the Rothwarf-Taylor model for both weak and strong bottleneck cases [78]. Together with the observed delayed pair-breaking, discussed above, the data seem to be generally consistent with the strong boson-bottleneck scenario. On the other hand, some doping-dependent studies of carrier relaxation dynamics suggest that the character of relaxation may change between the under-doped and the over-doped regions of the phase diagram.
In particular, Gedik et al. [52] observed that the relaxation time is temperature and excitation density dependent in the under-doped case, while in over-doped at low temperatures the gap recovery dynamics is independent on excitation density and temperature. Indeed, the observation in Bi 2 Sr 2 Ca 1−y Dy y Cu 2 O 8 series [52] is consistent with the observation on YBa 2 Cu 3 O 7− samples, where in the over-doped range the relaxation dynamics at low temperatures seems to be independent on temperature and excitation density [14]. On the other hand, systematic doping dependence studies on single-layer cuprates do not seem to show this disparity [99,106].

Energetics of the Gap Suppression in Cuprate Superconductors
There are two important aspects of dynamics in cuprates that are revealed by studying the excitation density dependent suppression of superconductivity. First of all, in the low excitation density limit, the photo-induced changes in the QP density in cuprates scale linearly with the excitation density. Such an observation is straightforward for the case of a superconductor with an isotropic gap. Here the number of broken pairs per absorbed photon with energy ℏ is simply given by n pe = ℏ ∕2Δ (or a fraction thereof-see Sect. 2.1.1). But, for an anisotropic gap with nodes this observation is unexpected [15]. Here, the quasiparticle density of states for energies ≪ Δ a , where Δ a is the maximal gap in the anti-nodal direction, is given by N( ) = N(0) Δ a (power depends on the topology of the nodes on the Fermi surface). Therefore, at low temperatures, a sub-linear behavior of the photo-induced QP density with excitation intensity is expected [15], with n pe ∝ A ( +1)∕( +2) . For quasi 2D cuprates, with = 1 , this would suggest n pe ∝ A 2∕3 dependence at low temperatures, which was never observed. In fact, even the temperature dependence of n pe extracted from time-resolved experiments, seems to be at odds with the behavior expected for a d-wave superconductor [15,107]. A possible origin of this apparent contradiction may be in a slow relaxation of hot QPs towards the nodes due to kinematic constraints (i.e. the energy and momentum conservation) [108].
Another interesting aspect related to the excitation dependent studies of suppression of superconductivity in cuprates is related to the energy required to deplete superconductivity. As shown in Sect. 2.1.2, the absorbed energy density required to suppress superconductivity roughly matches the condensation energy in NbN [22]. Systematic studies on cuprates [52,60,61], however, show that the energy required to suppress superconductivity exceeds the condensation energy by about one order of magnitude. Figure 11 presents the fluence dependence of the superconducting state depletion in La 1.84 Sr 0.16 CuO 4 using near-infrared pump pulses in two configurations. The saturation of the photo-induced change in reflectivity or THz conductivity at 4 ps time delay is observed at the fluence ≈ 7 μJ/cm 2 , which corresponds to the absorbed energy density of 2.4 k B K per Cu atom. This energy is about 8 times higher than the condensation energy ≈ 0.3 K/Cu [109]. One possible reason for this discrepancy is a large magnitude of the superconducting gap in La 2−x Sr x CuO 4 , compared to the Debye energy. In other words, in La 2−x Sr x CuO 4 there are several optical phonon branches with � < 2Δ a ( Δ a is the gap energy in anti-nodal direction), thus a large amount of the absorbed optical energy can be transferred to phonons with ℏ < 2Δ a , on the 100 fs timescale. Since pair-breaking by phonons with ℏ < 2Δ a is strongly suppressed, these could serve as an energy sink during the initial avalanche relaxation of hot carriers towards the gap edge.
Recently, a similar study in the electron-doped cuprate superconductor Pr 2−x Ce x CuO 4 was performed [61]. Here melting of the superconducting state with both near-infrared and THz pulses was studied. Using THz pulses, tuned just above the gap, the absorbed energy density required to deplete SC was found to match the condensation energy. With near-infrared excitation, however, the absorbed energy density was again a factor of 5-6 higher than the condensation energy-despite the fact that in Pr 2−x Ce x CuO 4 the gap, 2Δ , is small compared to all relevant bosonic energy scales [61]. The data imply that following optical pumping a rapid electron-boson energy transfer takes place, yet only selected bosonic modes (e.g. antiferromagnetic fluctuations, or specific lattice modes) couple to the condensate. Further systematic studies, with selective (tuned) pumping are clearly required here.

Open Questions
Most of the observations obtained by systematic temperature and excitation dependent studies of the dynamics in the superconducting state show great similarities between cuprate superconductors and conventional BCS superconductors. Both pair-breaking and the superconducting state recovery dynamics are found to be consistent with a boson-mediated pairing in cuprates. The nature of the pairing "glue" is, however, still an open question. On the other hand, unlike in BCS superconductors, photoexcitation does not seem to give rise to a reduction of the gap magnitude (as seen for NbN in Fig. 3), but rather to filling-up of the gap, i.e. to an increase of the spectral weight inside the gap [68].

THz Driven Enhancement of Superconductivity in BCS Superconductors
In Sect. 1.2 we addressed the early studies of enhancement of superconductivity with sub-gap excitation. Most of these experiments were conducted using continuous excitation with sub-gap radiation, and the effects (increase of the superconducting gap, critical current and transition temperature) were typically of the order of 1%, limited by the QP scattering rate and the radiation induced sample heating. Using the novel pulsed sources of intense infrared radiation one may expect to observe transient effects, whose magnitude would exceed magnitudes observed in early studies. Here we briefly review a study, which suggest superconductivity enhancement effects in a dirty-limit superconductor NbN [110] on a 100 ps timescale, when the system is driven by narrow-band THz pulses tuned close to the gap energy.
In the study, summarized in Fig. 13, narrow band THz pulses ( ≈ 30 GHz bandwidth, 20 ps duration) were used to both excite NbN and probe the resulting changes in the transmission. The studies using near-infrared excitation [22] demonstrated that the time evolution of the NbN superconductor can be, for time delays beyond the pair-breaking process, described by the time evolution of effective temperature T* of thermalized condensate, QPs and HFPs. As the transmission through the sample is strongly affected by the opening of the superconducting gap, one has-in principle-access to the time evolution of the gap, without having spectral resolution. With these assumption, two types of experiments were performed: above gap pumping at 2.08 THz, with traces shown in panel (a), and pumping with 1.26 THz pulses (which is below 2Δ for T < 11 K-see inset to Fig. 13a) with traces at different temperatures shown in panels (c1)-(c6).
Two interesting observations were made, both suggesting an additional mechanism present at temperatures close to T c , competing with the expected light induced suppression of superconductivity.
In experiments with above-gap excitation, the absorbed energy density required to deplete superconductivity, A dep , was determined at different base temperatures, similar to experiments performed with near-infrared pumping (see Fig. 6). The surprising observation is that A dep is increasing with increasing temperature (Fig. 13b) opposite to the experiments with near-infrared (1.55 eV) pumping. This is a first  [110]. a Dynamics of the relative change in transmission Tr (normalized to Tr 4.3K -Tr 16K ) at 4.3 K when pumping with 2.08 THz pulses of different absorbed energy density A (in mJ/cm 3 ). The shaded region presents the pump-probe intensity auto-correlation. Inset shows the T-dependence of gap with the two frequencies used in narrow band THz pump-probe experiments by the dashed lines. b Comparison of the absorbed energy density required to deplete superconductivity, measured with near-infrared pumping at 1.55 eV [22] and THz pumping at 8.3 meV (2.08 THz). The dashed line represents the condensation energy. (c1-c6) Recorded traces of differential transmission for pumping at 1.26 THz at different temperatures (black solid lines). The shaded region presents the pump-probe intensity auto-correlation. For red dashed lines present the behavior in the case of the so called T * -model, which adequately described the dynamics for near-infrared excitation [22]. At each temperature the maximum relative change in the gap Δ∕Δ is estimated from the measured Tr∕Tr (Color figure online) indication that an additional process may be competing with the common lightdriven suppression of superconductivity.
Even stronger evidence of a competing process is provided by pumping at 1.26 THz. Panels (c1)-(c6) present the recorded traces of the relative transmission change, compared to the simulations assuming the T * model, i.e. the assumption that the superconducting state at ∼ 30 ps after photo-excitation can be described by Δ * = Δ(T * ). The simulations are given by the dashed red and dotted blue lines. The non-monotonous behavior is a result of the fact that at 1.26 THz the transmission is a non-monotonous function of temperature [110]. For time delays less than ≈ 50 ps, the relative transmission changes systematically shows the opposite direction to the one predicted by the T * -model. This is especially pronounced for near-resonant pumping, when 2Δ∕h ≈ 1.26 THz, around 10 K. This observation suggests a transient gap enhancement (and the corresponding increase in condensate density) under THz pumping for timescales less than ≈ 50 ps.
Given the fact that these effects are pronounced in the vicinity of critical temperature, when the density of thermally excited QPs becomes prominent, the data are interpreted within the Eliashberg scenario [36,110]. Here, narrow-band excitation gives rise to a highly non-thermal QP distribution resulting in amplification of SC that would be washed out under continuous excitation. Taking the data at face value, the enhancement effects seem to be observed even in the case of above gap excitation, where QP excitation competes with pair-breaking. One could imagine that such experiments, if performed with pulses at even lower carrier frequencies and closer to the critical temperature, may even lead to stronger effects as discussed above.

Photo-Induced Superconductivity in Unconventional Superconductors
Following the ideas of the pioneering works on radiation-enhanced superconductivity addressed in Sect. 1.2, we discussed first experiments utilizing intense narrow band pulses, aiming to transiently enhance superconductivity in conventional BCS superconductors. We should, however, also briefly address more recent studies reporting superconducting-like transient states (with the typical duration of 1 ps) observed in cuprate superconductors [111][112][113][114], fullerides [115][116][117][118] and recently also on organic charge-transfer salts [119] when photoexciting at base temperatures far above the corresponding equilibrium critical temperatures in these systems. The leading idea in these experiments is to excite specific lattice modes/distortions by utilizing nonlinear phononics [120,121]. This way, pulsed excitation can either suppress a competing (stripe) order [111,122], or induce superconductivity by the lattice-driven modification of the pairing interaction [111][112][113][114][115]. Figure 14 presents the experimental approach and the results obtained in one of the first studies of this type, the study of light-induced melting of the competing stripe order in La 1.675 Eu 0.2 Sr 0.125 CuO 4 ( LESCO 1∕8 ) [111]. The non-superconducting LESCO 1∕8 is known to possess a static stripe order. Upon excitation with infrared pulses, the induced changes in reflectance (Fig. 14c) display a shoulder-like structure around 60 cm −1 , which was interpreted as a photoinduced c axis Josephson plasma resonance (JPR) [111]. The JPR is a general feature of layered cuprate superconductors, explained by the Josephson coupling between the adjacent superconducting ( CuO 2 ) layers and can be observed in equilibrium studies of c axis optical conductivity, as shown in Fig. 14b. The fact that the induced changes in reflectance, Δr∕r , are small (the step in Δr∕r , attributed to the JPR is only about 0.1%, to be compared to the ≈ 40and the reflectivity of 30% effect in the equilibrium), could be attributed to a large mismatch of the penetration depths. Indeed, the pump field penetration depth at 16 μm is about 200 nm while that of the THz probe field is about 10 μm [111]. The comparably large, frequency independent offset in Δr∕r was however attributed to the photoinduced changes in the high(er) frequency phonons [111]. By processing the raw data, taking into account the large penetration depth mismatch and subtracting a frequency independent component originating from perturbed higher-frequency modes, the authors showed that 2 from the photo-excited layer displayed the upturn at low frequencies, resembling the response of the superconductor [111]. The fact that the effect was peaked around the pump-photon energy of 80 meV ( 15 μm ) leads to the interpretation of photoinduced superconductivity via melting of the underlying static stripe order [111]. The result is particularly surprising, considering the excitation densities in these studies are of the order of 1 mJ/cm 2 , nearly two orders of magnitude higher than excitation densities required to suppress superconductivity in an optimally doped LSCO, when excited by NIR pulses polarized in the a-b plane (see Fig. 11). Following this pioneering work, the approach has been extended to the superconducting cuprates. Using the same approach, it was shown that the transient state with the low-frequency divergence in 2 ( ) , reminiscent of 1∕ inductive response of the condensate, could be realized over large range of doping and temperatures, approaching room temperature in strongly underdoped cuprates [112,114,121]. The interpretation of these data is, however, still under intense discussion, ranging from aspects of the data analysis to the different interpretations of the origin of the lightinduced modifications of low frequency ( ) [123][124][125][126][127][128][129][130][131]. Recently, several groups performed similar studies of the c axis THz response using both mid-infrared as well as near-infrared pulses for photoexcitation [122,125,126,[132][133][134]. While it was demonstrated that near-infrared excitation results in similar changes in the low-frequency conductivity as with mid-infrared excitation [132,134], the consensus on the nature of the photoinduced state is still lacking. In the superconducting state of the underdoped LSCO Niwa et al. [125] observed a photoinduced red shift of the JPR, consistent with the photoinduced suppression of superconductivity. It was argued that transient heating linked to the inhomogeneous excitation profile could give rise to potential artefacts in the analysis [125]. Studies in the superconducting state of La 1.905 Ba 0.095 CuO 4 [132], La 1.885 Ba 0.115 CuO 4 [122], and Pr 0.88 LaCe 0.12 CuO 4 [133], however reported a photoinduced red shift of the JPR and the appearance of an additional mode at higher frequencies. While at temperatures above T c no such photoinduced mode was observed in La 1.905 Ba 0.095 CuO 4 [132], it has been observed in La 1.885 Ba 0.115 CuO 4 , albeit at very high excitation density of 9 mJ/cm 2 [122]. Finally, following up on the first reports, suggesting signatures of transient (inhomogeneous) superconductivity up to room temperature in underdoped YBa 2 Cu 3 O 6+x [112,114], similar studies have been performed by both mid-infrared and near-infrared excitation [126]. While, similarly to earlier reports [112,114], a photoinduced increase of THz reflectivity was observed, the authors argue against the photoinduced (inhomogeneous) superconductivity [126]. From the above, it is clear that further systematic studies are required. Ideally, the experiments should be performed on a axis oriented thin films in transmission geometry, to insure homogeneous excitation profile, and to enhance the experimental sensitivity (as shown in Sect. 2, superconductivity results in vehement changes in thin film transmission).
We note that similar studies to those in cuprates were performed also on K 3 C 60 [115][116][117][118] and recently also on the organic charge-transfer salt -(BEDT-TTF) 2 Cu[N(CN) 2 ]Br [119]. In both materials, upon photoexcitation with intense mid-infrared pulses, spectral signatures of superconductivity were reported for base temperatures well above the corresponding equilibrium critical temperatures. In K 3 C 60 the mismatch between the pump and the probe penetration depths was still large [115], requiring multilayer modelling (as in cuprates). This does not seem to be the case for the organic charge-transfer salt [119]. At face value, the extracted optical conductivities in the transient states indeed seem like to result from photoinduced superconductivity, i.e. showing gapping (reduction) of 1 and 2 ( ) ∝ 1∕ . In K 3 C 60 , where the recent study suggests the photoinduced state to be long-lived with lifetime well in the nanosecond range [118], the optical conductivity data is supplemented by the DC transport studies, demonstrating the photoinduced drop in resistivity [118].
Despite the apparent agreement with the proposed scenario of light-enhanced/ induced superconductivity, where mid-infrared pulses induce large-amplitude structural distortions and thereby modulate the pairing strength [115,118], there are still several open questions. For example, experiments on K 3 C 60 performed at temperatures below T c show a photoinduced suppression of superconductivity [115] instead of the expected gap enhancement. Apart from the proposed exciton-cooling model [129] this observation seems at odds with the majority of the proposed scenarios [115,127,[135][136][137]. Secondly, the excitation densities in most of these experiments are extremely high. For example, the recent study reporting the generation of a long-lived (nanoseconds) metastable superconducting state in K 3 C 60 [118] suggests a fluence threshold for the effect of about 20 mJ/cm 2 . With the pump penetration depth of 200 nm, and the reflectivity of 30% [115,118] one obtains the absorbed energy density of the order of 700 J/cm 3 . The resulting temperature of the excited volume, assuming the system is quasi-thermalized on the 1 ns timescale, can be estimated from the temperature dependence of the total specific heat of K 3 C 60 [138], and amounts to about 500 K! Here we should mention a recently proposed alternative scenario, where the experimental data are accounted for by considering that the pulse, resonant with a phonon mode, creates a (nonsuperconducting) nonequilibrium state in which the linear response conductivity becomes negative [130,131].

Summary and Perspectives
In this brief review, we addressed some aspect of superconductors, driven out of equilibrium using femtosecond and picosecond optical pulses. The real-time experimental approaches and theoretical models presented here, building on the pioneering works from the early era of non-equilibrium superconductivity (Sects. 1.1, 1.2), provide new insights into this fascinating field of condensed matter physics.
This review is, however, far from comprehensive in view of addressing all aspects related to light-driven out-of-equilibrium phenomena in superconductors. We briefly discussed novel approaches, aiming to impulsively modulate the pairing strength and thereby enhance or induce superconductivity at temperatures far above the superconducting critical temperature (Sect. 3.1.1). In addition, we should also mention the recent efforts in using intense THz pulses to probe collective modes of superconductors (for classification of collective modes in superconductors see, e.g. [139]).
The amplitude mode of the superconducting order parameter, describing variations of the amplitude of the order parameter with the frequency of 2Δ∕h , commonly referred to as the Higgs mode, is a scalar excitation of the order parameter [140]. As it does not carry spin or charge, the Higgs mode does not couple linearly to the electromagnetic field. Only when superconductivity coexists with a charge density wave (CDW) order, as in 2H-NbSe 2 [141][142][143], has Higgs mode been detected using linear spectroscopy [141,143]. In this case, the amplitude mode of the CDW order couples to the Higgs mode via modulation of the density of states at the Fermi level and thereby perturbing the superconducting condensate.
There were numerous theoretical reports, suggesting that the excitation of Higgs mode could be realized by quenching the superconducting order [144] on a timescale shorter than the mode frequency [145,146]. These were followed by experimental studies using intense THz pulses, utilizing the THz pump-THz probe approach [147,148] and the third harmonic generation [148,149]. These results suggest the Higgs mode to be observed through its nonlinear coupling to the intense THz light fields. On the other hand, it has been argued that the third harmonic generation (THG) can also be attributed to the nonlinear optical excitation of collective lattice-charge fluctuations [150]. It has been argued that the THG signal due to density fluctuations should in fact dominate the response in a BCS superconductor [151]. While the question remains open [148,[151][152][153], we should mention yet another approach, where the Higgs mode is observed in the linear optical conductivity spectrum under injection of supercurrrent [154,155]. Finally, in addition to the Higgs mode, unconventional superconductors with reduced dimensionality and multi-band superconductors may host a much richer spectrum of collective excitations [139,156,157]. The first experiments, suggesting the observation of the collective oscillation of the relative phase between the order parameters in a two-band superconductor MgB 2 , the Leggett mode, have just been reported [158].
Following the short summary, where we briefly introduced another field of studies of superconductors that has been made possible by the advancement in generation of intense electromagnetic pulses-the studies of the collective modes, we should reiterate some of the open questions and the possible ways of addressing them.
We have shown that the dynamics of the superconducting state of conventional BCS superconductors, excited by femtosecond optical pulses, can be extremely well described by the phenomenological Rothwarf-Taylor model [29,78]. Surprisingly, the same can be said about experiments in the weak perturbation regime in unconventional d-wave superconductors, despite the fact that they are characterized by the anisotropic gap with nodes. Both pair-breaking and the superconducting state recovery dynamics are found to be qualitatively consistent with a boson-mediated pairing in cuprates. The question about the nature of the pairing "glue" is, however, still open. Here, systematic studies, by varying the energy of excitation photons, could provide further clues. Recent studies on the optimally doped PCCO suggest that pairing can be mediated either by bosons or antiferromagnetic fluctuations [61]. Thus, doping-dependent studies may be able to distinguish between the two scenarios.
Moreover, photoexcitation of cuprate superconductors seems to give rise to filling-up of the gap instead of the suppression of its magnitude [68]. Here a major progress in understanding can be expected from systematic tr-ARPES studies with high energy resolution and large momentum coverage [159][160][161].
Another rapidly developing field of research focuses on amplification of superconductivity. Here, studies in both conventional [110] and unconventional superconductors [111-114, 121, 122] require further investigations. In NbN [110], narrow band THz excitation studies suggest superconducting state enhancement to compete Journal of Low Temperature Physics (2020) 201:676-709 with pair-breaking, yet the spectrally resolved study is still lacking. Such studies could provide further evidence of non-thermal quasiparticle distribution, supporting the proposed scenario of Eliashberg enhancement of superconductivity [110]. Moreover, by tuning the excitation photon energy, detailed information of the underlying processes could be obtained, and compared to different theoretical proposals [41][42]. As far as the studies of photoinduced enhancement of superconductivity in unconventional superconductors are concerned, the main advancement could be expected by studying optically thin films. This way, one could insure a homogeneous excitation profile, and dramatically enhance the sensitivity.
Finally, in many of the unconventional superconductors, their ground state properties can be controlled by doping, pressure or strain. Here, a combination of pressure-tuning and impulsive excitation should provide a new playground to generate and investigate emergent (metastable) states, especially when pressure-tuning the materials into the vicinity of (quantum) critical points in their phase diagrams.