• Chinese Optics Letters
  • Vol. 20, Issue 9, 093701 (2022)
Tingyuan Jia1、2、3, Shaoming Xie1、2, Zeyu Zhang1、2、3、4、*, Qinxue Yin1、3, Chunwei Wang1、2、3、5, Chenjing Quan1、2, Xiao Xing1, Juan Du1、2、3、4、**, and Yuxin Leng1、2、3、5、***
Author Affiliations
  • 1State Key Laboratory of High Field Laser Physics and CAS Center for Excellence in Ultra-intense Laser Science, Shanghai Institute of Optics and Fine Mechanics (SIOM), Chinese Academy of Sciences (CAS), Shanghai 201800, China
  • 2Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
  • 3Hangzhou Institute for Advanced Study, University of Chinese Academy of Sciences, Hangzhou 310024, China
  • 4School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
  • 5School of Physical Science and Technology, ShanghaiTech University, Shanghai 201210, China
  • show less
    DOI: 10.3788/COL202220.093701 Cite this Article Set citation alerts
    Tingyuan Jia, Shaoming Xie, Zeyu Zhang, Qinxue Yin, Chunwei Wang, Chenjing Quan, Xiao Xing, Juan Du, Yuxin Leng. Role of the interlayer interactions in ultrafast terahertz thermal dynamics of bilayer graphene[J]. Chinese Optics Letters, 2022, 20(9): 093701 Copy Citation Text show less

    Abstract

    Bilayer graphene, which is highly promising for electronic and optoelectronic applications because of its strong coupling of the Dirac–Fermions, has been studied extensively for the emergent correlated phenomena with magic-angle manipulation. Due to the low energy linear type band gap dispersion relationship, graphene has drawn an amount of optoelectronic devices applications in the terahertz region. However, the strong interlayer interactions modulated electron-electron and electron-phonon coupling, and their dynamics in bilayer graphene have been rarely studied by terahertz spectroscopy. In this study, the interlayer interaction influence on the electron-electron and the electron-phonon coupling has been assigned with the interaction between the two graphene layers. In the ultrafast cooling process in bilayer graphene, the interlayer interaction could boost the electron-phonon coupling process and oppositely reduce the electron-electron coupling process, which led to the less efficient thermalization process. Furthermore, the electron-electron coupling process is shown to be related with the electron momentum scattering time, which increased vividly in bilayer graphene. Our work could provide new insights into the ultrafast dynamics in bilayer graphene, which is of crucial importance for designing multi-layer graphene-based optoelectronic devices.

    1. Introduction

    Bilayer graphene (BLG) is a kind of graphene composed of two coupled honeycomb-like carbon layers[1]. The strong interaction of the Dirac–Fermions in BLG has driven new insight into the strong correlation phenomena such as superconductor and correlated insulator phases in the magic-angle twisted BLG systems[26]. During the phase diagram investigating the Bardeen–Cooper–Schrieffer (BCS) state and Bose–Einstein condensate (BEC) crossover in the magic-angle graphene system, the large critical temperature (Tc) and Fermi temperature (TF) make the weak electron-phonon (e-p) coupling system in BCS theory not applicable to describe the emergent phenomena in the twist BLG (TBG)[2]. Therefore, it is crucial to investigate the effect of the layer-layer interaction on the electron-electron (e-e) and e-p coupling in BLG in the ultrafast time scale. Kar et al. studied the response of the photoexcited hot carriers to terahertz (THz) fields in low and moderately doped BLG deposited on a quartz substrate, which shows that the relaxation dynamics is dominated by carrier disorder interaction along with carrier acoustic phonon interactions, as compared to charge-impurity Coulomb interaction or carrier surface phonon interaction[7]. Recently, by means of high-pressure treatment, the ultrafast hot electron cooling process in BLG was interpreted as a strong interlayer shear phonon coupling, which could accelerate the relaxation dynamics[8]. However, the impact of strong interlayer coupling on the THz response, which is more sensitive to the electrons near the Fermi energy, had been seldom discussed[9].

    In the case of doped monolayer graphene (MLG), the ultrafast THz photoconductivity is reduced by the intraband carrier-carrier thermalization near the Fermi energy, which is different from the interband transition origin compared with the optical pump optical probe and angle-resolved photoelectron spectroscopy (ARPES)[8,10]. Similarly, THz is sensitive to the strong coupling induced flat band and the possible gap opening process near the Dirac cone, which is associated with the interlayer e-e and e-p interactions[6,11,12]. Therefore, to disentangle the electron dynamics near the Fermi surface of graphene, the visible-infrared pump THz probe spectroscopy (VIPTP) is introduced to investigate the ultrafast thermalization and hot carrier cooling process in both MLG and BLG. The strong interlayer interaction is proved to be able to accelerate the ultrafast cooling dynamics in BLG. However, the subpicosecond thermalization efficiency is found to be reduced in BLG. Furthermore, the carrier scattering time obtained by the Drude fitting indicates that the strong layer-layer coupling in the BLG is probably correlated with the band renormalization process. Therefore, the interlayer strong coupling plays different roles in the thermalization and hot carrier cooling processes in BLG. Specifically, for thermalization, the interlayer interaction-induced less efficient thermal process could be attributed to the carrier-carrier scattering enhancement rather than the interlayer e-p coupling cooling mechanism, which is possible due to the A-A stacking (Turbostratic stacking)[13]. In contrast, in the case of cooling process, the interlayer phonon scattering, which is related to the A-B stacking (Bernal stacking), the interlayer interaction could help to boost the cooling process[8]. Our investigation would not only propose the different physical origin between the interlayer interaction and the intralayer e-p coupling, but also give new insight into the ultrafast investigation of the strong coupling in the BLG system.

    2. Experimental Setup

    The VIPTP is driven by a 1 kHz Ti:sapphire regenerative amplifier with 800 nm central wavelength and 35 fs pulse duration and then through an optical parameter amplifier to change the center wavelength to 190–2600 nm with the beam size of 0.25cm2. The THz pulses are generated by optical rectification and detected by electro-optic sampling in a pair of 1 mm thick, (110)-oriented ZnTe crystals. The MLG and BLG were grown on a SiO2 substrate by chemical vapor deposition (CVD), respectively. As shown in Fig. S1 of Supplementary Material, the static THz transmission spectrum is shown to be 6% absorption for MLG and 8% absorption for BLG. Figure 1(c) shows the typical Raman spectra of MLG and BLG excited by a 532 nm laser, respectively. The 2D peak and G peak are related to the energy band structure of graphene[14]. The G peak position of both samples indicates that MLG and BLG have approximately the same Fermi level[13]. Compared with the MLG, the electronic band structure splits, both the conduction band and valence band are composed of two parabolas, and there are four possible double resonance scattering processes in BLG of the 2D peak, which have the FWHM=66cm1 in BLG and FWHM=45cm1 in MLG[1416]. As shown in Fig. 1(c), the Raman spectrum of MLG features a much higher 2D peak intensity than that of the G peak (I2D/IG2)[17]. Moreover, the peak width of the 2D peak increases, and the intensity ratio I2D/IG decreases in the case of BLG, which reconfirms the bilayer crystal structure.

    (a) Schematic diagram of MLG on SiO2 substrate; (b) schematic diagram of BLG film on SiO2 substrate; (c) Raman spectra of MLG and BLG.

    Figure 1.(a) Schematic diagram of MLG on SiO2 substrate; (b) schematic diagram of BLG film on SiO2 substrate; (c) Raman spectra of MLG and BLG.

    3. Results and Discussion

    The pump-induced photoconductivity (Δσ) is probed by a THz pulse by monitoring the change in the transmitted electrical field (ΔE=EpumpE0) following the optical pump as a function of pump-probe delay. The measurement is based on the principle that ΔE is proportional to Δσ=1+nZ0(11+ΔT/T01)[9], where, n=1.95 is the refractive index of the fused silica substrate, and Z0=377 is the free space impedance. The time-resolved signals Δσ of MLG and BLG films dependent on the different pump fluence under pump photon energy of 0.78 eV, 1.13 eV, and 1.91 eV are shown in Fig. S2 (in Supplementary Material). We use the deconvolution equation of the e-e interaction in BLG with the deconvolution fitting, |Δσ|=A×Conv(etτdecay,τpulse)+B×Conv(etτdecay,τslowrise),where Conv(etτdecay,τpulse) is the convolution of etτdecay with a Gaussian pulse of pulse width τ, τpulse is the pulse width used, and the calculated effective conductivity drop time can be obtained with τrise=(A+B)/[(Aτpulse)+(Bτslowrise)]; this is how long it takes for the initial light excitation to produce a peak negative THz conductivity. As shown in Fig. 2(a), with the Gaussian deconvolution fitting of the thermalization process in BLG and MLG, the e-e thermalization time in BLG is about 0.75 ps; in contrast, the e-e process takes about 0.35 ps, which indicates a higher electron thermalization in BLG with the help of interlayer interaction. For the cooling process, BLG is faster than MLG. The photoexcited THz conductivity increases with increasing pump fluence and photon energy due to the presence of the carrier multiplication process[18]. In contrast to the monolayer one, the carrier-carrier scattering-induced transient conductivity of BLG is reduced, which could be modulated either by the carrier mobility or the e-p coupling efficiency.

    (a) Photoinduced THz conductivity of MLG and BLG as a function of the rising time at 1.91 eV, 160 µJ/cm2; (b) photoinduced THz conductivity of MLG and BLG as a function of the delay time at 1.91 eV, 160 µJ/cm2.

    Figure 2.(a) Photoinduced THz conductivity of MLG and BLG as a function of the rising time at 1.91 eV, 160 µJ/cm2; (b) photoinduced THz conductivity of MLG and BLG as a function of the delay time at 1.91 eV, 160 µJ/cm2.

    In order to disentangle the e-p and optical phonon-acoustic phonon interaction effect between the graphene layers, the mono-exponential fitting is presented in Fig. 3[19]: Δσ(Δt)=A·etτdecay+C,where τdecay is the relaxation dynamics of the lattice temperature. C represents the longitude thermal conductivity diffusion process, which could be regarded as a plateau in our condition. As shown in Figs. 3(a)3(c), the relaxation lifetime of MLG increases with the pump fluence increasing from 4μJ/cm2 to 160μJ/cm2, while in BLG the relaxation time of the hot carrier cooling process increases with the pump fluence at first and then reaches saturation after 40μJ/cm2. For MLG, cooling is slower because electrons first cool via optical phonons and then there can be a hot-phonon bottleneck due to optical-to-acoustic phonon cooling[19]. However, the bottleneck is less severe because there is more phonon coupling in BLG due to the interlayer interaction[20]. As shown in Fig. 3(d), the cooling process of BLG is significantly more efficient than in MLG[11], which indicates that an extra layer-layer interaction-induced cooling channel is possible to affect the cooling process in BLG[21]. The interlayer phonon scattering, which is caused by the interlayer interaction, could help boost the cooling process. Considering that the Fermi energy and the degree of the disorder do not change obviously between these two samples, as represented by the Raman spectra shown in Fig. 1(c)[15], the intrinsic origin such as the carrier-carrier strong interaction, the phonon-phonon interaction assisted interband transition, and the interlayer shear mode acoustic phonon coupling should be taken into account[2,8,10,11,22].

    (a)–(c) Relaxation process of MLG (black) and BLG (red) as a function of pump fluence at different photon energies; (d) relaxation process of MLG (black) and BLG (red) as a function of photon energy at the same pump fluence.

    Figure 3.(a)–(c) Relaxation process of MLG (black) and BLG (red) as a function of pump fluence at different photon energies; (d) relaxation process of MLG (black) and BLG (red) as a function of photon energy at the same pump fluence.

    For the primary concern, the negative THz onductivity arises from the carrier-carrier thermalization process, which reduced the conductivity of the carriers near the Fermi surface; therefore, the measurement of the thermalization efficiency could represent the carrier-carrier Coulomb interaction strength in MLG and BLG[18]. Therefore, we introduce a parameter to describe the degree of sub-linearity. Figures 4(a)4(c) show the THz conductivity peak value |Δσ| of MLG (black) and BLG (red) at 0.78 eV, 1.13 eV, and 1.91 eV, respectively. The peak THz conductivities of MLG and BLG are fitted by the power exponent function. The degree of deviation of the thermalizing sub-linearity η can be fitted from the following formula[23]: Δσ(peak)=AFη,where Δσ(peak) is the peak value of the photo-induced THz conductivity at different pump wavelengths (photo energy), F is the pump fluence, η describes the degree of sub-linearity of the thermalization, and η increases as the thermalization efficiency decreases. The degrees of deviation of the thermalizing sub-linearity are summarized in Table 1.

     η at 0.78 eVη at 1.13 eVη at 1.91 eV
    MLG0.161±0.0430.198±0.0980.199±0.051
    BLG0.256±0.0680.320±0.0950.379±0.053

    Table 1. The Heating Parameter η

    (a)–(c) Optical-induced THz conductivity peak value |Δσ| of MLG (black) and BLG (red) under different photon energies as a function of pump fluence; (d) ultrafast hot carrier cooling process of MLG and ultrafast heating and hot carrier cooling process of BLG. Diagrams of ultrafast processes and relaxation dynamics involving optical pumping (straight arrows), electron scattering (curled arrows), and optical phonon scattering (vertical blue wiggled arrows). Filled (open) circles signify electrons (holes). (e) Schematic diagram of the ultrafast dynamics in graphene after photoexcitation.

    Figure 4.(a)–(c) Optical-induced THz conductivity peak value |Δσ| of MLG (black) and BLG (red) under different photon energies as a function of pump fluence; (d) ultrafast hot carrier cooling process of MLG and ultrafast heating and hot carrier cooling process of BLG. Diagrams of ultrafast processes and relaxation dynamics involving optical pumping (straight arrows), electron scattering (curled arrows), and optical phonon scattering (vertical blue wiggled arrows). Filled (open) circles signify electrons (holes). (e) Schematic diagram of the ultrafast dynamics in graphene after photoexcitation.

    The degree of deviation of the thermalizing sub-linearity η of the hot electron in ideal MLG is predicted to be 0.5 due to the electron capacity C being in direct proportion to the electron temperature T[24]. In the case of CVD graphene with high mobility, the thermal efficiency is found to be dependent on the doping level and thus the Fermi–Dirac distribution near the Fermi surface[25]. As shown in Figs. 4(a)4(c) and Table 1, the sub-linearity index η(BLG) is larger than η(MLG) and monotonically increases with the pump photon energy, which indicates less photo-induced thermal efficiency in BLG. As described in Figs. 4(d) and 4(e), for BLG, the strong interaction between layers may affect the band structure near the Dirac point, resulting in the decrease of the interband transition-induced intraband thermalization efficiency. Therefore, it is reasonable to conclude that the layer-layer interaction would reduce the thermal efficiency and thus the carrier-carrier interaction in the initial subpicosecond thermal process[26]. Moreover, the sublinear photon energy dependence of the sub-linearity index η in graphene demonstrates that the thermal efficiency is still related to the electron-optical phonon coupling, while the thermal efficiency decrease indicates that the origins for the thermal and cooling processes in both graphene samples are different[7,2729].

    Despite the e-p coupling, the strong layer interaction affected layer electron momentum scattering time should be investigated. We measured the 2D photoconductivity of two samples at the pump-probe delay time of 0 and 2 ps, respectively, as shown in Figs. 5(a)5(d). To understand the conductivity dispersion relation for different samples, we fit the experimental conductivity by the Drude model, which is often employed for the transport of charge carriers in graphene[9], Δσ˜(ω)=iΔDπ(ω+iΓ),where Γ is the average scattering rate for momentum changing collisions of charge carriers. The Drude weight ΔD characterizes the magnitude of the reduced conductivity. The carrier scattering times at different pump-probe delays are 0.023±0.059ps and 0.094±0.057ps in MLG, 0.477±0.069ps and 0.835±0.133ps in BLG, respectively. Due to the relationship below, σ=n|e|μe, μe=eτm*, one can find that the Drude conductivity is in direct proportion to the scattering time and the carrier density. As the Fermi energy and the pump fluence are fixed, the time dependent conductivity change should come from the carrier momentum scattering time τ, which is related to the electronic band structure near the Dirac cone[10,26]. Therefore, the thermalization of the Drude-like Dirac–Fermion should be modulated by the A-A stacking of the interlayer tunneling potential and thus the band structure[3,13]. In addition, the e-e scattering time increase could possibly be induced by the asymmetric band structure of BLG[30,31]. The momentum scattering time difference modulated by the interlayer interaction further confirms the different origins of the strong interlayer coupling modulation of the thermalization and the hot carrier cooling process in BLG.

    (a), (b) Extracted frequency dependence of the THz conductivity at delay times of 0 ps and 2 ps for MLG, solid lines show the fit of the complex conductivity to a Drude model; (c), (d) extracted frequency dependence of the THz conductivity at delay times of 0 ps and 2 ps for BLG, solid lines show the fit of the complex conductivity to a Drude model.

    Figure 5.(a), (b) Extracted frequency dependence of the THz conductivity at delay times of 0 ps and 2 ps for MLG, solid lines show the fit of the complex conductivity to a Drude model; (c), (d) extracted frequency dependence of the THz conductivity at delay times of 0 ps and 2 ps for BLG, solid lines show the fit of the complex conductivity to a Drude model.

    4. Conclusions

    In conclusion, the time-resolved THz conductivity dynamics of MLG and BLG have been studied by optical pump THz probe spectroscopy with different photon energies and fluences. By analyzing the thermal and cooling dynamics of the optical-induced negative THz conductivity, we find that the interlayer strong coupling process plays different roles in the ultrafast thermal and relaxation process of THz conductivity. In addition, the interlayer interaction of BLG could accelerate the hot carrier cooling process by means of the phonon-phonon coupling process due to the A-B stacking clusters in the illumination zone of BLG, and the carrier-carrier scattering time of the Dirac–Fermion reduction could be assigned to the band structure modulation induced by the A-A stacking interlayer tunneling process. Finally, the different carrier scattering times changed with the pump delay in both the MLG and BLG are obtained by utilizing the Drude model, and the momentum scattering time in BLG is much larger than that in MLG, indicating that the band structure was modulated by the layer-layer interaction. Our work will provide new insights for the application of ultrafast nonequilibrium heating and cooling pathways, and the improvement of multi-layer graphene-based nano optoelectronic devices[32].

    References

    [1] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, A. K. Geim. The electronic properties of graphene. Rev. Mod. Phys., 81, 109(2009).

    [2] Y. Cao, V. Fatemi, A. Demir, S. Fang, S. L. Tomarken, J. Y. Luo, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras, R. C. Ashoori, P. Jarillo-Herrero. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature, 556, 80(2018).

    [3] J. Gonzalez, T. Stauber. Kohn-Luttinger superconductivity in twisted bilayer graphene. Phys. Rev. Lett., 122, 026801(2019).

    [4] I. Das, X. Lu, J. Herzog-Arbeitman, Z.-D. Song, K. Watanabe, T. Taniguchi, B. A. Bernevig, D. K. Efetov. Symmetry-broken Chern insulators and Rashba-like Landau-level crossings in magic-angle bilayer graphene. Nat. Phys., 17, 710(2021).

    [5] J. M. Park, Y. Cao, K. Watanabe, T. Taniguchi, P. Jarillo-Herrero. Tunable strongly coupled superconductivity in magic-angle twisted trilayer graphene. Nature, 590, 249(2021).

    [6] J. Wu, Y. Zheng, Z. Zeng, R. Li. High-order harmonic generation from zigzag graphene nanoribbons. Chin. Opt. Lett., 18, 103201(2020).

    [7] S. Kar, V. L. Nguyen, D. R. Mohapatra, Y. H. Lee, A. K. Sood. Ultrafast spectral photoresponse of bilayer graphene: optical pump-terahertz probe spectroscopy. ACS Nano, 12, 1785(2018).

    [8] K. Ni, J. Du, J. Yang, S. Xu, X. Cong, N. Shu, K. Zhang, A. Wang, F. Wang, L. Ge, J. Zhao, Y. Qu, K. S. Novoselov, P. Tan, F. Su, Y. Zhu. Stronger interlayer interactions contribute to faster hot carrier cooling of bilayer graphene under pressure. Phys. Rev. Lett., 126, 027402(2021).

    [9] Z. Zhang, T. Lin, X. Xing, X. Lin, X. Meng, Z. Cheng, Z. Jin, G. Ma. Photoexcited terahertz conductivity dynamics of graphene tuned by oxygen-adsorption. Appl. Phys. Lett., 110, 111108(2017).

    [10] S. Ulstrup, J. C. Johannsen, F. Cilento, J. A. Miwa, A. Crepaldi, M. Zacchigna, C. Cacho, R. Chapman, E. Springate, S. Mammadov, F. Fromm, C. Raidel, T. Seyller, F. Parmigiani, M. Grioni, P. D. King, P. Hofmann. Ultrafast dynamics of massive Dirac Fermions in bilayer graphene. Phys. Rev. Lett., 112, 257401(2014).

    [11] M. T. Mihnev, J. R. Tolsma, C. J. Divin, D. Sun, R. Asgari, M. Polini, C. Berger, W. A. de Heer, A. H. MacDonald, T. B. Norris. Electronic cooling via interlayer Coulomb coupling in multilayer epitaxial graphene. Nat. Commun., 6, 8105(2015).

    [12] J. Huang, T. Fu, H. Li, Z. Shou, X. Gao. A reconfigurable terahertz polarization converter based on metal–graphene hybrid metasurface. Chin. Opt. Lett., 18, 013102(2020).

    [13] S. Carr, D. Massatt, S. Fang, P. Cazeaux, M. Luskin, E. Kaxiras. Twistronics: manipulating the electronic properties of two-dimensional layered structures through their twist angle. Phys. Rev. B, 95, 075420(2017).

    [14] J. Wu, H. Xu, J. Zhang. Raman spectroscopy of graphene. Acta Chim. Sin., 72, 301(2014).

    [15] Z. Tao, J. Du, Z. Qi, K. Ni, S. Jiang, Y. Zhu. Raman spectroscopy study of sp2 to sp3 transition in bilayer graphene under high pressures. Appl. Phys. Lett., 116, 133101(2020).

    [16] D. Sun, C. Divin, C. Berger, W. A. de Heer, P. N. First, T. B. Norris. Spectroscopic measurement of interlayer screening in multilayer epitaxial graphene. Phys. Rev. Lett., 104, 136802(2010).

    [17] A. Das, S. Pisana, B. Chakraborty, S. Piscanec, S. K. Saha, U. V. Waghmare, K. S. Novoselov, H. R. Krishnamurthy, A. K. Geim, A. C. Ferrari, A. K. Sood. Monitoring dopants by Raman scattering in an electrochemically top-gated graphene transistor. Nat. Nanotechnol., 3, 210(2008).

    [18] K. J. Tielrooij, J. C. W. Song, S. A. Jensen, A. Centeno, A. Pesquera, A. Zurutuza Elorza, M. Bonn, L. S. Levitov, F. H. L. Koppens. Photoexcitation cascade and multiple hot-carrier generation in graphene. Nat. Phys., 9, 248(2013).

    [19] T. Jia, W. Zhang, Z. Zhang, Z. Zhan, G. Ma, J. Du, Y. Leng. Role of the optical–acoustic phonon interaction in the ultrafast cooling process of CVD graphene. J. Phys. Chem. C, 125, 27283(2021).

    [20] E. A. A. Pogna, X. Jia, A. Principi, A. Block, L. Banszerus, J. Zhang, X. Liu, T. Sohier, S. Forti, K. Soundarapandian, B. Terrés, J. D. Mehew, C. Trovatello, C. Coletti, F. H. L. Koppens, M. Bonn, H. I. Wang, N. van Hulst, M. J. Verstraete, H. Peng, Z. Liu, C. Stampfer, G. Cerullo, K.-J. Tielrooij. Hot-carrier cooling in high-quality graphene is intrinsically limited by optical phonons. ACS Nano, 15, 11285(2021).

    [21] Z. Zhang, M. Hu, T. Jia, J. Du, C. Chen, C. Wang, Z. Liu, T. Shi, J. Tang, Y. Leng. Suppressing the trapping process by interfacial charge extraction in antimony selenide heterojunctions. ACS Energy Lett., 6, 1740(2021).

    [22] P. H. Tan, W. P. Han, W. J. Zhao, Z. H. Wu, K. Chang, H. Wang, Y. F. Wang, N. Bonini, N. Marzari, N. Pugno, G. Savini, A. Lombardo, A. C. Ferrari. The shear mode of multilayer graphene. Nat. Mater., 11, 294(2012).

    [23] H. A. Hafez, S. Kovalev, K. J. Tielrooij, M. Bonn, M. Gensch, D. Turchinovich. Terahertz nonlinear optics of graphene: from saturable absorption to high-harmonics generation. Adv. Opt. Mater., 8, 1900771(2019).

    [24] M. W. Graham, S.-F. Shi, D. C. Ralph, J. Park, P. L. McEuen. Photocurrent measurements of supercollision cooling in graphene. Nat. Phys., 9, 103(2012).

    [25] G. Jnawali, Y. Rao, H. Yan, T. F. Heinz. Observation of a transient decrease in terahertz conductivity of single-layer graphene induced by ultrafast optical excitation. Nano Lett., 13, 524(2013).

    [26] K. J. Tielrooij, L. Piatkowski, M. Massicotte, A. Woessner, Q. Ma, Y. Lee, K. S. Myhro, C. N. Lau, P. Jarillo-Herrero, N. F. van Hulst, F. H. Koppens. Generation of photovoltage in graphene on a femtosecond timescale through efficient carrier heating. Nat. Nanotechnol., 10, 437(2015).

    [27] A. Tomadin, S. M. Hornett, H. I. Wang, E. M. Alexeev, A. Candini, C. Coletti, D. Turchinovich, M. Klaui, M. Bonn, F. H. L. Koppens, E. Hendry, M. Polini, K. J. Tielrooij. The ultrafast dynamics and conductivity of photoexcited graphene at different Fermi energies. Sci. Adv., 4, eaar5313(2018).

    [28] D. Sun, C. Divin, M. Mihnev, T. Winzer, E. Malic, A. Knorr, J. E. Sipe, C. Berger, W. A. de Heer, P. N. First, T. B. Norris. Current relaxation due to hot carrier scattering in graphene. New J. Phys., 14, 105012(2012).

    [29] M. Mittendorff, T. Winzer, E. Malic, A. Knorr, C. Berger, W. A. de Heer, H. Schneider, M. Helm, S. Winnerl. Anisotropy of excitation and relaxation of photogenerated charge carriers in graphene. Nano Lett., 14, 1504(2014).

    [30] B. L. Huang, C. P. Chuu, M. F. Lin. Asymmetry-enriched electronic and optical properties of bilayer graphene. Sci. Rep., 9, 859(2019).

    [31] Z. Q. Li, E. A. Henriksen, Z. Jiang, Z. Hao, M. C. Martin, P. Kim, H. L. Stormer, D. N. Basov. Band structure asymmetry of bilayer graphene revealed by infrared spectroscopy. Phys. Rev. Lett., 102, 037403(2009).

    [32] W. Aroua. Metallic nanoparticles/graphene-molecules hybrid system-based active biosensor. Chin. Opt. Lett., 19, 123603(2021).

    Tingyuan Jia, Shaoming Xie, Zeyu Zhang, Qinxue Yin, Chunwei Wang, Chenjing Quan, Xiao Xing, Juan Du, Yuxin Leng. Role of the interlayer interactions in ultrafast terahertz thermal dynamics of bilayer graphene[J]. Chinese Optics Letters, 2022, 20(9): 093701
    Download Citation