New first order Raman-active modes in few layered transition metal dichalcogenides.

Although the main Raman features of semiconducting transition metal dichalcogenides are well known for the monolayer and bulk, there are important differences exhibited by few layered systems which have not been fully addressed. WSe2 samples were synthesized and ab-initio calculations carried out. We calculated phonon dispersions and Raman-active modes in layered systems: WSe2, MoSe2, WS2 and MoS2 ranging from monolayers to five-layers and the bulk. First, we confirmed that as the number of layers increase, the E', E″ and E2g modes shift to lower frequencies, and the A'1 and A1g modes shift to higher frequencies. Second, new high frequency first order A'1 and A1g modes appear, explaining recently reported experimental data for WSe2, MoSe2 and MoS2. Third, splitting of modes around A'1 and A1g is found which explains those observed in MoSe2. Finally, exterior and interior layers possess different vibrational frequencies. Therefore, it is now possible to precisely identify few-layered STMD.

The field of 2-Dimensional semiconducting transition metal dichalcogenides (STMDs) is growing very fast mainly due to: a) their moderate band gap which can be used for fabricating electronic devices such as field effect transistors (FETs)1234567; b) the presence of a direct band gap in monolayered STMDs, ideal for optoelectronic devices891011, and c) their unique valley polarization properties121314. The direct band gap of monolayered SMTDs has been measured and it results in an intense photoluminescence (PL) signal911 which can be enhanced at the edges possibly due to edge passivation15. The non-destructive technique of Raman spectroscopy has been intensively used to characterize monolayered STMDs, and other than the first-order modes, double resonant and triple resonant processes have been recently identified1617. Recently, important efforts in understanding the Raman spectroscopy of few layered systems have been carried out in WSe21819. However, a careful, systematic and deep study of the different Raman modes in few layered STMDs (FLSTMD) is needed since previous works use either the monolayered or the bulk systems to explain Raman spectra of FLSTMDs. Unfortunately, this kind of analysis might lead to an incorrect assignment of the active Raman modes20 since FLSTMDs exhibit different symmetries depending on the number of layers. In this manuscript, using first principles density functional perturbation theory (DFPT), we shed light on this subject and find that the particular symmetry in FLSTMD (e.g. MoS2, WS2, WSe2, MoSe2) results in new Raman-active modes that are neither present in the monolayer nor in the bulk spectra. In particular, we predict the existence of a new out-of-plane Raman-active mode, which is expected to be intense in the spectra of FLSTMDs such as WSe2, MoSe2 and MoS2. This mode has already been experimentally reported212223, but neither fully understood nor studied in detail in different FLSTMDs. We present here a polarized Raman study of this new mode, on synthesized chemical vapor transport samples242526, which supports our prediction. In our calculations, we additionally found a splitting in all Raman-active modes which explains the recently observed double and triple modes in MoSe2 when L = 3, L = 4 and L = 5 for the A′1 and A1g22 (see table 1), in addition to other modes related to in-plane vibrations that differentiate between exterior (surface) and interior (bulk) layer vibrations.

Table 1

First order Raman-active modes in cm−1 of few layered WSe2, MoSe2, MoS2 and WS2 (from L = 1 to L = 5 layers) and the bulk, calculated using DFPT, and experimental results published in the literature (in yellow). The point symmetry of the unit cell is indicated in the first row by its Shöenflies symbol. An asterisk * indicates the more intense calculated signals, “ir” indicates that the mode is also infrared-active, “in” indicates that just the interior layers vibrate and “ex” represents that just the exterior layers vibrate. The lower frequency modes have not been included since some of them do not appear in the most common Raman configurations

Structure	L = 1(cm−1)D3h	L = 2 (cm−1)D3d	L = 3(cm−1)D3h	L = 4(cm−1)D3d	L = 5(cm−1)D3h	Bulk(cm−1)D6h
WSe2	175.75(E″)	176.20(Eg2)	175.01(E″2)	175.22(Eg2)	174.69(E″2)	175.01(E1g)
249.36(E′)*ir	176.26(Eg2)19	175.58(E′2)ir	175.52(Eg2)19	175.07(E′2)ir	177.5(E1g)46
250.23(A′1*)	248.52(Eg1)*	175.84(E″2)19	176.38(Eg2)	175.60(E″2)	247.78(E2g)*
249.5(E′)21	250.83(A1g)*	176.33(E″2)	247.66(Eg1)*In	176.01(E′2)ir	251.61(A1g)*
249.5(A′1)21	306.97(A1g2)*	247.83(E′)*ir,In	248.35(Eg1)*Ex	176.23(E″2)	248.0(E2g)22
	309.00(A1g2)22	248.32(E′)*ir,Ex	249.57(A1g)	247.55(E′1)*ir,In	250.8(A1g)22
