Stimulated emission of Cooper pairs in a high-temperature cuprate superconductor

Stimulated emission of Cooper pairs in a high-temperature cuprate superconductor

Wentao Zhang Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA Department of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China    Tristan Miller Department of Physics, University of California, Berkeley, California 94720, USA Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA    Christopher L. Smallwood Department of Physics, University of California, Berkeley, California 94720, USA Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA    Yoshiyuki Yoshida    Hiroshi Eisaki Electronics and Photonics Research Institute, National Institute of Advanced Industrial Science and Technology, Ibaraki 305-8568, Japan    R. A. Kaindl Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA    Dung-Hai Lee    Alessandra Lanzara Department of Physics, University of California, Berkeley, California 94720, USA Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
July 29, 2019

The concept of stimulated emission of bosons has played an important role in modern science and technology, and constitutes the working principle for lasers. In a stimulated emission process, an incoming photon enhances the probability that an excited atomic state will transition to a lower energy state and generate a second photon of the same energy. It is expected, but not experimentally shown, that stimulated emission contributes significantly to the zero resistance current in a superconductor by enhancing the probability that scattered Cooper pairs will return to the macroscopically occupied condensate instead of entering any other state. Here, we use time- and angle-resolved photoemission spectroscopy to study the initial rise of the non-equilibrium quasiparticle population in a BiSrCaCuO cuprate superconductor induced by an ultrashort laser pulse. Our finding reveals significantly slower buildup of quasiparticles in the superconducting state than in the normal state. The slower buildup only occurs when the pump pulse is too weak to deplete the superconducting condensate, and for cuts inside the Fermi arc region. We propose this is a manifestation of stimulated recombination of broken Cooper pairs, and signals an important momentum space dichotomy in the formation of Cooper pairs inside and outside the Fermi arc region.

Numerous studies have utilized ultrafast infrared pulses to break Cooper pairs and create non-equilibrium quasiparticles (NEQPs) in conventional and high-temperature superconductorsGedik2004 (); Kaindl2005 (); Coslovich2011a (); Smallwood2012 (); Zhang2013 (). So far, most pump-probe spectroscopic studies on superconductors have focused on the recombination dynamics of NEQPs on picosecond time scaleGedik2005 (); Perfetti2007 (); Graf2011 (); Cortes2011 (); Smallwood2012 (). Such phenomena is usually understood within the phenomenological Rothwarf-Taylor (RT) model, whose basic ingredient is simple bimolecular recombinationRothwarf1967 (); Gedik2004 (); Smallwood2012 (). On the contrary, processes like stimulated emission play a greater role in the faster femtosecond buildup dynamic and become negligent at picosecond time scale. Important clues into what governs the short time scale dynamic may be provided by the emerging technique of time and angle resolved photoemission spectroscopy (trARPES), which adds simultaneous measurement of the superconducting gap recovery as well as momentum-dependent information about NEQP dynamicsSmallwood2012a (); Smallwood2014 ().

Here we report a systematic trARPES study of the temperature, momentum, doping, and density dependence of the initial buildup of NEQPs in high temperature superconductor BiSrCaCuO (Bi2212). In the superconducting state, and at low pump fluence, it takes ps for the NEQP population to reach the maximum. In contrast, in the normal state, or when the pump fluence is sufficiently strong to deplete the entire superconducting condensate, the maximum in the NEQP population is reached immediately (within our time resolution 300 fs) after the arrival of the pump pulse. The NEQP population residing outside the Fermi arc also builds up almost instantly. Along the Fermi arc, the NEQP buildup time shortens as the energy gap increases. We propose that these findings provide evidence of stimulated recombination of quasiparticles into the superconducting condensate.


Buildup of non-equilibrium states. Figure 1 shows the buildup of the nodal NEQP population for an underdoped Bi2212 superconductor ( = 78 K) above and below the superconducting transition temperature. Fig.1a shows the equilibrium photoemission intensity map (before the arrival of the pump pulse) in the superconducting state and the pump-induced change of the same map for different delay times ( represents the arrival time of the pulse). Here the pump fluence (2.4 Jcm) is far below the critical fluence ( Jcm) needed to fully deplete the superconducting condensateZhang2013 (). We observe a clear decrease in the photoemission intensity below the Fermi level, and an increase above the Fermi level, signifying the excitation of NEQPs. However, the change in intensity is greater at ps than at , indicating that the NEQP population requires a finite time to build up, and it can be evidenced in the integrated intensity as a function of energy shown in the left panel of Fig.1b.

Fig.1c shows the integrated change in the intensity above the Fermi level, defined as Perfetti2007 (); Cortes2011 (); Graf2011 (); Smallwood2012 (),


