Light–matter interactions in two-dimensional layered WSe2 for gauging evolution of phonon dynamics

  1. 1,2 ORCID Logo ,
  2. 3 and
  3. 1,2,3
1Department of Electrical Engineering, University of North Texas, Denton, TX 76203, United States
2Department of Materials Science and Engineering; PACCAR Technology Institute; University of North Texas, Denton, TX 76203, United States
3Department of Electrical and Computer Engineering, University of Texas, El Paso, TX 79968, United States
  1. Corresponding author email
Associate Editor: P. Leiderer
Beilstein J. Nanotechnol. 2020, 11, 782–797.
Received 22 Oct 2019, Accepted 23 Apr 2020, Published 12 May 2020
Full Research Paper
cc by logo


Phonon dynamics is explored in mechanically exfoliated two-dimensional WSe2 using temperature-dependent and laser-power-dependent Raman and photoluminescence (PL) spectroscopy. From this analysis, phonon lifetime in the Raman active modes and phonon concentration, as correlated to the energy parameter E0, were calculated as a function of the laser power, P, and substrate temperature, T. For monolayer WSe2, from the power dependence it was determined that the phonon lifetime for the in-plane vibrational mode was twice that of the out-of-plane vibrational mode for P in the range from 0.308 mW up to 3.35 mW. On the other hand, the corresponding relationship for the temperature analysis showed that the phonon lifetime for the in-plane vibrational mode lies within 1.42× to 1.90× that of the out-of-plane vibrational mode over T = 79 K up to 523 K. To provide energy from external stimuli, as T and P were increased, peak broadening in the PL spectra of the A-exciton was observed. From this, a phonon concentration was tabulated using the Urbach formulism, which increased with increasing T and P; consequently, the phonon lifetime was found to decrease. Although phonon lifetime decreased with increasing temperature for all thicknesses, the decay rate in the phonon lifetime in the monolayer (1L) material was found to be 2× lower compared to the bulk. We invoke a harmonic oscillator model to explain the damping mechanism in WSe2. From this it was determined that the damping coefficient increases with the number of layers. The work reported here sheds fundamental insights into the evolution of phonon dynamics in WSe2 and should help pave the way for designing high-performance electronic, optoelectronic and thermoelectric devices in the future.


Since the discovery of graphene, atomically thin two-dimensional layered materials have drawn intense attention due to their unique physical properties [1,2]. Two-dimensional (2D) layered materials beyond graphene, such as transitional metal dichalcogenides (TMDCs) [3], black phosphorus (BP) [4], and other families of layered materials [5], can be mechanically exfoliated or, in some cases, grown from bottom-up processes akin to graphene. While graphene is comprised of a single element on the periodic table, i.e., carbon, TMDCs are binary compounds which makes their lattice dynamics more complex compared to multilayer (ML) graphene [6]. The symmetry, force constants, and frequency variation with geometrical confinement in some TMDCs has been studied recently [7]. Monolayer (1L) TMDCs consist of a plane of a transition metal, M, sandwiched by chalcogenides, X, on either side to yield the stoichiometry MX2 [8]. The interlayer bonding in most ML TMDCs is through the weak van der Waals interaction while the intra-layer bonding is via the strong covalent interaction. This makes them inherently flexible and good candidates for flexible electronics [9], optoelectronics [10], and other related applications [11,12]. Amongst the TMDCs, WSe2 offers unique advantages for device applications, which includes its high mobility of ≈500 cm2/V·s at room temperature, and a strong spin–orbit coupling [3,13,14]. Thus, it is not surprising that a rich variety of electronic and optoelectronic devices have already been demonstrated using 1L WSe2 which harnesses its exceptional properties [13,15,16].

It is well-understood that the underlying factors governing the optical, electronic and thermal properties of solid-state materials are strongly influenced by phonons and their spatio-temporal response toward external stimuli. Raman and photoluminescence (PL) spectroscopy has been a remarkable tool to gauge phonon dynamics for a broad range of materials in the past, including nanocarbons [17]. Phonon dynamics in 2D TMDCs, just as in other materials, includes discerning factors such as phonon lifetime τ and the change in phonon concentration as determined from the characteristic energy parameter E0, calculated from the slope of the low-energy edge of the excitonic mode of the PL spectra. The theory of spectral line shape in the Raman spectrum predicts a Lorentzian distribution of a collected signal in a dispersive medium, where the full-width-half-maximum (FWHM) scales as 1/τ, and not surprisingly, τ is influenced by damping mechanisms. The FWHM is expected to be infinitesimally small for activated phonons in a dissipationless medium, and the crystal elastic waves of the harmonic oscillator model for the allowable phonon modes would thus yield an exceptionally large τ. However, natural systems inherently exhibit damping, and thus the FWHM of the Raman peaks have a finite width, indicating the presence of decay channels that reduce τ. In general, the phonon linewidths contain contributions arising from several scattering mechanisms such as the electron–electron interaction, i.e., Coulombic scattering, or the electron–phonon interaction, i.e., scattering of electrons from defects. Chakrabarty et al. reported that the linewidth of the A1g peak in single-layer MoS2 that was subsequently used in transistors, broadened due to n-type doping where the phonon linewidth renormalized under the presence of an electric filed [18]. Similarly, the Raman linewidths in graphene are found to increase with defects resulting from electron–impurity and electron–phonon scattering [19]. Moreover, the Raman linewidth broadening is also attributed to the confinement of the optical phonons. Specifically, in the case of low-dimensionality nanocrystallites, the wave function of the optical phonons no longer remains a continuous plane wave and thus the localization of the wave function leads to a relaxation in the conservation of the wave vector selection rules. The phonons with a nonzero wave vector also take part in the Raman scattering process along with the phonons with zero wave vector, which results in broadening of the phonon linewidths. Another important parameter in this analysis is the actual position of the allowable Raman modes, which are typically the higher energy optical phonons, and how these modes interact with external stimuli. For example, external radiation could be in the form of heat or optical energy, which also directly influences properties such as the electronic and optoelectronic transport and the thermal conductivity of the material.

In this work, we have conducted an in-depth analysis of the phonon dynamics in WSe2, where τ and the change in phonon concentration deciphered from an energy parameter E0 were quantified as a function of external stimuli, specifically temperature, T, of the WSe2 and the Raman laser power, P, in 1L, ML and bulk WSe2 samples using Raman and PL spectroscopy. The temperature dependence of the Raman shifts in 2D TMDCs such as MoS2 [20-24], and WS2 [25,26] have been extensively studied over a wide temperature range from which properties such as thermal conductivity was deciphered [23,27]. On the contrary, the temperature-dependent Raman analysis of WSe2 is rather limited with only one prior report discussing the thermal coefficients of various Raman modes in WSe2 where a comparative analysis of temperature-dependent Raman modes in WSe2 and MoSe2 was presented [28]. Moreover, the study related to the phonon lifetime in WSe2 is also limited, among which one prior work reported τ over a range T = 4.4–300 K in naturally abundant and isotopically pure WSe2 grown by the chemical vapor deposition (CVD) method [29]. In one of our recent works, exciton dynamics and phonon lifetimes in CVD-grown and mechanically exfoliated 1L WSe2 nanosheets were analyzed and τ, obtained from temperature-dependent Raman measurements in the [Graphic 1] mode, was found to decrease at high temperature due to increased phonon-induced scattering events which eventually also reduces the exciton density in WSe2 [30]. In this present work, we studied the phonon lifetime, τ, and phonon concentration in 1L, ML and bulk WSe2 using Raman and PL spectroscopy respectively where both the temperature as well as the laser excitation power were varied to better understand the phonon dynamics in 2D WSe2 as they are vitally necessary to truly harness its intriguing properties for devices. For example, analysis of the phonon dynamics in 2D WSe2 will shed insights on the impact of self-heating effects in WSe2 to illustrate its utility in electronic, optoelectronic and thermoelectric device platforms in the future.

In this work, we demonstrate that exposure to heat on the WSe2 crystallites as generated via external stimuli such as T and P causes a red-shift in both the [Graphic 2] and A1g Raman-active modes. Similarly, the FWHM of the Raman peaks is examined which represents the anharmonic terms in the lattice potential energy, and a shift in the peak position and peak broadening effects are noted here as a function of increasing temperature, which is directly related to the phonon damping mechanisms. Finally, we propose a mechanical model to help explain the effect of damping or the rate of decrease of τ as a function of T for quantum-confined 1L, as well as ML and bulk WSe2.