310.00(A1g2)21	248.31(E″1)*Ex	251.30(A1g)*	247.65(E′1)ir,In	247.00(E2g)30
	249.16(A′1)	305.57(A1g2)	248.55(E′1)*ir,Ex	251.00(A1g)30
251.23(A′1)*	306.62(A1g2)*	247.59(E″1)In	250(E2g)46
306.57(A′12)*	309.00(A1g2) 22	248.54(E″1) Ex	253(A1g)46
309.00(A′12)22	310.00(A1g2)21	248.83(A′1)
310.00(A′12)21		250.15(A′1)
	251.21(A′1)*
305.29(A′12)
306.84(A′12)*
309.00(A′12)22
310.00(A′12)21
MoSe2	167.60(E″)	167.72(Eg2)	166.55(E″2)	166.96(Eg2)	166.37(E″2)	166.81(E1g)
240.27(A′1)*	284.03(Eg1)*	167.51(E′2)ir	167.94(Eg2)	166.80(E′2)ir	168.00(E1g)30
240.5(A′1)22	285.90(Eg1)22	168.00(E″2)	239.49(A1g)	167.35(E″2)	241.84(A1g)*
241.2(A′1)34	240.68(A1g)*	238.74(A′1)*	241.23(A1g)*	167.84(E′2)ir	242.5(A1g)22
285.32(E′)*ir	348.04(A1g2)*	241.14(A′1)*	282.64(Eg1)*In	168.10(E″2)	243.00(A1g)34
287.2(E′)22	353.00(A1g2)22	282.69(E′1)*ir,In	284.06(Eg1)*Ex	238.42(A′1)	283.11(E2g)*
287.3(E′)34		284.18(E′1)*ir,Ex	345.72(A1g2)	238.03(A′1)22	283.7(E2g)34
	284.17(E″1)Ex	347.92(A1g2)*	240.02(A′1)
348.07(A′12)*	353.00(A1g2)22	240.12(A′1)22
353.00(A′12)22		241.40(A′1)*
	242.04(A′1)22
282.63(E′1)*ir,In
282.81(E′1)ir,In
284.07(E′1)*ir,Ex
282.73(E″2)In
284.10(E″2)Ex
345.54(A′12)
347.90(A′12)*
353.00(A′12)22
WS2	298.54(E″)	299.08(Eg2)	297.22(E″2)	297.73(Eg2)	297.10(E″2)	299.20(E1g)
359.24(E′)ir	358.26(Eg1)*	298.22(E′2)ir	298.94(Eg2)	297.83(E′2)ir	357.21(E2g)*
356.00(E′)16	419.39(A1g)*	298.92(E″2)	357.13(Eg1)*In	298.55(E″2)	356.5(E2g)21
357.50(E′)21	437.61(A1g2)	357.17(E′1)*ir,In	358.06(Eg1)*Ex	299.36(E′2)ir	420.27(A1g)*
418.7(A′1)*		358.13(E′1)*ir,Ex	417.69(A1g)	299.47(E″2)	421.00(A1g)21
417.50(A′1)16	358.13(E″2)Ex	419.81(A1g)*	357.09(E′1)ir,In
418.0(A′1)21	416.87(A′1)	436.14(A1g2)	357.21(E′1)*ir,In
	419.63(A′1)*	437.28(A1g2)	358.64(E′1)*ir,Ex
