Research ArticleCONDENSED MATTER PHYSICS

Moiré potential impedes interlayer exciton diffusion in van der Waals heterostructures

See allHide authors and affiliations

Science Advances  23 Sep 2020:
Vol. 6, no. 39, eaba8866
DOI: 10.1126/sciadv.aba8866

Abstract

The properties of van der Waals heterostructures are drastically altered by a tunable moiré superlattice arising from periodically varying atomic alignment between the layers. Exciton diffusion represents an important channel of energy transport in transition metal dichalcogenides (TMDs). While early studies performed on TMD heterobilayers suggested that carriers and excitons exhibit long diffusion, a rich variety of scenarios can exist. In a moiré crystal with a large supercell and deep potential, interlayer excitons may be completely localized. As the moiré period reduces at a larger twist angle, excitons can tunnel between supercells and diffuse over a longer lifetime. The diffusion should be the longest in commensurate heterostructures where the moiré superlattice is completely absent. Here, we experimentally demonstrate the rich phenomena of interlayer exciton diffusion in WSe2/MoSe2 heterostructures by comparing several samples prepared with chemical vapor deposition and mechanical stacking with accurately controlled twist angles.

INTRODUCTION

van der Waals (vdW) materials provide exciting opportunities to create new heterostructures with vastly expanded choices of materials, as the strict requirement of lattice matching between adjacent layers has been lifted (1). Most intriguingly, a moiré superlattice emerges because of the periodic variations in atomic alignment between adjacent layers. The period of the moiré superlattice can be readily controlled by choosing materials with different lattice constants or adjusting the twist angle between the layers, leading to a new paradigm in engineering quantum materials (24). For example, graphene bilayers can exhibit either superconducting or insulating phases driven by electron correlation controllable via the twist angle (58). Such a strategy has yet to be carefully explored in controlling semiconductor heterostructures built with transition metal dichalcogenides (TMDs) (9, 10).

A type II band alignment is typically found in a TMD heterobilayer (hBL), leading to rapid electron transfer to one layer and hole transfer to the other layer and the formation of interlayer excitons (IXs) (1113). These IX resonances are the lowest energy optical excitations and shifted from the intralayer excitons mainly by the band alignments (1419). The diffusion and mobility of various carriers including electrons, holes, and IXs are important properties because they determine the charge and energy transport processes in hBLs, which, in turn, are critical for building transistors and photovoltaic devices. Several experiments performed on WS2/WSe2, MoS2/WSe2, and MoSe2/WSe2 hBLs have reported long free carrier and IX diffusion length and electric field control of this transport (16, 17, 2022). Other recent experiments, on the other hand, suggest that IXs may be localized by the moiré potential on the order of 100 to 200 meV (18, 19, 23). These seemingly conflicting results suggest rich opportunities to control the optical properties of TMD heterostructures.

Here, we investigate the influence of the moiré potential on IX diffusion in WSe2/MoSe2 hBLs using spatially and time-resolved photoluminescence (PL) measurements. We carefully compare experiments performed on two types of samples either grown by chemical vapor deposition (CVD) or created by the mechanical exfoliation and transfer (MET) with a well-controlled twist angle. IX diffusion length in the rotationally aligned CVD-grown hBL exceeds the size of the hBL and is estimated to be a few microns. In contrast, no IX diffusion beyond the excitation laser spot is observed in the stacked hBL with ~1° twist angle, suggesting that IXs may be localized by the moiré potential. In another MET sample with ~3.5° twist angle, IX diffusion is observed but with a short diffusion length of ~1 μm despite nearly ~100-fold longer lifetime than that in the CVD-grown sample. These experiments suggest that IX diffusion is highly controllable via engineering the moiré superlattice and provides valuable guidance to extensive effort in searching for IX Bose-Einstein condensate (BEC) (24, 25) or localized quantum emitters in the TMD heterostructures (3, 4, 18, 19, 26, 27).

RESULTS

