arXiv:2010.05019v1 [cond-mat.supr-con] 10 Oct 2020 clrectto ihu hreo electric or is charge fluctuations without amplitude excitation gap. the scalar superconducting to the corresponding of mode phase Higgs fluctuatio and by amplitude characterized the is of state superconducting a ken, uecnutrNN[ NbN superconductor ssc, antdrcl opet xenlpoe nthe in probes external to couple [ directly regime linear-response cannot it such, As uny lhuhteefaue eeosre nthe in observed the were twice features fre- when these gap signal Although superconducting THG quency. the the matches in frequency peak pump-pulse a to rise give should smr ieyarsl fcag est utain to due fluctuations [ density breaking charge Cooper-pair of driven THz result the a likely more is asdmnn,a one u yfrhrtertclstudie theoretical further by al- out not pointed however, as is, THG dominant, the ways to fluctuations density charge of a supercon- to the addition below In signal temperature transition THG mode. ducting the Higgs of driven enhancement the major attribu to signal unambiguously THG to the straightforward to not leading however, pulse, [ is, pump (THG) the the generation of to third-harmonic frequency coupled the the further twice is with pulse mode pump parame Higgs order driven coher- superconducting The the the ter. trigger of and oscillations Higgs equilibrium, ent of out superconductor a 4 os[ tors ob ne an to ml steHgsapiuemd nasprodco [ superconductor a in mode amplitude Higgs prominent the A is often ample states. can non-equilibrium states, driven equilibrium in in visualized spectroscopy linear the to .Wt the With ]. eaet Tz amncgnrto a ensuggested been has generation harmonic (THz) Terahertz unu hnmn,wihmyntb accessible be not may which phenomena, quantum Novel 5 – 10 ffi .Hr,a nes arwadTzplecndrive can pulse THz narrowband intense an Here, ]. in rb orva ig oei superconduc- in mode Higgs reveal to probe cient ihtesprodcigtasto temperature transition superconducting the with below oevr 2 for Moreover, .,ad05Tz eprtr-eedn vlto fteTH the of 2 evolution where temperature-dependent temperature THz, 0.5 and 0.4, oteiprac fsatrn o bevto fteHiggs the of observation for scattering of importance driven the the to of contribution predominant the suggest results f Wang, U 3 alH .vnLoosdrecht, van M. H. Paul nentoa etrfrQatmMtras colo Phys of School Materials, Quantum for Center International = erpr ntm-eovdlna n olna eaet s terahertz nonlinear and linear time-resolved on report 1 oainlsmer pnaeul bro- spontaneously symmetry rotational (1) egyKovalev, Sergey 0 adSlcieTidHroi eeaini Superconduc in Generation Third-Harmonic Band-Selective 1 . nttt fRdainPyis emot-etu Dresde Helmholtz-Zentrum Physics, Radiation of Institute n . H ih2 with THz 0.7 and 6 T 5 c u Wang, Yue ydiigtesse ihitnenrobn H uss Fo pulses. THz narrowband intense with system the driving by 5 ,i a rudta h H signal THG the that argued was it ], 5 4 2 ple uecnutvt etradSaeKyLaboratory Key State and Center Superconductivity Applied ]. nttt fPyis oansGtnegUiest Main Gutenberg-University Johannes Physics, of Institute 7 auta hsk ehiceUiesta otud 4422 Dortmund, Universit¨at Technische Fakult¨at Physik, f vdnefrHgsApiueMd nteDryLimit Dirty the in Mode Amplitude Higgs for Evidence 4 < nttt fPyisI,Uiest fClge 03 Colog 50937 Cologne, of University II, Physics of Institute T 6 2 olbrtv noainCne fQatmMte,Beijin Matter, Quantum of Center Innovation Collaborative f c 1 ∆ h ig oeexcitation mode Higgs the , 5 = π a Dong, Tao iZa Gan, Zi-Zhao colo hsc,Pkn nvriy ejn 081 Chin 100871, Beijing University, Peking Physics, of School 2 h H sfudnal stoi ihrsett h pump-pu the to respect with isotropic nearly found is THG the ∆ π 11 f ( 4 T .Eprmnal it Experimentally ]. > .Tecontribution The ]. ,frtedrylmtsprodciggp2 gap superconducting dirty-limit the for ), / 2 4 antcdipoles. magnetic ∆ ogWu, Dong 2,3, π h H nest nrae oooial ihdecreasin with monotonically increases intensity THG the , ∗ 5 iY Shi, Li- eynGermanskiy, Semyon Dtd coe 3 2020) 13, October (Dated: s -wave 3 a-i Wang, Nan-Lin T 1 be ns - c – s - 3 ≈ hi Reinho Chris 6K hr-amncgnrto TG sobserved is (THG) generation Third-harmonic . 36 ig mltd oeo h it-ii a,pointing gap, dirty-limit the of mode amplitude Higgs [ us ihfrequency with pulse frequency n a enpeitdt lyadcsv oei h nonlinea the [ in superconductors role decisive the a of play response to predicted been has ing egt hs oe eutn rmteitrbn coupling inter-band so-calle the the from [ resulting also but mode, ampli- phase bands, Higgs Leggett individual the in only fluctuations not tude addressed superconductors band adydtcal naTGeprmn [ and experiment THG weaker a much in are detectable superconductor hardly fluc two-band phase a Leggett in the tuations while tha dirty fluctuations, predominant the density are charge in fluctuations the Thus, amplitude Higgs the scattering. limit, breaking impurity pair by by enhanced dominated icantly be to found [ was THG the limit electromagnetic the to couple charge linearly also [ field not is does mode, and Leggett neutral the parameters, order coupled hc sepce ob aietdb eoac ekin peak resonance a by manifested [ THG be the to expected is which ob nstoi [ anisotropic expected is be mode, fluctuations to density Higgs charge driven from the THG to the THG whereas the signature of a assignment field, with the driving signal porting the THG of isotropic polarization the an crystalli to exhibit to respect to respect found with was NbN field THz axis. driving the of larization di two the suggest Moreover, scenarios driv signal. competing THG the the to dominate will due mode contribution Higgs the pronounced, are scattering 11 13 10 oei superconductors. in mode etocp ftetobn uecnutrMgB superconductor two-band the of pectroscopy eety xesv hoeia tde fteTGi two- in THG the of studies theoretical extensive Recently, , – , inl xiisarsnnemxmma the at maximum resonance a exhibits signals G 17 12 15 c,Pkn nvriy ejn 081 China 100871, Beijing University, Peking ics, 13 -osnof 12 rse,Germany Dresden, 01328 n-Rossendorf, 4 ,6 3, .Hr,flcutoso eaiepaebtentetwo- the between phase relative of fluctuations Here, ]. .Epcal hncreaine correlation when Especially ]. ,tecnrbto fteHgsmd a esignif- be can mode Higgs the of contribution the ], .Smlrt h ig oe,teLgetmd with mode Leggett the modes, Higgs the to Similar ]. a-hitp Deinert, Jan-Christoph ueDemsar, Jure ff 14 ω , L ,518Miz Germany Mainz, 55128 z, – h upplefrequencies pump-pulse the r otud Germany Dortmund, 1 ∆ a ersnnl xie yamliyl THz multicycle a by excited resonantly be can 17 4 o eocpcPhysics, Mesoscopic for π i-unXu, Tie-Quan .Mroe,prmgei muiyscatter- impurity paramagnetic Moreover, ]. = 7 1 ]. e Germany ne, . 3Tza .I otat for contrast, In K. 4 at THz 03 ,China g, f 2 aifigtecniino 2 of condition the satisfying a n h Wang and s oaiain These polarization. lse igMgB ting 5 ff 1 Hong-Zhang rn eedne ntepo- the on dependences erent grIlyakov, Igor 17 temperature. g – 19 ff ,4 1, 4, 7, .Wiei h clean the in While ]. 17 f csand ects 2 = : ]. 0 † . 3, 1 2 / rimpurity or f = sup- ω ne L d n - r , 2

(a) (c) 13 nm thick MgB films grown on 5 × 5 mm MgO substrate MgB 2 Mg 2

2 of 0.33 mm thickness, with Tc ≈ 36K[20, 21]. Linear time-

B 36K domain THz spectroscopy was performed using GaAs photo- )/Tr w 1 conductive antenna as THz source, driven by Ti:sapphire os- Tr( 4 K

26 K

12 K cillator generating 35 fs pulses at 800 nm central frequency

c

36 K 20 K with the repetition rate of 80 MHz [23]. For the THG mea- 0

(d) b a surements, broadband THz pulses were generated via tilted s s

1 2 - 1 )