WSe2 nanosheets, including all 1L, ML and bulk samples, were mechanically exfoliated on top of SiO2/Si substrates (SiO2 thickness = 270 nm) using the scotch tape method [1,30,31]. The mechanical exfoliation process used in this study is similar to the one reported in our earlier work [30]. Prior to the exfoliation procedure, the substrates were cleaned using acetone, isopropyl alcohol and methanol solution in an ultrasonicator, and rinsed with DI water. Nitrogen blow drying and substrate heating on a hotplate at T = 110 °C for 5 min allowed any residual moisture from the substrate surface to be removed. The micro-Raman and PL measurements were conducted using a LabRAM HR Evolution NIR (HORIBA Scientific) device equipped with a 532 nm laser for excitation. The micro-Raman and PL spectra were collected by the CCD detector with a diffraction grating of 1800 gr/mm. The spectral resolution for Raman and PL measurements was 0.09 cm−1. The Raman band for Si at 520.7 cm−1 was used as a reference to calibrate the spectrometer. The laser spot size and spatial resolution were calculated to be ≈2.6 µm and 1.3 µm respectively and the supporting calculations are outlined in Section 1 of the Supporting Information File 1. The WSe2 nanosheets were viewed, and the data were collected using a 10× objective (NA = 0.25) of the Raman microscope with working distance 10.6 mm. Due to the instrument limitation, the 10× objective lens could only be used for the temperature-dependent measurement and hence the 10× objective lens was used to collect the data throughout this study. The temperature-dependent Raman and PL measurements were conducted using liquid nitrogen (LN2) in a THMS600 Linkam temperature cell at ambient pressure. During the temperature-dependent measurements, the sample was first cooled gradually using LN2 from room temperature, i.e., T = 298 K to T = 79 K after which it was warmed again to T = 298 K. Finally heating from room temperature to T = 523 K was conducted in a ceramic crucible inside the Linkam cell where a temperature controller enabled the stage to heat the sample.

Results and Discussion

Phonon modes in 2D WSe2

Mechanically exfoliated high-quality WSe2 nanomembranes were observed under an optical microscope through optical contrast differences on the SiO2/Si substrates. The optical image of monolayer (1L), ML, and bulk WSe2 is shown in Figure 1a. Monolayer TMDCs such as WSe2 have a ground state structure with D3h symmetry and this phase is called 1H. In bulk structures, the stacking of individual layers results in an alternating rotated sequence which is called the hexagonal symmetric 2H phase that belongs to the inversion-symmetric D6h point group. Additionally, WSe2 has a lattice constant a = 3.28 Å [32] (Figure 1b), and a van der Waals gap g = 3.36 Å [33] (Figure 1c, left). The crystal structure of WSe2 (0001) consists of a repetition of Se–W–Se trilayers, as depicted in Figure 1b.


Figure 1: (a) Optical image of monolayer (1L), multilayer (ML), and bulk WSe2 nanomembranes. Raman and PL measurements were done on three different sections in the regions marked a, b and c. (b) Schematic representation (top and side views) of the layered structure of (0001) WSe2. Schematic representation of frequency evolutions of the (c) [Graphic 3] and (d) A1g modes for 1L, 2L and 4L of WSe2. The layered structure of WSe2 on top of a representative substrate is shown in (c)-left, and the application of external stimuli such as laser power P and substrate temperature T. The interlayer gap g is shown as ≈3.36 Å. (e) The variation of the Raman spectra for the [Graphic 4] and A1g modes for mechanically exfoliated WSe2 for 1L, ML, and bulk. The [Graphic 5] mode exhibits a red-shift Δω ≈ 0.47 cm−1 while the A1g mode reveals a blue-shift Δω ≈ 1.4 cm−1 with increasing layer number. The corresponding error bars calculated from the standard deviation (1σ) of the peak positions in [Graphic 6] and A1g are shown in Supporting Information File 1, Figure S1a. (f) The PL spectra for 1L, ML, and bulk WSe2 nanomembranes, where the excitonic A-peak represents direct-gap optical transitions, while the I-peak is characteristic of indirect band gap optical transitions. The corresponding error bars for the PL A- and I-peak shifts are shown in Figure S1b.

The irreducible representation of the phonon modes at the center of the Brillouin zone in WSe2 can be described by the following in Equation 1 [7],


among which A1g, E1g, [Graphic 7] [Graphic 8] are the Raman active optical modes, [Graphic 9] [Graphic 10] are the acoustic phonon modes, [Graphic 11] [Graphic 12] are the infrared active modes, and E2u, B1u, [Graphic 13] [Graphic 14] are the four optically inactive modes; however the [Graphic 15] mode is active for bilayer and thicker WSe2. In this study we focus on the Raman active modes that are dominant and show a thickness dependence for the [Graphic 16] and A1g modes.

The respective atomic vibrations related to these two Raman-active modes for 1L and ML WSe2 are represented schematically in Figure 1c (right) for the in-plane [Graphic 17] mode and (d) the out-of-plane A1g mode, respectively. Additionally, in Figure 1c (left) the layer-wise arrangement of WSe2 is shown schematically where the two external stimuli that were used in our experiments, i.e., P and T are also illustrated.

Layer-dependent Raman and PL analysis in WSe2

As stated previously, the WSe2 nanomembranes were mechanically exfoliated from which 1L, ML and bulk regions were identified as shown by the representative case in Figure 1a. Raman and PL spectroscopy was then performed using a Raman laser excitation wavelength of 532 nm on three spatially uniform sections designated 1L, ML, and bulk WSe2 with regions a, b and c on the same section marked, as indicated in Figure 1a. Data obtained from repeated measurements in regions a, b, c, served as a means to gauge the standard deviation and the repeatability of the measurements. The error bars thus are calculated from the standard deviation (1σ) of the Raman and PL measurements done on regions a, b and c on the sample for 1L, ML and bulk WSe2 which is discussed in detail in Section 2 of Supporting Information File 1. For a fixed laser power (P = 3.35 mW) and room temperature (T = 298 K), the two Raman-active [Graphic 18] and A1g modes were found at ≈247 cm−1 and ≈256 cm−1, (Figure 1e) for 1L WSe2 which is consistent with the previous reports [34-36]. As the thickness of the WSe2 increased, the [Graphic 19] and A1g modes shifted in opposite frequency directions; the [Graphic 20] mode experienced a red-shift while the A1g mode blue-shifted as thickness increased (Figure 1e). The blue-shift in the A1g mode as thickness increased is explained on the basis of the interlayer interaction of the Se atoms in the neighboring planes. From the nearest-neighbor interaction model it is expected, to first-order, that as film thickness increases, a greater restoring force from overlying layers will be present, where the equilibrium average lattice vibrational amplitude is reduced; hence the frequencies of the modes here in the out-of-plane direction will blue-shift as thickness increases [37].

The red-shift of the [Graphic 21] mode as thickness increases is attributed to dielectric screening effects of the long-range Coulomb interaction, where the effective charges resulting from the relative displacement between the W and Se atoms is reduced; this causes the coulombic force to decrease, and hence the energy also decreases [38]. We note from Figure 1e that the blue-shift Δω(A1g) ≈ 1.4 cm−1 is more than three times larger compared to the red-shift [Graphic 22] ≈ 0.47 cm−1 when the sample thickness increased from 1L to bulk. Another important observation from Figure 1e is the increase in intensity of the A1g peak as the number of layers increases, which is likely due to positive reinforcement arising from interference effects for the out-of-plane modes [39]. Prior reports have shown that the difference in peak position, as well as the relative intensities of the [Graphic 23] and A1g peaks, can serve as a fingerprint for identifying 1L membranes and appears to be a viable platform to decipher material thickness [40,41]. It was found from the data in Figure 1e that the Δω was 7.5 cm−1, 8.4 cm−1 and 9.1 cm−1 for 1L, ML and bulk WSe2, respectively. Thus, the Raman peaks exhibit a larger separation as thickness increases.

Monolayer WSe2 undergoes a transition from direct to indirect bandgap (Eg) as the number of layers increases [42,43]. We corroborate this here as well with our PL analysis of WSe2 for 1L, ML and bulk samples, where Figure 1f depicts a single excitonic A-peak at 1.61 eV for 1L which is consistent with the prior reports [35,44]. As the number of layers increases, an indirect peak I emerges, which is clearly seen from the PL spectra of the ML and bulk samples. In contrast, the I peak is absent in the 1L PL spectra, demonstrating that Eg undergoes an indirect to direct evolution for the ML to 1L case, respectively. It should also be noted that the I peak undergoes a red-shift from the ML to the bulk, and consequently, Eg was found to be 1.42 and 1.37 eV, respectively, for the two cases.

