Abstract
We used pump-probe spectroscopy to study the dependence of optical pump-induced change of graphene extinction coefficient on polarization. We find that the optical pump-induced change of graphene extinction coefficient is polarization-dependent for light with energy identical to that of the excitation light; this dependence decreases rapidly after optical excitation. However, the change of extinction coefficient is isotropic for light with energy much smaller than that of the excitation light. The spatial distribution of momenta of non-equilibrium carriers accounts for this polarization dependence. These results provide a potential way for manipulating graphene’s transient optical response characteristics for applications in ultrafast photonics.
© 2015 Optical Society of America
1. Introduction
Graphene, one-atom-thick semiconductor with unique linear band structure near Dirac point, exhibits peculiar optical and electrical properties and has recently attracted significant attention [1–8 ]. The fascinating optical properties of graphene (such as its universal optical absorbance [9] and electrostatic gating-modified optical absorbance [10]) make it an ideal material for many photonics applications, including optical sensors [11, 12 ], optical modulators [13–15 ], and photodetectors [16, 17 ]. In addition, extinction coefficient k and refractive index n (i.e., optical constants) are two important and useful optical parameters for accurately describing the optical characteristics of graphene in these photonics devices/structures [11, 13, 14, 16–18 ]. Extinction coefficient k is especially important because most graphene-based optical devices utilize the graphene’s optical absorption response properties [11–13, 16–18 ]. Precise estimation of graphene’s extinction coefficient k is important for designing and optimizing graphene-based photonics and optoelectronics devices.
The extinction coefficient/absorbance of graphene, arising from the optical transition within the linear energy band structure, is related to the energy state occupancy corresponding to the optical transition [1, 4, 19 ]. The optical constants of graphene have been extensively characterized by several groups using spectroscopic ellipsometry [2–4, 19, 20 ]. These reports demonstrated that the extinction coefficient is isotropic in the graphene’s plane [3, 4, 16 ]. However, previous studies concentrated on the steady-state extinction coefficient. It is known that the energy state occupancy is modified in optically excited graphene owing to the injection of photo-excited carriers and subsequent carrier relaxation [7, 21–26 ], changing the extinction coefficient. Recent pump-probe studies reported by Mittendorff and Malic et al show that the initial spatial occupancy of momentum states by photo-excited carriers is anisotropic at optically pumped state [7, 21 ], and this occupancy immediately evolves to fully isotropic with carrier relaxation [7, 21, 27 ]. We can deduce that the transient change of extinction coefficient induced by optical pumping is polarization dependent. However, the polarization dependence of optical pump-induced change of graphene extinction coefficient has not been unambiguously investigated. A clear understanding of optical pump-induced changes of extinction coefficient and the dependence of this change on polarization can help predict graphene’s transient optical response characteristics (such as, transient reflection, transient transmission and transient absorption). This knowledge is likely to be important for developing graphene-based high-speed photonic devices, such as saturable absorbers and optical modulators.
In this work, we performed degenerate and nondegenerate pump-probe measurements to explicitly study the dependence of optical pump-induced changes of graphene extinction coefficient on polarization of incident light. The degenerate pump-probe experiment presented here is similar to that reported by Mittendorff et al, which was used to study distribution of photo-excited carriers in momentum space [7]. During the initial stage of photo-excited carrier relaxation, the dependence of pump-induced change of extinction coefficient on polarization was observed for light with photon energies identical to that of the excitation light. Spatial distributions of non-equilibrium carriers’ momenta determined the dependence of extinction coefficient change on polarization. Our results suggest an efficient way for controlling graphene’s transient extinction coefficient by manipulating the distribution of hot carriers. This work is closely related and in certain respect repeats the work of Mittendorff et al [7].
1.A Relationship between Δk and ΔT/T
For beams that are normally incident onto graphene placed on a substrate, the normalized optical transmission can be written as [28]:
Where n sub is the substrate’s refractive index, N is the number of graphene layers and σ is the graphene’s optical conductivity. The transmittance can also be related to the optical absorbance πα of a single graphene layer [28] as:According to the work of Bruna et al. [1], the extinction coefficient k can be related to πα as:Where n = 3.0 is the graphene’s refractive index and d 1 = 0.34 nm is the graphene’s thickness. From Eqs. (2) and (3) , the normalized optical transmittance can be expressed as:From this equation, the differential transmittance can be calculated as:For a single graphene layer (N = 1) on quartz substrate and incident wavelength of λ = 800 nm, the change Δk of the extinction coefficient can be directly related to the differential transmittance by . Therefore, a fundamental material property of graphene can be obtained from the differential transmittance. We performed pump-probe measurements to study the dependence of optical pump-induced change of Δk on polarization. In this method, the change Δkp-p(tdelay) (obtained from the pump-probe differential signal ΔT(t)/T using Eq. (5)) equals to the convolution of Δk(t) and Gaussian intensity distribution of the probe pulse, i.e., . The parameter τFWHM is the pulse width of the probe beam. Clearly, the dependence of Δkp-p(tdelay) on polarization is determined by Δk(t). Thus, in what follows we study the dependence of Δkp-p(tdelay) on polarization.
2. Experiment and analysis
2.A. Experimental setup
A Ti: sapphire mode-locked laser provided the laser pulses (repetition rate of 1 KHz, central wavelength of 800 nm). The pump-probe configuration used in our experiments is shown in Fig. 1(a) . To reduce the noise arising owing to the scattering of the pump beam, the pump and the probe beams were incident onto the sample at different angles. After passing through a λ/2 plate-Glan Taylor prism-λ/2 plate combination, the 800 nm linearly polarized probe beam was incident normally onto the graphene sample (the probe beam was incident normally onto the sample for avoiding any potential influence of oblique incidence of the probe beam on the Δk(tdelay) measurement [19, 20 ]). The probe beam’s linear polarization orientation (Fig. 1(b)) was controlled by the corresponding λ/2 plate of the combination. For the data presented here, the pump beam was incident onto the sample at an angle of 35° relative to the sample surface normal. It should be noted the measured polarization dependence of Δk is independent on the pump beam’s incident angle. i.e., the is ∼1.4 and 1 for degenerate and non-degenerate measurement, respectively. The pump beam’s linear polarization was controlled by using another λ/2 plate-Glan Taylor prism-λ/2 plate combination.
Two combinations of pump/probe frequency had been used for studying the changes of extinction coefficient induced by optical pumping. The 800 nm and 400 nm pulses were used separately for optical pumping. To reveal the relationship between the dependence of Δk(t) on polarization and the carrier distribution, the dependence of on polarization was compared for two types of optical excitation (λpump = 800 and 400 nm). For pump beams of 400 nm the pulses were obtained by frequency doubling of 800 nm pulses in a BBO crystal. Pump-induced transmittance change was measured using a balanced detector in conjunction with a lock-in amplifier referenced to a 383 Hz mechanically chopped pump. Pulse durations τFWHM at the sample were 300 ± 30 fs and 230 ± 30 fs, obtained by two-photon absorption autocorrelation in ZnSe for 800 nm pulses [29] and in CS2 for 400 nm pulses [30], respectively. The spot size of the probe beam at the sample was approximately 35 μm, the probe beam’s fluence was ∼0.13 mJ/cm2 (the density of carriers was ∼1.2 × 1013 cm−2). The spot size of the pump beam at the sample was ∼80 μm. The measured ΔT/T was converted into Δkp-p using Eq. (5).
We studied a chemically vapor-deposited (CVD) grown graphene monolayer that was transferred onto a quartz substrate. Figure 1(c) shows the typical Raman spectrum at 514 nm, which is similar to graphene monolayers spectra reported in literature (the 2D-peak is located at ∼2695 cm−1 and the G-peak is located at ∼1581 cm−1) [31]. This spectral distribution confirms that the sample studied by us was a graphene monolayer. Measurements performed at different locations on the sample yielded similar results.
2.B. Experimental results
First, we present the experimental results obtained from degenerate pump-probe measurements. Figures 2(a) and 2(b) show the probe polarization dependent Δkp-p time trace for pumping with s- and p-polarized pulses, respectively. Owing to optical absorption bleaching caused by Pauli blocking, optical excitation yielded strongly reduced extinction coefficient around zero-delay time [21, 22, 25, 26, 32 ]. Immediately following the photo-excitation, the change of extinction coefficient, Δk, recovered to zero owing to carrier thermalization and subsequent carrier cooling [21].
As Fig. 2 shows, the Δkp-p around the delay of 0 (i.e., the photo-excitation process and the initial stage of carrier relaxation) is a function of probe polarization. For either S or P pump polarizations the peak Δkp-p is maximal when the probe polarization is parallel to the pump polarization and monotonically decreases as the probe polarization is varied to cross with pump polarization. In other words, the pump-induced Δk is maximal for the light that has the same polarization as the excitation light. With carrier relaxation, this dependence of Δk on polarization decreases rapidly, and no dependence on polarization is observed for long delays (above 0.7 ps).
According to previous reports, the photo-excited carriers are primarily located in the excited state, and the spatial distribution of momenta of these carriers is maximized/vanishes in the direction normal/parallel to the polarization of the excitation light [7, 21, 22, 33 ]. Immediately following the photo-excitation, non-equilibrium carriers are scattered into energetically higher and lower states through two competing processes: carrier-carrier scattering and carrier-phonon scattering. This reduces the excitation state occupancy by the carriers [7, 21, 22, 34 ]. At the same time, this anisotropic distribution centered at the excitation state rapidly evolves into fully isotropic in the entire energy band [22, 27 ]. In energy states away from the excitation state, the distribution of hot carriers is fully isotropic throughout the entire non-equilibrium process [22, 27 ]. Because the optical pump-induced reduction of extinction coefficient indicates the inter-band transition bleaching caused by Pauli blocking, the observed dependence of Δkp-p on polarization originates from anisotropic distribution of non-equilibrium carriers (Fig. 3(a) ). If the distribution of carriers in the optical coupling state is isotropic, the pump-induced Δk should be polarization-independent (Fig. 3(b)).
2.C. Analysis of the dependence of Δk and max(Δkp-p) on polarization
In what follows we analyze the dependence of Δk on polarization, based on the distribution of carriers as reported in literature [7, 21, 22 ]. S polarization direction was set as the angular coordinate in the momentum space. According to the works of Mittendorff and Malic et al. [7, 21 ], the spatial occupancy of momentum states by the carriers photo-excited by light with polarization angle of is proportional to , where is the angle between the momentum vector of the photo-excited carriers and the polarization of excitation light. Clearly, the initial occupancy of the photo-excited carriers is anisotropic and peaks in the direction normal to the polarization of excitation light (Fig. 3(a)) [22]. As discussed above, the change of extinction coefficient caused by Pauli blocking arises owing to the interaction between the non-equilibrium carriers from the excitation light and the carriers that are photo-generated by the studied incident light [7]. If the momentum spatial occupation of nonequilibrium carriers (photo-generated by the excitation light) in the energy state corresponding to the studied wavelength is denoted by , the resulting Δk can be written as:
Only if is a function of the quantity depends on the polarization of light . If the distribution of carriers could be described as (e.g., the initial distribution of carriers photo-excited by the excitation light with polarization angle of ), one would obtain [7].
Because the laser pulses used in the present study were not extremely narrow, the distribution of carriers changes during probing owing to carrier relaxation [22]. As a result, the value of Δk is different for different parts of the probe pulse (Fig. 4 ). In what follows, we use the methods of differential calculus to derive the dependence of optical pump-induced Δkp-p on polarization, given the degenerate pump-probe measurements. As shown in Fig. 4, the probe pulse lags the pump pulse by . We assume the pump polarization and probe polarization are characterized by the angles and , respectively, relative to the S polarization direction. For randomly chosen segments of the probe and pump pulses (closed segments in Fig. 4), the momentum spatial occupation numbers of the carriers generated by the two segments are:
where and and denote the peak radiation light intensity (here, we take only the temporal envelope into account) of the pump and probe pulses, respectively. When these non-equilibrium carriers (photo-excited by the pump pulse segment at , i.e. ) interact with the probe pulse segment, the carriers undergo relaxation for the duration of .Former pump-probe experiments showed that the relaxation of the population of carriers could be satisfactorily described by a bi-exponential function [8, 34 ]. We used the functional form of for describing the population relaxation of those carriers photo-excited by the pump pulse segment. Here, τ1 = 110 ± 40 fs and τ2 = 900 ± 200 fs are the two decay times [8, 23 ] and A is a free parameter used to describe the weight of the two relaxation pathways. Same as for silicon [35], the relaxation of anisotropy in the distribution of carriers could be assumed to be described by a single exponential function . The parameter τaniso is the anisotropy decay time. Because this relaxation function is used for describing the ratio distribution of non-equilibrium carriers for different momenta vectors, the term is used for guaranteeing that the integration over does not change after carrier relaxation. Thus, when these carriers interact with the probe pulse segment at time , the distribution of the carriers in the excitation state can be written as:
According to Eq. (6), Δk is proportional to . By integrating over all segments of pump and probe pulses, the measured Δk p-p can be written as:
where:Thus, the measured peak Δkp-p (i.e., Δkp-p(tdelay = 0)) can be described as a cosine function of the probe polarization (angle ). The measured dependence of peak Δkp-p on polarization, induced by s- and p-polarized pulse pumping, is shown in Figs. 2(c) and 2(d), respectively. As the solid lines demonstrate, this dependence is satisfactorily described by a cosine function. From this fit, we estimated that ≈45 fs. This value broadly agrees with the reported anisotropy relaxation time in ΔT/T (i.e., 150 fs, corresponding to ≈30 fs, exp(−150/30)≈0) [7]. Clearly, the time is much smaller than the pulse width used in our experiment. Figure 3 and the integration in Eq. (9) can help explain why a clear dependence of Δk p-p on polarization could be observed when pulses with duration longer than were used. We suggest that the anisotropy in the distribution of carriers photo-excited by the pump pulse segment could be monitored by the probe pulse segment which is temporally very close to (and possibly overlaps with) this pump pulse segment.
According to Eq. (9), the ratio is always above 1 owing to the non-zero value of N. In principle, the dependence of max(Δkp-p) on polarization could be observed for any pulse. Detailed calculations show that decreases with pulse width and increases with (data not shown). This implies that the effective Δk decreases with the pulse width owing to the rapid relaxation of anisotropy in the carrier distribution. The analysis presented above is only valid for Δk corresponding to the optically pumped states. For Δk corresponding to other energy states, obtaining an explicit expression on the dependence of Δk on polarization is challenged by the complicate carrier relaxation process in graphene band structure [21].
To further study the optical pump-induced Δk and confirm the relationship between the dependence of optical pump-induced Δk on polarization and distribution of non-equilibrium carriers, we performed the measurements with the pump of 400 nm and probe of 800 nm [36]. The non-equilibrium carriers in the optically probed state (0.775 eV) were created by scattering of photo-excited carriers in the energy state of 1.55 eV (pump: 400 nm) [37]. As discussed above, the distribution of these non-equilibrium carriers was isotropic during the non-equilibrium process because the optically probed state was away from the excitation state [22, 27 ], and no dependence of Δk on polarization was expected. The obtained Δkp-p is shown in Fig. 5 vs. the delay time.
Following the photo-excitation, the extinction coefficient decreased owing to carrier thermalization, and Δkp-p reached a peak value when a distribution of hot carriers was established [21, 22 ]. As expected, the Δkp-p time traces for different probe beam polarizations were the same.
Figure 3 shows the schematic of the optical transition and carrier distribution. If the distribution of non-equilibrium carriers in the energy state matching the incident light energy is anisotropic (Fig. 3(a)), the optical pump-induced Δk is a function of light polarization. Conversely, if the distribution is isotropic, the optical pump-induced Δk is polarization-independent. Thus, for avoiding anisotropy of transient absorption in the spectral range, graphene can be excited using shorter-wavelength light.
3. Summary
The dependence of optical pump-induced change of graphene extinction coefficient on polarization was determined by the spatial distribution of momenta of non-equilibrium carriers. For the light with energy identical to that of the excitation light, the optical pump-induced change of graphene extinction coefficient was maximal/minimal for the light with polarization parallel/normal to the excitation light. This polarization dependence decreased rapidly with increasing carrier relaxation. For the light with energy away from that of the excitation light, the optical pump-induced change of graphene extinction coefficient was polarization-independent. Manipulating the excitation light is an efficient way for controlling the polarization characteristics of graphene extinction coefficient. This, in turn, allows manipulating transient optical response characteristics of graphene.
Acknowledgments
The authors would like to acknowledge helpful discussions with Dr. Xiang-Tian Kong. This work was supported by the Chinese National Key Basic Research Special Fund (2011CB922003), National Natural Science Foundation of China (11174159, 11374164, 11304166) and the Fundamental Research Funds for the Central Universities (65145005).
References and links
1. M. Bruna and S. Borini, “Optical constants of graphene layers in the visible range,” Appl. Phys. Lett. 94(3), 031901 (2009). [CrossRef]
2. X. Wang, M. Zhao, and D. D. Nolte, “Optical contrast and clarity of graphene on an arbitrary substrate,” Appl. Phys. Lett. 95(8), 081102 (2009). [CrossRef]
3. J. W. Weber, V. E. Calado, and M. C. M. van de Sanden, “Optical constants of graphene measured by spectroscopic ellipsometry,” Appl. Phys. Lett. 97(9), 091904 (2010). [CrossRef]
4. A. C. Diebold, F. J. Nelson, and V. K. Kamineni, “Spectroscopic Ellipsometry of Nanoscale Materials for Semiconductor Device Applications,” in Ellipsometry at the Nanoscale (Springer, 2013), pp. 557–581.
5. G. Teo, H. Wang, Y. Wu, Z. Guo, J. Zhang, Z. Ni, and Z. Shen, “Visibility study of graphene multilayer structures,” J. Appl. Phys. 103(12), 124302 (2008). [CrossRef]
6. T. Stauber, N. M. R. Peres, and A. K. Geim, “Optical conductivity of graphene in the visible region of the spectrum,” Phys. Rev. B 78(8), 085432 (2008). [CrossRef]
7. M. Mittendorff, T. Winzer, E. Malic, A. Knorr, C. Berger, W. A. de Heer, H. Schneider, M. Helm, and S. Winnerl, “Anisotropy of excitation and relaxation of photogenerated charge carriers in graphene,” Nano Lett. 14(3), 1504–1507 (2014). [CrossRef] [PubMed]
8. I. Gierz, J. C. Petersen, M. Mitrano, C. Cacho, I. C. E. Turcu, E. Springate, A. Stöhr, A. Köhler, U. Starke, and A. Cavalleri, “Snapshots of non-equilibrium Dirac carrier distributions in graphene,” Nat. Mater. 12(12), 1119–1124 (2013). [CrossRef] [PubMed]
9. K. F. Mak, M. Y. Sfeir, Y. Wu, C. H. Lui, J. A. Misewich, and T. F. Heinz, “Measurement of the optical conductivity of graphene,” Phys. Rev. Lett. 101(19), 196405 (2008). [CrossRef] [PubMed]
10. K. F. Mak, L. Ju, F. Wang, and T. F. Heinz, “Optical spectroscopy of graphene: From the far infrared to the ultraviolet,” Solid State Commun. 152(15), 1341–1349 (2012). [CrossRef]
11. S. H. Choi, Y. L. Kim, and K. M. Byun, “Graphene-on-silver substrates for sensitive surface plasmon resonance imaging biosensors,” Opt. Express 19(2), 458–466 (2011). [CrossRef] [PubMed]
12. F. Xing, Z.-B. Liu, Z.-C. Deng, X.-T. Kong, X.-Q. Yan, X.-D. Chen, Q. Ye, C.-P. Zhang, Y.-S. Chen, and J.-G. Tian, “Sensitive real-time monitoring of refractive indexes using a novel graphene-based optical sensor,” Sci Rep 2, 908 (2012). [CrossRef] [PubMed]
13. J.-T. Liu, N.-H. Liu, J. Li, X. Jing Li, and J.-H. Huang, “Enhanced absorption of graphene with one-dimensional photonic crystal,” Appl. Phys. Lett. 101(5), 052104 (2012). [CrossRef]
14. X. Zhu, L. Shi, M. S. Schmidt, A. Boisen, O. Hansen, J. Zi, S. Xiao, and N. A. Mortensen, “Enhanced light-matter interactions in graphene-covered gold nanovoid arrays,” Nano Lett. 13(10), 4690–4696 (2013). [CrossRef] [PubMed]
15. M. Liu, X. Yin, and X. Zhang, “Double-layer graphene optical modulator,” Nano Lett. 12(3), 1482–1485 (2012). [CrossRef] [PubMed]
16. A. Pospischil, M. Humer, M. M. Furchi, D. Bachmann, R. Guider, T. Fromherz, and T. Mueller, “CMOS-compatible graphene photodetector covering all optical communication bands,” Nat. Photonics 7(11), 892–896 (2013). [CrossRef]
17. M. Furchi, A. Urich, A. Pospischil, G. Lilley, K. Unterrainer, H. Detz, P. Klang, A. M. Andrews, W. Schrenk, G. Strasser, and T. Mueller, “Microcavity-integrated graphene photodetector,” Nano Lett. 12(6), 2773–2777 (2012). [CrossRef] [PubMed]
18. J. Zhu, Q. H. Liu, and T. Lin, “Manipulating light absorption of graphene using plasmonic nanoparticles,” Nanoscale 5(17), 7785–7789 (2013). [CrossRef] [PubMed]
19. F. J. Nelson, V. K. Kamineni, T. Zhang, E. S. Comfort, J. U. Lee, and A. C. Diebold, “Optical properties of large-area polycrystalline chemical vapor deposited graphene by spectroscopic ellipsometry,” Appl. Phys. Lett. 97(25), 253110 (2010). [CrossRef]
20. V. G. Kravets, A. N. Grigorenko, R. R. Nair, P. Blake, S. Anissimova, K. S. Novoselov, and A. K. Geim, “Spectroscopic ellipsometry of graphene and an exciton-shifted van Hove peak in absorption,” Phys. Rev. B 81(15), 155413 (2010). [CrossRef]
21. E. Malic, T. Winzer, E. Bobkin, and A. Knorr, “Microscopic theory of absorption and ultrafast many-particle kinetics in graphene,” Phys. Rev. B 84(20), 205406 (2011). [CrossRef]
22. E. Malic, T. Winzer, and A. Knorr, “Efficient orientational carrier relaxation in optically excited graphene,” Appl. Phys. Lett. 101(21), 213110 (2012). [CrossRef]
23. B. A. Ruzicka, S. Wang, J. W. Liu, K. P. Loh, J. Z. Wu, and H. Zhao, “Spatially resolved pump-probe study of single-layer graphene produced by chemical vapor deposition,” Opt. Mater. Express 2(6), 708–716 (2012). [CrossRef]
24. J. Shang, S. Yan, C. Cong, H.-S. Tan, T. Yu, and G. G. Gurzadyan, “Probing near Dirac point electron-phonon interaction in graphene,” Opt. Mater. Express 2(12), 1713–1722 (2012). [CrossRef]
25. R. W. Newson, J. Dean, B. Schmidt, and H. M. van Driel, “Ultrafast carrier kinetics in exfoliated graphene and thin graphite films,” Opt. Express 17(4), 2326–2333 (2009). [CrossRef] [PubMed]
26. G. Xing, H. Guo, X. Zhang, T. C. Sum, and C. H. A. Huan, “The Physics of ultrafast saturable absorption in graphene,” Opt. Express 18(5), 4564–4573 (2010). [CrossRef] [PubMed]
27. X.-Q. Yan, J. Yao, Z.-B. Liu, X. Zhao, X.-D. Chen, C. Gao, W. Xin, Y. Chen, and J.-G. Tian, “Evolution of anisotropic-to-isotropic photoexcited carrier distribution in graphene,” Phys. Rev. B 90(13), 134308 (2014). [CrossRef]
28. J. M. Dawlaty, S. Shivaraman, J. Strait, P. George, M. Chandrashekhar, F. Rana, M. G. Spencer, D. Veksler, and Y. Chen, “Measurement of the optical absorption spectra of epitaxial graphene from terahertz to visible,” Appl. Phys. Lett. 93(13), 131905 (2008). [CrossRef]
29. T. D. Krauss and F. W. Wise, “Femtosecond measurement of nonlinear absorption and refraction in CdS, ZnSe, and ZnS,” Appl. Phys. Lett. 65(14), 1739–1741 (1994). [CrossRef]
30. L. D. Ziegler and X. J. Jordanides, “Two-photon absorption resonance effects in the third-order responses of transparent liquids,” Chem. Phys. Lett. 352(3-4), 270–280 (2002). [CrossRef]
31. Z. Ni, Y. Wang, T. Yu, and Z. Shen, “Raman spectroscopy and imaging of graphene,” Nano Res. 1(4), 273–291 (2008). [CrossRef]
32. J. Shang, T. Yu, J. Lin, and G. G. Gurzadyan, “Ultrafast Electron-Optical Phonon Scattering and Quasiparticle Lifetime in CVD-Grown Graphene,” ACS Nano 5(4), 3278–3283 (2011). [CrossRef] [PubMed]
33. B. Y. Sun and M. W. Wu, “Microscopic theory of ultrafast dynamics of carriers photoexcited by THz and near-infrared linearly polarized laser pulses in graphene,” New J. Phys. 15(8), 083038 (2013). [CrossRef]
34. M. Breusing, S. Kuehn, T. Winzer, E. Malic, F. Milde, N. Severin, J. P. Rabe, C. Ropers, A. Knorr, and T. Elsaesser, “Ultrafast nonequilibrium carrier dynamics in a single graphene layer,” Phys. Rev. B 83(15), 153410 (2011). [CrossRef]
35. A. J. Sabbah and D. M. Riffe, “Femtosecond pump-probe reflectivity study of silicon carrier dynamics,” Phys. Rev. B 66(16), 165217 (2002). [CrossRef]
36. H. Yang, X. Feng, Q. Wang, H. Huang, W. Chen, A. T. S. Wee, and W. Ji, “Giant two-photon absorption in bilayer graphene,” Nano Lett. 11(7), 2622–2627 (2011). [CrossRef] [PubMed]
37. W.-T. Liu, S. W. Wu, P. J. Schuck, M. Salmeron, Y. R. Shen, and F. Wang, “Nonlinear broadband photoluminescence of graphene induced by femtosecond laser irradiation,” Phys. Rev. B 82(8), 081408 (2010). [CrossRef]