## Abstract

For many of the envisioned optoelectronic applications of graphene, it is crucial to understand the subpicosecond carrier dynamics immediately following photoexcitation and the effect of photoexcitation on the electrical conductivity—the photoconductivity. Whereas these topics have been studied using various ultrafast experiments and theoretical approaches, controversial and incomplete explanations concerning the sign of the photoconductivity, the occurrence and significance of the creation of additional electron-hole pairs, and, in particular, how the relevant processes depend on Fermi energy have been put forward. We present a unified and intuitive physical picture of the ultrafast carrier dynamics and the photoconductivity, combining optical pump–terahertz probe measurements on a gate-tunable graphene device, with numerical calculations using the Boltzmann equation. We distinguish two types of ultrafast photo-induced carrier heating processes: At low (equilibrium) Fermi energy (*E*_{F} ≲ 0.1 eV for our experiments), broadening of the carrier distribution involves interband transitions (interband heating). At higher Fermi energy (*E*_{F} ≳ 0.15 eV), broadening of the carrier distribution involves intraband transitions (intraband heating). Under certain conditions, additional electron-hole pairs can be created [carrier multiplication (CM)] for low *E*_{F}, and hot carriers (hot-CM) for higher *E*_{F}. The resultant photoconductivity is positive (negative) for low (high) *E*_{F}, which in our physical picture, is explained using solely electronic effects: It follows from the effect of the heated carrier distributions on the screening of impurities, consistent with the DC conductivity being mostly due to impurity scattering. The importance of these insights is highlighted by a discussion of the implications for graphene photodetector applications.

## INTRODUCTION

The potential for technological breakthroughs using graphene is particularly high in the field of optoelectronics, given graphene’s promising properties for terahertz (THz) technologies (*1*, *2*), photodetection (*3*), plasmonics (*4*), light harvesting (*5*), data communication (*6*, *7*), ultrafast laser technologies (*8*, *9*), and more. For all of these applications, it is of primary importance to accurately understand how photoexcitation and the subsequent ultrafast carrier dynamics occur and affect graphene’s conductivity. Here, the ultrafast carrier dynamics correspond to the initial energy relaxation of photoexcited electron-hole pairs, which occurs on a subpicosecond time scale, whereas the graphene DC photoconductivity is related to how momentum scattering takes place in the quasi-equilibrium state after initial energy relaxation.

The ultrafast energy relaxation dynamics of graphene have been studied experimentally using optical pump–optical probe spectroscopy (*10*–*15*), time-resolved angle-resolved photoemission spectroscopy (ARPES) measurements (*16*, *17*), time-resolved photocurrent scanning microscopy (*18*, *19*), ultrafast time-resolved multiphoton photoemission spectroscopy (*20*), optical pump–THz probe spectroscopy (*21*–*33*), and high-field THz spectroscopy (*34*, *35*). These time-resolved studies have identified that the energy relaxation dynamics of the photoexcited carriers consist mainly of carrier-carrier scattering and coupling to optical, acoustic, and remote substrate phonons. Carrier-carrier scattering leads to a very rapid (*t* ≲ 50 fs) broadening of the electron distribution as a function of energy, associated with thermalization within the electronic system, that is, carrier heating (*10*–*12*, *16*, *17*, *21*, *25*, *36*). Some theoretical studies (*37*–*41*) and combined experimental-theoretical work (*12*, *13*, *17*, *36*) have discussed the interband scattering processes that lead to carrier multiplication (CM), where the carrier density in the conduction band exceeds the carrier density immediately after photoexcitation. Other studies have discussed the occurrence of intraband scattering processes that lead to hot-CM (*16*, *25*, *30*, *42*–*44*), where the hot-carrier density (carriers with energy above the chemical potential) exceeds the carrier density immediately after photoexcitation. Important questions concerning the ultrafast carrier dynamics, which have not been addressed, pertain to how carrier heating processes change with Fermi energy, what the different conditions are for the occurrence of CM and hot-CM, and how the ultrafast carrier dynamics correlate with the THz photoconductivity.