Power-dependent Raman and PL analysis in WSe2

Besides the above analysis of layer number on the [Graphic 24] and A1g modes, phonon dynamics are also strongly influenced by external stimuli. In this work, we explore two such external stimuli, specifically the Raman laser power P and substrate temperature T and comment on the light–matter interactions that evolve here in the context of phonon dynamics.

Power-dependent Raman shifts in WSe2

Starting with the discussion on the laser power, the dependence of the [Graphic 25] and A1g mode in WSe2 was analyzed as a function of P. Najmaei et al. [45] have analyzed the laser-induced thermal effects in 1L MoS2 and with increasing thickness approaching the bulk by means of Raman spectroscopy. Here, for the first time, we report on the variation of the [Graphic 26] and A1g modes for 1L, ML and bulk WSe2 crystallites as a function of P at T = 298 K. In order to minimize the laser-induced effects on the crystallinity and structural changes in the WSe2, a laser P ≤ 3.35 mW was used. It was also found that the minimum P required to get a good signal-to-noise ratio for the Raman and PL data was when P ≥ 0.31 mW. Therefore, the only four powers that could be used in this study were 0.31 mW, 1 mW, 1.76 mW and 3.35 mW, i.e., 0.31 mW ≤ P ≤ 3.35 mW. This is further explained in detail in Section 3 of Supporting Information File 1. It is also important to note that the P was tuned at 0.308 mW first and then it increased gradually to 3.35 mW while the acquisition time was fixed to as low as 0.1 s to further minimize any local heating that could be caused by the laser-induced effects. We termed the change in frequency of the peaks ([Graphic 27] and A1g) with P as Δωp. For the 1L case, the red-shift, i.e., Δωp = ωp (P = 3.35 mW) − ωp (P = 0.31 mW) for both modes was observed, where for the [Graphic 28], Δωp ≈ 0.27 cm−1, and for the A1g, Δωp ≈ 0.52 cm−1 as P increased (Figure 2a) for T = 298 K. Figure 2b clearly depicts the peak shift Δωp for the [Graphic 29] mode as P increased from 0.31 mW to 3.35 mW, while the inset shows the corresponding data for the A1g mode. The Δωp is characterized by Equation 2 below [27],


where χP is denoted as the power-dependent Raman coefficient and is obtained from the slope of the first-order fit to the Δωp versus P plot (Figure 2b). The fitted coefficients χP for the [Graphic 30] and A1g modes were computed to be −0.02173 cm−1/mW and −0.02558 cm−1/mW, respectively. The variation of χP for 1L, ML and bulk WSe2 is plotted in Figure 2c, which clearly shows χP for the [Graphic 31] mode increases while it decreases for the A1g mode.


Figure 2: Raman and PL analysis as a function of laser power P at T = 298 K. (a) Raman spectroscopy of 1L WSe2 obtained for P in the range of 0.31 mW to 3.35 mW. (b) The Δωp of the [Graphic 32] and A1g (inset) Raman active modes for 1L WSe2 with increasing P. The corresponding error bars calculated from the standard deviation (1σ) of the Raman shifts in [Graphic 33] and A1g are shown in Supporting Information File 1, Figure S2a. The results were fit linearly from which the slope χP was computed to be ≈0.02173 cm−1/mW and 0.02558 cm−1/mW for the [Graphic 34] and A1g modes, respectively. (c) The variation χP for the [Graphic 35] and A1g modes for 1L, ML and bulk WSe2. (d) Phonon lifetime τ for the [Graphic 36] and A1g modes as a function of P where [Graphic 37] ≈ 2× τ(A1g) for all P. The corresponding error bars calculated for the FWHMs of [Graphic 38] and A1g modes are shown in Figure S2b. (e) The normalized PL spectra of the excitonic A-peak as a function of P. (f) The excitonic A-peak position variation as a function of P (left-axis) and E0 (right-axis) representing the phonon concentration that is computed from the low-energy edge of the PL spectra in (e) using Urbach formulism [49] for 1L WSe2. The value of ηEP, calculated from Equation 6, was found to be ≈1.6 for 1L WSe2. The corresponding error bars calculated for the PL A-peak shift 1/E0 are shown in Figure S2c.

Phonon lifetime from power-dependent Raman analysis in WSe2

As mentioned earlier, phonon dynamics in 2D TMDCs also involves the determination of the phonon lifetime τ. The phonon lifetime τ is determined from the FWHM, i.e., phonon linewidth (Γ) of the [Graphic 39] and A1g peaks in the Raman spectra. The measured Raman linewidths are a convolution of effects of both the actual Lorentzian vibrational distribution of the phonons and the instrument-induced broadening which is typically assumed to have a Gaussian response and is provided by the Voigt profile [46]. The actual Γ was therefore determined from the Voigt profile fitted to the experimental data, illustrated in Figure S4 for 1L WSe2 at T = 298 K in Section 4 of Supporting Information File 1.

The instrumental broadening Δυ, calculated to be 0.613 cm−1, is lower compared to the measured phonon linewidths Γ in the Raman active [Graphic 40] and A1g modes from our analysis and hence, the influence of instrumental broadening on the calculation of phonon lifetimes from the phonon linewidth of the Raman active modes can be considered as negligible. The approach used to calculate the instrumental broadening is discussed in detail in Section 4 of Supporting Information File 1. As P increased, Γ was seen to increase based on the broadening of the Raman peaks for dispersive or dissipative media. As P increased, the incoming optical radiation thermalizes the lattice vibrations in both modes. Given the asymmetric nature of the lattice potential energy plot, a larger mean spacing between atoms arises at the higher effective temperature which increases the amplitude of the phonon vibrations as well as their potential for scattering events [20,47].

By using the energy uncertainty relationship with the phonon linewidth [46] we calculated τ using,


where τ is in picoseconds, and ℏ is the modified Planck’s constant (5.3 ps·cm−1). The τ for the [Graphic 41] and A1g modes is depicted in Figure 2d, where τ decreased as P increased due to more damping in the phonon vibrations caused by the increased likelihood of scattering events. We define a parameter, the power-dependent phonon lifetime ratio ητP as given by Equation 4,


For the 1L case ητP[Graphic 42] and ητP(A1g) was determined to be 0.954 and 0.945, respectively. Another artifact from this data in Figure 2d is the important observation that [Graphic 43] ≈ 2 × τ(A1g) for all of the values of P measured. This appears to corroborate that the coupling in the A1g mode between layers is rather weak given the larger separation; however, for the in-plane vibrations the phonon modes are well coupled and hence the phonon lifetime in the in-plane direction is larger when compared to the out-of-plane direction.

Phonon concentration from power-dependent PL analysis in WSe2

The other aspect of phonon dynamics is the phonon concentration and the change in phonon concentration in WSe2 due to change in excitation energy, which has been analyzed in this study. The PL measurement in 1L WSe2 as P increased at T = 298 K (Figure 2e) depicts the variation of the excitonic A-peak, representing direct-gap transitions, where Eg (1L) = 1.61 eV. The shift of the A-peak towards lower energies as P increased is further illustrated in Figure 2f (left axis), while a broadening of the lower energy side of the PL peak in Figure 2e was also observed as P increased. Urbach [48] discussed this dependence of spectrum in the context of various types of excitation sources used which formed the basis for the so-called Urbach formulism. Ramos and Luzzi [49] used this Urbach formulism to explain the behavior displayed in the radiation emission band of semiconductors at high-excitation levels, and the slope of the low-energy edge of the spectrum was characterized with an empirical parameter E0 on a semi-logarithmic plot. This parameter E0 has a dimension of energy and varies with excitation energy [49]. Ko et al. [50] also showed that the value of E0 increases with P and T in 1L MoS2 as a result of excess phonons generated when the system is perturbed beyond equilibrium. Similarly, in this study, from this low-energy edge due to excitons [51] and using Urbach formulism, the slope of the low-energy edge of the excitonic A-peak (Figure 2e) is calculated by the following Equation 5 [49],


where I(ω) is the intensity, and E0 is in units of meV and is the so-called characteristic energy parameter [50], which is correlated to the phonon concentration. The variation of E0 (right-axis) as a function of P is shown in Figure 2f, which indicates that the phonon concentration increases with increasing P. This is likely due to the excess phonons generated above equilibrium due to the higher effective temperature as P increases [50]. Once again, we define a parameter, the power-dependent energy parameter ratio ηEP, given by Equation 6,


