Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Interlayer exciton mediated second harmonic generation in bilayer MoS2

## Abstract

Second-harmonic generation (SHG) is a non-linear optical process, where two photons coherently combine into one photon of twice their energy. Efficient SHG occurs for crystals with broken inversion symmetry, such as transition metal dichalcogenide monolayers. Here we show tuning of non-linear optical processes in an inversion symmetric crystal. This tunability is based on the unique properties of bilayer MoS2, that shows strong optical oscillator strength for the intra- but also interlayer exciton resonances. As we tune the SHG signal onto these resonances by varying the laser energy, the SHG amplitude is enhanced by several orders of magnitude. In the resonant case the bilayer SHG signal reaches amplitudes comparable to the off-resonant signal from a monolayer. In applied electric fields the interlayer exciton energies can be tuned due to their in-built electric dipole via the Stark effect. As a result the interlayer exciton degeneracy is lifted and the bilayer SHG response is further enhanced by an additional two orders of magnitude, well reproduced by our model calculations. Since interlayer exciton transitions are highly tunable also by choosing twist angle and material combination our results open up new approaches for designing the SHG response of layered materials.

## Introduction

Nonlinear optical spectroscopy of atomically thin crystals gives crucial insights into their physical properties and investigates potential applications1,2. Nonlinear optics is based on the interaction of photons within the crystal. This can give rise, among other processes, to second-harmonic generation (SHG) where two photons of identical energy combine into one photon of twice this energy3,4. Currently existing applications in photonics and laser physics make use of this process in bulk crystals for example using potassium dihydrogen phosphate (KDP) or lithium niobate (LiNbO3)5.

SHG is vital both as a spectroscopic tool6,7,8 and also for applications of a large class of materials such as II-VI and III-V semiconductors9,10,11,12, nanotubes13, magnetic- and nonmagnetic layered materials14,15,16,17,18,19,20, perovskites and antiferromagnetic oxides21,22.

The occurrence of SHG in a crystal is directly linked to its symmetry, where crystals with broken inversion symmetry can exhibit SHG and crystals which are inversion symmetric can typically not5. This dependence on symmetry makes nonlinear optics in atomically thin crystals such as graphene a very rich field of research23,24,25,26,27,28,29. Although pristine graphene is an inversion-symmetric crystal, the crystal symmetry or symmetry of the electronic states is potentially very sensitive to charges or single molecules in the layer vicinity30, an imbalance in population of valley states31, or electric currents32. However, graphene is a gapless material, so light-matter interaction cannot be controlled to the same extent as in atom-thin semiconductors based on transition metal dichalcogenides (TMD).

TMD monolayers such as MoS2 and WSe2 have broken inversion symmetry and show intrinsically strong SHG signals1,33,34. Excitons are strongly bound in these materials35 and govern optical processes also at room temperature, contrary to the model semiconductor GaAs10. This is crucial for nonlinear optics as excitons resonantly enhance light-matter interaction by several orders of magnitude35. As a result, in monolayers the excitonic contribution to SHG can be orders of magnitude higher than the intrinsic contribution from the crystal which exists also off-resonance17,18,36,37,38. An intrinsic challenge for exploiting SHG signals from nanostructures is the overall small signal, despite the giant response of a monolayer per unit thickness34,39,40. Hence, in the interest of both fundamental physics and applications, we need to uncover the exact origins of SHG in layered semiconducting materials to enable further tuning and amplification.

In this work we demonstrate strong, tunable SHG based on exciton resonances in nominally inversion-symmetric MoS2 bilayers with 2H stacking: Varying the laser excitation energy allows us to address resonantly not only intralayer and but also interlayer exciton states with large oscillator strength41, which makes MoS2 bilayers particularly interesting. At these specific energies we generate SHG signals in MoS2 bilayers which are (i) comparable in amplitude to off-resonance SHG monolayer signals and (ii) several orders of magnitude larger than reported in the literature for TMD bilayers with 2H stacking33,38,42,43. The prominent role of interlayer excitons in MoS2 bilayers allows for an efficient SHG tuning mechanism: in applied electric fields, we lift the degeneracy of the interlayer exciton states through the Stark effect. This results in an additional enhancement of the SHG signal at the interlayer exciton resonance by two orders of magnitude, surpassing the response of an ungated monolayer at this energy. We uncover the excitonic origin of the highly tunable SHG response through the linear polarization dependence of the SHG signal with respect to the armchair direction of the sample. We compare bilayer results with measurements in monolayers. We uncover variable SHG enhancement in trilayers located in the same van der Waals stacks. We present a model analysis of the emergence of the strong SHG in the bilayer, taking into account the impact of electric fields perpendicular to the bilayer as well as different sources of asymmetry.

## Results

### SHG in MoS2 bilayers on SiO2/Si

Very recently, details of the rich excitonic resonances in bilayer MoS2 have been uncovered41. We aim to distinguish the possible SHG generation that comes from a weak, intrinsic, crystal-related contribution from the possible excitonic contribution by using a tunable laser source. This allows comparing off-resonance with on-resonance experiments.

For the experimental results shown in Fig. 1, we investigate a 2H MoS2 bilayer deposited on 90 nm of SiO2 on top of a Si substrate, a very typical sample configuration (see Mtethods and Supplementary Information for details in the fabrication process). The bilayer sample (as well as the trilayer sample below) has been directly exfoliated from a 2H bulk crystal, i.e., without manual stacking individual layers. We ensure therefore 60° twist angle (and not 0° as for 3R stacking) between adjacent layers, as confirmed by measuring the energy separation of the A–B excitons and observing the interlayer exciton44. We use a pulsed Ti:Sapphire laser source coupled to an optical parametric oscillator (OPO, using 5 mW average power and 1 ps pulse duration of 80 MHz repetition rate correspond to 60 W peak power) and we scan twice the laser energy (i.e. 2 × EL) over an energy range that corresponds to the main optical transitions in the MoS2 bilayers, namely the interlayer and intralayer exciton resonances45,46,47,48.

