Abstract
We develop the model for the terahertz (THz) and infrared (IR) photoconductivity of graphene layers (GLs) at room temperature. The model accounts for the linear GL energy spectrum and the features of the energy relaxation and generation-recombination mechanisms inherent at room temperature, namely, the optical phonon absorption and emission and the Auger interband processes. Using the developed model, we calculate the spectral dependences of the THz and IR photoconductivity of the GLs. We show that the GL photoconductivity can change sign depending on the photon frequency, the GL doping and the dominant mechanism of the carrier momentum relaxation. We also evaluate the responsivity of the THz and IR photodetectors using the GL photoconductivity. The obtained results along with the relevant experimental data might reveal the microscopic processes in GLs, and the developed model could be used for the optimization of the GL-based photodetectors.
© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
The gapless energy carrier spectrum and high carrier mobility in graphene layers (GLs) [1] enable unique optical and transport properties. The GL conductivity dependence on the carrier effective temperature has been analyzed and reported in a number of publications (for example, [2–9]). Due to the specifics of the carrier scattering in the GLs, the variation of the GL DC conductivity under the terahertz (THz) and infrared (IR) irradiation (THz and IR photoconductivity) can exhibit interesting behavior [6, 10]. The response of the carrier system in the GLs caused by the THz and IR radiation enables their utilization in the bolometric photodetectors [11–15]. Extensive studies of the carrier dynamics in the GLs in the response to ultrashort optical pulses (for example, [16–24]) stimulate the creation of novel ultrafast optoelectronic devices. The GL ac (dynamic) conductivity under the optical pumping can become negative in the THz range due to the population inversion [25–28] (see also [16,19]).
The THz and IR photoconductivity in GLs is associated with the interplay of the variations of the carrier effective temperature and density (the quasi-Fermi energy). Depending on the dominating scattering mechanisms and the photon energies, the effective temperature can either exceed or be smaller than the lattice temperature. This can be accompanied by a splitting of the quasi-Fermi levels [25,29]. The interplay of the intra- and interband absorption of the incident THz and IR radiation [30] adds complexity to the processes governing the GL conductivity.
The GL photoconductivity was experimentally studied in a number of works demonstrating the effect of negative photoconductivity [10,31,32]. However, it was also shown that the sign of the THz photoconductivity of the GLs changes with the radiation frequency depending on the environment [32].
Despite a variety of the works aimed to describe and explain the features of the THz and IR photoconductivity, there is still no transparent model of the phenomena. The previous studies mainly were focused on the deeply cooled GLs and related bolometric photodetectors. Due to this, the main energy relaxation and generation-recombination mechanisms were associated with the carrier interaction with acoustic phonons and ambient thermal radiation. However, in the uncooled GLs, the carrier intraband and interband transitions due to the recombination and emission of the optical phonons [33] (see also [10]) and the carrier-carrier interaction (Auger processes) [34, 35] can dominate. The inclusion of these mechanisms is necessary for the explanation of the nontrivial spectral characteristics of the GL photoconductivity. Namely, the observed inversion of the photoconductivity sign requires an adequate but transparent model. In this paper, we develop such a model, which accounts not only for the variations of the carrier temperature but also for the possible splitting of the quasi-Fermi levels associated with the THz and IR radiation absorption. This approach allows us to consider the effect of the THz and IR photoconductivity at room temperatures in more details. We also evaluate the application of this effect for uncooled THz and IR photodetectors. The obtained results can be used for the optimization of these photodetectors and promote a deeper insight into the GL microscopic properties.
We focus on both the intrinsic (or compensated) and p-doped GLs at relatively elevated temperatures (such as room temperature). We assume that the carrier momentum relaxation is associated with impurity scattering, scattering on defects, and and on the acoustic phonons. The carrier scattering on the clustered impurities inside the GLs and the substrate inclusions in the decorated GLs is accounted for as well. The intraband and interband transitions of the carriers due to the emission and absorption of the optical phonons are considered as the crucial mechanism of the energy relaxation and the generation-recombination. Despite the specifics of the Auger generation-recombination in the GLs associated with the carrier dispersion law [34] and the carrier collinear scattering (see [35] and references therein), these generation-recombination mechanisms are also included into the model under consideration
2. GL conductivity dependence on the carrier temperature
Due to relatively high carrier densities in GLs at the room temperature, their pair collisions lead to the establishment of the quasi-Fermi energy distributions, fe(ε) and fh(ε), with the same effective temperature T. Therefore, we use the following formula for the real part of the intraband conductivity of the GL- channel σΩ at the frequency Ω (see [2–4]):
Here it is assumed that the momentum relaxation time as a function of the momentum p is equal to τp = τ0(pvW/T0)l = τ0(T/T0)lξl, where τ0 is its characteristric value, ξ = pvW /T0 is the normalized carrier energy in the GL, vW ≃ 108 cm/s is the characteristic velocity in GLs, and ħ is the Planck constant.If the scattering on the weakly screened charged objects (weakly screened Coulomb, long-scale scattering) prevails, τp ∝ p. This case corresponds to the superscript l in Eq. (1) equal to unity (l = 1). In the case of scattering on the neutral impurities and defects, the strongly screened charged impurities, and the acoustic phonons, τ ∝ p−1 (i.e., l = −1). If several mechanisms contribute to the carrier scattering, τp can be described by a combination the linear and the inverse dependences. For a simplified description, in some papers τp = const is used.
The scattering processes can depend on the presence of large-scale charged defects, clustered impurities, [36,37] and various particles in the GL or near the GL interface with the substrate. The latter can be important in the so-called decorated GLs [38,39].
When the carrier momentum relaxation is determined by the short–range scattering, Eq. (1) yields for the pertinent conductivity σΩ is given by
Here σ00 = (e2T0τ0/πħ2) is the characteristic intraband conductivity (it is equal to the low electric-field conductivity in the case of short-range scatterers). In the intrinsic (or compensated) GLs, the electron-hole symmetry leads to the equality of the pertinent quasi-Fermi energies μe = μh = μ (this does not generally imply that the positions of the quasi-Fermi levels and coincide; such a coincidence only occurs if μ = 0). In this case, Eq. (2) transforms to the following: In the case of the short–range scattering, one can find for the characteristic time τ from Eq. (1) . If the long-scale momentum relaxation is the dominant (l = 1), we obtain from Eq. (1) where is the Fermi-Dirac integral. Equations (3) and (4) show that the variations of both the quasi-Fermi energy μ and the effective temperature T caused by the radiation absorption determine the variation of the conductivity, i.e., the effect of photoconductivity. As demonstrated in the following, the quasi-Fermi energy changes not only due to the carrier density change but due to the heating.3. Generation-recombination and energy balance equations
The equation governing the interband balance of the carriers can be generally presented as
The terms in Eq. (5) correspond to the interband Auger generation-recombination processes and the processes associated with optical-phonon, acoustic–phonon, and radiative transitions (in particular, indirect transitions, for which some selection restrictions are lifted). In the situations under consideration (in particular in the temperature range in question), the characteristic times of the processes associated with acoustic phonons and radiative transitions (as was mentioned in the Introduction) are much longer than those related to the optical-phonon and Auger mechanisms [6,9,33,35]. Due to this, we disregard the terms GAc and GRad in Eq. (5).Due to the fast processes of the optical phonon decay into acoustic phonons and their effective removal, confirmed by high values of the G-layer thermal conductivity [40–42], the optical phonon system in the GLs is assumed to be in equilibrium with the lattice having the temperature T0.
For the rate, GOpt, of the interband and intraband transitions assisted by the optical phonons emission and absorption one can use rather general equations accounting for the transitions matrix elements, the features of carrier energy spectrum in the GLs, and the conservation of the carrier momentum and energy [33]. As shown [29] (see also [11, 43]), at ħω0 ≫ μ, T, T0, where ħω0 ≃ 200 meV is the energy of optical phonons in the GL, the general equations for GOpt [33] can be simplified by separating of the exponential factors and relatively slow varying factors depending on μ and T0. Similar procedure leads to a simplified formula for GAuger. As the result, in the situation under consideration, we present Eq. (5) describing the carrier balance in the conduction and valence bands of the GL in the following form:
The equation governing the energy electron and hole balance is given by:
Here τOpt and are the characteristic recombination and intraband relaxation times associated with the carrier interaction with the optical phonons, τAuger is the Auger recombination time, [2] (for a small μ) and are the GL interband and intraband absorption coefficients (determined by the interband and intraband ac conductivities), respectively, Σ0 is the characteristic carrier density, α = (e2/cħ) ≃ 1/137 is the fine structure constant, c is the speed of light in vacuum, and κ is the background dielectric constant. The factors Σ0/τOpt, , and Σ0/τAuger are weak functions of μ and T (mentioned above, see, for example, [29]). These factor can be considered as phenomenological parameters of the model.The exponential terms in the left-hand sides of Eqs. (6) and (7) are proportional to the rate of the generation-recombination and the rate of the energy transfer into the lattice, respectively. The terms in the right-hand sides of these equations correspond to the carrier photogeneration and to the power obtained by the electron-hole system due to the radiation absorption. The latter is proportional to (βintra + βinter)ħΩ IΩ.
Limiting our consideration by low radiation intensities IΩ, Eqs. (6) and (7) are presented as
where, due to a smallness of IΩ, we put and and . Here D0 = (4T0τ0/πħ) is the Drude factor. At T0 = 25 meV and τ0 = (10−12 – 10−13) s, one obtains D0 ≃ 5.3 – 52.7.4. GL photoconductivity
Using Eqs.(8) and (9), we express the normalized carrier quasi-Fermi energy μ/T and temperature variation (T − T0)/T0 via the radiation intensity IΩ, its frequency Ω, and the characteristic times τOpt, , and τAuger as
Here is the characteristic relaxation-recombination time, which accounts for the interband transitions associated with the optical phonons and the Auger generation recombination processes, as well as the intraband energy relaxation on the optical phonons, and b = τOpt/τAuger. The parameters a and b can be estimated using the data found previously [29,33] and [33,35], respectively. Equations (10) and (11) can be rewritten as where The first term in the square brackets of Eq. (12) is negative when Ω < ω0/(1 + b). In this case, the interband transitions promote a decrease in the carrier temperature. The point is that every act of the photon absorption brings the energy ħΩ to the carrier system, while the interband emission of optical phonon decreases the carrier system energy by the value ħω0, which is a smaller value. If b ≪ 1, these processes lead to the carrier heating. The Auger recombination (in the gapless GLs) does not directly change the carrier energy, but decreases the cooling role of the optical photon emission.In the limit of very strong Auger processes (τAuger ≪ τOpt, i.e., b ≫ 1), Eq. (12) yields
In the latter case, the radiative interband processes supply the energy, while the interband optical phonon emission does not work (only the intraband optical phonon processes contribute to the carrier energy relaxation). In this limit, we obtain from Eq. (13) μ = 0. The latter implies that the electron and hole quasi-Fermi levels coincide ().Figure 1 shows the carrier temperature variation (T − T0)/T0 and the quasi-Fermi energy μ/T normalized by IΩ/ĪΩ (i.e., the quantities [(T − T0)/T0](ĪΩ/IΩ) and (μ/T)(ĪΩ/IΩ)), versus the photon energy ħΩ in the GLs with dominating short-range scattering calculated using Eqs. (12) and (13) for different b (i.e., different τAuger) and different τ0. We set ħω0 = 200 meV, T0 = 25 meV, a ≃ π2(T0/ħω0)2(1 + 2.19T0/ħω0) ≃ 0.2 [24].
As seen from Fig. 1 that the photon energies where μ and (T − T0) change their signs do not coincide.
Using Eqs. (3) and (13) when |μ| ≪ 2T, for the difference ΔσΩ = σΩ→0 − σ00 between the dc conductivity under irradiation and the dc dark conductivity we obtain
so thatFigure 2 shows the GL photoconductivity (ΔσΩ/σ00) normalized by a factor IΩ/ĪΩ) calculated using Eq. (16) for different parameters. The GL photoconductivity as a function of the photon energy changes its sign twice at ħΩ = ħΩ0 and ħΩ = ħΩ1 which correspond to the points, where μ/T = εe = εh = 0. One can see that the absolute value of the GL photoconductivity steeply decreases with increasing b, i.e., with the decreasing Auger recombination time τAuger (and, hence, decreasing τR). If τAuger tends to zero, the quantities μ/T and Δσ0 also tend to zero, despite (T −T0) ≠ 0, as follows from Eqs. (13) and (16). All this is because at very intensive Auger processes the splitting of the quasi-Fermi levels vanishes, and the factor in the square brackets in the right-side of Eq. (2) is equal to unity. Hence, the conductivity becomes independent of the carrier temperature.
As seen from Fig. 2, the GL photoconductivity is negative ΔσΩ < 0 in the photon energy ranges 0 < ħΩ < ħΩ0 and ħΩ1 > ħΩ. In the intermediate range ħΩ0 < ħΩ < ħΩ1, the GL photoconductivity is positive.
Taking into account that in reality D0 > 4T0(1 + a)/ħω0, i.e., ω0τ0 > π(1 + a), from Eq. (16) we obtain
andFor τ0 = (10−12 – 10−13) s−1 and a ≃ 0.2, from Eqs. (17) and (18) one obtains Ω0/2π ≃ (1.35 – 3.16) THz and Ω1/2π ≃ (59.29 – 60.92) THz. Thus, the frequencies, at which the GL photoconductivity changes sign, correspond to the THz range (Ω0/2π) and to the Mid-IR range (Ω1/2π).
It should be noted that if the momentum relaxation time τp is independent of the momentum p (as was assumed in some GL conductivity models [10, 44, 45]), the conductivty variation Δσ0/σ00 ≃ μ/4T [46] (compare with Eq. (15). Therefore, the results for the case of τp = const are qualitatively close to those obtained in this section.
5. Clustered impurities and decorated GLs - long-range scattering
In this case, using Eq. (4), at relatively low intensities of the incident radiation IΩ when the deviation of μ from zero is small, we obtain
Hence, The latter is quite different from Eq. (15), which is valid in the case of short-range momentum relaxation. Equation (19) with an equation similar to Eq. (11) results in (compare with Eq. (14))Figure 3 shows the normalized responsivity of the photodetectors based on the GLs with the dominant long-range scattering calculated using Eq. (21).
As follows from the comparison of Eqs. (16) and (21) (see also Figs. 2 and 3), the photoconductivities involving the short- and long-range scattering mechanisms exhibit fairly different spectral dependences: the photoconductivity at long-range scattering is positive at low and high photon energy and negative in the range of the intermediate photon energies. This is because in the case of the dominant short-range scattering, the variation of the dc conductivity stimulated by the interband and intraband radiative transitions is associated primarily with the variation of μ (more exactly with the variation of μ/T), i.e., with the shifts of the electron and hole quasi-Fermi levels. The latter is in contrast with the case of the dominant long-scale scattering in which the variation of the dc conductivity is due to the variation of the carrier effective temperature. This can be seen from the comparison of Eqs. (3) and (19) describing a decrease in the conductivity with increasing −μ/T and an increase with increasing T, respectively.
6. Responsivity of the GL-based photodetectors
Using the obtained values of the photoconductivity, ΔσΩ, of the GLs one can find the responsivity RΩ = |ΔσΩ|E/ħΩ IΩ of the photodetectors based on undoped GL. Using Eqs. (16) and (21), for the cases of the dominant short- and long-scale scattering we arrive at
respectively, where Here E is the electric field along the GL plane and L is the spacing between the contacts. For example, setting τ0 = 10−12 s, τOpt = 10−12 s, κ = 4, L = 10−4 cm, and E = 100 V/cm, from Eq. (24) we obtain R ≃ 2.47 × 10−4 A/W. One needs to note that the expression for RΩ includes a large parameter D0. In particular, for the responsivity of the photodetector using high quality GLs (τ0 = 10−11 s) at ħΩ = 10 meV (i.e., Ω/2π ≃ 2.5 THz) Eq. (21) yields RΩ ≃ 1.8 A/W.Figures 4 and 5 show the RΩ/R versus ħΩ dependences calculated using Eqs. (12) and (22) for the GL-photodetectors with different τ0 and b at a = 0.2. As the GL photoconductivity, the detector responsivity RΩ also turns to zero at certain frequencies.
The responsivity of the GL-based photodetectors is proportional to the fundamental constant πα ≃ 0.023, which determines the coupling of the carriers and the incident radiation. The inclusion of the plasmonic effects associated with the excitation of the plasmons in the optical couplers (metal gratings and so on) and in the GLs can lead to a substantial reinforcement of the carrier-radiation coupling (see, for example [47] and references therein). As a result, the GL-detector responsivity can be markedly enhanced. The responsivity spectral characteristic can exhibit pronounced plasmonic resonances superimposed on the spectral dependences obtained above.
7. Heavily doped GLs
Since the pristine GLs are usually p-type, we will consider the GLs with sufficiently high acceptor and hole densities ΣA. In the limit of short Auger recombination time, μe = −μh and the variations of the carrier temperature do not lead to the variations of the GL conductivity and the effect of the photoconductivity vanishes. In a more realistic situation one can to assume the long recombination time τAuger, therefore, the term with τAuger can be omitted. In this case, we can use Eqs. (6) and (7) but with 2μ replaced by (μe + μh). Apart from this, the quantities and should be replaced by
and respectively, and instead of Σ0 one needs to put Σ = ΣA ≫ Σ0. As a result, we arrive at the following equations: From Eqs. (25) and (26), we obtainThe hole Fermi energies μe and μh are related to the effective temperature T as
Here U ∝ ΣA = const.Taking into account that at μ ≫ T, ℱ1(x) ≃ (x2/2 + π2/6) + x ln(1 + exp(−|x|), at T ∼ T0 from Eq. (28) for μh = μ we obtain
The screening length of the charged scatterers in the doped GLs with μ ≫ T, is much shorter than that in the intrinsic GLs. Therefore, it is natural to assume that in the former GLs the short-range scattering dominates. Considering Eq. (2) with the carrier temperature and the quasi-Fermi energies given by Eqs. (27), (28), and (30), we get
If μ0 > ħΩ1/2 ≃ (1 + a)ħω0/2, Eq. (15) yields ΔσΩ < 0 at all the photon energies. This is in line with the experimental results for the p-GLs with ΣA = (2 – 11) × 1012 cm−2, i.e., with μ0 = 192 – 367 meV [10].Equation (31) resembles Eq. (16). However, the latter includes a small factor [exp(−μ0/T0)/ΣA], different frequency-dependent factor of the interband absorption coefficient, and much shorter . Due to the former, the responsivity of the photodetectors based on doped GLs is smaller than that of the intrinsic GLs. The factor [exp(−μ0/T0)/ΣA] reflects the facts that the GL conductivity is a rather weak function of the carrier temperature at μ0 ≫ T0 and that heat received by the carrier system from the absorbed radiation is distributed among a larger number of the carriers (ΣA ≫ Σ0).
8. Conclusions
We developed an analytical model describing the THz and IR photoconductivity in the GLs at room temperature. Using this model we showed that the energy relaxation and the generation-recombination associated with the optical phonons and the Auger processes are crucial for determining the behavior of GL photoconductivity. We found that in the GLs with dominating short-range scattering mechanism of the carrier momentum relaxation the photoconductivity of the intrinsic GLs changes sign at a certain photon frequency Ω0 in the THz range. Depending on the GL parameters (the transition from the negative to positive photoconductivity) occurs at a certain photon frequency Ω1 in the mid-IR range (the photoconductivity is negative at Ω < Ω0 and Ω > Ω1, while it is positive in the range Ω0 < Ω < Ω1. In contrast, in the case of the long-range scattering dominances, the photoconductivity is positive at Ω < Ω0 and Ω > Ω1 and negative at Ω0 < Ω1. The change in the GL environment might result in the change from the short-range to long-range scattering. The evaluated responsivity of the GL-based THz and IR photodetectors indicates their prospects for different applications.
Funding
Japan Society for Promotion of Science (16H06361); Russian Science Foundation (14-29-00277); Russian Foundation for Basic Research (16-37-60110, 18-07-01379); RIEC Nation-Wide Collaborative Research Project, Japan; Office of Naval Research (Project Monitor Dr. Paul Maki).
Acknowledgments
The authors are grateful to A. Arsenin, V. E. Karasik, V. G. Leiman, and P. P. Maltsev for useful discussions.
Disclosures
The authors declare that there are no conflicts of interest related to this article.
References
1. A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, “The electronic properties of graphene,” Rev. Mod. Phys. 81(1), 109–162 (2009). [CrossRef]
2. L. A. Falkovsky and A. A. Varlamov, “Space-time dispersion of graphene conductivity,” The European Phys. J. B 56(4), 281–284 (2007). [CrossRef]
3. V. T. Vasko and V. Ryzhii, “Voltage and temperature dependence of conductivity in gated graphene,” Phys. Rev. B 76(23), 233404 (2007). [CrossRef]
4. E. H. Hwang, S. Adam, and S. D. Sarma, “Carrier transport in two-dimensional graphene layers,” Phys. Rev. Lett. 98(18), 186806 (2007). [CrossRef] [PubMed]
5. V. Vyurkov and V. Ryzhii, “Effect of Coulomb scattering on graphene conductivity,” JETP Lett. 88(5), 370–373 (2008). [CrossRef]
6. V. T. Vasko and V. Ryzhii, “Photoconductivity of intrinsic graphene,” Phys. Rev. B 77(19), 195433 (2008). [CrossRef]
7. E. H. Hwang and S. das Sarma, “Acoustic phonon scattering limited carrier mobility in two-dimensional extrinsic graphene,” Phys. Rev. B 77(11), 115449 (2008) [CrossRef]
8. E. H. Hwang and S. das Sarma, “Screening induced temperature dependent transport in 2D graphene,” Phys. Rev. B 79(16), 165404 (2009). [CrossRef]
9. O. G. Balev, V. T. Vasko, and V. Ryzhii, “Carrier heating in intrinsic graphene by a strong dc electric field,” Phys. Rev. B 79(16), 165432 (2009). [CrossRef]
10. J. N. Heyman, J. D. Stein, Z. S. Kaminski, A. R. Banman, A. M. Massari, and J. T. Robinson, “Carrier heating and negative photoconductivity in graphene,” J. Appl. Phys. 117(1), 015101 (2015). [CrossRef]
11. V. Ryzhii, T. Otsuji, M. Ryzhii, N. Ryabova, S. O. Yurchenko, V. Mitin, and M. S. Shur, “Graphene terahertz uncooled bolometers,” J. Phys. D: Appl. Phys. 46, 065102 (2013). [CrossRef]
12. Xu Du, D. E. Prober, H. Vora, and C. Mckitterick, “Graphene-based bolometers,” 2D Mater. 1(1), 1–22 (2014).
13. Qi Han, T. Gao, R. Zhang, Yi Chen, J. Chen, G. Liu, Y. Zhang, Z. Liu, X. Wu, and D. Yu, “Highly sensitive hot electron bolometer based on disordered graphene,” Sci Rep. 3, 3533 (2013). [CrossRef] [PubMed]
14. G. Skoblin, J. Sun, and A. Yurgens, “Graphene bolometer with thermoelectric readout and capacitive coupling to an antenna,” Appl. Phys. Lett. 112(6), 063501 (2018). [CrossRef]
15. Y. Wang, W. Yin, Q. Han, X. Yang, H. Ye, Q. Lv, and D. Yin, “Bolometric effect in a waveguide-integrated graphene photodetector,” Chin. Phys. B 25(11), 118103 (2016). [CrossRef]
16. T. Li, L. Luo, M. Hupalo, J. Zhang, M. C. Tringides, J. Schmalian, and J. Wang, “Femtosecond population inversion and stimulated emission of dense Dirac fermions in graphene,” Phys. Rev. Lett. 108(16), 167401 (2012). [CrossRef] [PubMed]
17. S. Boubanga-Tombet, S. Chan, T. Watanabe, A. Satou, V. Ryzhii, and T. Otsuji, “Ultrafast carrier dynamics and terahertz emission in optically pumped graphene at room temperature,” Phys. Rev. B 85(3), 035443 (2012). [CrossRef]
18. I. Gierz, J. C. Petersen, M. Mitrano, C. Cacho, I. E. Turcu, E. Springate, A. Stohr, A. Kohler, U. Starke, and A. Cavalleri, “Snapshots of nonequilibrium Dirac carrier distributions in graphene,” Nat. Mater. 12(12), 1119–1124 (2013). [CrossRef] [PubMed]
19. T. Watanabe, T. Fukushima, Y. Yabe, S. A. Boubanga-Tombet, A. Satou, A. A. Dubinov, V. Ya. Aleshkin, V. Mitin, V. Ryzhii, and T. Otsuji, “The gain enhancement effect of surface plasmonpolaritons on terahertz stimulated emission in optically pumped monolayer graphene,“ New J. Phys. 15(7), 07503 (2013). [CrossRef]
20. K. J. Tielrooij, J.C.W. Song, S. A. Jensen, A. Centeno, A. Pesquera, A. Z. Elorza, M. Bonn, L. S. Levitov, and F. H.L. Koppens, “Photoexcitation cascade and multiple hot-carrier generation in graphene,” Nat. Phys. 9(4), 248–252 (2013). [CrossRef]
21. E. Gruber, R. A. Wilhelm, R. Petuya, V. Smejkal, R. Kozubek, A. Hierzenberger, B. C. Bayer, I. Aldazabal, A. K. Kazansky, F. Libish, A. V. Krasheninnikov, M. Schleberger, S. Facsko, A. G. Borisov, A. Arnau, and F. Aumayr, “Ultrafast electronic response of graphene to a strong and localized electric field,” Nat. Comm. 7, 13948 (2016). [CrossRef]
22. G. X. Ni, L. Wang, M. D. Goldflam, M. Wagner, Z. Fei, A. S. McLeod, M. K. Liu, F. Keilmann, B. Ozyilmaz, A. H. Castro Neto, J. Hone, M. M. Fogler, and D. N. Basov, “Ultrafast optical switching of infrared plasmon polaritons in high-mobility graphene,” Nat. Photon. 10, 244–247 (2016). [CrossRef]
23. A. Mousavian, B. Lee, A. D. Stickel, and Y.-S. Lee, “Ultrafast photocarrier dynamics in single-layer graphene driven by strong terahertz pulses,” J. Opt. Soc. Am. B 35(6), 1255–1259 (2018). [CrossRef]
24. M. Baudisch, A. Marini, J. D. Cox, T. Zhu, F. Silva, S. Teichmann, M. Massicotte, F. Koppens, L. S. Levitov, F. J. G. de Abajo, and J. Biegert, “Ultrafast nonlinear optical response of Dirac fermions in graphene,” Nat. Commun. 9, 1018 (2018). [CrossRef] [PubMed]
25. V. Ryzhii, M. Ryzhii, and T. Otsuji, “Negative dynamic conductivity of graphene with optical pumping,” J. Appl. Phys. 101(8), 083114 (2007). [CrossRef]
26. V. Ryzhii, M. Ryzhii, A. Satou, T. Otsuji, A. A. Dubinov, and V. Y. Aleshkin, “Feasibility of terahertz lasing in optically pumped epitaxial multiple graphene layer structures,” J. Appl. Phys. 106(8), 084507 (2009). [CrossRef]
27. D. Svintsov, V. Ryzhii, A. Satou, T. Otsuji, and V. Vyurkov, “Carrier-carrier scattering and negative dynamic conductivity in pumped graphene,” Opt. Expess 22(17), 19873–19686 (2014). [CrossRef]
28. D. Svintsov, V. Ryzhii, and T. Otsuji, “Negative dynamic Drude conductivity in pumped graphene,” Appl. Phys. Exp. 7, 115101 (2014). [CrossRef]
29. V. Ryzhii, M. Ryzhii, V. Mitin, A. Satou, and T. Otsuji, “Effect of heating and cooling of photogenerated electron-hole plasma in optically pumped graphene on population inversion,” Jpn. J. Appl. Phys. 50(9), 094001 (2011). [CrossRef]
30. F. T. Vasko, V.V. Mitin, V. Ryzhii, and T. Otsuji, “Interplay of intra- and interband absorption in disordered graphene,” Phys. Rev. B 86(23), 235424 (2012). [CrossRef]
31. S. Zhuang, Y. Chen, Y. Xia, N. Tang, X. Xu, J. Hu, and Z. Chen, “Coexistence of negative photoconductivity and hysteresis in semiconducting graphene,” AIP Advances 6(4), 045214 (2016). [CrossRef]
32. C. J. Docherty, C. T. Lin, H. J. Joyce, R. J. Nicholas, L. M. Hertz, L. J. Li, and M. B. Johnston, “Extreme sensitivity of graphene photoconductivity to environmental gases,” Nat. Comm .3, 1228 (2012). [CrossRef]
33. F. Rana, P. A. George, J. H. Strait, S. Sharavaraman, M. Charasheyhar, and M. G. Spencer, “Carrier recombination and generation rates for intravalley and intervalley phonon scattering in graphene,” Phys. Rev. B 79(11)(11), 115447 (2009). [CrossRef]
34. M. S. Foster and I. L. Aleiner, “Slow imbalance relaxation and thermoelectric transport in graphene, ” Phys. Rev. B 79(8), 085415 (2009). [CrossRef]
35. G. Alymov, V. Vyurkov, V. Ryzhii, A. Satou, and D. Svintsov, “Auger recombination in Dirac materials: A tangle of many-body effects,” Phys. Rev. B 97(20), 205411 (2018). [CrossRef]
36. K. M. McCreary, K. Pi, A. G. Swartz, Wei Han, W. Bao, C. N. Lau, F. Guinea, M. I. Katsnelson, and R. K. Kawakami, “Effect of cluster formation on graphene mobility,” Phys. Rev. B 81(11), 115453 (2010). [CrossRef]
37. N. Sule, S. C. Hagness, and I. Knezevic, “Clustered impurities and carrier transport in supported graphene,” Phys. Rev. B 89(16), 165402 (2014). [CrossRef]
38. T. Stauber, G. Gomez-Santos, and F. Javier Garcia de Abajo, “Extraordinary absorption of decorated undoped graphene,” Phys. Rev. Lett. 112(7), 077401 (2014), [CrossRef]
39. V. Ryzhii, T. Otsuji, M. Ryzhii, V. Mitin, and M.S. Shur, “Effect of indirect interband transitions on terahertz conductivity in “decorated ” graphene bilayer heterostructures,” Lithuanian J. Phys. 55(4), 243–248 (2015).
40. A. A. Balandin, S. Ghosh, W. Bao, I. Calizo, D. Teweldebrhan, F. Miao, and C. N. Lau, “Superior thermal conductivity of single-layer graphene,” Nano Lett. , 8(3), 902–907 (2008). [CrossRef] [PubMed]
41. S. Ghosh, I. Calizo, D. Teweldebrhan, E. P. Pokatilov, D. L Nika, A. A. Balandin, W. Bao, F. Miao, and C. N. Lau, “Extremely high thermal conductivity of graphene: prospects for thermal management applications in nano-electronic circuits,” Appl. Phys. Lett. 92(15), 151911 (2008). [CrossRef]
42. A. A. Balandin, “Thermal properties of graphene and nanostructured carbon materials,” Nat. Mat. 10, 569–581 (2011). [CrossRef]
43. V. Ryzhii, M. Ryzhii, V. Mitin, and T. Otsuji, “Toward the creation of terahertz graphene injection laser,” J. Appl. Phys. 110(9), 094503 (2011). [CrossRef]
44. K. F. Mak, M. Y. Sfeir, Y. Wu, Ch. H. Lui, J. A. Misewich, and T. F. Heinz, “Measurement of the optical conductivity of graphene,” Phys. Rev. Lett. 101(19), 196405 (2008). [CrossRef]
45. H. M. Dong, W. Xu, and F. M. Peters, “Electrical generation of terahertz blackbody radiation from graphene,” Opt. Express 26(19), 24621–24626 (2018). [CrossRef]
46. V. Ryzhii, M. Ryzhii, D. S. Ponomarev, V. G. Leiman, V. Mitin, M. S. Shur, and T. Otsuji, “Negative photoconductivity and hot-carrier bolometric detection of terahertz radiation in graphene-phosphorene hybrid structures,” J. Appl. Phys.124, in press; arXiv:808.09083 [cond-mat.mes-hall].
47. B. Zhao, J. M. Zhao, and Z. M. Zhang, “Resonance enhanced absorption in a graphene monolayer using deep metal gratings,” J. Opt. Soc. Am. B 32(6), 1176–1185 (2018). [CrossRef]