The value of ηEP was calculated to be 1.6 which implies that E0 at P = 3.35 mW is 1.6 × E0 at P = 0.31 mW for 1L WSe2. Again, this increment in E0 is understandable due to the increase in phonon concentration as a result of increasing laser excitation energy.

Temperature-dependent Raman and PL analysis in WSe2

The temperature-dependent Raman analysis was carried out using a temperature cell where LN2 was used for cooling down to ≈79 K and a heater within the temperature cell allowed the temperature to reach 523 K. During the cooling process, it was found that the signal to noise ratio was very poor when both the T and P became very low, particularly in case of PL measurements. The best signal to noise ratio in both the Raman and PL measurements was found for P = 3.35 mW even when T was as low as ≈79 K and therefore the power was fixed at 3.35 mW throughout the temperature-dependent measurements in this study. It is also important to mention here that during the heating process, the intensity of the Raman signal started decreasing and above T > 523 K, no signal was detected, which may be attributed to material degradation arising in WSe2 through the possibility of an increased point defect density at higher temperatures, and/or through oxidative effects. From our preliminary checks on the samples, the Raman spectrum does not recover once the temperature is ramped down to 298 K. This indicates irreversible changes that may have transpired in the WSe2 as a result of the thermal cycle. In contrast, in the cooling down cycle from 298 to 79 K, and the subsequent warming back up cycle to 298 K, the Raman data was analyzed for any hysteresis effects. This data is shown in Section 5 of the Supporting Information File 1, Figure S5a and S5b for the [Graphic 44] and A1g modes in 1L WSe2, respectively in Section 5 of Supporting Information File 1. The largest difference in the Raman shift arising from hysteresis was found to be ≈0.21 cm−1 (T = 173 K) and 0.15 cm−1 (T = 248 K) for [Graphic 45] and A1g modes, respectively, suggesting that hysteresis is negligible for measurements at 298 K and below.

Temperature-dependent Raman shifts in WSe2

Figure 3a shows the temperature-dependent Raman spectra for 1L WSe2 for the [Graphic 46] and A1g modes. The observed changes in phonon frequencies with T in 1L WSe2 is attributed to the asymmetry in the interatomic potential versus displacement function [47], which leads to a larger average equilibrium lattice spacing. As the lattice expands or contracts with T, the equilibrium position of atoms and consequently the interatomic equilibrium separation changes, which induces a shift in the phonon energies. This behavior of the Raman peak frequencies with T is seen in many materials as T changes, from which a thermal expansion coefficient can be deduced [52,53]. Similar to the P dependence, here we termed the change in frequency of the peaks with T as ΔωT. As seen from Figure 3b, the red-shift of ΔωT[Graphic 47] was ≈6 cm−1 as T increased from 79 to 523 K, while its intensity remains constant. Similarly, the red-shift of ΔωT(A1g) was ≈4 cm−1; however the intensity of the A1g mode increased significantly with increasing T. This is likely due to the out-of-plane vibrations of the Se atoms that may be less constrained by the substrate for the A1g mode. In contrast, the in-plane vibrations of the [Graphic 48] mode are likely to be more constrained by the substrate at the Se-substrate interface even though temperature increased [54], and so the intensity of the [Graphic 49] mode did not increase significantly. In one of our prior reports [30], similar trends were observed in the intensities of [Graphic 50] and A1g mode in CVD-grown 1L WSe2. However, the ΔωT in [Graphic 51] and A1g modes were found to be less (ΔωT[Graphic 52] ≈ 5.10 cm−1, ΔωT(A1g) ≈ 3.54 cm−1) in the case of CVD grown 1L WSe2 [30], which could be attributed to the different film preparation methods used in both the studies [24].


Figure 3: Raman and PL analysis as a function of temperature T for P = 3.35 mW. (a) Raman spectroscopy of 1L WSe2 obtained for temperature in the range from 79 to 523 K. The corresponding error bars calculated from the standard deviation (1σ) of the Raman shifts and FWHM in [Graphic 53] and A1g modes are shown in Supporting Information File 1, Figure S2d and Figure S2e, respectively. (b) The net change in ΔωT and intensity of the [Graphic 54] and A1g modes from T = 79 K to T = 523 K. The net red-shift ΔωT was found to be ≈6 cm−1 and 4 cm−1 for the [Graphic 55] and A1g modes, respectively. (c) The frequency shift, ΔωT of the [Graphic 56] and A1g Raman active modes for 1L WSe2 with increasing T. The results were fit linearly using Equation 7 from which the slope χT was computed to be about −0.0145 cm−1/K and −0.0168 cm−1/K for the [Graphic 57] and A1g modes, respectively. The nonlinear perturbation of Raman shift for the [Graphic 58] mode mainly originates from the anharmonic effect of three-phonon processes, whereas for the A1g mode it is determined mainly by the three-phonon process in the low-T regime and in the high-T regime it is dominated by the thermal expansion effect. (d) The variation of χT for the [Graphic 59] and A1g modes with different layers of WSe2. (e) The normalized PL spectra of the excitonic A-peak as a function of T in 1L WSe2. (f) The excitonic A-peak position was used to calculate E0 and represents the phonon concentration that is computed from the low-energy edge of the PL spectra in (e) using Urbach’s formulism [49] for 1L WSe2. The phonon concentration E0 is shown for 1L, ML and bulk WSe2. The values of ηET, calculated using Equation 9, were found to be 8.31, 8.95 and 13.45 for 1L, ML and bulk WSe2, respectively. The corresponding error bars calculated for the PL A-peak shift 1/E0 are shown in Figure S2f.

The Se–substrate interface constraint may be assumed to be negligible in bulk WSe2. However, we noticed that the intensity of [Graphic 60] mode in bulk WSe2 (Supporting Information File 1, Figure S6) also remained less sensitive with temperature. This is likely due to the fact that the Se–substrate interface may not be the only constraining factor and there may be some other underlying mechanisms operative here, such as the the thermal properties of the substrate itself, or the sample preparation approach used, which could influence the vibrational properties and consquently the intensity of the in-plane [Graphic 61] mode in WSe2. Su et al. reported on the influence of the substrate, film morphology, and film preparation method on the Raman active modes in 1L, 2L and bulk MoS2 and showed that the [Graphic 62] mode in MoS2 is modified more strongly if an epitaxial sapphire substrate is used compared to SiO2/Si substrates [24]. As the SiO2/Si susbstrate and the preparation method (mechanical exfoliation) remained the same for 1L and bulk WSe2 used in our study, this may be the reason why the intensity of the [Graphic 63] Raman mode remained unchanged even in bulk WSe2. Moreover, in this temperature-dependent study, thickness was confirmed on representative 1L and bulk WSe2, as shown by the data in Figure S7 of Supporting Information File 1. The change in peak position of the [Graphic 64] and A1g modes with increasing T is plotted in Figure 3c. The results were fit to the Grüneisen model which was originally formulated for diamond [55], and is also applicable here, as given by Equation 7 below,


where ω0 is the frequency of vibration of the [Graphic 65] and A1g modes at T = 0 K and χT is the first-order Raman temperature coefficient calculated for both the [Graphic 66] and A1g modes (just as in the computation of χP, the slope of the linear fit represents the value of χT). The values for χP and χT for ML and bulk samples of WSe2 are listed in Table 1. The χT for the [Graphic 67] and A1g modes was computed to be −0.0145 cm−1/K and −0.0168 cm−1/K, respectively, which appear to be comparable to those reported for WSe2 [28] and other exfoliated TMDCs, such as MoS2 [22,23,56], MoSe2 [28], and WS2 [25]. The thickness dependence of χT for both the [Graphic 68] and A1g modes is plotted in Figure 3d, and it was found that χT for the A1g mode decreased, while on the other hand, it increased for the [Graphic 69] mode as the layer number increased. This dependence was similar to the variation noted for χP with layer thickness, as discussed earlier.

Table 1: Magnitudes of T- and P-dependent coefficients for the [Graphic 70] and A1g Raman-active modes and comparison to prior results reported for other TMDCs. Here χT[Graphic 71] represents the first-order temperature-coefficient for [Graphic 72]; χT(A1g) represents the first-order temperature-coefficient for A1g; χP[Graphic 73] represents the first-order power-coefficient for [Graphic 74]; and χP(A1g) represents the first-order power-coefficient for A1g. The χT and χP increased for the [Graphic 75] mode as the layer number increased, while a decrease for the A1g mode was noted.

