Abstract
Atomically thin layered van der Waals heterostructures feature exotic and emergent optoelectronic properties. With growing interest in these novel quantum materials, the microscopic understanding of fundamental interfacial coupling mechanisms is of capital importance. Here, using multidimensional photoemission spectroscopy, we provide a layer- and momentum-resolved view on ultrafast interlayer electron and energy transfer in a monolayer-WSe2/graphene heterostructure. Depending on the nature of the optically prepared state, we find the different dominating transfer mechanisms: while electron injection from graphene to WSe2 is observed after photoexcitation of quasi-free hot carriers in the graphene layer, we establish an interfacial Meitner-Auger energy transfer process following the excitation of excitons in WSe2. By analysing the time-energy-momentum distributions of excited-state carriers with a rate-equation model, we distinguish these two types of interfacial dynamics and identify the ultrafast conversion of excitons in WSe2 to valence band transitions in graphene. Microscopic calculations find interfacial dipole-monopole coupling underlying the Meitner-Auger energy transfer to dominate over conventional Förster- and Dexter-type interactions, in agreement with the experimental observations. The energy transfer mechanism revealed here might enable new hot-carrier-based device concepts with van der Waals heterostructures.
Similar content being viewed by others
Introduction
The unique physical properties of atomically thin two-dimensional (2D) materials1,2 and constantly improving fabrication methods3,4 have led to a great interest in novel quantum materials based on van der Waals (vdW) heterostructures5. By stacking 2D materials, vdW heterostructures inherit the properties from individual constituents, and exotic physical phenomena may emerge due to the interfacial interaction5,6,7. An emblematic example is the emergence of superconductivity in twisted bilayer graphene when stacked at the so-called âmagic angleâ8. As another example, interlayer excitons, which are spatially separated yet Coulomb-bound electron-hole pairs in semiconducting transition metal dichalcogenide (TMDC) heterostructures allow exceptional control of optoelectronic properties9,10,11. Out of the vdW heterostructure library, a basic optoelectronic building block is a monolayer (ML) semiconducting TMDC in contact with graphene12. This hybrid structure represents a model system as it combines the strong light-matter coupling of TMDCs and the high mobility of massless Dirac carriers of graphene13. The gapless electronic structure of graphene allows for harvesting low-energy photons, extending the spectral range covered by conventional photodetectors to the near-infrared wavelength, which is highly beneficial for photovoltaic applications14.
Optoelectronic functionality in vdW heterostructures arises from careful design and control of optical transitions and interfacial transfer processes. Particularly, interfacial charge (ICT) and energy transfer (IET) are key processes that have triggered extensive experimental and theoretical efforts15,16,17,18,19,20 Using time-resolved optical spectroscopies, a strong reduction of the exciton lifetime21 and optically active charge-transfer excitations of TMDC/graphene heterostructures have been observed22,23, suggesting strong interlayer coupling and the underlying mechanisms have been discussed24,25,26. Moreover, the efficiency of IET processes like Förster-type coupling (based on electronic dipole-dipole interaction) has recently been investigated theoretically, pointing out the importance of energy-momentum conservation between participating quasiparticles15. These studies provide our current understanding of the mechanisms of interfacial interactions. However, it is still challenging to clearly distinguish and unravel the involved interlayer charge and energy transfer in vdW heterostructures based on optical spectroscopies, primarily transient absorption/reflection and terahertz spectroscopy, which are inherently sensitive to the selective spectral range and limited momentum accessible. Therefore, a momentum-resolved probe is required to monitor the dynamics directly and achieve a complete picture of interfacial charge and energy transfer processes, including those involving momentum-forbidden dark states.
Here, we use time- and angle-resolved photoemission spectroscopy (trARPES) to investigate ultrafast interlayer carrier interactions in an epitaxially grown ML-WSe2/graphene heterostructure. Our trARPES setup combines a high-repetition-rate (500âkHz) femtosecond extreme ultraviolet (XUV) source27 coupled to a time-of-flight momentum microscope28 (see Methods). It allows the measurement of the four-dimensional (4D) photoemission intensity I(Ekin, kx, ky, Ît), where Ekin is the outgoing photoelectron kinetic energy, kx,ky are the in-plane momenta and Ît is the pump-probe delay, as shown in Fig. 1a, b. The probe photon energy of 21.7âeV allows accessing the entire Brillouin zone of the heterostructure and the variable pump wavelength allows us to photoexcite the heterostructure in a state-resolved manner. In the following, we present a time-, energy-, and momentum-resolved study on the excited-state dynamics in the heterostructure with two different pump photon energies: below the optical bandgap of WSe2 (1.2âeV) and in resonance with its first excitonic transition (1.55âeV).
Results
Interlayer quasi-free carrier transfer
First, we photoexcite the heterostructure with the pump photon energy centered at âÏpump=1.2âeV (pump pulse duration 200âfs FWHM), well below the optical bandgap of WSe229. The NIR-pump/XUV-probe experiments were performed with a pump fluence of Fâ=â5.3 mJ/cm2 and at room temperature. Figure 2aâd shows energy-resolved photoemission signals along the \({{{{{{{\rm{{K}}}}}}}^{{\prime} }}}-{{{{{{{\rm{K}}}}}}}}\) cut of the Brillouin zone, at selected time delays. The band mappings are contrast-enhanced using a multidimensional extension of the contrast limited adaptive histogram equalization (MCLAHE)30,31 for better visualization of the band structure. The momentum distributions above EF within the first 400âfs reveal that the excited states are localized in three different types of valleys: the Dirac cones of graphene at its K points (KGr) and the K and Q valleys of WSe2 (\({{{{{{{\rm{{K}}}}}}}_{WS{e}_{2}},{{{{{{{{\rm{Q}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}}}\)), as shown in Fig. 2e. The \({{{{{{{{\rm{Q}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\) valley localizes between the \({{{{{{{{\rm{K}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\) valley and the Î point. By performing energy-momentum integration in selected regions of interest (ROIs), we extracted excited-state dynamics within these three valleys (Fig. 2f). Upon arrival of the pump pulses, the excited-state population rapidly builds up at KGr (black curve) and decays with a time scale of ~200âfs. Strikingly, the conduction band minima (CBMs) at \({{{{{{{{\rm{K}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\) (red curve) and \({{{{{{{{\rm{Q}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\) valleys (green curve) are also being populated, however, with a delay of Îtâ=â51â±â9âfs (see SI) compared to the rise of hot-carrier population in graphene. Since below-bandgap pump photon energy does not allow the direct photoexcitation of WSe2, the delayed electron populations in the conduction bands arise through charge transfer from graphene to WSe2. Two/multiple photon excitation can safely be ruled out (details see SI). The excited-state population of the \({{{{{{{{\rm{Q}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\) valleys (Fig. 2f) could be raised via ICT from the graphene layer and the intervalley scattering from the \({{{{{{{{\rm{K}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\) valleys32,33.
These observations support the following picture of the underlying processes with a below-bandgap excitation: light is absorbed by graphene and populates unoccupied states at \({E}_{Gr}^{el}={E}_{D}+\hslash {\omega }_{{{{{\rm{pump}}}}}}/2\), leaving holes at \({E}_{Gr}^{h}={E}_{D}-\hslash {\omega }_{{{{{\rm{pump}}}}}}/2\) (Dirac energy EDâ>â0 for a p-doped system or EDâ<â0 for an n-doped system). The energy position of the Dirac point in our heterostructure is estimated to be ~ â0.1eV below the Fermi level, obtained from the conical crossing34,35 (see SI). The photoexcited carriers quickly reach a quasi-thermalized states in ~10âfs36 and could further increase their energy via intraband electron-electron scattering and interband Auger recombination in few tens of femtoseconds37,38. Once electrons gained a sufficient amount of energy to overcome the energy barrier, they scatter to WSe2 via a phonon-assisted tunneling process, filling the single-particle CBMs at \({{{{{{{{\rm{K}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\) and \({{{{{{{{\rm{Q}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\). This ICT mechanism is called interlayer hot-carrier injection, and is schematically illustrated in Fig. 2g. The excited electrons in WSe2 may subsequently scatter back to graphene and relax down towards the Fermi energy (EF). Based on the observed carrier dynamics, we performed microscopic calculations of the phonon-assisted interlayer tunneling process, allowing us to estimate the electronic wavefunction overlap between the involved conduction bands of WSe2 and graphene to be ~4% (see SI for details).
Interlayer energy transfer
Next, we select a pump photon energy of âÏpumpâ=â1.55âeV (pump pulse duration: 35âfs FWHM, pump fluence: Fâ=â1.7âmJ/cm2), near-resonant to the A-excitonic transition of WSe2. In this case, the pump photon energy allows both the WSe2 and the graphene layer to be simultaneously photoexcited. One striking observation is that the energy distribution of excited carriers at the \({{{{{{{{\rm{K}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\) valleys is centered at 0.63âeV (Fig. 3a), ~100âmeV lower than with below-bandgap excitation (Fig. 3b), as apparent from the energy distribution curves (EDCs) (first 100âfs). As discussed above, with 1.2âeV excitation, the \({{{{{{{{\rm{K}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\) valleys are filled with quasi-free electrons that have tunneled from the graphene layer. Therefore, this ~100âmeV energy difference is a direct photoemission signature of exciton formation, when near-resonantly pumping using 1.55âeV photons39: the bound electron-hole (el-h) pair reduces the quasi-free particle bandgap by the exciton binding energy. In addition to this excitonic feature, we also observe a transient shift of WSe2 valence bands. In Fig. 3d, EDCs at \({{{{{{{{\rm{K}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\) are shown at Îtâ=â0âfs (red) and Îtâ=âââ200âfs (black), in which the top two valence bands, VB1 and VB2, are fitted using Gaussian lineshape functions (see SI). The peak position of VB1 shifts towards the conduction band within the first 100âfs, transiently shrinking the electronic bandgap. This is due to the arrival of ICT-induced charge carriers from the graphene layer. With near-resonantly pumping the A-exciton, the occurrence of ICT and injection of quasi-free carriers from graphene to WSe2 is expected, similar to the case of below-bandgap excitation. This could lead to dynamical screening effect and the observed bandgap renormalization, as reported in highly-excited or doped ML TMDC materials40,41,42,43,44. As the magnitude of such a transient bandgap renormalization has been shown to scale with the excited charge carrier density42,45, we utilize the VB shift in the following as a measure of the ICT transferred carriers dynamics from graphene layer.
In addition to the excited-state dynamics in WSe2, important insight can be drawn from the energy-momentum distribution of hot carriers in graphene. As shown in the early-time 2D differential spectrum ÎI(E, k, Îtâ=â0âfs) (Fig. 3c), obtained by subtracting the spectrum at the negative time, hot carriers distribute in a broad energy range. The momentum-integrated spectrum along the linearly dispersing band in Fig. 3e clearly features the energy distribution of net electron gain (positive; red area) and loss (negative; blue area) following near-resonant photoexcitation. Remarkably, besides the modification of the distribution function near the Fermi level, we notice a strong negative peak at EâEFâ=ââ1.8âeV. As noted earlier, for direct photoexcitation in graphene the photoexcited carriers are expected to be spread ±0.77eV (âÏpump/2) around the Dirac point and quickly relax back to the Fermi level. Thus, this simple excitation mechanism cannot explain this peculiar feature in the valence band spectrum. The electron-electron scattering and Auger recombination could lead to a transient broadening of the momentum-space carrier distribution, but without any preferential energy localization38,46,47. hole transfer can also be ruled out, as the top valence band of WSe2 lies at EâEFâ=ââ1.0âeV. It would require a multi-phonon absorption to populate the hole-states localized deeply in the valence band, taking the typical phonon energy of ~0.17âeV in graphene48, a process of very low probability. However, the energy difference of deep-lying valence holes (EâEFâ=ââ1.8eV) and states near EF (EâEFâ=ââ0.2âeV) in graphene well matches the energy of the A-exciton in WSe2 (Eexâ~â1.6eV). Combined with the fast depletion of exciton population shown in Fig. 4a (black curve) extracted from the excited state of WSe2 (ROI1 in Fig. 3c), this brings about the following scenario for the excitation of these carriers: annihilation of excitons in WSe2 drives the intraband excitation of deep-lying valence electrons in graphene into empty hole states below the Dirac point. In more detail, this exciton energy transfer process, which we term Meitner-Auger energy transfer49,50, considers recombination of excitons in WSe2 with center-of-mass (COM) momentum Q and exciton energy Eex. The photoexcitation prepares the required hot hole vacancy below EF in graphene, thus enabling the intraband excitation. The photo-generated hole density plays an important role in the MA-type IET process (see the discussion of pump fluence dependence in SI). Besides the observation of the deep-lying hot holes, we also identify a substantial suppression of hole-like spectral weight (Meitner-Auger type IET-induced hot electrons) below the Fermi level with near-resonant excitation, supporting the occurrence of intraband transition in the graphene layer (details see SI, section Meitner-Auger type IET-induced hot electrons near the Fermi level). The momentum of the valence electron-hole pair kGr is determined by the Fermi velocity of the graphene bands and the transition energy EGr. This required momentum is provided by the optically pumped excitons which gain finite COM momenta during the population formation process via phonon-mediated dephasing and intravalley thermalization51,52,53,54 (see the discussion in SI). The highly efficient IET of the excitons and intraband electron-hole pairs is thus possible under the conservation of energy and momentum, i.e., Eexâ=âEGr and Qâ=âkGr. In a similar trARPES study of a ML WS2/graphene heterostructure, dominating interfacial charge transfer has been observed17. Compared with our study, the different charge transfer rates could be raised from the different band structure alignment near the interface and the density of defect sites26. While the additional exciton energy transfer was not excluded, its relative efficiency might be reduced due to the larger COM momentum required at the larger A-exciton energy of WS2 and the energy level alignment of these specific samples.
In order to gain information on the time scales of the energy and charge transfer processes, next we analyze the dynamics of excited-state populations extracted from the ROIs shown in Fig. 3c, including the excited-state carriers in WSe2 (ROI1), VB1 shifting (ROI2), hot electrons in graphene (ROI3) and IET-driven deep valence band holes (ROI4). The time trace of hot carriers in the CBM of WSe2 (black curve in Fig. 4a) contains two types of quasiparticles dynamics: the photo-generated excitons \({N}_{T}^{ex}\) and the ICT-induced quasi-free electrons \({N}_{T}^{el}\). The decay of excitons excite the valence band electrons in graphene via IET with a transfer time of ÏIET (Fig. 4f). On the other hand, the arrival of ICT-induced electrons transiently shifts the VBs of WSe2 (green curve in Fig. 4a) which therefore represents the dynamics of \({N}_{T}^{el}\) as discussed before. We assume VB1 and VB2 shift in the same way (fitting details see SI). The VB1 shifting shows a time delay of ~65âfs compared to the CB signal, evidencing the occurrence of interlayer hot electron injection after photoexcitation. The population of \({N}_{T}^{el}\) subsequently relaxes back to KGr, refilling the excited states of graphene (Fig. 4h). From the graphene side, the photoexcited hot electrons \({N}_{Gr}^{el}\) (red curve in Fig. 4b) could either scatter to conduction bands of WSe2 or relax by interband decay channels in graphene. Therefore, the relaxation of \({N}_{Gr}^{el}\) could be characterized with the charge transfer time of ÏICT and a decay time of \({\tau }_{Gr}^{el}\). The deep valence band holes \({N}_{Gr}^{h}\) (blue curve in Fig. 4b) are populated by exciton energy transfer on a time scale of ÏIET, which would relax back to the Fermi level with a lifetime of \({\tau }_{Gr}^{h}\).
The complete dynamics across the interface can be described with a set of coupled rate equations based on a multi-level scheme (details see SI). By numerically solving the rate-equation model, we disentangle the dynamics of IET and ICT. Our global fit describes the data well and yields the transfer times of ÏIETâ=â67â±â7âfs and ÏICTâ=â118â±â18fs. The lifetimes of electrons and IET-populated hot holes in graphene are simultaneously extracted as \({\tau }_{Gr}^{el}=84\pm 7\,{{{{{{{\rm{fs}}}}}}}}\) and \({\tau }_{Gr}^{h}=7\pm 4\,{{{{{{{\rm{fs}}}}}}}}\). Combining all our observations and analysis of the energy-momentum dynamics in WSe2 and graphene, we summarize the interfacial phenomena governing the non-equilibrium behavior of our heterostructure: first, the optical pump generates excitons in WSe2 and quasi-free carriers in graphene (Fig. 4e). Following photoexcitation, the exciton annihilation excites deep valence electrons in graphene via an IET process (Fig. 4f, g). Simultaneously, hot electrons in graphene are injected to the conduction bands of WSe2 via ICT which transiently shift the valence bands of WSe2 (Fig. 4h).
Discussion
To elucidate the interfacial coupling mechanism at play in our experiment, in particular the observed ultrafast energy transfer rate, we perform microscopic calculations of three types of IET mechanisms: Meitner-Auger, Förster, and Dexter energy transfer. The interlayer MA process is described by the dipole-monopole energy transfer from excitons to valence band excitation, schematically shown in Fig. 4e. The photoexcited hot holes in graphene quickly relax and distribute below EF near a transient chemical potential \({\mu }_{Gr}^{h*}\). This allows an MA-type transition from the deep valence band to the hot hole vacancy by absorbing the exciton energy. The microscopically calculated transfer rate is plotted as a function of Q in Fig. 4c with different transient chemical potentials for the hole distributions \({\mu }_{Gr}^{h*}\). When the hole vacancy is located around \({\mu }_{Gr}^{h*}=-0.3\,{{{{{{{\rm{eV}}}}}}}}\), the maximum transfer rate reaches ÎIETâ=â2.4âmeV, corresponding to a ÏIETâ=â270âfs transfer time. The MA-type IET process could describe the observed energy-momentum distribution of intraband transition of valence electrons in a reasonable quantitative agreement with the extracted transfer rate. As the transient chemical potential is subject to the doping level of graphene, we also calculate the MA-type IET with n-doped graphene by artificially increasing the Fermi energy. The IET-induced intraband transition in the valence bands is suppressed with decreased photo-generated hole vacancies. However, the intrinsically doped electrons above the Dirac point enable the MA-type IET in the conduction bands (details see SI).
Another IET mechanism is Förster energy transfer (Fig. 4f). The energy of the exciton excites an interband transition from valence bands to above Dirac point via the dipole-dipole coupling55. In contrast to the MA-type IET process, the interband excitation via Förster-type energy transfer populates the conduction bands of graphene above the Fermi level, independent of the photon-induced hot carriers distribution. The coupling strength is explicitly evaluated (for derivation, see SI) and determined by the momentum Q and interlayer distance d. The strong exciton oscillator strength and intrinsic in-plane exciton dipole moment in many 2D materials favor the Förster-type IET56. However, the calculated transfer rate is only 0.08âmeV (a transfer time of ~8.1âps), even assuming a tightly stacked heterostructure with interlayer distance of dâ=â0nm (Fig. 4d). Our calculations reveal that the IET process preferably excites an intraband rather than an interband transition. The experimentally observed energy-momentum distribution of excited-state hot holes supports this conclusion. To further distinguish the MA- and Förster-type IET, we calculate the transfer rates of these two mechanisms as a function of layer distance, and identify the distinct layer distance dependence (details see SI). In addition, we also performed calculations of Dexter-type IET (Fig. 4g), in which scenario the electron and hole components of excitons in WSe2 scatter to the graphene layer simultaneously. However, due to the small wavefunction overlap and the finite momentum distance between \({{{{{{{{\rm{K}}}}}}}}}_{{{{{{{{{\rm{WSe}}}}}}}}}_{{{{{{{{\rm{2}}}}}}}}}}\) and KGr, we found a very weak Dexter-type interlayer coupling strength, more than three orders of magnitude smaller compared to the other two mechanisms (see SI). Compared with Förster- and Dexter-type IET, the calculated transfer time of MA-type energy coupling is the closest to our experimental results. We can thus identify the MA-type conversion of excitons in WSe2 to intraband excitations in graphene as the dominant IET mechanism.
In this work, we provide a detailed microscopic picture of interfacial charge and energy transfer processes in photoexcited ML-WSe2/graphene heterostructures. Optical excitation of electrons in graphene leads to interlayer charge transfer of quasi-free electrons from the graphene layer to the K and Q valleys of the semiconductorâs conduction bands on a time scale of ~50âfs. In contrast, excitons in WSe2 decay through an interfacial Meitner-Auger energy transfer process with a time constant of ~70âfs. This previously unidentified process is governed by interlayer dipole-monopole interactions leading to annihilation of an exciton in WSe2 and non-vertical intraband excitations in graphene. The momentum of the electron-hole pair in graphene originates from the finite center of mass momentum of the hot excitons in WSe2. The interfacial Meitner-Auger mechanism is found to dominate the energy transfer process over established mechanisms like Förster- and Dexter-type transfer. This mechanism results in transient hole distributions as low as 2âeV below the Dirac points. These observations enrich the physical toolbox for designing van der Waals heterostructures and might be utilized in hot-carrier photovoltaic device concepts to harness the ultrafast and efficient carrier transfer processes at interfaces57.
Methods
Time- and angle-resolved photoemission spectroscopy
We used a 500âkHz tabletop femtosecond optical parametric chirped pulse amplification (OPCPA) laser system operating at a center wavelength of 800ânm and delivering average power up to 15âW. The high harmonic generation is produced in a vacuum chamber by tight focusing (10âμm) the second harmonic (400ânm) of the OPCPA fundamental on a thin and dense argon gas jet. We select the photons ~21.7âeV (110âmeV FWHM bandwidth) as the probe arm for trARPES experiment27. Concerning the pump arm, we used two different beams for this study. One pump beam is directly obtained from the OPCPA (800ânm, FWHMâ=â35âfs) and another one is the residual power of the compressed fiber amplifier (1030ânm, FWHMâ=â200âfs). The pump and probe beams are coupled into an ultra-high-vacuum (UHV) chamber and spatially overlapped at the sample position which is controlled by a six-axis manipulator (Carving, SPECS GmbH). The main UHV chamber is equipped with a unique combination of a hemispherical electron energy analyzer (PHOIBOS150, SPECS GmbH) and time-of-flight (ToF) momentum microscope (METIS1000, SPECS GmbH)28. On the one hand, the hemispherical analyzer, which can work in a multi-electrons per laser shot regime, provides high statistic energy/momentum cuts along a given momentum direction, as shown in Fig. 3. On the other hand, the momentum microscope allows for efficient, parallel, momentum-resolved detection of the full photoemission horizon from the surface as shown in Fig. 1b and Fig. 2aâe. All the experiments are performed at room temperature.
ML-WSe2/ML-graphene vdW heterostructure fabrication
Monolayer graphene on SiC (Si-terminated surface) was grown using the well-established recipe of sublimation growth at elevated temperatures in an argon atmosphere4. Note that, on SiC58, the graphene monolayer resides on top of a \((6\sqrt{3}\times 6\sqrt{3})\)R30â reconstructed carbon buffer layer that is covalently bound to the SiC substrate. WSe2 films were grown on the thus prepared MLG/SiC substrates via hybrid-pulsed-laser deposition (hPLD) in ultra-high vacuum35. Pure tungsten (99.99%) was ablated using a pulsed KrF excimer laser (248ânm) with a repetition rate of 10 Hz, while pure selenium (99.999%) was evaporated from a Knudsen cell at a flux rate of around 1.5âà /s as monitored by a quartz crystal microbalance. The deposition was carried out at 450â°C for 6 h, followed by two-step annealing at 640â°C and 400â°C for 1âh each.
DFT band structure calculations
We performed density functional theory (DFT) calculation of suspended ML WSe2 and graphene with the projector augmented wave code GPAW59 using GLLB-SC xc-functional, separately. The GLLB-SC is an orbital-dependent exact exchange-based functional including the spin-orbital coupling60. The relaxed lattice constant of WSe2 is aâ=â3.25âà . We sample the Brillouin Zone with a (15âÃâ15âÃâ1) k-point mesh, and set the cutoff energy for the plane-wave expansion at 600âeV. The bandgap is adjusted to fit our data. The calculated band structures of both materials are superimposed on each other and shown in Fig. 2a.
Data availability
The trARPES data generated in this study have been deposited in the Zenodo database under accession code https://doi.org/10.5281/zenodo.8210835.
Code availability
The source code for the trARPES data analysis is available on GitHub (https://github.com/mpes-kit). The data processing details are well described in ref. 61.
References
Novoselov, K. S., Mishchenko, A., Carvalho, A. & Neto, A. H. 2D materials and van der Waals heterostructures. Science 353, 6298 (2016).
Akinwande, D. et al. Graphene and two-dimensional materials for silicon technology. Nature 573, 507â518 (2019).
Liu, F. et al. Disassembling 2D van der Waals crystals into macroscopic monolayers and reassembling into artificial lattices. Science 367, 903â906 (2020).
Emtsev, K. V. et al. Towards wafer-size graphene layers by atmospheric pressure graphitization of silicon carbide. Nat. Mater. 8, 203â207 (2009).
Jin, C. et al. Ultrafast dynamics in van der Waals heterostructures. Nat. Nanotechnol. 13, 994â1003 (2018).
Lee, C.-H. et al. Atomically thin pân junctions with van der Waals heterointerfaces. Nat. Nanotechnol. 9, 676 (2014).
McGilly, L. J. et al. Visualization of moiré superlattices. Nat. Nanotechnol. 15, 580â584 (2020).
Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43â50 (2018).
Kunstmann, J. et al. Momentum-space indirect interlayer excitons in transition-metal dichalcogenide van der Waals heterostructures. Nat. Phys. 14, 801â805 (2018).
Rivera, P. et al. Observation of long-lived interlayer excitons in monolayer MoSe2âWSe2 heterostructures. Nat. Commun. 6, 6242 (2015).
Binder, J. et al. Upconverted electroluminescence via auger scattering of interlayer excitons in van der Waals heterostructures. Nat. Commun. 10, 2335 (2019).
Avsar, A. et al. Spinâorbit proximity effect in graphene. Nat. Commun. 5, 4875 (2014).
Garcia, J. H., Vila, M., Cummings, A. W. & Roche, S. Spin transport in graphene/transition metal dichalcogenide heterostructures. Chem. Soc. Rev. 47, 3359â3379 (2018).
Massicotte, M. et al. Photo-thermionic effect in vertical graphene heterostructures. Nat. Commun. 7, 12174 (2016).
Selig, M., Malic, E., Ahn, K. J., Koch, N. & Knorr, A. Theory of optically induced Förster coupling in van der Waals coupled heterostructures. Phys. Rev. B 99, 035420 (2019).
Froehlicher, G., Lorchat, E. & Berciaud, S. Charge versus energy transfer in atomically thin graphene-transition metal dichalcogenide van der Waals heterostructures. Phys. Rev. X 8, 011007 (2018).
Aeschlimann, S. et al. Direct evidence for efficient ultrafast charge separation in epitaxial WS2/graphene heterostructures. Sci. Adv. 6, 0761 (2020).
Hill, H. M. et al. Exciton broadening in WS2/graphene heterostructures. Phys. Rev. B 96, 205401 (2017).
Ferrante, C. et al. Picosecond energy transfer in a transition metal dichalcogenide-graphene heterostructure revealed by transient Raman spectroscopy. PNAS 119, 15 (2022).
Luo, D. et al. Twist-angle-dependent ultrafast charge transfer in MoS2-graphene van der Waals heterostructures. Nano Lett. 21, 8051â8057 (2021).
He, J. et al. Electron transfer and coupling in grapheneâtungsten disulfide van der Waals heterostructures. Nat. Commun. 5, 5622 (2014).
Krause, R. et al. Ultrafast charge separation in bilayer WS2/graphene heterostructure revealed by time-and angle-resolved photoemission spectroscopy. Front. Phys. 2021, 668149 (2021).
Zhou, H., Chen, Y. & Zhu, H. Deciphering asymmetric charge transfer at transition metal dichalcogenideâgraphene interface by helicity-resolved ultrafast spectroscopy. Sci. Adv. 7, 34 (2021).
Yuan, L. et al. Photocarrier generation from interlayer charge-transfer transitions in WS2-graphene heterostructures. Sci. Adv. 4, 1700324 (2018).
Fu, S. et al. Long-lived charge separation following pump-wavelengthâdependent ultrafast charge transfer in graphene/WS2 heterostructures. Sci. Adv. 7, 9 (2021).
Krause, R. et al. Microscopic understanding of ultrafast charge transfer in van-der-Waals heterostructures. Phys. Rev. Lett. 127, 276401 (2021).
Puppin, M. et al. Time-and angle-resolved photoemission spectroscopy of solids in the extreme ultraviolet at 500âkHz repetition rate. Rev. Sci. Instrum. 90, 023104 (2019).
Maklar, J. et al. A quantitative comparison of time-of-flight momentum microscopes and hemispherical analyzers for time-resolved ARPES experiments. Rev. Sci. Instrum. 91, 123112 (2020).
Li, Y. et al. Measurement of the optical dielectric function of monolayer transition-metal dichalcogenides: MoS2, MoSe2, WS2, and WSe2. Phys. Rev. B 90, 205422 (2014).
Stimper, V., Bauer, S., Ernstorfer, R., Schölkopf, B. & Xian, R. P. Multidimensional contrast limited adaptive histogram equalization. IEEE Access 7, 165437â165447 (2019).
Xian, R. P. et al. A machine learning route between band mapping and band structure. https://arxiv.org/abs/2005.10210 (2020).
Wallauer, R. et al. Momentum-resolved observation of exciton formation dynamics in monolayer WS2. Nano Lett. 21, 5867â5873 (2021).
Madéo, J. et al. Directly visualizing the momentum-forbidden dark excitons and their dynamics in atomically thin semiconductors. Science 370, 1199â1204 (2020).
Bostwick, A., Ohta, T., Seyller, T., Horn, K. & Rotenberg, E. Quasiparticle dynamics in graphene. Nat. Phys. 3, 36â40 (2007).
Nakamura, H. et al. Spin splitting and strain in epitaxial monolayer WSe2 on graphene. Phys. Rev. B 101, 165103 (2020).
Baudisch, M. et al. Ultrafast nonlinear optical response of Dirac fermions in graphene. Nat. Commun. 9, 1018 (2018).
Breusing, M., Ropers, C. & Elsaesser, T. Ultrafast carrier dynamics in graphite. Phys. Rev. Lett. 102, 086809 (2009).
Chen, Y., Li, Y., Zhao, Y., Zhou, H. & Zhu, H. Highly efficient hot electron harvesting from graphene before electron-hole thermalization. Sci. Adv. 5, 9958 (2019).
Dong, S. et al. Direct measurement of key exciton properties: energy, dynamics, and spatial distribution of the wave function. Nat. Sci. 1, e10010 (2021).
Chernikov, A., Ruppert, C., Hill, H. M., Rigosi, A. F. & Heinz, T. F. Population inversion and giant bandgap renormalization in atomically thin WS2 layers. Nat. Photonics 9, 466â470 (2015).
Liu, F., Ziffer, M. E., Hansen, K. R., Wang, J. & Zhu, X. Direct determination of band-gap renormalization in the photoexcited monolayer MoS2. Phys. Rev. Lett. 122, 246803 (2019).
Gao, S., Liang, Y., Spataru, C. D. & Yang, L. Dynamical excitonic effects in doped two-dimensional semiconductors. Nano Lett. 16, 5568â5573 (2016).
Dendzik, M. et al. Observation of an excitonic Mott transition through ultrafast core-cum-conduction photoemission spectroscopy. Phys. Rev. Lett. 125, 096401 (2020).
Steinhoff, A., Rosner, M., Jahnke, F., Wehling, T. O. & Gies, C. Influence of excited carriers on the optical and electronic properties of MoS2. Nano Lett. 14, 3743â3748 (2014).
Liang, Y. & Yang, L. Carrier plasmon induced nonlinear band gap renormalization in two-dimensional semiconductors. Phys. Rev. Lett. 114, 063001 (2015).
Tomadin, A., Brida, D., Cerullo, G., Ferrari, A. C. & Polini, M. Nonequilibrium dynamics of photoexcited electrons in graphene: Collinear scattering, Auger processes, and the impact of screening. Phys. Rev. B 88, 035430 (2013).
Gierz, I., Link, S., Starke, U. & Cavalleri, A. Non-equilibrium Dirac carrier dynamics in graphene investigated with time-and angle-resolved photoemission spectroscopy. Faraday Discuss. 171, 311â321 (2014).
Na, M. et al. Direct determination of mode-projected electron-phonon coupling in the time domain. Science 366, 1231â1236 (2019).
Meitner, L. On the origin of the β-ray spectra of radioactive substances. Z. Phys. 9, 131â144 (1922).
Auger, P. Sur les rayons β secondaires produits dans un gaz par des rayons X. CR Acad. Sci. (F.) 177, 169 (1923).
Selig, M. et al. Excitonic linewidth and coherence lifetime in monolayer transition metal dichalcogenides. Nat. Commun. 7, 13279 (2016).
Christiansen, D., Selig, M., Malic, E., Ernstorfer, R. & Knorr, A. Theory of exciton dynamics in time-resolved ARPES: intra-and intervalley scattering in two-dimensional semiconductors. Phys. Rev. B 100, 205401 (2019).
Selig, M. et al. Ultrafast dynamics in monolayer transition metal dichalcogenides: interplay of dark excitons, phonons, and intervalley exchange. Phys. Rev. Res. 1, 022007 (2019).
Pöllmann, C. et al. Resonant internal quantum transitions and femtosecond radiative decay of excitons in monolayer WSe2. Nat. Mater. 14, 889â893 (2015).
Förster, T. Intermolecular energy migration and fluorescence. Ann. Phys. 437, 55â75 (1948).
Kozawa, D. et al. Evidence for fast interlayer energy transfer in MoSe2/WS2 heterostructures. Nano Lett. 16, 4087â4093 (2016).
Paul, K. K., Kim, J.-H. & Lee, Y. H. Hot carrier photovoltaics in van der Waals heterostructures. Nat. Rev. Phys. 3, 178â192 (2021).
Riedl, C., Coletti, C. & Starke, U. Structural and electronic properties of epitaxial graphene on SiC (0 0 0 1): a review of growth, characterization, transfer doping and hydrogen intercalation. J. Phys. D: Appl. Phys. 43, 374009 (2010).
Mortensen, J. J., Hansen, L. B. & Jacobsen, K. W. Real-space grid implementation of the projector augmented wave method. Phys. Rev. B 71, 035109 (2005).
Gritsenko, O., van Leeuwen, R., van Lenthe, E. & Baerends, E. J. Self-consistent approximation to the kohn-sham exchange potential. Phys. Rev. A 51, 1944 (1995).
Xian, R. P. et al. An open-source, end-to-end workflow for multidimensional photoemission spectroscopy. Sci. Data 7, 442 (2020).
Acknowledgements
This work was funded by the Max Planck Society, the European Research Council (ERC) under the European Unionâs Horizon 2020 research and innovation program (Grant No. ERC-2015-CoG-682843), the German Research Foundation (DFG) within the Emmy Noether program (Grant No. RE 3977/1), through Projektnummer 18208777-SFB 951 âHybrid Inorganic/Organic Systems for Opto-Electronics (HIOS)â (CRC 951 project B12, M.S., D.C., A.K.), and the SFB/TRR 227 âUltrafast Spin Dynamicsâ (projects B07, project-ID: 328545488), and the Program DFG SPP2244 (project-ID: 443366970). S.B. acknowledges financial support from the NSERC-Banting Postdoctoral Fellowships Program. M.D. acknowledges financial support from the Göran Gustafsson Foundation and the Swedish Research Council under Grant No: 2022-03813. A.K. acknowledges financial support from DFG Projekt KN 427/14-1. A.C. and J.D.Z. acknowledge the financial support by the DFG SPP2244 (Project-ID: 443405595) and the Würzburg-Dresden Cluster of Excellence on Complexity and Topology in Quantum Matter (ct.qmat) (EXC 2147, Project-ID 390858490). K.W. and T.T. acknowledge support from the JSPS KAKENHI (Grant Numbers 21H05233 and 23H02052) and World Premier International Research Center Initiative (WPI), MEXT, Japan.
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Contributions
S.D., S.B., T.P., M.D., J.M., A.N., and L.R. performed the trARPES measurement. S.D. analyzed the data and wrote the first draft of the manuscript. R.E., L.R., and M.W. were responsible for developing all the experimental infrastructures. M.S. and D.C. performed the microscopic calculation with the guidance of A.K. R.P.X. and developed the 4D data processing code. P.R. and H.N. provided the epitaxially grown heterostructure, with support from U.S. and H.T. A.M., A.S., and M.S. conducted Raman and photoluminescence measurements, with guidance from M.J. and P.M. J.D.Z. and A.C. prepared the exfoliated ML sample with the hBN substrate provided by K.W. and T.T. All authors contributed to the final version of the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Hai Wang, and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Dong, S., Beaulieu, S., Selig, M. et al. Observation of ultrafast interfacial Meitner-Auger energy transfer in a Van der Waals heterostructure. Nat Commun 14, 5057 (2023). https://doi.org/10.1038/s41467-023-40815-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-40815-8
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.