Abstract
Being able to control the neutral excitonic flux is a mandatory step for the development of future room-temperature two-dimensional excitonic devices. Semiconducting Monolayer Transition Metal Dichalcogenides (TMD-ML) with extremely robust and mobile excitons are highly attractive in this regard. However, generating an efficient and controlled exciton transport over long distances is a very challenging task. Here we demonstrate that an atomically sharp TMD-ML lateral heterostructure (MoSe2-WSe2) transforms the isotropic exciton diffusion into a unidirectional excitonic flow through the junction. Using tip-enhanced photoluminescence spectroscopy (TEPL) and a modified exciton transfer model, we show a discontinuity of the exciton density distribution on each side of the interface. We introduce the concept of exciton Kapitza resistance, by analogy with the interfacial thermal resistance referred to as Kapitza resistance. By comparing different heterostructures with or without top hexagonal boron nitride (hBN) layer, we deduce that the transport properties can be controlled, over distances far greater than the junction width, by the exciton density through near-field engineering and/or laser power density. This work provides a new approach for controlling the neutral exciton flow, which is key toward the conception of excitonic devices.
Similar content being viewed by others
Introduction
Electronics relies on the control of the motion of charge carriers to process information. The losses caused by the charged particle current and the resulting need to improve the power efficiency have fueled lots of interest in recent years1,2. As it is based on the control of electrically neutral quasi-particles (excitons), insensitive to long-range Coulomb scattering mechanisms, excitronics is by nature much more power efficient as it presents only negligible ohmic losses3,4. Nevertheless, developing excitronic devices is challenging, as it requires the ability to control the neutral exciton properties, such as the recombination rates, diffusion length (LD) and propagation direction, in an optically active medium at room temperature, without the help of external electric or magnetic fields3,5.
Owing to their promising optical properties, Transition Metal Dichalcognenide monolayers emerged as a highly versatile platform for excitonics system at the nanoscale6,7,8,9,10,11. In particular, their large exciton binding energy allows operating at room temperature. Due to their unique band properties, excitonic transport in TMD-MLs has led to the discovery of new fundamental phenomena such as a valley hall effect12,13,14, or the observation of nonlinear behavior such as a halo in the spatial profile15 or negative effective diffusion16. With the increasing maturity of the field, several basic components necessary to control exciton information have been developed, like room temperature excitonic transistor in a Van der Waals vertical heterostructure17, or an excitonic diode able to filter excitons in lateral heterostructures18,19,20. The physical mechanisms involved in the exciton transfer processes through lateral heterojunctions, as well as their influence on the exciton distribution, dynamics and the resulting photoluminescence observed in previous work18 are still unclear. This is mainly due to: (i) the extreme sharpness of lateral heterojunctions making complex its optical characterization, (ii) the strong influence of external parameters like defects, strain and dielectric environment on exciton transport, making the fabrication of high-quality samples with controlled dielectric environment crucial to extract the intrinsic junction properties, (iii) the lack of a complete theoretical description able to describe the asymmetric exciton transfer through the junction. Exciton diffusion at sharp interfaces have been previously studied in other materials like organic semiconductors in so-called bulk-heterojunction21,22. In the context of solar cells, efficient fast exciton diffusion and harvesting have stimulated strong theoretical and experimental efforts. Nevertheless, in such system the interface is designed to dissociate the exciton into free charges and not to transfer the entire quasi-particle through the junction21,22. Up to now, in TMD-ML the exciton diffusion is mainly driven either by strain23,24 or dielectric gradient engineering in TMD based vertical heterostructure19. Both approaches require complex architectures, with nanometric precision of the strain or the electrostatic potential over large distances (micrometers). The fabrication of those excitonic guides as well as their coupling with other circuit elements is thus very challenging.
In this paper, we investigate experimentally and theoretically the effect of an atomically sharp MoSe2-WSe2 lateral heterostructure (LH) on exciton diffusion and distribution. We performed tip-enhanced photoluminescence (TEPL) spectroscopy experiments, allowing sub-wavelength spatial resolution down to 30 nm, and developed an exciton transfer model. We show that the difference in the energy gap at the LH generates a discontinuity in the exciton density distribution analog to the temperature discontinuity found at interfaces presenting thermal resistance (Kapitza resistance). In steady state conditions, the presence of this discontinuity results in unique non-reciprocal exciton transport properties, taking place over distances far greater (two orders of magnitude) than the junction width and experimentally evidenced by: a highly asymmetric photoluminescence (PL) profile, the quenching of the WSe2-related PL and an enhancement of the MoSe2-related PL. Furthermore, by comparing the diffusion properties of fully hBN-encapsulated LH and hBN supported LH (without top hBN layer), we demonstrate that the diffusion properties of the LH can be tuned, by the generated exciton population either by increasing the laser power density or by modifying the optical near-field configuration (at constant laser power).
Experimental results
Sample preparation and characterization
The high quality monolayer MoSe2-WSe2 LH is grown using a modified CVD method described in Ref. 25. We then use water-assisted deterministic transfer to pick up the LH from the as-grown substrate and transfer it on a supporting flake of exfoliated hBN on SiO2/Si substrate. Finally, a second exfoliated hBN flake is transferred to cover the structure partially26. As a result, we obtain two distinct areas, as shown in Fig. 1a, a fully encapsulated area hBN/WSe2-MoSe2/hBN/SiO2/Si (e-LH), and an uncapped area WSe2-MoSe2/hBN/SiO2/Si (un-LH). The dashed white and yellow line in Fig. 1a highlight the boundaries of the bottom and top hBN flakes, respectively.
Figure 1c, f displays the room temperature μ-PL spectra measured on WSe2 and MoSe2, far from the junction on the encapsulated (e-LH) and uncapped (un-LH) regions, respectively. For the e-LH zone, the PL spectrum measured on MoSe2 ML exhibits a pronounced peak centered around 1.575 eV with a full width at half maximum (FWHM) of 50 meV that, in agreement with previous room temperature measurements27, can be assigned to the neutral exciton (A\({}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\)). In the case of WSe2, the PL spectrum is asymmetric, due to an intense peak linked to the A\({}_{1s}^{{{{{{\rm{WS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\) transition located around 1.65 eV (FWHMâ~â30 meV) and a second weaker peakâ~â40 meV below A1s, attributed to the spin-forbidden dark exciton (XD)27. For the un-LH region, all previously described PL features appear. Furthermore, we observe a slight increase of the broadening of all transitions, of the order of 5%. One can also notice that all PL spectra are more asymmetric. This could be attributed to an increased contribution of the charged excitons (trions), as the unprotected sample can get chemically charged when exposed to ambient air. In any case, one can notice that the integrated PL intensity \({I}_{{{{{{\rm{PL}}}}}}}^{\infty }\)(\({A}_{1s}^{{{{{{\rm{WS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\)) measured far from the interface is â3 times stronger than \({I}_{{{{{{\rm{PL}}}}}}}^{\infty }\)(\({A}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\)). Figure 1d, g is spectrally integrated μ-PL intensity color maps. The red (respectively blue) color intensity is the value of integrated PL (shaded areas in c and f) from 1.620 eV to 1.650 eV (respectively from 1.550 eV to 1.580 eV). Figure 1d (g) corresponds to the mapping of the encapsulated (uncapped) region. One can clearly see that the MoSe2 layer is organized in 6 folded stars surrounded by WSe2. The white regions, signature of low emission regions, have different origin. They correspond to bilayer inclusions when located at the center of MoSe2 stars, but also reveal cracks or/and inclusions that are related to the transfer process. We used the far-field PL color maps in combination with Raman mappings to select interfaces away from any visible defects.
Near-field studies of the lateral heterojunction
In this type of LH, the junction between the two materials WSe2 and MoSe2 is extremely sharp, down to ~3 nm as measured by electron microscopy18. Therefore, in order to gain insight into the transport properties around the junction, we use a sub-wavelength resolution tool, TEPL imaging and spectroscopy. Figure 2a displays a schematic of the experimental set-up, showing the linearly polarized laser excitation of 633 nm wavelength (â1.96eV) focused onto the apex of an atomic force microscope (AFM) silver-coated tip. The exciton generation profile is then linked to the electric field exaltation under the tip, a Gaussian profile with a full width at half maximum of ~30 nm, and a tip position (xtip) controlled with AFM resolution. As illustrated in the Fig. 2a excitons diffuse away from the excitation spot over long distances before recombining. The collection of emitted photons is ensured by a long working distance high numerical aperture (0.7 NA) microscope objective, with a fully open collection aperture. The position of the junction is determined very precisely (~30 nm) using tip-enhanced Raman spectroscopy (TERS), as described in ref. 18 and supplementary information (Fig. 6).
TEPL line scans with 30 nm step size are obtained by scanning the tip along a line perpendicularly to the LH junction. Figure 2b, d displays typical TEPL spectra measured for three tip positions (xtip) along the measured line for both encapsulated and uncapped samples. The spectra (1) and (3) recorded far from the interface (500 nm inside the MoSe2 region and 1.25 μm inside the WSe2 region) are similar to the ones described in Fig. 1. The spectra of WSe2 in locations (2) are recorded at 100 nm (respectively 300 nm) from the interface for the e-LH (resp. un-LH) sample. They show contributions of both MoSe2 and WSe2.
All TEPL spectra can be fitted using three Lorentzian peaks to account for the previously described relative contributions of A\({}_{1s}^{{{{{{\rm{WS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\) (dashed red line), X\({}_{{{{{{\rm{WS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}^{D}\) (dashed green line) and A\({}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\) (dashed blue line) as illustrated by the diagram in Fig. 2a. Figure 2c, e displays the result of this fitting procedure. The peak positions and FWHMs (shaded area) of the bright excitons A\({}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\) and A\({}_{1s}^{{{{{{\rm{WS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\) are shown as a function of the tip position xtip in the upper panels, while the integrated PL intensity (IPL) appears in the lower panels. For clarity, the contribution of X\({}_{{{{{{\rm{WS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}^{D}\) is not displayed (see supplementary information). First, we point out that the A\({}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\) signature (blue star) is clearly observed while the excitation is taking place inside WSe2, (xtipâ>â0, right side of the solid gray line in Fig. 2c, e). A contrario, no signature of A\({}_{1s}^{{{{{{\rm{WS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\) is observed when the excitation takes place in MoSe2 (xtipâ<â0, left side of the solid gray line). This reveals that nonreciprocal filtering is taking place at the interface, allowing the excitons to cross the junction from WSe2 to MoSe2. The other transport direction being forbidden by the junction, in agreement with18,28. The difference in the local dielectric environment between e-LH and un-LH has no impact on the diode-like effect. There is no visible influence on both the energy and FWHM of the PL spectra after crossing or being blocked at the junction. This suggests that there is no alteration of the nature of the excitons in each TMD-ML near the vicinity of the junction.
It is very interesting to notice that, in both e-LH and un-LH, the TEPL intensity originating from WSe2 (red stars) strongly decreases as the excitation takes place closer to the interface (xtipâââ0). As a consequence, the resulting signal measured near the LH interface is strongly dominated by A\({}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\) (see Fig. 2c, e bottom panels). In other words, even if the tip is located inside WSe2, most of the generated excitons migrate through the LH junction into MoSe2 before recombining. This behavior can only be explained by a strongly anisotropic transport, with the diffusion toward the junction becoming more efficient than the other directions. The second effect resulting from such an efficient exciton transfer A\({}_{1s}^{{{{{{\rm{WS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\,\)âââA\({}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\) is that as the tip approaches the interface from the WSe2 side the MoSe2 PL intensity IPL(A\({}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\)) increases and nearly reaches the values of the integrated PL intensity of WSe2 far from the interface (I\({}_{PL}^{\infty }\)(A\({}_{1s}^{{{{{{\rm{WS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\))). Far from the junction, MoSe2 is three times less bright than WSe2, meaning that the junction is enhancing the MoSe2 flake brightness at its vicinity.
Finally, a pronounced difference can be seen between the two systems, with A\({}_{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}^{1s}\) signature extending ~1.2 μm away from the interface in un-LH versus ~400 nm for e-LH. Both distances being considerably larger than the physical junction width (â3 nm18)
Modified exciton transfer model
In order to get a better understanding of the underlying physical phenomena taking place close to the interface, we developed an exciton transfer model, which is detailed in supplementary information. We theoretically investigate, in the low exciton density regime (typically 1011 cmâ2), the variation of neutral exciton density versus tip position, to be compared to the TEPL results. To do so, we analytically solved the linear 1D steady-state diffusion equation in each material i (1 = WSe2, 2 = MoSe2), given by :
n(x,âxtip) is the exciton density at x position with the excitation taking place at xtip, the x-axis origin being placed at the LH interface, Di the effective diffusion coefficient and Ïi the effective lifetime, radiative (\({\tau }_{i}^{r}\)) and non-radiative (\({\tau }_{i}^{nr}\)) (1/Ïi = 1/\({\tau }_{i}^{r}\)â+â1/\({\tau }_{i}^{nr}\)). \({L}_{i}=\sqrt{{D}_{i}{\tau }_{i}}\) represents the effective diffusion length in the material i. We use \({\Gamma }_{i}(x,\, {x}_{{{{{{\rm{tip}}}}}}})={\Gamma }_{{0}_{i}}(P,\, {\alpha }_{i},\, \nu ){{{{{{\rm{e}}}}}}}^{-\frac{{(x-{x}_{{{{{{\rm{tip}}}}}}})}^{2}}{{w}^{2}}}\) to simulate the exciton generation under the tip, centered at xtip. The width w corresponds to the tip diameter (~30 nm), P is the enhanced laser power under the tip and αi the absorption coefficient of the material i at the laser energy hν. We set boundary conditions to be n(xâââ±ââ,âxtip)â=â0. In agreement with experimental results, we impose the continuity of the integrated PL intensity for all values of xtip. Finally, we considered that bright excitons are not interacting with other types of excitons (dark excitons, B-excitons, momentum-forbidden excitons).
As illustrated in Fig. 3a, we model the junction as an ideally thin interface of fixed width (ϵ) much smaller than the excitonic diffusion length inside the barrier (âªLD), where no electron-hole pair recombination can occur. To model the asymmetric local effective drift of neutral exciton through the junction, we introduce a local uniform force field:
As a result, the junction imposes the continuity of the excitonic flux density on both sides of the interface by a local constant flux density \({\overrightarrow{j}}_{n}={\mu }_{b}\overrightarrow{F}n-{D}_{b}\overrightarrow{\nabla }n\), where μb and Db are the exciton mobility and diffusion coefficient inside the interface.
Figure 3b displays the calculated n(x,âxtip) for two different tip positions. We emphasize first on the lateral heterostructure specificity: the exciton density is discontinuous at the interface (see continuous lines). As a reference, the continuous exciton density calculated for classical diffusion (Fâ=â0) is display in dashed lines. When the tip excitation takes place inside MoSe2, the excitons are blocked by the junction, and the exciton density on the other side of the junction stays close to zero (see continuous green line). Excitons interacting with the junction loose a velocity v*â=âμbF blocking the diffusion process. This blocking phenomenon causes a redistribution of the exciton that strongly accumulate at the interface. On the other hand, when the excitation takes place inside WSe2, the excitons cross the junction, and are accelerated, hence acquiring an extra velocity v*. One can clearly see that this drives those excitons further inside the MoSe2 layer, while bringing the excitonic density down to zero near the interface on the WSe2 side (see continuous black line).
To facilitate further developments, and better characterize the exciton discontinuity at junction we define a partition coefficient linking the exciton densities n(xâ=â0+) and n(xâ=â0â) on both sides of the interface. Because of the asymmetrical behavior of the junction, this partition coefficient κâ=ân(xâ=â0+)/n(xâ=â0â) can be written, depending on the direction of excitonic flux, as follows:
with L1, L2 the diffusion lengths of both materials. κ quantifies the discontinuity of the exciton density at the junction, with κ+â=âκ(xtipâ>â0) describing the exciton diffusion from WSe2 to MoSe2 and κââ=âκ(xtipâ<â0) the diffusion from MoSe2 to WSe2. In our case, as the local equilibrium is established in the steady state, κ+â~âκâ (see Table 1), to simplify the discussion we will therefore refer simply to κ independently of the direction of excitonic flux. When κâ=â1, the exciton distribution, displayed as a dashed line in Fig. 3b, is continuous as observed in classical diffusion. When κââªâ1, corresponding to the present case, the exciton distribution is strongly asymmetric (see continuous line in Fig. 3b). Our finding can be interpreted as an extrapolation of the interfacial thermal resistance so-called Kapitza resistance, which describes the temperature discontinuity at atomically flat interface between two materials29,30. For high-quality LHs, the absolute value of the exciton Kapitza resistance can be defined as a function of the partition coefficient as follows: (See more details in supplementary information).
To compare directly the prediction of our near-field model with the experimental results, we calculate IPL(A\({}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\)) and IPL(A\({}_{1s}^{WS{e}_{2}}\)), the integrated PL intensities of each material as a function of the tip position. In the linear regime, the normalized PL intensity of MoSe2 and WSe2 can be written as (see supplementary information for details):
with \({I}_{{{{{{\rm{PL}}}}}}}^{\infty }({A}_{1s}^{{{{{{\rm{WS}}}}}}{{{{{{\rm{e}}}}}}}_{2}})\) used for normalization. The experimental PL is fitted using four fitting parameters: Aκ and Bκ, amplitude parameters, giving us access to κ, the partition coefficient and L1, L2. Figure 3c, d shows the experimental normalized PL intensities of MoSe2 (blue dots) and WSe2 (red dots) for both e-LH and un-LH configurations, respectively. The black solid lines represent the fits using Eqs. (5) and (6). The resulting fitting parameters are displayed in Table 1. In both systems, the κââªâ1 values for e-LH and un-LH confirm the strong asymmetry of the junction and the resulting strong Kapitza exciton resistance.
DISCUSSION
One can see that the model describes extremely well the strong quenching of the WSe2-related PL and the appearance of MoSe2-related PL when approaching the junction from the right (xtipâ>â0), shedding light on the experimental results on unidirectional excitonic transport across the junction. It is also able to quantitatively describe the enhancement of IPL(A\({}_{1s}^{{{{{{\rm{MoS}}}}}}{e}_{2}}\)) at the junction and its progressive decrease to the typically lower value \({I}_{PL}^{\infty }\)(A\({}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\)) of monolayer MoSe2 as the tip is moved away from the junction (xtipâ<â0). We would like to point out that this observation is far from trivial, it indicates that to conserve the experimentally observed continuity of the PL intensity, the exciton density at the interface (n(xâ=â0±)) is enhanced. We believe that this could be linked to the drastic change in the exciton velocity at the interface. Indeed, when the tip is inside MoSe2, excitons diffusing toward WSe2 are abruptly stopped at the junction, their average velocity going to zero. The non-radiative lifetime (\({\tau }_{2}^{nr}\)), being in part linked to the probability for the exciton to encounter non-radiative traps during its lifetime, is by consequence increased31, resulting in an increased population at the interface. This would explain why the size of PL enhancement area (where IPL(A\({}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\)) â¥â\({I}_{{{{{{\rm{PL}}}}}}}^{\infty }\)(A\({}_{1s}^{{{{{{\rm{MoS}}}}}}{{{{{{\rm{e}}}}}}}_{2}}\))) is a function of the diffusion length (see Fig. 3), as it only occurs when excitons start accumulating at the junction. Other phenomenon may contribute to this variation as a modification of the non-radiative Auger recombination due to the local enhancement of the exciton density.
Finally, the model reveals that the important difference observed in the PL profile between e-LH and un-LH is linked to drastically different values of the diffusion length, much shorter in e-LH than in un-LH (see experimental results in Fig. 2 and theoretical results in Table 1). To understand why \({L}_{D}^{e-LH} \, < \, {L}_{D}^{un-LH}\) we need to consider that even if the TEPL experiments on e-LH and un-LH were performed using the same laser power (400μW), the presence or absence of the top hBN layer, strongly impacts the near-field configuration under the tip. In first approximation, and considering a spherical tip of 30nm diameter, we can estimate the near-field intensity ratio on the TMD layer as:
ZeâLH and ZunâLH are the distances from the tip center to the TMD-ML in e-LH and un-LH areas, respectively. nBN and nair are hBN and air refractive index, respectively. Considering a 3 nm-thick hBN top layer (measured by AFM), we can estimate the enhancement ratio to be ~23. We can then expect the exciton density generated under the tip in un-LH to be more than one order of magnitude larger than in e-LH. To characterize the exciton transport properties of the WSe2 layer versus exciton density we use spatio-temporal PL signal in a confocal microscope. A local pulsed laser excitation (λ= 692 nm, pulse duration 1.5 ps), generates excitons, and we record the time evolution of the PL profile. We then extract the time evolution of the squared width w2(t)â=âw2(0)â+âÎw2(t) which is used to determine an effective diffusion coefficient in the 80 first picosecond where 90% of the PL signal originates : Îw2(t)â=â4Defft. An important point is to be able to decorrelate the effective lifetime Ïeff (linked to the PL decay) and the diffusion coefficients, from which we deduce the effective diffusion length \({L}_{{{{{{\rm{eff}}}}}}}=\sqrt{{D}_{{{{{{\rm{eff}}}}}}}{\tau }_{{{{{{\rm{eff}}}}}}}}\). Experimental conditions and detailed results are given in the supplementary information (Figs. 2â4). The excitation power density was varied over several order of magnitude and the results compared between the e-LH and un-LH samples. Figure 4 shows that, for both samples, Leff increases with the excitation power density (i.e. the excitonic density) which is mainly due to the increase of the effective diffusion coefficient (See supplementary information). This trend, attributed to a weak Auger contribution has been previously observed in WS231,32 and explain the large difference observed in the diffusion lengths we observed in e-LH and un-LH. The experimental results have been modeled by using a Auger coefficient of 0.210â6cm2/s (See supplementary information). This value is two orders of magnitude lower than the one estimated on non encapsulated WS2 (without top and bottom hBN) where a large difference of the transport properties values with a fully encapsulated layer was observed31. In summary, we show that both samples exhibit similar exciton transport properties, indicating the importance of the hBN bottom layer that prevents strong Auger scattering effect.
One can notice here that, according to Eq. (7), the enhancement ratio is dominated by the tip-TMD layer distance. This offers a unique opportunity, as using plasmonic tips or nanoantennas in combination with hBN encapsulation allows to control the diffusion length in this type of structures. As a matter of fact, controlling the thickness of the top hBN layer from a single layer to 20 layers, for example, would change the enhancement factor by a factor 10, modifying the diffusion length by the same amount, all the while avoiding the flaws, lower optical quality and exposition to environment, that are observed in the uncapped LH.
In summary, we have performed a detailed tip-enhanced spectroscopy study of a MoSe2-WSe2 lateral heterostructure, and have developed a model to render the observation of the unilateral transport across the junction observed in near-field PL experiment. It accounts for a discontinuity of the exciton density at the interface. We have thus shown that the exciton diffusion properties follow a semi-classical process, due to the difference in the energy gap between the two materials: the usual isotropic in-plane diffusion of the excitons is frustrated by the junction, leading to an asymmetric diffusion, in which all generated excitons move away from the high bandgap WSe2 layer to recombine in the low bandgap MoSe2 layer. This transfer causes near the interface the quenching of WSe2-related PL and the enhancement of the MoSe2-related PL, well above the one observed in âbulkâ MoSe2 far from the junction. Interestingly, we observe similar asymmetric diffusion property for samples with or without top hBN. Both samples present similar intrinsic transport properties, probably linked to an efficient screening of the dielectric disorder by the bottom encapsulation by hBN layer. Finally, we have shown that the diffusion length in WSe2 is strongly dependent on the exciton density. This offers a new degree of freedom, as changing the laser power density or the near-field enhancement (for instance using optically resonant nanoantennas instead of the plasmonic tip) would allow tuning the diffusion length to any wanted value from tens of nanometer up to few micrometers. TMD-based lateral heterostructure is a rapidly evolving research field that could offer to combine TMDs with very different bandgaps, thus allowing partition zone engineering, or creating more complex designs using three or more TMDs (lateral excitonic quantum wells). This work offers both the theoretical and experimental tools to predict and control the new diffusion properties that will be at the origin of new excitronic devices.
Methods
We use water-assisted deterministic transfer to pick up the chemical vapor deposition LH from the as-grown substrate using polydimethylsiloxane (PDMS) and transfer it on a supporting flake of exfoliated hBN on SiO2/Si substrate26. TERS and TEPL are carried out with state-of-the-art commercial system (Trios OmegaScope-R coupled with LabRAM spectrometer, Horiba Scientific). Silver-coated tips with an apex radius of 15nm were used for tip-enhanced measurements. The Time evolution of the PL profile experiment is based on a diffraction-limited laser excitation that induces lateral diffusion of the photogenerated excitonic species. We used a Streak camera system to record the time evolution of the PL spatial profile IPLâ=âI(x,ât) with a time resolution of 5.5 ps. The Ti:Sa laser excitation is set to Eexâ=â1.79 eV, with a 80 MHz repetition frequency, 1.5 ps pulse width and we vary the excitation power from 10 μW to 1 mW.
Data availability
The data for all figures in the main text are available in the Source Data file. The data presented in Supplementary information are available upon request due to their large file size. Source data are provided with this paper.
References
Butov, L. Excitonic devices. Superlattices Microstruct. 108, 2 (2017).
Miller, D. Rationale and challenges for optical interconnects to electronic chips. Proc. IEEE 88, 728 (2000).
High, A. A., Novitskaya, E. E., Butov, L. V., Hanson, M. & Gossard, A. C. Control of exciton fluxes in an excitonic integrated circuit. Science 321, 229 (2008).
Peng, R. et al. Long-range transport of 2D excitons with acoustic waves. Nat. Commun. 13, 1334 (2022).
High, A. A., Hammack, A. T., Butov, L. V., Hanson, M. & Gossard, A. C. Exciton optoelectronic transistor. Opt. Lett. 32, 2466 (2007).
Wang, G. et al. Colloquium : Excitons in atomically thin transition metal dichalcogenides. Rev. Mod. Phys. 90, 021001 (2018).
Wang, J., Verzhbitskiy, I. & Eda, G. Electroluminescent devices based on 2D semiconducting transition metal dichalcogenides. Adv. Mater. 30, 1802687 (2018).
Novoselov, K. S., Mishchenko, A., Carvalho, A. & Castro Neto, A. H. 2D materials and van der Waals heterostructures. Science 353, aac9439 (2016).
Cheiwchanchamnangij, T. & Lambrecht, W. R. L. Quasiparticle band structure calculation of monolayer, bilayer, and bulk MoS 2. Phys. Rev. B 85, 205302 (2012).
Yuan, L., Wang, T., Zhu, T., Zhou, M. & Huang, L. Exciton dynamics, transport, and annihilation in atomically thin two-dimensional semiconductors. J. Phys. Chem. Lett. 8, 3371 (2017).
Rosati, R. et al. Interface engineering of charge-transfer excitons in 2D lateral heterostructures. Nat. Commun. 14, 2438 (2023).
Mak, K. F., McGill, K. L., Park, J. & McEuen, P. L. The valley Hall effect in MoS2 transistors. Science 344, 1489 (2014).
Ubrig, N. et al. Microscopic origin of the valley Hall effect in transition metal dichalcogenides revealed by wavelength-dependent mapping. Nano Lett. 17, 5719 (2017).
Zhou, B. T., Taguchi, K., Kawaguchi, Y., Tanaka, Y. & Law, K. T. Spin-orbit coupling induced valley Hall effects in transition-metal dichalcogenides. Commun. Phys. 2, 26 (2019).
Kulig, M. et al. Exciton diffusion and halo effects in monolayer semiconductors. Phys. Rev. Lett. 120, 207401 (2018).
Rosati, R., Perea-CausÃn, R., Brem, S. & Malic, E. Negative effective excitonic diffusion in monolayer transition metal dichalcogenides. Nanoscale 12, 356 (2020).
Unuchek, D. et al. Room-temperature electrical control of exciton flux in a van der Waals heterostructure. Nature 560, 340 (2018).
Beret, D. et al. Exciton spectroscopy and unidirectional transport in MoSe2-WSe2 lateral heterostructures encapsulated in hexagonal boron nitride. npj 2D Mater. Appl. 6, 84 (2022).
Shanks, D. N. et al. Interlayer exciton diode and transistor. Nano Lett. 22, 6599 (2022).
Bellus, M. Z. et al. Photocarrier transfer across monolayer MoS2-MoSe2 lateral heterojunctions. ACS Nano 12, 7086 (2018).
Firdaus, Y. et al. Long-range exciton diffusion in molecular non-fullerene acceptors. Nat. Commun. 11, 5220 (2020).
Siegmund, B. et al. Exciton diffusion length and charge extraction yield in organic bilayer solar cells. Adv. Mater. 29, 1604424 (2017).
Dirnberger, F. et al. Quasi-1D exciton channels in strain-engineered 2D materials. Sci. Adv. 7, eabj3066 (2021).
Cordovilla Leon, D. F., Li, Z., Jang, S. W., Cheng, C.-H. & Deotare, P. B. Exciton transport in strained monolayer WSe2. Appl. Phys. Lett. 113, 252101 (2018).
Najafidehaghani, E. et al. 1D p-n junction electronic and optoelectronic devices from transition metal dichalcogenide lateral heterostructures grown by one-pot chemical vapor deposition synthesis. Adv. Funct. Mater. 31, 2101086 (2021).
Paradisanos, I. et al. Controlling interlayer excitons in MoS2 layers grown by chemical vapor deposition. Nat. Commun. 11, 2391 (2020).
Poumirol, J.-M. et al. Unveiling the optical emission channels of monolayer semiconductors coupled to silicon nanoantennas. ACS Photonics 7, 3106 (2020).
Shimasaki, M. et al. Directional exciton-energy transport in a lateral heteromonolayer of WSe2-MoSe2. ACS Nano 16, 8205 (2022).
Alosious, S., Kannam, S. K., Sathian, S. P. & Todd, B. D. Kapitza resistance at water-graphene interfaces. J. Chem. Phys. 152, 224703 (2020).
Pollack, G. L. Kapitza resistance. Rev. Mod. Phys. 41, 48 (1969).
Zipfel, J. et al. Exciton diffusion in monolayer semiconductors with suppressed disorder. Phys. Rev. B 101, 115430 (2020).
Uddin, S. Z. et al. Enhanced neutral exciton diffusion in monolayer WS2 by exciton-exciton annihilation. ACS Nano 16, 8005 (2022).
Acknowledgements
Toulouse acknowledges partial funding from ANR IXTASE (ANR-20-CE30-0032), ANR HiLight (ANR-19-CE24-0020-01), ANR Ti-P (ANR-21-CE30-0042), NanoX project 2DLight (ANR-17-EURE-0009), the Institut Universitaire de France, and the EUR grant ATRAP-2D NanoX ANR-17-EURE-0009 in the framework of the âProgramme des Investissements dâAvenirâ, the Institute of quantum technology in Occitanie IQO and a UPS excellence PhD grant. Growth of hexagonal boron nitride crystals was supported by JSPS KAKENHI (Grants No. 19H05790, No. 20H00354 and No. 21H05233). The Jena group received financial support of the Deutsche Forschungsgemeinschaft (DFG) through a research infrastructure grant INST 275/257-1 FUGG, CRC 1375 NOA (Project B2), SPP2244 (Project TU149/13-1) as well as DFG grant TU149/16-1. This project has also received funding from the joint European Unionâs Horizon 2020 and DFG research and innovation program FLAG-ERA under grant TU149/9-1.
Author information
Authors and Affiliations
Contributions
Z.G., A.G., and A.T. developed the growth method and fabricated the lateral heterostructures. T.T. and K.W. grew the hBN crystals. I.P. carried out hBN encapsulation. D.B., L.L., P.R, D.L, and X.M. developed, performed and analyzed effective diffusion coefficient and exciton lifetime experiments. H.L. and J.-M.P. performed TEPL and TERS. V.P. and J.-M.P. optimized the TEPL and TERS set-up. H.L., J.-M.P., N.C, and V.P developed the analytical modeling. J-M.P., V.P., L.L., B.U. wrote the manuscript with inputs from all the authors and supervised the project.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewers 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
Source data
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
Lamsaadi, H., Beret, D., Paradisanos, I. et al. Kapitza-resistance-like exciton dynamics in atomically flat MoSe2-WSe2 lateral heterojunction. Nat Commun 14, 5881 (2023). https://doi.org/10.1038/s41467-023-41538-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-023-41538-6
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.