Alloying as a Route to Monolayer Transition Metal Dichalcogenides

Apr 26, 2018 - Tan, S. J. R.; Abdelwahab, I.; Ding, Z.; Zhao, X.; Yang, T.; Loke, G. Z. J.; Lin, H.; Verzhbitskiy, I.; Poh, S. M.; Xu, H.; Nai, C. T.;...
0 downloads 4 Views 4MB Size
Article Cite This: ACS Appl. Energy Mater. XXXX, XXX, XXX−XXX

www.acsaem.org

Alloying as a Route to Monolayer Transition Metal Dichalcogenides with Improved Optoelectronic Performance: Mo(S1−xSex)2 and Mo1−yWyS2 Zhiming Shi, Qingyun Zhang, and Udo Schwingenschlögl* Physical Science and Engineering Division, King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia S Supporting Information *

ABSTRACT: On the basis of first-principles and cluster expansion calculations, we propose an effective approach to realize monolayer transition metal dichalcogenides with sizable band gaps and improved optoelectronic performance. We show that monolayer Mo(S1−xSex)2 and Mo1−yWyS2 with x = 1/3, 2/3 and y = 1/3, 1/2, 2/3 are stable according to phonon calculations and realize 1T′ or 1T″ phases. The transition barriers from the 2H phase are lower than for monolayer MoS2, implying that the 1T′ or 1T″ phases can be achieved experimentally. Furthermore, it turns out that the 1T″ monolayer alloys with x = 1/3, 2/3 and y = 1/3, 2/3 are semiconductors with band gaps larger than 1 eV, due to trimerization. The visible light absorption and carrier mobility are strongly improved as compared to 2H monolayer MoS2, MoSe2, and WS2. Thus, the 1T″ monolayer alloys have the potential to expand the applications of transition metal dichalcogenides, for example, in solar cells. KEYWORDS: transition metal dichalcogenide, monolayer, trimerization, band gap, light absorption



intercalated.19,20 Restacked 1T phases of LiMoS2 and KMoS2 have been confirmed by electron diffraction, 21,22 and stabilization by substitutional doping of Re, Tc, Mn, and Sn atoms, which serve as electron donors, has been observed.23,24 Phase engineering of ultrathin VIB TMD nanosheets from the 2H to the 1T phase turns out to be an effective method to improve the conductivity25 and therefore to enhance the material performance in electrocatalytic hydrogen evolution,26 low-resistance contact transistors,27 and electrochemical supercapacitors.28 Stability of the 1T′ and 1T″ phases under ambient conditions has been demonstrated by many techniques, including scanning transmission electron microscopy, electron diffraction, Raman, and X-ray photoelectron spectroscopy,29,30 in agreement with theoretical results.31,32 The 1T′ phase synthesized by intercalating hydrogenated Li turns out to be stable in air for more than three months.33 For nanoelectronics and photovoltaics, it is of basic importance to realize semiconductor nanostructures with continuously tunable band gaps. Alloying is an efficient approach to achieve this goal. For example, the band gap and light emission of low-dimensional ternary semiconductors can be modified gradually by changing the constituent stoichiometries.34,35 On the other hand, calculations have shown that monolayer alloys of the type Mo(S,Se,Te)2 are thermodynami-

INTRODUCTION In recent years, access to more and more two-dimensional materials, such as graphene,1 hexagonal boron nitride,2,3 transition metal dichalcogenides (TMDs),4 and black phosphorene,5,6 has led to the discovery of fundamentally new physics and to technological interest for a variety of applications. Belonging to the most studied two-dimensional materials, most TMDs undergo a transition from indirect to direct band gap semiconductors by exfoliation from bulk to monolayer, which results in great potential in electronic and optoelectronic devices.7,8 One special feature of MoS2 is polymorphism, with phases of distinct electronic characteristics. Given by the arrangement of the S atoms, monolayer MoS2 appears in 2H (trigonal prismatic D3h) and 1T (octahedral Oh) symmetries, the 2H phase being semiconducting and the 1T phase being metallic.9,10 The phases can be easily transformed into each other by intralayer atomic plane gliding, i.e., a transversal displacement of one of the S planes. For pristine monolayer MoS2, MoSe2, and WS2 the 2H phase is energetically favorable and semiconducting, while the 1T phase is metastable and metallic.11−13 Pristine monolayer PtSe2 shows the opposite properties.14 Distorted 1T phases with dimerization (1T′) and trimerization (1T″) of the metal atoms recently have been reported for the group VIB TMDs, giving rise to interesting phenomena (quantum spin Hall effect,15 ferroelectricity,16 Weyl semimetallicity17,18). The 2H phase of MoS2 has been reported to convert into a 1T phase when alkali atoms, such as Li and K, are © XXXX American Chemical Society

Received: February 24, 2018 Accepted: April 19, 2018

A

DOI: 10.1021/acsaem.8b00288 ACS Appl. Energy Mater. XXXX, XXX, XXX−XXX

Article

ACS Applied Energy Materials

Figure 1. Mixing energy of 2H bulk (a) Mo(S1−xSex)2 and (d) Mo1−yWyS2. (b, e) 2H monolayer structures obtained from the bulk ground state structures. (c, f) Distorted (1T′ or 1T″) monolayer structures obtained from the 2H monolayer structures. Yellow, cyan, orange, and blue spheres represent S, Mo, Se, and W atoms, respectively.

