The photoluminescence (PL) in monolayer transition metal dichalcogenides (TMDs) is dominated by the recombination of electrons in the conduction band with holes in the spin-orbit split valence bands, and there are two distinct emission features referred to as the A-peak (ground state exciton) and B-peak (higher spin-orbit split state). The intensity ratio of these two features varies widely, and several contradictory interpretations have been reported. In this work, we analyze the room temperature PL from MoS2, MoSe2, WS2, and WSe2 monolayers and identify the underlying cause of observed variations in emission profile. We determine that PL variations arise from differences in the non-radiative recombination associated with defect densities. Therefore, the relative intensities of the A- and B-emission features can be used to qualitatively assess the non-radiative recombination and a low B/A ratio is indicative of low defect density and high sample quality. We also performed polarization-resolved PL measurements. Emission from TMD monolayers is governed by unique optical selection rules which make them promising materials for valleytronic operations. We observe a notably higher valley polarization in the B-exciton relative to the A-exciton. The high polarization is a consequence of the shorter B-exciton lifetime resulting from rapid relaxation of excitons from the B-exciton to the A-exciton of the valence band. Our work clarifies disparities reported in the literature relating to the emission profile and provides a straightforward means to assess sample quality.

Monolayer transition metal dichalcogenides (TMDs) are a new class of materials that hold promise for electronic and optoelectronic applications.1,2 In the past several years, extensive experimental and theoretical research has provided a growing foundation of knowledge for these materials, led by the discovery that many TMDs (e.g., MoS2, MoSe2, WS2, and WSe2) transition from indirect- to direct-gap semiconductors at the monolayer limit.3,4 Structurally, they are composed of a plane of metal atoms positioned between a top and bottom chalcogen layer, arranged in a hexagonal lattice as viewed normal to the surface [Figs. 1(a) and 1(b)].

FIG. 1.

Monolayer transition metal dichalcogenides. (a) The side view and (b) plan view of a monolayer TMD. A plane of metal atoms is positioned between a top and bottom chalcogen layer and arranged in a hexagonal lattice as viewed normal to the surface. (c) A schematic diagram of the single-particle electronic band structure highlighting the spin-split A- and B-excitonic emission states and the valley dependent optical selection rules. (d) An optical micrograph of a representative monolayer MoS2 sample.

FIG. 1.

Monolayer transition metal dichalcogenides. (a) The side view and (b) plan view of a monolayer TMD. A plane of metal atoms is positioned between a top and bottom chalcogen layer and arranged in a hexagonal lattice as viewed normal to the surface. (c) A schematic diagram of the single-particle electronic band structure highlighting the spin-split A- and B-excitonic emission states and the valley dependent optical selection rules. (d) An optical micrograph of a representative monolayer MoS2 sample.

Close modal

Despite rapid progress, there are many fundamental optoelectronic aspects of monolayer TMDs that are not fully understood. The photoluminescence (PL) emission character and interpretation thereof have varied widely. As a case in point: in monolayer MoS2, some studies report only a single emission feature associated with optical transitions between the highest valence band and the conduction band at the K-points, often referred to as the ground state A-exciton emission.3,5 However, others observe a strong A-emission feature accompanied by a weaker second emission peak at 100-200 meV higher energy.6–8 This second peak, referred to as the B-exciton peak, is associated with transitions between the spin-orbit split valence band and the conduction band. Finally, a number of groups identify two distinct emission peaks of comparable intensities, associated with both A- and B-emissions.4,9,10

In addition to disparities in the number of spectral components observed, published interpretations of the relative intensities of these two features are varied and contradictory. The presence of B-peak emission has been associated with high quality samples by some,11 whereas others claim the opposite.12 It has been suggested that B-peak emission can only occur in W-based TMDs13 or when a sample is optically excited below the electronic bandgap.14 

Similar issues exist for the degree of polarization observed in the A- and B-exciton emission from these materials.15–17 The isolated monolayers have two inequivalent K-valleys at the edges of the Brillouin zone, labeled K and K′. The valence band maxima at K (K′) are populated by spin up (spin down) holes, leading to valley dependent optical selection rules15,16 [Fig. 1(c)]. A high degree of valley polarized emission is expected in the isolated monolayers.15 However, experimental values are scattered and typically much lower than predicted.18–20 

Here we address the discrepancies in interpretation of these PL emission features by analyzing a large number of different monolayer TMDs (MoS2, MoSe2, WS2, and WSe2) to better understand the conditions responsible for various emission characteristics and valley polarizations. We find that both A- and B-emission intensities can vary widely from sample-to-sample, consistent with other reports, leading to a variety of emission profiles as well as B/A intensity ratios. We show that these observed variations arise from differences in the non-radiative recombination associated with the defect density in a given sample. This relationship between PL profile and exciton dynamics provides a facile method to assess sample quality: a low B/A ratio indicates low defect density and high sample quality, whereas a large B/A ratio signals a high defect density and poor-quality material. By comparing the degree of valley polarization from A- and B-excitons for a given sample, we find a notably higher valley polarization in the B-exciton. The high polarization is a consequence of the shorter B-exciton lifetime resulting from rapid relaxation of excitons from the higher energy spin-orbit split state (B-exciton) to the ground state (A-exciton) of the valence band. This additional relaxation pathway shortens the effective lifetime for B-excitons, subsequently reducing the opportunity for intervalley scattering and enhancing the degree of valley polarization.

