Raman scattering from the bulk inactive out–of–plane B2g1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\bf{B}}}_{{\bf{2}}{\bf{g}}}^{{\bf{1}}}$$\end{document} mode in few–layer MoTe2

We report a study of Raman scattering in few-layer MoTe2 focused on high-frequency out-of-plane vibrational modes near 291 cm−1 which are associated with the bulk-inactive B2g1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\rm{B}}}_{2{\rm{g}}}^{1}$$\end{document} mode. Our temperature-dependent measurements reveal a double peak structure of the feature related to these modes in the Raman scattering spectra of 4- and 5-layer MoTe2. In accordance with literature data, the doublet’s lower- and higher-energy components are ascribed to the Raman-active A1g/A′1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\bf{A}}{\boldsymbol{^{\prime} }}}_{{\bf{1}}}$$\end{document} vibrations involving, respectively, only the inner and surface layers. We demonstrate a strong enhancement of the inner mode’s intensity at low temperature for 1.91 eV and 1.96 eV laser light excitation which suggests a resonant character of the Raman scattering processes probed under such conditions. A resonance of the laser light with a singularity of the electronic density of states at the M point of the MoTe2 Brillouin zone is proposed to be responsible for the observed effects.

Almost a decade ago layered transition metal dichalcogenides (TMDs), known in solid-state physics since the late 1960s 1 , were rediscovered as a vast group of materials which in the mono-and few-layer form can complement graphene in the quest for alternatives to electronics based on silicon and other conventional semiconductors. Since publication of seminal works of K. F. Mak et al. 2 and A. Splendiani et al. 3 dedicated to molybdenum disulfide, a huge amount of effort has been put into understanding physical and chemical properties of different members of the TMD family. When subject to a more detailed study, thin-film TMDs turned out to offer a number of fascinating possibilities not only in the field of electronics 4,5 , but also optoelectronics 6 , photovoltaics 7,8 , and thermoelectricity 9 , just to mention a few domains of the utmost importance from the research and development point of view. Moreover, they proved themselves to constitute an interesting platform for fundamental investigations. Indeed, due to specific crystal structure, the energy bands and resulting properties of TMDs, composed of layers of covalently bound metal and chalcogen atoms which are kept together by relatively weak van der Waals forces, can be conveniently tuned by simply changing the number of layers (N). This is how, in particular, the band gap of semiconducting TMDs is transformed from indirect in the three-dimensional crystals to direct in the monolayer limit. Along with light emission, a strong thickness dependence in TMD materials is exhibited also by their lattice dynamics, commonly probed with the use of Raman scattering (RS) spectroscopy 10 . A pronounced influence of the number of layers on the RS spectrum is observed for both low-energy shear and breathing modes 11,12 and higher-energy modes of in-and out-of-plane oscillations [13][14][15][16] . The impact of sample's thickness on the lattice vibrations becomes even stronger when resonant excitation of RS is employed. Under such conditions, the resulting change in the electronic configuration of investigated TMD material leads to several effects which cannot be studied with non-resonant RS measurements, such as the enhancement of particular features in the RS spectrum 17 related to multiphonon processes 18 or quenching of Raman-active modes due to quantum interference 19,20 . The observation of Davydov splitting of phonon modes in few-layer TMD films, which arises from interlayer interactions, is also facilitated by resonant excitation of RS spectra 21,22 . It has been reported for all semiconducting few-layer TMD materials: MoS 2 23 , MoSe 2 24 , WS 2 25,26 , WSe 2 27 and MoTe 2 28-30 . Last but not least, resonant RS can also unveil those vibrational modes which are symmetry allowed, but barely visible under non-resonant excitation. A prominent example of such "brightening" is the observation of high-frequency out-of-plane modes of A 1g / ′ A 1 symmetry, which in N-layer TMDs originate from the bulk-inactive B 2g 1 mode 28,[31][32][33] . Surprisingly, in contrast to extensively studied low-and mid-frequency oscillations of few-layer TMD crystals, these modes have not received so far much attention. Strong enhancement of this phonon peak was first reported by Yamamoto et al. 31 , where the appearnce in the RS spectrum was atributed to translation symmetry breaking in thin MoTe 2 layers. The work of Froehlicher et al. 28 exposed the fine doublet structure of these out-of-plane modes at room temperature, which occurs due to surface effects. However, the splitting presented in previous works was rather weak and the additional features were not very distinct. This suggests the sensitivity of the studied system to the resonant conditions is high and should be appropriately adjusted for specific vibrations.
To this aim we employ various excitation wavelengths to measure the RS in few-layer MoTe 2 and we follow temperature evolution of the RS spectra in 2-to 5-layer-thick MoTe 2 studied for 1.91 eV and 1.96 eV laser light excitation. A doublet structure of the out-of-plane mode in question, which can be hardly noticed at room temperature, becomes clearly observable at low temperature in tetra-and pentalayer samples. This results from a substantial increase of its lower-energy component's intensity with decreasing temperature. The experimental results are discussed with respect to a previously reported temperature evolution of the A 1g / ′ A 1 mode in MoTe 2 19 . We conclude that the observed temperature evolution of the investigated out-of-plane mode is related to resonant character of the RS for the employed light excitation. We propose that the resonance of the laser light with a singularity of the electronic density of states at the M point of the Brilloiun zone in MoTe 2 might be responsible for the observed anomalous increase of the RS features.