cally stable at room temperature for all compositions.36 Twodimensional ternary TMDs, including Mo(S1−xSex)2 and Mo1−yWyS2,37−44 are receiving much attention due to tunable band gaps38−42 and high catalytic activity for hydrogen generation.43,44 In addition, excellent contact and defect properties give rise to great potential in photovoltaic applications.45−50 The most widely used approaches to prepare 2H mono- and few-layer TMDs are mechanical exfoliation,51 chemical vapor deposition,52,53 and chemical liquid exfoliation,2,54,55 the latter providing high productivity. Regrettably, the consequences of the transition from the 2H into the 1T′ or 1T″ phase are not known. For this reason, we employ systematic first-principles calculations to investigate the structural and electronic properties of monolayer Mo(S1−xSex)2 and Mo1−yWyS2. We show that a semiconducting behavior can been realized in 1T″ monolayer alloys and that the optoelectronic properties are improved with respect to 2H monolayer MoS2, MoSe2, and WS2.

thicker than 15 Å are applied to avoid interlayer interaction. We use 0.3 Å−1 Monkhorst−Pack k-grids in the structural relaxations and 21 k-points between high symmetry points of the Brillouin zone in the band structure calculations. The force and total energy convergence thresholds are set to 0.01 eV/Å and 10−6 eV, respectively. Phase transition barriers are obtained by the climbing image nudged elastic band method60 with 10 images along the pathway from the 2H structure to the 1T′ or 1T″ structure. The forces are converged to 0.02 eV/Å. G0W0 calculations are performed on 9 × 9 × 1 Monkhorst−Pack kgrids using 450 empty bands. These k-grids have the same point density as a 15 × 15 × 1 k-grid in the case of the primitive unit cell of MoS2, for which the band gap is found to deviate only by 0.01 eV from the value obtained for a 19 × 19 × 1 k-grid. This indicates that the k-grids of the G0W0 calculations are wellconverged.





RESULTS AND DISCUSSION Cluster Expansion Calculations. Experimentally, 1T′ and 1T″ monolayer TMDs can be obtained by Li-assisted chemical exfoliation from 2H bulk TMDs. The primitive unit cell of the 2H bulk phase contains two monolayers. We built 2H bulk Mo(S1−xSex)2 and Mo1−yWyS2 alloys by cluster expansion and determine the thermodynamic ground state structures, which are displayed in Figure S1. The mixing energy is given in Figure 1. Two ground states are obtained for 2H bulk Mo(S1−xSex)2, for x = 1/3 and 2/3, see the convex hull (red lines) in Figure 1a. We note that these structures are composed of two identical monolayers with 33% and 67% Se atoms, respectively, which previously has been predicted by cluster expansion for 2H monolayer Mo(S1−xSex)2.61 Eight ground states are obtained for 2H bulk Mo1−yWyS2, for y = 1/8, 1/4, 1/3, 5/12, 7/12, 2/3, 3/4, and 7/8, see the convex hull (red lines) in Figure 1d. For y = 1/8, 5 /12, 7/12, and 7/8, the structures are composed of two different monolayers, as shown in Figure S1. Separation and subsequent relaxation result in a total of seven monolayer alloys: Mo(S1−xSex)2 with x = 1/3 and 2/3, Mo1−yWyS2 with y = 1/4, 1 /3, 1/2, 2/3, and 3/4. The 1T′ or 1T″ monolayers are obtained by transversal displacement of one of the S/Se planes. Imaginary frequencies in the phonon spectra, see the path M (1/2, 0, 0) to Γ (0, 0, 0) to K (1/3, 1/3, 0) in Figure S2, demonstrate that Mo0.25W0.75S2 and Mo0.75W0.25S2 are unstable, while the other five monolayer alloys show no imaginary frequencies (except for tiny regions around the Γ point due to artifacts of the numerical procedure used) and therefore are

METHODOLOGY Supercells of the alloys under investigation are constructed and sampled by the alloy-theoretic automated toolkit.56 To measure the quality of the cluster expansion, we use the cross-validation score57 1 M

M

∑ (EiDFT − EiCE)2 i

where EiCE is the energy predicted by cluster expansion based on M − 1 structures, and EiDFT is the corresponding energy from first-principles calculations. We use M = 55 for Mo(S1−xSex)2 and M = 61 for Mo1−yWyS2. The supercells include at most 24 atoms for Mo(S1−xSex)2 and 36 atoms for Mo1−yWyS2, with cross-validation scores of 3.4 and 0.4 meV/ atom, respectively, which means that the cluster expansions are well-converged. We find for both Mo(S1−xSex)2 and Mo1−yWyS2 that the effective cluster interaction is converged when the cluster diameter exceeds 5.5 Å. Total energies are calculated by first-principles calculations using the generalized gradient approximation in the Perdew− Burke−Ernzerhof (PBE) flavor and projector-augmented wave potentials, as implemented in the Vienna Ab-initio Simulation Package.58 The cutoff energy determining the number of plane wave basis functions is set to high a value of 500 eV. Band structures are calculated using the Heyd−Scuseria−Ernzerhof (HSE06) hybrid functional.59 In all calculations, vacuum slabs B

DOI: 10.1021/acsaem.8b00288 ACS Appl. Energy Mater. XXXX, XXX, XXX−XXX

Article

ACS Applied Energy Materials