Monolayer transition metal dichalcogenides are synthesized via chemical vapor deposition (CVD) on SiO2/Si substrates (275 nm oxide) as detailed at the end of the manuscript. Figure 1(d) displays an optical image of MoS2 following a growth. The nucleation density and lateral size of randomly oriented MoS2 flakes vary across the substrate, providing regions of isolated monolayer islands as well as continuous film regions. Triangular growth is typical, although we have also observed hexagonal, star-like, and rounded flakes.

The as-grown TMDs are characterized at room temperature under ambient conditions. Raman, PL, and reflectance measurements are acquired using a commercial system equipped with a 50× objective so that the area analyzed is approximately a 2 µm diameter circle. The Raman spectra [Fig. 2(a)] of a representative MoS2 sample display the in-plane E12g mode at 384 cm−1 and out-of-plane A1g mode at 404 cm−1, consistent with monolayer MoS2.21,22 The normalized differential reflection spectrum (1 − R/R0, where R0 is the background signal from the bare substrate) is measured at the same location and displays clear peaks in Fig. 2(b) at 1.86 eV and 2.01 eV, corresponding to the A- and B-excitonic transitions, respectively. By contrast, only the A-peak is easily discernible in the photoluminescence spectra for this particular MoS2 monolayer [Fig. 2(c)], corresponding to radiative recombination via the lowest energy channel. Upon closer inspection and re-plotting the PL intensity on a log-scale [Fig. 2(d)], the B-peak is also evident in PL emission, albeit approximately two orders of magnitude lower than the dominant A-peak.

FIG. 2.

Characterization of monolayer MoS2. (a) Raman spectra exhibit in-plane and out-of-plane Raman modes with a peak separation Δk = 20 cm−1, characteristic of monolayer MoS2. (b) Normalized differential reflection spectra (1 − R/R0) find the A- and B-exciton at 1.86 and 2.01 eV, respectively. Photoluminescence spectra are plotted on a (c) linear and (d) logarithmic scale. Emission from both A- and B-excitons becomes clear when viewed on the log scale. The vertical dashed lines in (b) and (d) are at 1.86 and 2.01 eV. A laser excitation of 532 nm is used for Raman and PL measurements.

FIG. 2.

Characterization of monolayer MoS2. (a) Raman spectra exhibit in-plane and out-of-plane Raman modes with a peak separation Δk = 20 cm−1, characteristic of monolayer MoS2. (b) Normalized differential reflection spectra (1 − R/R0) find the A- and B-exciton at 1.86 and 2.01 eV, respectively. Photoluminescence spectra are plotted on a (c) linear and (d) logarithmic scale. Emission from both A- and B-excitons becomes clear when viewed on the log scale. The vertical dashed lines in (b) and (d) are at 1.86 and 2.01 eV. A laser excitation of 532 nm is used for Raman and PL measurements.

Close modal

Significant qualitative differences in the PL emission are reported in the literature for single monolayers, particularly regarding the presence or absence of the B-peak. To gain further insight into these reported differences and to understand the sources and implications of these variations, we examined the PL emission from additional as-grown MoS2 monolayer samples. Figure 3(a) displays spectra obtained under identical conditions from five different monolayer samples synthesized in three growth runs. The PL intensity of the A-exciton varies over an order of magnitude, from 0.113 ct/ms [Fig. 3(a), red] to 2.315 ct/ms [Fig. 3(a), light green]. Small variations in peak emission energy are also observed, typically within ∼10 meV of 1.84 eV, which are likely due to inhomogeneous distributions of strain, doping, and sample-to-substrate distances.8,23 Figure 3(b) displays the same spectra as Fig. 3(a) after normalizing to the A-peak intensity and emission energy (EA). It is apparent that the relative contribution of the B-peak changes significantly from sample-to-sample. Furthermore, there appears to be a trend between the intensity of the prominent A-peak and the appearance of a B-peak; samples having lower intensity A-peak emission exhibit a more noticeable B-peak [e.g., Figs. 3(a) and 3(b) red curve].

FIG. 3.

Photoluminescence from Mo-based monolayer TMDs. Photoluminescence from representative MoS2 monolayers is plotted in (a) ct/ms and (b) normalized to the maximum emission intensity. All spectra exhibit emission from both A- and B-excitons. The intensity for each peak is determined by 2-Lorentzian fitting, as shown in the inset of (a). (c) The B-intensity vs. A-intensity is plotted for all 24 MoS2 samples measured and is well described by the dashed line; y = 0.008 + 0.009x. (d) The B/A ratio vs. A-intensity is plotted for low intensity MoS2 samples. The horizontal line indicates 1%. B/A varies from 29% to less than 1%. [(e) and (f)] PL from representative MoSe2 is displayed in ct/ms and normalized, respectively. (g) A linear relationship is also observed in MoSe2 B-intensity vs. A-intensity and the (h) B/A ratio vs. A-intensity varies.