Experiments are performed at a temperature of T = 4 K in vacuum in a confocal microscope (excitation/detection spot diameter of the order of the wavelength), see supplement and7 for details. We plot in Fig. 1b a series of SHG spectra from the bilayer for a number of selected laser energies between EL = 0.9 and 1.13 eV in steps of ≈6 meV. Each SHG spectrum is a single peak that shifts in energy as we vary the excitation laser energy EL. The spectral width of the SHG signal is limited by the laser pulse duration (ps).

We start with the off-resonance case: we do not detect any SHG signal within our typical integration time (2 minutes and excitation power of a few mW) when twice the laser energy (2 × EL) is tuned below the intralayer A energy (i.e. the lowest lying direct optical transition with large oscillator strength). But surprisingly for a crystal with inversion center, as we increase the laser energy, we see clear resonances in the SHG response, i.e. a strong variation in SHG intensity as a function of laser energy EL in Fig. 1b. In Fig. 1c, we plot the differential white light reflectivity spectrum of this sample on the same energy scale.

Comparing SHG spectra and white light reflectivity, we deduce that the maxima in the SHG intensity occur when 2 × EL is in resonance with the intralayer (A and B) and interlayer (IE) exciton energies44.

The surprisingly strong SHG signal could suggest a symmetry breaking in our nonencapsulated sample. A possible reason can be a drastically different dielectric environment experienced by top and bottom MoS2 layers from Fig. 1b.

### SHG in MoS2 bilayers encapsulated in hBN

To verify the impact of the dielectric environment, we have fabricated MoS2 bilayer samples encapsulated in high quality hexagonal BN49. This more symmetric dielectric environment50,51,52 eliminates potential sources of symmetry breaking. In addition, we also exclude the possibility of adsorbates or molecules at the sample surface in encapsulated samples which might impact the measurements as in surface SHG studies53 and suggested in graphene24,30. In Fig. 1e we plot the SHG spectra as a function of twice the laser energy 2 × EL using a smaller step size of ≈3 meV between two adjacent SHG spectra. The SHG signal coming from hBN is several orders of magnitude weaker than from the MoS2 bilayer, see33 and our data in the supplement. Our experiments in Fig. 1e show that also hBN-encapsulated MoS2 bilayers in a far more symmetric dielectric environment show SHG at the exciton resonances. Comparing the signal strength for the A-exciton resonance, we record an SHG signal about an order of magnitude weaker for the hBN-encapsulated sample (Fig. 1e) as compared to the nonencapsulated sample on SiO2 (Fig. 1b). This indicates that indeed a sizeable contribution to the SHG in bilayers results from a very different dielectric environment for the top layer as compared to the bottom layer. The resonances are spectrally narrower in the encapsulated samples, see Fig. 1e, consistent with TMD exciton transitions in high quality hBN-encapsulated layers54. This spectral narrowing in combination with cavity effects in encapsulated samples will modulate the wavelength dependence of SHG intensity. Comparison with white light reflectivity data in Fig. 1f allows attributing the different interlayer and intralayer transitions. The slight overall shift in energy of the exciton transitions for encapsulated versus nonencapsulated samples comes mainly from renormalization of all the Coulomb energies as the effective dielectric constant is different for the two bilayer samples.

Our surprising finding on interlayer exciton-mediated SHG merits further investigation. In Fig. 1d, we plot the polarization dependence of the SHG response at the interlayer exciton resonance. We excite with linearly polarized light and the strength of the SHG signal collected in the same polarization depends on how the crystallographic axes are aligned with respect to the laser polarization. We clearly observe a 6-fold rotational symmetry expected for the space group of a 2H bilayer33, see discussion below. This polarization dependence is a strong indication that the SHG signal is due to intrinsic effects linked to crystal and exciton symmetry. We have performed measurements at additional exciton resonances that give the same polarization dependence.

### Comparison of SHG in hBN-encapsulated mono-, bi-, and trilayers

We provide in Fig. 2a–c results for hBN-encapsulated mono-, bi-, and trilayers, for the same sample used in ref. 41. Our aim is to compare the surprising SHG amplitude for bilayers with the monolayer and trilayers, for which crystal inversion symmetry is broken and as a consequence strong SHG is expected33. Strikingly, for all the experiments on the mono-, bi-, and trilayers, we see that the SHG signal is orders of magnitudes enhanced when twice the laser energy 2 × EL is in resonance with an excitonic transition, as compared with a nonresonant situation.

Previous reports on SHG in MoS2 bilayers did not focus on exciton resonances1,33,38,42,43,55 and hence signals from monolayers were 3 orders of magnitude higher than for bilayers in the off-resonant case. From our measurements in panels Fig. 2a, b we deduce that EL = 0.95 eV (i.e. SHG energy of 1.9 eV) is below the A-intralayer exciton resonance for monolayers and bilayers. We directly compare the measured SHG signal at this laser energy for monolayers and bilayers, see peaks in Fig. 2e, and find that the monolayer signal is indeed 2 orders of magnitude larger than the bilayer signal. The situation is drastically different as we change laser energy: At EL = 1.005 eV twice the laser energy 2 × EL is in resonance with the interlayer IE transition of the bilayer resulting in a strong amplification of the SHG signal while being nonresonant for a monolayer, see comparison in Fig. 2d. These different situations with respect to the exciton resonance result in a comparable overall SHG amplitude for mono- and bilayers.