the 1T than for the 2H phases and obtain the following transition barriers (per formula unit): 1.57 eV for MoS2, 1.35 eV for MoSe2, 1.69 eV for WS2, 1.29 eV for Mo(S0.67Se0.33)2, 1.37 eV for Mo(S0.33Se0.67)2, 1.43 eV for Mo0.67W0.33S2, 1.56 eV for Mo0.5W0.5S2, and 1.41 eV for Mo0.33W0.67S2. Thus, the transition barriers of the monolayer alloys are lower than that of monolayer MoS2; i.e., phase transition is easier than for MoS2 (for which it has been achieved experimentally). We also find that the recovery barrier back to the 2H phase is at least 0.793 eV for the considered alloys, i.e., higher than for pristine MoSe2 (0.638 eV) and WS2 (0.792 eV), indicating that spontaneous recovery is very difficult under ambient conditions. Because the phase transitions in MoS2 happen during the process of alkali atom intercalation (charge injection from the alkali atoms into the sample), we address in Figure 2b the transition barrier for negatively charged Mo(S0.67Se0.33)2 as an example (1, 2, or 3 electrons per formula unit). Results for the other monolayer alloys are qualitative similar, see Figure S3. The transition barrier decreases significantly when the negative charge increases. The fact that the 1T″ phase of Mo(S0.67Se0.33)2 is more stable than the 2H phase in the −3 charge state implies that charge injection supports the phase transition, similar to the case of monolayer MoS2.63,64 Semiconducting Monolayer Alloys. The band structures in Figure 3 display semiconducting states for 1T″ Mo(S 0 . 6 7 Se 0 . 3 3 ) 2 , Mo(S 0 . 3 3 Se 0 . 6 7 ) 2 , Mo 0 . 6 7 W 0 . 3 3 S 2 , and Mo0.33W0.67S2 with direct band gaps at the Γ point. According to Table 1, see also the band structures in Figure S4, the band

stable. In the case of Mo0.5W0.5S2, chains of MoS2 and WS2 alternate along the zigzag direction as in the 1T′ phase.15 The structures for x = 1/3 and 2/3 (y = 1/3 and 2/3) contain one Se and S (W and Mo) atom, respectively, per √3 × √3 supercell of MoS2 and MoSe2 (MoS2 and WS2). We observe that these atoms realize √3 × √3 patterns, implying that alloying provides an approach to realize 1T″ phases. The 1T″ MoS2 has been predicted to have 22 meV per formula unit lower energy than the 1T′ structure.62 We also find that relaxation of 1T′ structures results in 1T″ structures, which therefore are more stable. In addition, relaxation of the 2H and 1T″ structures of Mo(S0.67Se0.33)2 and Mo(S0.33Se0.67)2 (Figure 1b,c) shows that clustering of the S and Se atoms is energetically unfavorable. In all 2H monolayer alloys (Figure 1b,e) the bond lengths resemble those of the pristine monolayers (2.41, 2.53, and 2.41 Å for the Mo−S, Mo−Se, and W−S bonds, respectively), reflecting negligible distortions, while the 1T′ and 1T″ monolayer alloys show significant distortions (Figure 1c,f). Ab-initio molecular dynamics simulations (10 ps) have been performed at 300 K for each 1T″ monolayer alloy, employing a 4√3 × 4√3 supercell with 144 atoms, to confirm thermodynamic stability. 2H−1T Phase Transition. To evaluate whether the 1T′ and 1T″ monolayer alloys can be achieved experimentally, we calculate the phase transition barriers and compare with the pristine monolayers, see Figure 2a. We find higher energy for

Table 1. PBE and HSE06 Band Gaps (in eV) 1T′ or 1T″ phase Mo(S0.67Se0.333)2 Mo(S0.33Se0.67)2 Mo0.67W0.333S2 Mo0.5W0.5S2 Mo0.33W0.67S2 MoS2 MoSe2 WS2

Figure 2. Energy along the minimum barrier pathways from the 2H phase into the 1T′ or 1T″ phase for (a) the uncharged monolayers and (b) charged monolayer Mo(S0.67Se0.33)2.

2H phase

PBE

HSE06+SOC

PBE

0.73 0.67 0.80 0 0.82 0 0 0

1.28 1.16 1.25 0 1.25 0 0 0

1.64 1.57 1.66 1.74 1.71 1.74 1.50 1.89

Figure 3. Band structures of the monolayer alloys. Solid lines refer to the PBE level and dashed lines to the HSE06+SOC level. C

DOI: 10.1021/acsaem.8b00288 ACS Appl. Energy Mater. XXXX, XXX, XXX−XXX

Article

ACS Applied Energy Materials gaps at the PBE level are less than half as compared to the 2H phase. Since the generalized gradient approximation typically underestimates the band gaps of semiconductors, and the spin− orbit coupling (SOC) cannot be neglected, we have repeated the band structure calculations with the HSE06 hybrid functional. It turns out that the conduction bands shift to higher energy, and the band gaps grow significantly, see Table 1, being comparable to widely used silicon (1.1 eV). The effect of the SOC is weak at the valence band edges but leads to band splitting at the conduction band edges. Our results are similar to previous reports on Rashba splitting in pristine 1T″ TMDs.16,62 For 1T′ Mo0.5W0.5S2 we observe band crossing at the Fermi level between the Y and Γ points, similar to 1T′ MoS2, for which a quantum spin Hall state with large gap has been predicted.15 The mechanism of band gap opening in the 1T″ monolayer alloys turns out to be intriguing and gives rise to a new strategy for materials design. We use 1T″ Mo(S0.67Se0.33)2 as an example to reveal the mechanism. According to Figure 4, in unrelaxed

Figure 5. Densities of states projected on the Mo 4d orbitals in (a, c, e) unrelaxed 1T Mo(S0.67Se0.33)2 and (b, d, f) 1T″ Mo(S0.67Se0.33)2 with the Se atoms at the positions indicated by orange in the insets.