FIG. 3.

Photoluminescence from Mo-based monolayer TMDs. Photoluminescence from representative MoS2 monolayers is plotted in (a) ct/ms and (b) normalized to the maximum emission intensity. All spectra exhibit emission from both A- and B-excitons. The intensity for each peak is determined by 2-Lorentzian fitting, as shown in the inset of (a). (c) The B-intensity vs. A-intensity is plotted for all 24 MoS2 samples measured and is well described by the dashed line; y = 0.008 + 0.009x. (d) The B/A ratio vs. A-intensity is plotted for low intensity MoS2 samples. The horizontal line indicates 1%. B/A varies from 29% to less than 1%. [(e) and (f)] PL from representative MoSe2 is displayed in ct/ms and normalized, respectively. (g) A linear relationship is also observed in MoSe2 B-intensity vs. A-intensity and the (h) B/A ratio vs. A-intensity varies.

Close modal

To quantify individual contributions from A- and B- emissions, each PL spectrum is fit with the sum of two independent Lorentzian curves [the inset of Fig. 3(a)]. In addition to the five samples presented in Figs. 3(a) and 3(b), 19 other MoS2 monolayers, synthesized in 4 separate growth runs, were measured. Interestingly, we observe a non-zero B-peak intensity in all 24 samples. Furthermore, there is a monotonic relationship between the maximum peak intensity for the A- and B-emissions I(A) and I(B) [Fig. 3(c)] which is well-fit with a simple linear relationship: I(B) = 0.008 + 0.009*I(A). We note that, surprisingly, the fit does not pass through the origin (0,0). A straightforward extrapolation of the observed linear relationship suggests that there may be a non-zero B-peak intensity even as the A-emission becomes vanishingly small.

In the majority of samples measured, the B-emission appears to be only a small fraction of the A-intensity and one might overlook the emission from this higher energy feature. However, for samples exhibiting low-intensity emission from the A-peak, the B-peak becomes especially evident, as shown by the red and dark-blue curves in Figs. 3(a) and 3(b) corresponding to samples with the lowest values of I(A). To quantify the contribution of the B-emission feature to the photoluminescence spectra, we calculate the B/A ratio, I(B)/I(A), for each of the 24 MoS2 samples [Fig. 3(d)]. Notably, the B/A ratio exhibits considerable variation across the samples investigated and monotonically increases with decreasing I(A). We observe a ratio of 29% for the sample exhibiting the lowest A-emission intensity [Fig. 3(a), red curve]. Samples with I(A) greater than 8 ct/ms exhibit B/A <1% (horizontal line) and are not shown for clarity. The observed behavior suggests that the A- and B-emission intensities may become comparable (B/A = 1) in samples that exhibit even lower intensity A-exciton emission, and that the B-exciton emission will dominate (B/A > 1) when the A-emission is below 0.008 ct/ms for our measurement conditions.

Photoluminescence measurements were also performed on monolayers of the closely related TMD MoSe2. Emission from the dominant A-exciton occurs at lower energy (relative to MoS2 monolayers) and is observed near 1.52 eV [Fig. 3(e)]. Careful inspection reveals a small PL emission peak identified as the B-exciton ∼190 meV above the A-exciton [Fig. 3(f)], consistent with the expected valence band splitting and previous observations.24 Behavior in the selenium-based material is nearly identical to that observed in MoS2 monolayers, with all samples exhibiting a non-zero B-peak emission, a linear relationship between the A- and B-peak intensities with a non-zero intercept [Fig. 3(g)], and a variety of B/A ratios [Fig. 3(h)]. We note that a similar behavior is observed when using a different laser excitation of 488 nm (2.54 eV). Such excitation is expected to be well above the measured electronic bandgap of ∼2.2 eV for MoSe2.25 

For completeness, we also investigate the tungsten-based TMDs and present the results in Fig. 4. In both WS2 and WSe2, the general behaviors observed in MoX2 are repeated; all samples exhibit a measurable B-exciton emission with a linear correlation between A- and B-emission intensities. Additionally, the B/A ratio varies across the measured samples and is found to monotonically decrease as the A-intensity increases [Figs. 4(c) and 4(f)]. We therefore conclude that these are general characteristics for TMD monolayers.

FIG. 4.

Photoluminescence from W-based monolayer TMDs. [(a) and (d)] Photoluminescence from WS2 and WSe2 monolayers, respectively. The insets show emission from the B-exciton. For the case of WS2, Raman modes are also evident from SiO2 (black curve) as well as WS2 above 2.4 eV. [(b) and (e)] Plots detail the B-intensity vs. A-intensity for WS2 and WSe2 emission, respectively, and find a linear relationship in both instances. [(c) and (f)] Plots detail the B/A ratio vs. A-intensity for WS2 and WSe2 respectively.