We briefly describe the samples investigated here. We focus on WSe2/MoSe2 hBLs with R-stacking style or near this stacking style, where a deep moiré potential is expected (28, 29). The stacking style is confirmed by transmission electron microscope (TEM) images or second harmonic generation (SHG) (see the Supplementary Materials). The optical microscope images and structural illustrations for the CVD (sample A) and MET (samples B and C) methods are shown in Fig. 1 (A, B, D, and E, respectively). The lattice constants of MoSe2 and WSe2 are very close to each other. In samples A-1 and A-2 grown by the CVD method, the lattice in each layer is slightly distorted to form a commensurate and rotationally aligned (i.e., 0° twist angle) hBL as illustrated in Fig. 1C. The hBL regions are limited to the bright inner triangle with a typical lateral size of ~10 μm (Fig. 1A). In the MET samples, a twist angle θ = 1.1° ± 0.3° (sample B-1), θ = 1.3° ± 0.3° (sample B-2), and θ = 3.5° ± 0.3° (sample C) from the R-stacking leads to the formation of a moiré superlattice (Fig. 1F). The twist angle is determined from a high-resolution optical image, SHG (see the Supplementary Materials), and further confirmed by IX lifetime measurements presented later. The presence of the moiré potential and its periodicity plays a major role in determining exciton diffusion as illustrated in Fig. 1 (C and F) and demonstrated in the rest of the paper. Both intralayer excitons and IXs are observed (PL spectra included in the Supplementary Materials), while we focus on the properties of the IXs in this paper.

Fig. 1 Two types of WSe2/MoSe2 heterostructures and the absence/presence of moiré superlattice.

(A) Optical images of CVD-grown sample A-1. (B) Illustrated CVD-grown hBL, which consists of a monolayer WSe2 covering on top of a monolayer MoSe2-WSe2 lateral heterostructure. (C) Schematics of exciton diffusion in a CVD-grown WSe2/MoSe2 bilayer with a commensurate structure. (D) Optical image of mechanically stacked WSe2/MoSe2 hBL sample B-1. (E) Illustration of the hBN-encapsulated WSe2/MoSe2 hBL. (F) impeded exciton diffusion by the moiré potential in a MET hBL. The moiré supercell size decreases with an increasing twist angle θ between the two layers.

IX diffusion is investigated using spatially and spectrally resolved PL in all samples as shown in Fig. 2 (A to F). We examine the spatial PL line profiles at a few selected energies at the center, lower, and higher energy side of the broad IX spectral feature. No pronounced energy dependence is observed in any of these samples. However, the spatial IX diffusion is remarkably different between these samples. The diffusion lengths in the CVD samples A-1 and A-2 are the longest and limited by the physical size of the hBL region. In contrast, no IX diffusion beyond the excitation laser spot size is observed in the MET samples B-1 and B-2 with θ = ~1°. The full width at half maximum (FWHM) of a Gaussian function fitted to the excitation laser spot yields ~1 μm. The size of the IX PL spot appears slightly larger than the excitation laser spot size in Fig. 2E, which is caused by the imaging optics (e.g., microscope objective and curved mirrors inside the spectrograph) at different wavelengths (i.e., IXs near 900 nm and the excitation laser at 660 nm). Interestingly, the IX diffusion is found to depend on the moiré period and the IX lifetimes. When the moiré period is reduced from 20 nm in sample B-1 to 5.7 nm in another MET sample C with θ = ~3.5°, IX diffusion beyond the excitation spot size becomes observable as shown in Fig. 2 (C and F).

Fig. 2 Spatially resolved PL images.

Two-dimensional PL images and spectra of the IXs with the spatial (energy) coordinate were shown in the horizontal (vertical) axis for (A) CVD-grown sample A-1, (B) MET sample B-1 with θ = 1.1° ± 0.3°, and (C) MET sample C with θ = 3.5° ± 0.3°. A much longer IX diffusion length in sample A-1 is clearly observed. The corresponding PL line profiles were taken from (D) sample A-1, (E) sample B-1, and (F) sample C, respectively. Note that the PL profiles from sample A-1 are truncated by the boundaries of the hBL region as indicated by the white dashed lines in (A). In addition, line profiles of excitation laser spots (bottom gray lines) are displayed. The 900-nm laser profile (top gray shaded) is displayed in (E) to illustrate that the PL image is limited by excitation spot size and imaging optics in sample B-1. Measured data points are fitted by a Gaussian function. a.u., arbitrary units.

