We review our recent time- and angle-resolved photoemission spectroscopy experiments, which measure the transient electronic structure of optically driven graphene. For pump photon energies in the near infrared (), we have discovered the formation of a population-inverted state near the Dirac point, which may be of interest for the design of THz lasing devices and optical amplifiers. At lower pump photon energies (), for which interband absorption is not possible in doped samples, we find evidence for free carrier absorption. In addition, when mid-infrared pulses are made resonant with an infrared-active in-plane phonon of bilayer graphene (), a transient enhancement of the electron-phonon coupling constant is observed, providing interesting perspective for experiments that report light-enhanced superconductivity in doped fullerites in which a similar lattice mode was excited. All the studies reviewed here have important implications for applications of graphene in optoelectronic devices and for the dynamical engineering of electronic properties with light.
INTRODUCTION
Dynamical engineering of the electronic structure of solids with light is emerging as an important new area of research, which complements more traditional routes of control using chemistry, pressure, or magnetic fields. Especially interesting are new techniques that can drive materials close to lattice and other resonances at mid-infrared and THz frequencies. Indeed, in these driven solids, new properties emerge away from equilibrium. These result naturally from modulation with a strong, periodic field that creates new effective Hamiltonians with new eigenstates.1
Resonant excitation of the crystal lattice has, for example, been used to induce superconductivity far above the equilibrium critical temperature in cuprates2–5 and K3C60,6 to induce insulator-to-metal phase transitions,7–9 and to melt magnetic or orbital order.10–12 Recently, the formation of photon-dressed Floquet-Bloch states upon mid-infrared pumping has been shown to result in a topological phase transition in Bi2Se3.13,14 Although these techniques have been primarily applied to strongly correlated materials, the underlying mechanism is more general and should be studied in simpler conditions.
Here, we set out to study the dynamical electronic properties of monolayer and bilayer graphene in three different excitation regimes. (1) For pump photon energies bigger than twice the chemical potential, electron-hole pairs are generated. Photo-excited electrons are found to rapidly accumulate above the Dirac point at the bottom of the conduction band where they form a short-lived population-inverted state (Fig. 1).15 By improving the temporal resolution of the experiment to better than 10 fs, we are able to resolve the primary scattering events that are responsible for the formation of the population inversion (Fig. 2).16,17 (2) For lower pump photon energies, free carriers in the vicinity of the Fermi level are accelerated by the applied laser field and then heated via inelastic scattering. For this excitation regime, we observe a thermal electronic distribution at all times within the experimental energy and time resolution (Fig. 3).15 (3) When the pump frequency is tuned to be resonant to the infrared-active in-plane E1u phonon in bilayer graphene at Eph = 200 meV, in addition to carrier heating via free carrier absorption, the electronic band structure is modulated by the coherent oscillation of the carbon atoms. In this case, we observe anomalous carrier dynamics that we attribute to a transient enhancement of the electron-phonon coupling constant (Fig. 4).18,19
METHODS
Samples
For the present study, we used epitaxial monolayer and bilayer graphene samples on hydrogen-terminated silicon carbide crystals. The 6H-SiC(0001) substrates were hydrogen-etched and annealed in an Argon atmosphere, resulting in carbon coverages of one and two monolayers, respectively. These carbon layers were then decoupled from the underlying substrate by hydrogen intercalation, giving quasi-freestanding monolayer and bilayer graphene.20 The equilibrium band structure of these samples can be found in Ref. 18. For all samples, the top of the valence band is unoccupied, indicating hole-doping with the chemical potential at about −200 meV below the Dirac point. The samples were exposed to air and reinserted into an ultra-high vacuum at the Artemis user facility in Harwell, United Kingdom, where the tr-ARPES experiments were performed. The original band structure was recovered by a mild annealing.
Setup
Tr-ARPES measurements have been performed at the Artemis user facility at Harwell, United Kingdom. Femtosecond laser pulses tunable between 800 nm and 15 μm wavelength are typically used to excite the sample. At a variable delay from this excitation, the sample is illuminated with a second extreme ultra-violet (XUV) probe pulse that ejects photoelectrons from the sample. These photoelectrons are measured with a hemispherical analyzer, yielding the photocurrent as a function of kinetic energy, Ekin, and emission angle, θ. From the kinetic energy, the binding energy, EB, of the electrons in the solid can be obtained, and the measured emission angle can be converted into in-plane momentum, , so that the image on the two-dimensional detector directly maps the band structure along a particular cut in momentum space.
The setup is based on a 30 fs/780 nm/1 kHz Titanium:Sapphire (Ti:Sa) laser. Pump pulses are generated with a commercial optical parametric amplifier (1 to 2 μm) with subsequent difference frequency generation (4 to 15 μm). About 1 mJ of laser energy is used for high-order harmonics generation in an Argon gas jet, generating XUV photons in the energy range between 10 and 40 eV. A single-grating time-preserving monochromator is used to select a single harmonic at , yielding <30 fs long XUV pulses with a spectral width of about 300 meV.21 The XUV probe is focused onto the sample with a toroidal mirror resulting in a typical spot size of 200 μm.
In order to improve the temporal resolution of the setup to below 10 fs, the output spectrum of the Ti:Sa laser was broadened in a Neon filled fiber and recompressed using chirped mirrors. The compressed pulses were used both for high order harmonics generation and to pump the sample.16
The photocurrent was measured along the ΓK-direction in the vicinity of the K-point as a function of momentum, energy, and pump-probe delay. Electron distribution functions have been obtained by integrating the data over momentum.15,18,22
RESULTS
Interband excitation
The interband excitation regime is discussed in detail in Fig. 1. In Fig. 1(a), we show a schematic of the excitation mechanism. Measured electron distributions for various pump-probe delays before, during, and after excitation are shown in Fig. 1(b).15 At negative time delays, before the arrival of the pump pulse, the measured electron distributions (blue data points in Fig. 1(b)) follow a Fermi-Dirac distribution (continuous black lines). Upon arrival of the pump pulse, however, clear deviations from a single Fermi-Dirac function develop (red data points in Fig. 1(b)). In this case, the data can be fitted with two Fermi functions, one with the chemical potential inside the valence band and another one with the chemical potential inside the conduction band (dashed black lines in Fig. 1(b)). The two chemical potentials merge within ∼130 fs and a single Fermi-Dirac distribution is recovered (yellow data points in Fig. 1(b)). We attribute the transient distribution of electrons described by two separate chemical potentials located in the valence and conduction bands, respectively, to a population-inverted state.15
(a) Schematic of direct interband excitation. (b) Momentum-integrated photocurrent in the vicinity of the K-point in monolayer graphene as a function of delay (left) together with distribution functions at selected pump-probe delays (right) for .
(a) Schematic of direct interband excitation. (b) Momentum-integrated photocurrent in the vicinity of the K-point in monolayer graphene as a function of delay (left) together with distribution functions at selected pump-probe delays (right) for .
The measurements displayed in Fig. 1 show a thermal electron distribution at all times, either described by a single or two separate Fermi-Dirac distributions for valence and conduction bands. This indicates that the temporal resolution of these measurements (σ = 35 fs, FWHM = 80 fs, determined from the width of the rising edge of the pump probe signal) is too bad to allow access to the primary thermalization events.
Extreme timescales
The ultrafast <30 fs dynamics after photo-excitation are predicted to be dominated by different Auger processes as sketched in Figs. 2(a) and 2(b).23,24 A photo-excited electron-hole pair can recombine and transfer its energy to a second electron inside the conduction band (Fig. 2(a)). This process (Auger recombination) is dominant in conventional semiconductors with a parabolic dispersion.25–28 In graphene, due to the conical dispersion and the absence of a band gap, this process is strongly suppressed by the lack of occupied (empty) states at the bottom (top) of the conduction (valence) band. Instead, carrier relaxation is expected to be dominated by impact ionization (Fig. 2(b)), where the excess energy of an electron high up in the conduction band is used to generate a second electron-hole pair.23,24,29
Primary thermalization events tracked down with sub 10 fs pulses. (a) Schematic of Auger heating. This process is expected to be strongly suppressed in graphene due to the lack of empty (occupied) states at the top (bottom) of the valence (conduction) band. (b) Schematic of impact ionization. The presence of occupied (empty) states at the top (bottom) of the valence (conduction) band results in a large phase space for impact ionization. (c) Comparison of the total number of carriers inside the conduction band NCB and their average kinetic energy ECB/NCB as a function of delay. For about 25 fs after photoexcitation, NCB keeps increasing while ECB/NCB already decreases, indicating impact ionization.
Primary thermalization events tracked down with sub 10 fs pulses. (a) Schematic of Auger heating. This process is expected to be strongly suppressed in graphene due to the lack of empty (occupied) states at the top (bottom) of the valence (conduction) band. (b) Schematic of impact ionization. The presence of occupied (empty) states at the top (bottom) of the valence (conduction) band results in a large phase space for impact ionization. (c) Comparison of the total number of carriers inside the conduction band NCB and their average kinetic energy ECB/NCB as a function of delay. For about 25 fs after photoexcitation, NCB keeps increasing while ECB/NCB already decreases, indicating impact ionization.
These two processes can be easily identified by measuring the number of carriers inside the conduction band, NCB, as well as their average kinetic energy, ECB/NCB. In the case of Auger recombination, NCB decreases while ECB/NCB increases and vice versa for impact ionization. These numbers are directly accessible in a tr-ARPES experiment, as the photocurrent is proportional to the number of electrons at any given energy.16
In order to resolve these primary scattering events, we repeated the experiment with a temporal resolution of σ = 9 fs (FWHM = 21 fs) using a pump photon energy of . Figure 2(c) shows both NCB as well as ECB/NCB as a function of pump-probe delay. We find that, during the first ∼25 fs after photo-excitation, NCB keeps increasing while ECB/NCB already decreases, clearly indicating impact ionization.16
Free carrier absorption
For experiments in which the pump photon energy was below the threshold for interband excitation, metallic carriers in our hole-doped samples absorb energy from the pump by phonon-assisted intraband absorption (Fig. 3(a)). This excitation mechanism dominates the pump-probe signal in our samples if . In Fig. 3(b), we show electron distributions for different delays after excitation at . We find that a thermal electronic distribution is maintained at all times within the experimental resolution. The only effect of the pump pulse is a broadening of the distribution, indicating an elevated electronic temperature.15
(a) Schematic of free carrier absorption. (b) Momentum-integrated photocurrent in the vicinity of the K-point in monolayer graphene as a function of delay (left) together with distribution functions at selected pump-probe delays (right) for .
(a) Schematic of free carrier absorption. (b) Momentum-integrated photocurrent in the vicinity of the K-point in monolayer graphene as a function of delay (left) together with distribution functions at selected pump-probe delays (right) for .
Phonon pumping
The excitation schemes discussed above only redistribute the electrons inside the Dirac cone and leave both the single particle dispersion as well as the coupling strength of various many-body interactions unaffected. Therefore, in order to manipulate the electronic properties with light, one needs to apply different driving schemes. In this section, we choose the pump frequency to be resonant with the infrared-active in-plane bond stretching E1u phonon in bilayer graphene at Eph = 200 meV. (Monolayer graphene does not possess any infrared-active modes and cannot be used for the present purpose.) In this case, in addition to carrier heating via free carrier absorption (Fig. 4(a)), the crystal lattice and thus the band structure is periodically modulated (Fig. 4(b)). Both the electronic peak temperature (Fig. 4(c)) and the fast cooling time (Fig. 4(d)), typically interpreted in terms of optical phonon emission, are found to decrease when the pump pulse is tuned to the phonon resonance. We attribute this to a transient increase of the electron-phonon coupling constant.18,19
Low-energy excitations of bilayer graphene. (a) Schematic of free carrier absorption in bilayer graphene. This process dominates carrier heating if . (b) Resonant excitation of the in-plane bond-stretching phonon at periodically splits the Dirac cone along the ΓK-direction. (c) Peak electronic temperature as a function of pump wavelength. (d) Initial fast cooling time attributed to the emission of strongly coupled optical phonons as a function of pump-probe delay. Green (red) data points are for on- (off-) resonance excitation. Measurements in (c) and (d) were performed for a constant fluence of 0.26 mJ/cm2. Dashed blue lines in (c) and (d) show the optical conductivity associated with the driven phonon from Ref. 52.
Low-energy excitations of bilayer graphene. (a) Schematic of free carrier absorption in bilayer graphene. This process dominates carrier heating if . (b) Resonant excitation of the in-plane bond-stretching phonon at periodically splits the Dirac cone along the ΓK-direction. (c) Peak electronic temperature as a function of pump wavelength. (d) Initial fast cooling time attributed to the emission of strongly coupled optical phonons as a function of pump-probe delay. Green (red) data points are for on- (off-) resonance excitation. Measurements in (c) and (d) were performed for a constant fluence of 0.26 mJ/cm2. Dashed blue lines in (c) and (d) show the optical conductivity associated with the driven phonon from Ref. 52.
DISCUSSION
Most of the pump-probe data on graphene is interpreted in terms of a three-temperature model (see, e.g., Refs. 30 and 31). On time scales short compared to the typical experimental resolution the photo-excited electrons thermalize into a Fermi-Dirac distribution at an elevated electronic temperature, Te.32,33 In a second step, on a timescale on the order of 200 fs, the electrons cool down by the emission of a subset of strongly coupled optical phonons, mainly the E2g mode at q = Γ and the mode at q = K.34–39 Once the electrons and the optical phonons have reached the same temperature, further cooling of the electron-optical-phonon system occurs via emission of acoustic phonons within several picoseconds.30,34–37,40 Similar models have been applied to describe the cooling dynamics of hot carriers within many different sample systems such as metals,41,42 semimetals,43 semiconductors,44 topological insulators,45 and high-temperature superconductors.46
In this discussion, we want to focus on the peculiarities of photo-excited graphene that go beyond the well-known conclusions obtained from a three-temperature model, i.e., (1) the observation of a transient population inversion, (2) pre-thermal scattering events dominated by impact ionization, and (3) anomalous dynamics due to resonant excitation of the crystal lattice.
(1) and (2) are related to the conical band structure of graphene and the vanishing density of states at the Dirac point. While the occurrence of a population-inverted state after strong optical excitation is commonly observed in conventional semiconductors with a sizable band gap, it is not a priori clear whether such a state might exist in graphene. Ryzhii et al.47 predicted that photo-excited electron-hole pairs in graphene relax by optical phonon emission, resulting in an accumulation of electrons (holes) at the bottom (top) of the conduction (valence) band and a negative optical conductivity in the THz regime. The effect has been demonstrated by Li et al.48 and recently confirmed by more detailed experimental15 and theoretical investigations.49,50 Aside from the interaction with optical phonons, several theoretical23,24 and experimental papers29 have stressed the importance of different Auger scattering channels in photo-excited graphene.
In our tr-ARPES experiments performed in different excitation regimes and with different temporal resolutions, we are able to observe the complete dynamics during the buildup and decay of the population inversion.17 From Fig. 2, we know that impact ionization rapidly accumulates the photo-excited electrons at the bottom of the conduction band, assisting the buildup of the population-inversion via optical phonon emission.16 While two separate Fermi-Dirac distributions can be clearly distinguished for about 130 fs in Fig. 1(b),15 complete thermalization of the carriers into a single Fermi-Dirac distribution takes about 250 fs.17 The decay of the population inversion within hundreds of femtoseconds is dominated by Auger recombination with a minor contribution of optical phonon emission.49
The observed anomalous carrier dynamics that occur when the E1u phonon is resonantly excited can be understood in terms of a transiently enhanced electron-phonon coupling constant.18,19 In Figs. 4(c) and 4(d), we observe lower peak electronic temperatures as well as faster cooling rates at resonance with the E1u mode.18 A more detailed analysis of the energy-dependence of the scattering rates shows that resonant excitation of the phonon results in systematically higher scattering rates at all energies compared to off-resonance excitation.19 In Ref. 19, we simulated the coupled electron-phonon dynamics using Boltzmann equations and showed that neither the difference in the measured electron distribution nor an enhanced occupancy of the driven phonon mode can reproduce our data, leaving an enhanced electron-phonon coupling constant as the only plausible explanation for our findings.
This can also be understood from intuitive arguments without the need to invoke theoretical models. An increase in electronic temperature increases the phase space for optical phonon emission and is expected to result in faster cooling of the electronic system.51 We, however, observe a faster cooling rate for smaller electronic temperatures (Figs. 4(c), 4(d) and Ref. 18). Further, an increased occupation of the driven phonon will enhance the probability for both stimulated phonon emission and phonon absorption. Phonon emission (carrier cooling) only dominates over phonon absorption (carrier heating) in the presence of a population-inverted electronic distribution. We observe a thermal electronic distribution at all times before, during, and after excitation in the mid-infrared.15,18 Finally, in the framework of a two-temperature model, an increased cooling rate at smaller temperature difference between electrons and phonons can only come from an increased coupling strength between the two systems.
CONCLUSION
We have reviewed the electronic-structural dynamics of optically excited graphene at different pump wavelengths throughout the near- and mid-infrared. These dynamics were tracked with time- and angle-resolved photoemission spectroscopy. We find that free-carrier absorption efficiently heats the electrons to high temperatures but maintains a thermal electronic distribution at all times. After interband excitation, carriers rapidly accumulate at the bottom of the conduction band, resulting in a short-lived population-inverted state. Using sub 10 fs pulses, we were able to show that the predominant scattering channel during the first 25 fs after interband excitation is impact ionization. Finally, resonant excitation of the E1u mode in bilayer graphene leads to faster relaxation at lower peak electronic temperatures that we attribute to a transient enhancement of the electron-phonon coupling constant.
While the observed short-lived population inversion may in principle be used to amplify TeraHertz pulses, the presence of impact ionization most likely cannot be exploited for efficient light harvesting as the absence of a band gap in graphene makes electron-hole separation difficult.
The observed transient enhancement of the electron-phonon coupling strength presents a significant step towards dynamical band structure engineering with important consequences for the occurrence and the understanding of light-induced superconductivity in graphene and other carbon compounds, such as K3C60.6
ACKNOWLEDGMENTS
We thank F. Calegari, S. Aeschlimann, M. Chávez-Cervantes, M. Mitrano, C. Cacho, E. Springate, J. Petersen, C. R. Ast, and U. Starke for their contributions to a successful outcome of the experiments and B. Sachs, T. O. Wehling, S. Lichtenstein, M. Eckstein, and M. A. Sentef for theoretical support. The research presented here received funding from Laserlab Europe, STFC, and the German Science Foundation (SPP 1459 and SFB 925).