FIG. 4.

Photoluminescence from W-based monolayer TMDs. [(a) and (d)] Photoluminescence from WS2 and WSe2 monolayers, respectively. The insets show emission from the B-exciton. For the case of WS2, Raman modes are also evident from SiO2 (black curve) as well as WS2 above 2.4 eV. [(b) and (e)] Plots detail the B-intensity vs. A-intensity for WS2 and WSe2 emission, respectively, and find a linear relationship in both instances. [(c) and (f)] Plots detail the B/A ratio vs. A-intensity for WS2 and WSe2 respectively.

Close modal

In order to explain the general behavior we observe in our PL data of monolayers, various radiative and non-radiative recombination pathways must be considered. The exciton lifetime (τE) is sensitive to both radiative recombination time (τR) and non-radiative recombination time (τNR) through the relationship26 

(1)

This is applicable for the A-exciton lifetime, τE,A, as well as the B-exciton lifetime, τE,B. Each material system is investigated using steady-state (cw) excitation conditions which should generate similar initial exciton populations in each sample. Therefore, the sample-to-sample differences in PL intensity (IPL) observed for a given monolayer material can be related to variations in the exciton lifetime as IPLτEτR.26 The radiative recombination time is an intrinsic property, and unlikely to vary from sample to sample in a particular TMD material at a given temperature.27 Non-radiative recombination, however, depends on a variety of factors and can vary widely. In particular, progressively shorter non-radiative recombination times are expected as the defect density increases, providing more non-radiative channels.28–30 Thus the PL intensity of both A- and B-peaks will be sensitive to changes in the density of defects mediating non-radiative recombination, with the emission intensity decreasing as τNR becomes shorter.

Excitons in the B-band have an additional available pathway, in which they scatter to the lower energy A-band. This energetically favored rapid relaxation process, having recombination time τB−A, is known to occur on the sub-picosecond timescale31 for MoS2. The additional relaxation pathway will modify Eq. (1) to

(2)

for the B-exciton. While relevant experiments in other TMD monolayers are unavailable, we surmise that the associated τB−A is comparable to or faster than in MoS2, due to the larger valence band splitting in the other TMDs. This relaxation will reduce the available exciton population in the B-channel while simultaneously increasing the A-band population, resulting in significantly different exciton lifetimes, τE,A and τE,B, and corresponding emission intensities IPLτEτR. Following Eq. (2), we see that τB−A provides an upper bound of ∼1 ps for the B-exciton lifetime, which is significantly shorter than the expected lifetime of ∼800 ps for the A-exciton in MoS2 at room temperature.27 We note that the radiative lifetime is strongly temperature dependent in monolayer TMDs, with lifetime decreasing as temperature decreases.32,33 In this manuscript, all measurements and discussion focus on room temperature conditions.

In the ideal case, where only radiative recombination and the τB−A recombination occur, a significantly higher PL intensity is expected for A-excitons than B-excitons. However, in reality, non-radiative pathways are common, arising from factors such as defects, charge trapping, and exciton-exciton anihiliation,34–36 and will modify both τE,A and τE,B. The subsequent impact of these additional non-radiative pathways is very different for A- and B-excitons, in that the effect on B-excitons will be minimized due to the very rapid intraband relaxation provided by the τB−A pathway. According to Eq. (1), an A-exciton having radiative recombination time of 800 ps and non-radiative recombination of 10 ps would result in an effective A-exciton lifetime, τE,A of ∼10 ps. By contrast, a B-exciton radiative lifetime of 800 ps combined with the τB−A = 1 ps would first result in an effective B-exciton lifetime of ∼1 ps, which is then only moderately impacted to τB ∼ 900 fs by the addition of the 10 ps non-radiative recombination.

Importantly, this indicates that the relative intensities of the A- and B-exciton emission features of monolayers can be used to qualitatively assess the non-radiative recombination, and thus the quality of the sample. A clearly dominant A-emission indicates a long A-exciton lifetime, relative to the B-exciton lifetime, and a low density of non-radiative defects. By contrast, comparable intensities of A-and B-peaks indicate that non-radiative recombination due to a high density of defects has significantly reduced the effective A-exciton lifetime such that it is nearly as rapid as that of the B-to-A relaxation pathway, τB−A. Thus, the B/A peak ratio directly reflects the quality of the sample, with low B/A indicating a low density of non-radiative defects and long A-exciton lifetime. Note that this information can be obtained using cw lasers and steady state conditions, avoiding the more complex apparatus of time-resolved measurement as well as the use of pulsed lasers, which are known to alter the emission spectra and potentially damage the sample due to transient high power densities.37 Additional Raman spectra of representative MoS2 and WS2 samples are presented in the supplementary material and provide an independent assessment of sample quality. The in-plane E12g Raman mode is associated with crystalline quality, with sharper FWHM indicative of lower defect density.7,9,38 Raman spectra exhibit a considerably narrower E12g mode in samples exhibiting lower B/A ratio, providing further support for our assertions that the smaller B/A ratio indicates a lower defect density.