The understanding of the conductivity of unexcited graphene comes mainly from transport measurements and theoretical modeling of momentum scattering processes. Typically, in graphene devices with a mobility below 10,000 cm^{2}/Vs, that is, most substrate-supported graphene devices, the conductivity is limited by long-range scattering, including, for example, scattering between carriers and Coulomb impurities in the substrate (*45*–*47*) or at corrugations (*48*). Only in the case of ultraclean, suspended, or hexagonal boron nitride–encapsulated graphene are much higher mobilities achieved because impurity scattering is suppressed. In that case, phonon scattering has to be taken into account as well to explain the observed conductivity (*49*). The conductivity of graphene after photoexcitation has been mainly addressed by optical pump–THz probe measurements (*21*–*33*). This technique consists of photoexcitation using an optical pump pulse that excites one electron per absorbed photon from valence to conduction band, followed by measurement of the transmission of a THz probe pulse. Given the THz photon energy of a few meV, smaller than the Fermi energy of most graphene samples, THz absorption takes place because of intraband electronic transitions and therefore provides a direct, contact-free measurement of graphene’s low-frequency (nearly DC) conductivity. The change in conductivity at THz frequencies following photoexcitation is referred to as the THz photoconductivity. Recent experiments on the (Fermi energy–dependent) THz photoconductivity have indicated that the sign of the photoconductivity can be both positive and negative. Explanations of this surprising effect include the modified carrier distribution due to carrier heating (*21*, *22*, *25*, *30*), increased momentum scattering with optical phonons (*27*, *32*), THz emission (*24*), and different combinations of these effects (*23*, *26*, *28*, *29*). Although all these studies have been performed on substrate-supported graphene devices, the effect of photoexcitation on the dominant momentum scattering mechanism—long-range Coulomb impurity scattering—has so far not been studied nor has the importance of screening been properly taken into account.

Here, we address the key issues in the observed ultrafast carrier dynamics and photoconductivity of graphene, for which so far controversial and partially contradicting explanations have been put forward. We investigate these issues by measuring the THz photoconductivity in a single, well-characterized, substrate-supported graphene device with tunable carrier density. Moreover, we vary the photon energy of the pump pulse, arguably the most revealing parameter for understanding photodynamical processes, such as carrier heating (*19*, *25*) and (hot-)CM processes (*50*). We complement the experimental results with theoretical calculations that take into account all relevant energy relaxation processes and consider that the conductivity is governed by momentum scattering through long-range interactions with impurities, thus only considering electronic effects. This accurately reproduces the observed experimental trends without freely adjustable parameters. Our theory is based on the massless Dirac fermions model for graphene carriers (*51*), which is a low-energy approximation valid around the Dirac point, that is, the crossing of two bands arising from the π orbitals of the carbon atoms. The electronic dispersion around the Dirac point is conical and particle-hole symmetrical, with one conduction and one valence band above and below the Dirac point, respectively (*51*). For convenience, the Fermi energy is measured from the Dirac point in this approximation.

Our combined experimental and theoretical results lead to an intuitive physical picture based on two distinct regimes of carrier heating. (i) In the case of low carrier density (*E*_{F} ≲ 0.1 eV), energy relaxation of the initially excited electron-hole pairs occurs through electronic transitions from valence to conduction band, leading to a broadened carrier distribution, that is, interband heating (here and in the following, the terms interband and intraband refer to the two bands of the massless Dirac fermions model mentioned above). Under certain conditions, this interband heating corresponds to the creation of secondary carriers in the conduction band and thus CM, although this process competes with relaxation of the photoexcited carriers from conduction to valence band due to phonon emission and other carrier-carrier scattering channels. Here, photoexcitation and the subsequent ultrafast dynamics lead to an increase in THz conductivity (positive THz photoconductivity), which can be partially understood from the additional free carriers that are created. (ii) For higher initial carrier densities (*E*_{F} ≳ 0.15 eV), interband heating is suppressed because of the different scattering phase space for interband and intraband processes in graphene (*52*). The most effective ultrafast energy relaxation channel consists of electronic transitions from below the chemical potential to above the chemical potential, that is, intraband heating. Under the conditions studied here, this intraband heating corresponds to the creation of hot carriers in the conduction band and thus hot-CM. Here, photoexcitation and the subsequent ultrafast dynamics lead to a reduction of the conductivity (negative THz photoconductivity) mainly due to reduced screening of the impurities, leading to increased momentum scattering. In short, we find that the transition from interband to intraband heating and the transition from positive to negative THz photoconductivity are determined by the ultrafast carrier dynamics following photoexcitation and are thus naturally correlated. Finally, we show that energy is efficiently transduced from incident light to graphene carriers, regardless of the Fermi energy, meaning that both interband and intraband heating are efficient. Thus, our results show that a large collection of experimental observations related to photoexcited graphene can be explained by a simple and intuitive picture, where energy relaxation occurs through ultrafast and efficient carrier heating, and the photoconductivity is the result of long-range Coulomb scattering, where screening plays a crucial role both for the ultrafast dynamics and for the sign of the photoconductivity.