relative shifts of the Mo atoms along the x, y, or z axis. What is the consequence of such a shift? For instance, for a shift along the x axis the overlap between the Mo dxy, dxz orbitals and the S/Se 3p orbitals increases; i.e., the dxy and dxz states shift to higher energy with respect to the dyz states, and a band gap is opened. Indeed, the partial densities of states in Figure 5b,d,e show that the conduction band minimum is composed of dyz, dxz states, dxy, dyz states, or dxy, dxz states depending on the direction in which the respective Mo atom is shifted. We have calculated 200 random structures (2√3 × 2√3 supercells) for each Mo(S0.67Se0.33)2, Mo(S0.33Se0.67)2, Mo0.67W0.33S2, and Mo0.33W0.67S2 to confirm that they undergo trimerization and exhibit band gaps similar to our results for the √3 × √3 supercells, see Figure S5. Carrier Mobility and Absorption Properties. The sizable band gaps of the semiconducting 1T″ monolayer alloys call for evaluation of the potential in solar cells, the absorption coefficient and carrier mobility being the most important parameters. To this aim, we employ the G0W0 method combined with the Bethe−Salpeter equation to consider exciton effects and to obtain accurate frequency-dependent dielectric tensors. The absorption coefficients are shown in Figure 6a (solid lines) and are compared to monolayer MoS2, MoSe2, and WS2 (dotted lines). The visible light range (1.62− 3.11 eV) is highlighted. The first absorption peaks of monolayer MoS2, MoSe2, and WS2 are found at energies of 1.97, 1.75, and 2.20 eV, respectively, while they appear at energies of 1.03, 0.92, 1.03, and 1.05 eV for 1T″ Mo(S 0 . 6 7 Se 0 . 3 3 ) 2 , Mo(S 0 . 3 3 Se 0 . 6 7 ) 2 , Mo 0 . 6 7 W 0 . 3 3 S 2 , and Mo0.33W0.67S2. We obtain for the monolayer alloys stronger absorption than for the pristine monolayers in the low energy range from 1.62 to 2.42 eV and in the infrared range. However, in the high energy range (2.42−3.11 eV) monolayer MoS2 and WS2 show the strongest absorption. The HSE06+SOC band gaps are used to construct the band alignment diagram of Figure 6b. It turns out that a type-II band alignment (for the separation of electrons and holes) can be formed by any two monolayer alloys. The electrical transport is key for solar cell devices. Table 2 shows that the hole effective masses of the monolayer alloys are comparable to that of monolayer MoS2, whereas the electron

Figure 4. Densities of states of unrelaxed 1T MoS2 and Mo(S0.67Se0.33)2 as well as the semiconducting 1T″ monolayer alloys.

1T MoS2 the electronic states near the Fermi level originate from hybridized Mo 4d and S 3p orbitals. Both unrelaxed 1T MoS2 and Mo(S0.67Se0.33)2 show metallicity, i.e.; alloying with Se without structural relaxation does not open a band gap. However, it enhances the Mo 4d contributions with respect to the S/Se 3p contributions at the Fermi level. In unrelaxed 1T MoS2 the Mo atoms are octahedrally coordinated by S atoms so that the five 4d orbitals split into dxy, dyz, and dxz states (occupied by a total of two electrons, which leads to metallicity) and d3z2−r2, dx2−y2 states (unoccupied). Structural relaxation opens a substantial and unexpected band gap, see Figure 4. To investigate its origin, we plot orbitally decomposed densities of states in Figure 5. In the case of unrelaxed 1T Mo(S0.67Se0.33)2 the states near the Fermi level are still due to almost degenerate dxy, dyz, and dxz states, see Figure 5a,c,e. Interestingly, it turns out that the structural distortions developing around the three nonequivalent Mo atoms per supercell are distinctively different, see the insets in Figure 5b,d,f. We find that the distortions are well-approximated by D

DOI: 10.1021/acsaem.8b00288 ACS Appl. Energy Mater. XXXX, XXX, XXX−XXX

Article

ACS Applied Energy Materials

Figure 6. (a) Absorption coefficients and (b) band alignment diagram (HSE06+SOC level) with band gaps.

Table 2. Carrier Effective Masses, Deformation Potentials, Elastic Moduli, and Carrier Mobilities along the Armchair (x) and Zigzag (y) Directions at T = 300 K (PBE Level) structure

electron/hole

mx*/m0

my*/m0

Elx (eV)

Ely (eV)

Cx (J/m2)

Cy (J/m2)

μx (cm2/(V s))

μy (cm2/(V s))

Mo(S0.67Se0.33)2

electron hole electron hole electron hole electron hole

5.3 0.6 3.7 0.7 3.0 0.7 1.6 0.6

6.4 0.6 14.5 0.6 2.4 0.8 1.6 0.7

−2.5 −4.1 −2.3 −0.9 −1.6 −3.6 −0.5 −3.8

−2.2 −3.2 −2.5 −7.2 −1.9 −2.8 −2.9 −2.7

101 101 107 107 115 115 117 117

102 102 107 107 117 117 119 119

11 396 17 5311 120 344 4421 474

12 649 3 101 108 569 118 799

Mo(S0.33Se0.67)2 Mo0.67W0.33S2 Mo0.33W0.67S2

effective masses are much higher.65 The carrier mobilities can be estimated from the expression66

CONCLUSION



ASSOCIATED CONTENT

We have shown that alloying provides a feasible route to realizing semiconducting monolayer TMDs with sizable band gaps and improved optoelectronic performance. Combination of first-principles and cluster expansion calculations has enabled us to identify stable 1T′ and 1T″ monolayer alloys: Mo(S0.67Se0.33)2, Mo(S0.33Se0.67)2, Mo0.67W0.33S2, Mo0.5W0.5S2, and Mo0.33W0.67S2. The fact that the transition barriers from the 2H phase are lower than in the case of monolayer MoS2 implies that these alloys can be achieved experimentally. Most interestingly, we observe that 1T″ Mo(S0.67Se0.33)2, Mo(S0.33Se0.67)2, Mo0.67W0.33S2, and Mo0.33W0.67S2 exhibit sizable band gaps, induced by trimerization of the metal atoms, amounting to 1.23, 1.16, 1.24, and 1.16 eV at the HSE06+SOC level, respectively. Their behavior thus is fundamentally different from metallic 1T′ monolayer MoS2, MoSe2, and WS2. In addition, as compared to 2H monolayer MoS2, MoSe2, and WS2, the absorption of light is enhanced in the low energy visible region, and 1T″ Mo(S0.33Se0.67)2 and Mo0.33W0.67S2 provide very high carrier mobilities along the armchair direction. Thus, the predicted 1T″ monolayer alloys are of great interest for the fields of two-dimensional electronics and optoelectronics.