The short lifetime of the B-exciton is ideal for valley polarization. The degree of circularly polarized emission, Pcirc, is related to both exciton lifetime and valley lifetime through the relationship17,20,39

(3)

where P0 is the initial polarization and τs is the valley relaxation time. As evident in Eq. (3), for a fixed P0, higher Pcirc can be obtained either by increasing valley lifetime or by decreasing exciton lifetime. In a TMD monolayer such as WS2 where the A-emission is dominant, τE,B is shorter than τE,A, as discussed previously. Consequently, emission from the B-exciton should exhibit a higher degree of polarization and provides an internal check of our arguments above.

We test our assertions in monolayer WS2 measured at room temperature. The normalized differential reflection spectrum identifies the A- and B-excitonic features at 1.96 eV and 2.35 eV, respectively [Fig. 5(a), dashed lines]. The degree of valley polarization is measured using 588 nm (2.11 eV) excitation for the A-peak and 488 nm (2.54 eV) excitation for the B-peak measurement to be near resonance for each feature, with the excitation energies indicated by the orange and green lines [Fig. 5(a)], respectively. These excitation conditions ensure that the energy separation between laser excitation and PL emission energy is nearly equal for the A- and B-peak measurements, suppressing any dependence on excitation energy.18,20,28

FIG. 5.

Valley polarization from A- and B-peak emissions in WS2 monolayers. (a) Normalized differential reflection spectrum of a monolayer WS2 region. The dashed lines indicate the A- and B-excitonic features, and the colored lines indicate the laser excitation used for polarization measurements of the A-peak (588 nm, orange line) and B-peak (488 nm, green line). (b) σ+ (red) and σ (blue) PL components for 588 nm excitation. (c) σ+ and σ PL components for 488 nm excitation. The inset of (c) highlights emission from the B-peak. The narrow peaks above 2.4 eV are Raman peaks arising from SiO2 as well as WS2. The valley polarization from the B-peak emission is more than a factor of three greater than the valley polarization from the A-peak measured at the the same location.

FIG. 5.

Valley polarization from A- and B-peak emissions in WS2 monolayers. (a) Normalized differential reflection spectrum of a monolayer WS2 region. The dashed lines indicate the A- and B-excitonic features, and the colored lines indicate the laser excitation used for polarization measurements of the A-peak (588 nm, orange line) and B-peak (488 nm, green line). (b) σ+ (red) and σ (blue) PL components for 588 nm excitation. (c) σ+ and σ PL components for 488 nm excitation. The inset of (c) highlights emission from the B-peak. The narrow peaks above 2.4 eV are Raman peaks arising from SiO2 as well as WS2. The valley polarization from the B-peak emission is more than a factor of three greater than the valley polarization from the A-peak measured at the the same location.

Close modal

The sample is excited with σ+ helicity light, and the emission is analyzed for σ+ and σ helicity. The subsequent degree of circular polarization is computed as17 

(4)

where I(σ+/−) are the polarization resolved PL intensities. As evident in Fig. 5(b), a Pcirc of 7% is measured at the A-emission peak for 588 nm excitation, comparable to previously reported values under similar conditions.18,39 However, the circularly polarized emission at the B-peak on the same sample [Fig. 5(c), inset 2.36 eV, dashed line] is significantly enhanced to a value of 25% for 488 nm excitation. The narrow peaks present above 2.4 eV arise from Raman features and do not affect the measured valley polarization at the B- emission energy. The considerably higher Pcirc for B-emission than A-emission can be attributed to the shorter τE,B compared to τE,A as discussed previously in our model.

We conclude that the ratio of the B- and A-exciton intensities reflects the density of non-radiative defects, thereby providing a qualitative measure of sample quality. This information helps clarify why significant variations in these PL components have been reported in the literature and enhances our fundamental understanding of excitonic dynamics and valley polarization for both A- and B-emissions in monolayer transition metal dichalcogenides. We identify non-zero PL emission from the B-exciton in all TMD monolayers investigated (MoS2, MoSe2, WS2, and WSe2). Sample-to-sample differences in non-radiative recombination, stemming from differences in defect densities, lead to variations in the A- and B-emission intensities and emission profiles. The variations in the photoluminescence (in the form of B/A intensity ratio) can be utilized to assess sample quality, with a large B/A ratio indicative of high defect density and low sample quality. While all samples exhibit B-exciton emission, the intensity can be orders of magnitude lower than the ground-state A-emission, making identification challenging in high-quality samples. Emission from the B-exciton exhibits an increased degree of valley polarization, relative to the A-exciton. The enhanced polarization results from the presence of an additional relaxation pathway, τB−A, which is available only to the B-exciton population and shortens the lifetime τΕ,B, subsequently reducing the opportunity for intervalley scattering.