IX diffusion can be quantitatively analyzed using the following equation (30, 31)n(r,t)t=·(DIXn(r,t)n(r,t)μu0n(r,t))n(r,t)τIX+G(r,t)(1)where n(r, t) is the IX density, DIX is the diffusion coefficient, μ is the exciton mobility, u0 is the IX interaction energy, τIX is the IX lifetime, and G(r, t) is the exciton generation term, respectively (more details in the Supplementary Materials). To avoid the complications arising from exciton-exciton annihilation and repulsive interaction at high excitation density (20, 32, 33), all experiments presented here are performed at a low excitation power of 1 μW, corresponding to an estimated IX density of ~109/cm2. We present the power-dependent studies for all three samples in the Supplementary Materials. Making the assumption that the interaction effect can be neglected at low density, the steady-state IX density follows the analytical solution (32)n(r)K0(rLIX)e(rr')2/w2dr'(2)

This solution is a convolution between the Gaussian profile of incident laser and the modified Bessel function of the second kind K0 with the diffusion length LIX=DIXτIX, where τIX and DIX are the IX lifetime and diffusion coefficient, respectively. A series of solutions are generated by varying LIX as shown in Fig. 3A. Comparing the 1/e width of the simulated curves and the experimental data, we extracted a range of diffusion lengths for CVD sample A-1 (A-2) to be 3 to 16 (2.4 to 13.5) μm and for sample C to be ~1 μm. The uncertainty in diffusion lengths in CVD-grown samples partially arises from the truncated data due to the limited size of the flake. We used several different models to analyze LIX and the range of values as shown in Fig. 3B (see the Supplementary Materials for more details).

Fig. 3 Determination of IXs diffusion length (LIX) and measured IX lifetimes.

(A) Calculated spatial distribution of IX density, in which n(r) is convoluted by the Gaussian laser profile with ~1-μm FWHM spot size. The measurements are compared to the simulated density distribution function to obtain the diffusion length for each sample. (B) IX diffusion length LIX from samples A-1, A-2, and C. (C) Time-resolved PL (TRPL) spectra from samples A-1, B-1, and C, respectively. The TRPL spectrum was taken at the center of the IX resonances. The TRPL spectrum is fitted by a single- or double-exponential decay function. The lifetimes extracted from fittings are summarized in (D) for all five samples. Note that the error bars of lifetimes are smaller than the dots.

We further analyze the diffusion coefficients of IXs by measuring IX lifetimes using time-resolved PL (TRPL) measurements. The data in Fig. 3C are fitted with either a biexponential or single-exponential decay function, and extracted lifetimes are shown in Fig. 3D. In sample A-1 (A-2), a biexponential fit leads to a fast and slow component of τ1A = 1.2 (1.0) ns and τ2A = 5.5 (5.9) ns and a weighted average lifetime of 1.7 (1.8) ns. This IX lifetime is similar to τB = ~1 ns for samples B-1 and B-2 but markedly different from τC = ~100 ns for sample C. The lifetimes measured from both types of hBLs are consistent with the previous experiments (14, 19, 34, 35). The remarkably longer lifetime found in sample C is related to the momentum space shift between the conduction and valence bands in twisted bilayers, leading to an indirect transition in both real and momentum space (36). This angle-dependent lifetime is confirmed in many of our hBLs prepared by the MET method. Using DIX = LIX2IX, an effective diffusion coefficient is estimated to be ~50 (30) cm2/s for sample A-1 (A-2), larger than the typical value reported for incommensurate hBLs (16, 20, 21). The extracted value for sample C is ~0.1 cm2/s, consistent with previous studies (see the Supplementary Materials for details).

The absence of the moiré pattern in sample A is directly proven by a high-resolution scanning TEM (STEM) image as shown in Fig. 4A. The brightest points identify the positions of Mo atoms, the second brightest points correspond to W atoms, and the dimmest ones are from the Se atoms. The commensurate atomic registry is proven by the uniform pattern observed in the STEM image. The schematic in Fig. 4B illustrates the 3R-like stacking of this atomic registry, in which the Se atoms are directly positioned on top of the Mo atoms. A MET hBL cannot be easily transferred and identified on a standard TEM grid. Thus, we perform atomic force microscopy (AFM) measurements to examine the atomic alignment in a mechanically stacked hBL. We follow the same procedure to prepare a sample without the top hexagonal boron nitride (hBN) to expose the WSe2/MoSe2 hBL, which enables AFM imaging (see Materials and Methods for details). Despite some contaminants on the sample, we observe a disordered hexagonal pattern with a period of ~20 nm, suggesting a small twist angle of ~1° in this particular sample (Fig. 4C).

Fig. 4 Crystal structures of CVD-grown and MET hBLs.

(A) A STEM image of a CVD-grown hBL similar to sample A shows a commensurate atomic registry, where the Bragg-filtered image is shown in the inset for comparison. Scale bar, 1 nm. (B) Schematic showing the ideal atomic registry with top (side) view in the top (bottom) panel. The STEM contrast originates from the difference of total atomic number in each atomic column. The different atomic columns are indicated by black dashed lines: column 1, one Mo atom and two Se atoms (blue triangle); column 2, one W atom (red triangle); column 3, two Se atoms (orange triangle). (C) Filtered AFM phase image of an hBL prepared similar to samples B and C, showing a hexagonal pattern consistent with the moiré superlattice. The bright spots at the upper left corner are from surface contaminations.

DISCUSSION

Recognizing that the quality of hBLs can vary, we repeated these measurements on different samples and on different locations of the same sample (more data in the Supplementary Materials). All experiments yielded consistent observations. The longest diffusion length is always observed in the CVD-grown commensurate WSe2/MoSe2 hBLs. A finite diffusion length of ~1 μm is observed in a stacked hBL with a relatively large twist angle, while no diffusion beyond the laser spot is observed in stacked hBLs with a small twist angle (θ = ~1°).

The strikingly different IX diffusion in CVD and MET hBLs is most naturally explained by taking into account the absence/presence of the moiré potential. While there are other differences between CVD-grown and MET hBLs—e.g., substrates, defect density, and interface quality—these differences cannot account for the different IX diffusions and the controllable changes on the stacked samples. The CVD samples are not grown on an atomically smooth substrate and likely exhibit higher defect density. Both factors should lead to a shorter diffusion length than the MET samples (33). Furthermore, the notably longer IX lifetime in the stacked hBL with θ = ~3.5° (sample C) should lead to a longer diffusion length than the CVD sample A if their diffusion coefficients were similar.

The difference between two MET hBLs arises from the expected changes in moiré periods and measured exciton lifetimes. For a moiré crystal with a small periodicity, excitons may tunnel between different supercells, and diffusion is observable over a long IX lifetime. For a moiré crystal with a relatively large periodicity, the tunneling between different supercells is exponentially reduced. IXs become localized within a moiré supercell due to a reduced lifetime and decreased velocity. Our conclusion is qualitatively consistent with previous experiments (16, 17, 20, 21) that reported finite IX or carrier diffusion in hBLs with a small moiré period in the sense exciton diffusion beyond the diffraction limit can be observed. Although some quantitative variations in diffusion lengths and diffusion coefficients exist in the literature, they mostly arise from the different diffusion models and analysis (see the Supplementary Materials) and, to a less extent, the preparation procedure of MET bilayers. To establish a more detailed correlation between the moiré supercell size and exciton diffusion, future experimental techniques with higher spatial and/or temporal resolutions are required.

Moiré potential has been demonstrated to have a profound influence on the electronic and optical properties of vdW heterostructures. However, its role in exciton diffusion has been largely ignored in previous studies on hBLs consisting of TMD monolayers with rather different lattice constants or with larger twist angles (16, 17, 20). Our experiments present a complementary view of IX diffusion from previous studies and highlight the influence of the moiré potential on exciton diffusion in hBLs. These studies provide critical guidance for intensive efforts on searching for IX BEC or quantum emitters in TMD heterostructures. IX BEC only occurs when the exciton density exceeds a threshold value. In this case, mobile excitons are preferred and are more likely found in either CVD-grown samples or MET samples with small moiré periods. On the other hand, localized excitons leading to a regular array of quantum emitters are more likely hosted in stacked hBLs with a small twist angle and large moiré period.

MATERIALS AND METHODS

Optical measurements

For the CVD-grown sample A, PL measurements were taken at T = 4 K inside a cryogen-free low-vibration cryostat equipped with a three-axis piezo-positioner and an objective lens (numerical aperture, 0.82) mounted inside the cryostat. The excitation sources at 632 nm (750 nm) were provided by a He-Ne (continuous wave Ti:Sapphire) laser. The PL signal was then sent to a 0.75-m monochromator and detected by a nitrogen-cooled charge-coupled device camera. For TRPL measurements, a 635-nm pulsed laser diode (~90-ps width and 40-MHz repetition rate) was used as the excitation source, and the PL was detected by a fast avalanche diode. The PL decays were recorded using the time-correlated single-photon counting technique with a temporal resolution of ~100 ps. Slightly different experimental conditions for measurements performed on samples B and C can be found in the Supplementary Materials.

CVD growth process and STEM characterizations

Single-crystal WSe2/MoSe2 hBLs were grown on sapphire substrates by CVD following a previously developed one-step growth process (37, 38). High-purity MoO2 (99%; Sigma-Aldrich), WO3, and Se powders (99.5%; Alfa Aesar) were used as the initial reactants. The hBLs were grown at 880°C in H2/Ar flowing gas at low pressure (5 to 40 torr). During the growth, the flow rates for H2/Ar gas were set to 6/60 standard cubic centimeters per minute. For STEM measurements, the samples were capped with poly(methylmethacrylate) (PMMA) (950K A4) by spin-coating and then baked at 100°C for 60 min. The WSe2/MoSe2 hBLs were then transferred onto a TEM grid by PMMA-based wet transfer procedures. Annular dark-field STEM images were obtained using a spherical aberration-corrected TEM (JEOL-2100F). The detailed growth and characterization processes can be found in (29).

AFM measurements

We obtained the phase image using the tapping mode of AFM (NX10, Park Systems) with a conductive cantilever (PPP-EFM) at ~28°C and ~68% relative humidity. To obtain the image, the tip was brought to the sample as close as possible, i.e., if the setpoint was lowered further, then a large oscillating noise appeared due to the tip-sample interaction. The hexagonal pattern was only observable in selected areas of the sample due to surface contaminants. The inverse fast Fourier transform of an AFM phase image was obtained using XEI (Park Systems) software.

Mechanically stacked bilayer preparation

The MoSe2, WSe2 monolayers (2D semiconductors), and the top hBN capping layer were mechanically exfoliated onto a polydimethylsiloxane (PDMS) sheet. The bottom hBN layer was directly exfoliated onto a SiO2/Si substrate. The crystal orientations of MoSe2 and WSe2 monolayers were identified by polarization-resolved SHG. Using a hemisphere PDMS dot with polycarbonate (PC) film–based dry-transfer technique (39), the top hBN layer was picked up and used to stack the MoSe2 and WSe2 bilayer by picking up each layer sequentially. The rotational alignment of the crystal orientations of the monolayers was guided by high-resolution optical microscope images. The stacked hBN/MoSe2/WSe2 layers were transferred on the bottom hBN layer, and then the heterostructure was dipped in an acetone bath to clean off PC residue. Thermal annealing in a high vacuum (~10−7 mbar) at 200°C for 4 hours was performed to improve the interfaces. Polarization-resolved SHG measurements on the heterostructure region were performed to confirm the twist angle and distinguish between R- and H-stacking orders (40).

SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/39/eaba8866/DC1

https://creativecommons.org/licenses/by-nc/4.0/

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.

REFERENCES AND NOTES

Acknowledgments: Funding: The spectroscopic experiments performed by J.C. at UT-Austin and X. Li were supported by the Department of Energy, Basic Energy Science program via grant DE-SC0019398. Partial support for K.T. was provided by the NSF MRSEC program DMR-1720595, which also facilitated the collaboration between the group of C.-K.S. and X.L. L.S. and C.-K.S. were supported by the Welch Foundation via grant F-1662 and F-1672. M.S. and K.J. were supported by NSF EFMA-1542747. C.-K.S. acknowledges support from the U.S. Air Force via grant FA2386-18-1-4097. W.-H.C. acknowledges the support from the Ministry of Science and Technology (MOST) of Taiwan (105-2119-M-009-014-MY3 and 107-2112-M-009-024-MY3). W.-T.H. acknowledges the support from the MOST of Taiwan (MOST-107-2917-I-564-010). M.-H.L. and M.-W.C. acknowledge the support from MOST of Taiwan. S.G. and C.-Y.W. acknowledge the support from the MOST of Taiwan (MOST 108-2119-M-007-008). J.Q. acknowledges the support from the China Scholarship Council (grant no. 201706050068). K.W. and T.T. acknowledge support from the Elemental Strategy Initiative conducted by the MEXT, Japan and the CREST (JPMJCR15F3), JST. S.K. was financially supported by the National Research Foundation (NRF) of Korea grant funded by the Korea Government (2017R1D1B04036381). The collaboration between National Tsing-Hua University and The University of Texas at Austin is facilitated by the Global Networking Talent (NT 3.0) Program, Ministry of Education in Taiwan. Author contributions: J.C. and W.-T.H. led the optical experiments. L.S., H.-Y.C., and C.-Y.W. assisted the experiments. J.C. and L.-S.L. led the sample preparations. K.T., M.S., and K.J. assisted in making the samples. T.T. and K.W. provided hBN samples. M.-H.L., W.-T.H., and L.-S.L. performed STEM measurements and analysis. J.Q. and S.K. performed AFM measurements and analysis. J.C., W.-T.H., X.L., and W.-H.C. wrote the manuscript. X.L., M.-W.C., S.G., S.K., C.-K.S., and W.-H.C. supervised the project. All authors discussed the results and commented on the manuscript at all stages. 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. Additional data related to this paper may be requested from the authors.
View Abstract

Stay Connected to Science Advances

Navigate This Article