as a function of time for T = 20 K T and T = 90 K T. Consistent with Fig. 1a, the maximum change in intensity occurs at a much later time than . We will define the time at which is near maximum, as the “build-up” time. Later, we will describe a procedure to extract . For the decay rate of NEQP population above is much faster than that at below , consistent with other reports in the literatureGedik2004 (); Gedik2005 (); Smallwood2012 (). From now on, we focus our attention on for the usually overlooked . Strikingly, in the superconducting state, the maximum intensity change () is reached at 1 ps, while in the non-superconducting state the maximum intensity change is reached instantly. If is integrated only above the maximum superconducting gap (), no delay in is observed. In supplementary Figure 1 we demonstrate that the finite cannot be due to the limitation of our time resolution.

Pump intensity dependence of the buildup of non-equilibrium states. Figure 2a shows the low temperature () for different pump fluences. Very different buildup behaviors are observed for fluences below and above the critical fluence ( Jcm). Specifically, at pump fluence equal to 24 Jcm, reaches its maximum instantly within our time resolution. This fast buildup is similar to the normal state data in Fig. 1c. On the other hand at lower pump fluence, appears as large as 0.9 ps. In marked contrast with the low temperature data, for the buildup of non-equilibrium electrons does not show a significant delay for any pump fluence studied (Fig. 2b).

To precisely define the buildup time, , we fit using a resolution-convolved function given by


where is the long-time NEQP population decay rate. The Gaussian accounts for the experimental time resolution, and its full width at half maximum (FWHM) is 300 fs. Equation 2 captures the dynamics of the initial buildup of non-equilibrium quasiparticles quite well (see solid lines in Fig. 2a). Fig. 2c shows the extracted as a function of pump fluence. At low temperature, is 400 fs at the lowest fluence, and decreases as the fluence increases until it saturates at the normal state value of 100 fs. The saturation occurs at the critical fluence required to deplete the superconducting condensateZhang2013 (), but note that the superconducting gap does not close until about 0.6 ps laterSmallwood2014 (). The absence of a buildup time above or above the critical fluence (Fig. 2a) suggests a connection between the slow buildup and the presence of the superconducting condensate.

Momentum dependence of the buildup of non-equilibrium states. In momentum space, similar delayed buildup is observed in throughout the Fermi arc region. Fig. 3a shows the momentum-dependent buildup time in an optimally doped sample in response to a weak (2.4 Jcm) pump. As previously reportedSmallwood2012 (), the long-time decay rate () of the NEQP population increases from node to off-node. Inside the Fermi arc region, there is a delay in buildup, which decreases away from the node. In contrast, for a momentum cut outside the Fermi arc region, there is no delay, and drops to the normal state value of about 100 fs. We have also compared the buildup time for the NEQP population in samples with different doping level (overdoped sample, Fig. 3b), and the same momentum dependence of is observed: the build up time is non zero only on the Fermi arc. For the overdoped sample, all the momentum cuts go through the Fermi arc region and hence they all show an apparent delay in the buildup of the NEQP population. However, compared with the optimally- and under-doped samples, the delay is less pronounced (Fig. 3c). The momentum dependence of the build up rate for the NEQP population is shown in Fig. 3d. For each momentum we extract the corresponding gap value. The NEQP is then plotted as a function of gap size. Gap values are extracted using a standard procedure of fitting the symmetrized energy distribution curves at the Fermi momentum to phenomenological single particle spectral functionsNorman1998 () (supplementary Figure 2 and supplementary Discussion 2). It is clear that there is a distinct mechanism for the buildup time for the states beyond the Fermi arc region.

Note that the difference cannot be attributed to the lack of a gap inside the Fermi arc region, since the Fermi arc itself does not develop until about 0.6 ps after pumpingSmallwood2014 (). Interestingly, the build up time of the NEQP population is identical to that of the photoinduced shift of the superconducting gap at the same momentum (Fig. 3e). The similar timescale further suggests that the superconducting condensate is a key ingredient in establishing the build up time.


In the following we propose an explanation of these observations that takes into account the quantum coherence of the Cooper pair condensate. When the pump pulse photoexcites the material, the input energy is partially stored in the electronic degrees of freedom and partially in the lattice degrees of freedom. Where there is a superconducting condensate, those NEQPs whose time reversal partner is also excited can quickly recombine into Cooper pairs due to a process of coherent recombination stimulated by the existing pair condensate. This process serves to reduce the proportion of input energy initially stored in electronic degrees of freedom, hence causing the observed delay in the buildup of the NEQP population. Such a stimulated recombination process is absent when the superconducting condensate is depleted by either the temperature or the high pump fluence. After , what are left are “unpaired” excited quasiparticles, meaning the time reversed partners of these quasiparticles are not simultaneously excited. At a longer time scale, due to momentum-changing impurity, some subset of those “unpaired” excited quasiparticles will change into “paired” ones and quickly recombine. The bottleneck of the recombination process is the conversion of unpaired NEQPs into paired ones, which is considerably slowerRothwarf1967 (); Gedik2004 (); Smallwood2012 (). This process is depicted as real space diffusion triggered “bimolecular” recombination.