t -1 4 K (b) pulse-front scheme utilizing a LiNbO3 crystal [24, 25], driven 30

4 36 K cm -1 by a Ti:sapphire amplifier generating 100 fs pulses at 800 nm 4 K, MgO W

5

4 K, MgB /MgO 2 with pulseenergyof 1.5 mJ and repetition rateof 1 kHz. Spec-

2 20 (10 36 K, MgB /MgO

2

1,2 tral distribution of the generated THz pulses was optimized s for around 0.7 THz central frequency. Narrow-band multicy-

0

10 (e) cle THz pulses were produced using bandpass filters with 20%

x30

2D p bandwidth. The THz radiation was detected by electro-optic

1 0 sampling at a ZnTe crystal. Temperature dependent measure- Electric field (a.u) field Electric 1,36K ments were performed using helium-flow cryostats. (w)/s

-10 1

s The two-band BCS-type superconductor MgB2 crystal- lizes in a hexagonal structure with the space group P6/mmm

0

1 2 -1 0 1 2 3 [Fig. 1(a)] [26–31]. Two superconducting gaps of 2∆π = 1.1

Time delay (ps) Frequency (THz) and 2∆σ = 3.6 THz were foundin the π- and σ-bands, respec- tively [28]. Figure 1(b) shows THz electric field transients Figure 1. (a) Layered crystallographic structure of MgB2. THz elec- transmitted through the MgB2/MgO sample at 4 and 36 K, tromagnetic waves propagated along the c-axis. (b) Transmitted THz compared to that through the reference bare MgO substrate electric field as a function of time delay, for MgB2/MgO at 4 and at 4 K. Since the transmission through samples is much lower, 36 K and for the bare MgO substrate at 4 K. The data for MgB2/MgO at 4 and 36 K is magnified by a factor of 30 for clarity. (c) Ratio of the corresponding time-domain traces are multiplied by 30 for transmission for various temperatures with respect to 36 K. The ar- clarity. By Fourier transformation of the time-domain signals, rows indicate the frequencies at which the ratio curves cross one. (d) one can obtain transmission and phase shift in frequency do- Real σ1 and imaginary σ2 parts of optical conductivity at 36 and main, from which complex optical conductivity can be - 4 K. The arrow indicates the crosspoint of the σ1 and σ2 curves at rived. −1 τ = 1.34 THz for 36 K. (e) Ratio of σ1 for various temperatures Transmission ratio with respect to 36 K for various tem- with respect to 36 K. The dashed line indicates the excitation gap peratures below T is shown in Fig. 1(c) as a function of fre- 2∆ = 1.03 THz at4 K. c π quency. For all the temperatures the transmission ratio in- creases first with decreasing frequency, and then drops con- tinuously down to a value smaller than one. The curves of Motivated by these theoretical advances, we performed the different temperatures below Tc intersect with each other, THz spectroscopic measurements in the two-band supercon- forming an isosbestic point [32–34]. As indicated by the ductor MgB2 in the linear- and nonlinear-response regimes. dashed line, the isosbestic point around 1.03 THz coincides ∆ The linear response reveals the smaller gap 2 π to be in the with the onset of optical conductivity σ1 at4 K[seeFig. 1(d)], dirty superconductor limit. On the other hand, no spectro- which reflects formation of the superconducting gap in the π scopic signature is evident in the vicinity of the larger gap band as discussed below [28–30, 35]. 2∆σ in the accessible spectral range, suggesting it may be in Figure 1(d) shows the extracted complex optical conductiv- the clean-limit. The THz third-harmonic generation is stud- ity in the normalstate (36 K) and at 4 K, clearly demonstrating ied as a function of temperature for various driving frequen- opening of the gap below 36 K. In the normal state, the data cies and as a function of polarization with respect to the crys- are consistent with the free carrier Drude response. The lim- talline sample orientation. The temperature-dependent THG ited bandwidth prevents us from performing detailed analysis display peaks at temperatures, where 2 f matches 2∆π, while using two Drude components with reliable accuracy. Instead no anomalies are observed for 2 f matching either the fre- we estimate the scattering rate within the simple Drude sce- quency of the larger gap or of the Leggett mode. Together nario as the crossing point between σ1 and σ2. The extracted with the nearly polarization-independentTHG intensity, these value τ−1 = 1.34 THz is indicated by the arrow in Fig. 1(d). results provide strong evidence for the Higgs mode of the gap The estimated scattering rate in our sample is comparable to in the π-band. the values reported in Ref. [36, 37] but considerably lower High-quality single-crystalline MgB2 thin films with the c- than the values of 4.5 and 9 THz reported in Ref. [38] and axis epitaxy were grown on MgO (111) substrate by using a [35], respectively. hybrid physical-chemical vapor deposition method, and char- In the superconducting state the spectral weight of σ1 at acterized by X-ray diffraction and charge transport measure- the finite frequency is depleted due to the opening of the su- ments [20–22]. Optical measurements were carried out on a perconducting gap as shown in Figs. 1(d) and (e). At 4 K, 3

(b) σ1 drops almost to zero below the superconducting gap at (a) STM, 2 D

f = 0.5 THz s

3.5

80 2∆π = 1.03 THz ∼ 4.26 meV, while σ2 exhibits a 1/ω in- STM, 2 D p

Raman, w

L

3.0 THz, 2 D ductive response. We also note the upturn in σ1 at frequencies p Theory, w

L ∼ 40 below 0.7 THz, which is consistently observed in MgB2 2.5 at low frequencies [35, 38]. This could be caused by an ad- 2.0 0 THz) (

ditional absorption of unpaired carriers, which seems to be D

1.5 2 Electric field (kV/cm) field Electric intrinsic in nature [38]. -40 The observed superconducting gap and temperature depen- 1.0

-80 dence of the optical conductivity are comparable with that re- 0.5

0 10 20 30 ported for thin-film samples in Ref. [35, 38]. Quantitatively, 0.0 2 f = 2 D

Time delay (ps) broad-band normal state optical response was modelled also p

(d) (c)

4

f using two Drude terms, from which different scattering rates 1 = 0 . 3 T H z

for the σ- and π-bands were derived [39, 40]. Similar to the 0 . 4 T H z − 1 3

∆ ∆ (norm.) case in Ref. [36, 37], our data suggest that 2 π < τ < 2 σ. f 3 I Note that the scattering rates extracted in literature are gener- 0 . 5 T H z ally higher in the π-band than in the σ-band [39, 40]. Thus, 2 (arb. units) (arb.

0 . 6 T H z 4 p D

T /

c our extracted value likely corresponds to the scattering rate of 2 |

f = 0 . 7 T H z

(3) 1

the σ-band. This suggests that superconductivity can be con- c THG intensity intensity THG | sidered to be in the dirty limit for the π-band and in the clean f = 0.3 THz

0 limit for the σ-band. The contrast between the two bands ac- 0

0.5 1.0 1.5 2.0 count for the fact that only the smaller gap is clearly observed 10 20 30 40

2f/2 D (T) Temperature (K) using optical spectroscopy, and seems crucial also for the un- p derstanding of the nonlinear response, as presented in the fol- Figure 2. (a) Electric field of f = 0.5 THz pump pulse in time lowing. domain. (b) Temperature dependence of superconducting gaps as Due to the low-frequency residual conductivity observed determined by tunnelling spectroscopy (STM) for the σ and π even at 4 K [Fig. 1(e)], the temperature dependence of the bands (squares and circles) [28], and as estimated from Fig. 1(c) superconducting gap 2∆π(T) cannot be straightforwardly de- (spheres). Leggett mode versus temperature as observed in Raman rived, e.g. by fitting the frequency-dependent optical conduc- spectroscopy (left-triangles) [41, 42] and obtained from theoretical tivity according to the Mattis-Bardeen theory [35]. To ex- calculations (right-triangles) [43]. The horizontal solid lines indicate 2 f for f = 0.3, 0.4, 0.5, 0.6, and 0.7 THz, respectively. (c) The tract qualitatively the temperature dependence of the gap in corresponding normalized THG intensity I3 f as a function of tem- the π−band, we approximate 2∆π(T) by the point where the perature, measured through two band-pass filters of 3 f . The vertical normalised transmission crosses the unity. As indicated by the dashed lines indicate the peak positions in the temperature-dependent arrows in Fig. 1(c), the crossing point shifts continuously to- curves, where the corresponding resonance condition 2 f = 2∆π wards lower frequency with increasing temperature. By mul- is fulfilled simultaneously. (d) The normalized third-order suscep- (3) 2 4 tiplying the extracted crossing point frequencies with the ratio tibility |χ | /∆π versus normalized frequency 2 f /2∆π(T) for f = of 2∆ = 1.03 THz and the crossing point frequency at 4 K 0.3 THz at various temperatures below Tc. Only one band-pass filter π of 3 f was utilized. The temperature dependent gap 2∆ (T) is param- ∆ T π we obtain 2 π( ) as shown in Fig. 2(c) by the solid symbols. eterized using the experimental values as shown by the dash-dotted The result agrees very well with the tunneling data [28]. line in (b). The dashed lines provide a guide to the eye. Third-harmonic radiation is observed by pumping the sys- tem with multicycle THz pulses. Figure 2(a) displays the time-domain waveform of the f = 0.5 THz pump pulse. Driven by this pulse, the generated radiation of the sample nor for the Leggett mode that was previously determined from is recorded after a 3 f = 1.5 THz bandpass filter in time do- Raman response [41–43], despite that such resonance peaks main [see Fig. 3(a) for a illustrative sketch of the experimental were predicted in the clean limit [14, 15]. setup]. Starting from 5 K, the observed THG intensity in- Appearance of a maximum in an experimental I3 f (T) curve creases continuously with increasing temperature until about can also be coincident, because while the superconductinggap 19 K [Fig. 2(c)], which is followed by a continuous decrease increases with decreasing temperature, the effective driving at higher temperatures towards Tc. Above Tc, the THG sig- field inside the thin film E f (T) decreases due to the increasing nal is dramatically reduced. The enhancement of the third- of the reflectivity. To rule out such a coincidence, we per- harmonic generation below Tc clearly reflects the nonlinear formed THG measurements for various pump-pulse frequen- response of the superconducting state. Moreover, the intensity cies, i.e. f = 0.3, 0.4, 0.6, and 0.7 THz. As presented in maximum for 3 f = 1.5 THz occurs right at a resonance condi- Fig. 2(c), for f = 0.3 and 0.4 THz, i.e. 2 f < 2∆π(T = 0), tion where 2 f = 2∆π, as indicated by the vertical dashed line the temperature-dependent THG curves exhibit maxima at in Figs. 2(b) and 2(c). However, such resonant behaviour is 30 and 28 K, respectively, matching the resonance condi- neither observed in the vicinity of temperatures where the res- tion 2 f = 2∆π (indicated by the arrows). Together with onance condition would be fulfilled for the gap in the σ-band, the 0.5 THz results, the systematic shifting of the maximum 4

matching the resonance condition [Fig. 2(b)(c)] strongly evi- (a) q q q

2 3 1 MgB /MgO dence the driven Higgs oscillations of the π band. 2 Based on the Anderson’s pseudospin formalism [5, 6], it was predicted for a single-band superconductor that the THG THz intensity I3 f due to the THz driven Higgs oscillations is pro-

6 4 −1 c portional to E f ∆ |2 f − 2∆| , where E f is the amplitude of

BPF- f the driving THz electric field and ∆(T) is the superconducting W GP1 W GP2 W GP3 BPF-3 f

(b) gap. This means that the normalized third-order susceptibil- (c) (3) 2 4 −1 ity |χ | /∆ ∝ |2 f − 2∆| . Hence, in Fig. 2(d) we present f = 0.4 THz

3f = 1.2 THz, 20 K | (3)|2 ∆4 ∆ 40 χ / π versus normalized frequency 2 f /2 π(T) in MgB2 at 3f = 1.5 THz, 4.6 K (3) various temperatures below Tc for f = 0.3 THz. The χ (T) 2 20 is estimated by a simultaneous measurement of the THG and

the transmission of the fundamental frequency through one 0 band-pass filter of 3 f . Thegap2∆π(T) is parameterized using 1 the experimental values as shown by the dash-dotted line in -20 Electric field (kV/cm) field Electric Fig. 2(b). As shown in Fig. 2(d), a peak appears very close to (norm.) intensity THG

-40 the position of resonance 2 f /2∆π(T) = 1 in good agreement with the theoretical prediction of the Higgs amplitude mode. 0

0 10 20 0 30 60 90 For f = 0.6 and 0.7 THz, the THG signal decreases mono- Time delay (ps) Polarization q ( deg. ) tonically with increasing temperature. While no resonance is 2 expected for the gap in the π-band as 2 f > 2∆π, there is nei- ther peak-like enhancement at temperatures where 2 f = 2∆ . Figure 3. (a) Schematic illustration of the setup for polarization de- σ pendent THG measurements. The THz electromagnetic waves prop- These observationssuggest a much weaker THG for the Higgs agate along the sample c-axis. BPF-f and BPF-3f denote band- mode of the gap in the σ-band that is in the clean-limit. In- pass filters for the fundamental frequency and the THG, respec- stead, the THG seems to be dominated by the coupling to the tively. WGP’s denote wird-grid polarizers. (b) Waveform of the −1 Higgs mode of the π-band, where τ > 2∆π in the dirty limit. f = 0.4 THz pump pulse for the polarization-dependent measure- Together, these observations indicate that the observed THG ments. (c) Normalized THG intensity as a function of the pump- pulse polarization, varied by tuning the angle θ of WGP2, which is signals for 2 f < 2∆π(T = 0) is mainly due to the Higgs mode 2 of the smaller gap in the π−band rather than charge density obtained for f = 0.4 and 0.5 THz at 20 and 4.6 K, respectively. fluctuations, in line with recent theoretical studies [17]. Further support for the THG being a result of coupling to the Higgs mode can be provided by polarization dependent sample is presented in Fig. 3(c) for 3 f = 1.2 and 1.5 THz at study. Theoretical studies on s- and d-wave superconductors 20 and 4.6 K, respectively. The peak fields of the correspond- showed that THG from the fully symmetric A1g Higgs am- ing pump pulses were set at 47 and 54 kV/cm for f = 0.4 and plitude mode is isotropic, while the contribution of charge 0.5 THz, respectively. These experimental settings were cho- density fluctuations is strongly dependent on the pump-pulse sen in a way that 2 f < 2∆π and the THG response is in the polarization [7, 10, 15, 16]. It is worth noting that the single- perturbative regime. In the presented THG intensity, a fac- 2 crystalline thin film, as used in this experiment, is crucial for tor, cos (θ2 − θ3), due to the polarization difference between such study. To discriminate the contributions of charge den- WGP2 and WGP3 has been taken into account. Figure 3(c) sity fluctuations and the A1g Higgs mode to the observed THG clearly shows that the THG intensity is nearly independent signals, we performed polarization-dependent measurements. on the pump-pulse polarization for both the f = 0.4 and the As illustrated in Fig. 3(a), the polarization of the pump pulse 0.5 THz pump pulses. As the observed THG is dominated is tuned by the wire-grid polarizer (WGP2) directly in front by the response of the π-band, the observed independence on of the sample. To assure a constant pump fluence at the sam- the polarization provides further evidence for the Higgs mode ple position, another WGP1 is placed between WGP2 and the being responsible for the THG. bandpass filter (BPF-f ) of the fundamental frequency. After To conclude, linear and nonlinear response of the super- the sample, an additional wire-grid polarizer (WGP3) is in- conducting state in the two-band superconductor MgB2 has serted with a fixed polarization direction at 45◦ with respect been studied by time-resolved THz spectroscopy. The third- to the horizontal direction. With this polarizer, the electric- harmonic generation is observed to be resonantly enhanced optics sampling at the ZnTe crystal (not shown) was optimized when the lower-energygap 2∆π matches twice the pump-pulse for detection. Before measuring the polarization dependence frequency. Such resonant behaviour is absent when such a of the THG signals, the setup was calibrated by adjusting θ1 condition is met for the higher-energy gap in the σ−band and while setting θ2 fordifferent polarizations, to ensure a constant for the Leggett phase mode. Intensity of the third-harmonic pump fluence at the sample position. radiation is found nearly isotropic with varying the pump- Figure 3(b) shows waveform of the pump pulse of f = pulse polarization. Together, these results provide evidence 0.4 THz. The polarization-dependent THG intensity from the for the nonlinear response of the π-band Higgs amplitude 5 mode, in agreement with theoretical predictions [17]. Com- [21] C. Zhang, Y. Wang, D. Wang, Y. Zhang, Z. H. Liu, Q. R. Feng, plementary pump-probe approaches [43, 44] are expected to and Z. Z. Gan, . Appl. Phys. 114, 023903 (2013). address further aspects of the nonlinear response of the multi- [22] C. G. Zhuang, S. Meng, C. Y. Zhang, Q. R. Feng, Z. Z. band superconductors [4]. Gan, H. Yang, Y. Jia, H. H. Wen, and X. X. Xi, ZW gratefully acknowledges G. Blumberg, B. B. Jin, J. Appl. Phys. 104, 013924 (2008). [23] L. Y. Shi, Y. Q. Liu, T. Lin, M. Y. Zhang, S. J. I. I. Mazin, R. Shimano, and N. Tsuji for illuminating dis- Zhang, L. Wang, Y. G. Shi, T. Dong, and N. L. Wang, cussions and/or communications. T.D. acknowledges finan- Phys. Rev. B 98, 094414 (2018). cial support by the Alexander von Humboldt foundation. The [24] K.-L. Yeh, M. C. Hoffmann, J. Hebling, and K. A. Nelson, work in Mainz was partially supported by the DFG (Ger- Appl. Phys. Lett. 90, 171121 (2007). man Research Foundation) in the framework of the Collabo- [25] S. Kovalev, R. M. A. Dantas, S. Germanskiy, J.-C. Deinert, rative Research Centre SFB TRR 173, “Spin + X”. The work B. Green, I. Ilyakov, N. Awari, M. Chen, M. Bawatna, J. Ling, in Beijing was supported by the National Science Founda- F. Xiu, P. H. M. van Loosdrecht, P. Sur´owka, T. Oka, and Z. Wang, Nature Commun. 11, 2451 (2020). tion of China (No.11888101) and the National Key Research [26] J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, and and Development Program of China (No.2016YFA0300902, J. Akimitsu, Nature 410, 63 (2001). 2017YFA0302904). The work in Cologne was partially sup- [27] J. Kortus, I. I. Mazin, K. D. Belashchenko, V. P. Antropov, and ported by the DFG via Project No. 277146847 — Collabora- L. L. Boyer, Phys. Rev. Lett. 86, 4656 (2001). tive Research Center 1238: Control and Dynamics of Quan- [28] M. Iavarone, G. Karapetrov, A. E. Koshelev, W. K. Kwok, G. W. tum Materials (Subproject No. B05). Crabtree, D. G. Hinks, W. N. Kang, E.-M. Choi, H. J. Kim, H.- S.K. and T.D. contributed equally to this work. J. Kim, and S. I. Lee, Phys. Rev. Lett. 89, 187002 (2002). [29] S. Souma, Y. Machida, T. Sato, T. Takahashi, H. Matsui, S. C. Wang, H. Ding, A. Kaminski, J. C. Campuzano, S. Sasaki, and K. Kadowaki, Nature 423, 65 (2003). [30] P. C. Canfield and G. W. Crabtree, Physics Today 56, 34 (2003). ∗ [email protected] [31] X. X. Xi, Rep. Prog. Phys. 71, 116501 (2008). † [email protected] [32] D. Vollhardt, Phys. Rev. Lett. 78, 1307 (1997). [1] P. W. Anderson, Phys. Rev. 112, 1900 (1958). [33] M. Greger, M. Kollar, and D. Vollhardt, [2] P. W. Anderson, Phys. Rev. 130, 439 (1963). Phys. Rev. B 87, 195140 (2013). [3] D. Pekker and C. M. Varma, [34] Z. Wang, M. Schmidt, J. Fischer, V. Tsurkan, Annu. Rev. Condens. Matter Phys. 6, 269 (2015). M. Greger, D. Vollhardt, A. Loidl, and J. Deisenhofer, [4] R. Shimano and N. Tsuji, Nat. Commun. 5, 3202 (2014). Annu. Rev. Condens. Matter Phys. 11, 103 (2020). [35] R. A. Kaindl, M. A. Carnahan, J. Orenstein, D. S. Chemla, [5] R. Matsunaga, N. Tsuji, H. Fujita, A. Sugioka, K. Makise, H. M. Christen, H.-Y. Zhai, M. Paranthaman, and D. H. Lown- Y. Uzawa, H. Terai, Z. Wang, H. Aoki, and R. Shimano, des, Phys. Rev. Lett. 88, 027003 (2002). Science 345, 1145 (2014). [36] J. J. Tu, G. L. Carr, V. Perebeinos, C. C. Homes, M. Strongin, [6] N. Tsuji and H. Aoki, Phys. Rev. B 92, 064508 (2015). P. B. Allen, W. N. Kang, E.-M. Choi, H.-J. Kim, and S.-I. Lee, [7] R. Matsunaga, N. Tsuji, K. Makise, H. Terai, H. Aoki, and Phys. Rev. Lett. 87, 277001 (2001). R. Shimano, Phys. Rev. B 96, 020505 (2017). [37] B. B. Jin, P. Kuzel, F. Kadlec, T. Dahm, J. M. Red- [8] X. Yang, C. Vaswani, C. Sundahl, M. Mootz, L. Luo, wing, A. V. Pogrebnyakov, X. X. Xi, and N. Klein, J. H. Kang, I. E. Perakis, C. B. Eom, and J. Wang, Appl. Phys. Lett. 87, 092503 (2005). Nature Photonics 13, 707 (2019). [38] A. V. Pronin, A. Pimenov, A. Loidl, and S. I. Krasnosvobodt- [9] H. Chu et al, Nature Commun. 11, 1793 (2020). sev, Phys. Rev. Lett. 87, 097003 (2001). [10] L. Schwarz and D. Manske, Phys. Rev. B 101, 184519 (2020). [39] V. Guritanu, A. B. Kuzmenko, D. van der Marel, [11] T. Cea, C. Castellani, and L. Benfatto, S. M. Kazakov, N. D. Zhigadlo, and J. Karpinski, Phys. Rev. B 93, 180507 (2016). Phys. Rev. B 73, 104509 (2006). [12] N. Tsuji, Y. Murakami, and H. Aoki, [40] M. Ortolani, P. Dore, D. Di Castro, A. Perucchi, S. Lupi, V. Fer- Phys. Rev. B 94, 224519 (2016). rando, M. Putti, I. Pallecchi, C. Ferdeghini, and X. X. Xi, [13] H. Krull, N. Bittner, G. S. Uhrig, D. Manske, and A. P. Schny- Phys. Rev. B 77, 100507 (2008). der, Nature Communications 7, 11921 (2016). [41] G. Blumberg, A. Mialitsin, B. S. Dennis, M. V. [14] T. Cea and L. Benfatto, Phys. Rev. B 94, 064512 (2016). Klein, N. D. Zhigadlo, and J. Karpinski, [15] Y. Murotani, N. Tsuji, and H. Aoki, Phys. Rev. Lett. 99, 227002 (2007). Phys. Rev. B 95, 104503 (2017). [42] G. Blumberg, A. Mialitsin, B. S. Dennis, N. D. Zhigadlo, and [16] T. Cea, P. Barone, C. Castellani, and L. Benfatto, J. Karpinski, Physica C 46, 75 (2007). Phys. Rev. B 97, 094516 (2018). [43] F. Giorgianni, T. Cea, C. Vicario, C. P. Hauri, W. K. Withanage, [17] Y. Murotani and R. Shimano, Phys. Rev. B 99, 224510 (2019). X. Xi, and L. Benfatto, Nature Physics 15, 341 (2019). [18] T. Yu and M. W. Wu, Phys. Rev. B 96, 155311 (2017). [44] K. Katsumi, N. Tsuji, Y. I. Hamada, R. Matsunaga, J. Schnee- [19] M. Silaev, Phys. Rev. B 99, 224511 (2019). loch, R. D. Zhong, G. D. Gu, H. Aoki, Y. Gallais, and R. Shi- [20] C. Zhang, Y. Wang, D. Wang, Y. Zhang, mano, Phys. Rev. Lett. 120, 117001 (2018). Q. R. Feng, and Z. Z. Gan, IEEE Transactions on Applied Superconductivity 23, 7500204 (2013).