Results
Molybdenum ditelluride (MoTe 2 ) is a semiconducting TMD with a relatively small band-gap energy close to 1.3 eV 34 which makes it well-suited for resonant RS study with the use of commonly available light sources like He-Ne (1.96 eV) or Kr-Ar (1.91 eV) lasers. The former laser energy coincides, in particular, with substantial density of states in few-layer MoTe 2 , which occurs at the M and K points of the Brillouin zone 32 . Similar to other semiconducting tungsten-and molybdenum-based TMDs, MoTe 2 crystallizes in a trigonal prismatic (2H) structure, which symmetries in the limit of bulk material belong to the D 6h point group. As a result, with six atoms in the unit cell, multilayer MoTe 2 possesses 18 modes of normal vibrations out of which 7 are Raman-active (doubly degenerate in-plane: E 1g , E 2g 1 and E 2g 1 , and non-degenerate out-of-plane A 1g ), 6 are infrared (IR)-active (doubly degenerate in-plane E 1u 1 and E 1u 2 , and non-degenerate out-of-plane A 2u 1 and A 2u 2 ), and 5 are optically inactive (doubly degenerate in-plane E 2u , and non-degenerate out-of-plane: B 1g , B 2g 1 and B 1u ) 35 . On the other hand, the point groups describing symmetries of few-layer MoTe 2 composed of even and odd number of layers are D 3d and D 3h , respectively. This leads to different symmetry labelling of vibrational modes depending on the parity of N. In a structure with odd (even) N, the Raman-active modes are denoted by E′, E″ and ′ A 1 (E g and A 1g ). The difference in the labelling is presented in Fig. 1(a) for the group of the out-of-plane B 2g 1 -related modes. Figure 1(b) shows microscopic images of the thin MoTe 2 layers under study. All phonon modes observed in the RS spectra of few-layer MoTe 2 are presented in Fig. 2, which demonstrates the results obtained for monolayer (1 L) to pentalayer (5 L) excited with several different laser wavelenghts. It can be seen that both the relative intensities and spectral lineshapes of individual peaks for a particular thickness strongly depend on the excitation energy.
Let us focus first on the peaks due to the out-of-plane A 1g / ′ A 1 vibrations at ~174 cm −1 , which correspond to the Raman-active A 1g mode in bulk MoTe 2 28-30 . Their characteristic pattern originating from the Davydov splitting can be clearly observed in 3 L, 4 L and 5 L MoTe 2 for 1.98 eV, 1.96 eV or 1.91 eV, while it remains hardly distinguishable for 1.58 eV excitation (see Fig. 2). The relative intensities of the lower-and higher-energy components of the ′ A 1 mode in 3 L MoTe 2 also depend on the excitation energy. The higher-energy component of the ′ A 1 mode in 3 L MoTe 2 , which is due to vibrations occuring with the same phase in all three covalently bound MoTe 2 layers is of higher intensity than its lower-energy counterpart in the RS spectrum excited with 2.41 eV light. On the contrary, the lower-energy component due to the mode in which the tellurium atoms in the middle MoTe 2 layer vibrate out-of-phase as compared to vibrations in the outer layers, is of higher intensity in the RS spectrum excited with 1.96 eV or 1.91 eV light. Finally, the Davydov-split components of the ′ A 1 mode in 3 L MoTe 2 are of similar intensity in the RS spectrum excited with 1.98 eV light. As compared to the in-plane mode, the relative intensities of spectral features originating from the out-of-plane mode also strongly depend on the excitation energy. For 1 L MoTe 2 the out-of-plane mode dominates the RS spectrum excited with 1.98 eV, 1.96 eV or 1.91 eV light, whereas the A 1g / ′ A 1 peak is much weaker than the one due to in-plane vibrations corresponding to the E 2g 1 mode in the bulk (observed at ~235 cm −1 ). For 5 L MoTe 2 the out-of-plane mode is weaker than its in-plane counterpart for 2.41 eV and 2.33 eV excitation while the opposite happens when 1.58 eV excitation is employed. This behaviour indicates a resonant RS of light on the out-of-plane modes, which occurs when the laser energy coincides with a maximum of electronic density of states in the crystal lattice 32 . This effect was previously reported for MoTe 2 29,30,36 and explained in terms of quantum interference of contrubutions from different points of the Brillouin zone and the electron-phonon coupling 20 . A possible role of resonant vs non-resonant contributions to the RS efficency was also suggested in ref. 19 .
A focal point of this paper is the feature apparent at ~291 cm −1 in the RS spectra of flakes composed of at least two MoTe 2 layers. Because of symmetries identifing this mode it is labelled as A 1g / ′ A 1 for few-layer thick structures. The Mulliken notation is exactly the same as for the peaks at ~174 cm −1 described in the preceding paragraph. Therefore, in order to avoid any confusion, unless otherwise stated, the A 1g / ′ A 1 label will always refer in what follows to the feature at ~291 cm −1 . It corresponds to vibrations in which the tellurium atoms in each layer oscillate out-of-phase with respect to the molybdenum atoms (see Fig. 1). There are two such modes in bulk MoTe 2 . The atoms in neighboring layers vibrate out-of-phase in the optically inactive B 2g 1 mode while in the infrared-active (A 2u ) mode the in-phase vibrations occur. As a result none of the modes contribute to the RS in bulk MoTe 2 . Raman-inactive (IR-active) is also the corresponding ″ A 2 mode in monolayer MoTe 2 . Additional symmetry elements, like the inversion symmetry center for even or the mirror plane for odd number of layers, activate the phonon mode in the material consisting of at least two MoTe 2 layers. Two out-of-plane modes in 2 L MoTe 2 are similar to those present in the bulk. However, due to the symmetry, the A 1g mode, in which the atoms in both layers move out-of-phase with respect to each other becomes Raman-active (see Fig. 1). The other mode, A 2u , in which the vibrations in both layers occur in-phase, is IR-active. A similar line of reasoning may be applied to any other N-layer structure in which N is an even number. The number of Raman-active A 1g modes and the number of IR-active A 2u modes in such cases equals N 2 28 . If N has an odd value, the Raman-active modes must hold the mirror symmetry with respect to the plane crossing the Mo atoms in the central layer (see Fig. 1). The number of these ′ A 1 modes is equal to − N 1 2 and the number of corresponding IR-active ″ A 2 modes amounts to + N 1 2 28 . The above considerations lead, in particular, to a conclusion that there are: one Raman-active ′ A 1 mode and two ″ A 2 IR-active modes in 3 L MoTe 2 , two Raman-active A 1g modes and two A 2u IR-active modes in 4 L MoTe 2 , and two Raman-active ′ A 1 modes and three IR-active ″ A 2 modes in 5 L MoTe 2 . Corresponding normal mode displacements are schematically presented in Fig. 1. Two different modes of Raman-active vibrations in 4 L and 5 L samples are referred to as the surface mode (s) and inner mode (i) 28,33 . The atomic displacements in the corresponding normal modes occur in the outside layers (s) and in central layers (i), respectively. Based on theoretical calculations 28 the energy splitting between the modes has been estimated to be ~1.2 cm −1 . Previous works attributed the splitting mainly to surface effects. The Davydov splitting, which is due to van der Waals interactions between neighboring layers (being of the order of ~0.1 cm −1 ) is negligible in this case. Experimental study carried out using 1.96 eV excitation show that the inner mode appears in the RS spectra of N ≥ 4 L MoTe 2 as a faint shoulder on the peak related to the surface mode 28 . Our results confirm that observation. However, the lower-energy peak related to the bulk inactive B 2g 1 mode can be observed in the spectra measured at room temperature just as a weak feature (see Fig. 3). In order to improve its visibility by changing the excitation conditions, we have studied the RS as a function of temperature. It is well known that the semiconductor energy band structure changes with temperature with the band gap usually widening when the temperature is lowered. Temperature modulation therefore allows to form variable resonance conditions without changing the laser excitation energy 17,19 . The effect of temperature on the investigated RS can be appreciated in Figs 3 and 4 showing a series of RS spectra excited with 1.96 eV and 1.91 eV light and measured as a function of temperature ranging from 5 K to 300 K on 2-5 L MoTe 2 .
It can be clearly seen in Figs 3 and 4 that the lineshape of the investigated spectral features substantially changes with decreasing temperature. Let us start the analysis of the results with the most staggering effects observed for 4 L and 5 L MoTe 2 . For both excitation energies, the lower-energy component of the doublet significantly gains the intensity with decreasing temperature. The effect can be appreciated for both MoTe 2 structures, however, at T = 5 K it is stronger for 4 L than for 5 L flakes. It is also stronger for 1.91 eV than for 1.96 eV excitation. At T = 5 K the lower-energy component of the A 1g / ′ A 1 modes dominates the high-energy part of the RS spectrum excited with 1.91 eV. Following the theoretical disscusion presented in ref. 28 we attribute the higher-energy component of the investigated peak in 4 L and 5 L to the surface mode (s) of the out-of-plane A 1g or ′ A 1 vibrations, respectively, which involve atomic dispacements in the outer MoTe 2 layers (see Fig. 1). The lower-energy component is the Raman-active A 1g / ′ A 1 inner mode (i) in 4 L/5 L MoTe 2 . The energy splitting between the components matches the value predicted by a linear chain model. This points at surface effects as the main factors which affect the splitting between the surface and the inner mode in the sample. An observation, which goes beyond the previously reported data, concerns the full width at half maximum (FWHM) of both components (see Table 1).
The lower-energy components of the investigated doublets are broader than their higher-energy counterparts. This may suggest that the former features correspond to more than one vibrational mode. Additional vibrations which may contribute to the lower energy peak are most likely the IR-active modes. They are expected to be quasi-resonant with the inner modes. Although they should not appear in the RS spectrum, it is known that resonant excitation and/or disorder in the structure can make them Raman-active 25,33,37 . This conclusion is supported by the RS from 2 L and 3 L MoTe 2 (see Figs 3 and 4). In both cases there is just one Raman-active mode of the out-of-plane vibrations, which is actually the surface mode. However, at low temperature, the lower-energy peak related to the out-of-plane vibrations becomes apparent in the spectrum. Its intensity is lower than the intensity of its higher-energy counterpart but its presence in the RS spectra remains beyond any doubts. As in the case of 4 L and 5 L MoTe 2 , also for 2 L and 3 L MoTe 2 the peaks can be seen more distinctively for 1.91 eV excitation.   In order to highlight the resonant character of changes in the RS spectra of few-layer MoTe 2 described above, the best way is to plot the intensities of individual spectral features as a function of temperature. For better comparison, it is convenient to relate the intensities of the out-of-plane modes to the intensity of the in-plane E g /E′ mode, as this mode is not expected to exhibit any resonant behaviour (see the grey points in Fig. 5). Such an approach accounts for small variations of the excitation laser intensity 17 .
The temperature dependences of the relative intensities of the A 1g / ′ A 1 related peaks observed in the RS spectra of 2-5 L MoTe 2 samples under study are presented in Fig. 5. As previously stated, the most pronunced is the effect of temperature on the RS spectrum excited with 1.91 eV light. The peak due to the A 1g (s) mode in 4 L MoTe 2 looses its intensity with decreasing temperature. This behaviour is opposite to the temperature evolution of the peak related to the A 1g (i) mode (with a possible contribution from the IR-active out-of-plane mode). This peak is almost seven times more intense at T = 5 K than at room temperature. A similar effect of temperature can be noticed for spectral features of 5 L MoTe 2 . In this case, the ′ A 1 (i) mode is three times stronger at T = 5 K than at room temperature. We interpret the observed effect of temperature on the out-of-plane vibrational modes in terms of the transformation of the MoTe 2 band structure induced by changing the temperature, which we discuss in more detail in the following section.

