Resonant enhancement of photo-induced superconductivity in K3C60

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. There is evidence that K3C60 can host a photo-induced superconducting state. Now, resonant excitation at low frequencies allows this phenomenon at room temperature and low pumping fluence.

Typical experimental results reported until now are displayed in Fig. 2c.K 3 C 60 powders were held at a base temperature T = 100 K ≫ T c = 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 ≪ T c = 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 σ 1 (ω) and an imaginary conductivity σ 2 (ω) that diverges towards low frequencies as σ 2 (ω) ∝ 1/ω.The divergent σ 2 (ω) implies (through Kramers-Kronig relations) the presence of a peak in σ 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 Article https://doi.org/10.1038/s41567-023-02235-9has 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 σ 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 nearinfrared signal beams of two phase-locked optical parametric amplifiers 23 .This source, illustrated schematically in Fig. 3a and described in detail in Supplementary Section 4, was used to generate narrowbandwidth 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 Δτ pump ≈ 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(ω) and σ 1 (ω) (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 superconductingstate 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 σ 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 buckyball 25 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 calculations 26,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 to the terahertz response, the superconducting-like nature of the transient state is independent of the specific choice of assumption 21 .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(ω, τ)), σ 1 (ω, τ) and σ 2 (ω, τ) and plotted as a function of pump-probe time delay τ in Fig. 2e, showing the time evolution of the system.The first two quantities are the frequency-averaged values of R(ω) and σ 1 (ω) 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 σ 2 (ω).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 superconductinglike 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 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 states 26,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 mechanism [28][29][30] .However, it is not yet clear exactly how this excitation is transformed in the presence of tunnelling between neighbouring C 60 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 T 1u 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 C 60 3− 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 t 1u 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 C 60 3− structure, these molecular levels are degenerate.Applying a distortion along a T 1u coordinate lifts this degeneracy, leaving a doubly degenerate t 1u 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

Article
https://doi.org/10.1038/s41567-023-02235-9 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 T 1u 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 K 3 C 60 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 T c .
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.b, Frequency dependence of the photo-susceptibility of K 3 C 60 defined as the gradient of the lost spectral weight in σ 1 in the low-fluence limit (Supplementary Section 8) measured 10 and 50 ps after photo-excitation.These measurements were carried out at a base temperature T = 100 K.These data were obtained by accounting for the inhomogeneous excitation of the probed volume under the assumption of a square-root fluence dependence of the photo-induced changes in the complex refractive index of the material.The solid curves are guides to the eye.

Fig. 1 |
Fig. 1 | Crystal structure and phase diagram of K 3 C 60 .a, Crystal structure of the organic molecular solid K 3 C 60 .C 60 molecules are situated at the vertices of a face-centred-cubic lattice.Potassium atoms (red) occupy the interstitial voids.b, Pressure-temperature phase diagram of the fcc-A 3 C 60 alkali-doped fulleride family of compounds showing the presence of superconducting, metallic, paramagnetic insulating (PI) and anti-ferromagnetic insulating (AFI) phases.Physical pressure tunes the spacing between the C 60 molecules.The grey line indicates the boundary between the insulating and metallic/superconducting compounds.The blue shaded area indicates where superconductivity is observed at equilibrium.The star indicates the K 3 C 60 compound investigated in this work, which superconducts below the critical temperature (T c ) of 20 K.

Fig. 2 |
Fig. 2 | Photo-induced metastable superconductivity in K 3 C 60 generated with intense 170 meV excitation pulses.a, Schematic of the experimental setup.Pump pulses with 170 meV (41 THz) photon energy (close to resonance with the illustrated phonon mode) were generated by means of white light generation (WLG), followed by optical parametric amplification (OPA) and subsequent difference frequency generation (DFG) of the signal and idler beams.These pulses were stretched to a duration of Δτ pump ≈ 1 ps by propagation in a highly dispersive CaF 2 rod.The photo-induced changes in the far-infrared optical properties of K 3 C 60 were detected with phase-sensitive transient terahertz (THz) time-domain spectroscopy.b, Reflectivity (sample-diamond interface) and real and imaginary parts of the optical conductivity of K 3 C 60 measured after cooling across the equilibrium superconducting transition.c, Same quantities measured at equilibrium (red lines) and 10 ps after excitation (symbols).The data in filled blue circles, open blue circles and open brown triangles were obtained accounting for the inhomogeneous excitation of the probed volume under the assumption of a linear, square-root and saturating (respectively) fluence dependence of the photo-induced changes in the complex refractive index of the material (Supplementary Section 6).These data were acquired at a base temperature T = 100 K with an excitation fluence of approx.18 mJ cm −2 .d, Same quantities as in c but measured at a base temperature T = 295 K with an excitation fluence of approx.15 mJ cm −2 .e, Time dependence of the average value (indicated by *) of reflectivity and real part of the optical conductivity σ 1 (ω) in the 5-10 meV spectral range, and light-induced 'superfluid density' (n sf ) extracted from a two-fluid model fit and expressed as a fraction of the total charge carrier density (n tot ).Filled and open symbols indicate the results of the linear and square-root reconstruction models, respectively.The red dotted lines indicate the value of the corresponding quantity in equilibrium.These data were acquired at a base temperature T = 100 K with a fluence of approximately 18 mJ cm −2 and a pump-pulse duration of Δτ pump ≈ 1 ps.Transient optical spectra corresponding to these measurements are reported in Supplementary Fig.6.3.

Fig. 4 |
Photosusceptibility(a.u.) Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.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/.