Arxiv:2010.05019V1 [Cond-Mat.Supr-Con] 10 Oct 2020 Clrectto Ihu Hreo Electric a Or Is Charge fluctuations Without the Amplitude Excitation Gap
Total Page:16
File Type:pdf, Size:1020Kb
Band-Selective Third-Harmonic Generation in Superconducting MgB2: Evidence for Higgs Amplitude Mode in the Dirty Limit Sergey Kovalev,1 Tao Dong,2,3, ∗ Li-Yu Shi,3 Chris Reinhoffer,4 Tie-Quan Xu,5 Hong-Zhang Wang,5 Yue Wang,5 Zi-Zhao Gan,5 Semyon Germanskiy,4 Jan-Christoph Deinert,1 Igor Ilyakov,1 Paul H. M. van Loosdrecht,4 Dong Wu,3 Nan-Lin Wang,3, 6 Jure Demsar,2 and Zhe Wang7, 4, 1, † 1Institute of Radiation Physics, Helmholtz-Zentrum Dresden-Rossendorf, 01328 Dresden, Germany 2Institute of Physics, Johannes Gutenberg-University Mainz, 55128 Mainz, Germany 3International Center for Quantum Materials, School of Physics, Peking University, Beijing 100871, China 4Institute of Physics II, University of Cologne, 50937 Cologne, Germany 5Applied Superconductivity Center and State Key Laboratory for Mesoscopic Physics, School of Physics, Peking University, Beijing 100871, China 6Collaborative Innovation Center of Quantum Matter, Beijing, China 7Fakult¨at Physik, Technische Universit¨at Dortmund, 44221 Dortmund, Germany (Dated: October 13, 2020) We report on time-resolved linear and nonlinear terahertz spectroscopy of the two-band superconductor MgB2 with the superconducting transition temperature Tc ≈ 36 K. Third-harmonic generation (THG) is observed below Tc by driving the system with intense narrowband THz pulses. For the pump-pulse frequencies f = 0.3, 0.4, and 0.5 THz, temperature-dependent evolution of the THG signals exhibits a resonance maximum at the temperature where 2 f = 2∆π(T), for the dirty-limit superconducting gap 2∆π = 1.03 THz at 4 K. In contrast, for f = 0.6 and 0.7 THz with 2 f > 2∆π, the THG intensity increases monotonically with decreasing temperature. Moreover, for 2 f < 2∆π the THG is found nearly isotropic with respect to the pump-pulse polarization. These results suggest the predominant contribution of the driven Higgs amplitude mode of the dirty-limit gap, pointing to the importance of scattering for observation of the Higgs mode in superconductors. Novel quantum phenomena, which may not be accessible [10, 12]. Especially when correlation effects and/or impurity to the linear spectroscopy in equilibrium states, can often be scattering are pronounced, the contribution due to the driven visualized in driven non-equilibrium states. A prominent ex- Higgs mode will dominate the THG signal. Moreover,the two ample is the Higgs amplitude mode in a superconductor [1– competing scenarios suggest different dependences on the po- 4]. With the U(1) rotational symmetry spontaneously bro- larization of the driving THz field with respect to crystalline ken, a superconducting state is characterized by fluctuations axis. NbN was found to exhibit an isotropic THG signal with of the amplitude and phase of the superconducting gap. The respect to the polarization of the driving field, a signature sup- Higgs mode corresponding to the amplitude fluctuations is a porting the assignment of the THG to the driven Higgs mode, scalar excitation without charge or electric/magnetic dipoles. whereas the THG from charge density fluctuations is expected As such, it cannot directly couple to external probes in the to be anisotropic [7]. linear-response regime [4]. Recently, extensive theoretical studies of the THG in two- Terahertz (THz) harmonic generation has been suggested band superconductors addressed not only the Higgs ampli- tobeanefficient probe to reveal Higgs mode in superconduc- tude fluctuations in individual bands, but also the so-called tors [5–10]. Here, an intense narrowband THz pulse can drive Leggett phase mode, resulting from the inter-band coupling a superconductor out of equilibrium, and trigger the coher- [13–17]. Here, fluctuations of relative phase between the two- ent Higgs oscillations of the superconducting order parame- coupled order parameters, the Leggett mode, is also charge ter. The driven Higgs mode with twice the frequency of the neutral and does not linearly couple to the electromagnetic pump pulse is further coupled to the pump pulse, leading to field [13]. Similar to the Higgs modes, the Leggett mode with the third-harmonic generation (THG) [4]. Experimentally it frequency ωL can be resonantly excited by a multicycle THz = arXiv:2010.05019v1 [cond-mat.supr-con] 10 Oct 2020 is, however, not straightforward to unambiguously attribute pulse with frequency f satisfying the condition of 2 f ωL, the THG signal to the driven Higgs mode. In addition to a which is expected to be manifested by a resonance peak in major enhancement of the THG signal below the supercon- the THG [14–17]. Moreover, paramagnetic impurity scatter- ducting transition temperature Tc, the Higgs mode excitation ing has been predicted to play a decisive role in the nonlinear should give rise to a peak in the THG signal when twice the response of the superconductors [17–19]. While in the clean pump-pulse frequency matches the superconducting gap fre- limit the THG was found to be dominated by pair breaking quency. Although these features were observed in the s-wave [11, 15], the contribution of the Higgs mode can be signif- superconductor NbN [5], it was argued that the THG signal icantly enhanced by impurity scattering. Thus, in the dirty is more likely a result of charge density fluctuations due to limit, the Higgs amplitude fluctuations are predominant than the THz driven Cooper-pair breaking [11]. The contribution the charge density fluctuations, while the Leggett phase fluc- of charge density fluctuations to the THG is, however, not al- tuations in a two-band superconductor are much weaker and ways dominant, as pointed out by further theoretical studies hardly detectable in a THG experiment [17]. 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 de- 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].