Discussion
It has previously been noted that the resonance of the exciting light with the electronic transitions in a semiconductor structure can result in strong enhancement of particular RS modes 38 . This effect was observed in several semiconductor TMDs including bulk MoS 2 17 . It was shown that modification of the MoS 2 band structure caused by either the temperature or pressure can lead to significant increase of the intensity of the RS peak due to the out-of-plane Raman-active A 1g mode. We propose a similar explanation for the observed enhancement of the A 1g (i) mode in 4 L MoTe 2 reported in this work.
In our case, however, a two-dimensional space of tuning parameters which affect the MoTe 2 band structure is spanned by temperature (T) and the number of layers in the sample (N). Symbolically, it can be expressed by a formula E i (k, T, N), where E stands for the energy, i is the band index, and k indicates the wave vector. Following this formalism we can say that the results shown in the bottom panel of Fig. 5 clearly demonstrate that the excitation energy of 1.91 eV employed in our experiment corresponds to such a point in the E i (k, T, N) space, where in the vicinity of a particular k, T = 5 K and N = 4 a substantial density of electronic states exists in the MoTe 2 bands linked by optically active transitions. As also shown in Fig. 5, tuning E, N, or T by as little as 50 meV, plus/minus one layer, and 100 K, respectively, while keeping the other two parameters constant, can efficiently quench the intensity of the A 1g -mode-related peak in the RS spectrum of few-layer MoTe 2 . It shows how sensitive the resonant RS is as a tool to probe the local density of electronic states in a semiconductor system.
We suggest that the strong enhancement of the A 1g (i) mode's intensity in the RS spectrum of 4 L MoTe 2 is associated with electronic transitions taking place at the so-called M point of the Brillouin zone of the hexagonal crystal lattice. Indeed, as demonstrated by transmission measurements performed on bulk 1,39 and few-layer 2H-MoTe 2 30,34 , in the vicinity of 1.91 eV there exists a substantial density of electronic states which gives rise to a doublet of excitonic transitions known as A′ and B′. Their labelling intentionally resembles that of the optical transitions at the K points of the Brillouin zone as in both cases the local minima and maxima of the conduction and valence band, respectively, are split into pairs of subbands by spin-orbit interaction. Essentially independent of the sample thickness, the A′ transition occurs at room temperature at about 1.73 eV whereas the B′ is observed at about 1.98 eV. With decreasing temperature, due to thermal shrinkage of the crystal lattice and resulting effective repulsion of the energy bands, the energies of all excitonic transitions experience a blueshift. For the lowest-energy K-point transition between the spin-orbit-split subbands in the conduction and valence bands of 1 L MoTe 2 it amounts to about 0.1 eV as inferred from the comparison of optical data recorded at 4.5 K 34 and room temperature 40 . Worth noting here is that the blueshift of fundamental band-gap transitions is about 0.1 eV when going from room to liquid helium temperature is characteristic of all monolayer TMDs, as shown by temperature-dependent reflectance contrast measurements performed on WSe 2 41 , MoSe 2 42 , and WS 2 43 . To the first approximation, the same temperature evolution as for the K points can be assumed also for transitions taking place at other points of the Brillouin zone, as the uniform change in the equilibrium distance between atoms in the MoTe 2 crystal lattice should have a homogeneous impact on the whole band structure. We should then expect the A′ and B′ transitions to occur at 4.5 K at about 1.83 eV and 2.08 eV. Although it is difficult to precisely estimate the broadening of these transitions based on the experimental results, presented in refs 1,30,34,39 , it must amount to at least a few tens of meV. It means that with 1.91 eV excitation we most probably probed the high-energy shoulder of the maximum in the local density of electronic states which gives rise to the A′ transition. This conjecture also explains the sensitivity of our results to the excitation energy since changes in the resonance conditions should be much more dynamic for the slope than for the close-to-maximum part of the peak in the local density of states. In this respect it would be desirable to extend our study to excitation energies falling in the range from about 1.75 eV to 1.91 eV.
A question arises why only the inner modes are in resonance with 1.91 eV laser excitation, while the surface modes react to it relatively weakly. Definitely, a more formal theoretical approach is necessary to account for this observation. We suppose, however, that it may be related to the effect of particular electronic excitation on the polarizability of bonds involved in different vibrational modes. For example, as it was demonstrated in ref. 20 , for 3 L MoTe 2 the contributions from different regions of the Brillouin zone add with particular signs (plus or minus). Consequently, this may lead to a partial quenching of some features observed in the RS spectrum. This effect can be easily recognized for the high-energy component of the out-of-plane ′ A 1 vibrations in MoTe 2 , where the contributions from the vicinity of the M point partially cancel each other out for E L = 1.96 eV excitation 19,20 . This is not the case of the low-energy part of the mode. We believe that a similar effect is responsible for our observations. Temperature modulation of the MoTe 2 band structure modifies resonant excitation conditions and may result in different effects on individual vibrational modes.