437.09(A′12)	357.23(E″1)In
358.66(E″1)Ex
416.73(A′1)
418.49(A′1)
420.14(A′1)*
435.75(A′12)
437.57(A′12)
MoS2	284.00(E″)	283.78(Eg2)	281.91(E″2)	281.79(Eg2)	281.84(E″2)	282.05(E1g)
	383.61(E′)*ir	381.65(Eg1)*	283.24(E′2)ir	283.28(Eg2)	282.29(E′2)ir	286.00(E1g)47
	384.7(E′)34	402.86(A1g)*	283.83(E″2)	379.57(Eg1)*In	282.95(E″2)	287.00(E1g)48
	384.0(E′)29	463.84(A1g2)*	379.89(E′1)*ir,In	381.54(Eg1)*Ex	283.57(E′2)ir	380.35(E2g)*
	402.65(A′1)*		381.75(E′1)*ir,Ex	400.67(A1g)	284.01((E′2)ir	383.00(E2g)48
	406.1(A′1)34		381.74(E″1)Ex	402.80(A1g)*	380.14(E′1)*ir,In	383.50(E2g)22
	403.00(A′1)29		400.00(A′1)	461.12(A1g2)	380.50(E′1)ir,In	383.7(E2g)34
			403.05(A′1)*	463.39(A1g2)*	381.46(E′1)ir,Ex	405.02(A1g)*
			463.62(A′12)*		380.33(E″1)In	408.6(A1g)22
					381.47(E″1)*Ex	408.7(A1g)34
					399.73(A′1)
					401.43(A′1)
					403.44(A′1)*
					461.28(A′12)
					463.18(A′12)*
					463.00(A′12)23

Bulk trigonal prismatic STMDs belong to P63/mmc (D6h4), which is a nonsymmorphic space group, and thus, in order to analyze their vibrations in the unit cell, we need to consider that their unit cell exhibits D6h point group symmetry; this is not the case for the monolayered or FLSTMDs, which are considered quasi 2-D systems in the “Z” direction. The monolayer in STMDs is non-centrosymmetric (no inversion symmetry) and its vibrational modes belong to the D3h point group. Therefore, in agreement with the irreducible representations of this point group, one must use the symbols A′1 and E′, for the out of plane and in-plane vibrations, respectively. For understanding the case of FLSTMDs (L = 2,3,4,5…), their particular symmetries need to be analyzed since the symmetry changes when varying the number of layers from even to odd numbers. Systems with an odd number of layers (L = 1,3,5,…) exhibit D3h point-group symmetry, do not possess a center of inversion, and the out-of-plane Raman-active modes belong to the A′1 irreducible representation. For systems with an even number of layers (L = 2,4,…), the modes belong to the D3d point group which exhibits a center of inversion, and the out-of-plane Raman-active modes possess the A1g irreducible representation. In this context, it is important to emphasize that the absence of inversion symmetry plays a crucial role in nonlinear optical properties such as the second harmonic generation2728, and valley polarization121314.

Results

We synthesized crystals of the trigonal prismatic phase of WSe2 by chemical vapor transport techniques242526 (see methods section), and by adhesive tape exfoliation we were able to produce few layered WSe2 (See figure 1). Raman spectra were attained using four different laser lines (488 nm, 514.5 nm, 633 nm and 647 nm) to study their vibrational properties depending on the number of layers (See figure 2a–2d), and photoluminescence studies were carried out to confirm the presence of the direct optical band gap at around 750 nm (1.65 eV) for the monolayer22. It is noticed that as the number of layers increases, the PL intensity decreases, appearing two signals which shift to higher wavelengths (lower energies) agreeing with the work of Tonndorf et al22 (see figure 2e).

Figure 1

Optical and atomic force microscope images of exfoliated WSe2 crystals a) Monolayer; b) Bilayer; c)Tetra-layer; d) Penta-layer.

Figure 2

Raman spectra of few layered WSe2 with different excitation wavelengths; (a) 488 nm; (b) 514.5 nm; (c) 633 nm; (d) 647 nm; (e) Photoluminescence (PL) spectra obtained with the 488 nm.

We first considered the case of monolayered WSe2, in which the experimental Raman spectrum exhibits the presence of the perpendicular mode A′1 and the in-plane E′ almost degenerate at around 250 cm−1 192122; according to our calculations A′1 is at 250.23 cm−1 and the E′ at 249.36 cm−1 (see table 1). By adding layers, we observed that the out of plane modes A′1 (A1g) shift to higher frequencies and the E′ (Eg) displace to lower frequencies (see table 1), a behavior that has been reported experimentally by different authors in WSe2 and other STMDs151921222329. Experimentally, these modes in the bulk 3-D crystals are associated with the A1g (located at 251 cm−1) and the E2g (around 247 cm−1) irreducible representations of the D6h point group, exhibiting inversion symmetry (see table 1)2230.