Material T| [Graphic 76] (cm−1/K) T| (A1g) (cm−1/K) P| [Graphic 77](cm−1/mW) P| (A1g) (cm−1/mW) Ref.
WSe2 (1L) 0.0145 0.0168 0.02173 0.02558 this work
WSe2 (ML) 0.02015 0.01345 0.0247 0.02271 this work
WSe2 (bulk) 0.0214 0.0122 0.02671 0.02032 this work
WSe2 (1L) 0.0048 0.0032 [26]
MoSe2 (1L) 0.0054 [28]
WS2 (1L) 0.006 0.006 [25]
MoS2 (1L) 0.0124 0.0143 [22]
MoS2 (ML) 0.0132 0.0123 [23]
MoS2 (bulk) 0.0147 0.0123 [56]

Incidentally, the temperature-dependent Raman active [Graphic 78] and A1g modes actually exhibited some nonlinearity in certain regions that are described as nonlinear perturbations in Figure 3c. It was observed that the nonlinear perturbations of the Raman shift in the A1g mode were present over the entire temperature regime, whereas for the [Graphic 79] mode, a distinctively linear region is presented in the high-T regime above ≈350 K. To analyze the physical origin of these nonlinear perturbation, a physical model [26] is invoked as expressed by Equation 8,


where ΔωE is the Raman shift change induced by lattice thermal expansion, which leads to a red-shift as discussed earlier. Additionally, the ΔωA is the Raman shift attributed to anharmonic effects of the three- and four-phonon process. A light scattering process can be assumed to be comprised of: the absorption of a photon, the emission of a photon, and the creation of an optical phonon which then decays anharmonically into two phonons, three phonons, etc. The production of two and three phonons is known as a three-phonon and four-phonon process, respectively. The nonlinear perturbations in the [Graphic 80] mode mainly originate from the anharmonic effect of three-phonon processes (Figure 3c), whereas for the A1g mode, there is a contribution from the thermal expansion effect as well as in the high-T regime given the out-of-plane A1g vibrations that are less constrained by the substrate. However, the three-phonon anharmonic effect still presides in the low-T range for the A1g mode [26]. As a result of this the nonlinear perturbations were found to be persistent in the A1g mode even at high temperature. Additionally, the laser absorbance efficiency in WSe2 will change with temperature, which could also affect the phonon shifts in the [Graphic 81] and A1g modes, as found in Figure 3c.

Phonon concentration from temperature-dependent PL analysis in WSe2

Photoluminescence spectroscopy was carried out similarly as the temperature was varied at fixed P = 3.35 mW for 1L, ML and bulk WSe2 samples. The effect of temperature on the broadening and red-shift of the direct-band excitonic A-peak in 1L WSe2 is shown in Figure 3e. At high temperature, the phonon concentration increases, and the phonon scattering results in the broadening of exciton radiation [56]. This increased phonon scattering also causes the phonon lifetime to decay at high temperature, which is discussed in the next section. Additionally, the intensity of the PL A-peak also decreases at high temperature due to the temperature-induced enhancement of the exciton–phonon coupling and nonradiative recombination [57]. The change in phonon concentration was measured from E0 according to Equation 5, and the variation of E0 with T for 1L, ML and bulk samples is shown in Figure 3f. The energy parameter, E0, increases exponentially with T and its value is much higher for ML and bulk samples when compared with 1L samples across the entire temperature range shown in Figure 3f. Similar to ηEP, we define a parameter, the temperature-dependent energy parameter ratio ηEP, given by Equation 9,


The values of ηEP (1L), ηET (ML), ηET (bulk) were found to be 8.31, 8.95 and 13.45, respectively, which confirms the fact that the phonon concentration is higher in the bulk compared to the 1L case as T increased from 79 to 523 K. It is not surprising that in increasing temperature, more phonon modes are excited, which in turn increases the phonon concentration and hence the value of E0. The value of ηET (1L) > 5× ηEP (1L), which is possibly due to the more direct coupling of energy into the system by thermal means, in contrast to the local variation caused by an increase in the laser power where the energy coupling and the optothermal transduction may not be as effective. The following fit in Equation 10 was applied to the E0(T) data as a function of T,


where A1, B1 are constants having dimension of energy, and R0 is the exponential increment rate of E0. From the fitting, the value of R0 (1L) for WSe2 was found to be ≈0.0017/K, which is almost 3× < R0 (ML) and R0 (bulk), as shown by the comparative values in Table 2. This indicates the rate of increment in phonon concentration is much smaller in 1L WSe2 when compared with ML and bulk samples as temperature increases. This is again consistent with intuitive considerations where a larger number of phonon modes are available and excited in the 3D bulk state since phonon density increases as T3, while in a quantum confined system such as the 1L case, the phonon density is not as strong a function of T. This also has a direct influence on the thermal conductivity of quantum-confined 1D structures where the thermal conductivity should be high and not vary substantially with T, as has already been deduced through experimental verifications of thermal conductivity in MoS2 [23,27] and other layered materials [58].

Table 2: Determination of R0 and D0 which correspond to the increment rate of phonon concentration and the decay rate in phonon lifetime, respectively, with increasing T. The R0 and D0 values decreased by 3× and 2×, respectively for 1L when compared to bulk WSe2.

WSe2 layers R0 (/K) D0 (/K) [Graphic 82]
1L 0.0017 0.0031
ML 0.00487 0.0047
bulk 0.00557 0.0063

Phonon lifetime from temperature-dependent Raman analysis in WSe2

The observed broadening of the Raman peak in WSe2 with increasing T arises from damping of the excited optical phonon, and the line width Γ is inversely proportional to the phonon lifetime τ which is calculated using Equation 3. Figure 4a and 4b show the variation of τ for 1L, ML and bulk samples with increasing T for the [Graphic 83] and A1g modes, respectively. It is clear from this data that τ decreases for all sample thicknesses as T increases, similar to the dependence of τ on P, as shown by the data in Figure 2d. This decrease occurs as the rate of phonon–phonon scattering events increases with T due to the associated increase in the phonon thermal occupancy and hence, their interaction [47]. The phonon lifetime for 1L WSe2 was found to be ≈1.19 ps for the [Graphic 84] mode (Figure 4a) at T = 298 K which is consistent with the value of τ(≈1.2 ps) in [Graphic 85] mode in ME 1L WSe2 reported in our prior work [30]. It has been also found that τ is higher in sapphire substrate (≈1.4 ps) due to the decreased interfacial defect density and built-in strain on sapphire that reduces the electron–phonon interaction in the material [30]. The value of τ obtained in our study is also comparable to what has been reported for τ in 1L CVD grown WSe2 (≈0.76 ps) [30], bilayer CVD grown WSe2 (≈2.4 ps) [29], MoS2 nanoparticle (≈1.4 ps) [59], and 1L graphene (≈1.2 ps) [60]. We compared our results with the phonon lifetimes of conventional bulk 3D materials such as GaN [38,61], AlN [61], ZnO [61] and the values are provided in Table 3. Similar to τP, here we also define a define a parameter, the temperature-dependent phonon lifetime ratio ητT as given by Equation 11 below,


The values of ητT for 1L WSe2 were found to be 0.66 and 0.55 for the [Graphic 86] and A1g modes, respectively. The value of ητT is less for both Raman modes when compared with ητP, which is likely due to more scattering occurring over the wider temperature range, when compared with the laser power variation.


Figure 4: The phonon lifetime τ variation with T. Temperature-dependent phonon lifetime variation in 1L, ML and bulk WSe2 nanomembranes for the (a) [Graphic 87] and (b) A1g modes. The value of ητT, calculated using Equation 11 for 1L, ML and bulk WSe2 samples yielded 0.66, 0.51 and 0.34 for [Graphic 88] mode and 0.55, 0.50 and 0.33 for A1g mode, respectively. (c) Phonon lifetime τ for the [Graphic 89] and A1g modes as a function of T. It was found that 1.42 × τ(A1g) < [Graphic 90] < 1.90 × τ(A1g) by taking the maximum and minimum ratio of [Graphic 91] over τ(A1g). (d) The proposed 3D schematic representation of our analysis which shows the effect of external excitation, in this case the incoming laser power, in the vibrational modes in 2D layered WSe2. (e) The rate of change in τ with respect to T for the [Graphic 92] mode in 1L, ML and bulk WSe2. (f) A damped harmonic oscillator model is proposed to qualitatively represent the low, moderate and high damping scenarios in the WSe2 samples, where ξ1L, ξML and ξbulk are the damping coefficients associated with 1L, ML and bulk WSe2, respectively. From the phonon concentration and decay rate data as a function of WSe2 thickness, it was concluded that ξ1L < ξML < ξbulk as shown by the three scenarios: ξ1L(left), ξML(middle), ξbulk(right).

