Abstract
Photo-excitation at terahertz and mid-infrared frequencies has emerged as an effective way to manipulate functionalities in quantum materials, in some cases creating non-equilibrium phases that have no equilibrium analogue. In K3C60, a metastable zero-resistance phase was observed that has optical properties, nonlinear electrical transport and pressure dependencies compatible with non-equilibrium high-temperature superconductivity. Here we demonstrate a two-orders-of-magnitude increase in photo-susceptibility near 10 THz excitation frequency. At these drive frequencies, a metastable superconducting-like phase is observed up to room temperature. The discovery of a dominant frequency scale sheds light on the microscopic mechanism underlying photo-induced superconductivity. It also indicates a path towards steady-state operation, limited at present by the availability of a suitable high-repetition-rate optical source at these frequencies.
Similar content being viewed by others
Main
The search for new non-equilibrium functional phases in quantum materials, such as optically induced ferroelectricity1,2, magnetism3,4,5, charge density wave order6,7, non-trivial topology8,9 and superconductivity10,11,12,13,14,15,16,17,18, has become a central research theme in condensed-matter physics. In the case of K3C60 (Fig. 1a), mid-infrared optical pulses have been extensively documented to yield an unconventional non-equilibrium phase that exhibits metastable zero resistance14, extraordinarily high mobility and a superconducting-like gap in the optical conductivity12,14 that reduce with applied pressure13, and nonlinear current–voltage characteristics19. All these observations are indicative of non-equilibrium high-temperature superconductivity, observed at base temperatures far exceeding the highest equilibrium superconducting critical temperature of any alkali-doped fulleride (Fig. 1b).
Typical experimental results reported until now are displayed in Fig. 2c. K3C60 powders were held at a base temperature T = 100 K ≫ Tc = 20 K and irradiated with 1-ps-long pulses with 170 meV photon energy (wavelength λ ≈ 7.3 µm, frequency f ≈ 41 THz) at a fluence of 18 mJ cm−². This strong excitation regime yielded a long-lived transient state with pronounced changes in both the real and imaginary parts of the optical conductivity, measured using phase-sensitive terahertz time-domain spectroscopy. The transient optical properties displayed in Fig. 2c are reminiscent of those of the equilibrium superconducting state measured in the same material at T ≪ Tc = 20 K (cf. Fig. 2b) and are suggestive of transient high-temperature superconductivity. These signatures consist of perfect reflectivity, a gap in the real part of the optical conductivity \({\sigma }_{1}\left(\omega \right)\) and an imaginary conductivity \({\sigma }_{2}\left(\omega \right)\) that diverges towards low frequencies as \({\sigma }_{2}\left(\omega \right)\propto 1/\omega\). The divergent \({\sigma }_{2}(\omega )\) implies (through Kramers–Kronig relations) the presence of a peak in \({\sigma }_{1}\) centred at zero frequency, with a width limited by the lifetime of the state that also determines the carrier mobility.
These data were obtained by accounting for the inhomogeneous excitation of the probed volume using a multilayer model. Here we show the results of this reconstruction under the assumption of a linear (filled blue circles) and sublinear (open symbols)20 dependence of the photo-induced changes in the terahertz refractive index on the mid-infrared pump fluence. Two sublinear models are shown, with an assumed square-root fluence dependence in blue circles and saturating fluence dependence in brown triangles, as detailed in Supplementary Section 6. Allowing for a finite-temperature superconductor in which a varying density of uncondensed quasiparticles also contributes to the terahertz response, the superconducting-like nature of the transient state is independent of the specific choice of assumption21. Only quantitative differences associated with the relative densities of induced superfluid and heated quasiparticles, which can be extracted by fitting with a two-fluid model, emerge. This modelling also reveals that the enhancement of conductivity observed in these experiments is connected not to an increase in the carrier density, which remains constant upon photo-excitation, but rather to a giant increase in the carrier mobility.
Three spectrally integrated figures of merit are extracted from the snapshots of reflectivity (R\(\left(\omega ,\tau \right)\)), \({\sigma }_{1}\left(\omega ,\tau \right)\) and \({\sigma }_{2}\left(\omega ,\tau \right)\) and plotted as a function of pump-probe time delay \(\tau\) in Fig. 2e, showing the time evolution of the system. The first two quantities are the frequency-averaged values of R(ω) and \({\sigma }_{1}(\omega )\) from 5 to 10 meV, a frequency range that lies below the photo-induced energy gap, for which a zero-temperature superconductor with infinite lifetime would give values of 1 and 0, respectively. The third figure of merit is the fractional superfluid density, which is proportional to the divergence of \({\sigma }_{2}(\omega )\). This is determined by fitting the photo-excited optical properties with a two-fluid model where one fluid represents the remaining normal carriers with a finite scattering rate and the other has zero scattering rate, giving a superconducting-like contribution. Details of this fitting procedure are given in Supplementary Section 7.
For low excitation fluences, the system becomes superconducting-like after photo-excitation and relaxes on a time scale of a few picoseconds. As already seen in the spectrally resolved measurements, for high excitation fluences, the system enters a metastable regime in which the superconducting-like optical properties persist for much longer times, up to several nanoseconds.
We note that the temperature dependence reported in ref. 12 shows transient superconducting-like optical properties up to a temperature of 150–200 K. For higher temperatures, the gapping and extracted superfluid density are severely reduced. Examples of such spectra measured at room temperature are shown in Fig. 2d. Nevertheless, the pressure scaling reported in ref. 13 suggests that traces of non-equilibrium superconductivity may survive up to higher temperatures, raising the prospect that with more effective driving, a full manifestation of the metastable superconducting-like state may be possible at 300 K.
Until now, these experiments have been limited to excitation photon energies between 80 and 165 meV (20–40 THz), such that a more comprehensive search for a dominant excitation frequency scale has remained out of reach. Many potentially important resonances at lower frequencies (hν < 80 meV) have remained unexplored, primarily due to the lack of a suitable high-intensity pump source that operates in this range. In the present work, we explore excitation at energies between 24 and 80 meV (6–20 THz). This energy range hosts many excitations, both vibrational (phonons) and electronic in nature, including a broad polaronic peak seen in \({\sigma }_{1}\) centred at about 60 meV (15 THz). The possible relevance of this excitation has been highlighted in ref. 22, although this prediction could not be tested to date.
To achieve wide tuneability, we used a terahertz source based on chirped-pulse difference frequency generation, mixing the near-infrared signal beams of two phase-locked optical parametric amplifiers23. This source, illustrated schematically in Fig. 3a and described in detail in Supplementary Section 4, was used to generate narrow-bandwidth pulses with photon energies spanning the range from 24 to 145 meV (6–35 THz). All measurements reported here were carried out with an excitation bandwidth of hΔνpump ≈ 4 meV (1 THz) and \(\Delta {{\tau }}_{{\rm{pump}}}\approx 600\) fs pulse duration. The same probing protocol as that reported in Fig. 2 was used here to detect changes in the complex optical properties for probe energies spanning 4–18 meV (1–4.5 THz).
Figure 3b–d shows reflectivity and complex conductivity spectra measured after photo-excitation with pulses tuned to 41 meV photon energy (λ ≈ 30 µm, f ≈ 10 THz) at base temperatures of 100 K and room temperature, respectively. The response is very similar to the case reported in Fig. 2 for 170 meV (41 THz) excitation, manifested on metastable timescales but persisting here up to room temperature despite an excitation fluence almost two orders of magnitude weaker. Figure 3e displays the time evolution of the optical properties with a transient amplifying state at early pump-probe delays, as already reported in ref. 24, which relaxes into a metastable state with superconducting-like optical properties.
Figure 4a shows the scaling with fluence of \(R(\omega )\) and \({\sigma }_{1}(\omega )\) (averaged in the 5–10 meV range), as well as the fractional superfluid density in response to photo-excitation at 170 meV (41 THz) and 41 meV (10 THz). These measurements were carried out at a pump-probe time delay of 10 ps and thus refer to the metastable phase. The figure shows how all figures of merit approach their equilibrium superconducting-state values as the fluence increases, with the fluence required being roughly 50 times less for 41 meV (10 THz) compared to 170 meV (41 THz) excitation.
Similar fluence dependence measurements were carried out by varying the photon energy of the pump and maintaining a constant 4 meV (1 THz) bandwidth with 600 fs pulse duration. For all excitation photon energies between 24 meV (6 THz) and 145 meV (35 THz), the photo-induced changes in the optical properties were qualitatively similar to those shown in Figs. 2 and 3, with only the size of the response for a given fluence differing. From each fluence dependence, we extracted a figure of merit for the photo-susceptibility, defined as the rate of depletion of spectral weight in \({\sigma }_{1}\) with excitation fluence, in the limit of low fluence (Supplementary Section 8). Plots of the pump-frequency-dependent photo-susceptibility are shown in Fig. 4b for both 10 ps and 50 ps pump-probe time delay. A peak centred at 41 meV (10 THz) with about 16 meV FWHM width is observed in these measurements. Next, we will discuss three distinct energy scales that coincide with this resonance; sequentially, they relate to ‘on-ball’ orbital excitations, phonons and excitons.
Superconductivity in alkali-doped fullerides is believed to be mediated by a dynamical Jahn–Teller distortion, which leads to an effective negative Hund’s coupling for the orbitals of a single buckyball25 and to a low-spin S = 1/2 state. A theoretical model based on this assumption has been successful at providing a quantitatively correct phase diagram for fulleride superconductors from ab initio calculations26,27. In this model, the local ground state of the system is a sixfold degenerate low-spin state, which features intraorbital pairs that delocalize over two molecular orbitals. As detailed in Supplementary Section 10, a first set of local excited states also features such pairs, albeit with a different angular momentum (that is, a different interorbital phase for the delocalized pair). Ab initio calculations predict an energy splitting of 37 meV between these two sets of states26,27. The observed resonance may therefore be related to the creation of interorbital pairs with local angular momentum, which may also contribute to superconductivity, as suggested in the Suhl–Kondo mechanism28,29,30. However, it is not yet clear exactly how this excitation is transformed in the presence of tunnelling between neighbouring C60 molecules and why the creation of such pairs may support metastable superconductivity at such high temperatures. Furthermore, as the local parity of this excited state would be different from that of the ground state, condensation in this configuration may give rise to a superconductor with different symmetry. This possibility, whilst tantalizing, remains speculative and should be tested with more comprehensive ultrafast probing methods.
Turning to phonon excitations, we also note that the 41 meV resonance frequency identified here coincides with an infrared-active T1u phonon, which predominantly consists of intramolecular motion of the C atoms. Whereas the atomic motions of the 170 meV molecular mode discussed previously in ref. 12 are directed along the tangential directions of the C603− molecule, those of the 41 meV mode are predominantly along the radial directions (Supplementary Section 9). By performing frozen-phonon calculations using density functional theory (DFT), we evaluated the different effect of these distortions on the three t1u molecular levels at the Fermi energy, which we map out from DFT wave functions as maximally localized Wannier functions (Supplementary Section 9). In the undistorted C603− structure, these molecular levels are degenerate. Applying a distortion along a T1u coordinate lifts this degeneracy, leaving a doubly degenerate t1u orbital lowered in energy. This electronic configuration is prone to developing a Jahn–Teller distortion that may lead to an enhanced negative Hund’s coupling, possibly facilitating the onset of superconductivity at higher temperatures. The strength of the induced splitting is quadratic in the phonon coordinate and is enhanced when driving the 41 meV mode compared to the 170 meV one, suggesting that the observed resonance may arise from a more efficient manipulation of the electronic degrees of freedom when driving the 41 meV T1u mode.
Finally, we address the electronic excitations discussed in ref. 22, where an enhanced effect on tuning the drive to lower frequencies was already predicted. This work proposed the existence of a dressed exciton in the excitation spectrum of K3C60 at the same energy scale as the resonance reported here. Excitation of the system at this frequency would generate excitons that would act as a reservoir able to incoherently cool the quasiparticles, resulting in the reemergence of superconductivity at base temperatures in excess of the equilibrium Tc.
We expect the importance of this discovery to be capitalized on in future work. The extreme efficiency improvement due to resonant enhancement, nearing two orders of magnitude, is expected to also dramatically reduce unwanted dissipation. This, taken in conjunction with the observed nanosecond-long lifetime, suggests that excitation of the sample with a train of pulses of only 400 μJ cm−2 delivered at a 100 MHz repetition rate—as determined by the inverse lifetime of this state—may yield continuous wave operation. Because this effect is documented here to persist up to room temperature, continuous wave operation would likely have important practical implications. To make this regime experimentally accessible, single-order-of-magnitude improvements in the efficiency of the process or the light–matter coupling strength, combined with suitable developments in high-repetition-rate terahertz sources, would be required.
Methods
All methods can be found in the Supplementary Information.
Data availability
Source data are provided with this paper. All other data that support the plots in this paper and other findings of this study are available from the corresponding authors upon reasonable request.
References
Nova, T. F., Disa, A. S., Fechner, M. & Cavalleri, A. Metastable ferroelectricity in optically strained SrTiO3. Science 364, 1075–1079 (2019).
Li, X. et al. Terahertz field–induced ferroelectricity in quantum paraelectric SrTiO3. Science 364, 1079–1082 (2019).
Disa, A. S. et al. Polarizing an antiferromagnet by optical engineering of the crystal field. Nat. Phys. 16, 937–941 (2020).
Wang, X. et al. Light-induced ferromagnetism in moiré superlattices. Nature 604, 468–473 (2022).
Radu, I. et al. Transient ferromagnetic-like state mediating ultrafast reversal of antiferromagnetically coupled spins. Nature 472, 205–208 (2011).
Kogar, A. et al. Light-induced charge density wave in LaTe3. Nat. Phys. 16, 159–163 (2020).
Wandel, S. et al. Enhanced charge density wave coherence in a light-quenched, high-temperature superconductor. Science 376, 860–864 (2022).
Wang, Y. H., Steinberg, H., Jarillo-Herrero, P. & Gedik, N. Observation of Floquet-Bloch states on the surface of a topological insulator. Science 342, 453 (2013).
McIver, J. W. et al. Light-induced anomalous Hall effect in graphene. Nat. Phys. 16, 38–41 (2020).
Fausti, D. et al. Light-induced superconductivity in a stripe-ordered cuprate. Science 331, 189–191 (2011).
Hu, W. et al. Optically enhanced coherent transport in YBa2Cu3O6.5 by ultrafast redistribution of interlayer coupling. Nat. Mater. 13, 705–711 (2014).
Mitrano, M. et al. Possible light-induced superconductivity in K3C60 at high temperature. Nature 530, 461–464 (2016).
Cantaluppi, A. et al. Pressure tuning of light-induced superconductivity in K3C60. Nat. Phys. 14, 837–841 (2018).
Budden, M. et al. Evidence for metastable photo-induced superconductivity in K3C60. Nat. Phys. 17, 611–618 (2021).
Buzzi, M. et al. Photo-molecular high temperature superconductivity. Phys. Rev. 10, 031028 (2020).
Buzzi, M. et al. Phase diagram for light-induced superconductivity in κ-(ET)2-X. Phys. Rev. Lett. 127, 197002 (2021).
Cremin, K. A. et al. Photoenhanced metastable c-axis electrodynamics in stripe-ordered cuprate La1.885Ba0.115CuO4. Proc. Natl Acad. Sci. 116, 19875 (2019).
Isoyama, K. et al. Light-induced enhancement of superconductivity in iron-based superconductor FeSe0.5Te0.5. Commun. Phys. 4, 160 (2021).
Wang, E. et al. Nonlinear transport in a photo-induced superconductor. Preprint at https://arxiv.org/abs/2301.06425 (2023).
Dodge, J. S., Lopez, L. & Sahota, D. G. Optical saturation produces spurious evidence for photoinduced superconductivity in K3C60. Phys. Rev. Lett. 130, 146002 (2023).
Buzzi, M., Nicoletti, D., Rowe, E., Wang, E. & Cavalleri, A. Comment on arXiv:2210.01114: Optical saturation produces spurious evidence for photoinduced superconductivity in K3C60. Preprint at https://arxiv.org/abs/2303.10169 (2023).
Nava, A., Giannetti, C., Georges, A., Tosatti, E. & Fabrizio, M. Cooling quasiparticles in A3C60 fullerides by excitonic mid-infrared absorption. Nat. Phys. 14, 154–159 (2018).
Liu, B. et al. Generation of narrowband, high-intensity, carrier-envelope phase-stable pulses tunable between 4 and 18 THz. Opt. Lett. 42, 129–131 (2017).
Buzzi, M. et al. Higgs-mediated optical amplification in a nonequilibrium superconductor. Phys. Rev. 11, 011055 (2021).
Capone, M., Fabrizio, M., Castellani, C. & Tosatti, E. Strongly sorrelated superconductivity. Science 296, 2364–2366 (2002).
Nomura, Y., Sakai, S., Capone, M. & Arita, R. Unified understanding of superconductivity and Mott transition in alkali-doped fullerides from first principles. Sci. Adv. 1, e1500568 (2015).
Nomura, Y., Sakai, S., Capone, M. & Arita, R. Exotics-wave superconductivity in alkali-doped fullerides. J. Phys. Condens. Matter 28, 153001 (2016).
Rice, M. J., Choi, H. Y. & Wang, Y. R. Three-band superconductivity in K3C60 and Rb3C60. Phys. Rev. B 44, 10414–10416 (1991).
Kondo, J. Superconductivity in transition metals. Prog. Theor. Phys. 29, 1–9 (1963).
Suhl, H., Matthias, B. T. & Walker, L. R. Bardeen-Cooper-Schrieffer theory of superconductivity in the case of overlapping bands. Phys. Rev. Lett. 3, 552–554 (1959).
Acknowledgements
The research leading to these results received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007-2013)/ERC Grant Agreement No. 319286 (QMAC, A.C.). We acknowledge support from the Deutsche Forschungsgemeinschaft (DFG) through the Cluster of Excellence ‘The Hamburg Centre for Ultrafast Imaging’ (EXC 1074 – project ID 194651731, A.C.). We thank M. Volkmann and P. Licht for their technical assistance. We are also grateful to B. Fiedler and B. Höhling for their support in the fabrication of the electronic devices used on the measurement setup and to J. Harms for assistance with graphics.
Funding
Open access funding provided by Max Planck Society.
Author information
Authors and Affiliations
Contributions
E.R., B.L., M. Först, M.B. and A.C. conceived the experiment. A.C. supervised the project. The setup shown in Fig. 2 was built and related measurements were performed by M.B. and G.J. The setup shown in Fig. 3 was built by B.L. with support from M. Först, and related measurements were performed by E.R., B.Y. and Y.Z. Data analysis was performed by E.R. and M.B. K3C60 samples were provided by D.P. and M.R. DFT frozen-phonon calculations were performed by M. Fechner. The single-molecule Hamiltonian treatment was developed by G.J. The manuscript was written by E.R., M.B., G.J., M. Fechner and A.C. with contributions from all other authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Physics thanks the anonymous reviewers for their contribution to the peer review of this work.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Supplementary Information
Supplementary Discussion Sections 1–10, Figs. 1–11 and Tables 1 and 2.
Source data
Source Data Fig. 2
Source data for optical spectra and time-dependent optical properties plotted in Fig. 2.
Source Data Fig. 3
Source data for optical spectra and time-dependent optical properties plotted in Fig. 3.
Source Data Fig. 4
Source data for fluence-dependent optical properties and frequency-dependent photosusceptibility plotted in Fig. 4.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Rowe, E., Yuan, B., Buzzi, M. et al. Resonant enhancement of photo-induced superconductivity in K3C60. Nat. Phys. 19, 1821–1826 (2023). https://doi.org/10.1038/s41567-023-02235-9
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41567-023-02235-9
This article is cited by
-
Dynamical onset of light-induced unconventional superconductivity—a Yukawa-Sachdev-Ye-Kitaev study
Communications Physics (2024)
-
Theory of resonantly enhanced photo-induced superconductivity
Nature Communications (2024)
-
Magnetic field expulsion in optically driven YBa2Cu3O6.48
Nature (2024)
-
Stronger pairs with resonant excitation
Nature Physics (2023)
-
Superconducting nonlinear transport in optically driven high-temperature K3C60
Nature Communications (2023)