Although SHG signals in trilayers higher than in bilayers have been reported at fixed laser energy, very little is known of the energy dependence of the nonlinear susceptibility of the trilayer. The trilayer crystal has no inversion center, as the monolayer, and sizeable nonresonant SHG can be expected. In addition, there exist several intralayer and interlayer exciton resonances41 that can potentially enhance SHG. For the measurements in Fig. 2c, we find an SHG signal with strong amplitude variations over the investigated energy range 1.87 and 2.15 eV. We identify the two intralayer transitions for excitons in the central layer (AL2) and in the outer layers (AL1 and AL3) respectively, between 1.9 and 1.92 eV. Local maxima around 1.98 eV can be attributed to interlayer excitons, with B-exciton contributions around 2.1 eV. From our measurements, we conclude that the trilayer combines the high intrinsic SHG response of the monolayer with the tunability (see below) of the interlayer exciton-mediated SHG. The exciton transitions for mono- and trilayers are identified and labeled in more detail by comparing with white light reflectivity41 shown in the supplement.

The A-exciton resonance in trilayers for the middle layer L2 is redshifted compared to the A-exciton resonance for the outer layers L1 and L3. Because of this shift, these contributions do not cancel each other as could be expected from symmetry if they would be at the same energy. The A-exciton of the monolayer gives about 1500 counts (Fig. 2a), the A-exciton of L2 gives also about 1500 counts and at the energy of A-exciton for L1 and L3 we get about 3000 counts (Fig. 2c). If for the trilayer the A-excitons in all 3 layers would be degenerate, we would expect to measure roughly only 3000 − 1500 = 1500 overall counts for the A-excitons. But in an MoS2 trilayer this degeneracy among the A-excitons is lifted and we record in total 4500 counts in Fig. 2c.

In addition to measurements at T = 4 K, we have carried out SHG spectroscopy for mono-, bi-, and trilayer MoS2 also at room temperature, see supplement. We observe the typical broadening of the excitonic transitions and the bandgap shift to lower energy7. It has been reported that the excitonic enhancement of the SHG signal also applies to room temperature experiments on monolayers56,57. Here we show that this is also true for MoS2 bilayers and trilayers when twice the laser energy is tuned into resonance with interlayer and intralayer excitons, underlining the importance of exciton-mediated SHG. This is of practical relevance as changing the laser energy in order to enhance the SHG signal is simpler in some experiments than changing the sample or device temperature.

In agreement with earlier estimations58,59, we find a power of the SHG emitted by the monolayers of the order of a few to 10 fW for 3 mW laser excitation power. For the bilayer, the maximum measured SHG power is of the order of 1 fW for 5 mW excitation. For the trilayer we find similar power as for the monolayer. We discuss the sheet polarizability and the conversion efficiency in more detail in the Supplementary Information.

### Interlayer exciton-mediated SHG tuning

We now show that by applying an electric field normal to the bilayer, we can further increase the SHG signal specifically at the interlayer exciton resonance. The interlayer exciton has an in-built, static electric dipole. In the absence of electric fields, there are two degenerate interlayer exciton configurations, with a hole delocalized over both layers and an electron in the top or bottom layer, respectively, as shown in Fig. 3d. This strong absorption feature is highly tunable in energy through the Stark effect41,48,60. In Fig. 3, we show that the SHG signal is also tunable in amplitude and also spectrally when an external electric field is applied to the bilayer. In previous works, the effect of the electric field application on SHG response of 2H bilayers42,43 has been addressed only in terms of the impact of doping and crystal symmetry breaking, while interlayer exciton tuning as in our resonant SHG experiment has not been observed.

In Fig. 3a, we compare the SHG signal of a bilayer with and without an applied electric field. For this gated sample, the IE SHG signal reaches 5 counts/s at Fz = 0. As we apply a field of Fz = 0.17 MV/cm this signal increases by a factor of 25 to about 125 counts/s at the IE resonance maximum. To study this tuning in more detail, we plot in Fig. 3b the SHG amplitude as a function of the applied electric field for six separate experiments. We see a quadratic increase of the IE SHG signal as a function of applied electric field Fz. This quadratic increase is due to mixing with intralayer excitons, as we show below.

In contrast, for the intralayer excitons, we do not record any significant increase of the SHG amplitude as Fz is increased. At Fz = 0.17 MV/cm the SHG resonance shows a slight red-shift (Fig. 3a), possibly indicating charging effect, i.e., a shift toward trion transitions. Strong charging effects as a possible source of symmetry breaking in bilayer SHG response are discussed for WSe243. In our experiments, the target is to apply a static electric field and not to generate charging effects. As we trace the evolution of the neutral interlayer exciton as a function of the applied field, we rely on the Stark effect and not charging effects for the increase in SHG amplitude.

We now turn our attention to the spectral range over which the SHG enhancement occurs. In a 2H homobilayer at Fz = 0 two degenerate IE configurations exist: with a hole delocalized over both layers bound to an electron either in the bottom (IE1) or top (IE2) layer, see sketch in Fig. 3d. Application of a nonzero Fz lifts the degeneracy of IE1 and IE2 as their static, permanent dipoles point in opposite directions. This leads to a Stark shift to lower and higher energy, respectively41. This lifting of the degeneracy can be seen in Fig. 3c in white light reflectivity experiments. We extract a splitting between IE1 and IE2 of roughly 10 meV for Fz = 0.13 MV/cm. In Fig. 3e, we plot SHG spectroscopy results at the IE resonance for Fz = 0, 0.07 and 0.13 MV/cm, respectively. In addition to the increase in SHG amplitude with Fz, we therefore observe that the spectral range for which amplification is observed is widened. This broadening is the consequence of the energy splitting between IE1 and IE2 transitions in an applied electric Fz. The broadening can also be seen directly in Fig. 3a. Tuning both amplitude and spectral range of the SHG response of the bilayer is a direct consequence of exploiting the properties of interlayer excitons applicable to many van der Waals bilayer systems61,62,63,64,65,66 and will be further discussed below.