3

μ=



eℏ C kBTme*md (Eli)2

where m*e denotes the carrier effective mass in the transport direction, and md = mx*m*y , with x and y being the armchair and zigzag directions, respectively. Furthermore, Eil = ΔVi/(Δl/ l0) is the deformation potential of the valence band maximum (holes) or conduction band minimum (electrons), where ΔVi is the corresponding energy shift when the system is compressed or dilated by a distance Δl in the transport direction (equilibrium distance l0). Finally, C is the elastic modulus related to the longitudinal acoustic wave in the transport direction, derived from 2(E − E0)/S0 = C(Δl/l0)2, where E is the total energy, and S0 is the equilibrium area of the monolayer. We use for Δl/l0 five points in the range from −0.5% to 0.5% to obtain C and Eil by fitting. The calculated carrier mobilities in Table 2 reflect for 1T″ Mo(S0.33Se0.67)2 and Mo0.33W0.67S2 strong directional anisotropy. In addition, along the armchair direction, both alloys are found to outperform the carrier mobilities of monolayer MoS2 (340 cm2/(V s)), MoSe2 (240 cm2/(V s)), and WS2 (1103 cm2/(V s)),65 as well as that of phosphorene (2300 cm2/(V s)),67 which can be attributed to small deformation potentials (weak response of the states at the band edge to phonons), i.e., to reduced scattering of the carriers by phonons.

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsaem.8b00288. E

DOI: 10.1021/acsaem.8b00288 ACS Appl. Energy Mater. XXXX, XXX, XXX−XXX

Article

ACS Applied Energy Materials



Liu, Z. L.; Sun, J. T.; Lei, T.; Liu, C.; Wang, J. O.; Ibrahim, K.; Leonard, D. N.; Zhou, W.; Guo, H. M.; Wang, Y. L.; Du, S. X.; Pantelides, S. T.; Gao, H.-J. Intrinsically Patterned Two-Dimensional Materials for Selective Adsorption of Molecules and Nanoclusters. Nat. Mater. 2017, 16, 717−721. (15) Qian, X. F.; Liu, J. W.; Fu, L.; Li, J. Quantum Spin Hall Effect in Two-Dimensional Transition Metal Dichalcogenides. Science 2014, 346, 1344−1347. (16) Shirodkar, S. N.; Waghmare, U. V. Emergence of Ferroelectricity at a Metal-Semiconductor Transition in a 1T Monolayer of MoS2. Phys. Rev. Lett. 2014, 112, 157601. (17) Sun, Y.; Wu, S.-C.; Ali, M. N.; Felser, C.; Yan, B. Prediction of Weyl Semimetal in Orthorhombic MoTe2. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, 161107. (18) Qi, Y.; Naumov, P. G.; Ali, M. N.; Rajamathi, C. R.; Schnelle, W.; Barkalov, O.; Hanfland, M.; Wu, S.-C.; Shekhar, C.; Sun, Y.; Süß, V.; Schmidt, M.; Schwarz, U.; Pippel, E.; Werner, P.; Hillebrand, R.; Förster, T.; Kampert, E.; Parkin, S.; Cava, R. J.; Felser, C.; Yan, B.; Medvedev, S. A. Superconductivity in Weyl Semimetal Candidate MoTe2. Nat. Commun. 2016, 7, 11038. (19) Mattheiss, L. F. Band Structures of Transition-MetalDichalcogenide Layer Compounds. Phys. Rev. B: Condens. Matter Mater. Phys. 1973, 8, 3719−3740. (20) Py, M. A.; Haering, R. R. Structural Destabilization Induced by Lithium Intercalation in MoS2 and Related-Compounds. Can. J. Phys. 1983, 61, 76−84. (21) Heising, J.; Kanatzidis, M. G. Structure of Restacked MoS2 and WS2 Elucidated by Electron Crystallography. J. Am. Chem. Soc. 1999, 121, 638−643. (22) Wypych, F.; Solenthaler, C.; Prins, R.; Weber, T. Electron Diffraction Study of Intercalation Compounds Derived from 1T-MoS2. J. Solid State Chem. 1999, 144, 430−436. (23) Enyashin, A. N.; Yadgarov, L.; Houben, L.; Popov, I.; Weidenbach, M.; Tenne, R.; Bar-Sadan, M.; Seifert, G. New Route for Stabilization of 1T-WS2 and MoS2 Phases. J. Phys. Chem. C 2011, 115, 24586−24591. (24) Raffone, F.; Ataca, C.; Grossman, J. C.; Cicero, G. MoS2 Enhanced T-Phase Stabilization and Tunability Through Alloying. J. Phys. Chem. Lett. 2016, 7, 2304−2309. (25) Voiry, D.; Mohite, A.; Chhowalla, M. Phase Engineering of Transition Metal Dichalcogenides. Chem. Soc. Rev. 2015, 44, 2702− 2712. (26) Lukowski, M. A.; Daniel, A. S.; Meng, F.; Forticaux, A.; Li, L. S.; Jin, S. Enhanced Hydrogen Evolution Catalysis from Chemically Exfoliated Metallic MoS2 Nanosheets. J. Am. Chem. Soc. 2013, 135, 10274−10277. (27) Cho, S.; Kim, S.; Kim, J. H.; Zhao, J.; Seok, J.; Keum, D. H.; Baik, J.; Choe, D. H.; Chang, K. J.; Suenaga, K.; Kim, S. W.; Lee, Y. H.; Yang, H. Phase Patterning for Ohmic Homojunction Contact in MoTe2. Science 2015, 349, 625−628. (28) Acerce, M.; Voiry, D.; Chhowalla, M. Metallic 1T Phase MoS2 Nanosheets as Supercapacitor Electrode Materials. Nat. Nanotechnol. 2015, 10, 313−318. (29) Eda, G.; Fujita, T.; Yamaguchi, H.; Voiry, D.; Chen, M.; Chhowalla, M. Coherent Atomic and Electronic Heterostructures of Single-Layer MoS2. ACS Nano 2012, 6, 7311−7317. (30) Voiry, D.; Salehi, M.; Silva, R.; Fujita, T.; Chen, M. W.; Asefa, T.; Shenoy, V. B.; Eda, G.; Chhowalla, M. Conducting MoS2 Nanosheets as Catalysts for Hydrogen Evolution Reaction. Nano Lett. 2013, 13, 6222−6227. (31) Enyashin, A. N.; Seifert, G. Density-Functional Study of LixMoS2 Intercalates (0 ≤ x ≤ 1). Comput. Theor. Chem. 2012, 999, 13−20. (32) Kan, M.; Wang, J. Y.; Li, X. W.; Zhang, S. H.; Li, Y. W.; Kawazoe, Y.; Sun, Q.; Jena, P. Structures and Phase Transition of a MoS2 Monolayer. J. Phys. Chem. C 2014, 118, 1515−1522. (33) Tan, S. J. R.; Abdelwahab, I.; Ding, Z.; Zhao, X.; Yang, T.; Loke, G. Z. J.; Lin, H.; Verzhbitskiy, I.; Poh, S. M.; Xu, H.; Nai, C. T.; Zhou, W.; Eda, G.; Jia, B.; Loh, K. P. Chemical Stabilization of 1T′ Phase