Table 3: Comparison of τ in 2D TMDCs and conventional bulk 3D materials. From our work, the τ for the [Graphic 93] mode was found to be almost 2× higher compared to the A1g mode for 1L, ML and bulk WSe2.

Material Mode Phonon lifetime (τ) (ps) Ref.
WSe2 (1L) [Graphic 94], A1g 1.19, 0.62 this work
WSe2 (ML) [Graphic 95], A1g 1.32, 0.65 this work
WSe2 (bulk) [Graphic 96], A1g 1.38, 0.74 this work
WSe2 (1L) [Graphic 97](SiO2/Si, sapphire) 1.2, 1.4 [30]
WSe2 (2L) E2g,[Graphic 98] 2.4, 1.3 [29]
MoS2 (nanoparticle) E1u 1.4 [59]
graphene G-Band 1.2 [60]
graphite G-Band 2.4 [60]
GaN (bulk) [Graphic 99] 2.56 [44]
GaN (bulk) [Graphic 100] 10.1, 1.4 [61]
AlN (bulk) [Graphic 101] 4.4, 0.83 [61]
ZnO (bulk) [Graphic 102] 5.9, 0.9 [61]

The values of ητT for both the [Graphic 103] and A1g modes in ML and bulk samples were found to be very similar, as is evident from the data in Figure 4a and 4b. The variation of τ with T in 1L WSe2 for both [Graphic 104] and A1g modes is depicted in Figure 4c where it was found that [Graphic 105] > τ(A1g) for all temperatures, which is quite similar to what we found in Figure 2d where we showed [Graphic 106] ≈ 2× τ(A1g) for all P. Specifically, the upper and lower bounds for τ were determined to be 1.42 × τ(A1g) < [Graphic 107] < 1.90 × τ(A1g) for all T. The decay in τ caused by the thermal excitation energy is shown schematically in Figure 4d; here the effect of external stimuli, such as the energy of incoming photons from the laser beam, on the vibrational modes of the atoms in WSe2 is illustrated. The external energy causes damping in the in-plane [Graphic 108] mode and the out-of-plane A1g mode within the WSe2 crystal, which affects the phonon frequencies and subsequently their lifetime and concentration.

Phonon lifetime and decay mechanisms in WSe2

Another perspective that elucidates the dependence of τ with T is the rate at which τ decays with T. As seen from Figure 4a and 4b, the rate of decay in τ in 1L, ML and bulk WSe2 varies widely in low T, but remains almost constant at high T. This is because the decay in τ is governed by various scattering mechanisms which are caused by the interaction of phonons with defects, boundaries, other phonons, etc. It is important to analyze the scattering as a whole through the analysis of the individual scattering events themselves, as deduced from Matthiessen’s rule [38,46,62]. As discussed earlier, the phonon linewidth broadening arises from the scattering of phonons with defects, doping, electrons, etc., and consequently, the resulting τ is an average of the contributions from each of those scattering sources. Hence, the phonon concentration, defect density, doping concentrations, etc. determine the effective strength of a scattering source and disparities between them cause the differences in the measured lifetimes in 1L, ML and bulk WSe2 which we observed in Figure 4a and 4b. Interestingly, the defect-induced scattering events will not dominate the overall scattering mechanisms as the defect concentration will not change significantly due to the stable crystal arrangement in WSe2 over the course of the measurements. However, this is not the case for other scattering events such as phonon–phonon scattering. According to the Bose–Einstein distribution, the population of phonons will change with T, which will in turn influence the number of phonon–phonon scattering events [46]. Thus, the phonon concentration increases with T (Figure 3f) and consequently τ will be dominated by phonon–phonon scattering at high T. Figures 4a and 4b also show that τ decays similarly in 1L, ML and bulk WSe2 at higher T which may be explained through a dominant source common to all samples. To confirm this, we examined the rate of change in τ with respect to T, i.e., [Graphic 109] for the Raman active [Graphic 110] mode in 1L, ML and bulk WSe2 which is shown in Figure 4e. It was found that at high T there was a convergence of [Graphic 111], indicating that the phonon–phonon scattering was dominant in all the specimens which governed the decay mechanism at high T. In contrast, at low T, τ (Figure 4a and 4b) and consequently [Graphic 112] (Figure 4e) vary widely in 1L, ML, and bulk WSe2, as nanostructure differences become the more dominant scattering source at low T. These dissimilarities in the decay in 1L, ML and bulk WSe2 at low T are characterized by fitting the results of Figure 4a, i.e., τ(T) for the [Graphic 113] mode, to the exponential decay function as noted in Equation 12,


where A2, B2 are constants having units of time and D0 is the decay rate in τ with units of K. The value of D0 (1L) was ≈0.0031/K while the D0 (bulk) was ≈0.0063/K which is almost 2× higher; the data are once again summarized in Table 2.

We conceptualize these results using a mechanical harmonic oscillator model to qualitatively compare the decay in τ with T due to damping mechanisms governed by phonon–phonon scattering as discussed above in WSe2 as a function of sample thickness. Figure 4f depicts this visually, where three different damping scenarios are illustrated. Based on our data in Figure 4a and Figure 3f, the damping is low, moderate and high in 1L (top schematic in Figure 4f), ML (middle schematic in Figure 4f) and bulk samples (bottom schematic in Figure 4f), respectively, since the rate of decay of τ increases from 1L to bulk samples. This further confirmed the result we observed earlier from the PL analysis of WSe2 with varying T (Figure 3f) where it was found that the rate of increment R0 in phonon concentration with increasing T is lowest for 1L and highest for bulk samples, suggesting that the mean free path for scattering between phonons decreases, which increases the probability of collisions when the phonon concentration increases, eventually decreasing phonon lifetime. We compared our phonon lifetime τ and temperature T to a proposed model of the standard harmonic oscillator exhibiting damping with, [Graphic 114] where x represents displacement of the oscillator, t is the time and ξ is the damping coefficient which we figuratively relate to the three scenarios for 1L, ML and bulk as ξ1L, ξML and ξbulk, respectively. Thus, from the damping in the tabulated τ values, the three damping ratios, ξ1L, ξML and ξbulk are shown in Figure 4f. It is clear from our analysis and from the values of D0 (Table 2) that ξ1L < ξML < ξbulk, i.e., the damping is least for the 1L case and it is the highest for bulk WSe2. Understanding the phonon dynamics and factors such as phonon lifetime, decay rate and nonradiative decay pathways, such as defects, will prove to be vital for truly harnessing the unique properties of WSe2 and designing high-performance electronic, optoelectronic, thermal, and thermoelectric devices from WSe2 in the future.


In conclusion, we report, the dependence of the Raman and PL of 1L, ML and bulk WSe2 by external stimuli such as laser power and temperature, from which we analyze the phonon lifetime and change in phonon concentration in WSe2. Our results show that both temperature and laser power affect the phonon dynamics in WSe2. It was found that the phonon concentration increases, followed by a decrease in the phonon lifetime, for increasing thickness of WSe2 with increasing energy of the external stimuli (temperature and laser power). We also analyzed the first-order temperature- and power-dependent [Graphic 115] and A1g coefficients and found that the coefficients increase for the [Graphic 116] mode and decrease for the A1g mode as the number of layers increases. A comparative analysis of the phonon frequencies and their relative changes from 1L to bulk WSe2 samples was conducted, which showed a frequency and intensity shift for 1L, ML, and bulk samples. The mechanical model we proposed provides a framework in which to explain the damping mechanisms in WSe2 with varying thickness and the increase in energy from external stimuli. The insights from this work reveal the importance of light–matter interactions in 2D WSe2 to alter the phonon spectrum toward the quantum-confined limit of monolayers. These findings can be broadly applied to other layered materials to help guide the design of high-performance electronic, optoelectronic, thermal, and thermoelectric devices based on WSe2 in the future.

Supporting Information

The supporting information includes: Section 1: Calculation of laser spot size, Section 2: Calculation of error bars for Raman & PL measurements, Section 3: Optimization of laser powers used in Raman and PL measurements, Section 4: Calculation of instrumental broadening and its effect on phonon lifetime analysis, Section 5: Calculation of hysteresis during cooling down and warming up processes, Section 6: Temperature-dependent Raman spectra in bulk WSe2, and Section 7: AFM, Raman and PL characterization of 1L and bulk WSe2.