Conclusions
In conclusion, we have investigated Raman scattering from the out-of-plane vibrational modes (A 1g / ′ A 1 ) in few-layer MoTe 2 of thickness ranging from one to five layers. Our temperature-dependent measurements performed with the use of several excitation energies clearly demonstrate a doublet structure of one of the features corresponding to these modes in the RS spectrum at about 291 cm −1 , especially pronounced in 4 L-and 5 L-thick samples excited with 1.91 eV and 1.96 eV laser light. The lower-and higher-energy components of the doublet have been identified as the inner and surface modes of the out-of-plane oscillations, respectively. We have shown that for 1.91 eV excitation a strong enhancement of the inner mode's contribution to the RS spectrum is observed for 4 L and 5 L MoTe 2 . Based on temperature evolution of this effect, we qualitatively interpret it as being associated with the resonance between the A 1g (i) mode and the A′ excitonic transition at the M point of the MoTe 2 Brillouin zone, which at 4.5 K is expected to occur at about 1.83 eV. Taking into account the finite broadening of the A′ transition (at the level of few tens of meV), the enhancement of the A 1g (i)-mode-related feature's intensity in the RS spectrum we have demonstrated for 1.91 eV excitation can be understood as originating from the onset of resonant conditions (from the high-energy side) reaching their maximum at slightly lower energies, not yet available in the experimental set-up exploited for the measurements reported in this paper. A striking observation about the resonance mentioned above is its extreme sensitivity to the sample thickness which cannot result from the MoTe 2 band structure, which at the M point of the Brillouin zone barely changes with the number of layers in the sample. It may hint at important details we have not considered in our simplified interpretation, like, for instance, the thickness and temperature dependence of electron-phonon interactions. Definitely, a much more strict and formal analysis of our experimental data is necessary to fully account for all the results reported in this manuscript, not only qualitatively but also at the quantitative level. We do hope that the present work will stimulate a theoretical study which will shed more light on many aspects of resonant Raman scattering in TMD systems which have not yet been thoroughly understood.