From the Raman spectra we noticed that with the laser line of 514.5 nm a new peak located at 310.4 cm−1 appears for the bilayer, tri-layer, tetra-layer and penta-layer (see figure 2a–d). Other authors have identified this peak at 308 cm−1 21, 309 cm−1 22 and 310 cm−1 19. One report assigned this peak to a combination of shear modes and the E2g mode21, and another publication related this peak to the inactive B2g mode of bulk WSe222. Our calculations on FLSTMD predicted Raman-active modes around 306 cm−1 (See table 1) that are close in frequency to the infrared-active A2u (305 cm−1)31 and the Raman-inactive B2g mode of the bulk phase (301 cm−1)30. These new Raman-active modes involve the vibrations of the chalcogen and the transition metal atoms, belonging to the A′1 or A1g irreducible representations (for odd and even number of layers, respectively), and do not appear in the monolayer. We propose to label these new modes A′12 and A1g2 (See table 1 and figure 3). For the bi-layer (306.97 cm−1) and tri-layer (306.57 cm−1), only one Raman-active mode appears around this frequency, but for the tetra-layer (305.57 cm−1 and 306.62 cm−1) and penta-layer (305.29 cm−1 and 306.84 cm−1), two modes appear that are very close in frequency, corresponding to exterior and interior layer vibrations (See table 1 and figure 4).

Figure 3

Experimental Raman spectra of WSe2 taken with a 514 nm laser under conventional geometry (VV) and cross polarization (VH).

(a) Monolayer (L = 1). (b) Bilayer (L = 2). (c) Trilayer (L = 3). Insets show models of the out of plane vibrational modes present in each case. (d) Calculated DFPT phonon dispersion of the WSe2 monolayer.

Figure 4

Models showing the vibrations of A′1, E′ and E″ Raman-active modes for L = 5 WSe2.

A′1(1) at 248.83 cm−1; A′1(2) at 250.15 cm−1; A′1(3) at 251.21 cm−1; A′12(1) at 305.29 cm−1; A′12(2) at 306.84 cm−1;E′1(1) at 247.55 cm−1; E′1(2) at 247.65 cm−1; E′1(3) at 248.55 cm−1; E″1(1) at 247.59 cm−1 and E″1(2) at 248.54 cm−1.

In order to experimentally demonstrate that these first order high frequency modes belong to vibrations perpendicular to the lattice plane, parallel and cross polarized Raman spectra were acquired using a 514.5 nm laser, confirming that both the A1g and A′1 signals for systems with different number of layers vanish completely in the cross-polarized Raman configuration (VH), as expected by group theory (See figure 3). In figure 3 it is also shown a signal around 260 cm−1 for monolayer, bilayer and trilayer which under no strain could be regarded as a second order peak caused by a double resonance effect involving the LA(M) phonon located around 130 cm−1 (See figure 3d) which will produce a peak circa 260 cm−1; double resonance effects involving the LA(M) phonon have already been reported in WS215.

Interestingly, these new modes have also been reported experimentally in MoSe2 at ca. 353 cm−1 22, and appear in our calculations between 345.54 cm−1 and 348.07 cm−1, depending on the number of layers (See table 1). For few-layered MoS2 and WS2, the high frequency out of plane modes, according to DFPT, are around 463 cm−1 and 437 cm−1, respectively. Unfortunately, these modes have not been properly identified experimentally for several reasons. First, for MoS2 the mode at 464 cm−1 is close to the 2LA(M) mode located at 455 cm−1 and it is often confused with it, and sometimes assigned to the Raman inactive mode in the bulk crystal, the A2u mode, but no further discussion has been made in the literature2332. Second, for few layered WS2 the mode centered at 437 cm−1 exhibits a very small intensity, probably due to the small polarizability of this mode (the components of the Raman tensor are very small), thus making it difficult to identify this mode experimentally.