The final paragraphs detail the growth techniques and characterization conditions used in this study. Synthesis of monolayer TMDs is performed at ambient pressure on SiO2/Si (275 nm oxide) substrates. A dedicated 2″ quartz tube furnace is used for each material to prevent cross contamination. Prior to use, all SiO2/Si substrates undergo a standard cleaning procedure consisting of: (i) ultrasonic in acetone; (ii) ultrasonic in isopropyl alcohol; (iii) submersion in Piranha etch (3:1 mixture of H2SO4:H2O2); and (iv) thorough rinsing in de-ionized water. Although exact precursor amounts, gas flow, and growth temperatures are optimized for each material, a similar methodology is utilized to grow the MoS2, MoSe2, WS2, and WSe2 monolayers. A quartz boat containing metal oxide powder (Alfa Aesar 10812 MoO3 or Alfa Aesar 13398 WO3) is positioned at the center of the tube furnace. Two SiO2/Si wafers are positioned face-down, directly above the oxide precursor. The upstream wafer contains perylene-3,4,9,10-tetracarboxylic acid tetrapotassium salt (PTAS) seeding molecules, while the downstream substrate is untreated. The hexagonal PTAS molecules are carried downstream to the untreated substrate and promote lateral growth of the TMD materials.40 Growth occurs on both upstream (PTAS coated) and downstream (untreated) substrates. A separate quartz boat containing the chalcogen powder (Alfa Aesar 36208 Se or Alfa Aesar 10755 S) is placed upstream, outside the furnace-heating zone. High purity argon (99.999%, Airgas) flows through the furnace as the temperature is ramped from room temerature to the target temperature. Upon reaching the target temperature, H2 (99.999%, AirGas) is added to the Ar flow for WS2, WSe2, and MoSe2 growth and maintained throughout the 10 min soak and subsequent cooling. During MoS2 growth, Ar is the only carrier gas used and continues to flow during the 10 min soak and cooling. Growth temperatures are 825 °C, 825 °C, 625 °C, and 700 °C for WS2, WSe2, MoS2, and MoSe2, respectively. Additional growth details can be found in our previous studies.41,42

Optical spectroscopy measurements are performed at room temperature in an atmosphere using a Horiba LabRam confocal Raman/PL microscope system. A 50× objective is used for all Raman and PL measurements, which focuses the laser to a diameter of ∼2 µm at the sample surface. The laser power is measured at the sample position and is optimized for each material system to achieve high signal-to-noise ratio. Raman and PL spectra presented in Fig. 2 are acquired using 532 nm excitation of 20 µW. Photoluminescence spectra presented in Figs. 3 and 4 are acquired using 488 nm excitation with 3 µW power for WS2, 488 nm excitation with 6 µW power for WSe2, 532 nm excitation with 20 µW power for MoS2, and 532 nm excitation with 37 µW power for MoSe2. Photoluminescence spectra presented in Figs. 5(b) and 5(c) are acquired using 588 nm excitation with 10 µW and 488 nm excitation with 10 µW, respectively.

See supplementary material for additional Raman and PL characterization.

This research was performed while S.V.S. held a National Research Council fellowship at NRL. This work was supported by core programs at NRL and the NRL Nanoscience Institute.