## Discussion

We first focus on the model considerations for SHG mediated by the interlayer excitons in bilayer MoS2. Here the SHG signal is very sensitive to the electric field, being one of the most striking experimental observations. Further we briefly address other contributions to the SHG. The detailed theory is presented in the supplement.

A TMD monolayer is described by the noncentrosymmetric point group D3h which lacks an inversion center and allows for the SHG. An ideal homobilayer in absence of external fields has an inversion symmetry and is described by the D3d point group where the SHG is forbidden. Qualitatively, this is because the constituent monolayers in the bilayer are rotated by an angle of π with respect to each other in 2H stacking67,68, and the electric field at a double frequency generated in one monolayer is compensated by the contribution of the other, Fig. 1a.

An electric field F perpendicular to the monolayers breaks the symmetry and enables the SHG described by a phenomenological relation

$${P}_{2\omega ,x}=\chi ({E}_{\omega ,x}^{2}-{E}_{\omega ,y}^{2}),\quad {P}_{2\omega ,y}=-2\chi {E}_{\omega ,x}{E}_{\omega ,y},$$
(1)

where Eω is the incident field and P2ω is the polarization induced at a double frequency, x and y are the main in-plane axes of the structure. Since the bilayer is centrosymmetric, the nonlinear susceptibility χ in Eq. (1) changes sign at the space inversion, particularly, χ(−Fz) = −χ(Fz). This relation reflects the three-fold symmetry of the bilayer where the quadratic components of the vector component $${E}_{x}^{2}-{E}_{y}^{2}$$, −2ExEy transform as vector components Px, Py, see ref. 37 for details. Polarization dependence of the SHG predicted by Eq. (1) is observed experimentally, Fig. 1(d), showing sixfold symmetry of the crystal structure.

Intralayer excitons show negligible Stark shifts in applied fields41, as their static dipole moment is negligible, as for excitons in separate monolayers69,70. Therefore, due to their strong in-built electric dipole, the main effect of the field Fz is on the indirect excitons IE1 and IE2 whose degeneracy is lifted yielding $${E}_{{{{{{{{{\rm{IE}}}}}}}}}_{i}}={E}_{{{{{{{{\rm{IE}}}}}}}}}\mp {F}_{z}\mu$$, where ±μ is the static dipole moment of the corresponding IE, Fig. 3(d). As a result of the energy splitting (Fig. 3c), the SHG contributions from IE1 and IE2 do not cancel, resulting in a strong SHG signal as the energy separation between IE1 and IE2 increases (Fig. 3a). The second-order susceptibility in the vicinity of the IE resonances in the linear-in-Fz regime can be calculated taking into account two-photon excitation and coherent one-photon emission from the exciton states37 as (see supplement for details):

$${\chi }_{{{{{{{{\rm{IE}}}}}}}}}=-2\mu {F}_{z}\frac{| T{| }^{2}}{{{{\Delta }}}^{2}}\frac{{C}_{2}| {{{\Phi }}}_{{{{{{{{\rm{B:1s}}}}}}}}}(0){| }^{2}}{{\left(2\hslash \omega -{E}_{{{{{{{{\rm{IE}}}}}}}}}+\frac{| T{| }^{2}}{{{\Delta }}}+{{{{{{{\rm{i}}}}}}}}{{{\Gamma }}}_{{{{{{{{\rm{IE}}}}}}}}}\right)}^{2}}.$$
(2)

Here C2 is the combination of the interband momentum matrix element, electron charge, free-electron mass and the band gap, ΦB:1s(0) is the envelope function of the B:1s exciton at the coinciding electron and hole coordinates, T is the hole tunneling matrix element between the layers, Δ = EB − EIE is the energy splitting between the IE and B-excitons (which are mixed due to the hole tunneling45,47), and ΓIE is the damping of the IE resonance. As demonstrated in previous works41,60, the IE exciton optical activity is mainly due to the mixing with B-excitons, thus, in Eq. (2) the ratio T2 accounts for the IE-B exciton mixing in the two-photon excitation and single-photon emission channel. Equation (2) holds provided μFz ΓIE, otherwise two peaks in the SHG spectra are expected to split by 2μFz (see supplement for details).

Two important conclusions can be drawn from Eqs. (1) and (2): First, P2ω is maximized when twice the laser energy is resonant with an indirect exciton state as 2 × EL = 2ω = EIE. Second, the intensity of the SHG scales as $${F}_{z}^{2}$$. Experimental SHG intensities in Fig. 3a, e as well as in Figs. 1b, e and 2b are strongly enhanced at the IE resonance. In the measurements in Fig. 3b the SHG amplitude grows quadratically as a function of the applied Fz. So both experimental observations strongly support our analysis.

Importantly, the experiment demonstrates also the SHG on other excitonic species in bilayers, including weakly field-dependent effect on intralayer A-excitons, Fig. 3(a, b), and also the SHG at zero field on intralayer A:2s and B-excitons. A weaker effect of the electric field F on the A-excitons can be readily understood taking into account the fact that the A-excitons are mixed, e.g., with the interlayer B-excitons via hole tunneling. Corresponding energy distance $${{\Delta }}^{\prime} ={E}_{{{{{{{{\rm{IE(B)}}}}}}}}}-{E}_{{{{{{{{\rm{A}}}}}}}}}\gg {{\Delta }}$$ resulting in a weaker susceptibility to Fz.