In addition, our calculations in FLSTMDs reveal Raman-active splitting modes, both out of plane and in-plane; this can be seen for instance from the phonon dispersions of WSe2 at the Γ point (see figure 1S in supplementary information). While mono-layered and bi-layered WSe2 exhibit only one perpendicular mode at 250.23 cm−1 and 250.83 cm−1, respectively, the tri-layer (249.17 cm−1 and 251.23 cm−1) and tetra-layer (249.58 cm−1 and 251.3 cm−1) possess two branches, and the penta-layer reveals three (248.83 cm−1, 250.15 cm−1 and 251.21 cm−1) (See table 1 and figure 4). Recently, these three out of plane modes have been reported experimentally in five layered MoSe2 at 238.03 cm−1, 240.12 cm−1 and 242.04 cm−1 22; in our calculations, these modes appear at 238.42 cm−1, 240.02 cm−1 and 241.4 cm−1 and are in very good agreement with the experimental data (See table 1). The main reason why we cannot experimentally resolve them in WSe2 is because the in-plane modes are very close in frequency to the perpendicular modes, thus making this task very difficult and dependent on the resolution of the instrument. In other FLSTMDs, such as WS2 and MoS2, the in-plane and out of plane signals are separated by several wave numbers, around 60 cm−1 for WS2 and 20 cm−1 for MoS2, thus making possible to identify the splitting when having a high resolution instrument working at low temperatures. For WS2, our calculations depict that there should be two modes for both the tri-layer and tetra-layer between 416.87 cm−1 and 419.81 cm−1 (see table 1), and for the penta-layer there should be three modes between 416.73 cm−1 and 420.14 cm−1. In the case of the MoS2 tri-layer and tetra-layer, the two splitting modes appear between 400 cm−1 and 403 cm−1, and for the penta-layer, three modes appear between 399.73 cm−1 and 403.44 cm−1 (See table 1). On the other hand, it might be interesting to try finding the splitting in other 2D materials.

Our calculations indicate that the splitting of A′1 (A1g), starting from L = 3, involves two main types of mode families: One, in which all the Se or S atoms compress or expand at the same time, thus shifting to higher frequencies as the number of layers grows (reaching the maximum frequency in bulk), and another, in which the chalcogen vibrations in each layer exhibit opposite directions in an alternative way, resulting in lower frequencies (See table 1 and figure 4).

It is also noteworthy that the new in-plane modes can differentiate between interior and exterior layers, and belong to the Eg irreducible representation for even numbers of layers, and to the E′ and E″ for odd numbers of layers. In WSe2 the in plane modes are around 176 cm−1 and 248 cm−1 agreeing with the calculations presented by Luo and co-workers for L = 2 to L = 419 (see table 1). In our Raman spectra, despite using different laser lines, it is difficult to observe a signal around 176 cm−1 due to the small intensity of these modes, however, around 248 cm−1 we observe the signal of these modes which is mixed with the out of plane modes (see figure 2). For the case of MoSe2 the in-plane modes are located circa 167 cm−1 and 283 cm−1, in WS2 at about 298 cm−1 and 357 cm−1, and in MoS2 are nearby 282 cm−1 and 380 cm−1 (to differentiate between these in plane modes a super index has been used in table 1). For example, in WSe2, two in-plane modes appear for the tri-layer (247.83 cm−1 and 248.33 cm−1) and tetra-layer (247.66 cm−1 and 248.35 cm−1), and three modes appeared for the penta-layer (247.55 cm−1, 247.65 cm−1 and 248.55 cm−1. See figure 4). Further, it is worth noticing that Raman-active E″ modes are also active for the tri-layer at 248.32 cm−1 and the penta-layer at 248.54 cm−1, and only the exterior layers are involved (See figure 4). In addition, the E′ in-plane mode has its maximum frequency value for the monolayer, and when more layers are added, a splitting takes place in which the interior layers (bulk-like) exhibit lower frequencies than the exterior ones (surface layers). Therefore, when reaching the bulk phase (no exterior layers), the E1g mode preserves the lower frequency values of the interior layers, thus explaining the experimentally observed shift to lower frequencies as the number of layers is increased (See figure 4). All the E′ modes in odd numbers of layers are also infrared-active (See table 1), while for even numbers of layers the E1g modes, as expected from the principle of mutual exclusion, are not infrared-active since they are symmmetric with respect to inversion. The calculated phonon dispersions of WSe2 for L = 1 to L = 5 and bulk are shown in figure S1 (supplementary information) where the new modes appear as new branches at the Γ point. These results can be extended to other FLSTMDs with trigonal prismatic layers (See table 1).

We also expect that the new modes mentioned above will appear in hetero-bilayered systems formed by stacking different layers of STMDs. For example, for a bilayered system made of one layer of WSe2 and another layer of WS2, our calculations demonstrate that, besides the A1 (C3v symmetry) vibrational modes at 251.7 cm−1 for WSe2 and at 415.4 cm−1 for WS2, two other perpendicular Raman-active modes appear at different frequencies, one at 311.1 cm−1 for WSe2 and another at 427.9 cm−1 for WS2. Note that these modes are close in frequency to those present in the pure bilayered systems. However, the shift here is mainly due to an induced strain causing the WS2 to expand and the WSe2 to contract to preserve an energy local minimum33. Since our calculations consider free standing layers of STMDs, the effect of the substrate has not been considered and could affect mainly the layer to which it has interaction depending on the particular substrate: Several experimental works, including ours, use SiO2/Si substrates162122232934, and the Raman modes exhibit little change from those calculated here and by other authors when no strain is considered181935, indicating a weak interaction with the substrate. In addition, the photoluminescence results (See figure 2) show a clear signal at 750 nm which corresponds to the direct optical band gap of monolayer WSe2, thus no evident strain from the substrate or from the sample itself is detected since strain would alter the optical band gap353637. Nevertheless, a careful study with different substrates is required to address the role of a particular substrate with FLSTMDs.