Top and side views of ground state structures, phonon spectra of the ground state structures, minimum energy pathways, PBE band structures, and band gaps (PDF)

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Udo Schwingenschlögl: 0000-0003-4179-7231 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The research reported in this publication was supported by funding from King Abdullah University of Science and Technology (KAUST).



REFERENCES

(1) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666−669. (2) Paine, R. T.; Narula, C. K. Synthetic Routes to Boron-Nitride. Chem. Rev. 1990, 90, 73−91. (3) Coleman, J. N.; Lotya, M.; O’Neill, A.; Bergin, S. D.; King, P. J.; Khan, U.; Young, K.; Gaucher, A.; De, S.; Smith, R. J.; Shvets, I. V.; Arora, S. K.; Stanton, G.; Kim, H. Y.; Lee, K.; Kim, G. T.; Duesberg, G. S.; Hallam, T.; Boland, J. J.; Wang, J. J.; Donegan, J. F.; Grunlan, J. C.; Moriarty, G.; Shmeliov, A.; Nicholls, R. J.; Perkins, J. M.; Grieveson, E. M.; Theuwissen, K.; McComb, D. W.; Nellist, P. D.; Nicolosi, V. TwoDimensional Nanosheets Produced by Liquid Exfoliation of Layered Materials. Science 2011, 331, 568−571. (4) Novoselov, K. S.; Jiang, D.; Schedin, F.; Booth, T. J.; Khotkevich, V. V.; Morozov, S. V.; Geim, A. K. Two-Dimensional Atomic Crystals. Proc. Natl. Acad. Sci. U. S. A. 2005, 102, 10451−10453. (5) Li, L. K.; Yu, Y. J.; Ye, G. J.; Ge, Q. Q.; Ou, X. D.; Wu, H.; Feng, D. L.; Chen, X. H.; Zhang, Y. B. Black Phosphorus Field-Effect Transistors. Nat. Nanotechnol. 2014, 9, 372−377. (6) Liu, H.; Neal, A. T.; Zhu, Z.; Luo, Z.; Xu, X. F.; Tomanek, D.; Ye, P. D. Phosphorene: An Unexplored 2D Semiconductor with a High Hole Mobility. ACS Nano 2014, 8, 4033−4041. (7) Radisavljevic, B.; Radenovic, A.; Brivio, J.; Giacometti, V.; Kis, A. Single-Layer MoS2 Transistors. Nat. Nanotechnol. 2011, 6, 147−150. (8) Wang, Q. H.; Kalantar-Zadeh, K.; Kis, A.; Coleman, J. N.; Strano, M. S. Electronics and Optoelectronics of Two-Dimensional Transition Metal Dichalcogenides. Nat. Nanotechnol. 2012, 7, 699−712. (9) Wypych, F.; Schö llhorn, R. 1T-MoS2, a New Metallic Modification of Molybdenum-Disulfide. J. Chem. Soc., Chem. Commun. 1992, 19, 1386−1388. (10) Bissessur, R.; Kanatzidis, M. G.; Schindler, J. L.; Kannewurf, C. R. Encapsulation of Polymers Into MoS2 and Metal to Insulator Transition in Metastable MoS2. J. Chem. Soc., Chem. Commun. 1993, 20, 1582−1585. (11) Gupta, U.; Naidu, B. S.; Maitra, U.; Singh, A.; Shirodkar, S. N.; Waghmare, U. V.; Rao, C. N. R. Characterization of Few-Layer 1TMoSe2 and Its Superior Performance in the Visible-Light Induced Hydrogen Evolution Reaction. APL Mater. 2014, 2, 092802. (12) Chou, S. S.; Huang, Y.-K.; Kim, J.; Kaehr, B.; Foley, B. M.; Lu, P.; Dykstra, C.; Hopkins, P. E.; Brinker, C. J.; Huang, J.; Dravid, V. P. Controlling the Metal to Semiconductor Transition of MoS2 and WS2 in Solution. J. Am. Chem. Soc. 2015, 137, 1742−1745. (13) Ambrosi, A.; Sofer, Z.; Pumera, M. 2H → 1T Phase Transition and Hydrogen Evolution Activity of MoS2, MoSe2, WS2 and WSe2 Strongly Depends on the MX2 Xomposition. Chem. Commun. 2015, 51, 8450−8453. (14) Lin, X.; Lu, J. C.; Shao, Y.; Zhang, Y. Y.; Wu, X.; Pan, J. B.; Gao, L.; Zhu, S. Y.; Qian, K.; Zhang, Y. F.; Bao, D. L.; Li, L. F.; Wang, Y. Q.; F