Supporting Information File 1: Additional material.
Format: PDF Size: 1.0 MB Download


We greatly appreciate the funding support from the Air Force Office of Scientific Research (grant number FA9550-15-1-0200) who provided funding support that allowed us to pursue this work. A.B.K is grateful to the support received from the PACCAR Technology Institute and Endowed Professorship at the University of North Texas, Denton.


  1. Novoselov, K. S.; Jiang, D.; Schedin, F.; Booth, T. J.; Khotkevich, V. V.; Morozov, S. V.; Geim, A. K. Proc. Natl. Acad. Sci. U. S. A. 2005, 102, 10451–10453. doi:10.1073/pnas.0502848102
    Return to citation in text: [1] [2]
  2. Butler, S. Z.; Hollen, S. M.; Cao, L.; Cui, Y.; Gupta, J. A.; Gutiérrez, H. R.; Heinz, T. F.; Hong, S. S.; Huang, J.; Ismach, A. F.; Johnston-Halperin, E.; Kuno, M.; Plashnitsa, V. V.; Robinson, R. D.; Ruoff, R. S.; Salahuddin, S.; Shan, J.; Shi, L.; Spencer, M. G.; Terrones, M.; Windl, W.; Goldberger, J. E. ACS Nano 2013, 7, 2898–2926. doi:10.1021/nn400280c
    Return to citation in text: [1]
  3. Wang, Q. H.; Kalantar-Zadeh, K.; Kis, A.; Coleman, J. N.; Strano, M. S. Nat. Nanotechnol. 2012, 7, 699–712. doi:10.1038/nnano.2012.193
    Return to citation in text: [1] [2]
  4. Li, L.; Yu, Y.; Ye, G. J.; Ge, Q.; Ou, X.; Wu, H.; Feng, D.; Chen, X. H.; Zhang, Y. Nat. Nanotechnol. 2014, 9, 372–377. doi:10.1038/nnano.2014.35
    Return to citation in text: [1]
  5. Tao, L.; Cinquanta, E.; Chiappe, D.; Grazianetti, C.; Fanciulli, M.; Dubey, M.; Molle, A.; Akinwande, D. Nat. Nanotechnol. 2015, 10, 227–231. doi:10.1038/nnano.2014.325
    Return to citation in text: [1]
  6. Tan, P. H.; Han, W. P.; Zhao, W. J.; Wu, Z. H.; Chang, K.; Wang, H.; Wang, Y. F.; Bonini, N.; Marzari, N.; Pugno, N.; Savini, G.; Lombardo, A.; Ferrari, A. C. Nat. Mater. 2012, 11, 294–300. doi:10.1038/nmat3245
    Return to citation in text: [1]
  7. Zhao, Y.; Luo, X.; Li, H.; Zhang, J.; Araujo, P. T.; Gan, C. K.; Wu, J.; Zhang, H.; Quek, S. Y.; Dresselhaus, M. S.; Xiong, Q. Nano Lett. 2013, 13, 1007–1015. doi:10.1021/nl304169w
    Return to citation in text: [1] [2]
  8. Lv, R.; Robinson, J. A.; Schaak, R. E.; Sun, D.; Sun, Y.; Mallouk, T. E.; Terrones, M. Acc. Chem. Res. 2015, 48, 56–64. doi:10.1021/ar5002846
    Return to citation in text: [1]
  9. Michel, M.; Biswas, C.; Tiwary, C. S.; Saenz, G. A.; Hossain, R. F.; Ajayan, P.; Kaul, A. B. 2D Mater. 2017, 4, 025076. doi:10.1088/2053-1583/aa66ff
    Return to citation in text: [1]
  10. Saenz, G. A.; Karapetrov, G.; Curtis, J.; Kaul, A. B. Sci. Rep. 2018, 8, 1276. doi:10.1038/s41598-018-19367-1
    Return to citation in text: [1]
  11. Fadil, D.; Hossain, R. F.; Saenz, G. A.; Kaul, A. B. J. Mater. Chem. C 2017, 5, 5323–5333. doi:10.1039/c7tc01001j
    Return to citation in text: [1]
  12. Hossain, R. F.; Deaguero, I. G.; Boland, T.; Kaul, A. B. npj 2D Mater. Appl. 2017, 1, 28. doi:10.1038/s41699-017-0034-2
    Return to citation in text: [1]
  13. Kaul, A. B. J. Mater. Res. 2014, 29, 348–361. doi:10.1557/jmr.2014.6
    Return to citation in text: [1] [2]
  14. Yankowitz, M.; McKenzie, D.; LeRoy, B. J. Phys. Rev. Lett. 2015, 115, 136803. doi:10.1103/physrevlett.115.136803
    Return to citation in text: [1]
  15. Fang, H.; Chuang, S.; Chang, T. C.; Takei, K.; Takahashi, T.; Javey, A. Nano Lett. 2012, 12, 3788–3792. doi:10.1021/nl301702r
    Return to citation in text: [1]
  16. Huang, C.; Wu, S.; Sanchez, A. M.; Peters, J. J. P.; Beanland, R.; Ross, J. S.; Rivera, P.; Yao, W.; Cobden, D. H.; Xu, X. Nat. Mater. 2014, 13, 1096–1101. doi:10.1038/nmat4064
    Return to citation in text: [1]
  17. Dresselhaus, M. S.; Jorio, A.; Hofmann, M.; Dresselhaus, G.; Saito, R. Nano Lett. 2010, 10, 751–758. doi:10.1021/nl904286r
    Return to citation in text: [1]
  18. Neumann, C.; Reichardt, S.; Venezuela, P.; Drögeler, M.; Banszerus, L.; Schmitz, M.; Watanabe, K.; Taniguchi, T.; Mauri, F.; Beschoten, B.; Rotkin, S. V.; Stampfer, C. Nat. Commun. 2015, 6, 8429. doi:10.1038/ncomms9429
    Return to citation in text: [1]
  19. Venezuela, P.; Lazzeri, M.; Mauri, F. Phys. Rev. B 2011, 84, 035433. doi:10.1103/physrevb.84.035433
    Return to citation in text: [1]
  20. Lanzillo, N. A.; Glen Birdwell, A.; Amani, M.; Crowne, F. J.; Shah, P. B.; Najmaei, S.; Liu, Z.; Ajayan, P. M.; Lou, J.; Dubey, M.; Nayak, S. K.; O'Regan, T. P. Appl. Phys. Lett. 2013, 103, 093102. doi:10.1063/1.4819337
    Return to citation in text: [1] [2]
  21. Najmaei, S.; Ajayan, P. M.; Lou, J. Nanoscale 2013, 5, 9758–9763. doi:10.1039/c3nr02567e
    Return to citation in text: [1]
  22. Taube, A.; Judek, J.; Jastrzębski, C.; Duzynska, A.; Świtkowski, K.; Zdrojek, M. ACS Appl. Mater. Interfaces 2014, 6, 8959–8963. doi:10.1021/am502359k
    Return to citation in text: [1] [2] [3]
  23. Sahoo, S.; Gaur, A. P. S.; Ahmadi, M.; Guinel, M. J.-F.; Katiyar, R. S. J. Phys. Chem. C 2013, 117, 9042–9047. doi:10.1021/jp402509w
    Return to citation in text: [1] [2] [3] [4] [5]
  24. Su, L.; Zhang, Y.; Yu, Y.; Cao, L. Nanoscale 2014, 6, 4920–4927. doi:10.1039/c3nr06462j
    Return to citation in text: [1] [2] [3]
  25. M, T.; Late, D. J. ACS Appl. Mater. Interfaces 2014, 6, 1158–1163. doi:10.1021/am404847d
    Return to citation in text: [1] [2] [3]
  26. Huang, X.; Gao, Y.; Yang, T.; Ren, W.; Cheng, H.-M.; Lai, T. Sci. Rep. 2016, 6, 32236. doi:10.1038/srep32236
    Return to citation in text: [1] [2] [3] [4]
  27. Yan, R.; Simpson, J. R.; Bertolazzi, S.; Brivio, J.; Watson, M.; Wu, X.; Kis, A.; Luo, T.; Hight Walker, A. R.; Xing, H. G. ACS Nano 2014, 8, 986–993. doi:10.1021/nn405826k
    Return to citation in text: [1] [2] [3]
  28. Late, D. J.; Shirodkar, S. N.; Waghmare, U. V.; Dravid, V. P.; Rao, C. N. R. ChemPhysChem 2014, 15, 1592–1598. doi:10.1002/cphc.201400020
    Return to citation in text: [1] [2] [3] [4]
  29. Wu, W.; Morales-Acosta, M. D.; Wang, Y.; Pettes, M. T. Nano Lett. 2019, 19, 1527–1533. doi:10.1021/acs.nanolett.8b04269
    Return to citation in text: [1] [2] [3]
  30. Bandyopadhyay, A. S.; Adhikari, N.; Kaul, A. B. Chem. Mater. 2019, 31, 9861–9874. doi:10.1021/acs.chemmater.9b04086
    Return to citation in text: [1] [2] [3] [4] [5] [6] [7] [8] [9]
  31. Castellanos-Gomez, A.; Buscema, M.; Molenaar, R.; Singh, V.; Janssen, L.; van der Zant, H. S. J.; Steele, G. A. 2D Mater. 2014, 1, 011002. doi:10.1088/2053-1583/1/1/011002
    Return to citation in text: [1]
  32. Brixner, L. H. J. Inorg. Nucl. Chem. 1962, 24, 257–263. doi:10.1016/0022-1902(62)80178-x
    Return to citation in text: [1]
  33. Rasmussen, F. A.; Thygesen, K. S. J. Phys. Chem. C 2015, 119, 13169–13183. doi:10.1021/acs.jpcc.5b02950
    Return to citation in text: [1]
  34. Chen, L.; Liu, B.; Ge, M.; Ma, Y.; Abbas, A. N.; Zhou, C. ACS Nano 2015, 9, 8368–8375. doi:10.1021/acsnano.5b03043
    Return to citation in text: [1]
  35. Liu, B.; Fathi, M.; Chen, L.; Abbas, A.; Ma, Y.; Zhou, C. ACS Nano 2015, 9, 6119–6127. doi:10.1021/acsnano.5b01301
    Return to citation in text: [1] [2]
  36. Browning, R.; Kuperman, N.; Solanki, R.; Kanzyuba, V.; Rouvimov, S. Semicond. Sci. Technol. 2016, 31, 095002. doi:10.1088/0268-1242/31/9/095002
    Return to citation in text: [1]
  37. Luo, X.; Zhao, Y.; Zhang, J.; Toh, M.; Kloc, C.; Xiong, Q.; Quek, S. Y. Phys. Rev. B 2013, 88, 195313. doi:10.1103/physrevb.88.195313
    Return to citation in text: [1]
  38. Molina-Sánchez, A.; Wirtz, L. Phys. Rev. B 2011, 84, 155413. doi:10.1103/physrevb.84.155413
    Return to citation in text: [1] [2] [3]
  39. O’Brien, M.; McEvoy, N.; Hanlon, D.; Hallam, T.; Coleman, J. N.; Duesberg, G. S. Sci. Rep. 2016, 6, 19476. doi:10.1038/srep19476
    Return to citation in text: [1]
  40. Li, H.; Zhang, Q.; Yap, C. C. R.; Tay, B. K.; Edwin, T. H. T.; Olivier, A.; Baillargeat, D. Adv. Funct. Mater. 2012, 22, 1385–1390. doi:10.1002/adfm.201102111
    Return to citation in text: [1]
  41. Lee, C.; Yan, H.; Brus, L. E.; Heinz, T. F.; Hone, J.; Ryu, S. ACS Nano 2010, 5, 2695–2700. doi:10.1021/nn1003937
    Return to citation in text: [1]
  42. Xu, K.; Wang, Z.; Du, X.; Safdar, M.; Jiang, C.; He, J. Nanotechnology 2013, 24, 465705. doi:10.1088/0957-4484/24/46/465705
    Return to citation in text: [1]
  43. Zhao, W.; Ribeiro, R. M.; Toh, M.; Carvalho, A.; Kloc, C.; Castro Neto, A. H.; Eda, G. Nano Lett. 2013, 13, 5627–5634. doi:10.1021/nl403270k
    Return to citation in text: [1]
  44. Huang, J.; Yang, L.; Liu, D.; Chen, J.; Fu, Q.; Xiong, Y.; Lin, F.; Xiang, B. Nanoscale 2015, 7, 4193–4198. doi:10.1039/c4nr07045c
    Return to citation in text: [1] [2]
  45. Najmaei, S.; Liu, Z.; Ajayan, P. M.; Lou, J. Appl. Phys. Lett. 2012, 100, 013106. doi:10.1063/1.3673907
    Return to citation in text: [1]
  46. Beechem, T.; Graham, S. J. Appl. Phys. 2008, 103, 093507. doi:10.1063/1.2912819
    Return to citation in text: [1] [2] [3] [4]
  47. Cowley, R. A. Adv. Phys. 1963, 12, 421–480. doi:10.1080/00018736300101333
    Return to citation in text: [1] [2] [3]
  48. Urbach, F. Phys. Rev. 1953, 92, 1324. doi:10.1103/physrev.92.1324
    Return to citation in text: [1]
  49. Ramos, J. G.; Luzzi, R. Solid State Commun. 1974, 14, 1275–1278. doi:10.1016/0038-1098(74)90646-2
    Return to citation in text: [1] [2] [3] [4] [5]
  50. Ko, P. J.; Abderrahmane, A.; Thu, T. V.; Ortega, D.; Takamura, T.; Sandhu, A. J. Nanosci. Nanotechnol. 2015, 15, 6843–6846. doi:10.1166/jnn.2015.10348
    Return to citation in text: [1] [2] [3]
  51. Thomas, D. G.; Hopfield, J. J.; Power, M. Phys. Rev. 1960, 119, 570–574. doi:10.1103/physrev.119.570
    Return to citation in text: [1]
  52. Hart, T. R.; Aggarwal, R. L.; Lax, B. Phys. Rev. B 1970, 1, 638–642. doi:10.1103/physrevb.1.638
    Return to citation in text: [1]
  53. Beechem, T.; Graham, S.; Kearney, S. P.; Phinney, L. M.; Serrano, J. R. Rev. Sci. Instrum. 2007, 78, 061301. doi:10.1063/1.2738946
    Return to citation in text: [1]
  54. Su, L.; Yu, Y.; Cao, L.; Zhang, Y. Nano Res. 2015, 8, 2686–2697. doi:10.1007/s12274-015-0775-1
    Return to citation in text: [1]
  55. Zouboulis, E. S.; Grimsditch, M. Phys. Rev. B 1991, 43, 12490–12493. doi:10.1103/physrevb.43.12490
    Return to citation in text: [1]
  56. Wilson, J. A.; Yoffe, A. D. Adv. Phys. 1969, 18, 193–335. doi:10.1080/00018736900101307
    Return to citation in text: [1] [2] [3]
  57. Liu, Y.; Zhai, W.; Gu, P.; Zhang, T.; Wang, P.; Liu, W.; Cui, T.; Wang, Y.; Zhang, Y. Optik (Munich, Ger.) 2013, 124, 3059–3062. doi:10.1016/j.ijleo.2012.09.032
    Return to citation in text: [1]
  58. Xu, X.; Chen, J.; Li, B. J. Phys.: Condens. Matter 2016, 28, 483001. doi:10.1088/0953-8984/28/48/483001
    Return to citation in text: [1]
  59. Sun, Q.-C.; Mazumdar, D.; Yadgarov, L.; Rosentsveig, R.; Tenne, R.; Musfeldt, J. L. Nano Lett. 2013, 13, 2803–2808. doi:10.1021/nl401066e
    Return to citation in text: [1] [2]
  60. Kang, K.; Abdula, D.; Cahill, D. G.; Shim, M. Phys. Rev. B 2010, 81, 165405. doi:10.1103/physrevb.81.165405
    Return to citation in text: [1] [2] [3]
  61. Bergman, L.; Alexson, D.; Murphy, P. L.; Nemanich, R. J.; Dutta, M.; Stroscio, M. A.; Balkas, C.; Shin, H.; Davis, R. F. Phys. Rev. B 1999, 59, 12977–12982. doi:10.1103/physrevb.59.12977
    Return to citation in text: [1] [2] [3] [4] [5] [6]
  62. Gökden, S.; Tülek, R.; Teke, A.; Leach, J. H.; Fan, Q.; Xie, J.; Özgür, Ü.; Morkoç, H.; Lisesivdin, S. B.; Özbay, E. Semicond. Sci. Technol. 2010, 25, 045024. doi:10.1088/0268-1242/25/4/045024
    Return to citation in text: [1]
Other Beilstein-Institut Open Science Activities