## RESULTS

### Photoconductivity versus carrier density

We use optical pump–THz probe spectroscopy, as illustrated in Fig. 1A. A visible pump pulse with a duration of ~100 fs, a spot size of ~1 cm, and a photon energy *E*_{ph} typically between 1.0 and 2.5 eV excites electrons from valence band states with energy − *E*_{ph}/2 to conduction band states with energy *E*_{ph}/2. After a time delay *t* ~ 300 fs (see fig. S1A), we probe the pump-induced change in the transmission of a quasi-single cycle THz pulse (center frequency, ~0.6 THz). The time delay is much longer than the thermalization time of the electronic distribution, so the electronic system is in a quasi-equilibrium state (*53*). The probe yields the time-dependent THz photoconductivity Δσ_{THz} = σ_{THz}(*t*) − σ_{THz}(0), where σ_{THz}(0) is the THz conductivity at equilibrium, that is, before the pump pulse. The device in Fig. 1A (see Methods for details) consists of chemical vapor depostion (CVD) graphene on quartz with a drop cast ionic gate deposited on top to control the equilibrium carrier density or, equivalently, the Fermi energy. Through resistance measurements using the source and drain contacts, we extract the equilibrium DC conductivity σ_{0} (see also fig. S1B). We infer a mobility of ~1000 cm^{2}/Vs by combining the transport measurements with Raman measurements to determine the carrier density (see fig. S1C). This mobility confirms that the conductivity is indeed governed by impurity scattering (*45*–*47*). It also indicates a momentum scattering time on the order of tens of femtoseconds, meaning that the THz (photo)conductivity is mainly real and has very little frequency dependence in our experimental THz window (~0.4 to 1.8 THz). Figure 1B shows the optically measured THz photoconductivity Δσ_{THz} as a function of the equilibrium DC conductivity σ_{0}, which is obtained simultaneously using electrical measurements (see Methods for details). We use a fixed fluence of 8 μJ/cm^{2} and a pump wavelength of 800 nm and observe a transition from a positive signal for low steady-state conductivity (low carrier density) to a negative signal for higher steady-state conductivity (higher carrier density). A similar sign change of the THz photoconductivity was observed in the studies of Shi *et al.* (*28*), Frenzel *et al.* (*29*), Jensen *et al.* (*30*), and Wang *et al.* (*33*).

To understand the physics behind these experimental results, we show in Fig. 1C the results of our theoretical calculation of the THz photoconductivity. To model the system, we first solve a semiclassical Boltzmann equation that takes into account carrier-carrier and carrier–optical phonon scattering in graphene. Our method has been described in detail in the studies of Brida *et al.* (*12*) and Tomadin *et al.* (*53*). The calculation is nonperturbative and yields the full nonequilibrium electron and optical phonon distributions in time following photoexcitation. We do not integrate the Boltzmann equation during the photoexcitation process but use a Fermi-Dirac distribution with an additional Gaussian population of photoexcited electrons in the conduction band (centered around *E*_{ph}/2) and a corresponding population of holes in the valence band (centered around − *E*_{ph}/2) as the initial condition for the electron distribution immediately after photoexcitation. For the sake of comparison with experiment, the Fermi-Dirac distribution always has a finite carrier density in the conduction band, that is, positive Fermi energy. We verify that, after a short initial transient *t* ≲ 20 fs, the quasi-equilibrium electron distribution consists of two Fermi-Dirac distributions, in valence and conduction bands, with the same temperature but different chemical potentials. This is in agreement with time- and angle-resolved photoemission spectroscopy, where the carrier distribution is probed as a function of time (*16*, *17*, *36*). The values of the temperature and the chemical potentials, obtained by fitting the electron distribution, are then used to calculate the photoconductivity using the semiclassical Boltzmann equation linearized around the quasi-equilibrium distribution, taking into account scattering between carriers and long-range Coulomb impurities (see Methods for details). The measured device mobility of ~1000 cm^{2}/Vs dictates a surface impurity density of 5 × 10^{12}/cm^{2}.