A cartoon of of this physical scenario is shown in Fig. 4a. Note that this explanation is similar to that of earlier optical spectroscopy studies, in that the finite buildup time arises because much of the initial photoexcitation energy is stored in the lattice degrees of freedomDemsar2003 (); Kusar2008 (); Beck2011 (). However, the main difference is that we attribute this initial energy specifically to the presence of the superconducting condensate. This is a key aspect to account for the temperature, energy and fluence dependence here reported, which otherwise cannot be explained.

The importance of including stimulated recombination into the RT equations is summarized in Fig. 4b, where the delay dependence of the number of non-equilibrium electrons is shown. Without superconducting condensate (black curve) no buildup time is observed. However, when the RT equations are augmented to include the stimulated recombination channel (see Supplementary Figure 3 and Discussion 3), the results are qualitatively consistent with the experimental observations (see red curve in Fig. 4b) . Note that due to the electronic inhomogeneity of Bi-2212 the momentum conservation rules are relaxed, i.e. when “paired” NEQP recombine they can emit phonons (bosons) with a sufficiently wide range momentum. Within this scenario, the momentum dependence of the buildup time can be easily accounted for by the momentum dependence of the recombination rate (), that increases as we move away from the nodeSmallwood2012 (). Additionally, the absence of a buildup time beyond the Fermi arc tip, put a strong constraint on the presence of a condensate in this momentum region, another important manifestation of the nodal off nodal dichotomyFu2006 ().

One alternate interpretation of our results is that initial NEQPs far from the node scatter towards the node, allowing the NEQP population near the node to build up over a longer period of time. Such scattering would be induced by the anisotropic form of the superconducting gap. However, this interpretation does not account for the absence of a buildup time at high fluences or high temperatures. Above the critical fluence, the superconducting gap is not suppressed until about 0.6 ps after the pump pulseSmallwood2014 (), some time after the NEQP population has reached its maximum, and above , there is still an anisotropic pseudogap.

To conclude, our main finding is that upon photoexcitation the delay in the buildup of non-equilibrium quasiparticle population in higher temperature superconductors is sensitive to whether there is a superconducting condensate. We propose that our observation reveals the stimulated recombination of photoexcited quasiparticles.


Measurements were taken on our home-built trARPES system. The experimental setup is the same as that reported previouslyGraf2011 (); Smallwood2012 (); Smallwood2012a (); Zhang2013 (). Transient states in the solids are created by pumping with infrared laser pulses ( eV), then they are probed by a following ultra-violet pulse ( eV) using ARPES setup equipped with a SPECS PHOIBOS 150 electron spectrometer. Time resolution is achieved by varying the delay time between the probe and pump pulses, with a total resolution of 300 fs determined by a cross-correlation of pump and probe pulses measured on hot electrons. The total energy resolution is 23 meV characterized by fitting the Fermi edge to a resolution-convolved Fermi-Dirac distribution function. Samples with superconducting transition of underdoped ( K, UD78K) and optimally doped ( K, OP91K) and overdoped ( K, OD78K) were studied. The samples were cleaved in situ at 20 K in ultra high vacuum with a base pressure lower than mbar.


  • (1) Gedik, N. et al. Single-quasiparticle stability and quasiparticle-pair decay in . Phys. Rev. B 70, 014504– (2004).
  • (2) Kaindl, R. A., Carnahan, M. A., Chemla, D. S., Oh, S. & Eckstein, J. N. Dynamics of cooper pair formation in . Phys. Rev. B 72, 060510– (2005).
  • (3) Coslovich, G. et al. Evidence for a photoinduced nonthermal superconducting-to-normal-state phase transition in overdoped . Phys. Rev. B 83, 064519– (2011).
  • (4) Smallwood, C. L. et al. Tracking cooper pairs in a cuprate superconductor by ultrafast angle-resolved photoemission. Science 336, 1137–1139 (2012).
  • (5) Zhang, W. T. et al. Signatures of superconductivity and pseudogap formation in nonequilibrium nodal quasiparticles revealed by ultrafast angle-resolved photoemission. Phys. Rev. B 88, 245132– (2013).
  • (6) Gedik, N. et al. Abrupt transition in quasiparticle dynamics at optimal doping in a cuprate superconductor system. Phys. Rev. Lett. 95, 117005– (2005).
  • (7) Perfetti, L. et al. Ultrafast electron relaxation in superconducting by time-resolved photoelectron spectroscopy. Phys. Rev. Lett. 99, 197001– (2007).
  • (8) Graf, J. et al. Nodal quasiparticle meltdown in ultrahigh-resolution pump-probe angle-resolved photoemission. Nat Phys 7, 805–809 (2011).
  • (9) Cortés, R. et al. Momentum-resolved ultrafast electron dynamics in superconducting . Phys. Rev. Lett. 107, 097002 (2011).
  • (10) Rothwarf, A. & Taylor, B. N. Measurement of recombination lifetimes in superconductors. Phys. Rev. Lett. 19, 27–30 (1967).
  • (11) Smallwood, C. L., Jozwiak, C., Zhang, W. T. & Lanzara, A. An ultrafast angle-resolved photoemission apparatus for measuring complex materials. Review of Scientific Instruments 83, 123904 (2012).
  • (12) Smallwood, C. L. et al. Time- and momentum-resolved gap dynamics in . Phys. Rev. B 89, 115126– (2014).
  • (13) Norman, M. R., Randeria, M., Ding, H. & Campuzano, J. C. Phenomenology of the low-energy spectral function in high- superconductors. Phys. Rev. B 57, R11093–R11096 (1998).
  • (14) Demsar, J. et al. Pair-breaking and superconducting state recovery dynamics in . Phys. Rev. Lett. 91, 267002– (2003).
  • (15) Kusar, P. et al. Controlled vaporization of the superconducting condensate in cuprate superconductors by femtosecond photoexcitation. Phys. Rev. Lett. 101, 227001– (2008).
  • (16) Beck, M. et al. Energy-gap dynamics of superconducting thin films studied by time-resolved terahertz spectroscopy. Phys. Rev. Lett. 107, 177007 (2011).
  • (17) Fu, H. & Lee, D.-H. Dichotomy between the nodal and antinodal excitations in high-temperature superconductors. PHYSICAL REVIEW B 74 (2006).