DOI: 10.1021/acsaem.8b00288 ACS Appl. Energy Mater. XXXX, XXX, XXX−XXX

Article

ACS Applied Energy Materials Transition Metal Dichalcogenides with Giant Optical Kerr Nonlinearity. J. Am. Chem. Soc. 2017, 139, 2504−2511. (34) Bailey, R. E.; Nie, S. M. Alloyed Semiconductor Quantum Dots: Tuning the Optical Properties Without Changing the Particle Size. J. Am. Chem. Soc. 2003, 125, 7100−7106. (35) Liu, Y. K.; Zapien, J. A.; Shan, Y. Y.; Geng, C. Y.; Lee, C. S.; Lee, S. T. Wavelength-Controlled Lasing in ZnxCd1‑xS Single-Crystal Nanoribbons. Adv. Mater. 2005, 17, 1372−1377. (36) Komsa, H. P.; Krasheninnikov, A. V. Two-Dimensional Transition Metal Dichalcogenide Alloys: Stability and Electronic Properties. J. Phys. Chem. Lett. 2012, 3, 3652−3656. (37) Dumcenco, D. O.; Kobayashi, H.; Liu, Z.; Huang, Y.-S.; Suenaga, K. Visualization and Quantification of Transition Metal Atomic Mixing in Mo1‑xWxS2 Single Layers. Nat. Commun. 2013, 4, 1351. (38) Chen, Y. F.; Xi, J. Y.; Dumcenco, D. O.; Liu, Z.; Suenaga, K.; Wang, D.; Shuai, Z. G.; Huang, Y. S.; Xie, L. M. Tunable Band Gap Photoluminescence from Atomically Thin Transition-Metal Dichalcogenide Alloys. ACS Nano 2013, 7, 4610−4616. (39) Gong, Y. J.; Liu, Z.; Lupini, A. R.; Shi, G.; Lin, J. H.; Najmaei, S.; Lin, Z.; Elias, A. L.; Berkdemir, A.; You, G.; Terrones, H.; Terrones, M.; Vajtai, R.; Pantelides, S. T.; Pennycook, S. J.; Lou, J.; Zhou, W.; Ajayan, P. M. Band Gap Engineering and Layer-by-Layer Mapping of Selenium-Doped Molybdenum Disulfide. Nano Lett. 2014, 14, 442− 449. (40) Li, H.; Duan, X.; Wu, X.; Zhuang, X.; Zhou, H.; Zhang, Q.; Zhu, X.; Hu, W.; Ren, P.; Guo, P. Growth of Alloy MoS2xSe2(1−x) Nanosheets with Fully Tunable Chemical Compositions and Optical Properties. J. Am. Chem. Soc. 2014, 136, 3756−3759. (41) Mann, J.; Ma, Q.; Odenthal, P. M.; Isarraraz, M.; Le, D.; Preciado, E.; Barroso, D.; Yamaguchi, K.; von Son Palacio, G.; Nguyen, A. 2-Dimensional Transition Metal Dichalcogenides with Tunable Direct Band Gaps: MoS2(1‑x)Se2x Monolayers. Adv. Mater. 2014, 26, 1399−1404. (42) Feng, Q.; Zhu, Y.; Hong, J.; Zhang, M.; Duan, W.; Mao, N.; Wu, J.; Xu, H.; Dong, F.; Lin, F. Growth of Large-Area 2D MoS2(1‑x)Se2x Semiconductor Alloys. Adv. Mater. 2014, 26, 2648−2653. (43) Kiran, V.; Mukherjee, D.; Jenjeti, R. N.; Sampath, S. Active Guests in the MoS2/MoSe2 Host Lattice: Efficient Hydrogen Evolution Using Few-Layer Alloys of MoS2(1‑x)Se2x. Nanoscale 2014, 6, 12856−12863. (44) Gong, Q.; Cheng, L.; Liu, C.; Zhang, M.; Feng, Q.; Ye, H.; Zeng, M.; Xie, L.; Liu, Z.; Li, Y. Ultrathin MoS2(1‑x)Se2x Alloy Nanoflakes for Electrocatalytic Hydrogen Evolution Reaction. ACS Catal. 2015, 5, 2213−2219. (45) Wu, Y.-B.; Yang, W.; Wang, T.-B.; Deng, X.-H.; Liu, J.-T. Broadband Perfect Light Trapping in the Thinnest Monolayer Graphene-MoS 2 Photovoltaic Cell: The New Application of Spectrum-Splitting Structure. Sci. Rep. 2016, 6, 20955. (46) Fontana, M.; Deppe, T.; Boyd, A. K.; Rinzan, M.; Liu, A. Y.; Paranjape, M.; Barbara, P. Electron-Hole Transport and Photovoltaic Effect in Gated MoS2 Schottky Junctions. Sci. Rep. 2013, 3, 1634. (47) Pezeshki, A.; Shokouh, S. H. H.; Nazari, T.; Oh, K.; Im, S. Electric and Photovoltaic Behavior of a Few-Layer α-MoTe2/MoS2 Dichalcogenide Heterojunction. Adv. Mater. 2016, 28, 3216−3222. (48) Hao, L.; Liu, Y.; Gao, W.; Han, Z.; Xue, Q.; Zeng, H.; Wu, Z.; Zhu, J.; Zhang, W. Electrical and Photovoltaic Characteristics of MoS2/Si p−n junctions. J. Appl. Phys. 2015, 117, 114502. (49) Chow, P. K.; Jacobs-Gedrim, R. B.; Gao, J.; Lu, T.-M.; Yu, B.; Terrones, H.; Koratkar, N. Defect-Induced Photoluminescence in Monolayer Semiconducting Transition Metal Dichalcogenides. ACS Nano 2015, 9, 1520−1527. (50) Yao, J. D.; Zheng, Z. Q.; Yang, G. W. Promoting Photosensitivity and Detectivity of the Bi/Si Heterojunction Photodetector by Inserting a WS2 Layer. ACS Appl. Mater. Interfaces 2016, 8, 12915− 12924. (51) Mak, K. F.; Lee, C.; Hone, J.; Shan, J.; Heinz, T. F. Atomically Thin MoS2: A New Direct-Gap Semiconductor. Phys. Rev. Lett. 2010, 105, 136805.