We now discuss the enhancement of the SHG at the intralayer exciton resonances in the nominally centrosymmetric situation when Fz = 0. Possible origin of SHG in this case could be the quadrupolar or magneto-dipolar SHG, where akin to graphene case23,27, the effect is related to the light wavevector qz at the normal incidence. The symmetry analysis again yields Eq. (1) with χqz. Particularly, a (tiny) phase difference of ϕ ~ qza where a is the interlayer distance in the bilayer yields an imbalance of the contributions of individual monolayers and results in the resonant contributions at A- and B-excitons similar to calculated in37 for monolayers, but being by a factor of ϕ ~ 10−2…10−3 smaller (see supplement). Our estimates show, however, that this effect is too small to account on its own for the surprising experimental findings which demonstrate that the resonant susceptibilities of the bilayer are roughly one order of magnitude smaller than those of the monolayer.

Thus, we arrive at the conclusion that, despite relatively symmetric environment of our bilayer sample, the structure lacks an inversion center. Possible options are (i) small built-in electric fields and (ii) inequivalence of the intralayer excitons in monolayers forming a bilayer. Option (i) results in the replacement of Fz by Fz + F0 where F0 is the normal component of the build-in field. In this case at Fz = −F0 the effectively symmetric situation can be realized. While the measured dependence of the SHG intensity for IE, Fig. 3, does not contradict this scenario, note the minimum at Fz ≈ − 0.035 MV/cm, the SHG intensity measured at the A-exciton resonance does not significantly drop in this electric field range. Alternatively, we may suppose that the build-in field is inhomogeneous in the sample plane within our detection spot and the contributions of the A-exciton and IE states come from slightly different nanoscopic regions.

Option (ii) implies that the energies of A- and B-excitons, their nonradiative broadening, or their oscillator strengths are inequivalent in the two monolayers that form the bilayer. This may be due to the dielectric disorder, which, although suppressed, can still be present in state-of-the-art, hBN-encapsulated samples51. In this situation, the IE is mainly activated by the electric field, while intralayer excitons are less sensitive to Fz (see supplement). Inhomogeneous broadening effects due to disorder, impurities, etc. are inevitable and inhomogeneities on the order of one micrometer can occur even in hBN-encapsulated samples51,52. For the in-plane disorder the crucial parameter is the ratio between the size of the inhomogeneity, l, and the laser spot diameter, δ. For lδ, inhomogeneous broadening effects will effectively cancel out and the structure will effectively maintain the space inversion. However, for lδ, the inhomogeneities will lift the inversion symmetry, thus a contribution of one monolayer will dominate over the other and SHG occurs. It is important to note that, the asymmetry of the hBN environment due to the different thicknesses of the top and bottom layers and ‘vertical’ disorder caused, e.g., by the fluctuations of the top layer thickness, can result in the SHG at Fz = 0 at the A-exciton, e.g., due to the difference of the intralayer A-exciton energies and their nonradiative dampings. We stress that the six-fold symmetry of the SHG signal (see supplement) clearly rules out, e.g., an in-plane asymmetry due to the strain or in-plane fields20.

In summary, by varying excitation laser energy, dielectric environment, and applied electric fields we show strong and tunable exciton-mediated SHG in 2H MoS2 bilayers that can surpass the off-resonance SHG signal in monolayers. Drastic enhancement of the bilayer SHG amplitude is observed when twice the laser energy is in resonance with the excitonic transitions. The SHG signal outside exciton resonances remains orders of magnitude smaller than the monolayer signal, indicating a very small intrinsic χ value from broken inversion symmetry. At the interlayer exciton resonance, we tune the SHG signal by over an order of magnitude in electric fields applied perpendicular to the layer and demonstrate that the spectral width of the SHG resonance increases. With our model calculations, we relate the SHG in applied electric fields to the Stark splitting of the interlayer exciton and its mixing with intralayer excitons, which results in a quadratic dependence of the SHG amplitude on the applied electric field.

Our scheme for tuning SHG based on interlayer excitons with a permanent static electric dipole can be applied to a variety of other systems with strong interlayer exciton resonances, such as homobilayer MoSe266. Very importantly it has been shown recently that several heterobilayer systems host interlayer excitons with static electric dipoles and high optical oscillator strength, for example, MoSe2/WS2 with hybridized conduction states61,62, WSe2/WS2 with hybridized valence states63, MoS2/WS2 with electrically tunable valence state hybridization64 and similar predictions of electronic state hybridization for specific twist angles in MoTe2/MoSe265. It can be expected that also in these and other systems SHG can also be efficiently tuned via the Stark effect, as demonstrated here for homobilayer MoS2.

## Methods

### Sample fabrication

The 90-nm-thick SiO2/Si substrates were cleaned for 10 min in aceton and isopropanol using an ultrasonication bath and were subsequently exposed in oxygen plasma for 60 s. Bulk 2H MoS2 (2D Semiconductors) was first exfoliated on Nitto Denko tape and the exfoliated areas were attached on a polydimethylsiloxane (PDMS) stamp, supported by a microscope glass slide. Monolayers, 2H bilayers, and trilayers were identified based on the optical contrast under an optical microscope prior to transfer on the SiO2/Si substrate. For the hBN-encapsulated samples, hBN flakes were first exfoliated on a Nitto Denko tape from high-quality bulk crystals while the same PDMS-assisted transfer process on SiO2/Si substrates was followed. A staircase sample of monolayers, 2H bilayers and trilayers was susbequently transferred and capped in hBN. Between every transfer step, annealing at 150 °C was applied for 60 min. Optical images of hBN encapsulated but also bare MoS2 in SiO2 are shown in Supplementary Information. The thickness of the bottom hBN layers was selected to optimize the visibility of the interlayer exciton in the reflectivity spectra. For the electric field device, the same process was followed including the additional transfer of few-layered graphite (FLG) flakes. The stack, starting from bottom to top consists of Si, 90 nm SiO2, 130 nm hBN, FLG, few-layered (≈15 nm) hBN, 2L MoS2, few-layered hBN (≈20 nm), and FLG again. The bottom and top FLG are in contact with gold (Au) electrodes to apply a potential difference and generate an electric field perpendicular to the structure. The precise sequence of the complete stack from bottom to top includes hBN/FLG/hBN/2H-MoS2/hBN/FLG. A schematic representation of the device is shown in Supplementary Information.