In the following paragraphs, we will present experimental and theoretical results that elucidate the processes that dominate the ultrafast dynamics following photoexcitation of graphene. Subsequently, we will discuss in detail the reason for the transition from positive to negative photoconductivity. First, we point out that the results in Fig. 1C are not obtained by fitting our theoretical calculation to the experimental results in Fig. 1B, so the agreement between the two is remarkable. Our calculation of the conductivity uses the general Boltzmann equation, rather than simplified versions, such as the Drude equation with separate contributions from momentum scattering and from changes in the electron distribution, captured by the Drude weight (*26*–*28*), or the Sommerfeld expansion of the Boltzmann equation, which is only valid for high doping and low fluence, such that the electron temperature is lower than the Fermi temperature *T*_{F} = *E*_{F}/*k*_{B}, with *k*_{B} as the Boltzmann constant (*25*). Furthermore, we only consider the dominant momentum scattering process, namely, long-range impurities scattering, rather than assuming a modified momentum scattering rate due to enhanced phonon interactions, as in the studies of Shi *et al.* (*28*), Frenzel *et al.* (*29*), Hafez *et al.* (*31*), and Mihnev *et al.* (*32*). Thus, the experimental results are explained by our model that takes into account solely electronic effects, although including secondary effects would likely improve the quantitative agreement between the experimental and the theoretical results. Further improvement could come from including effects of electron-hole puddles (*54*), which could be substantial in substrate-supported CVD graphene devices. Finally, we point out that the photoconductivity is a nonequilibrium property and, as such, that its magnitude and sign change with time in a fashion that depends on the initial conditions (for example, the equilibrium carrier density) and on the details of the perturbation, such as the fluence [see, for example, the study of Jensen *et al.* (*30*)] and the frequency of the pump pulse.

### Interband-to-intraband heating transition

The calculation of the evolution of the electron distribution with time allows us to get insight into the ultrafast energy relaxation that takes place after photoexcitation. To this end, we calculate the carrier density in the conduction band before (*n*_{CB,0}) and after (*n*_{CB}) photoexcitation, according to , where ε is the carrier energy and the formulas for the density of states *ν*(ε) and the quasi-equilibrium carrier distribution *f*_{0}(ε; *t*) are given in Methods. We also calculate the hot-carrier density, that is, the density of electrons in the conduction band with an energy exceeding the chemical potential μ_{e}, before (*n*_{HC,0}) and after (*n*_{HC}) photoexcitation, according to , where μ_{e}(*t*) is defined in Methods. These quantities are shown in Fig. 1D as functions of the equilibrium DC conductivity. We find that at low equilibrium DC conductivity, that is, close to the charge neutrality point, photoexcitation leads to an increase in the carrier density in the conduction band: *n*_{CB} − *n*_{CB,0} > 0. This means that interband scattering processes, where electrons/holes are transferred between conduction/valence band, have occurred. At high conductivity, that is, away from the charge neutrality point, photoexcitation does not lead to a change in the carrier density in the conduction band (*n*_{CB} ≈ *n*_{CB,0}), whereas it leads to an increase in the hot-carrier density: *n*_{HC} − *n*_{HC,0} > 0. This means that intraband scattering events have occurred. Obviously, upon approaching the Dirac point, the hot-carrier density converges to the conduction band carrier density, as we see in Fig. 1D.

The physical interpretation of this interband-to-intraband scattering transition becomes clear when we examine the electron distribution before and after photoexcitation. We calculate these for an equilibrium Fermi energy of 0.05 eV (Fig. 2A) and 0.4 eV (Fig. 2B), corresponding to σ_{0} ≲ 5 *e*^{2}/h and σ_{0} ≃ 50 *e*^{2}/h, respectively (e is the electron charge and h is Planck’s constant). Before photoexcitation, at room temperature *T* = 300 K, the distributions are slightly broadened, and for both equilibrium Fermi energies, the hole density in the valence band is small compared to the electron density in the conduction band. At time *t* = 300 fs following photoexcitation with *E*_{ph} = 1.5 eV, the electron distributions for the two equilibrium Fermi energies have broadened substantially, indicating carrier heating. However, the broadened distributions are qualitatively different, owing to different relaxation processes dominating the dynamics depending on the equilibrium Fermi energy. At *E*_{F} ≃ 0.05 eV, the heated electron distribution contains a significant hole density in the valence band and an increased electron density in the conduction band (see cyan-shaded area in Fig. 2A). This indicates that electronic transitions have taken place from valence band to conduction band. By contrast, at *E*_{F} ≃ 0.4 eV, the heated electron distribution does not contain holes in the valence band. However, it contains an increased hot electron density (see orange-shaded area in Fig. 2B). Thus, we see that the occurrence of both interband and intraband transitions can be understood from heating of the graphene carriers. Furthermore, we see that the transition from positive to negative photoconductivity correlates with a transition from interband heating to dominant intraband heating.

Because our calculations show an increased (hot) carrier density for low (high) equilibrium Fermi energy, this prompts us to investigate whether (hot-)CM occurs. We study this by comparing the electron distributions before and after photoexcitation to the absorbed photon density, which equals the initial density of photoexcited carriers *n*_{exc}. CM takes place when the carrier density in the conduction band with respect to equilibrium, *n*_{CB} − *n*_{CB,0}, is larger than *n*_{exc}. Hot-CM takes place when the hot-carrier density with respect to equilibrium, *n*_{HC} − *n*_{HC,0}, is larger than *n*_{exc}. Our calculations show that CM can take place for low equilibrium DC conductivity, that is, close to the charge neutrality point (see Fig. 1D). Upon increasing the carrier density, CM disappears as a result of suppressed interband heating, whereas hot-CM increases with increasing carrier density, as result of intraband heating.

### Photoconductivity versus photon energy

Our optical pump–THz probe experiment (and the same holds for other all-optical experiments) does not give direct insight into (hot-)CM numbers because we do not harvest electrons through an electrical measurement, which would allow for counting electrons and comparing with absorbed photons. However, a key indication of the occurrence of efficient carrier heating, as well as (hot-)CM, is that a larger photon energy leads to a larger density of secondary carriers (hot carriers) (*19*, *25*, *42*). The reason for this is that additional energy from photons is transduced to additional excited carriers in the electronic system. Figures 3 and 4 explore this trend in our experiments. The results in Fig. 3 correspond to the device tuned away from the charge neutrality point, with an estimated Fermi energy |*E*_{F}| ≃ 0.25 eV (the ionic gate is at 0 V, where we measure σ_{0} ≈ 8 *e*^{2}/h). Here, interband heating is suppressed, whereas intraband heating takes place. Figure 4, instead, corresponds to 0.5 V on the ionic gate, where we measure σ_{0} ≈ 2 *e*^{2}/h. Here, the carrier density is low, and interband heating is present.

Figure 3A shows the measured (peak) THz photoconductivity as a function of pump pulse fluence, parametrized by the photoexcited carrier density, where wavelength-dependent absorption is taken into account. The results for two pump photon energies of 1.0 and 2.5 eV are shown. The larger photon energy leads to a larger negative photoconductivity at constant excitation density, implying a larger hot-carrier density due to efficient intraband scattering. These observations confirm earlier results which showed linear scaling with photon energy, at very low fluence (*25*). Figure 3B shows the measured and the theoretically calculated THz photoconductivity as a function of photon energy, at fixed fluence. The experimental and theoretical results are in good agreement (except for a constant scaling factor) and show a larger negative photoconductivity for larger photon energy. This effect is due to the higher electron temperature at larger photon energy. Because of saturation effects, the scaling with photon energy is not perfectly linear. This saturation likely comes from the reduced heating efficiency for increasing fluence and increasing photon energy, as has been observed in the study of Jensen *et al.* (*30*). We will address this in more detail in Fig. 5.

Figure 3C shows the calculated hot-CM factor, defined as HCM = (*n*_{HC} − *n*_{HC,0})/*n*_{exc}, as a function of photon energy, obtained from our theoretical approach. For the investigated conditions, calculations show that the HCM factor is larger than 1 and increases with photon energy. For a photon energy of 2.5 eV and a fluence corresponding to a photoexcited carrier density of 0.2 × 10^{12}/cm^{2}, we calculate that seven hot carriers are created on average per absorbed photon. This is fully compatible with the very recent experimental demonstration of five generated and subsequently collected hot carriers per absorbed photon in a photo-Nernst device (*44*). The plots of the electron distribution in Fig. 3D illustrate the microscopic picture explaining hot-CM: Excitation with a larger photon energy leads to a stronger broadening of the Fermi-Dirac distribution.

The panels in Fig. 4 show the same quantities as in the corresponding panels of Fig. 3 but closer to the charge neutrality point (see also fig. S2). Because of electron-hole puddles, the effective Fermi energy will not reach a value below 0.05 to 0.1 eV. In Fig. 4 (A and B), we show the THz photoconductivity, measured for three photon energies *E*_{ph} = 0.5, 1.0, and 2.5 eV, as we vary the fluence of the pump pulse, compared to the results of the theoretical calculation. In the calculation, we used *E*_{F} = 0.05 eV. The THz photoconductivity is now positive, with a larger signal for a larger photon energy. This suggests that interband heating is also an efficient process. We again observe saturation effects that are likely related to the decreased heating efficiency at higher fluence and photon energy [see the study of Jensen *et al.* (*30*) and Fig. 5]. An additional effect giving rise to saturation could be a reduction of conductivity (leading to reduced positive photoconductivity) due to populating optical phonons at increasing fluence and photon energy (*27*, *32*).

Analogous to the case of intraband heating and hot-CM, in Fig. 4C, we show the calculated CM factor CM = (*n*_{CB} − *n*_{CB,0})/*n*_{exc} as a function of the photon energy. Under the investigated conditions, we find that the CM factor is larger than 1 and increases with photon energy. For a photon energy of 2.5 eV and a fluence corresponding to a photoexcited carrier density of 0.2 × 10^{12}/cm^{2}, we calculate that >3 secondary carriers are created on average per absorbed photon. The plots of the carrier distribution in Fig. 4D show stronger broadening of the Fermi-Dirac distribution for higher photon energy, corresponding to more electronic transitions from valence to conduction band. Thus, although the occurrence of CM arises from specific microscopic interactions, such as collinear scattering events (*12*, *37*, *39*), the thermodynamic picture of Fig. 4D indicates that it can be regarded as interband heating with high efficiency.

### Heating efficiency

Ultimately, both interband and intraband heating are mechanisms that slow down the dissipation of the electron energy density by phonons, and in this sense, they have the same potential to increase the efficiency of a light-harvesting or photodetection device. To assess the ultimate potential of graphene-based photodetection, we characterize the efficiency of carrier heating, that is, the energy transfer from photons to electrons, by calculating the carrier temperature at time *t* = 0.1 ps and the ratio between the electron energy density at time *t* = 0.1 ps and the photon energy density of the pump pulse. We chose this delay time because it corresponds to the quasi-equilibrium state, where the initial ultrafast energy relaxation, including interband and intraband heating, has already occurred, whereas equilibration with the phonon system is far from complete, that is, the hot carrier distribution has been formed and has not cooled yet. Figure 5 shows these results for a photoexcited carrier density of 0.5 × 10^{12}/cm^{2} and 0.17 × 10^{12}/cm^{2} and 800 nm light. The calculations show that the quasi-equilibrium electron temperature is somewhat higher for lower carrier density (see Fig. 5A), mainly due to the smaller electronic heat capacity. We also see that the higher fluence corresponds to a lower heating efficiency (see Fig. 5B), in agreement with the saturation effects in our measured photoconductivity versus fluence (see Figs. 3 and 4) and with the optical pump–THz probe study of Jensen *et al.* (*30*). Furthermore, time-resolved ARPES (*16*, *17*) and ultrafast photoluminescence (*10*) studies showed less efficient heating because these experiments were done at a fluence that was typically two orders of magnitude larger than the fluence in our THz study. The reason for the reduced heating efficiency is that, at a higher fluence, the energy transfer to phonons occurs faster (see Fig. 5, C and D). This thermodynamic picture shows that the energy transfer from photons to electrons is exceptionally efficient (for the fluences used here) and depends only weakly on the equilibrium Fermi energy. This once more asserts the great intrinsic potential of graphene for photodetection applications, provided that hot carriers can be extracted before their cooling occurs (typically on a picosecond time scale at room temperature).

## DISCUSSION

We have observed a correlation between, on the one hand, the transition from positive to negative THz photoconductivity and, on the other hand, the transition from interband to intraband heating, respectively. In the following, we discuss the physical origins behind positive and negative THz photoconductivity, which in our physical picture can be traced back solely to electronic effects and, in particular, to impurity screening by many-body interactions. There are two reasons to expect an increase of conductivity following photoexcitation (positive THz photoconductivity). First, the pump laser increases the electron density in the conduction band and generates a finite hole density in the valence band. The free carrier density in the system is larger, and thus, the conductivity increases. Second, the generation of hot electrons entails an increase of the electron temperature. Picturing graphene as a zero-gap semiconductor, it is natural to expect that a higher temperature *T* ≃ *E*_{F}/*k*_{B} corresponds to a larger conductivity as well. The reduction of conductivity (negative THz photoconductivity) is more subtle because its origin resides in the reduction of electronic screening of the long-range Coulomb interaction between graphene’s carriers and charged impurities in the substrate. More precisely, the polarizability function of graphene’s carriers, evaluated at the Fermi wave vector where most scattering events take place, decreases with temperature and thus reduces the effectiveness of electronic screening (*55*). When screening is less effective, impurities in the substrate interact more strongly with graphene’s carriers, and the conductivity decreases because of enhanced scattering.

The experimentally observed correlation between the transition from interband to intraband heating and from positive to negative THz photoconductivity can thus be explained as follows: In the vicinity of the charge neutrality point, where interband heating takes place, a larger free carrier density is available for conduction, and this preempts the effect of the reduced screening at higher electron temperature. The THz photoconductivity is positive. At higher carrier density, where intraband heating takes place, the system evolves to a regime where screening of the impurities in the substrate is less effective and the THz conductivity is decreased with respect to its equilibrium value. The THz photoconductivity is negative. Thus, electronic screening plays a crucial role—first, during the ultrafast dynamics that lead to a thermalized electron distribution (*12*, *53*) and, second, for the resulting conductivity that corresponds to this thermalized distribution. Both CM and interband recombination due to carrier-carrier scattering vanish when a nonregularized dynamical screening model is used (*12*, *53*). Our calculations show that such a screening model yields positive THz photoconductivity at all equilibrium carrier densities, contrary to the experimental results, because it overestimates the free carrier density in the conduction band at a large equilibrium carrier density.

To put our results in perspective, we point out that, for semiconductors, the concept of CM is highly attractive because it allows to overcome the Shockley-Queisser limit of conventional solar cells (*56*) and photodetection with gain, where multiple electrons are detected per absorbed photon. Typical light harvesting and photodetection devices using semiconductors rely on the photovoltaic effect, where photoexcited electron-hole pairs are separated by the intrinsic electric field present at the junction between *p*-doped and *n*-doped regions. However, in graphene, the photoresponse at a *pn*-junction (for example, created by local gates) is typically dominated by the photo-thermoelectric (PTE) effect (*57*, *58*), at least under zero-bias voltage. Here, the different Seebeck coefficients in the *p* and *n* regions give rise to an electron-heat driven photoresponse. It turns out that, in this case, efficient heating (and hot-CM) is highly beneficial for creating a large PTE response, and indeed harvesting of multiple hot carriers per absorbed photon has been demonstrated very recently in a high-responsivity photodetector (*44*). Furthermore, in a recent study, an increased photoresponse at a graphene *pn*-junction photodetector (with a 3-dB bandwidth of 65 GHz) was observed upon increasing the bias voltage, which only occurred at low Fermi energy (*59*). This observation was interpreted as arising from a photoconductive contribution to the photoresponse, coming from additional carriers in the conduction band, thus increasing the device conductivity. Our results explain why this effect only occurs at low Fermi energy: Only in this regime does photoexcitation actually lead to additional carriers in the conduction band, that is, *n*_{CB} > *n*_{CB,0} (see Fig. 1D). By contrast, at higher Fermi energies *n*_{CB} ≈ *n*_{CB,0}, so that no photoconductive effect is expected. The insights from our work lead to direct input for optimizing graphene-based photodetector devices.

## METHODS

### Experimental design

For the gate-tunable device, we used CVD graphene (commercially grown by Graphene Supermarket), which was transferred onto a quartz substrate. Quartz was used as opposed to the more common silicon substrate for its negligible photoconductivity and high THz transmission. The ionic gate consisted of LiClO_{4} mixed with poly ethyl oxide in an 8:1 ratio by mass. This mixture was dissolved in methanol to apply it to the graphene. The source and drain contacts were made from 100 nm of Au and 5 nm of Cr. Turbostratic graphene was obtained on the C face of SiC through thermal decomposition in argon atmosphere (see fig. S2 and caption for details). The optical pump–THz probe setups that were used to measure THz photoconductivity were very similar to the one described in an earlier study (*25*), with the pump path containing an optical parametric amplifier and mixing stages (Light Conversion, TOPAS) to convert the incident 800 nm light to different wavelengths. We measured the pump-induced change in THz transmission at the peak of the pump-probe time trace (see fig. S1A), with the 800-nm sampling beam (which measures the THz pulse through electro-optic sampling) positioned at the peak of the quasi-single cycle THz pulse. In this configuration, the pump-induced change in THz absorption was measured, corresponding to the pump-induced change in the real part of the conductivity, as extracted using the thin-film approximation (*25*).

### Calculation of the THz conductivity

The expression for the intraband conductivity σ(*t*) in graphene at angular frequency ω reads (*55*)(1)where − *e* is the electron charge, *v*_{F} the Fermi velocity in graphene, *ν*(ε) = 2|ε|/[π(*ℏv*_{F})^{2}] is the density of states, λ = 0, 1 labels electron states in the conduction band or hole states in the valence band, and *f*_{λ}(ε; *t*) is the corresponding quasi-equilibrium distribution function at energy ε > 0 (measured from the Dirac point) and time *t*. The quantity *f*_{λ}(ε; 0) represents the equilibrium distribution before photoexcitation. The scattering mechanism responsible for the finite conductivity determines the transport scattering time τ(ε; *t*), which depends on the distribution function. Equation 1 is justified as long as *ℏ*ω ≪ 2*E*_{F}, such that the probe pulse does not produce interband electron transitions. The expression for the transport scattering time due to scattering with long-range Coulomb impurities reads (*60*)(2)where *n*_{i} is the impurity density, is the polar angle of the wave vector * k*′,

*q*is the modulus of the difference between

**′ and , and ε**

*k*_{k}=

*ℏv*

_{F}

*is the dispersion of electrons and holes. We denote by*

**k***W*

_{λ}(

*q*;

*t*) the screened Coulomb potential that reads(3)where is the average dielectric constant around the graphene plane, ϵ(

*q*, 0;

*t*) is the static electronic dielectric function in graphene, and the Coulomb impurities have been taken to sit at the graphene plane for simplicity. In the well-known random phase approximation (RPA), the dielectric function is related to the polarizability function χ

^{(0)}(

*q*, ω;

*t*) of the noninteracting electron system (that is, the Lindhard function) by ϵ(

*q*, ω;

*t*) = 1 −

*v*

_{q}χ

^{(0)}(

*q*, ω;

*t*). The expression for the Lindhard function of the thermalized photoexcited electron system in graphene has been obtained in the study of Tomadin

*et al.*(

*53*). The Lindhard function depends on the electron temperature

*T*(

*t*) and the chemical potentials μ

_{e}(

*t*), μ

_{h}(

*t*) for electrons in the conduction band and holes in the valence band at time

*t*. In our approach, we obtained these quantities by integrating in time the semiclassical Boltzmann equation for the coupled electron-optical phonon system, as detailed in the studies of Brida

*et al.*(

*12*) and Tomadin

*et al.*(

*53*). Then, at each time

*t*of interest, we calculated the transport scattering time on an energy mesh with Eq. 2 and, finally, the conductivity with Eq. 1. We remark that this approach takes screening into account, within the RPA approximation, at each instant in time, without making assumptions on the energy dependence of the transport scattering time.

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/4/5/eaar5313/DC1

fig. S1. Sample characterization.

fig. S2. Intrinsically undoped graphene.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is **not** for commercial advantage and provided the original work is properly cited.

## REFERENCES AND NOTES

**Acknowledgments:**We thank Z. Mics, J. Song, N. Richter, and P. Fantuzzi for fruitful discussions and F. Vialla for graphical support.

**Funding:**This work was supported by the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 696656 “Graphene Flagship,” the EU grants FP7-ICT-2013-613024-GRASP and MoQuaS FET-ICT-2013-10 610449, Fondazione Istituto Italiano di Tecnologia, the Spanish Ministry of Economy and Competitiveness through the Severo Ochoa Programme for Centres of Excellence in R&D (SEV-2015-0522), Fundacio Cellex Barcelona, the Mineco grants Ramon y Cajal (RYC-2012-12281), Plan Nacional (FIS2013-47161-P), the Government of Catalonia (2014-SGR-1535), the European Research Council StG CarbonLight (307806), the German Research Foundation (DFG Priority Program Graphene SPP 1459 and SFB TRR173), and the State Research Centre for Innovative and Emerging Materials and the Graduate School of Excellence Materials Science in Mainz (MAINZ) GSC 266. E.H. acknowledges support from Engineering and Physical Sciences Research Council fellowship EP/K041215/1. K.-J.T. acknowledges support through a Mineco Young Investigator Grant (FIS2014-59639-JIN).

**Author contributions:**A.T. and M.P. performed the calculations based on the Boltzmann equation. S.M.H., E.M.A., and E.H. prepared and characterized the gate-tunable graphene devices and performed the optical pump–THz probe measurements on this device. H.I.W., D.T., M.K., and M.B. performed the optical pump–probe spectroscopy measurements on the turbostratic graphene device. A.C. and C.C. provided the turbostratic graphene device. A.T., F.H.L.K., M.P., and K.-J.T. interpreted the experimental and theoretical results. A.T. and K.-J.T. wrote the paper with input from all authors.

**Competing interests:**The authors declare that they have no competing interests.

**Data and materials availability:**All data needed to evaluate the conclusions in the paper are present in the paper and the Supplementary Materials. Additional data related to this paper may be requested from the authors.

- Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).