(52) Lee, Y. H.; Zhang, X. Q.; Zhang, W. J.; Chang, M. T.; Lin, C. T.; Chang, K. D.; Yu, Y. C.; Wang, J. T. W.; Chang, C. S.; Li, L. J.; Lin, T. W. Synthesis of Large-Area MoS2 Atomic Layers with Chemical Vapor Deposition. Adv. Mater. 2012, 24, 2320−2325. (53) Najmaei, S.; Liu, Z.; Zhou, W.; Zou, X. L.; Shi, G.; Lei, S. D.; Yakobson, B. I.; Idrobo, J. C.; Ajayan, P. M.; Lou, J. Vapour Phase Growth and Grain Boundary Structure of Molybdenum Disulphide Atomic Layers. Nat. Mater. 2013, 12, 754−759. (54) Zeng, Z.; Yin, Z.; Huang, X.; Li, H.; He, Q.; Lu, G.; Boey, F.; Zhang, H. Single-Layer Semiconducting Nanosheets: High-Yield Preparation and Device Fabrication. Angew. Chem., Int. Ed. 2011, 50, 11093−11097. (55) Mahler, B.; Hoepfner, V.; Liao, K.; Ozin, G. A. Colloidal Synthesis of 1T-WS2 and 2H-WS2 Nanosheets: Applications for Photocatalytic Hydrogen Evolution. J. Am. Chem. Soc. 2014, 136, 14121−14127. (56) van de Walle, A.; Asta, M.; Ceder, G. The Alloy Theoretic Automated Toolkit: A User Guide. CALPHAD: Comput. Coupling Phase Diagrams Thermochem. 2002, 26, 539−553. (57) van de Walle, A.; Ceder, G. Automating First-Principles Phase Diagram Calculations. J. Phase Equilib. 2002, 23, 348−359. (58) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 1758−1775. (59) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207−8215. (60) Sheppard, D.; Xiao, P. H.; Chemelewski, W.; Johnson, D. D.; Henkelman, G. A Generalized Solid-State Nudged Elastic Band Method. J. Chem. Phys. 2012, 136, 074103. (61) Kang, J.; Tongay, S.; Li, J. B.; Wu, J. Q. Monolayer Semiconducting Transition Metal Dichalcogenide Alloys: Stability and Band Bowing. J. Appl. Phys. 2013, 113, 143703. (62) Bruyer, E.; Di Sante, D.; Barone, P.; Stroppa, A.; Whangbo, M.H.; Picozzi, S. Possibility of Combining Ferroelectricity and RashbaLike Spin Splitting in Monolayers of the 1T-Type Transition-Metal Dichalcogenides MX2 (M = Mo, W; X = S, Se, Te). Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 94, 195402. (63) Enyashin, A. N.; Yadgarov, L.; Houben, L.; Popov, I.; Weidenbach, M.; Tenne, R.; Bar-Sadan, M.; Seifert, G. New Route for Stabilization of 1T-WS2 and MoS2 Phases. J. Phys. Chem. C 2011, 115, 24586−24591. (64) Gao, G. P.; Jiao, Y.; Ma, F. X.; Jiao, Y. L.; Wacawik, E.; Du, A. J. Charge Mediated Semiconducting-to-Metallic Phase Transition in Molybdenum Disulfide Monolayer and Hydrogen Evolution Reaction in New 1T′ Phase. J. Phys. Chem. C 2015, 119, 13124−13128. (65) Zhang, W.; Huang, Z.; Zhang, W.; Li, Y. Two-Dimensional Semiconductors with Possible High Room Temperature Mobility. Nano Res. 2014, 7, 1731−1737. (66) Bruzzone, S.; Fiori, G. Ab-Initio Simulations of Deformation Potentials and Electron Mobility in Chemically Modified Graphene and Two-Dimensional Hexagonal Boron-Nitride. Appl. Phys. Lett. 2011, 99, 222108. (67) Fan, X. F.; Chang, C. H.; Zheng, W. T.; Kuo, J. L.; Singh, D. J. The Electronic Properties of Single-Layer and Multilayer MoS2 Under High Pressure. J. Phys. Chem. C 2015, 119, 10189−10196.

G

DOI: 10.1021/acsaem.8b00288 ACS Appl. Energy Mater. XXXX, XXX, XXX−XXX