1. Introduction
Two-dimensional (2D) materials exhibit unique properties that depend on various factors, including the material’s chemical composition, atomic arrangement, thickness, and interlayer interactions [
1,
2]. Modification of the intrinsic properties of crystalline materials is often necessary to achieve new fundamental effects or create favorable conditions for device fabrication. By carefully tuning these parameters, researchers can unlock new properties and behaviors in solid-state materials, leading to innovative technologies and scientific discoveries [
2,
3,
4]. Recent studies have explored the effects of mechanical deformation and stress on the properties of 2D materials [
5,
6,
7,
8,
9]. In particular, the ability of 2D materials to be easily strained is of great interest, as strain engineering can significantly impact their electronic and optical properties [
8,
10,
11,
12,
13,
14,
15,
16]. Indeed, many theoretical works have predicted that biaxial strain is particularly effective in tuning the band structure of transition metal dichalcogenides (TMDs) [
17,
18,
19,
20,
21,
22].
Atomically thin materials of the MoS
2-type family are particularly well-suited for studying the effects of mechanical deformation and stress because they can withstand extreme nonhomogeneous deformations before rupture [
9,
23,
24,
25]. The electronic and optical properties of MoS
2 monolayers are strongly coupled to the valley/spin/orbital degrees of freedom and the lattice structure, making them sensitive to mechanical deformation or stress [
5,
11,
25,
26,
27,
28]. To explore the effect of biaxial strain on TMDs, various methods have been developed, such as deposition on nanocones [
29] or pillars [
30] and epitaxial growth of superlattices [
31]. In these methods, the strain obtained is typically around 1-2%. However, values of strain of 5% are reached at the center of MoS
2 monolayer domes that are investigated in this work [
11,
32].
In this letter, we present a resonant Raman spectroscopy study of MoS
2 dome sample using 23 different laser excitation energies in the NIR and visible ranges. Our measurements allowed us to obtain the Raman excitation profiles (REPs) of the first-order E
’ and A
1’ Raman modes, as well as the REP of the second-order 2LA Raman band. Our results show that the three Raman bands are significantly enhanced at the indirect gap transition at 1.61 eV and by the C exciton at 2.72 eV, but the interaction of the first-order modes with the A and B excitons was shown to be very weak in the strained MoS
2 dome. It was also observed redshifts in the energies of the A and B excitons, in agreement with previous photoluminescence (PL) results [
11,
18,
21,
33], as well as in the energy of the C excitons, that was not yet reported in the literature. Our work provides important insights into the effects of biaxial strain on the excitonic and optical properties of 2D materials, which could have implications for the development of innovative technologies and scientific discoveries.
3. Results
We began our study by characterizing the Raman modes of a selected dome.
Figure 1 presents the optical image, the Raman spectra and intensity Raman maps acquired, and positions of A
’1 and E
’ modes from sample. The dashed circles around the dome delineate its edge.
Figure 1(a) shows that the semispherical dome has a diameter of approximately (2.4 ± 0.1)
m.
Figure 1(b) presents the Raman spectra recorded at the center of the dome (black line) and at the MoS
bulk substrate (red line) with 2.34 eV laser energy. The bulk’s Raman spectrum exhibits two modes, E
and A
1g, located around 384 cm
-1 and 409 cm
-1, respectively [
35,
36]. The Raman spectra at the dome’s center exhibit four modes, where the two weaker peaks come from the bulk’s substrate and the two more intense peaks, shifted to lower wavenumbers around 373 cm
-1 and 402 cm
-1, are associated with the E
’ and A
1’ modes of the monolayer dome, respectively. In the literature, the positions of the E
’ and A
1’ modes in unstrained monolayer MoS
are 385 cm
-1 and 404 cm
-1, respectively [
36,
37,
38]. Therefore, the redshifts in the mode positions of the strained monolayer MoS
studied in this work with respect to the unstrained monolayer MoS
are 12 cm
-1 and 2 cm
-1 for the E
’ and A
1’ modes, respectively [
38]. Previous studies have shown that the effect of the biaxial strain at the dome is more significant for the E
mode, with a displacement rate of 2.2 cm
/% [
5]. According to this relation, the shift of 12 cm
observed in our work for the E
’ mode corresponds to a value of ≈ 5.5% of strain at the dome’s center shown in
Figure 1(a).
Figures 1(c)-(f) show the intensity maps and positions of the dome peaks acquired with a laser energy of 2.34 eV. The intensity maps of the E’ and A’1 Raman peaks indicate an increase in the modes intensities at the dome’s region compared to the bulk spectra, as displayed in Figures 1(c) and 1(d), respectively. The positions of the E’ and A1’ Raman peaks at the monolayer dome’s center are red-shifted by about 11 cm-1 and 7 cm-1, respectively, when compared to the bulk’s positions, as shown in Figures 1(e) and 1(f), respectively.
Figure 2(a) shows the Raman spectra recorded at the dome’s center using five different laser excitation energies: 1.59 eV (780 nm), 1.68 eV (738 nm), 1.92 eV (646 nm), 2.38 eV (521 nm) and 2.73 eV (454 nm). The three dashed lines serve as guides to follow the positions of the in-plane E
’, the out-of-plane A
1’ and the 2LA modes. The Raman spectra of the bulk sample were subtracted from the spectra, so that it only displays the peaks of the dome’s center. In addition to the first order-modes, we can observe the second-order 2LA Raman bands [
39] centered around 440 cm
-1.
In the spectrum of the 1.59 eV excitation energy, we can observe an intense and sharp peak associated with the out-of-plane A
1’ mode and the broad 2LA band. Notice that the in-plane E
’ mode is absent in this spectrum. In the 1.68 eV excitation energy spectrum, the A
1’ mode is still the most intense, the 2LA band is observed and we can see the appearance of a broad band in the position of the E
’ mode. In the 1.92 eV excitation energy spectrum, we can observe that the broad 2LA band becomes the most intense, and the appearance of the E
’ mode is clearly observed. The broad band near the position of the E
’ mode is related to a double-resonance Raman process involving different valleys of MoS
[
40]. By increasing the excitation energy, we can observe in the 2.38 eV excitation energy spectrum that the E
’ mode becomes as intense as the A
1’ mode, and the 2LA band is absent in the spectrum. In the highest 2.73 eV excitation energy spectrum, we see that the E
’ mode becomes more intense than the A
1’ mode.
Figure 2(b) shows the result of a multiple excitation energy Raman map at the MoS
2 dome’s center obtained using different 23 laser lines with excitation energies ranging from 1.59 to 2.73 eV. The peak intensities in these figures were normalized by the quartz intensity, and are represented by the color bar. The energies of the incident photon are on the vertical scale, and the horizontal scale represents the Raman shifts. The blank gap in excitation energies between 1.69 and 1.88 eV in
Figure 2(b) is a region where we do not have available laser lines.
Figure 3 shows the Raman excitation profiles (REPs) of the A
’ , E
’ and 2LA bands, that is, the intensity of each feature as a function of the laser excitation energy. We will discuss separately the results of the first-order modes and the 2LA band. The REPs of the first-order A
1’ (open red squares) and E
’ (dark blue triangles) modes are shown in
Figure 3, and the lines represent the best fits of the experimental data considering the equation for Raman intensity as a function of the laser energy E
L for a first-order (one-phonon) process given by [
40]:
where the index
i denotes the electronic or excitonic transition, the index
k denotes the two phonons (A
1’ or E
’) and the three terms in the numerator represent the matrix elements of the electron-radiation (absorption of the incident photon), exciton-phonon, and electron-radiation (emission of the scattered photon) interactions, respectively. The two terms in the denominator give rise to the resonant enhancement of the Raman peaks when the incident or scattered photon energies match the exciton energy. The damping constant
is related to the finite lifetime of the exciton
i involved in the Raman process, E
is the energy of exciton
i and E
is the energy of phonon
k.
Figure 3(a) clearly shows three resonances in the REPs in the investigated range of energies. Notice that the out-of-plane A
’1 mode is enhanced by the three resonances, whereas the in-plane E
’ mode is weakly enhanced at lower energy resonance, but is more intense than the A
1’ mode at higher energies in
Figure 3(a).
In the fitting process, the same set of parameters of the excitons (energies and damping constants) was used to fit the REPs of the different phonon modes.
Table 1 shows the fitting parameters (exciton energies and damping constants) that provide the best fit of the experimental REP data and are represented by the solid curves in
Figure 3(a). The values for the first resonance (E
) are more accurate since we have more experimental data in this energy region. The lack of experimental points in the range 1.69-1.88 eV prevented us to measure with accuracy the second resonance (E
) energy and linewidth. The accuracy in the value of the third resonance energy (E
) is also poor since we do not have experimental points above 2.8 eV.
We will now discuss the resonance Raman behavior of the 2LA band, that comes from a second-order phonon process involving two phonons of the longitudinal acoustic (LA) branch near the Brillouin zone edge (K and M points)[
39].
Figure 3(b) shows the REP of the 2LA mode, where the blue circles represent the experimental data and the curves represent the best fit considering the expression for the intensity of a second-order Raman process as a function of the laser energy E
L given by following expression[
39]:
where the two middle terms in the numerator of Eq.
2 represent the exciton-phonon interactions involving two phonons with opposite momenta. The denominator shows three terms that gives the resonances of the incident photon and the scattered photon with one and by two phonons. The REP of the 2LA Raman band shown in
Figure 3(b) was fitted using the same energy values of the E
, E
and E
shown in
Table 1. We can observe that the 2LA band is significantly enhanced at lower energies (E
and E
) but only very weakly enhanced at the higher energy resonance (E
).
Figure 3(c) shows the photoluminescence spectra recorded at the dome center (upper spectrum), at the MoS
bulk substrate (middle spectrum), and the subtraction between the two first ones in the bottom, where we can only observe the contribution of the strained MoS
single layer.
4. Discussion
Let us now discuss the physical origin of the three resonances E, E and E observed in the resonace Raman results. For the first-order A1’ and E’, we have observed that the out-of-plane A1’ mode is strongly enhanced at E and E, whereas the in-plane E’ is only strongly enhanced at the higher resonance energy E. On the other hand, the second-order 2LA band is significantly enhanced at E and E and practically not observed at E.
The highest resonance energy E
can be attributed to exciton C of MoS
, whose energy in unstrained monolayer MoS
is around 2.9 eV as observed by resonance Raman spectroscopy [
41], and is far isolated from the other lower energy excitons in the optical spectrum. The observed value of the C-exciton energy at 2.72 eV allows us to conclude that is redshifted by about 0.18 eV in the strained domes with respect to the unstrained MoS
monolayer, similarly to the case of the energies of A and B excitons in the strained MoS
domes. Interestingly, the 2LA band is only very weakly enhanced at higher energies, revealing thus a weak interaction of the zone-edge LA phonons with the C exciton.
In order to assign the E
and E
resonances in the REP we acquired the photoluminescence (PL) spectrum on the same dome, whose REP curves are displayed in
Figure 3(a,b). The topmost trace in
Figure 3 (c) is the PL spectrum recorded with the laser spot centred on the top of the MoS
dome. The emission from the MoS
substrate (or bulk) adjacent to the dome is given by the middle trace in
Figure 3 (c). Finally, the bottommost trace in
Figure 3 (c) is obtained by subtracting the bulk contribution from the spectrum recorded on the dome. It must be emphasized that the spectrum intensities were not corrected by the spectral response of the spectrometer (detector and gratings), and the intensities at lower energies are underestimated. Nevertheless, in the PL spectrum two clear resonances are observed. The latter match very well the E
and E
REP resonances, which are indicated by the arrows in the lower part of
Figure 3 (c).
Previous PL studies on strained MoS
domes similar to those investigated here showed the presence of two main contributions to the emission spectrum acquired close to the top of the dome [
11] like in the present case. The lower and higher energy resonances in the PL spectrum were attributed to the indirect and direct exciton, respectively. As a matter of fact, for sufficiently high strain values (typically greater than 2%) the maximum of the valence band (VB) of the MoS
monolayer undergoes a transition from K to
(the same occurs in WS
and WSe
monolayers as reported in Ref. [
11]). At the same time, the minimum of the conduction band (CB) at K moves at lower energy while remaining the lowest state of the CB. Consequently, the exciton transition with lowest energy becomes indirect in character (
-K
), while the direct exciton transition (K
-K
) is at higher energy. It should be noted though that when the direct and indirect excitons are resonant they hybridise and their direct vs indirect character is smeared out. In particular, the indirect exciton may gain sufficiently oscillator strength and become bright despite its k-space indirect character. Recently, evidence of exciton hybridisation was also observed in the strain dependence of the exciton magnetic moment in WS
domes [
13]. By comparing the PL difference spectrum in
Figure 3 (c) with the results on MoS
domes reported in Ref. [
11], we may attribute the E
resonance in the REP to the indirect exciton and the E
resonance to the direct exciton. Indeed, in Ref. [
11], it was observed that the indirect exciton energy ranges between 1.65 eV and 1.62 eV depending on the position on the dome, in agreement with the E
resonance. The direct exciton was found to vary from 1.75 eV to 1.78 eV for increasing strain that suggests the E
resonance (E
=1.82 ± 0.05 eV) being associated to the direct (or A) exciton state. Nevertheless, due to the large energy uncertainty of E
and its low spectral weight, we cannot exclude a possible contribution from the B exciton.
Finally, the different relative weight of the indirect (corresponding to E
) and direct (corresponding to E
) excitons in the PL spectra and in the REP curves probably comes from the exciton-phonon (or electron-phonon) interaction. The PL process involves two optical transitions and the Raman process involves not only the two optical transitions, but also the exciton-phonon intereaction, whose matrix element is given by the middle term in the numerator of Eq.
1, which is specific for each exciton
i and phonon
k.