### Experimental setup

Optical spectroscopy is performed in a home-built micro-spectroscopy setup assembled around a closed-cycle, low vibration helium cryostat with a temperature controller (T = 4 K to 300 K). For SHG measurements, we use ps pulses, generated by a tunable OPO synchronously pumped by a mode-locked Ti:sapphire laser. SHG signal is collected in reflection geometry. A combination of linear polarizers and halfwave plates allows the control of excitation and detection polarization for the polarization-resolved measurements with the setup sketched in Supplementary Information. The light is focused onto the sample at T = 4.2 K using a microscope objective (NA = 0.8). The position of the sample with respect to the focus can be adjusted with cryogenic nanopositioners. The reflected light from the sample is sent to a spectrometer with a 150 grooves per millimeter grating. The spectra are recorded by a liquid-nitrogen cooled charged coupled device (CCD) array. For low temperature white light reflectance measurements a white light source; a halogen lamp is used with a stabilized power, focused initially on a pin-hole that is imaged on the sample. The emitted and/or reflected light was dispersed in a spectrometer and detected by a Si-CCD camera. The excitation/detection spot diameter is 1 μm, i.e., smaller than the typical size of the homobilayers. We obtained differential reflectivity from reflectivity spectra as (RML − Rsub)/Rsub, where RML is the intensity reflection coefficient of the sample with the MoS2 layer and Rsub is the reflection coefficient of the hBN/SiO2 stack.

## Data availability

The data that support the findings of this study are available from the corresponding authors upon request.

## References

1. 1.

Zhang, J., Zhao, W., Yu, P., Yang, G. & Liu, Z. Second harmonic generation in 2D layered materials. 2D Mater. 7, 042002 (2020).

2. 2.

Trovatello, C. et al. Optical parametric amplification by monolayer transition metal dichalcogenides. Nat. Photonics 15, 6–10 (2020).

3. 3.

Franken, P. A., Hill, A. E., Peters, C. W. & Weinreich, G. Generation of optical harmonics. Phys. Rev. Lett. 7, 118–119 (1961).

4. 4.

Terhune, R. W., Maker, P. D. & Savage, C. M. Optical harmonic generation in calcite. Phys. Rev. Lett. 8, 404–406 (1962).

5. 5.

Yariv, A. & Yeh, P. Photonics: Optical Electronics in Modern Communications, Vol. 6 (Oxford University Press, New York, 2007).

6. 6.

Heinz, T. F., Chen, C., Ricard, D. & Shen, Y. Spectroscopy of molecular monolayers by resonant second-harmonic generation. Phys. Rev. Lett. 48, 478 (1982).

7. 7.

Shree, S., Paradisanos, I., Marie, X., Robert, C. & Urbaszek, B. Guide to optical spectroscopy of layered semiconductors. Nat. Rev. Phys. 3, 39–54 (2021).

8. 8.

Fiebig, M., Pavlov, V. V. & Pisarev, R. V. Second-harmonic generation as a tool for studying electronic and magnetic structures of crystals. JOSA B 22, 96–118 (2005).

9. 9.

Lafrentz, M. et al. Second-harmonic generation spectroscopy of excitons in ZnO. Phys. Rev. B 88, 235207 (2013).

10. 10.

Bergfeld, S. & Daum, W. Second-harmonic generation in GaAs: experiment versus theoretical predictions of $${\chi }_{xyz}^{(2)}$$. Phys. Rev. Lett. 90, 036801 (2003).

11. 11.

Yakovlev, D. R. et al. Exciton spectroscopy of semiconductors by the method of optical harmonics generation (review). Phys. Solid State 60, 1471–1486 (2018).

12. 12.

Lu, X., Moille, G., Rao, A., Westly, D. A. & Srinivasan, K. Efficient photoinduced second-harmonic generation in silicon nitride photonics. Nat. Photonics 15, 131–136 (2021).

13. 13.

De Dominicis, L. et al. Second- and third- harmonic generation in single-walled carbon nanotubes at nanosecond time scale. Appl. Phys. Lett. 85, 1418–1420 (2004).

14. 14.

Sun, Z. et al. Giant nonreciprocal second-harmonic generation from antiferromagnetic bilayer CrI3. Nature 572, 497–501 (2019).

15. 15.

Löchner, F. J. et al. Controlling second-harmonic diffraction by nano-patterning MoS2 monolayers. Opt. Express 27, 35475–35484 (2019).

16. 16.

Cunha, R. et al. Second harmonic generation in defective hexagonal boron nitride. J. Phys.: Condens. Matter 32, 19LT01 (2020).

17. 17.

Seyler, K. L. et al. Electrical control of second-harmonic generation in a WSe2 monolayer transistor. Nat. Nanotechnol. 10, 407–411 (2015).

18. 18.

Wang, G. et al. Giant enhancement of the optical second-harmonic emission of WSe2 monolayers by laser excitation at exciton resonances. Phys. Rev. Lett. 114, 097403 (2015).

19. 19.

Lin, K.-Q., Bange, S. & Lupton, J. M. Quantum interference in second-harmonic generation from monolayer WSe2. Nat. Phys. 15, 242–246 (2019).

20. 20.

Mennel, L. et al. Optical imaging of strain in two-dimensional crystals. Nat. Commun. 9, 1–6 (2018).

21. 21.