1.
W.
Choi
,
N.
Choudhary
,
G. H.
Han
,
J.
Park
,
D.
Akinwande
, and
Y. H.
Lee
,
Mater. Today
20
,
116
(
2017
).
2.
Z.
Lin
,
A.
McCreary
,
N.
Briggs
,
S.
Subramanian
,
K.
Zhang
,
Y.
Sun
,
X.
Li
,
N. J.
Borys
,
H.
Yuan
,
S. K.
Fullerton-Shirey
,
A.
Chernikov
,
H.
Zhao
,
S.
McDonnell
,
A. M.
Lindenberg
,
K.
Xiao
,
B. J.
LeRoy
,
M.
Drndić
,
J. C. M.
Hwang
,
J.
Park
,
M.
Chhowalla
,
R. E.
Schaak
,
A.
Javey
,
M. C.
Hersam
,
J.
Robinson
, and
M.
Terrones
,
2D Mater.
3
,
042001
(
2016
).
3.
K. F.
Mak
,
C.
Lee
,
J.
Hone
,
J.
Shan
, and
T. F.
Heinz
,
Phys. Rev. Lett.
105
,
136805
(
2010
).
4.
A.
Splendiani
,
L.
Sun
,
Y.
Zhang
,
T.
Li
,
J.
Kim
,
C.-Y.
Chim
,
G.
Galli
, and
F.
Wang
,
Nano Lett.
10
,
1271
(
2010
).
5.
W.
Zhao
,
Z.
Ghorannevis
,
L.
Chu
,
M.
Toh
,
C.
Kloc
,
P.-H.
Tan
, and
G.
Eda
,
ACS Nano
7
,
791
(
2013
).
6.
X.
Dou
,
K.
Ding
,
D.
Jiang
,
X.
Fan
, and
B.
Sun
,
ACS Nano
10
,
1619
(
2016
).
7.
Y.
Xie
,
Z.
Wang
,
Y.
Zhan
,
P.
Zhang
,
R.
Wu
,
T.
Jiang
,
S.
Wu
,
H.
Wang
,
Y.
Zhao
,
N.
Tang
, and
X.
Ma
,
Nanotechnology
28
,
084001
(
2017
).
8.
Z.
Liu
,
M.
Amani
,
S.
Najmaei
,
Q.
Xu
,
X.
Zou
,
W.
Zhou
,
T.
Yu
,
C.
Qiu
,
A. G.
Birdwell
,
F. J.
Crowne
,
R.
Vajtai
,
B. I.
Yakobson
,
Z.
Xia
,
M.
Dubey
,
P. M.
Ajayan
, and
J.
Lou
,
Nat. Commun.
5
,
5246
(
2014
).
9.
Y.
Yu
,
C.
Li
,
Y.
Liu
,
L.
Su
,
Y.
Zhang
, and
L.
Cao
,
Sci. Rep.
3
,
1866
(
2013
).
10.
Y.-C.
Lin
,
W.
Zhang
,
J.-K.
Huang
,
K.-K.
Liu
,
Y.-H.
Lee
,
C.-T.
Liang
,
C.-W.
Chu
, and
L.-J.
Li
,
Nanoscale
4
,
6637
(
2012
).
11.
C.
Cong
,
J.
Shang
,
X.
Wu
,
B.
Cao
,
N.
Peimyoo
,
C.
Qiu
,
L.
Sun
, and
T.
Yu
,
Adv. Opt. Mater.
2
,
131
(
2014
).
12.
Q.
Ji
,
Y.
Zhang
,
T.
Gao
,
Y.
Zhang
,
D.
Ma
,
M.
Liu
,
Y.
Chen
,
X.
Qiao
,
P.-H.
Tan
,
M.
Kan
,
J.
Feng
,
Q.
Sun
, and
Z.
Liu
,
Nano Lett.
13
,
3870
(
2013
).
13.
H.
Zeng
,
G.-B.
Liu
,
J.
Dai
,
Y.
Yan
,
B.
Zhu
,
R.
He
,
L.
Xie
,
S.
Xu
,
X.
Chen
,
W.
Yao
, and
X.
Cui
,
Sci. Rep.
3
,
1608
(
2013
).
14.
A.
Steinhoff
,
J.-H.
Kim
,
F.
Jahnke
,
M.
Rösner
,
D.-S.
Kim
,
C.
Lee
,
G. H.
Han
,
M. S.
Jeong
,
T. O.
Wehling
, and
C.
Gies
,
Nano Lett.
15
,
6841
(
2015
).
15.
D.
Xiao
,
G.-B.
Liu
,
W.
Feng
,
X.
Xu
, and
W.
Yao
,
Phys. Rev. Lett.
108
,
196802
(
2012
).
16.
T.
Cao
,
G.
Wang
,
W.
Han
,
H.
Ye
,
C.
Zhu
,
J.
Shi
,
Q.
Niu
,
P.
Tan
,
E.
Wang
,
B.
Liu
, and
J.
Feng
,
Nat. Commun.
3
,
887
(
2012
).
17.
K. F.
Mak
,
K.
He
,
J.
Shan
, and
T. F.
Heinz
,
Nat. Nanotechnol.
7
,
494
(
2012
).
18.
A. T.
Hanbicki
,
G.
Kioseoglou
,
M.
Currie
,
C. S.
Hellberg
,
K. M.
McCreary
,
A. L.
Friedman
, and
B. T.
Jonker
,
Sci. Rep.
6
,
18885
(
2016
).
19.
G.
Wang
,
E.
Palleau
,
T.
Amand
,
S.
Tongay
,
X.
Marie
, and
B.
Urbaszek
,
Appl. Phys. Lett.
106
,
112101
(
2015
).
20.
G.
Kioseoglou
,
A. T.
Hanbicki
,
M.
Currie
,
A. L.
Friedman
,
D.
Gunlycke
, and
B. T.
Jonker
,
Appl. Phys. Lett.
101
,
221907
(
2012
).
21.
P.
Tonndorf
,
R.
Schmidt
,
P.
Böttger
,
X.
Zhang
,
J.
Börner
,
A.
Liebig
,
M.
Albrecht
,
C.
Kloc
,
O.
Gordan
,
D. R. T.
Zahn
,
S.
Michaelis de Vasconcellos
, and
R.
Bratschitsch
,
Opt. Express
21
,
4908
(
2013
).
22.
Y.
Zhan
,
Z.
Liu
,
S.
Najmaei
,
P. M.
Ajayan
, and
J.
Lou
,
Small
8
,
966
(
2012
).
23.
D.-H.
Lien
,
J. S.
Kang
,
M.
Amani
,
K.
Chen
,
M.
Tosun
,
H.-P.
Wang
,
T.
Roy
,
M. S.
Eggleston
,
M. C.
Wu
,
M.
Dubey
,
S.-C.
Lee
,
J.-H.
He
, and
A.
Javey
,
Nano Lett.
15
,
1356
(
2015
).
24.
X.
Lu
,
M. I. B.
Utama
,
J.
Lin
,
X.
Gong
,
J.
Zhang
,
Y.
Zhao
,
S. T.
Pantelides
,
J.
Wang
,
Z.
Dong
,
Z.
Liu
,
W.
Zhou
, and
Q.
Xiong
,
Nano Lett.
14
,
2419
(
2014
).
25.
M. M.
Ugeda
,
A. J.
Bradley
,
S.-F.
Shi
,
F. H.
da Jornada
,
Y.
Zhang
,
D. Y.
Qiu
,
W.
Ruan
,
S.-K.
Mo
,
Z.
Hussain
,
Z.-X.
Shen
,
F.
Wang
,
S. G.
Louie
, and
M. F.
Crommie
,
Nat. Mater.
13
,
1091
(
2014
).
26.
L.
Yuan
and
L.
Huang
,
Nanoscale
7
,
7402
(
2015
).
27.
M.
Palummo
,
M.
Bernardi
, and
J. C.
Grossman
,
Nano Lett.
15
,
2794
(
2015
).
28.
K. M.
McCreary
,
M.
Currie
,
A. T.
Hanbicki
,
H.-J.
Chuang
, and
B. T.
Jonker
,
ACS Nano
11
,
7988
(
2017
).
29.
M.
Amani
,
D.-H.
Lien
,
D.
Kiriya
,
J.
Xiao
,
A.
Azcatl
,
J.
Noh
,
S. R.
Madhvapathy
,
R.
Addou
,
S.
Kc
,
M.
Dubey
,
K.
Cho
,
R. M.
Wallace
,
S.-C.
Lee
,
J.-H.
He
,
J. W.
Ager
,
X.
Zhang
,
E.
Yablonovitch
, and
A.
Javey
,
Science
350
,
1065
(
2015
).
30.
M. R.
Rosenberger
,
H.-J.
Chuang
,
K. M.
McCreary
,
C. H.
Li
, and
B. T.
Jonker
,
ACS Nano
12
,
1793
(
2018
).
31.
S.
Cha
,
J. H.
Sung
,
S.
Sim
,
J.
Park
,
H.
Heo
,
M.-H.
Jo
, and
H.
Choi
,
Nat. Commun.
7
,
10768
(
2016
).
32.
G.
Wang
,
A.
Chernikov
,
M. M.
Glazov
,
T. F.
Heinz
,
X.
Marie
,
T.
Amand
, and
B.
Urbaszek
,
Rev. Mod. Phys.
90
,
021001
(
2018
).
33.
G.
Moody
,
C.
Kavir Dass
,
K.
Hao
,
C.-H.
Chen
,
L.-J.
Li
,
A.
Singh
,
K.
Tran
,
G.
Clark
,
X.
Xu
,
G.
Berghäuser
,
E.
Malic
,
A.
Knorr
, and
X.
Li
,
Nat. Commun.
6
,
8315
(
2015
).
34.
H.
Shi
,
R.
Yan
,
S.
Bertolazzi
,
J.
Brivio
,
B.
Gao
,
A.
Kis
,
D.
Jena
,
H. G.
Xing
, and
L.
Huang
,
ACS Nano
7
,
1072
(
2013
).
35.
P. D.
Cunningham
,
K. M.
McCreary
,
A. T.
Hanbicki
,
M.
Currie
,
B. T.
Jonker
, and
L. M.
Hayden
,
J. Phys. Chem. C
120
,
5819
(
2016
).
36.
P. D.
Cunningham
,
K. M.
McCreary
, and
B. T.
Jonker
,
J. Phys. Chem. Lett.
7
,
5242
(
2016
).
37.
A. T.
Hanbicki
,
M.
Currie
,
G.
Kioseoglou
,
C. S.
Hellberg
,
A. L.
Friedman
, and
B. T.
Jonker
,
Nanoscale
9
,
17422
(
2017
).
38.
M. S.
Kim
,
S. J.
Yun
,
Y.
Lee
,
C.
Seo
,
G. H.
Han
,
K. K.
Kim
,
Y. H.
Lee
, and
J.
Kim
,
ACS Nano
10
,
2399
(
2016
).
39.
B.
Zhu
,
H.
Zeng
,
J.
Dai
,
Z.
Gong
, and
X.
Cui
,
Proc. Natl. Acad. Sci. U. S. A.
111
,
11606
(
2014
).
40.
X.
Ling
,
Y.-H.
Lee
,
Y.
Lin
,
W.
Fang
,
L.
Yu
,
M. S.
Dresselhaus
, and
J.
Kong
,
Nano Lett.
14
,
464
(
2014
).
41.
O. M. J.
van’t Erve
,
A. T.
Hanbicki
,
A. L.
Friedman
,
K. M.
McCreary
,
E.
Cobas
,
C. H.
Li
,
J. T.
Robinson
, and
B. T.
Jonker
,
J. Mater. Res.
31
,
845
(
2016
).
42.
K. M.
McCreary
,
A. T.
Hanbicki
,
G. G.
Jernigan
,
J. C.
Culbertson
, and
B. T.
Jonker
,
Sci. Rep.
6
,
19159
(
2016
).

Supplementary Material