# Fileset

[s41598-021-95700-5.pdf](https://mdr.nims.go.jp/filesets/b0af6c63-d08e-4f7e-be2b-c2b3291c6721/download)

## Creator

M. Grzeszczyk, J. Szpakowski, A. O. Slobodeniuk, T. Kazimierczuk, M. Bhatnagar, [T. Taniguchi](https://orcid.org/0000-0002-1467-3105), [K. Watanabe](https://orcid.org/0000-0003-3701-8119), P. Kossacki, M. Potemski, A. Babiński, M. R. Molas

## Rights

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

## Other metadata

[The optical response of artificially twisted MoS$$_2$$ bilayers](https://mdr.nims.go.jp/datasets/f86b28c7-925d-4c07-8eda-13a45d01d8c5)

## Fulltext

The optical response of artificially twisted MoS bilayers1Vol.:(0123456789)Scientific Reports |        (2021) 11:17037  | https://doi.org/10.1038/s41598-021-95700-5www.nature.com/scientificreportsThe optical response of artificially twisted MoS2 bilayersM. Grzeszczyk  1*, J. Szpakowski1, A. O. Slobodeniuk  2, T. Kazimierczuk  1, M. Bhatnagar  1, T. Taniguchi  3, K. Watanabe  4, P. Kossacki  1, M. Potemski  1,5, A. Babiński  1 & M. R. Molas  1*Two-dimensional layered materials offer the possibility to create artificial vertically stacked structures possessing an additional degree of freedom—the interlayer twist. We present a comprehensive optical study of artificially stacked bilayers (BLs) MoS2 encapsulated in hexagonal BN with interlayer twist angle ranging from 0 ◦ to 60◦ using Raman scattering and photoluminescence spectroscopies. It is found that the strength of the interlayer coupling in the studied BLs can be estimated using the energy dependence of indirect emission versus the A 1g–E12g energy separation. Due to the hybridization of electronic states in the valence band, the emission line related to the interlayer exciton is apparent in both the natural (2H) and artificial (62◦ ) MoS2 BLs, while it is absent in the structures with other twist angles. The interlayer coupling energy is estimated to be of about 50 meV. The effect of temperature on energies and intensities of the direct and indirect emission lines in MoS2 BLs is also quantified.Two-dimensional (2D) van der Waals (vdW) crystals have emerged as a new generation of materials with extraor-dinary properties. For instance, widely studied semiconducting transition metal dichalcogenides (S-TMDs) transform from indirect- to direct-band gap, optically-bright semiconductors when thinned down to a mon-olayer (ML), which results in unique electronic structures and consequent optical properties 1–4. The family of 2D layered materials grows day by day, hugely expanding the scope of possible phenomena to be explored in two dimensions. Growing is also the number of possible vdW heterostructures that one can create. Such 2D materials currently cover a vast range of properties allowing potential applications in, i.a. spintronic devices 5, optoelectronics 6,7, tunnel field-effect transistors 8,9, single-photon sources 10,11, and quantum information process-ing 12,13. Rapid advances in fabrication methods, like chemical vapor deposition (CVD) growth and mechanical exfoliation techniques, have contributed to increased interest in artificial stacking of different layered materials on top of each other. The simplistic approach of producing vertical vdW heterostructures without the constraints of crystal lattice mismatch enables integrating various 2D materials to create diverse systems with new electronic properties that are not present in pristine components. In addition to the selection of compounds in terms of their properties, a new degree of freedom has emerged: the twist angle between stacked layers, which gives rise to the group of the so-called twistronic materials 14,15. The twist angle is responsible for the occurrence of moiré patterns, that leads to new and intriguing phenomena, like the formation of secondary Dirac points in graphene on hexagonal boron nitride (hBN) 16,17 or hybridized (moiré) excitons in vdW heterostructures formed by stacked two S-TMD MLs 18–23.Raman scattering (RS) and photoluminescence (PL) spectroscopies are extensively used experimental tech-niques to characterize layered materials. Particularly, they can be used to determine the thickness of each S-TMD thin layer due to the energy dependence of phonon vibrations 24–31 as well as to the direct-indirect band-gap transformation 1–4.In this work, we investigate interlayer interactions in high-quality artificially stacked twisted MoS2 bilayers (BLs) encapsulated in hBN flakes using the RS and PL spectroscopies. Our results indicate that the interlayer coupling can be determined by the comparison of the emission energies of the indirect transition versus the energy separation between two basic intralayer phonon modes (A1g and E 12g ). The origin of the apparent emis-sion due to the interlayer excitons in both the natural (2H) and artificial (62◦ ) MoS2 BLs and its absence for the BLs with other twist angles is associated with hybridization of electronic states in the valence band. We also OPEN1Institute of Experimental Physics, Faculty of Physics, University of Warsaw, ul. Pasteura 5, 02‑093  Warsaw, Poland. 2Department of Condensed Matter Physics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 5, 121 16  Prague 2, Czech Republic. 3International Center for Materials Nanoarchitectonics, National Institute for Materials Science, 1‑1 Namiki, Tsukuba 305‑0044, Japan. 4Research Center for Functional Materials, National Institute for Materials Science, 1‑1 Namiki, Tsukuba 305‑0044, Japan. 5Laboratoire National des Champs Magnétiques Intenses, CNRS-UGA-UPS-INSA-EMFL, 25, Avenue des Martyrs, 38042 Grenoble, France. *email: magdalena.grzeszczyk@fuw.edu.pl; maciej.molas@fuw.edu.plhttps://orcid.org/0000-0001-6861-3098https://orcid.org/0000-0001-5798-0431https://orcid.org/0000-0001-6545-4167https://orcid.org/0000-0003-3712-371Xhttps://orcid.org/0000-0002-1467-3105https://orcid.org/0000-0003-3701-8119https://orcid.org/0000-0002-7558-1044https://orcid.org/0000-0001-8881-6618https://orcid.org/0000-0002-5591-4825https://orcid.org/0000-0002-5516-9415http://crossmark.crossref.org/dialog/?doi=10.1038/s41598-021-95700-5&domain=pdf2Vol:.(1234567890)Scientific Reports |        (2021) 11:17037  | https://doi.org/10.1038/s41598-021-95700-5www.nature.com/scientificreports/investigate the evolution of energies and intensities of the emission lines due to direct and indirect transitions with temperature ranging from 5 to 300 K.ResultsAtomic structures of twisted bilayers MoS2.  MoS2 belongs to the family of S-TMDs with chemical formula MX2 where M = Mo or W and X = S, Se or Te, which the most common crystallographic structure is a hexagonal phase. In that case, X–Mo–X atoms in a monolayer (1 L) of MoX2 are arranged in a trigonal pris-matic structure, which does not exhibit inversion symmetry. A BL formed by the stacking of two MLs exhibits an additional degree of freedom—the twist angle between the layers. Due to the hexagonal symmetry of a MX2 ML, the different arrangements are characterized by twist angles ranging from 0 ◦ to 60◦ . The schematic illustra-tion of three patterns of MoS2 BLs with the twist angles of 60◦ , 15◦ , and 0 ◦ are presented in Fig. 1. The most stable structures, which exhibit the strongest coupling, are found for twist angles equal to 0 ◦ and 60◦ ascribed correspondingly to the 3R and 2H stackings32,33. We do not use the label 3R for the 0 ◦ BL (see Fig. 1), because the 3R unit cell involves the atoms from three consecutive layers (2H unit cell is composed of two layers). While both the 2H and 3R polytypes are accessible in natural and CVD-grown MoS2 multilayers 33, the BLs with other twist angles can only be constructed by the artificial stacking of individual monolayers. It is important to point out that the inversion symmetry is restored in the 2H BL, while the zero-twist angle BL has neither inversion nor in-plane mirror.Raman scattering spectroscopy.  In order to investigate the interlayer coupling in the twisted BLs and therefore to distinguish their pristine properties, we compare their RS and PL spectra with those obtained for 1 L and 2 L exfoliated from 2H bulk MoS2 . Note that the twist angle between MLs forming BLs was determined using the second-harmonic generation (SHG) technique, which permits to reveal the crystallographic orienta-tion of individual MLs34. The measured RS spectra of the selected BLs with twist angles of 62◦ , 18◦ , and 6 ◦ accompanied with the ones of 1 L and 2 L MoS2 are presented in Fig. 2a, b.First, we focus on the low-energy range of the RS spectra, presented in Fig. 2. There are no modes related to rigid vibrations in 1 L affirming their interlayer nature. For 2 L, two modes are observed at ∼ 23 cm−1 and ∼ 44 cm−1 , which are associated with the in-plane shear (SM) and out-of-plane breathing (BM) vibrations of rigid layers, respectively 33,35,36. In twisted BLs with angles of 62◦ and 6 ◦ , the SM and BM modes are also apparent. This indicates that the coupling between adjacent layers in these two BLs is relatively strong due to their resem-blance to the natural 2H and 3R stackings. Note that the relatively smaller intensity of the SM peak in the 6 ◦ BLs as compared with the 2H and 62◦ ones can be attributed to the smaller interlayer bond polarizability for the 3R polytype 37. This effect was applied to identify the stacking order of S-TMD layers 30,37. On the other hand, only the BM-related peak can be observed for the sample with the 18◦ angle as well as for other BLs with angles in the range from ∼ 10◦ to ∼56◦ (data not shown). The higher energy range of all RS spectra shown in Fig. 2b is dominated by two phonon modes ascribed to the in-plane ( E12g ) and out-of-plane ( A1g ) intralayer vibrations. In few layers 2H MoS2 , the E12g ( A1g ) phonon mode experiences a red (blue) shift with increasing the layer thick-ness. The energy separation between those modes, i.e. �ω = ωA1g − ωE12g , is a useful tool to determine the number of MoS2 layers 24. As can be seen in the Figure, the �ω energy difference varies for different twist angles ( � ), which is summarized for all samples in Fig. 2c. The �ω values are largest for samples with a twist angle of 6 ◦ and 62◦ which correspond well to natural 2H and 3R BLs 38. In the intermediate cases (6◦ < � < 60◦ ), the �ω falls lower than for 2 L MoS2 , but it is considerably larger than in 1 L. The �ω can be also used to characterize the effective interlayer mechanical coupling strength or, in other words, the distance between layers in twisted BLs 32. In such case, we can conclude that the structures with the twist angle equal to 6 ◦ and 62◦ are characterized by the strongest coupling (the smallest interlayer distance), while the coupling strength is weaker (the interlayer distance is bigger) for other BLs with the twist angles substantially different from 0 and 60 degrees.Note that our RS results are consistent with previous works devoted to the twisted MoS2 BLs exfoliated on Si/SiO2 substrates 32,39–42.Figure 1.   Atomic structure and stacking order of bilayers MoS2 for different twist angles: 60◦ , 15◦ and 0 ◦ . Note that the side view demonstrates the energetically favourable structures for angles of 60◦ and 0 ◦.3Vol.:(0123456789)Scientific Reports |        (2021) 11:17037  | https://doi.org/10.1038/s41598-021-95700-5www.nature.com/scientificreports/Photoluminescence experiment.  Previous investigation of PL spectra of artificially twisted BLs MoS2 were carried out at room temperature 32,39,42. The PL spectroscopy at low temperature may provide however a more accurate analysis, which is due to the much lower emission linewidths 4,43. The PL spectra of the selected BLs with twist angles of 62◦ , 18◦ , and 6 ◦ accompanied with those of 1 L and 2 L MoS2 measured at low tem-perature (T = 5 K) are shown in Fig. 3a. Similarly to the aforementioned analysis of the RS, we begin with an examination of the PL spectra of the natural 1 L and 2 L MoS2 . The 1 L spectrum consists of two narrow emission lines apparent in the vicinity of the optical band gap (so-called A exciton), which can be ascribed to the neutral and charged excitons in accordance with previous reports 44–46. In contrast, the PL spectrum of the 2 L MoS2 is composed of two distinct emission bands: (1) the transitions in the vicinity of the direct A (XA ) and B (XB ) exci-tons formed at the K ± points of the Brillouin zone (BZ), which are observed correspondingly at about 1.93 eV and 2.07 eV; (2) the significantly much intense transitions, denoted as I and I’, apparent at about 1.5 eV, which Figure 2.   Normalized room-temperature Raman scattering spectra of thin layers MoS2 : monolayer (1 L), bilayer (2 L) and selected homobilayers with twist angles of 6 ◦ , 18◦ and 62◦ in the energy range of (a) interlayer and (b) intralayer vibrations measured under excitation of 2.41 eV laser light. The spectra are normalized to the intensity of the A 1g peak. (c) The energy separation between the A 1g and E 12g peaks, i.e. �ω = ωA1g − ωE12g , as a function of twist angle ( θ ) for all the studied twisted homobilayers. The corresponding �ω for 1 L and 2 Ls, represented by solid symbols, are also shown.Figure 3.   (a) Normalized low-temperature ( T = 5 K) photoluminescence spectra of thin MoS2 layers: monolayer (1 L), bilayer (2 L) and selected homobilayers with twist angles of 6 ◦ , 18◦ and 62◦ measured under excitation of 2.41 eV laser light. The spectra are normalized to the intensity of the X A emission. (b) The emission energy of the direct (XA and X B ) and indirect (I) transitions as a function of twist angle ( θ ) for all the studied twisted homobilayers. The full and open black circles represent correspondingly the X A and X B , while the I emissions are indicated by the open black diamonds. (c) The evolution of the indirect emission energy, i.e. I, versus the energy separation between the A 1g and E 12g peaks, i.e. �ω = ωA1g − ωE12g . The corresponding emission energies for the 2 L are also indicated with solid orange points.4Vol:.(1234567890)Scientific Reports |        (2021) 11:17037  | https://doi.org/10.1038/s41598-021-95700-5www.nature.com/scientificreports/are ascribed correspondingly to an indirect recombination process between the � and K point in the conduction band (CB) and Ŵ points in the valence band (VB) of the BZ 32,43,47. The PL spectra of the artificial twisted BLs also comprise two direct- and indirect-related bands. While the energies of direct transitions are hardly affected by the twist angles, the energies of the indirect ones change significantly by about 150 meV. The low-intensity emission bands, denoted with * and apparent between the aforementioned I and X A lines, can be described as a remaining part of the defect-related emission, which is reported commonly for 1-L MoS2 exfoliated on Si/SiO2 substrates 44–46.To appreciate the effect of the twist angle on the PL spectra, the energy evolution of the X A , X B and I transi-tions for all the studied samples as a function of twist angle is presented in Fig. 3b. As can be seen in the Figure, there is no effect of the twist angle on the X A and X B energies. This reflects a pure two-dimensional character of both the A and B excitons, i.e. these complexes are distributed spatially within a single layer even in a bulk form of MoS2 48–50. For the indirect transitions, we focus only on the I emission lines as they are observed for all studied samples. The emission reaches the lowest energies for border cases (6◦ and 62◦ ), which are similar to the value obtained for the natural 2 L MoS2 . For BLs with other twist angles, the I energies are substantially higher by about 150 meV and they are almost independent of twist angle. This effect can be associated with changes in the interlayer distance between MLs forming the studied BLs as a function of twist angle, as it was reported previously in Refs. 32,39.Note that the values of I energy and �ω parameter for BL with twist angle of about 30◦ are considerably smaller as compared to the corresponding values obtained for other angles within the range 6 ◦ < � < 60◦ (see Figs. 2c, 3b). It is surprising as both the results indicate opposite cases: the weaker/stronger coupling between ML from the RS/PL results, but it requires more sophisticated theoretical analysis devoted to the effect of twist angles and is out of the scope of this work.As the energy separation between the A 1g and E 12g peaks ( �ω ) can be considered as a probe of the interlayer distance in the studied artificial MoS2 BLs, one can plot the energy dependence of the indirect transition I as a function of �ω in Fig. 3c. The presented results can be divided into two groups: (1) for structures with twist angle close to 0 ◦ and 60◦ , i.e. �ω ∼23 cm−1 , the I emission energy equals approx. 1.5 eV; (2) for samples with twist angle significantly different from 0 ◦ and 60◦ , i.e. �ω ∼22 cm−1 , the energy dispersion of the I line is of about 80 meV centered at ∼1.58 eV. These results are consistent with previously obtained for stacked triangle-shaped MoS2 BLs grown on Si/SiO2 substrate using CVD technique32. Summarizing, we can assume that the analysis of the indirect emission versus the �ω can be used as a useful tool to probe the interlayer distance in the twisted homo-BLs. This allows to distinguish between two cases of twist angles, which are in the vicinity of 0 ◦ and 60◦ or the twist angle is in-between.While the X A and X B lines are observed for all studied samples, there is an additional emission line, labelled IL, apparent between the former ones for the 2 L and 62◦ samples, which emission was not reported so far. According to recent results of reflectance contrast (RC) experiments carried out on natural 2H MoS2 BLs at low temperature51–53, that line can be associated with the recombination of the so-called interlayer A exciton. The IL line originates from the hybridization of electronic states in the VB of the 2H-stacked BL due to the interlayer coupling (see Ref. 51 for details). It is important to mention that the observation of the IL emission is not strictly limited to the 2H phase, reported so far in Refs.51–53. As the IL emission is also observed in PL spectra measured on the BL with 62◦ , it indicates that the IL exciton is also present in twisted bilayers with angles varying slightly from exact 60◦ . It is similar to the case of hybridized excitons reported in MoSe2/WS2 heterobilayers 23.The strength of the related VB coupling is described by the interlayer hopping term, t⊥ . Based on a kp model of 2H BLs in the vicinity of K± points 51, the t⊥ parameter is given by:where �2HA-B is the A-B energy difference for natural 2H BL and �2Hv  represents the spin-orbit VB splitting without interlayer coupling. Note that the detailed investigation of the band structure of the 2H and 0 ◦ BLs in the vicinity of the K ± point of the BZ is performed in Supporting Information. According to that analysis, we assume that: (1) the �2Hv  is an order of magnitude larger than its counterpart in the CB �2Hc  54; (2) the binding energies of the A and B excitons are comparable; (3) the �2Hv  can be roughly approximated by the A-B energy difference for the BL with twist angle of 0 ◦ . Consequently, the diagram of relevant subbands at the K + point of the BZ in the MoS2 BLs with twist angle of 60◦ (2H) and 0 ◦ is shown in Fig. 4. As can be appreciated from the Figure, the t⊥ parameter can be evaluated using the energy separation between the A and B excitons measured in the 60◦ (2H) and 0 ◦ BLs. The extracted �2HA-B ∼180 meV (the corresponding �62◦A-B is about 10 meV smaller, which suggests the lower strength of interlayer coupling in twisted BL) and �6◦A-B ∼150 meV are in very good agreement with the reported ones for natural BLs51–53. Using these values, the extracted interlayer hopping parameter ( t⊥ ) is found to be on the order of 50 meV, which agrees very well with the value obtained using the RC experiment performed on the MoS2 BLs grown using CVD technique and encapsulated in hBN flakes (49 meV) 53.The energy evolution of both the direct (XA , X B , and IL) and indirect (I and I’) lines as a function of tempera-ture measured on BL (2 L) and selected homobilayers with twist angles of 6 ◦ , 18◦ and 62◦ is shown in Fig. 5. As can be appreciated in the Figure, all lines within the direct- and indirect-related transitions located at ∼ 2.0 eV and ∼ 1.5 eV are characterized by the same type of evolution upon increasing temperature, i.e. redshift and blueshift, respectively. The X A , X B , and IL peaks redshift when temperature is increased from 5 to 300 K. This can be associated with the reduction of the direct band gap resulting from the temperature expansion of those layers in lateral directions. Consequently, these evolutions can described by the relation proposed by O’Donnell et al.55, which expresses the temperature dependence of the band gap in terms of the average energy of acoustic phonons involved in the electron-phonon interaction 〈�ω〉 . The relation reads E(T) = E0 − S��ω�[coth(��ω�/2kBT)− 1] , (1)t⊥ =12√(�2HA-B)2 − (�2Hv )2,5Vol.:(0123456789)Scientific Reports |        (2021) 11:17037  | https://doi.org/10.1038/s41598-021-95700-5www.nature.com/scientificreports/where E0 stands for the band gap at absolute zero temperature, S is the coupling constant, and kB denotes the Boltzmann constant. We found that 〈�ω〉 stayed on nearly the same level, ∼ 21 meV, for the X A , X B , and IL transitions in 2H BL 47. As a consequence, we kept it fixed during fittings of all the experimental data shown in Fig. 5. It can be seen that the fitted curves correctly reproduce the energies of excitonic lines (see Fig. 5), which suggests that binding energies of investigated complexes do not depend on temperature. Interestingly, the overall redshift of X A lines are of about 60–70 meV, while the corresponding shift for the MoS2 BL exfoliated on Si/SiO2 substrate was found to be of almost 90 meV 43. Simultaneously, the crystal expansion across the layers leads to a larger separation between the layers, which results in the blueshift of the indirect band gap. The I emission lines experience almost monotonic blueshifts in the range of about 20–40 meV, while the shift of ∼ 80 meV is reported in Ref. 43. Both these results indicate that the hBN encapsulation of the BLs induces much smaller expansion of Figure 4.   Diagram of relevant subbands in the CB and VB at the K + point of the Brillouin zone in the bilayers MoS2 with twist angle of (a) 60◦ (2H) and (b) 0 ◦ . The red and green wavy lines show the A and B exciton transitions. The t⊥ and �2Hv  represent the interlayer hopping term and the spin-orbit VB splitting without interlayer coupling, respectively.Figure 5.   Temperature evolution of the energies of the transitions extracted from PL spectra measured on MoS2 : (a) bilayer (2 L) and selected homobilayers with twist angles of (b) 62◦ , (c) 18◦ and (d) 6 ◦ . The solid black curves are fits to the data obtained with the aid of Odonnell’s formula.6Vol:.(1234567890)Scientific Reports |        (2021) 11:17037  | https://doi.org/10.1038/s41598-021-95700-5www.nature.com/scientificreports/the flake in both directions (in lateral directions and across the layers) with increasing temperature as compared with the BLs exfoliated on Si/SiO2 substrate 43.The temperature evolution of the integrated intensities of the direct- and indirect-related transitions is pre-sented in Fig. 6. It can be seen that the total intensities of all direct transitions mostly decrease monotonically with increasing temperature. The observed trends are in agreement with previous reports for direct emissions apparent in MoS2 ML 47,56 and probably result from the competition between the efficiencies of the radiative and non-radiative recombination channels. Surprisingly, the temperature dependencies of the corresponding indirect transition intensities are clearly non-monotonic. First, quick reductions of the related intensities is observed from 5 K to about 100–150 K, which can be associated with increase of kinetic energies of excitons and interplay between radiative and non-radiative processes. Then it is followed by slower increase of the intensities up to room temperature, which can be explained in terms of the increased population of phonons at higher temperatures.SummaryA systematic investigation of optical properties of hBN-encapsulated artificially stacked MoS2 BLs with vary-ing interlayer twist angle has been conducted. It has been shown that RS spectroscopy is a sensible tool for determining the interlayer coupling and spacing in the BL. The strongest coupling has been observed in MoS2 homobilayers with twist angles of 6 ◦ and 62◦ , which reproduce well the structure of 2H and 3R polytypes found in natural and CVD-grown MoS2 multilayers. Increasing the twist angle of homobilayers from 6 ◦ to ∼ 30◦ leads to an increase in the interlayer spacing, and thus a decrease in the interlayer coupling. This effect can be observed as the absence of the low-energy interlayer phonon modes in the RS spectrum as well as the energy shift of ∼ 140 meV for the indirect transition between K-Ŵ points of the Brillouin zone. The PL of artificially assembled homobilayers, in comparison with the natural 2 L MoS2 , features a single line associated with indirect emission, but its energy varies considerably depending on the interlayer twist. We conclude that gaining a new degree of freedom by inter-layer twisting of artificially assembled flakes permits the control over the energy of the indirect transition, This offers further insight in few-layer 2D systems and can be a useful tool for future device applications.MethodsThe investigated MoS2 BLs, MLs, and hBN flakes were fabricated by two-stage PDMS-based mechanical exfolia-tion of bulk crystals. An unoxidized Si wafer was used as a substrate. In order to ensure the best quality of the substrate surface, they were annealed at 200 ◦ C and kept on hot-plate until the first non-deterministic transfer of h-BN flakes. Subsequent layers were transferred deterministically, to reduce inhomogeneities between each transfer sample were annealed. The complete structures were annealed at 160 ◦ C for 1.5 h to ensure the best layer-to-layer and layer-to-substrate adhesion and to eliminate a substantial portion of air pockets on the interfaces between the constituent layers.The SHG measurements performed to identify the relative twist angle of the exfoliated S-TMD MLs form-ing homobilayers were taken at T = 300 K using a home-built setup with a femtosecond Ti:Sapphire laser with excitation at 800 nm (1.55 eV). For each measurement, the laser light had a typical incident power of 500 µ W, was Figure 6.   Temperature evolution of the integrated intensities of the (solid circular points) direct- and (open diamond points) indirect-related emissions extracted from PL spectra measured on MoS2 : (a) bilayer (2 L) and selected homobilayers with twist angles of (b) 62◦ , (c) 18◦ and (d) 6 ◦ . Note that the vertical scales are logarithmic.7Vol.:(0123456789)Scientific Reports |        (2021) 11:17037  | https://doi.org/10.1038/s41598-021-95700-5www.nature.com/scientificreports/linearly polarized, and focused to a spot size of 1 µ m by a 50x objective lens. A set of a motorized half-wave plate and a fixed linear polarizer were used to analyse an SHG signal, which was detected by a Si avalanche photodiode.The PL and RS measurements were performed using � = 515 nm (2.41 eV) radiation from continuous-wave Ar-ion or diode lasers. For the PL and RS experiments at room temperature ( T = 300 K), the excitation light was focused through a 100x long-working distance objective with a 0.55 numerical aperture (NA) producing a spot of about 1 µ m diameter. The signal was collected via the same microscope objective, sent through a 1 m monochromator, and then detected by using a liquid nitrogen-cooled charge-coupled device (CCD) camera. To detect low-energy RS up to about ±10 cm−1 from the laser line, a set of Bragg filters was implemented in both excitation and detection paths. The temperature-dependent PL measurements were performed using an analogous setup with small modifications (a 50× long-working distance objective and a 0.5 m monochroma-tor). Moreover, the studied samples were placed on a cold finger in a continuous flow cryostat mounted on x–y motorized positioners. The excitation power focused on the sample was kept at 200 µ W during all measurements to avoid local heating.Data availibilityThe datasets obtained during experiments and analysis in course of manuscript preparation are available from the corresponding author on reasonable request.Received: 27 May 2021; Accepted: 27 July 2021References  1.  Mak, K. F., Lee, C., Hone, J., Shan, J. & Heinz, T. F. Atomically thin MoS2 : A new direct-gap semiconductor. Phys. Rev. Lett. 105, 136805. https://​doi.​org/​10.​1103/​PhysR​evLett.​105.​136805 (2010).  2.  Arora, A. et al. Excitonic resonances in thin films of WSe2 : From monolayer to bulk material. Nanoscale 7, 10421. https://​doi.​org/​10.​1039/​C5NR0​1536G (2015).  3.  Arora, A., Nogajewski, K., Molas, M., Koperski, M. & Potemski, M. Exciton band structure in layered MoSe2 : From a monolayer to the bulk limit. Nanoscale 7, 20769. https://​doi.​org/​10.​1039/​C5NR0​6782K (2015).  4.  Molas, M. R. et al. The optical response of monolayer, few-layer and bulk tungsten disulfide. Nanoscale 9, 13128. https://​doi.​org/​10.​1039/​C7NR0​4672C (2017).  5.  Ciorciaro, L., Kroner, M., Watanabe, K., Taniguchi, T. & Imamoglu, A. Observation of magnetic proximity effect using resonant optical spectroscopy of an electrically tunable MoSe2 CrBr3 heterostructure. Phys. Rev. Lett. 124, 197401. https://​doi.​org/​10.​1103/​PhysR​evLett.​124.​197401 (2020).  6.  Sun, J. et al. Lateral 2d WSe2 p–n homojunction formed by efficient charge-carrier-type modulation for high-performance opto-electronics. Adv. Mater. 32, 1906499. https://​doi.​org/​10.​1002/​adma.​20190​6499 (2020).  7.  Shi, H. et al. Ultrafast electrochemical synthesis of defect-free In2 Se2 flakes for large-area optoelectronics. Adv. Mater. 32, 1907244. https://​doi.​org/​10.​1002/​adma.​20190​7244 (2020).  8.  Bandurin, D. A. et al. High electron mobility, quantum hall effect and anomalous optical response in atomically thin InSe. Nat. Nanotechnol. 12, 223–227. https://​doi.​org/​10.​1038/​nnano.​2016.​242 (2016).  9.  Britnell, L. et al. Resonant tunnelling and negative differential conductance in graphene transistors. Nat. Commun. 4, 1794. https://​doi.​org/​10.​1038/​ncomm​s2817 (2013). 10.  Koperski, M. et al. Single photon emitters in exfoliated WSe2 structures. Nat. Nanotechnol. 10, 503–506. https://​doi.​org/​10.​1038/​nnano.​2015.​67 (2015). 11.  Kern, J. et al. Single-photon emitters: Nanoscale positioning of single-photon emitters in atomically thin WSe2 (adv. mater. 33/2016). Adv. Mater. 28, 7032–7032. https://​doi.​org/​10.​1002/​adma.​20167​0228 (2016). 12.  Kumar, S., Kaczmarczyk, A. & Gerardot, B. D. Strain-induced spatial and spectral isolation of quantum emitters in mono- and bilayer WSe2 . Nano Lett. 15, 7567–7573. https://​doi.​org/​10.​1021/​acs.​nanol​ett.​5b033​12 (2015). 13.  Branny, A., Kumar, S., Proux, R. & Gerardot, B. D. Deterministic strain-induced arrays of quantum emitters in a two-dimensional semiconductor. Nat. Commun. 8, 15053. https://​doi.​org/​10.​1038/​ncomm​s15053 (2017). 14.  Hennighausen, Z. & Kar, S. Twistronics: A turning point in 2d quantum materials. Electron. Struct. 3(1), 014004 (2021). 15.  He, F. et al. Moiré patterns in 2d materials: A review. ACS Nanohttps://​doi.​org/​10.​1021/​acsna​no.​0c104​35 (2021) (PMID: 33769797). 16.  Ponomarenko, L. A. et al. Cloning of Dirac fermions in graphene superlattices. Nature 497, 594–597. https://​doi.​org/​10.​1038/​natur​e12187 (2013). 17.  Hunt, B. et al. Massive Dirac fermions and Hofstadter butterfly in a van der Waals heterostructure. Science 340, 1427–1430. https://​doi.​org/​10.​1126/​scien​ce.​12372​40 (2013). 18.  Rivera, P. et al. Observation of long-lived interlayer excitons in monolayer MoSe2 –WSe2 heterostructures. Nat. Commun. 6, 6242. https://​doi.​org/​10.​1038/​ncomm​s7242 (2015). 19.  Zhang, N. et al. Moiré intralayer excitons in a MoSe2/MoS2 heterostructure. Nano Lett. 18, 7651–7657. https://​doi.​org/​10.​1021/​acs.​nanol​ett.​8b032​66 (2018). 20.  Jin, C. et al. Observation of moiré excitons in WSe2/WS2 heterostructure superlattices. Nature 567, 76–80. https://​doi.​org/​10.​1038/​s41586-​019-​0976-y (2019). 21.  Tran, K. et al. Evidence for moiré excitons in van der Waals heterostructures. Nature 567, 71–75. https://​doi.​org/​10.​1038/​s41586-​019-​0975-z (2019). 22.  Seyler, K. L. et al. Signatures of moiré-trapped valley excitons in MoSe2/WSe2 heterobilayers. Nature 567, 66–70. https://​doi.​org/​10.​1038/​s41586-​019-​0957-1 (2019). 23.  Alexeev, E. M. et al. Resonantly hybridized excitons in moiré superlattices in van der Waals heterostructures. Nature 567, 81–86. https://​doi.​org/​10.​1038/​s41586-​019-​0986-9 (2019). 24.  Lee, C. et al. Anomalous lattice vibrations of single- and few-layer MoS2 . ACS Nano 4, 2695–2700. https://​doi.​org/​10.​1021/​nn100​3937 (2010). 25.  Gołasa, K. et al. Resonant Raman scattering in MoS2—From bulk to monolayer. Solid State Commun. 197, 53–56. https://​doi.​org/​10.​1016/j.​ssc.​2014.​08.​009 (2014). 26.  Yamamoto, M. et al. Strong enhancement of Raman scattering from a bulk-inactive vibrational mode in few-layer MoTe2 . ACS Nano 8, 3895–3903. https://​doi.​org/​10.​1021/​nn500​7607 (2014). 27.  Lui, C. H. et al. Observation of interlayer phonon modes in van der Waals heterostructures. Phys. Rev. B 91, 165403. https://​doi.​org/​10.​1103/​PhysR​evB.​91.​165403 (2015).https://doi.org/10.1103/PhysRevLett.105.136805https://doi.org/10.1039/C5NR01536Ghttps://doi.org/10.1039/C5NR01536Ghttps://doi.org/10.1039/C5NR06782Khttps://doi.org/10.1039/C7NR04672Chttps://doi.org/10.1039/C7NR04672Chttps://doi.org/10.1103/PhysRevLett.124.197401https://doi.org/10.1103/PhysRevLett.124.197401https://doi.org/10.1002/adma.201906499https://doi.org/10.1002/adma.201907244https://doi.org/10.1038/nnano.2016.242https://doi.org/10.1038/ncomms2817https://doi.org/10.1038/ncomms2817https://doi.org/10.1038/nnano.2015.67https://doi.org/10.1038/nnano.2015.67https://doi.org/10.1002/adma.201670228https://doi.org/10.1021/acs.nanolett.5b03312https://doi.org/10.1038/ncomms15053https://doi.org/10.1021/acsnano.0c10435https://doi.org/10.1038/nature12187https://doi.org/10.1038/nature12187https://doi.org/10.1126/science.1237240https://doi.org/10.1126/science.1237240https://doi.org/10.1038/ncomms7242https://doi.org/10.1021/acs.nanolett.8b03266https://doi.org/10.1021/acs.nanolett.8b03266https://doi.org/10.1038/s41586-019-0976-yhttps://doi.org/10.1038/s41586-019-0976-yhttps://doi.org/10.1038/s41586-019-0975-zhttps://doi.org/10.1038/s41586-019-0975-zhttps://doi.org/10.1038/s41586-019-0957-1https://doi.org/10.1038/s41586-019-0957-1https://doi.org/10.1038/s41586-019-0986-9https://doi.org/10.1021/nn1003937https://doi.org/10.1021/nn1003937https://doi.org/10.1016/j.ssc.2014.08.009https://doi.org/10.1016/j.ssc.2014.08.009https://doi.org/10.1021/nn5007607https://doi.org/10.1103/PhysRevB.91.165403https://doi.org/10.1103/PhysRevB.91.1654038Vol:.(1234567890)Scientific Reports |        (2021) 11:17037  | https://doi.org/10.1038/s41598-021-95700-5www.nature.com/scientificreports/ 28.  Grzeszczyk, M. et al. Raman scattering of few-layers MoTe2 . 2D Mater. 3, 025010. https://​doi.​org/​10.​1088/​2053-​1583/3/​2/​025010 (2016). 29.  Kipczak, L., Grzeszczyk, M., Olkowska-Pucko, K., Babiński, A. & Molas, M. R. The optical signature of few-layer ReSe2 . J. Appl. Phys. 128, 044302. https://​doi.​org/​10.​1063/5.​00152​89 (2020). 30.  Holler, J. et al. Low-frequency Raman scattering in WSe2—MoSe2 heterobilayers evidence for atomic reconstruction. Appl. Phys. Lett. 117, 013104. https://​doi.​org/​10.​1063/5.​00122​49 (2020). 31.  Grzeszczyk, M. et al. Breathing modes in few-layer MoTe2 activated by h-BN encapsulation. Appl. Phys. Lett. 116, 191601. https://​doi.​org/​10.​1063/1.​51280​48 (2020). 32.  Liu, K. et al. Evolution of interlayer coupling in twisted molybdenum disulfide bilayers. Nat. Commun. 5, 4966. https://​doi.​org/​10.​1038/​ncomm​s5966 (2014). 33.  van Baren, J. et al. Stacking-dependent interlayer phonons in 3R and 2H MoS2 . 2D Mater. 6, 025022. https://​doi.​org/​10.​1088/​2053-​1583/​ab0196 (2019). 34.  Hsu, W.-T. et al. Second harmonic generation from artificially stacked transition metal dichalcogenide twisted bilayers. ACS Nano 8, 2951–2958. https://​doi.​org/​10.​1021/​nn500​228r (2014). 35.  Zhang, X. et al. Raman spectroscopy of shear and layer breathing modes in multilayer MoS2 . Phys. Rev. B 87, 115413. https://​doi.​org/​10.​1103/​PhysR​evB.​87.​115413 (2013). 36.  Zhao, Y. et al. Interlayer breathing and shear modes in few-trilayer MoS2 and WSe2 . Nano Lett. 13, 1007–1015. https://​doi.​org/​10.​1021/​nl304​169w (2013). 37.  Puretzky, A. A. et al. Low-frequency Raman fingerprints of two-dimensional metal dichalcogenide layer stacking configurations. ACS Nano 9, 6333–6342. https://​doi.​org/​10.​1021/​acsna​no.​5b018​84 (2015). 38.  Lee, J.-U. et al. Raman signatures of polytypism in molybdenum disulfide. ACS Nano 10, 1948–1953. https://​doi.​org/​10.​1021/​acsna​no.​5b058​31 (2016). 39.  van der Zande, A. M. et al. Tailoring the electronic structure in bilayer molybdenum disulfide via interlayer twist. Nano Lett. 14, 3869. https://​doi.​org/​10.​1021/​nl501​077m (2014). 40.  Huang, S. et al. Low-frequency interlayer Raman modes to probe interface of twisted bilayer MoS2 . Nano Lett. 16, 1435. https://​doi.​org/​10.​1021/​acs.​nanol​ett.​5b050​15 (2016). 41.  Lin, M.-L. et al. Moiré phonons in twisted bilayer MoS2 . ACS Nano 12, 8770–8780. https://​doi.​org/​10.​1021/​acsna​no.​8b050​06 (2018). 42.  Liao, M. et al. Precise control of the interlayer twist angle in large scale MoS2 homostructures. Nat. Commun. 11, 2153. https://​doi.​org/​10.​1038/​s41467-​020-​16056-4 (2020). 43.  Du, L. et al. Temperature-driven evolution of critical points, interlayer coupling, and layer polarization in bilayer MoS2 . Phys. Rev. B 97, 165410. https://​doi.​org/​10.​1103/​physr​evb.​97.​165410 (2018). 44.  Cadiz, F. et al. Well separated trion and neutral excitons on superacid treated MoS2 monolayers. Appl. Phys. Lett. 108, 251106. https://​doi.​org/​10.​1063/1.​49548​37 (2016). 45.  Cadiz, F. et al. Excitonic linewidth approaching the homogeneous limit in MoS2-based van der Waals heterostructures. Phys. Rev. X 7, 021026. https://​doi.​org/​10.​1103/​PhysR​evX.7.​021026 (2017). 46.  Molas, M. R. et al. Tuning carrier concentration in a superacid treated MoS2 monolayer. Sci. Rep. 9, 1989. https://​doi.​org/​10.​1038/​s41598-​018-​38413-6 (2019). 47.  Tongay, S. et al. Thermally driven crossover from indirect toward direct bandgap in 2d semiconductors: MoSe2 versus MoS2 . Nano Lett. 12, 5576–5580. https://​doi.​org/​10.​1021/​nl302​584w (2012). 48.  Koperski, M. et al. Optical properties of atomically thin transition metal dichalcogenides: Observations and puzzles. Nanophotonics 6, 1289–1308. https://​doi.​org/​10.​1515/​nanoph-​2016-​0165 (2017). 49.  Arora, A. et al. Interlayer excitons in a bulk van der Waals semiconductor. Nat. Commun. 8, 639. https://​doi.​org/​10.​1038/​s41467-​017-​00691-5 (2017). 50.  Arora, A. et al. Valley-contrasting optics of interlayer excitons in Mo- and W-based bulk transition metal dichalcogenides. Nanoscale 10, 15571–15577. https://​doi.​org/​10.​1039/​C8NR0​3764G (2018). 51.  Slobodeniuk, A. O. et al. Fine structure of k-excitons in multilayers of transition metal dichalcogenides. 2D Mater. 6, 025026. https://​doi.​org/​10.​1088/​2053-​1583/​ab0776 (2019). 52.  Gerber, I. C. et al. Interlayer excitons in bilayer MoS2 with strong oscillator strength up to room temperature. Phys. Rev. B 99, 035443. https://​doi.​org/​10.​1103/​physr​evb.​99.​035443 (2019). 53.  Paradisanos, I. et al. Controlling interlayer excitons in MoS2 layers grown by chemical vapor deposition. Nat. Commun. 11, 2391. https://​doi.​org/​10.​1038/​s41467-​020-​16023-z (2020). 54.  Kormányos, A. et al.k · p theory for two-dimensional transition metal dichalcogenide semiconductors. 2D Mater. 2, 022001, https://​doi.​org/​10.​1088/​2053-​1583/2/​2/​022001 (2015). 55.  ODonnell, K. P. & Chen, X. Temperature dependence of semiconductor band gaps. Appl. Phys. Lett. 58, 2924–2926. https://​doi.​org/​10.​1063/1.​104723 (1991). 56.  Wang, G. et al. Spin-orbit engineering in transition metal dichalcogenide alloy monolayers. Nat. Commun. 6, 10110. https://​doi.​org/​10.​1038/​ncomm​s10110 (2015).AcknowledgementsThe work has been supported by the National Science Centre, Poland (Grant Nos. 2017/27/B/ST3/00205, 2017/27/N/ST3/01612, 2018/31/B/ST3/02111), EU Graphene Flagship project (No. 785219), the ATOMOPTO project (TEAM programme of the Foundation for Polish Science, co-financed by the EU within the ERD-Fund), and the CNRS via IRP “2DM” project. The Polish participation in EMFL is supported by the DIR/WK/2018/07 grant from Polish Ministry of Science and Higher Education. K. W. and T. T. acknowledge support from the Elemental Strategy Initiative conducted by the MEXT, Japan, (Grant No. JPMXP0112101001), JSPS KAKENHI (Grant No. JP20H00354), and the CREST (JPMJCR15F3), JST.Author contributionsM.G., J.S., M.B., and M.R.M carried out the RS and PL experiments. A.O.S. developed the theory. T.K., P.K., and M.R.M. performed the SHG experiments. K.W. and T.T. grew the hBN crystals. M.G. fabricated the investigated samples. M.G. and M.R.M. performed the data analysis. M.G., M.P., A.B., and M.R.M. wrote the manuscript with contribution from all other co-authors.Competing interest The authors declare no competing interests.https://doi.org/10.1088/2053-1583/3/2/025010https://doi.org/10.1063/5.0015289https://doi.org/10.1063/5.0012249https://doi.org/10.1063/1.5128048https://doi.org/10.1063/1.5128048https://doi.org/10.1038/ncomms5966https://doi.org/10.1038/ncomms5966https://doi.org/10.1088/2053-1583/ab0196https://doi.org/10.1088/2053-1583/ab0196https://doi.org/10.1021/nn500228rhttps://doi.org/10.1103/PhysRevB.87.115413https://doi.org/10.1103/PhysRevB.87.115413https://doi.org/10.1021/nl304169whttps://doi.org/10.1021/nl304169whttps://doi.org/10.1021/acsnano.5b01884https://doi.org/10.1021/acsnano.5b05831https://doi.org/10.1021/acsnano.5b05831https://doi.org/10.1021/nl501077mhttps://doi.org/10.1021/acs.nanolett.5b05015https://doi.org/10.1021/acs.nanolett.5b05015https://doi.org/10.1021/acsnano.8b05006https://doi.org/10.1038/s41467-020-16056-4https://doi.org/10.1038/s41467-020-16056-4https://doi.org/10.1103/physrevb.97.165410https://doi.org/10.1063/1.4954837https://doi.org/10.1103/PhysRevX.7.021026https://doi.org/10.1038/s41598-018-38413-6https://doi.org/10.1038/s41598-018-38413-6https://doi.org/10.1021/nl302584whttps://doi.org/10.1515/nanoph-2016-0165https://doi.org/10.1038/s41467-017-00691-5https://doi.org/10.1038/s41467-017-00691-5https://doi.org/10.1039/C8NR03764Ghttps://doi.org/10.1088/2053-1583/ab0776https://doi.org/10.1103/physrevb.99.035443https://doi.org/10.1038/s41467-020-16023-zhttps://doi.org/10.1088/2053-1583/2/2/022001https://doi.org/10.1088/2053-1583/2/2/022001https://doi.org/10.1063/1.104723https://doi.org/10.1063/1.104723https://doi.org/10.1038/ncomms10110https://doi.org/10.1038/ncomms101109Vol.:(0123456789)Scientific Reports |        (2021) 11:17037  | https://doi.org/10.1038/s41598-021-95700-5www.nature.com/scientificreports/Additional informationSupplementary Information The online version contains supplementary material available at https://​doi.​org/​10.​1038/​s41598-​021-​95700-5.Correspondence and requests for materials should be addressed to M.G. or M.R.M.Reprints and permissions information is available at www.nature.com/reprints.Publisher’s note  Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.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://​creat​iveco​mmons.​org/​licen​ses/​by/4.​0/.© The Author(s) 2021https://doi.org/10.1038/s41598-021-95700-5https://doi.org/10.1038/s41598-021-95700-5www.nature.com/reprintshttp://creativecommons.org/licenses/by/4.0/ The optical response of artificially twisted MoS bilayers Results Atomic structures of twisted bilayers MoS.  Raman scattering spectroscopy.  Photoluminescence experiment.  Summary Methods References Acknowledgements