In summary, we have addressed Raman spectroscopic modes of FLSTMDs by measuring directly few layered exfoliated WSe2 samples, and confirming our observations by performing DFPT first principles calculations. We find that there are novel, out of plane and in plane, Raman-active modes that are not present in the monolayer or in the bulk, and are capable of differentiating between the vibrations of the outer and interior layers. Our results also explain recent observations reported for few layered WSe2, MoSe2 and MoS2. These new modes reported here will be important when characterizing few-layered systems as well as to identify monolayers. In addition, the new modes might play an important role in double resonance Raman processes in FLSTMDs.

Methods

In order to study and understand the origin of the main features observed in Raman spectra of FLSTMDs, ab-initio density functional theory (DFT) and density functional perturbation theory (DFPT) calculations were carried out in monolayers and few-layered systems of WSe2, MoSe2, MoS2 and WS2 (from L = 1 to L = 5, and also the bulk), as implemented in the plane wave code CASTEP38 in the Materials Studio package. Trigonal unit cells were considered under the local density approximation (LDA) using the Ceperly-Alder-Perdew and Zunger (CA-PZ) functional3940 with 9 × 9 × 3 Monkhorst-Pack K-points and a plane waves cut off of 500 eV with a norm-conserving pseudopotential. All the structures were relaxed, including the cells, until the forces became smaller than 0.01 eV/Å and with self- consistent energy tolerances less than 5 × 10−7 eV/atom. A vacuum of 16 Å between the layers was considered. The electronic structure results obtained are in agreement with those reported by other groups within the LDA formalism4142. In order to provide a reasonable description of the van der Waals interaction, we have considered a dispersion correction for both LDA (LDA-D)43 and for general gradient approximation (GGA-PW91)44 finding that LDA provides a better description for the interlayer distance than LDA-D and corrected GGA-PW91, as described in a previous publication33. The phonon dispersions and Raman scattering modes were calculated with the above parameters, but extending the plane wave cut off to 720 eV with DFPT using the linear response methodology for insulators45.

The few layered WSe2 samples were obtained by mechanical exfoliation of WSe2 crystals produced using the chemical vapor transport method with iodine as transport agent242526. First, WSe2 powder was synthesized by heating a mixture containing stoichiometric amounts of tungsten (Acros Organics 99.9%) and selenium (Acros Organics 99.5+%) at 1000°C for 3 days in an evacuated and sealed quartz ampoule (10 mm ID, 12 mm OD, 150 mm length). The mixture was slowly heated from room temperature to 1000°C for 12 hours, in order to avoid any explosion due to the strong exothermic reaction. The WSe2 crystals were grown from that synthesized powder using chemical vapor transport with iodine (Sigma-Aldrich, 99.8+%) as the transport gas at 1.5 mg/cm3. The experiments ran for 10 days in an evacuated and sealed quartz ampoule (10 mm ID, 12 mm OD, 100 mm length). The source and growth zones were kept at 950°C and 816°C, respectively. The resulted crystals were washed with hexane and dried in vacuo to remove any residual iodine. Both WSe2 powder and crystals were analyzed using XRD, which confirmed that both WSe2 powders and crystals were pure 2H phase (trigonal prismatic). The micro-Raman measurements were performed in a triple monochromator spectrometer (DILOR XY) equipped with an Ar/Kr laser and in a Renishaw inVia confocal microscope-based Raman spectrometer using the 488 nm, 514.5 nm and 633 nm and 647 nm excitation wavelengths. A back scattering geometry at room temperature was used and the laser power was kept below 0.6 mW to avoid damage on the samples. The laser spot size using the 100× objective was about 4 μm2, always smaller that the dimensions of the flakes. A λ/2 polarizer was used to achieve the VV and VH configurations. The accumulations times varies between 5 and 10 minutes depending on the sample and the excitation wavelength.