Methods
MoTe 2 samples were prepared by an all-dry PDMS-based exfoliation technique 44 of a bulk crystal purchased from HQ graphene and deposited onto a Si/(90 nm) SiO 2 substrate. The thickness of studied MoTe 2 flakes was determined by optical contrast and Raman spectroscopy. The unpolarized RS measurements were carried out in the backscattering geometry. Raman spectra were taken at various laser excitation wavelenghts: 785 nm (1.58 eV), 625.4 nm (1.98 eV), 632.8 nm (1.96 eV), 647.1 nm (1.91 eV), 514.5 nm (2.41 eV), and 532 nm (2.33 eV). The Stokes scattering was measured on 2-5 L MoTe 2 in a temperature range from 5 K to 300 K. The investigated samples were placed on a cold finger of a continuous flow cryostat, mounted on the x-y motorized positioners. The excitation light was focused by means of a 50x magnification long working distance objective. The spot diameter of the focused beam was 1 μm. The excitation power focused on the sample was kept at 100 μW during all measurments to avoid local heating. The RS signal was collected via the same microscope objective, sent through a 0.75 m monochromator and detected with a liquid-nitrogen cooled CCD camera.

Data Availability
The datasets obtained during experiments and analysis in course of manuscript preparation are available from the corresponding author on reasonable request.