This work was supported by Berkeley Lab’s program on Ultrafast Materials Sciences, funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Materials Sciences and Engineering Division, under Contract No. DE-AC02-05CH11231.

Author contributions

W. Z. conceived the project and the experiments together with A. L. The experiments were taken by W. Z., T. M. and C. L. S. Samples were grown and characterized by Y. Y. and H. E. The manuscript was written by W. Z., A. L. and D. H. L. with input from all authors.

Additional information

The authors declare no competing financial interests.

Figure 1: Buildup time in the nodal non-equilibrium electron excitations. Measurements were taken along a nodal cut ((0, 0) – Y(, )) in an underdoped Bi2212 ( 78 K). a, Equilibrium photoemission intensity at delay time ps as a function of energy and momentum and pump-induced change of photoemission spectra as a function of energy and momentum at delay time 0 and 0.9 ps. The intensities as a function of energy by integrating the original photoemission spectra at delay time -1.2 ps, 0 ps and 0.9 ps are shown in the left panel. The measurement was taken at 20 K with a pump fluence of 2.4 Jcm. b, The integrated photoemission intensity across the arrowed momentum region for the three different delay time shown in (a). c, The spectral gain , by integrating the change of photoemission intensity above the Fermi level at 20 K, above the energy of maximum superconducting gap at 20 K, and above at 90 K as a function of delay time at pump fluence 2.4 Jcm. The width of the pump pulse in the cartoon is determined by frequency-resolved optical gating (FROG).
Figure 2: Pump fluence dependence of buildup time along nodal direction. a and b, as a function of delay time and pump fluences at 20 K and 90 K, respectively, for an underdoped sample ( K). The shaded line is a guide to the eyes. c, The extracted buildup time from a as a function of pump fluence.
Figure 3: Momentum dependence of the buildup time. The measurements were taken at 20 K on nearly optimally (91 K) and overdoped (78 K) Bi2212 samples with pump fluence of 2.4 Jcm. a and b, as a function of delay time measured inside (red curves) and outside (blue curve) Fermi arc for optimally doped sample (panel a) and inside the Fermi arc for an overdoped samples (panel b). c, Doping dependence of the buildup time at the node. d, The extracted buildup rate from a as a function of energy gap at the Fermi momentum for four cuts shown in inset of b. The green data point is a separate measurement at the Fermi surface angle of . The brown region e, Comparison between and the change of energy gap as a function of delay time at the Fermi surface angle of .
Figure 4: Response of a superconductor to pump photons. a, Cartoons to show the response of a superconductor to a pump laser for fluences lower and higher than the critical fluence. b, Simulations of the numbers of non-equilibrium electrons () as a function of delay time. A non-zero buildup time is observed only in the presence of a condensate. The curves are the numerical solutions of the equations in supplementary Discussion 3 using the same parameters as shown in supplementary Figure 3b but with for the absence of the stimulated emission of Cooper pairs into the condensate and for non-vanishing stimulated emission.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description