Abdelwahab, I. et al. Highly enhanced third-harmonic generation in 2D perovskites at excitonic resonances. ACS Nano 12, 644–650 (2018).

22. 22.

Fiebig, M., Frohlich, D., Leute, S. & Pisarev, R. V. Second harmonic spectroscopy and control of domain size in antiferromagnetic YMnO3. J. Appl. Phys. 83, 6560–6562 (1998).

23. 23.

Glazov, M. M. & Ganichev, S. D. High frequency electric field induced nonlinear effects in graphene. Phys. Rep. 535, 101–138 (2014).

24. 24.

Dean, J. J. & van Driel, H. M. Second harmonic generation from graphene and graphitic films. Appl. Phys. Lett. 95, 261910 (2009).

25. 25.

Hendry, E., Hale, P. J., Moger, J., Savchenko, A. & Mikhailov, S. A. Coherent nonlinear optical response of graphene. Phys. Rev. Lett. 105, 097401 (2010).

26. 26.

Mikhailov, S. A. Theory of the giant plasmon-enhanced second-harmonic generation in graphene and semiconductor two-dimensional electron systems. Phys. Rev. B 84, 045432 (2011).

27. 27.

Glazov, M. M. Second harmonic generation in graphene. JETP Lett. 93, 366–371 (2011).

28. 28.

Gullans, M., Chang, D., Koppens, F., de Abajo, F. G. & Lukin, M. D. Single-photon nonlinear optics with graphene plasmons. Phys. Rev. Lett. 111, 247401 (2013).

29. 29.

Soavi, G. et al. Hot electrons modulation of third-harmonic generation in graphene. ACS Photonics 6, 2841–2849 (2019).

30. 30.

Yu, R., Cox, J. D. & de Abajo, F. J. G. Nonlinear plasmonic sensing with nanographene. Phys. Rev. Lett. 117, 123904 (2016).

31. 31.

Golub, L. & Tarasenko, S. Valley polarization induced second harmonic generation in graphene. Phys. Rev. B 90, 201402 (2014).

32. 32.

Bykov, A. Y., Murzina, T. V., Rybin, M. G. & Obraztsova, E. D. Second harmonic generation in multilayer graphene induced by direct electric current. Phys. Rev. B 85, 121413 (2012).

33. 33.

Li, Y. et al. Probing symmetry properties of few-layer MoS2 and h-BN by optical second-harmonic generation. Nano Lett. 13, 3329–3333 (2013).

34. 34.

Malard, L. M., Alencar, T. V., Barboza, A. P. M., Mak, K. F. & de Paula, A. M. Observation of intense second harmonic generation from MoS2 atomic crystals. Phys. Rev. B 87, 201401 (2013).

35. 35.

Wang, G. et al. Colloquium: Excitons in atomically thin transition metal dichalcogenides. Rev. Mod. Phys. 90, 021001 (2018).

36. 36.

Wang, G. et al. Exciton states in monolayer MoSe2: impact on interband transitions. 2D Mater. 2, 045005 (2015).

37. 37.

Glazov, M. et al. Intrinsic exciton-state mixing and nonlinear optical properties in transition metal dichalcogenide monolayers. Phys. Rev. B 95, 035311 (2017).

38. 38.

Zhao, M. et al. Atomically phase-matched second-harmonic generation in a 2D crystal. Light. Sci. Appl. 5, e16131–e16131 (2016).

39. 39.

Wen, X., Gong, Z. & Li, D. Nonlinear optics of two-dimensional transition metal dichalcogenides. InfoMat 1, 317–337 (2019).

40. 40.

Weismann, M. & Panoiu, N. C. Theoretical and computational analysis of second-and third-harmonic generation in periodically patterned graphene and transition-metal dichalcogenide monolayers. Phys. Rev. B 94, 035435 (2016).

41. 41.

Leisgang, N., Shree, S. & Paradisanos, I. et al. Giant stark splitting of an exciton in bilayer MoS2. Nat. Nanotechnol. 15, 901–907 (2020).

42. 42.

Klein, J. et al. Electric-field switchable second-harmonic generation in bilayer MoS2 by inversion symmetry breaking. Nano Lett. 17, 392–398 (2017).

43. 43.

Yu, H., Talukdar, D., Xu, W., Khurgin, J. B. & Xiong, Q. Charge-induced second-harmonic generation in bilayer WSe2. Nano Lett. 15, 5653–5657 (2015).

44. 44.

Paradisanos, I. et al. Controlling interlayer excitons in MoS2 layers grown by chemical vapor deposition. Nat. Commun. 11, 1–7 (2020).

45. 45.

Deilmann, T. & Thygesen, K. S. Interlayer excitons with large optical amplitudes in layered van der Waals materials. Nano Lett. 18, 2984–2989 (2018).

46. 46.

Slobodeniuk, A. et al. Fine structure of k-excitons in multilayers of transition metal dichalcogenides. 2D Mater. 6, 025026 (2019).

47. 47.

Gerber, I. C. et al. Interlayer excitons in bilayer MoS2 with strong oscillator strength up to room temperature. Phys. Rev. B 99, 035443 (2019).

48. 48.

Peimyoo, N. et al. Electrical tuning of optically active interlayer excitons in bilayer MoS2. Nat. Nanotechnol. 16, 888–893 (2021).

49. 49.

Taniguchi, T. & Watanabe, K. Synthesis of high-purity boron nitride single crystals under high pressure by using Ba–BN solvent. J. Cryst. Growth 303, 525–529 (2007).

50. 50.

Stier, A. V., Wilson, N. P., Clark, G., Xu, X. & Crooker, S. A. Probing the influence of dielectric environment on excitons in monolayer WSe2: insight from high magnetic fields. Nano Lett. 16, 7054–7060 (2016).

