# Fileset

[sciadv.ade3759.pdf](https://mdr.nims.go.jp/filesets/3e63d077-999c-4cc1-8231-bc0aca4ac998/download)

## Creator

Jingda Wu, Dongyang Yang, Jing Liang, Max Werner, Evgeny Ostroumov, Yunhuan Xiao, [Kenji Watanabe](https://orcid.org/0000-0003-3701-8119), [Takashi Taniguchi](https://orcid.org/0000-0002-1467-3105), Jerry I. Dadap, David Jones, Ziliang Ye

## Rights

[Creative Commons BY Attribution 4.0 International](https://creativecommons.org/licenses/by/4.0/)

## Other metadata

[Ultrafast response of spontaneous photovoltaic effect in 3R-MoS            <sub>2</sub>            –based heterostructures](https://mdr.nims.go.jp/datasets/c8f76296-9bfe-47c0-a5b6-977db03a7fc0)

## Fulltext

Ultrafast response of spontaneous photovoltaic effect in 3R-MoS2–based heterostructuresA P P L I E D  S C I E N C E S  A N D  E N G I N E E R I N G  Ultrafast response of spontaneous photovoltaic effect in 3R-MoS2–based heterostructures Jingda Wu1,2†, Dongyang Yang1,2†, Jing Liang1,2†, Max Werner1,2, Evgeny Ostroumov1,2,  Yunhuan Xiao1,2, Kenji Watanabe3, Takashi Taniguchi4, Jerry I. Dadap1,2, David Jones1,2,  Ziliang Ye1,2* Rhombohedrally stacked MoS2 has been shown to exhibit spontaneous polarization down to the bilayer limit and can sustain a strong depolarization field when sandwiched between graphene. Such a field gives rise to a spontaneous photovoltaic effect without needing any p-n junction. In this work, we show that the photovoltaic effect has an external quantum efficiency of 10% for devices with only two atomic layers of MoS2 at low tem-peratures, and identify a picosecond-fast photocurrent response, which translates to an intrinsic device band-width at ∼100-GHz level. To this end, we have developed a nondegenerate pump-probe photocurrent spectroscopy technique to deconvolute the thermal and charge-transfer processes, thus successfully revealing the multicomponent nature of the photocurrent dynamics. The fast component approaches the limit of the charge-transfer speed at the graphene-MoS2 interface. The remarkable efficiency and ultrafast photoresponse in the graphene-3R-MoS2 devices support the use of ferroelectric van der Waals materials for future high-per-formance optoelectronic applications. Copyright © 2022  The Authors, some  rights reserved;  exclusive licensee  American Association  for the Advancement  of Science. No claim to  original U.S. Government  Works. Distributed  under a Creative  Commons Attribution  License 4.0 (CC BY).   INTRODUCTION When two layers of transition metal dichalcogenides (TMDs) are stacked in parallel, both the inversion symmetry and mirror symme-try are spontaneously broken, leading to an electric polarization along the out-of-plane direction (1). Under an electric field larger than the coercive field, the polarization direction can be flipped as one layer slides relative to the neighboring layer, which has been termed sliding ferroelectricity (2–5). While the artificially stacked bilayers (BLs) have limited domain sizes due to twist misalignment, homogeneous polarization in a scale as large as the entire flake has been observed in BLs directly exfoliated from a 3R-MoS2 bulk crystal (6), making it ideal for certain optoelectronic applications. In particular, when 3R-MoS2 is sandwiched between graphene (7), the spontaneous polarization gives rise to a photovoltaic effect, known as spontaneous photovoltaic effect (8), where photo-excited carriers in MoS2 transfer asymmetrically to the graphene under a largely unscreened depolarization field. An external quantum efficiency (EQE) of 16% has been observed in a device with 10-layer 3R-MoS2. The high EQE in the graphene/3R-MoS2/graphene heterostruc-tures suggests a potentially ultrafast photocurrent dynamics compa-rable to the charge-transfer process at the MoS2-graphene interface, which lasts approximately 1 ps (9–13). Such an ultrafast photocur-rent, if confirmed, can potentially be leveraged for high-speed optical communications. Although the MoS2 photocarrier dynam-ics generally can be measured by ultrafast techniques such as optical pump-probe spectroscopy (14), it remains a challenge to probe the intrinsic dynamics of the aforementioned photovoltaic effect, due to the existence of multiple photocurrent contributions. Because the spontaneous polarization naturally induces image charges of oppo-site polarities at the two electrodes, both the bolometric and the photo-thermoelectric (PTE) effects exist in the device (15–20). In addition, the asymmetric tunneling barrier is also known for gen-erating a photocurrent without external bias (21, 22). Here, we dis-entangle these electronic and thermal contributions and measure their corresponding response times by performing ultrafast photo-current autocorrelation and nondegenerate pump-probe photocur-rent measurements. We conclude that the photocurrent has approximately 2-ps fast dynamics for devices of various thicknesses and a slow dynamics that lasts for 25 ps or more. Ultrafast optical techniques are powerful tools for probing the photocarrier dynamics (9, 14, 23). Photocurrent autocorrelation technique has been previously used to study the photocurrent dy-namics in two-dimensional (2D) material–based photodetectors (24–28). In this approach, two equal-power laser pulses are split from a single femtosecond pulse and illuminate the heterostructure to generate a photocurrent (Fig. 1A). When these two pulses coin-cide on the device, a saturation in the photocurrent can be observed due to either electronic or thermal mechanisms. For example, elec-tronic saturation may arise from absorption saturation and in-creased interactions between carriers, whereas thermal saturation may be induced by the nonlinear temperature increase under illu-mination (29–32). By scanning the delay between these two pulses, a characteristic time associated with the saturation mechanism can be measured. When multiple mechanisms contribute to the photocur-rent, the intrinsic photocurrent time response can be obscured. RESULTS AND DISCUSSION We first perform the photocurrent autocorrelation measurement on a heterostructure composed of both BL and four-layer (4L) 3R- MoS2. The optical image of the heterostructure is shown in fig. 1Department of Physics and Astronomy, The University of British Columbia, 6224 Agricultural Road, Vancouver, BC V6T 1Z1, Canada. 2Stewart Blusson Quantum Matter Institute, The University of British Columbia, 2355 East Mall, Vancouver, BC V6T 1Z4, Canada. 3Research Center for Functional Materials, National Institute for Materials Science, 1-1 Namiki, Tsukuba 305-0044, Japan. 4International Center for Materials Nanoarchitectonics, National Institute for Materials Science, 1-1 Namiki, Tsukuba 305-0044, Japan. *Corresponding author. Email: zlye@phas.ubc.ca †These authors contributed equally to this work. Wu et al., Sci. Adv. 8, eade3759 (2022) 16 December 2022                                                                                                                                                        1 of 7  S C I E N C E  A D VA N C E S | R E S E A R C H  A R T I C L E  Downloaded from https://www.science.org at National Institute for Materials Science on February 04, 2023S1. From the current-voltage characteristics shown in Fig. 1B, a spontaneous photovoltaic behavior is clearly observed. The zero- bias photocurrent mapping shows a finite photoresponse in BL and 4L regions, and the photocurrent distribution agrees with the MoS2 thickness in the graphene overlapped area (Fig. 1C; photocur-rent maps for other devices can be found in fig. S2). This photocur-rent grows sublinearly with optical power (fig. S3), giving rise to a negative autocorrelation signal. A representative autocorrelation result is shown in Fig. 1D. A dip of the photocurrent strength is ob-served at zero delay, with symmetric recovery dynamics at both pos-itive and negative delays, as expected for conventional photocurrent autocorrelation signals (24–28). More quantitative analysis of the autocorrelation signal in the BL area is presented in Fig. 2. With a fluence of 140 μJ cm−2 for a spot size of approximately 1 μm in diameter at 770 nm (50-fs pulse width) for each pulse, the saturation dip at zero time delay is ap-proximately 15% of the steady-state signal. The recovery time is fitted to be 17 ps, more than one order of magnitude slower than the charge-transfer time at the MoS2-graphene interface (10, 12), suggesting that the saturation mechanism is likely not electronic. If such a slowing is caused by the charge transfer at the MoS2- MoS2 interface, a previous study suggests that an external electric field can speed up the process when the TMD is not monolayer (24). However, when we apply a bias voltage from −0.5 to 0.5 V across the device, no substantial change is observed in the recovery time (Fig. 2A, inset; the raw data are shown in fig. S4). Instead, we suggest that this photocurrent saturation originates from photo-thermal effects. As shown in Fig. 2C, the photoresponse of our device has a surprisingly strong temperature dependence. When the heterostructure is heated by 10 K above room tempera-ture (RT), the photocurrent drops by more than 25% (Fig. 2C, inset). More notably, the photoresponsivity increases by more than one order of magnitude after cooling the heterostructure from RT to 3 K (Fig. 2C, see also in another device in fig. S5, full temperature range in fig. S6A). At the A-exciton resonance, the EQE approaches 10%, which is remarkable for an atomically thin device. In addition, we also observe a pronounced temperature de-pendence in the total resistance of the device (fig. S6B), which mo-tivates us to develop a shunt-resistance model. In this model, within a few tens of kelvin above RT, which is relevant to our pump-probe photocurrent measurement, only the shunt resistance is expected to change markedly due to its tunneling current origin (33). By mod-eling the measured photocurrent to be proportional to the shunt resistance, which exponentially decreases with increasing tempera-ture, we are able to fit the temperature dependence of the photocur-rent (Fig. 2C, inset, details in the Supplementary Materials). Other factors that can contribute to this strong temperature dependence include the variation in spontaneous polarization (34, 35) and Fig. 1. Photocurrent generation in a Gr/3R-MoS2/Gr device. (A) Schematic of the device and ultrafast time-resolved photocurrent measurement. The heterostructure consists of a bilayer (BL) or few-layer 3R-MoS2 flake that is sandwiched between two graphene electrodes. The whole device is encapsulated with hexagonal boron nitride (hBN). The photocurrent is measured with a lock-in amplifier through graphene electrodes. The two ultrafast pulses are separated by a time delay Δτ and focused colli-nearly on the device. (B) Photocurrent-voltage characteristics on a BL device at room temperature (RT). The device shows photovoltaic effect with nonzero photocurrent without bias. (C) Photocurrent map of a device consists of both BL and 4L regions at zero bias. The resolution of the scan is 400 nm in both directions. The darker blue regions are also BL but with a different polarity. Scale bar, 2 μm. Both (B) and (C) are measured under 20-μW 532-nm continuous-wave (CW) laser illumination. (D) Autocorrelation signal from the BL and 4L region with 770-nm femtosecond pulse excitation, which is above the indirect bandgap of the material. Wu et al., Sci. Adv. 8, eade3759 (2022) 16 December 2022                                                                                                                                                        2 of 7  S C I E N C E  A D VA N C E S | R E S E A R C H  A R T I C L E  Downloaded from https://www.science.org at National Institute for Materials Science on February 04, 2023change in photocarrier recombination time or exciton linewidth, which affects the absorption of the heterostructure. As a result, when an ultrafast laser pulse is absorbed by the device, the transient temperature of the illuminated region quickly rises, which decreases the shunt resistance and subsequently decreases the photocurrent generated by the following pulse. In this case, the saturation recovery time corresponds to the device heat dissipation time. At low temperatures, the photocurrent becomes less temperature de-pendent (fig. S6), which is consistent with the smaller autocorrela-tion signal observed at 3 K, possibly due to the internal quantum efficiency limit. We confirm the photo-thermal origin of the saturation by mea-suring the photocurrent autocorrelation at cryogenic temperature (Fig. 2B). After most phonon modes are frozen at low temperature, the thermal conductivity of the 2D materials decreases substantially (36, 37), which should lead to a longer heat dissipation time. Exper-imentally, we find a recovery time that is approximately three times longer (55 ps) than that at RT, when the heterostructure is cooled to 3 K. All these evidences suggest that the dynamics observed in the autocorrelation measurements is limited by the device heat dissipation. In addition to this long recovery time, we also observe a transient middle peak in the autocorrelation measurement at 3 K. Exponential fitting of that small peak suggests a time constant of approximately 4 ps. Similar peaks have been observed in the 4L region (fig. S7A) and on other BL devices of varying amplitudes up to RT (fig. S7B). We interpret this feature as a result of the electronic temperature saturation in the graphene electrodes. As discussed in our previous work (7), the top and bottom graphene electrodes are naturally doped with charges of opposite polarities in these devices, which consequently acquire opposite Seebeck coefficients. Upon laser heating, a PTE current in a direction opposite to the photovol-taic current is generated, which reduces the net photocurrent. In the autocorrelation experiment, the two subsequent laser pulses with a small delay time saturate the graphene electronic temperature, which decreases the PTE effect and therefore increases the overall photocurrent. The characteristic time of the middle peak also agrees with graphene electronic temperature relaxation through supercollision (27, 28, 38). Because the photo-thermal effect of the laser pulse obscures the electronic response speed in autocorrelation measurements, an al-ternative method is needed to separate the thermal and electronic contributions. Here, we develop a nondegenerate pump-probe pho-tocurrent spectroscopy technique that permits the determination of the intrinsic speed of our devices (Fig. 3A). First, we use a strong sub-bandgap infrared pulse (1030 nm, 50 fs, 350 μJ cm−2; IR- pulse) to heat the graphene electrodes, which subsequently heats the MoS2 region. Second, we use a weak visible pulse (670 nm, 50 fs, 2 μJ cm−2; VIS-pulse) that is resonant with the A-exciton of MoS2 to generate a photocurrent pulse but without generating suffi-cient heat. At negative delays when the VIS-pulse arrives earlier than the IR-pulse, despite the fact that the IR-pulse generates a finite photo-current, the photo-thermal saturation effect should be negligible because little heat is deposited by the VIS-pulse (details in the Sup-plementary Materials). The photocurrent should decrease consider-ably immediately after the zero delay because the substantial heating induced by the IR-pulse would reduce the photocurrent generated by the VIS-pulse, in accordance with the results from our autocor-relation experiments. At a large positive delay, the photocurrent should be restored to the steady state as the heat is fully dissipated. Therefore, as we scan the time delay, we expect to see an asymmetric pump-probe photocurrent signal similar to photocarrier dynamics from an optical pump-probe measurement, which exhibits a quick rise followed by an exponential decrease in signal. Because the Fig. 2. Temperature-dependent measurements in the BL region. Autocorrela-tion signals of the BL region measured at (A) RT and (B) 3 K. The dotted lines are from experimental measurements, and black lines are exponential fittings to extract the time constants. Inset of (A): photocurrent time constant extracted from bias voltage–dependent autocorrelation measurements. The dashed line in-dicates a photocurrent time constant, τ, of 17 ps. (C) Energy-dependent photores-ponsivity of the BL region measured at RT and 3 K. A and B peaks correspond to the A-exciton and B-exciton absorption features in BL 3R-MoS2. The average power is kept at around 0.5 μW at all wavelengths to avoid saturation. Inset: photocurrent of the BL region at temperatures above RT with 532-nm CW laser excitation. The dots are experimental data, and the solid line is a fit based on the shunt-resis-tance model. Wu et al., Sci. Adv. 8, eade3759 (2022) 16 December 2022                                                                                                                                                        3 of 7  S C I E N C E  A D VA N C E S | R E S E A R C H  A R T I C L E  Downloaded from https://www.science.org at National Institute for Materials Science on February 04, 2023pump-probe photocurrent dynamics is approximately a cross-cor-relation between the photocurrent and device temperature dynam-ics, the decay of the signal is likely limited by the heat dissipation, while the rise of the signal should be determined by the photocur-rent decay and lattice heating. Experimentally, a highly asymmetric temporal response is ob-served in the pump-probe photocurrent measurement (Fig. 3B). The pump-probe photocurrent signal captures the photocurrent change due to thermal saturation. As expected, we find a very sharp drop of the photocurrent around the zero delay. The drop lasts for approximately 4 ps (Fig. 3B, inset), which gives an upper bound for the characteristic time of the transient photocurrent. It is followed by a 20-ps slow recovery dynamics similar to that ob-served in the autocorrelation experiment, corresponding to the cooling of the device. These pump-probe photocurrent results confirm the existence of a picosecond photocurrent response. In ad-dition, a small photocurrent saturation is observed ahead of the sharp drop before zero delay, suggesting that multiple dynamics could be associated with the photocurrent. The delayed IR-pulse can affect the photocurrent generated earlier by the VIS-pulse if the photocurrent lasts long enough to overlap with the heating pulse. With the fast photocurrent component determined from the quick drop, the saturation long before the zero delay indicates the existence of a slow photocurrent component. In addition, after symmetrizing the pump-probe photocurrent signal by adding itself to the time-reversed copy, the result shows similar dynamics as the photocurrent autocorrelation signal, except that the autocorrelation is affected by the PTE effect near zero delay (details in the Supple-mentary Materials). Since the signal in the pump-probe photocurrent experiment is much larger at positive delay than that at negative delay, we confirm that the slow decay in the auto-correlation signal is dominated by the photo-thermal satura-tion effect. To have a quantitative understanding of the photocurrent dy-namics, a phenomenological model is developed to model the joint response of the device temperature evolution (thermal pulse) and transient photocurrent ( photocurrent pulse). In this model, each VIS-pulse generates two concurrent photocurrent pulses in the device with independent decay constants I(t) = I1e−t/t1 + I2e−t/t2. Upon overlapping with the thermal pulse, the in-stantaneous photocurrent is reduced to αI, where α is a dimension-less saturation factor and its temperature dependence is defined in the Supplementary Materials. As the device temperature is a func-tion of time, α is consequently time dependent. After calculating the evolution of the transient lattice temperature of the device excited by an IR-pulse at a delay of Δτ through a two-temperature model (Supplementary Materials), we can obtain the average photocurrent through the cross-correlation relation IðDtÞ ¼ 1trepÐ trep0 aðtÞIðt þ DtÞdt, where I is the average photocur-rent and τrep is the repetition time of the laser pulse, which is much longer than the photocurrent response time. The calculated I fits our experimental results well (Fig. 3B), with the transient pho-tocurrent profile I(t) presented in Fig. 3C. The fast decay in the early stage has a time constant t1 of approximately 2 ps, which is similar to the graphene/TMD charge-transfer time. The measured time cons-tant is nearly independent of the power of IR-pulse or VIS-pulse. This characteristic time of the transient photocurrent approaches that of the charge transfer at the graphene-TMD interface [9–13], Fig. 3. Pump-probe photocurrent measurement. (A) Illustration of the heating dynamics of the device by a sub-bandgap pulse (IR-pulse). Graphene is first heated by the IR-pulse and then heats MoS2. The photocurrent is generated by the VIS-pulse, and it becomes weaker when the device is thermally saturated upon heating. (B) Pump-probe photocurrent measurement on a BL device. The red dotted line represents experimental results, and the solid black curve is the fitting from the numerical model. The measurements are carried out with an IR-pulse(VIS-pulse) fluence of 350(2) μJ cm−2. The pump-probe photocurrent signal is normalized for comparison with the numerical model. The 20-ps scan around zero delay is performed at a finer temporal resolution than the rest of the delay range. Inset on the left: a zoom-in semilog plot of the pump-probe photocurrent signal near zero delay. The black straight lines are for guidance. (C) Transient photocurrent responses for different processes. The total profile consists of two photocurrent generation processes, with the fast and slow processes having characteristic times of 2 and 25 ps, respectively.  Wu et al., Sci. Adv. 8, eade3759 (2022) 16 December 2022                                                                                                                                                        4 of 7  S C I E N C E  A D VA N C E S | R E S E A R C H  A R T I C L E  Downloaded from https://www.science.org at National Institute for Materials Science on February 04, 2023suggesting a likely faster interlayer charge-transfer process at the MoS2 interface, similar to those observed in TMD heterostructures (9). On the other hand, the slow component has a time constant t2 of approximately 25 ps. It was previously found that the two-com-ponent photocurrent can originate from defect-related processes (30, 32, 39), which agrees with the dynamics we observe. To further confirm the electronic and thermal contributions in the photocurrent dynamics, a comparison is made between the BL and 4L regions with the pump-probe photocurrent measurements. Similar to the BL region, an asymmetric time dependence before and after zero delay as well as a sharp drop in photocurrent around zero delay are observed in the 4L region (Fig. 4A, isolated curve in fig. S8), indicating a similarly fast dynamics at the gra-phene-MoS2 interface. Nevertheless, the pump-probe photocurrent signal before zero delay is much more prominent in the 4L than BL region, suggesting a much larger slow component contribution in the thick area. This contrast might result from a larger defect density in 2D or the extra MoS2-MoS2 interfaces in the 4L region. An exter-nal bias can be applied to tune the slow component of the photo-current, which is observed before zero delay (Fig. 4A). As fitted by the aforementioned model, the slow component becomes faster when the electric field is along the photocurrent direction and vice versa (Fig. 4B). On the other hand, the recovery time at positive delay is dominated by the device cooling dynamics and is therefore largely bias independent. These results further confirm the thermal saturation nature and the existence of a two-component photocur-rent response in the 3R-MoS2 photovoltaic device. In conclusion, we have observed a large temperature dependence in the spontaneous photovoltaic effect in heterostructures compris-ing atomically thin 3R-MoS2 and graphene. At low temperatures, the EQE can reach 10% even in the thinnest area with only two atomic layers of MoS2. Upon laser heating, such a temperature de-pendence leads to a strong thermal saturation of the photocurrent, which dominates in the commonly used photocurrent autocorrela-tion measurement and obscures the photocurrent dynamics. In ad-dition, we have developed a nondegenerate pump-probe photocurrent spectroscopy technique, which can distinguish the electronic and thermal dynamics. With this measurement tech-nique, we find that the transient photocurrent has two contribu-tions with distinct temporal responses. We quantify the fast photocurrent response to be on the picosecond level, which is similar to the charge-transfer process at the graphene-TMD inter-face, suggesting an intrinsic device bandwidth of hundreds of giga-hertz. At RT, the photoresponsivity and speed of our homobilayer device are comparable to those of TMD heterobilayer devices with type II band alignment (9, 40, 41). Our results may stimulate future uses of ferroelectric van der Waals devices in optoelectronic appli-cations that require high performance with low power consumption and built-in memory function, such as in high-speed optical com-munications and optical computing. Furthermore, the nondegener-ate pump-probe photocurrent spectroscopy technique can also be applied to study the interplay between various photo-excitation mechanisms in other types of photosensors. MATERIALS AND METHODS Sample fabrication The 3R-MoS2 flakes are exfoliated from a bulk crystal (HQ Gra-phene) on SiO2/Si substrates. Graphene electrodes are cut into spe-cific shapes with a femtosecond laser. The flake thickness is identified by an optical microscope and reflection contrast spectro-scopy. The 3R phase is confirmed by second-harmonic generation measurement (42). Exfoliated hexagonal boron nitride (hBN) is used for encapsulation. All devices are fabricated using dry transfer method under ambient conditions with a homebuilt transfer stage. Electrical contact is achieved by overlapping the graphene with gold electrodes prepatterned by optical lithography on heavily p-doped Si/SiO2 substrates. Photocurrent measurement For the autocorrelation measurements, a homebuilt wavelength- tunable optical parametric oscillator (OPO) (76-MHz repetition rate, 50-fs pulse width) is used to illuminate the sample, which sits in a continuous-flow optical cryostat (Oxford Microstat-He), through a long working distance 100× objective lens (Mitutoyo, nu-merical aperture = 0.5). The signal is detected through a lock-in am-plifier phase-locked to a mechanical chopper. In the pump-probe photocurrent measurement, we split some light from the pump laser (light conversion, 1030 nm) of the OPO and then combine it with the OPO beam for collinear measurements. The photocur-rent is measured by modulating the VIS-pulse from the OPO. For wavelength-dependent measurements, a supercontinuum white light source (YSL) (∼100-ps pulse width, 10 MHz) is wavelength- selected by an acousto-optic tunable filter to illuminate the Fig. 4. Pump-probe photocurrent measurement in the 4L region. (A) Bias-de-pendent pump-probe photocurrent measurement in the 4L region with bias volt-ages ranging from −0.25 to 0.5 V. The measurements are carried out with an IR- pulse(VIS-pulse) fluence of approximately 350(2) μJ cm−2. (B) Numerical modeling results with the response time of the slow component varying from 35 to 50 ps. Wu et al., Sci. Adv. 8, eade3759 (2022) 16 December 2022                                                                                                                                                        5 of 7  S C I E N C E  A D VA N C E S | R E S E A R C H  A R T I C L E  Downloaded from https://www.science.org at National Institute for Materials Science on February 04, 2023device. For device photocurrent mapping, current-voltage charac-teristics, and temperature dependence, a 532-nm continuous- wave laser is used. Supplementary Materials This PDF file includes: Figs. S1 to S14 Notes S1 to S7 References REFERENCES AND NOTES  1. J. Sung, Y. Zhou, G. Scuri, V. Zólyomi, T. I. Andersen, H. Yoo, D. S. Wild, A. Y. Joe, R. J. Gelly, H. Heo, S. J. Magorrian, D. Bèrubè, A. M. M. Valdivia, T. Taniguchi, K. Watanabe, M. D. Lukin, P. Kim, V. I. Fal’ko, H. Park, Broken mirror symmetry in excitonic response of reconstructed domains in twisted MoSe2/MoSe2 bilayers. Nat. Nanotechnol. 15, 750–754 (2020).  2. X. Wang, K. Yasuda, Y. Zhang, S. Liu, K. Watanabe, T. Taniguchi, J. Hone, L. Fu, P. Jarillo- Herrero, Interfacial ferroelectricity in rhombohedral-stacked bilayer transition metal di-chalcogenides. Nat. Nanotechnol. 17, 367–371 (2022).  3. M. Vizner Stern, Y. Waschitz, W. Cao, I. Nevo, K. Watanabe, T. Taniguchi, E. Sela, M. Urbakh, O. Hod, M. Ben Shalom, Interfacial ferroelectricity by van der Waals sliding. Science 372, 1462–1466 (2021).  4. M. Wu, J. Li, Sliding ferroelectricity in 2D van der Waals materials: Related physics and future opportunities. Proc. Natl. Acad. Sci. U.S.A. 118, e2115703118 (2021).  5. A. Weston, E. G. Castanon, V. Enaldiev, F. Ferreira, S. Bhattacharjee, S. Xu, H. Corte-León, Z. Wu, N. Clark, A. Summerfield, T. Hashimoto, Y. Gao, W. Wang, M. Hamer, H. Read, L. Fumagalli, A. V. Kretinin, S. J. Haigh, O. Kazakova, A. K. Geim, V. I. Fal’ko, R. Gorbachev, Interfacial ferroelectricity in marginally twisted 2D semiconductors. Nat. Nanotechnol. 17, 390–395 (2022).  6. J. Liang, D. Yang, J. Wu, J. I. Dadap, K. Watanabe, T. Taniguchi, Z. Ye, Optically probing the asymmetric interlayer coupling in rhombohedral-stacked MoS2 bilayer. Phys. Rev. X 12, 041005 (2022).  7. D. Yang, J. Wu, B. T. Zhou, J. Liang, T. Ideue, T. Siu, K. M. Awan, K. Watanabe, T. Taniguchi, Y. Iwasa, M. Franz, Z. Ye, Spontaneous-polarization-induced photovoltaic effect in rhom-bohedrally stacked MoS2. Nat. Photon. 16, 469–474 (2022).  8. T. Akamatsu, T. Ideue, L. Zhou, Y. Dong, S. Kitamura, M. Yoshii, D. Yang, M. Onga, Y. Nakagawa, K. Watanabe, T. Taniguchi, J. Laurienzo, J. Huang, Z. Ye, T. Morimoto, H. Yuan, Y. Iwasa, A van der Waals interface that creates in-plane polarization and a spontaneous photovoltaic effect. Science 372, 68–72 (2021).  9. X. Hong, J. Kim, S.-F. Shi, Y. Zhang, C. Jin, Y. Sun, S. Tongay, J. Wu, Y. Zhang, F. Wang, Ultrafast charge transfer in atomically thin MoS2/WS2 heterostructures. Nat. Nanotechnol. 9, 682–686 (2014). 10. L. Yuan, T.-F. Chung, A. Kuc, Y. Wan, Y. Xu, Y. P. Chen, T. Heine, L. Huang, Photocarrier generation from interlayer charge-transfer transitions in WS2-graphene heterostructures. Sci. Adv. 4, e1700324 (2018). 11. X. Zhang, D. He, L. Yi, S. Zhao, J. He, Y. Wang, H. Zhao, Electron dynamics in MoS2-graphite heterostructures. Nanoscale 9, 14533–14539 (2017). 12. D. Luo, J. Tang, X. Shen, F. Ji, J. Yang, S. Weathersby, M. E. Kozina, Z. Chen, J. Xiao, Y. Ye, T. Cao, G. Zhang, X. Wang, A. M. Lindenberg, Twist-angle-dependent ultrafast charge transfer in MoS2-graphene van der Waals heterostructures. Nano Lett. 21, 8051–8057 (2021). 13. Y. Chen, Y. Li, Y. Zhao, H. Zhou, H. Zhu, Highly efficient hot electron harvesting from gra-phene before electron-hole thermalization. Sci. Adv. 5, eaax9958 (2019). 14. C. Jin, E. Y. Ma, O. Karni, E. C. Regan, F. Wang, T. F. Heinz, Ultrafast dynamics in van der Waals heterostructures. Nat. Nanotechnol. 13, 994–1003 (2018). 15. J.-Y. Wu, Y. T. Chun, S. Li, T. Zhang, J. Wang, P. K. Shrestha, D. Chu, Broadband MoS2 field- effect phototransistors: Ultrasensitive visible-light photoresponse and negative infrared photoresponse. Adv. Mater. 30, 1705880 (2018). 16. W. Zhang, C.-P. Chuu, J.-K. Huang, C.-H. Chen, M.-L. Tsai, Y.-H. Chang, C.-T. Liang, Y.-Z. Chen, Y.-L. Chueh, J.-H. He, M.-Y. Chou, L.-J. Li, Ultrahigh-gain photodetectors based on atomically thin graphene-MoS2 heterostructures. Sci. Rep. 4, 3826 (2014). 17. M. Buscema, M. Barkelid, V. Zwiller, H. S. van der Zant, G. A. Steele, A. Castellanos-Gomez, Large and tunable photothermoelectric effect in single-layer MoS2. Nano Lett. 13, 358–363 (2013). 18. D. J. Groenendijk, M. Buscema, G. A. Steele, S. Michaelis de Vasconcellos, R. Bratschitsch, H. S. J. van der Zant, A. Castellanos-Gomez, Photovoltaic and photothermoelectric effect in a double-gated WSe2 device. Nano Lett. 14, 5846–5852 (2014). 19. N. M. Gabor, J. C. Song, Q. Ma, N. L. Nair, T. Taychatanapat, K. Watanabe, T. Taniguchi, L. S. Levitov, P. Jarillo-Herrero, Hot carrier–assisted intrinsic photoresponse in graphene. Science 334, 648–652 (2011). 20. M. Freitag, T. Low, F. Xia, P. Avouris, Photoconductivity of biased graphene. Nat. Photon. 7, 53–59 (2013). 21. W. J. Yu, Q. A. Vu, H. Oh, H. G. Nam, H. Zhou, S. Cha, J.-Y. Kim, A. Carvalho, M. Jeong, H. Choi, A. H. C. Neto, Y. H. Lee, X. Duan, Unusually efficient photocurrent extraction in monolayer van der Waals heterostructure by tunnelling through discretized barriers. Nat. Commun. 7, 13278 (2016). 22. L. Hou, Q. Zhang, M. Tweedie, V. Shautsova, Y. Sheng, Y. Zhou, H. Huang, T. Chen, J. H. Warner, Photocurrent direction control and increased photovoltaic effects in all-2D ultrathin vertical heterostructures using asymmetric h-BN tunneling barriers. ACS Appl. Mater. Interfaces 11, 40274–40282 (2019). 23. D. Daranciang, M. J. Highland, H. Wen, S. M. Young, N. C. Brandt, H. Y. Hwang, M. Vattilana, M. Nicoul, F. Quirin, J. Goodfellow, T. Qi, I. Grinberg, D. M. Fritz, M. Cammarata, D. Zhu, H. T. Lemke, D. A. Walko, E. M. Dufresne, Y. Li, J. Larsson, D. A. Reis, K. Sokolowski-Tinten, K. A. Nelson, A. M. Rappe, P. H. Fuoss, G. B. Stephenson, A. M. Lindenberg, Ultrafast pho-tovoltaic response in ferroelectric nanolayers. Phys. Rev. Lett. 108, 087601 (2012). 24. M. Massicotte, P. Schmidt, F. Vialla, K. G. Schädler, A. Reserbat-Plantey, K. Watanabe, T. Taniguchi, K. J. Tielrooij, F. H. L. Koppens, Picosecond photoresponse in van der Waals heterostructures. Nat. Nanotechnol. 11, 42–46 (2016). 25. M. Massicotte, P. Schmidt, F. Vialla, K. Watanabe, T. Taniguchi, K. J. Tielrooij, F. H. L. Koppens, Photo-thermionic effect in vertical graphene heterostructures. Nat. Commun. 7, 12174 (2016). 26. D. Sun, G. Aivazian, A. M. Jones, J. S. Ross, W. Yao, D. Cobden, X. Xu, Ultrafast hot-carrier- dominated photocurrent in graphene. Nat. Nanotechnol. 7, 114–118 (2012). 27. M. W. Graham, S.-F. Shi, D. C. Ralph, J. Park, P. L. McEuen, Photocurrent measurements of supercollision cooling in graphene. Nat. Phys. 9, 103–108 (2013). 28. A. C. Betz, S. H. Jhang, E. Pallecchi, R. Ferreira, G. Fève, J. M. Berroir, B. Plaçais, Supercollision cooling in undoped graphene. Nat. Phys. 9, 109–112 (2013). 29. Z. Nie, R. Long, L. Sun, C.-C. Huang, J. Zhang, Q. Xiong, D. W. Hewak, Z. Shen, O. V. Prezhdo, Z.-H. Loh, Ultrafast carrier thermalization and cooling dynamics in few-layer MoS2. ACS Nano 8, 10931–10940 (2014). 30. K. T. Vogt, S.-F. Shi, F. Wang, M. W. Graham, Ultrafast photocurrent and absorption mi-croscopy of few-layer transition metal dichalcogenide devices that isolate rate-limiting dynamics driving fast and efficient photoresponse. J. Phys. Chem. C 124, 15195–15204 (2020). 31. D. Sun, Y. Rao, G. A. Reider, G. Chen, Y. You, L. Brézin, A. R. Harutyunyan, T. F. Heinz, Ob-servation of rapid exciton–exciton annihilation in monolayer molybdenum disulfide. Nano Lett. 14, 5625–5629 (2014). 32. H. Wang, C. Zhang, F. Rana, Ultrafast dynamics of defect-assisted electron–hole recombi-nation in monolayer MoS2. Nano Lett. 15, 339–345 (2015). 33. T. Georgiou, R. Jalil, B. D. Belle, L. Britnell, R. V. Gorbachev, S. V. Morozov, Y.-J. Kim, A. Gholinia, S. J. Haigh, O. Makarovsky, L. Eaves, L. A. Ponomarenko, A. K. Geim, K. S. Novoselov, A. Mishchenko, Vertical field-effect transistor based on graphene–WS2 heterostructures for flexible and transparent electronics. Nat. Nanotechnol. 8, 100–103 (2013). 34. K. Yasuda, X. Wang, K. Watanabe, T. Taniguchi, P. Jarillo-Herrero, Stacking-engineered fer-roelectricity in bilayer boron nitride. Science 372, 1458–1462 (2021). 35. Y. Liu, S. Liu, B. Li, W. J. Yoo, J. Hone, Identifying the transition order in an artificial ferro-electric van der Waals heterostructure. Nano Lett. 22, 1265–1269 (2022). 36. Y. Liu, Z.-Y. Ong, J. Wu, Y. Zhao, K. Watanabe, T. Taniguchi, D. Chi, G. Zhang, J. T. Thong, C.- W. Qiu, K. Hippalgaonkar, Thermal conductance of the 2D MoS2/h-BN and graphene/h-BN interfaces. Sci. Rep. 7, 1–8 (2017). 37. A. A. Balandin, Thermal properties of graphene and nanostructured carbon materials. Nat. Mater. 10, 569–581 (2011). 38. J. C. Song, M. Y. Reizer, L. S. Levitov, Disorder-assisted electron-phonon scattering and cooling pathways in graphene. Phys. Rev. Lett. 109, 106602 (2012). 39. M. M. Furchi, D. K. Polyushkin, A. Pospischil, T. Mueller, Mechanisms of photoconductivity in atomically thin MoS2. Nano Lett. 14, 6165–6170 (2014). 40. C.-H. Lee, G.-H. Lee, A. M. Van Der Zande, W. Chen, Y. Li, M. Han, X. Cui, G. Arefe, C. Nuckolls, T. F. Heinz, J. Guo, J. Hone, P. Kim, Atomically thin p–n junctions with van der Waals het-erointerfaces. Nat. Nanotechnol. 9, 676–681 (2014). 41. E. Y. Ma, B. Guzelturk, G. Li, L. Cao, Z.-X. Shen, A. M. Lindenberg, T. F. Heinz, Recording in-terfacial currents on the subnanometer length and femtosecond time scale by terahertz emission. Sci. Adv. 5, eaau0073 (2019). 42. M. Zhao, Z. Ye, R. Suzuki, Y. Ye, H. Zhu, J. Xiao, Y. Wang, Y. Iwasa, X. Zhang, Atomically phase- matched second-harmonic generation in a 2D crystal. Light Sci. Appl. 5, e16131 (2016). 43. S. M. Sze, Y. Li, K. K. Ng, Physics of Semiconductor Devices (John Wiley & Sons, 2021). Wu et al., Sci. Adv. 8, eade3759 (2022) 16 December 2022                                                                                                                                                        6 of 7  S C I E N C E  A D VA N C E S | R E S E A R C H  A R T I C L E  Downloaded from https://www.science.org at National Institute for Materials Science on February 04, 202344. D. Moldovan, M. Andelkovic, F. Peeters, pybinding v0.9.5: A Python package for tight- binding calculations (2020). 45. E. Pop, V. Varshney, A. K. Roy, Thermal properties of graphene: Fundamentals and appli-cations. MRS Bull. 37, 1273–1281 (2012). 46. D. Saha, S. Mahapatra, Analytical insight into the lattice thermal conductivity and heat capacity of monolayer MoS2. Phys. E 83, 455–460 (2016). Acknowledgments: We thank M. Kuiri and R. Su for help on device fabrication and J. Young and M. Na for help on measurements and modeling, respectively. Funding: We acknowledge support from the Natural Sciences and Engineering Research Council of Canada, Canada Foundation for Innovation, New Frontiers in Research Fund, Canada First Research Excellence Fund, and Max Planck–UBC–UTokyo Center for Quantum Materials. Z.Y. is also supported by the Canada Research Chairs Program. K.W. and T.T. acknowledge support from JSPS KAKENHI (grant numbers 19H05790, 20H00354, and 21H05233). Author contributions: J.W., D.Y., and Z.Y. conceived this work. D.Y., J.L., and J.W. fabricated the samples. J.W. and D.Y. conducted the measurement with the help from J.I.D., Y.X., and J.L. J.W., D.Y., and Z.Y. analyzed the data. M.W., E.O., and D.J. developed the optical parametric oscillator. K.W. and T.T. provided the hBN crystal. Z.Y. supervised the project. J.W., Z.Y., and J.I.D. wrote the manuscript based on the input from all other authors. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials.  Submitted 15 August 2022 Accepted 14 November 2022 Published 16 December 2022 10.1126/sciadv.ade3759 Wu et al., Sci. Adv. 8, eade3759 (2022) 16 December 2022                                                                                                                                                        7 of 7  S C I E N C E  A D VA N C E S | R E S E A R C H  A R T I C L E  Downloaded from https://www.science.org at National Institute for Materials Science on February 04, 2023Use of this article is subject to the Terms of serviceScience Advances (ISSN ) is published by the American Association for the Advancement of Science. 1200 New York Avenue NW,Washington, DC 20005. The title Science Advances is a registered trademark of AAAS.Copyright © 2022 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claimto original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).Ultrafast response of spontaneous photovoltaic effect in 3R-MoS2–basedheterostructuresJingda Wu, Dongyang Yang, Jing Liang, Max Werner, Evgeny Ostroumov, Yunhuan Xiao, Kenji Watanabe, TakashiTaniguchi, Jerry I. Dadap, David Jones, and Ziliang YeSci. Adv., 8 (50), eade3759. DOI: 10.1126/sciadv.ade3759View the article onlinehttps://www.science.org/doi/10.1126/sciadv.ade3759Permissionshttps://www.science.org/help/reprints-and-permissionsDownloaded from https://www.science.org at National Institute for Materials Science on February 04, 2023https://www.science.org/content/page/terms-service INTRODUCTION RESULTS AND DISCUSSION MATERIALS AND METHODS Sample fabrication Photocurrent measurement Supplementary Materials This PDF file includes: REFERENCES AND NOTES Acknowledgments