Angle-tunable intersubband photoabsorption and enhanced photobleaching in twisted bilayer graphene
Eva A. A. Pogna, Xianchong Miao, Driele von Dreifus, Thonimar V. Alencar, Marcus V. O. Moutinho, Pedro Venezuela, Cristian Manzoni, Minbiao Ji, Giulio Cerullo, Ana Maria de Paula
Angle-tunable intersubband photoabsorption and enhanced photobleaching in twisted bilayer graphene
Eva A. A. Pogna , Xianchong Miao , Driele von Dreifus , Thonimar V. Alencar , Marcus V. O. Moutinho , Pedro Venezuela , Cristian Manzoni , Minbiao Ji , Giulio Cerullo , Ana Maria de Paula Istituto di Nanoscienze CNR-NANO, Laboratory NEST, Piazza San Silvestro 12, Pisa 56127, Italy Laboratory of Surface Physics and Department of Physics, Fudan University, Shanghai 200433, China Departamento de Física, Instituto de Ciências Exatas, Universidade Federal de Minas Gerais, Belo Horizonte-MG 31270-901, Brazil Departamento de Física, Instituto de Ciências Exatas e Biológicas, Universidade Federal de Ouro Preto, Ouro Preto-MG 35400-000, Brazil Núcleo Multidisciplinar de Pesquisas em Computação - NUMPEX-COMP, Campus Duque de Caxias, Universidade Federal do Rio de Janeiro, Duque de Caxias- RJ 25265-970, Brazil Instituto de Física, Universidade Federal Fluminense, UFF, Niterói-RJ 24210-346, Brazil IFN-CNR, Dipartimento di Fisica, Politecnico di Milano, Piazza L. da Vinci 32, Milano 20133, Italy
Abstract
Van der Waals heterostructures obtained by artificially stacking two-dimensional crystals represent the frontier of material engineering, demonstrating properties superior to those of the starting materials. Fine control of the interlayer twist angle has opened new possibilities for tailoring the optoelectronic properties of these heterostructures. Twisted bilayer graphene with a strong interlayer coupling is a prototype of twisted heterostructure inheriting the intriguing electronic properties of graphene. Understanding the effects of the twist angle on its out-of-equilibrium optical properties is crucial for devising optoelectronic applications. With this aim, we here combine excitation-resolved hot photoluminescence with femtosecond transient absorption microscopy. The hot charge carrier distribution induced by photo-excitation results in peaked absorption bleaching and photo-induced absorption bands, both with pronounced twist angle dependence. Theoretical simulations of the electronic band structure and of the joint density of states enable to assign these bands to the blocking of interband transitions at the van Hove singularities and to photo-activated intersubband transitions. The tens of picoseconds relaxation dynamics of the observed bands is attributed to the angle-dependence of electron and phonon heat capacities of twisted bilayer graphene.
Van der Waals heterostructures composed of stacked single-layer crystals represent a promising platform for novel electronic, photonic and spintronic devices [1-8]. Wide tunability of the optoelectronic properties is achieved by both selection of the stacked materials and control of the relative rotation between the crystal registries [9, 5,10-12], i.e. the twist angle. In the case of two layers of graphene, strong interlayer coupling drives the hybridization of the electronic states and can lead to the formation of moiré super-lattices for finite twist angle. The rotational misalignment of twisted bilayer graphene (tBLG) enables to explore also new properties as topological valley transport [13], superconductivity [14] at the magic twist angle θ ~ 1.1°, or the topologically protected 1D chiral states appearing at the moiré domain walls [15]. The Dirac cones associated to the two layers in tBLG are no longer aligned in the energy- momentum landscape and two saddle points appear at their low-energy intersections. The saddle points contribute to van Hove singularities (vHs) in the density of states, which occur at energies E vHs that depend on the twist angle [16]. The vHs are detected as a peak in the optical conductivity of tBLG [17-19] and they have been related to the increase of hot photoluminescence (PL) for resonant excitation [20]. These vHs appear at much lower energies than in single-layer graphene (SLG) [16,21-24] and allow angle controlled visible and near-infrared (NIR) absorption enhancement. This property has been exploited in a tBLG based photo-detector to selectively increase the photocurrent generation at the vHs [25]. The changes in crystal symmetry, electronic band structure and optical conductivity come along with modifications of the relaxation pathways available to charge carriers driven out of equilibrium. The carrier dynamics in SLG is well-known as it has been intensively studied by means of ultrafast transient absorption [26-31] (TA), time-resolved photocurrent [32,33] and photoemission spectroscopies [34-36]. Excitation with femtosecond light pulses produces out-of-equilibrium distributions of electrons and holes in the conduction and valence bands, which thermalize in ~ 50 fs via Coulomb two-body scattering to hot Fermi-Dirac (FD) distributions [30, 31]. The relaxation of the hot carriers proceeds through scattering with strongly coupled optical phonons (SCOP) till a common temperature is reached in about 200 fs. The cooling of the thermalized electron-phonon bath is finally driven by the lattice and achieved in less than 10 picoseconds via anharmonic decay of the hot phonons. In defected SLG, hot carriers can release their excess energy also by defect-mediated emission of acoustic phonons (supercollision mechanism), which occurs on a ps-time scale dependent on substrate and defect densities [31,33,37-43]. The hot carrier distribution inhibits interband absorption over a broad energy range due to Pauli blocking, so that in the TA spectra of SLG a decreased absorption, or photo-bleaching (PB), signal is detected. In Bernal stacked bilayer graphene, together with the PB band, Limmer et al. [44] observed a photoinduced absorption (PA) with a single decay time of 5 ps, about twice as long as the PB signal. The PA has been attributed to transitions between subbands formed within the conduction and the valence bands, named intersubband (ISB) transitions, which are absent in SLG since it lacks a subband structure. Also time- and angle-resolved photoemission spectroscopy (TARPES) has demonstrated long-living (up to 10 ps) population inversion in bilayer graphene [45,46] attributed to the presence of a small gap. Previous studies in tBLG by TA microscopy [47,48] and TARPES [49] have explored the twist-angle dependence and shown long decay times attributed to bound excitons and to states reached by two-photon absorption [48, 50]. Here we study the out-of-equilibrium optical properties of tBLG as a function of the twist angle by femtosecond TA microscopy with high sensitivity, which allows to explore low excitation regimes in which angle-dependent contribution from ISB transitions emerges. Specifically, we investigate heterostructures with the vHs located in the NIR and twist angle θ in the range of ~ 6°-8°. The tBLG samples are prepared by chemical vapor deposition as described in the Methods and they have tens of μm linear dimension. A lattice rotation in the real space by an angle θ corresponds to a rotation of θ in the phase-space of the hexagonal Brillouin zones (BZs), as sketched in Fig. 1(a-b). Accordingly, the Dirac points of the two layers (see Figure 1(c)), K A and K B , no longer coincide and are displaced by a momentum ΔK = 2K sin(θ/2) , where
K = 4π/3a and a ≈ st /2 nd layer [52]. In stack 1, the hexagons are rotated by an angle θ = 7°± 1°. The angle determination from the optical images is affected by a large uncertainty related to the definition of the single-crystals shape. Twist angle is usually measured by field high-resolution transmission electron microscopy [53,54], scanning tunneling microscopy [16] and lattice resolution atomic force microscopy [55]. Here we adopt an all-optical, far-field approach, based on measuring the hot-photoluminescence excitation (hot-PLE) spectrum under femtosecond laser excitation. As previously observed [20], due to the enhancement of pump absorption, a peak is expected in the bilayer emission when the laser is tuned at the vHs. In Fig. 1(e) the hot PL intensity of tBLG normalized to that of SLG ( I tBLG /I SLG ) is shown as function of laser excitation energy. The experimental data are fitted with a Gaussian function centered at E vHs = 1.3880 ± 0.007 eV. From the energy of the vHs, we extract a more accurate estimation of the twist angle considering the correlation θ [°] = A-√ (B-C E vHs [eV]), with A = 45.380, B = 2.030 × 10 , C = 4.186 × 10 extracted from analysis of the calculated joint density of states (JDOS) as a function of θ reported in the Fig. S-1 of the Supplementary Material (SM). We obtain the value θ = 7.31° ± 0.04° for stack 1, in good agreement with the evaluation from the optical image inspection. The hot PL map in Fig. 1(e) for laser excitation tuned at 1.4 eV, close to E vHs , reports indeed a ≈ -1 , the 2D peak at 2648 cm -1 and an additional peak at 1621 cm -1 are observed for tBLG, while for SLG these band are observed at 1588 cm -1 (G) and 2645 cm -1 (2D). In the bilayer, the G-band peak has a two times higher intensity, as shown in the inset of Fig. 1(f). The enhancement of the G-band and the extra peak, L a , testify the formation of a new band structure by interlayer coupling. It is thus an indication of the high crystalline quality of the samples and of a clean interface between the two monolayers forming the tBLG. The L a peak above the G-band position for the bilayer is the mode activated by intralayer electron-phonon scattering process from the unfolded longitudinal optical phonon (LO) [56] branch. The measured L a frequency of 1621 cm -1 is consistent with the theoretical predictions [56] for a tBLG with θ ~ 7°. The LO dispersion is expected to have a non-monotonic dependence on the twist angle, with a maximum frequency at θ ≈
10° but only slowly varying in the range of interest, thus making it unpractical to precisely define the twist angle only from the Raman data. We apply ultrafast TA microscopy to study the photophysics of tBLG samples as function of the twist angle. A NIR pump pulse with photon energy E pump is used to photo-excite the sample off-resonance with the vHs ( E pump = 0.8 eV stack 1 and E pump = 1.19 eV for stacks 2-3), while a second tunable pulse, the probe, monitors the relaxation of the induced hot carriers’ distributions as function of pump-probe time delay t . Pump and probe beams are collinearly coupled to a confocal microscope to achieve a ≈ 1 μm spatial resolution. The differential transmission ( ΔT/T ) is measured as a function of the probe photon energy E probe , using two distinct high sensitivity TA microscopes described in the Methods, guaranteeing a spectral coverage from 0.85 eV to 1.65 eV. Selected
ΔT/T maps at fixed E probe and time delay t = 150 fs are reported in Figure 2, for three tBLG stacks with different θ . As expected, the SLG regions exhibit a PB signal at all the investigated E probe , detected as a positive ΔT/T signal. The PB signal only shows a smooth monotone increase towards lower energy, following the momentum-spreading of the hot-carriers FD distribution. This behavior contrasts to that of tBLG, where a complex dependence on E probe emerges showing also a PA signal. Specifically, under low photo-excitation fluences, we identify two characteristic components in the transient response: a clear peak in the PB and the appearance of PA at lower energy. The PB peak corresponds to a factor 4-6 enhancement with respect to the SLG signal and it is centered in proximity to E vHs (see first column of maps of Figure 2). The PA is detected for E probe < E vHs . At lower probe energies, about 300 meV below E vHs , the ΔT/T signal of tBLG becomes barely distinguishable from that of SLG. While a similar behavior is observed in all the investigated samples, the energy position of PB and PA bands changes with θ . The larger is the twist angle, the higher is E probe at which similar tBLG/SLG contrast is observed. In stacks 2 and 3, bilayer regions with different twist angle are present, so that the signal intensity evolves differently according to their orientation. For all samples, the precise value of the twist angle is determined by recording the hot-PLE spectrum. The experimental ΔT/T spectra of monolayer and bilayer regions at t = 150 fs for several tBLG samples are reported in Figure 3(a), together with the PLE spectra used to extract E vHs and θ (blue symbols). As θ decreases, the overall TA spectrum experiences a red-shift. The PLE peak almost overlaps with the PB peak, with the larger mismatch appearing in the stack with the most intense PA band. The signal near the vHs, indeed, results from the interplay of PA and PB contributions whose relative intensity varies with the twist angle and excitation fluence. Remarkably, the PA emerges at low photo-excited carrier densities of the order of n = 1 × 10 cm -2 , as estimated from an incident fluence of 50 μJ cm -2 considering 4% absorption. At an incident fluence of 350 μJ cm -2 , corresponding to n = 7 × 10 cm -2 , the PA band disappears and a PB is observed over all the spectral range (see Fig. S-2 and S-3 of SM). The lack of previous evidences [47,48] of contribution from the ISB to the TA of tBLG could be traced back to the higher absorbed photon flux attained by pumping resonantly to the vHs. This dominance of the PB over the PA for high laser fluence was also observed for the intra- and inter-band transitions in suspended single layer graphene [31]. The dependence of the experimental TA spectrum on crystal orientation is traced back to the modifications of the phonon modes, electronic structure and of the JDOS, which we theoretically compute as a function of the twist angle and the hot-electron temperature T e (see details in Methods section) [57]. Figure 3(b) shows the calculated JDOS considering the interband and intersubband transitions between electronic states for the initial hot carrier temperature of 1500 K and same twist angles as in panel (a). The electronic states populated by the pump pulse allow for the ISB transitions. The temperature of the quasi-equilibrium hot electron distribution is obtained considering a three temperature model [58] (see Fig. S-4 and S-5 and details in the SM). The black (red) curves in Fig. 3(b) show the JDOS for the interlayer ISB (interband-IB) transitions as a function of the twist angle θ . The ISB curves are calculated for the electron and hole transitions as indicated in Fig. 1(c) by the solid and dashed red arrows, respectively. The higher energy peak in the JDOS for ISB transitions corresponds to holes’ excitations, while the lower-energy peak to the electrons. The difference in energy for the electron and hole ISB transitions is due to the small difference of their curvatures. Note that for all the twist angles, the electron ISB peak appears about 50 meV below the IB peak due to the vHs transitions, while the hole ISB peak is almost overlapped with the vHs transitions (see Fig. S-1 of SM). Figure 3(c) shows a comparison of both IB and ISB JDOS for carrier temperature of 1500 K and room temperature (300 K); the blue arrows indicate that the ISB JDOS intensity increases with the temperature while the IB JDOS intensity decreases. The intensity of the JDOS is related to the dynamical conductivity (or absorption). In this way, the difference of the JDOS at these two temperatures is proportional to the ΔT/T spectrum and qualitatively explains the opposite sign of the contributions from ISB and IB transitions to the TA signal. The increase of the JDOS associated to ISB transitions due to photo-excitation corresponds to the PA signal observed at energies just below E vHs in the experimental data of Figure 3(a). The ISB peak for holes is very close to the E vHs transition that dominates the ΔT/T signal at these energies. The strong angle dependence of the PA energy position results in ~ 250 meV tuning with only a ~ 1.26° twist angle change. Model predictions of the PA energy position based on the JDOS calculation, with the folding method, agree very well with the experimental data. In the SM, the expected values for a large range of twist angles (from 0° to 30°) are reported in Fig. S-1 and appear always slightly red-shifted with respect to the vHs transition. Figure 4 reports the
ΔT/T dynamics of stack 1. A similar relaxation scenario is observed for all the stacks, once the shift of energy position of the transient spectrum is taken into account (see Fig. S-6 in the SM). The recovery dynamics of tBLG is consistently different from that of SLG. The dynamics of tBLG (full colored dots) is compared to that of SLG (open black symbols) for E probe corresponding to the main features of the optical response: the PB peak at the vHs, the PA corresponding to ISB transitions and the low-energy PB band. The ΔT/T dynamics in SLG is energy-independent and it is described by a bi-exponential decay convoluted to an instrumental response function (IRF) of ≈ 120 fs. The extracted time constants are τ = 120 ± 20 fs and τ = 1980 ± 60 fs. These values are characteristic of SLG samples with low defect density [43]. The transient transmissivity depends superlinearly on T e through the FD distribution function, so that the time evolution of T e and ΔT/T are not correlated in time. However, while the fast relaxation component is strongly influenced by the superlinearity, the slow relaxation component of transient transmissivity is generally assumed as a reliable estimation of cooling time via anharmonic coupling of the thermalized bath of hot-electron and SCOPs with the acoustic phonons [59]. The
ΔT/T(t) dynamics of tBLG, on the contrary, strongly depends on the transition energy even though, away from the PA and the vHs, it is very similar to that of SLG. At the vHs, the PB not only has a much higher initial intensity than for SLG, but it decays on a longer time scale. The dynamics, after deconvolution with the IRF, is best fit by a biexponential decay with τ = 140 ± 20 fs and τ = 2500 ± 100 fs. The PA band is best described by a three-exponential decay with time constants τ = 120 ± 30 fs, τ = 1960 ± 40 fs and τ = 70 ± 2 ps with the additional decay component responsible for the tail extended to 100 ps in Figure 4(b). Such long lifetimes can be explained by the peculiar temperature dependence of the JDOS that emerges at finite twist angles. We model the decay curves by the ratio of the calculated JDOS as a function of the electron temperature T e , see Fig. 4(c-d). We calculate the time evolution of T e from a three temperature model, considering the electron energy relaxation through strongly coupled and weakly coupled phonon modes [58]. The specific heats for the electrons, strongly coupled phonons and weakly coupled phonons are obtained considering their density of states from tight-binding [60,57] and force constant [61] methods parameterized by first principles calculations. For tBLG, it is important to include the additional low frequency and high density phonon modes at the Γ point, that are not present in SLG, as representing the out-of-plane breathing of adjacent layers interacting via weak van der Waals forces (details in the SM). The calculated ΔT/T(t) dynamics, both at the vHs transition and at the ISB transition, reproduce very well the experimental results. Indicating that the density of states enhancement at the vHs and the changes in the electron and phonon heat capacities can explain the observed long decay of the PA and PB bands. To conclude, we demonstrate that at low excited carrier density the transient optical response of tBLG shows two peaks, one due to an photoinduced absorption from intersubband transitions and another one due to photobleaching of the interband transitions in resonance with the vHs ( E vHs ). Similarly, to the vHs, the energy positions of these peaks can be controlled with the twist angle. At the vHs, we observe a carrier cooling time which is about 30% longer than for SLG, while the photoinduced intersubband signal shows a long-living tail up to 100 ps. A simple model, based on the JDOS enhancement and on the evolution of the electronic temperature as the hot electrons equilibrate with the tBLG phonon modes, explains well the experimental TA signal, both for the peaks position and for the time decays, without recurring to more complex multiparticle excitations as bound excitons or non-linear processes such as two-photon absorption. Methods
Sample preparation and characterization
The twisted bilayer graphene single crystals were prepared by ambient pressure chemical vapor deposition (AP-CVD) synthesis on polycrystalline Cu foils. The nucleation conditions on Cu foils are controlled to allow both constituent layers to form as hexagonal single crystals, and each tBLG grain has a specific twist angle. The samples were transferred to a transparent amorphous quartz substrate for the optical images and spectroscopy measurements. More details about the sample preparation are presented in Ref. [62]. To select the samples for TA studies we first measured the sample twist angle by an analysis of the flake shapes on optical microscopy white light images. The optical images were obtained using an Olympus microscope (DP72) in reflection mode with a 100X objective lens. To obtain the twist angle we drew a hexagon for each layer and measured the angle between the edges of the two hexagons [52], as can be seen for the sample in Figure 1(d). We also characterized the sample twist angle by Raman spectroscopy using an alpha300R WITec confocal Raman spectrometer at excitation energy of 1.96 eV (633 nm). The Raman spectra were measured in backscattering geometry on a triple monochromator spectrometer equipped with a liquid-N2 CCD. All measurements were performed at room temperature.
For the PLE measurements a PL image is acquired for each excitation laser energy, that is tuned to cover the range of the vHs transition for each sample. The PL intensity is averaged at the bilayer region and divided by the PL at the monolayer region to obtain each point in the PLE spectra. The PL signal is measured in back scattering, separated from the excitation laser by a dichroic mirror (Semrock FF665-Di02). A band pass filter (75 nm centered at 610 nm) is used in front of the detector to further block any scattered laser light (Chroma 610/75M).
High-sensitivity TA microscopy
The TA data here reported have been acquired using two distinct microscopes working in complementary probe energy ranges, enabling an overall spectral coverage of 1 eV. The microscope used for the probe energy range of 0.8-1.4 eV is based on a Er-doped fiber laser (Toptica-Femtofiber pro) generating 300 mW, 150 fs pulses centered at 1550 nm with 40 MHz repetition rate. A portion of the output of the laser is used as the pump pulse tuned to 0.8 eV and is modulated with an acousto-optic modulator operating at 1 MHz. The probe pulse extending from 0.73 to 1.45 eV is obtained by super continuum (SC) generation from the laser fundamental focused in a Er-doped highly-nonlinear fiber. The low-energy branch of the SC and the fundamental are filtered out with a short pass filter cutting at 1500 nm. The pump and the probe beams are collinearly focused on the sample with an objective (Olympus-LCPLN-IR with magnification 100 X and NA = 0.85) to a spot size of about 1 μm. The probe transmitted by the sample is collected with an 8 mm focal length achromatic doublet, spectrally filtered with a monochromator with 5 nm resolution and detected by an InGaAs amplified photodiode with 4 MHz bandwidth. The TA signal is measured with a lock-in amplifier with 300 ms effective time constant, resulting in a
ΔT/T sensitivity of 10 -7 . TA dynamics is monitored by changing the pump-probe delay with an optical delay line, while the TA maps (images) at fixed probe photon energy and time delay are acquired by moving the sample with a motorized three-axis piezo-stack linear stage (Newport NPXYZ100). The microscope used for the probe energy range of 1.25-1.65 eV is based on a commercial femtosecond optical parametric oscillator (OPO, Insight DS+, Newport, CA) used as the laser source, with the fundamental 1040 nm (1.1923 eV) beam as the pump and the tunable output 750-990 nm (1.252-1.653 eV) as the probe. The pump beam is modulated at 20 MHz by an electro-optical modulator, collinearly combined with the probe beam and delivered into a laser-scanning microscope (FV1200, Olympus). The laser beams are focused onto the sample with an air objective (UPLSAPO 20 X, NA = 0.75, Olympus) to about 1.5 μm diameter, raster scanned by a pair of galvo mirrors, transmitted through the sample, passed through an optical filter to block the pump beam, and then directed onto a large- area silicon photodiode. Pump fluence is kept at approximately 50 μJ cm -2 (2 x 10 incident photons/cm ). The pump and probe beams are combined with a dichroic mirror (DMSP1000, Thorlabs) and aligned collinearly. The probe beam is optically filtered by two short-pass filters (FF011010/SP-25, Semrock), collected by a photodiode, and demodulated with a commercial lock-in amplifier (HF2LI, Zurich Instruments) to extract the TA signal ( ΔA ). A pixel dw ell time of 2 μs and image size of 320×320 pixels are used in the experiments. The cross-correlation between pump and probe pulses at the sample is measured to be ~ 0.3 ps, and the time delay interval between adjacent images is set to Δt = 133 fs. A single t scan of 20 ps takes about 1.5 min, and a complete spectral scan takes about 30 min. Theoretical model
The electronic structures were obtained by folding the SLG calculation results. For SLG, we follow the procedure given in Ref. [60, 57], where the electronic structure calculations are based on a fifth neighbors tight-binding approach, with one orthonormalized p z orbital per site, in which the hopping parameters are fitted to reproduce density functional theory (DFT) calculations with many-body corrections. The electronic transitions between bands in tBLG are captured considering two SLG BZs twisted by a certain angle θ . In this way, we considered only four electronic bands for the tBLG, one valence and one conduction band for each layer. For low temperatures, only the electronic IB transitions are possible. The magenta arrow in Fig. 1(c) shows the IB that coincides with the transition between the vHs, that gives rise to a peak in the JDOS at an energy E vHs which is a function of θ . When the temperature increases, ISB transitions between valence, or between conduction bands, are also possible (red arrows in Fig. 1(c)). The allowed IB and ISB transitions do not depend on the size of the minigap, which appears due to the hybridization of the eigenstates near to the band intersections. Our model cannot reproduce the interaction between layers and the atomic relaxation is not taken into account. However, Havener et al. [19] have shown that, although the interlayer interactions perturb the bands from each layer producing minigaps, which splits the valence and the conduction bands in tBLG, the difference in energy between the minigaps is generally too small to be seen directly in the PL and absorption spectra, contributing only to the peak broadening. Also, since the allowed electronic transitions do not depend on the size of the minigap and the PL peak is usually symmetric, which means that both transitions from minigaps to the conduction and the valence bands contribute equally, the energy position of the JDOS peak obtained by the folding approach is supposed to give a good estimation of the twisted angle. Therefore, the folding approach used in our simulations, which does not consider interaction between layers, is supposed to give good results [19, 56]. In our model, it is possible to analyze the JDOS for IB and ISB transitions separately, as a function of the electronic temperature. We compute the JDOS as 𝐽𝐽 ( 𝑇𝑇 𝑒𝑒 , 𝐸𝐸 ) = 𝛾𝛾𝜋𝜋𝑁𝑁 𝑘𝑘 � 𝑓𝑓�𝐸𝐸 𝑖𝑖 , 𝒌𝒌 , 𝑇𝑇 𝑒𝑒 � − 𝑓𝑓 ( 𝐸𝐸 𝑗𝑗 , 𝒌𝒌 , 𝑇𝑇 𝑒𝑒 ) �𝐸𝐸 𝑖𝑖 , 𝒌𝒌 − 𝐸𝐸 𝑗𝑗 , 𝒌𝒌 + 𝐸𝐸� + 𝛾𝛾 , 𝑗𝑗 , 𝒌𝒌 (1) where the sum is performed on a uniform grid of N k = 2400 × 2400 k-points over the two dimensional BZ and for the different pairs of bands i ; j while γ =1 × 10 -2 eV is a small energy broadening parameter related to the inverse of the electron-hole lifetime. The function 𝑓𝑓�𝐸𝐸 𝑖𝑖 , 𝒌𝒌 , 𝑇𝑇 𝑒𝑒 � = [1 + exp ( 𝐸𝐸 𝑖𝑖 , 𝒌𝒌 / 𝑇𝑇 𝑒𝑒 )] −1 is the Fermi function for a specific band energy i of the state k and temperature T e (the chemical potential μ is close to the Fermi level and both are set to zero). Using only IB transitions in the JDOS calculation we find the energy peaks corresponding to the E vHs for a specific twist angle θ . For the same values of E vHs obtained by the PLE peaks, we compute the normalized JDOS for IB and ISB transitions, as a function of the probe energy for T e = 1500 K, which is the estimated temperature for the hot electrons excited by the pump laser (see Fig. S-4 and S-5 in the SM). We see that the peaks corresponding to ISB transitions are always at lower energies than the IB transitions between the vHs, as shown in Fig. 3(c). An estimation of the ΔT/T dynamics can be obtained assuming that electron-radiation matrix element is constant for all transitions. In this way, the dynamical conductivity σ (or absorption α ) becomes proportional to the JDOS which we compute including both the IB and the ISB as function of the T e . Details of the T e time-evolution calculation are described in the SM. Acknowledgements
We acknowledge financial support from Graphene FET Flagship Core 3 Project, Grant No. 881603 and the Brazilian funding agencies: FAPERJ (grant number E-26/010.101126/2018), Fapemig, CNPq, Capes and INCT Carbon Nanomaterials. Prof. Po-Wen Chiu is gratefully acknowledged for preparing the tBLG samples.
Supplementary Material : twist angle estimation from E vHs , fluence dependence of TA spectra, electron temperature calculation, relaxation dynamics of tBLG stacks) is available. References [1] Geim, A. K.; Grigorieva, I. V. Van der Waals Heterostructures.
Nature , , 419. [2] Xia, F.; Wang, H.; Xiao, D.; Dubey, M.; Ramasubramaniam, A. Two-dimensional Material Nanophotonics. Nat. Photonics , , 899. [3] Ferrari, A. C. et al. Science and Technology Roadmap for Graphene, Related Two-dimensional Crystals, and Hybrid Systems. Nanoscale , , 4598. [4] Liu, Y.; Weiss, N. O.; Duan, X.; Cheng, H.-C.; Huang, Y.; Duan, X. Van der Waals Heterostructures and Devices. Nat. Rev. Mater. , , 16042. [5] Novoselov, K.; Mishchenko, A.; Carvalho, A.; Neto, A. C. 2D Materials and Van der Waals Heterostructures. Science , , 9439. [6] Zhong, D.; Seyler, K. L.; Linpeng, X.; Cheng, R.; Sivadas, N.; Huang, B.; Schmidgall, E.; Taniguchi, T.; Watanabe, K.; McGuire, M. A.; Yao, W.; Xiao, D.; Fu, K.-M. C.; Xu, X. Van der Waals Engineering of Ferromagnetic Semiconductor Heterostructures for Spin and Valleytronics. Sci. Adv . , , 1603113. [7] Unuchek, D.; Ciarrocchi, A.; Avsar, A.; Watanabe, K.; Taniguchi, T.; Kis, A. Roomtemperature Electrical Control of Exciton Flux in a Van der Waals Heterostructure. Nature , , 340. [8] Gibertini, M.; Koperski, M.; Morpurgo, A.; Novoselov, K. Magnetic 2D Materials and Heterostructures. Nat. Nanotech. , 14, 408. [9] Mishchenko, A. et al. Twist-controlled Resonant Tunnelling in Graphene/Boron Nitride/Graphene Heterostructures.
Nat. Nanotech. , , 808. [10] Ribeiro-Palau, R.; Zhang, C.; Watanabe, K.; Taniguchi, T.; Hone, J.; Dean, C. R. Twistable Electronics with Dynamically Rotatable Heterostructures. Science , , 690. [11] Kockum, A. F.; Miranowicz, A.; De Liberato, S.; Savasta, S.; Nori, F. Ultrastrong Coupling between Light and Matter. Nat. Rev. Phys . , , 19. [12] Dean, C.; Young, A.; Wang, L.; Meric, I.; Lee, G.-H.; Watanabe, K.; Taniguchi, T.; Shepard, K.; Kim, P.; Hone, J. Graphene Based Heterostructures. Solid State Commun. , , 1275. [13] Ju, L.; Shi, Z.; Nair, N.; Lv, Y.; Jin, C.; Velasco Jr, J.; Ojeda-Aristizabal, C.; Bechtel, H. A.; Martin, M. C.; Zettl, A.; Analytis, J.; Wang, F. Topological Valley Transport at Bilayer Graphene Domain Walls. Nature , , 650. [14] Cao, Y.; Fatemi, V.; Fang, S.;Watanabe, K.; Taniguchi, T.; Kaxiras, E.; Jarillo-Herrero, P. Unconventional Superconductivity in Magic-angle Graphene Superlattices. Nature , , 34. [15] Sunku, S. S.; Ni, G. X.; Jiang, B. Y.; Yoo, H.; Sternbach, A.; McLeod, A. S.; Stauber, T.; Xiong, L.; Taniguchi, T.; Watanabe, K.; Kim, P.; Fogler, M. M.; Basov, D. N. Photonic Crystals for Nano-light in Moiré Graphene Superlattices. Science , , 1153. [16] Li, G.; Luican, A.; Dos Santos, J. L.; Neto, A. C.; Reina, A.; Kong, J.; Andrei, E. Observation of Van Hove Singularities in Twisted Graphene Layers. Nat. Phys . , , 109. [17] Zou, X.; Shang, J.; Leaw, J. N.; Luo, Z.; Luo, L.; La-O Vorakiat, C.; Cheng, L.; Cheong, S. A.; Su, H.; Zhu, J.-X.; Liu, Y.; Loh, K. P.; Castro Neto, A. H.; Yu, T.; Chia, E. E. M. Terahertz Conductivity of Twisted Bilayer Graphene. Phys. Rev. Lett. , , 067401. [18] Wang, Y.; Ni, Z.; Liu, L.; Liu, Y.; Cong, C.; Yu, T.; Wang, X.; Shen, D.; Shen, Z. Stacking-Dependent Optical Conductivity of Bilayer Graphene. ACS Nano , , 4074. [19] Havener, R. W.; Liang, Y.; Brown, L.; Yang, L.; Park, J. Van Hove Singularities and Excitonic Effects in the Optical Conductivity of Twisted Bilayer Graphene. Nano Lett. , , 3353. [20] Alencar, T. V.; von Dreifus, D.; Moreira, M. G. C.; Eliel, G. S. N.; Yeh, C.-H.; Chiu, P.-W.; Pimenta, M. A.; Malard, L. M.; de Paula, A. M. Twisted Bilayer Graphene Photoluminescence Emission Peaks at Van Hove Singularities. J. Phys.
Condens. Mat. , , 175302. [21] Lopes dos Santos, J. M. B.; Peres, N. M. R.; Castro Neto, A. H. Graphene Bilayer with a Twist: Electronic Structure. Phys. Rev. Let t. , , 256802. [22] Tabert, C. J.; Nicol, E. J. Optical Conductivity of Twisted Bilayer Graphene. Phys. Rev.B , , 121402. [23] Moon, P.; Koshino, M. Optical Absorption in Twisted Bilayer Graphene. Phys. Rev. B , , 205404. [24] Rozhkov, A.; Sboychakov, A.; Rakhmanov, A.; Nori, F. Electronic Properties of Graphenebased Bilayer Systems. Phys. Rep. , , 1. [25] Yin, J.; Wang, H.; Peng, H.; Tan, Z.; Liao, L.; Lin, L.; Sun, X.; Koh, A. L.; Chen, Y.; Peng, H.; Liu, Z. Selectively Enhanced Photocurrent Generation in Twisted Bilayer Graphene with Van Hove Singularity. Nat. Commun. , , 10699. [26] Sun, D.; Wu, Z.-K.; Divin, C.; Li, X.; Berger, C.; de Heer, W. A.; First, P. N.; Norris, T. B. Ultrafast Relaxation of Excited Dirac Fermions in Epitaxial Graphene Using Optical Differential Transmission Spectroscopy. Phys. Rev. Lett. , , 157402. [27] Newson, R. W.; Dean, J.; Schmidt, B.; van Driel, H. M. Ultrafast Carrier Kinetics in Exfoliated Graphene and Thin Graphite Films. Opt. Express , 17, 2326. [28] Strait, J. H.;Wang, H.; Shivaraman, S.; Shields, V.; Spencer, M.; Rana, F. Very Slow Cooling Dynamics of Photoexcited Carriers in Graphene Observed by Optical-pump Terahertzprobe Spectroscopy.
Nano Lett. , , 4902. [29] Breusing, M.; Kuehn, S.; Winzer, T.; Malić, E.; Milde, F.; Severin, N.; Rabe, J. P.; Ropers, C.; Knorr, A.; Elsaesser, T.
Ultrafast Nonequilibrium Carrier Dynamics in a Single Graphene Layer.
Phys. Rev. B , , 153410. [30] Brida, D.; Tomadin, A.; Manzoni, C.; Kim, Y. J.; Lombardo, A.; Milana, S.; Nair, R. R.; Novoselov, K. S.; Ferrari, A. C.; Cerullo, G.; Polini, M. Ultrafast Collinear Scattering and Carrier Multiplication in Graphene. Nat. Commun. , , 1987. [31] Malard, L. M.; Mak, K. F.; Neto, A. H. C.; Peres, N. M. R.; Heinz, T. F. Observation of Intra- and Inter-band Transitions in the Transient Optical Response of Graphene. New J. Phys. , , 015009. [32] Song, J. C.; Rudner, M. S.; Marcus, C. M.; Levitov, L. S. Hot Carrier Transport and Photocurrent Response in Graphene. Nano Lett. , , 4688. [33] Graham, M. W.; Shi, S.-F.; Wang, Z.; Ralph, D. C.; Park, J.; McEuen, P. L. Transient Absorption and Photocurrent Microscopy Show that Hot Electron Supercollisions Describe the Rate-Limiting Relaxation Step in Graphene. Nano Lett. , , 5497. [34] Gilbertson, S.; Dakovski, G. L.; Durakiewicz, T.; Zhu, J.-X.; Dani, K. M.; Mohite, A. D.; Dattelbaum, A.; Rodriguez, G. Tracing Ultrafast Separation and Coalescence of Carrier Distributions in Graphene with Time-resolved Photoemission. J. Phys. Chem. Lett. , , 64. [35] Gierz, I.; Petersen, J. C.; Mitrano, M.; Cacho, C.; Turcu, I. E.; Springate, E.; Stöhr, A.; Köhler, A.; Starke, U.; Cavalleri, A. Snapshots of Non-equilibrium Dirac Carrier Distributions in Graphene . Nat. Mater. , , 1119. [36] Ulstrup, S.; Johannsen, J. C.; Crepaldi, A.; Cilento, F.; Zacchigna, M.; Cacho, C.; Chapman, R. T.; Springate, E.; Fromm, F.; Raidel, C. Ultrafast Electron Dynamics in Epitaxial Graphene Investigated with Time-and Angle-resolved Photoemission Spectroscopy. J. Condens. Matter Phys , , 164206. [37] Song, J. C. W.; Reizer, M. Y.; Levitov, L. S. Disorder-Assisted Electron-Phonon Scattering and Cooling Pathways in Graphene. Phys. Rev. Lett. , , 106602. [38] Song, J. C. W.; Tielrooij, K. J.; Koppens, F. H. L.; Levitov, L. S. Photoexcited Carrier Dynamics and Impact-excitation Cascade in Graphene. Phys. Rev. B , , 155429. [39] Betz, A.; Jhang, S. H.; Pallecchi, E.; Ferreira, R.; Fève, G.; Berroir, J.-M.; Plaçais, B. Supercollision Cooling in Undoped Graphene. Nat. Phys. , , 109. [40] Johannsen, J. C.; Ulstrup, S.; Cilento, F.; Crepaldi, A.; Zacchigna, M.; Cacho, C.; Turcu, I. C. E.; Springate, E.; Fromm, F.; Raidel, C.; Seyller, T.; Parmigiani, F.; Grioni, M.; Hofmann, P. Direct View of Hot Carrier Dynamics in Graphene. Phys. Rev. Lett. , 111, 027403. [41] Graham, M. W.; Shi, S.-F.; Ralph, D. C.; Park, J.; McEuen, P. L. Photocurrent Measurements of Supercollision Cooling in Graphene.
Nat. Phys. , , 103. [42] Tielrooij, K.; Song, J.; Jensen, S.; Centeno, A.; Pesquera, A.; Elorza, A. Z.; Bonn, M.; Levitov, L.; Koppens, F. Photoexcitation Cascade and Multiple Hot-carrier Generation in Graphene. Nat. Phys. , , 248. [43] Alencar, T. V.; Silva, M. G.; Malard, L. M.; de Paula, A. M. Defect-Induced Supercollision Cooling of Photoexcited Carriers in Graphene. Nano Lett. , , 5621. [44] Limmer, T.; Feldmann, J.; Da Como, E. Carrier Lifetime in Exfoliated Few-Layer Graphene Determined from Intersubband Optical Transitions. Phys. Rev. Lett. , , 217406. [45] Ulstrup, S. et al. Ultrafast Dynamics of Massive Dirac Fermions in Bilayer Graphene. Phys. Rev. Lett. , , 257401. [46] Gierz, I.; Mitrano, M.; Petersen, J. C.; Cacho, C.; Turcu, I. C. E.; Springate, E.; Sthr, A.; Khler, A.; Starke, U.; Cavalleri, A. Population Inversion in Monolayer and Bilayer Graphene. J. Phys. Condens. Mat. , , 164204. [47] Patel, H.; Havener, R. W.; Brown, L.; Liang, Y.; Yang, L.; Park, J.; Graham, M. W. Tunable Optical Excitations in Twisted Bilayer Graphene Form Strongly Bound Excitons. Nano Lett. , , 5932. [48] Patel, H.; Huang, L.; Kim, C.-J.; Park, J.; Graham, M. W. Stacking Angle-tunable Photoluminescence from Interlayer Exciton States in Twisted Bilayer Graphene. Nature Commun. , , 1445. [49] Fukumoto, K.; Boutchich, M.; Arezki, H.; Sakurai, K.; Felice, D. D.; Dappe, Y. J.; Onda, K.; Koshihara, S. Ultrafast Electron Dynamics in Twisted Graphene by Femtosecond Photoemission Electron Microscopy. Carbon , , 49. [50] Yang, H.; Feng, X.; Wang, Q.; Huang, H.; Chen, W.; Wee, A. T. S.; Ji, W. Giant Two-Photon Absorption in Bilayer Graphene. Nano Lett. , , 2622. [51] Chu, Z.-D.; He, W.-Y.; He, L. Coexistence of Van Hove Singularities and Superlattice Dirac Points in a Slightly Twisted Graphene Bilayer. Phys. Rev. B , , 155419. [52] H. Ribeiro, K. Sato, G. Eliel, E. de Souza, C.C. Lu, P.W. Chiu, R. Saito, M. Pimenta Carbon , , 138. [53] Havener, R. W.; Zhuang, H.; Brown, L.; Hennig, R. G.; Park, J. Angle-resolved Raman Imaging of Interlayer Rotations and Interactions in Twisted Bilayer Graphene. Nano Lett. , , 3162. [54] Latychevskaia, T.; Escher, C.; Fink, H.-W. Moiré Structures in Twisted Bilayer Graphene Studied by Transmission Electron Microscopy. Ultramicroscopy , , 46. [55] Schmidt, H.; Rode, J. C.; Smirnov, D.; Haug, R. J. Superlattice Structures in Twisted Bilayers of Folded Graphene. Nat. Commun. , , 5742. [56] Eliel, G. S. N.; Moutinho, M. V. O.; Gadelha, A. C.; Righi, A.; Campos, L. C.; Ribeiro, H. B.; Chiu, P.-W.; Watanabe, K.; Taniguchi, T.; Puech, P.; Paillet, M.; Michel, T.; Venezuela, P.; Pimenta, M. A. Intralayer and Interlayer Electron-phonon Interactions in Twisted Graphene Heterostructures. Nat. Commun. , , 1221. [57] Vela, A.; Moutinho, M.; Culchac, F.; Venezuela, P.; Capaz, R. B. Electronic Structure and Optical Properties of Twisted Multilayer Graphene. Phys. Rev. B , , 155135. [58] F. Caruso, F.; Novko, D; Draxl, C. Photoemission signatures of nonequilibrium carrier dynamics from first principles. Phys. Rev. B , , 035128. [59] Tielrooij, K. J. et al. Out-of-plane heat transfer in van der Waals stacks: electron hyperbolic phonon coupling . Nat. Nanotech. , , 41. [60] Venezuela, P.; Lazzeri, M; Mauri, F. Theory of double-resonant Raman spectra in graphene: Intensity and line shape of defect-induced and two-phonon bands. Phys. Rev. B , , 035433. [61] Cocemasov, A.I.; Nika, D.L.; Balandin, A. A. Phonons in twisted bilayer graphene. Phys. Rev. B , , 035428. [62] Lu, C.-C.; Lin, Y.-C.; Liu, Z.; Yeh, C.-H.; Suenaga, K.; Chiu, P.-W. Twisting Bilayer Graphene Superlattices. ACS Nano , 2587. Figure 1 (a) Sketch of the crystal structure of a tBLG with twist angle θ ; (b) the corresponding BZs of two rotated graphene layers as the red and blue hexagons; (c) energy structure indicating the optical transition at the vHs at photon energy E vHs (magenta line) and intersubband transitions from photo-excited FD of electrons (red solid lines) and holes (red dashed lines); ( d) bright field optical image of stack 1 with blue and red hexagons identifying a layer rotation angle of about 7°, scale bar 10 μm. (e) PL intensity map for laser excitation energy of 1.4 eV, in resonance with the E vHs identified by the peak in the PLE intensity of tBLG ( I tBLG ) divided by that of the corresponding SLG ( I SLG ); (f) Raman spectra of stack 1 at the bilayer (blue line) and monolayer (red line) regions normalized by the intensity of the G peak. The inset shows a zoom-in of the bilayer G and L a peaks divided by the monolayer G peak. Figure 2
Color map images of the transient transmissivity ΔT/T (10 -4 ) at time delay t = 0, for three tBLG stacks with twist angle increasing from the lower to the upper row of panels: θ = 7.31°, 7.62° and 8.57°. The probe energies, for the first column of maps are taken in resonance with the vHs at E probe =1.632, 1.512, 1.378 eV; second column at the maximum of the PA band at E probe = 1.512, 1.363, 1.240 eV; third column with low-energy PB band at E probe =1.319, 1.253, 0.886 eV, as indicated in the panels . Figure 3
Transient transmissivity spectra of SLG (red-dot and solid line) and tBLG (blackdot and solid line) at fixed time delay t = 150 fs for samples with different θ as indicated by the labels. The blue dots are the PLE intensity ratio of the tBLG to the SLG fitted with a Gaussian function (blue solid line) peaked at E vHs . (b) calculated JDOS for the ISB transition (black) and the E vHs transition (red) for a hot electron temperature of 1500 K. (c) A comparison of the IB and ISB JDOS for 1500 K and the room temperature 300 K. The difference of the JDOS at these two temperature is proportional to the ΔT/T spectrum.
Figure 4 (a) Transient transmissivity
ΔT/T dynamics of tBLG (full colored dots) and SLG (open black symbols) for the stack with θ = 7.31° at probe energy E probe = 1.378 eV corresponding to E vHs ; 1.240 eV corresponding to the ISB; 0.866 eV in the low-energy PB band E < E ISB . (b) Comparison of
ΔT/T(t) normalised to the maximum and with an extended time scale. (c) Theoretical estimated decay curves based on the electron temperature time decay. In the simulations we use E probe = 1.39, 1.24 and 0.70 eV to best reproduce the experimental decay behavior. (d) Normalized theoretical curves obtained in (c). Supplementary Material E vHs The electronic band structure and the corresponding density of states are evaluated for different twist angles. The energy position of the calculated vHs and ISB electronic transitions for holes and electrons are reported in Fig. S-1 as function of the twist angle. In analogy with Ref. 1 we find a convenient expression of θ fitting the theoretical data in the inset in the range 6.8° < θ < 8.6° with the function θ [°] = A-√ (B-C E vHs [eV]) , with A = 45.380, B = 2.030 × 10 , C = 4.186 × 10 . Figure S-1
Calculated energy of vHs for electrons 𝐸𝐸 𝑣𝑣𝑣𝑣𝑣𝑣− (red dots and line) and holes 𝐸𝐸 𝑣𝑣𝑣𝑣𝑣𝑣+ (magenta dots and lines) and of ISB transitions for electrons 𝐸𝐸 𝐼𝐼𝐼𝐼𝐼𝐼− (black dots and lines) and holes 𝐸𝐸 𝐼𝐼𝐼𝐼𝐼𝐼+ (blue dots and lines) as a function of the twist angle. The inset shows for the vHs energy in our experimental region of twist angles together with the best fit function (red solid line).
To calculate the uncertainty in the angle estimation we derive the expression: 𝛿𝛿𝛿𝛿 [°] =
𝐶𝐶2�𝐼𝐼−𝐶𝐶 𝐸𝐸 𝑣𝑣𝑣𝑣𝑣𝑣 [ 𝑒𝑒𝑒𝑒 ] 𝛿𝛿𝐸𝐸 , considering the 𝛿𝛿𝐸𝐸 as the PLE peak width. From the PLE results, the fitted Gaussian widths give 𝛿𝛿𝐸𝐸 ≤ E vHs = 1.388 ± 0.007 eV, for example, we find that θ = 7.312° ± 0.04°.
2. Fluence dependence of transient absorption spectra
The intersubband (ISB) transitions energetically overlap with interband optical transitions. The hot-electron distribution causes simultaneously the bleaching of interband absorption due to Pauli blocking and the activation of ISB absorption with contributions to the transient transmissivity of opposite signs. The interplay of these contributions depends on the fluence and, at high density of photo-excited carriers, we observe that the bleaching is dominating over the photo-induced absorption over the entire probe photon energy range. Accordingly, while the TA of SLG scales linearly with fluence, the spectra of tBLG change shape with fluence, as reported in Fig. S-2(a) for four different tBLG samples. The dynamics as function of pump-probe delay and probe wavelength is reported in Fig. S-2(b) at the two compared fluences for the stack with θ = 7.93°, taken as prototype of tBLG. While the photobleaching signal increases with excitation power, due to the increase in initial electron temperature and to the superlinear dependence of ΔT/T on T e , the instantaneous PA signal decreases. In Fig. S-3 we report the pump power dependence of relaxation dynamics of the stack with θ = 7.31° at probe energy E probe = 1.240 eV corresponding to the ISB, and E probe = 0.956 eV in the low-energy PB band. Figure S-2 (a)
ΔT/T spectra of tBLG at incident fluence of 350 μJ cm -2 (open symbols) and 50 μJ cm -2 (full symbols) for different twist angles θ ; (b) ΔT/T as a function of probe wavelength and time delay at low incident fluence 50 μJ cm -2 (upper panel) and hig h incident fluence 350 μ J cm -2 (lower panel) for the tBLG with θ = 7.93° . Figure S-3
Transient transmissivity dynamics at different excitation power at probe photon energy (a) E probe = 1.240 eV corresponding to the ISB and (b) E probe = 0.958 eV in the low-energy PB band. Normalised dynamics at E probe = 1.240 eV are reported in the inset of panel a.
3. Calculating the electron temperature
In order to obtain the electron and phonon temperatures as a function of time we adopt a three temperature model [2]. This model considers that after a transition to an out-of-equilibrium state the electrons rapidly exchange energy with phonons via electron-phonon interaction achieving a quasi-equilibrium Fermi distribution. Considering the occupation number of electrons and phonons as Fermi-Dirac and Bose-Einstein distributions respectively, we can associate an effective temperature for the electrons, T e , and for phonons, T ν which are separated into strongly coupled (ν = sc) and weakly coupled (ν = wc) phonons. Thus, a system of coupled equations for the three temperatures ( T e , T sc , T wc ) is used to simulate the dynamics of the electrons and phonons after the laser pulse excitation. For a Gaussian pulse I(t) with fluence F and pulse duration τ p , the temporal dependence of the effective temperatures is given by: 𝑑𝑑𝑇𝑇 𝑒𝑒 𝑑𝑑𝑑𝑑 = 𝐼𝐼 ( 𝑑𝑑 ) 𝛽𝛽𝑐𝑐 𝑒𝑒 − 𝑔𝑔 𝑣𝑣𝑠𝑠 𝑐𝑐 𝑒𝑒 ( 𝑇𝑇 𝑒𝑒 − 𝑇𝑇 𝑣𝑣𝑠𝑠 ) − 𝑔𝑔 𝑤𝑤𝑠𝑠 𝑐𝑐 𝑒𝑒 ( 𝑇𝑇 𝑒𝑒 − 𝑇𝑇 𝑤𝑤𝑠𝑠 ) 𝑑𝑑𝑇𝑇 𝑣𝑣𝑠𝑠 𝑑𝑑𝑑𝑑 = 𝑔𝑔 𝑣𝑣𝑠𝑠 𝑐𝑐 𝑣𝑣𝑠𝑠 ( 𝑇𝑇 𝑒𝑒 − 𝑇𝑇 𝑣𝑣𝑠𝑠 ) − 𝜏𝜏 𝑜𝑜 ( 𝑇𝑇 𝑣𝑣𝑠𝑠 − 𝑇𝑇 𝑤𝑤𝑠𝑠 ) 𝑑𝑑𝑇𝑇 𝑤𝑤𝑠𝑠 𝑑𝑑𝑑𝑑 = 𝑔𝑔 𝑤𝑤𝑠𝑠 𝑐𝑐 𝑤𝑤𝑠𝑠 ( 𝑇𝑇 𝑒𝑒 − 𝑇𝑇 𝑤𝑤𝑠𝑠 ) + 𝑐𝑐 𝑣𝑣𝑠𝑠 𝑐𝑐 𝑤𝑤𝑠𝑠 𝜏𝜏 𝑜𝑜 ( 𝑇𝑇 𝑣𝑣𝑠𝑠 − 𝑇𝑇 𝑤𝑤𝑠𝑠 ) (1) where β is a dimensionless parameter which determines the fraction of the pulse energy absorbed by the electrons and τ o = 2.5 ps is the optical phonon lifetime [3-4]. The respective specific heats for the electrons (c e ), strongly coupled phonons (c sc ) and weakly coupled phonons (c wc ) are obtained by the following integrals: 𝑐𝑐 𝑒𝑒 = � 𝜖𝜖𝑁𝑁 ( 𝜖𝜖 ) 𝜕𝜕𝑓𝑓 𝐷𝐷 ( 𝜖𝜖 , 𝑇𝑇 𝑒𝑒 ) 𝜕𝜕𝑇𝑇 𝑒𝑒∞−∞ 𝑑𝑑𝜖𝜖 𝑐𝑐 𝜈𝜈 = � 𝜔𝜔𝐹𝐹 𝜈𝜈 ( 𝜔𝜔 ) 𝜕𝜕𝑛𝑛 𝐼𝐼 ( 𝜔𝜔 , 𝑇𝑇 𝜈𝜈 ) 𝜕𝜕𝑇𝑇 𝜈𝜈∞0 𝑑𝑑𝜔𝜔 (2) (3) Here, N( 𝜖𝜖 ) and F ν (ω) are the density of states for electrons, strongly coupled (ν = sc) and weakly coupled (ν = wc) phonons obtained, for both SLG and tBLG, from tight-binding [5,6] and force constant [7] methods parameterized by first principles calculations, and 𝑓𝑓 𝐷𝐷 ( 𝑛𝑛 𝐼𝐼 ) is the Fermi-Dirac (Bose-Einstein) distribution. The difference between weakly and strongly coupled phonons in graphene is addressed considering that the optical A’ (TO) and E (LO) phonons in the vicinity of K and Γ points of the Brillouin zone are those with high electron-phonon coupling. In this way, in our model for SLG, only phonons with frequency close to (1350 ± 2) cm -1 and (1580 ± 2) cm -1 are considered in the calculation of the strongly coupled quantities, while the others modes are used for the weakly coupled ones. For tBLG, the additional low frequency and high density phonon modes close to ZA (88 ± 2) cm -1 at the Γ point are also considered as strongly coupled phonons. This mode is not present in the SLG because its polarization represents the out-of-plane breathing between two adjacent layers due to the weak van der Waals interaction. Figure S-4
Temperature dependence of the specific heat for electron and phonons and the electron-phonon coupling factors for SLG and tBLG. The indexes sc and wc stand for strongly coupled and weakly coupled phonons, respectively. The temperature dependent electron-phonon coupling factor can be calculated as in Ref. 8 by: 𝑔𝑔 𝜈𝜈 = − 𝜋𝜋ℏ𝑘𝑘 𝐼𝐼 𝜆𝜆 𝜈𝜈 〈𝜔𝜔 〉 𝜈𝜈 𝑁𝑁 ( 𝜖𝜖 𝐹𝐹 ) � 𝑁𝑁 ( 𝜖𝜖 ) 𝜕𝜕𝑓𝑓 𝐷𝐷 ( 𝜖𝜖 , 𝑇𝑇 𝑒𝑒 ) 𝜕𝜕𝜖𝜖 𝑑𝑑𝜖𝜖 ∞−∞ (4) Figure S-4 shows the temperature dependence of the specific heat for electron and phonons and the electron-phonon coupling factors for SLG and tBLG. The parameters used are 𝜆𝜆 𝑣𝑣𝑠𝑠 〈𝜔𝜔 〉 𝑣𝑣𝑠𝑠 = 566 meV and 𝜆𝜆 𝑤𝑤𝑠𝑠 〈𝜔𝜔 〉 𝑤𝑤𝑠𝑠 = 69 meV , that are the same values used for SLG [2]. Figure S-5
The temporal evolution of the electron and phonon effective temperatures after excitation for a laser fluence of 50 μJ cm -2 . Figure S-5 shows the evolution of the effective temperatures for an incident fluence F = 50 μJ cm -2 and a pulse duration of τ p = 200 fs. To consider the difference in optical absorption for SLG and tBLG, we used β = 50 for the former and β = 25 for the latter. These values were chosen to compare with the experimental results.
4. Relaxation dynamics of tBLG stacks
Similar relaxation dynamics is measured in all the probed tBLG stacks presented in Figure 2 and 3 as long as the probe energy is tuned to match the energy of the three different bands we have discussed. The dynamics of tBLG are reported in Fig. S-6 as colored dots compared to that of the monolayer, displayed as black symbols. The probe energies are specified in the graphs, and correspond to the vHs (pink dots), the peak of the PA band (red symbols) and the low energy tail.
Figure S-6
Transient transmissivity
ΔT/T dynamics of tBLG (full colored dots) and SLG (open black symbols) for the stacks with θ = 8.31°, 8.01°, 7.69°, 7.45°, at the probe energies indicated in the graphs . References [1] Havener, R. W.; Liang, Y.; Brown, L.; Yang, L.; Park, J. Van Hove Singularities and Excitonic Effects in the Optical Conductivity of Twisted Bilayer Graphene.
Nano Lett. , , 3353. [2] F. Caruso, F.; Novko, D; Draxl, C. Photoemission signatures of nonequilibrium carrier dynamics from first principles. Phys. Rev. B , , 035128. [3] Wang, H. et al.Ultrafast relaxation dynamics of hot optical phonons in graphene. App. Phys. Lett. , , 081917. [4] Rana, F.; George, P. A.; Strait, J. H.; Dawlaty, J.; Shivaraman, S.; Chandrashekhar, M.; Spencer, M. G. et al. Carrier recombination and generation rates for intravalley and intervalley phonon scattering in graphene. Phys. Rev. B , , 115447. [5] Venezuela, P.; Lazzeri, M; Mauri, F. Theory of double-resonant Raman spectra in graphene: Intensity and line shape of defect-induced and two-phonon bands. Phys. Rev.B , , 035433. [6] Vela, A.; Moutinho, M.; Culchac, F.; Venezuela, P.; Capaz, R. B. Electronic Structure and Optical Properties of Twisted Multilayer Graphene. Phys. Rev. B , , 155135. [7] Cocemasov, A.I.; Nika, D.L.; Balandin,A. A. Phonons in twisted bilayer graphene. Phys. Rev. B , , 035428. [8] Wang, X. Y.; Riffe, D. M.; Lee, Y.-S.; Downer, M. C. Time-resolved electron-temperature measurement in a highly excited gold target using femtosecond thermionic emission. Phys. Rev. B ,50