51. 51.

Raja, A. et al. Dielectric disorder in two-dimensional materials. Nat. Nanotechnol. 14, 832–837 (2019).

52. 52.

Rhodes, D., Chae, S. H., Ribeiro-Palau, R. & Hone, J. Disorder in van der Waals heterostructures of 2D materials. Nat. Mater. 18, 541 (2019).

53. 53.

Shen, Y. Surface properties probed by second-harmonic and sum-frequency generation. Nature 337, 519–525 (1989).

54. 54.

Cadiz, F. et al. Excitonic linewidth approaching the homogeneous limit in MoS2-based van der Waals heterostructures. Phys. Rev. X 7, 021026 (2017).

55. 55.

Stiehm, T. et al. Supercontinuum second harmonic generation spectroscopy of atomically thin semiconductors. Rev. Sci. Instrum. 90, 083102 (2019).

56. 56.

Lafeta, L. et al. Second-and third-order optical susceptibilities across excitons states in 2D monolayer transition metal dichalcogenides. 2D Mater. 8, 035010 (2021).

57. 57.

Khan, A. R. et al. Extraordinary temperature dependent second harmonic generation in atomically thin layers of transition-metal dichalcogenides. Adv. Optical Mater. 8, 2000441 (2020).

58. 58.

Klimmer, S. et al. All-optical polarization and amplitude modulation of second-harmonic generation in atomically thin semiconductors. Nat. Photonics 15, 837–842 (2021).

59. 59.

Woodward, R. et al. Characterization of the second-and third-order nonlinear optical susceptibilities of monolayer MoS2 using multiphoton microscopy. 2D Mater. 4, 011006 (2016).

60. 60.

Lorchat, E. et al. Excitons in bilayer MoS2 displaying a colossal electric field splitting and tunable magnetic response. Phys. Rev. Lett. 126, 037401 (2021).

61. 61.

Alexeev, E. M. et al. Resonantly hybridized excitons in moiré superlattices in van der Waals heterostructures. Nature 567, 81–86 (2019).

62. 62.

Zhang, L. et al. Twist-angle dependence of moiré excitons in WS2/MoSe2 heterobilayers. Nat. Commun. 11, 1–8 (2020).

63. 63.

Tang, Y. et al. Tuning layer-hybridized moiré excitons by the quantum-confined Stark effect. Nat. Nanotechnol. 16, 52–57 (2021).

64. 64.

Kiemle, J. et al. Control of the orbital character of indirect excitons in MoS2/WS2 heterobilayers. Phys. Rev. B 101, 121404 (2020).

65. 65.

Ruiz-Tijerina, D. A. & Fal’ko, V. I. Interlayer hybridization and moiré superlattice minibands for electrons and excitons in heterobilayers of transition-metal dichalcogenides. Phys. Rev. B 99, 125424 (2019).

66. 66.

Horng, J. et al. Observation of interlayer excitons in MoSe2 single crystals. Phys. Rev. B 97, 241404 (2018).

67. 67.

Gong, Z. et al. Magnetoelectric effects and valley-controlled spin quantum gates in transition metal dichalcogenide bilayers. Nat. Commun. 4, 2053 (2013).

68. 68.

Kim, W., Ahn, J. Y., Oh, J., Shim, J. H. & Ryu, S. Second-harmonic Young’s interference in atom-thin heterocrystals. Nano Lett. 20, 8825–8831 (2020).

69. 69.

Roch, J. G. et al. Quantum-confined stark effect in a MoS2 monolayer van der Waals heterostructure. Nano Lett. 18, 1070–1074 (2018).

70. 70.

Verzhbitskiy, I., Vella, D., Watanabe, K., Taniguchi, T. & Eda, G. Suppressed out-of-plane polarizability of free excitons in monolayer WSe2. ACS Nano 13, 3218–3224 (2019).

71. 71.

Goryca, M. et al. Revealing exciton masses and dielectric properties of monolayer semiconductors with high magnetic fields. Nat. Commun. 10, 1–12 (2019).

## Acknowledgements

Toulouse acknowledges funding from ANR 2D-vdW-Spin, ANR MagicValley, ANR IXTASE, ANR HiLight, ITN 4PHOTON Marie Sklodowska Curie Grant Agreement No. 721394 and the Institut Universitaire de France. Growth of hexagonal boron nitride crystals was supported by the Elemental Strategy Initiative conducted by the MEXT, Japan, Grant Number JPMXP0112101001, JSPS KAKENHI Grant Number JP20H00354 and the CREST(JPMJCR15F3), JST. M.M.G. and L.E.G. acknowledge the RFBR and CNRS joint project 20-52-16303. L.E.G. was supported by the Foundation for the Advancement of Theoretical Physics and Mathematics “BASIS”.

## Author information

Authors

### Contributions

S.S. and I.P. performed the measurements, the data analysis. M.M.G, L.E.G, and I.C.G. performed the model calculations. T.T. and K.W. grew the high-quality boron nitride. C.R., D.L., A.B., and L.L. mounted the optical setup. S.S., I.P., and C.R. fabricated the van der Waals stacks. X.M. and B.U. supervised the project and suggested the experiment. All authors discussed the results and commented on the manuscript at all stages.

### Corresponding authors

Correspondence to Mikhail M. Glazov, Bernhard Urbaszek or Ioannis Paradisanos.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

Peer review information Nature Communications thanks Zhipei Sun, Xuejin Zhang and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Shree, S., Lagarde, D., Lombez, L. et al. Interlayer exciton mediated second harmonic generation in bilayer MoS2. Nat Commun 12, 6894 (2021). https://doi.org/10.1038/s41467-021-27213-8

• Accepted:

• Published:

• DOI: https://doi.org/10.1038/s41467-021-27213-8