Oxidation-Induced Topological Phase Transition in Monolayer 1T

Aug 6, 2018 - This website uses cookies to improve your user experience. By continuing to use the site, you are accepting our use of cookies. Read the...
1 downloads 0 Views 3MB Size
Letter Cite This: J. Phys. Chem. Lett. 2018, 9, 4783−4788

pubs.acs.org/JPCL

Oxidation-Induced Topological Phase Transition in Monolayer 1T′WTe2 Jiali Yang,†,‡ Yuanjun Jin,‡ Wangping Xu,‡ Baobing Zheng,*,‡,§,∥ Rui Wang,*,†,‡ and Hu Xu*,‡ †

Institute for Structure and Function and Department of Physics, Chongqing University, Chongqing 400044, China Department of Physics and Institute for Quantum Science and Engineering, Southern University of Science and Technology, Shenzhen 518055, China § College of Physics and Optoelectronic Technology, Nonlinear Research Institute, Baoji University of Arts and Sciences, Baoji 721016, China ∥ School of Physics and Technology, Wuhan University, Wuhan 430072, China

J. Phys. Chem. Lett. 2018.9:4783-4788. Downloaded from pubs.acs.org by UNIV OF BRITISH COLUMBIA on 01/03/19. For personal use only.



S Supporting Information *

ABSTRACT: Monolayer (ML) tungsten ditelluride (WTe2) is a well-known quantum spin Hall (QSH) insulator with topologically protected gapless edge states, thus promising dissipationless electronic devices. However, experimental findings exhibit the fast oxidation of ML WTe2 in ambient conditions. To reveal the changes of topological properties of WTe2 arising from oxidation, we systematically study the surface oxidation reaction of ML 1T′-WTe2 using first-principles calculations. The calculated results indicate that the fast oxidation of WTe2 originates from the existence of H2O in air, which significantly promotes the oxidation of ML 1T′-WTe2. More importantly, this low-coverage oxidized WTe2 loses its topological features and is changed into a trivial insulator. Furthermore, we propose a fully oxidized ML WTe2 that can still possess the QSH insulator states. The topological phase transition induced by oxidation provides exotic insight into understanding the topological features of layered transition-metal dichalcogenide materials.

S

conditions because of their large exposure area.29 The surface of the WTe2 film is severely corroded with time when it is exposed to ambient air conditions, especially for single-layer WTe2.30−32 The fast degradation of the WTe2 film in air stems from the surface oxidation reaction only on top of the degraded surface.30 Once the surface layer of the WTe2 film is oxidized, the oxidation products effectively protect the inner layer of the WTe2 film from further atmospheric corrosion, resulting in the self-limiting and surface-saturating behavior of the WTe2 film. Fast surface oxidation often has detrimental effects on the electronic properties of the ML 1T′-WTe2. Therefore, a crucial issue is how to obtain its predicted electronic phases when we fabricate the corresponding nanoscale devices in a single layer. Thus, the role of surface oxidation on the topological properties of ML 1T′-WTe2 must be further elaborated. In this work, we focus on the influence of surface oxidation on the topological properties of ML 1T′-WTe2. Low-coverage surface oxidation is studied according to the reaction pathway calculations. Oxidation together with the H2O molecule is taken into account in light of the existence of water in air. We essentially elucidate the changes in the electronic properties of

ince the experimental discovery of graphene, many attempts1−6 have been made to investigate two-dimensional (2D) materials. Layered transition-metal dichalcogenide (TMD) materials,7−9 a class of typical 2D materials, have attracted extensive attention because of their unique sandwich structures and unexpected properties. In particular, their intriguing electronic properties with potential applications in practical devices have been intensively reported in recent years.10−15 The typical bulk phases of TMDs are referred to as 2H, 1T, and 1T′. Most TMD materials crystallize in a semiconducting 2H structure.16 Contrary to common 2H TMDs, the most energetically stable structure of bulk WTe2 is the 1T′ phase,16,17 which can be considered a distorted 1T phase. The intriguing physical properties of 1T′-WTe2 have prompted researchers to focus on this material extensively, such as in studies of its topological properties,18−22 superconductivity,23 extraordinary magnetoresistance effects,24,25 thermoelectric properties,26 etc. Bulk WTe2 is a type-II Weyl semimetal with topologically nontrivial Fermi arcs that connect the electron and the hole pockets.18 Theoretical and experimental studies demonstrate that monolayer (ML) WTe2 is a quantum spin Hall (QSH) insulator27,28 that has a sizable bulk band gap of 55 meV and conductive edge states. The topologically protected edge states make the 1T′-WTe2 a promising material for fabricating spintronic devices due to its low dissipation. However, fewlayer 1T′-WTe2 films are relatively unstable under ambient © 2018 American Chemical Society

Received: June 25, 2018 Accepted: August 6, 2018 Published: August 6, 2018 4783

DOI: 10.1021/acs.jpclett.8b01999 J. Phys. Chem. Lett. 2018, 9, 4783−4788

Letter

The Journal of Physical Chemistry Letters

To understand the effects of oxidation on topological properties and for the sake of simplicity, we start with a discussion of the low-coverage oxidation of ML 1T′-WTe2. Several possible dissociation configurations (see details in Figure S1 of the Supporting Information (SI)) of O2 on ML 1T′-WTe2 are considered, and the corresponding adsorption energies in the dissociation form are also listed. We find that the energetically favorable configuration (Figure S1f), in which two O atoms are completely inserted into the structure of ML WTe2 and each O atom coordinates with three adjacent atoms, has the lowest adsorption energy of −1.982 eV per O2. This low adsorption energy implies that 1T′-WTe2 is in favor of the surface oxidation, which supports fast degradation on top of the 1T′-WTe2 film observed in the experiment.31 To elucidate the oxidation of ML 1T′-WTe2, we chose the energetically favorable configuration to simulate surface oxidation using the CI-NEB method. The oxidative decomposition process on the ML WTe2 surface, together with the reaction energy barrier, is shown in Figure 2. During the

oxidized ML 1T′-WTe2 during the surface reaction, including the band structure, density of states (DOS), and differential charge density. In addition, to study the topological properties of ML 1T′-WTe2 in a fully oxidized condition, we propose an energetically favorable rectangular structure of 2D WTe2O. The calculated Z2 topological invariant and nontrivial edge states suggest that this 2D WTe2O can continue to possess the QSH insulator states. We performed first-principles calculations as implemented in the Vienna ab initio simulation package (VASP).33,34 The exchange−correlation energy was described by the generalized gradient approximation (GGA) with the Perdew−Burke− Ernzerhof functional.35 The Kohn−Sham equations were solved self-consistently within the projector augmented wave method.36,37 The plane wave basis set was truncated with a kinetic energy of 500 eV. We considered van der Waals (vdW) interactions by the DFT-D3 method.38,39 Spin polarization was included in the calculations. The 1T′-WTe2 in the bulk phase with optimized lattice constants of a = 6.313 Å and b = 3.490 Å were employed to construct the ML 1T′-WTe2. To avoid interlayer interactions, two neighboring slabs were separated by a perpendicular vacuum of 18 Å. All of the structures were fully relaxed by applying the conjugate gradient method until the total energy and Hellmann−Feynman force convergence values were less than 10−5 eV and 0.02 eV/Å, respectively. The 7 × 13 × 1 and 2 × 3 × 1 Monkhorst−Pack k-meshes were used for structural relaxation and calculations of reaction, respectively. Additionally, we used Bader charge analysis40 to calculate the amount of charge transfer. The climbing image nudged elastic band (CI-NEB) method41 was employed to study reaction pathways and energies for the surface oxidation of ML WTe2. Furthermore, ab initio molecular dynamics simulations were performed to study the dynamic behavior of O2 dissociation facilitated by a H2O molecule. The top and side views of the ML 1T′-WTe2 with a p(1 × 1) cell are shown in Figure 1a,b, respectively. W atoms are

Figure 2. Potential energy profile for O2 dissociation on ML 1T′WTe2. The potential energy barriers (Eb) are 0.66, 0.72, and 0.50 eV.

oxidation reaction, the O2 molecule initially dissociates into two individual O atoms by climbing over a reaction energy barrier of 0.66 eV. The relatively high energy barrier indicates that decomposition of the O2 molecule on the ML 1T′-WTe2 surface may take place at a slow rate. Following the above step, the two separated O atoms draw closer to the WTe2 substrate and bond to the W and Te atoms by overcoming potential barriers of 0.72 and 0.50 eV, which are marked in the potential energy profile (see Figure 2). Consequently, the oxidized ML 1T′-WTe2 reaches the most energetically favorable adsorption structure. The relatively high energy barrier to breaking the O−O covalent bond in the O2 molecule is unfavorable for the oxidative reaction of ML 1T′-WTe2, implying that the oxidative process of the 1T′-WTe2 surface may be slow if only O2 is present. Indeed, because H2O and O2 coexist in ambient conditions, H2O and O2 molecules are simultaneously adsorbed when we expose the WTe2 film in ambient conditions. As a result, the surface reaction of the ML WTe2 film in the coexistence of H2O and O2 molecules is much closer to actual conditions. Therefore, we will study the oxidation of ML 1T′-WTe2 combined with H2O. The calculated results show that the H2O molecule can easily move close to O2 by climbing over a very low energy barrier of 0.12 eV (see Figure 3a). More importantly, once H2O reaches the desirable adsorption site, O2 immediately dissociates into two isolated O atoms by

Figure 1. Top (a) and side (b) views of ML 1T′-WTe2.

sandwiched between two Te layers along the c axis, whose Te atoms are asymmetric due to structural distortion from the 1T phase to 1T′ phase. Each W atom is surrounded by six Te atoms, and the three nearest-neighbor W atoms and one Te atom form a triangular pyramidal building block. The interaction between the WTe2 and the adsorbed molecule can be quantitatively described as the adsorption energy (Ead), which is defined as Ead = Eadtot − E WTe2 − Eadmolecule

(1)

where Eadtot is the total energy of the supercell with the adsorbed molecule, EWTe2 is the energy of the pure (3 × 4) supercell WTe2, and Eadmolecule is the isolated molecule energy. 4784

DOI: 10.1021/acs.jpclett.8b01999 J. Phys. Chem. Lett. 2018, 9, 4783−4788

Letter

The Journal of Physical Chemistry Letters

To gain deeper insight into the influence of oxidation on the ML 1T′-WTe2, we analyzed the electronic properties of the oxidized WTe2 layer. The spin-polarized calculation suggests that the spin state of the adsorbed O2 molecule completely vanishes after this oxidation reaction, and thus, the oxidized WTe2 layer is nonmagnetic. The calculated band structure, DOS, and differential charge density are illustrated in Figure 4. The valence band maximum (VBM) is located along the A−X direction, while the conduction band minimum (CBM) is located at the Y point, suggesting that the low-coverage oxidized ML 1T′-WTe2 is an indirect band gap semiconductor. It is well-known that the pristine ML 1T′-WTe2 is a QSH insulator with a band gap of 55 meV.28 However, all of the crystalline symmetries are destroyed after partial oxidation. The space group of the low-coverage oxidized 1T′-WTe2 is C1; that is, there is no symmetry that can protect the topological features. Therefore, the oxidized ML 1T′-WTe2, whose band structure is no longer sensitive to spin-orbital coupling (SOC), is a trivial insulator with a band gap of 90 meV. The major contribution of the CBM and the VBM arises from the d orbitals of W, p orbitals of Te, and p orbitals of O, as shown in Figure 4a. Inspection of the DOS for Te and O further demonstrates that the formation of a Te−O bond is more likely than the formation of a W−O bond in oxidized ML 1T′WTe2, which is consistent with experimental findings that the bond energy of the O−Te bond is much stronger than those of O−W bonds on degraded WTe2.31 Charge transfer between the dissociated O atoms and the ML 1T′-WTe2 can provide a fundamental understanding of atomic bonding during oxidation. Thus, we calculated the differential charge density of the oxidized WTe2 layer. As shown in Figure 4b, because of the high electronegativity of O, the existence of dense charge acquisition around the O atoms is clearly visible, indicating that O gains electrons from the ML 1T′-WTe2 during oxidation. Meanwhile, charge depletion for the oxidized WTe2 layer is found near the O adsorption sites. Bader charge analysis shows that the O atoms capture almost one electron per O atom, and these electrons mainly originate from neighboring Te and W atoms. Surprisingly, the Te atoms tend to lose electrons more easily than the W atoms despite the high electronegativity of Te atoms. The total amount of electron transfer from the adjacent Te and W atoms is 1.51e and 0.75e, respectively.

Figure 3. (a) Diffusion behavior of H2O and O2 on ML 1T′-WTe2. (b) Depressed potential energy profile for O2 dissociation together with H2O on ML 1T′-WTe2.

overcoming a lower reaction barrier of only 0.3 eV (as shown in Figure 3b). Then, two dissociated O atoms move to Te atoms, forming Te−O bonds. Compared with the results in the absence of H2O, the presence of H2O effectively lowers the barrier of O2 dissociation and can significantly accelerate the oxidation of 1T′-WTe2. To further confirm the O2 dissociation promoted by H2O on the ML 1T′-WTe2 surface, we also perform molecular dynamics simulations to study the oxidation of ML 1T′-WTe2 together with H2O (see Figure S2 in the SI). Although the H2O molecule always retains its molecular form, the O2 molecule dissociates into oxygen atoms rapidly in view of the sharp increase in the distance between two oxygen at 8 ps. Combined with the reaction pathway obtained by using the CI-NEB approach, we can further conclude that the presence of H2O accelerates the dissociation of O2 molecules and may account for fast degradation of the WTe2 flake in air.

Figure 4. (a) Band structure and DOS of dissociated O2 on ML 1T′-WTe2. (b) Differential charge dissociation (Δρ = ρ(total) − ρ(WTe2) − ρ(2Oi)) adsorption of O2. The isosurface value is set to 0.007 e/bohr3 (violet: electron accumulation; green: electron depletion). The Fermi levels are marked by the blue dashed line. 4785

DOI: 10.1021/acs.jpclett.8b01999 J. Phys. Chem. Lett. 2018, 9, 4783−4788

Letter

The Journal of Physical Chemistry Letters

Figure 5. (a) Top and side views of ML WTe2O. (b) Phonon spectrum of ML WTe2O.

Figure 6. (a) Band structures of ML WTe2O without and with SOC. (b) Edge states of the ML WTe2O linearly cross at the X̅ point.

Dirac cone is protected by mirror-reflection symmetry. Two crossing bands belong to the opposite mirror eigenvalues ±1, and the effective model of low-energy excitations relative to the inverted center X point (kX) in the ML WTe2O can be described by

During the above discussions, we considered one only O2 molecule on a p(3 × 4) supercell to investigate the oxidation of ML 1T′-WTe2. This low-coverage oxidation significantly influences the electronic properties of ML 1T′-WTe2 and changes the topologically nontrivial ML 1T′-WTe2 into an indirect band gap semiconductor; that is, it undergoes a phase transition when oxidation occurs under low oxygen coverage. In practice, ML WTe2 is usually found in oxygen-rich conditions as it is exposed to air. Under the high oxygen coverage condition, ML 1T′-WTe2 is likely to be oxidized completely and then generate new oxidation products, for example, ML WTe2O. To investigate the fully oxidative WTe2, we predicted a 2D stoichiometric WTe2O with a rectangular lattice, which is the most energetically favorable structure. The crystal structure of the ML WTe2O is presented in Figure 5a. Clearly, this ML crystal structure is completely different from that of the pristine ML 1T′-WTe2, implying that oxidation may thoroughly change its structure. To verify dynamical stability, we performed the phonon calculation with the frozen-phonon approach, as implemented in the PHONOPY code;42 the corresponding result is shown in Figure 5b. There is no imaginary frequency in the Brillouin zone (BZ), suggesting that this structure is dynamically stable. As shown in Figure 6a, the electronic band structure without SOC show that the bands linearly cross along the X−S direction to form a Dirac cone. The ML WTe2O shares the time-reversal T and point group C2v, which contains inversion I, 2-fold rotational symmetry C2 along the z axis, and mirrorreflection symmetries mx and my. Without SOC, all the points along X−S share the unitary mirror symmetry Cs, and thus, the

H(q) = vxqx σ x +

2 vy ijjj D qy yzzz jjq − D zzσy 2 jj y qy zz k {

(2)

where σi are Pauli matrixes, q = k − kX are two components of the momentum vectors, and the two bands linearly cross at

( ,± q ), indicating the positions of the Dirac cone in the X− π a

D y

S driection. The parameters vx and vy are anisotropic Fermi velocities along x and y axes, respectively.43,44 Detailed proof of eq 2 is provided in the SI. Once SOC is considered, this Dirac cone is destroyed with a gap of 58 meV due to the strong SOC effect. Considering the inversion symmetry of the ML WTe2O, we calculated the Z2 topological invariant (ν) to reveal the band topology based on the parity eigenvalues of the occupied Bloch states at four time-reversal invariant momenta (TRIM) points. The products of parity eigenvalues at TRIM Γ, X, and Y are +1, whereas the product of parity eigenvalue at TRIM S is −1. Therefore, the Z2 topological invariant is ν = 1. This result indicates that our predicted ML WTe2O is a good 2D candidate to realize the QSH states. The other typical feature of QSH nature is the topologically protected gapless edge states that connect the valence and conduction bands. Thus, we employed a Wannier tight-binding Hamiltonian45,46 as implemented in WannierTools package47 to obtain the surface 4786

DOI: 10.1021/acs.jpclett.8b01999 J. Phys. Chem. Lett. 2018, 9, 4783−4788

Letter

The Journal of Physical Chemistry Letters

Nitride Nanomesh: Functionality from a Corrugated Monolayer. Angew. Chem., Int. Ed. 2007, 46, 5115−5119. (2) Kim, K. S.; Zhao, Y.; Jang, H.; Lee, S. Y.; Kim, J. M.; Kim, K. S.; Ahn, J.-H.; Kim, P.; Choi, J.-Y.; Hong, B. H. Large-Scale Pattern Growth of Graphene Films for Stretchable Transparent Electrodes. Nature 2009, 457, 706. (3) 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. (4) Vogt, P.; De Padova, P.; Quaresima, C.; Avila, J.; Frantzeskakis, E.; Asensio, M. C.; Resta, A.; Ealet, B.; Le Lay, G. Silicene: Compelling Experimental Evidence for Graphenelike Two-Dimensional Silicon. Phys. Rev. Lett. 2012, 108, 155501. (5) Liu, H.; Neal, A. T.; Zhu, Z.; Luo, Z.; Xu, X.; Tománek, D.; Ye, P. D. Phosphorene: An Unexplored 2D Semiconductor with a High Hole Mobility. ACS Nano 2014, 8, 4033−4041. (6) Anasori, B.; Lukatskaya, M. R.; Gogotsi, Y. 2D Metal Carbides and Nitrides (Mxenes) for Energy Storage. Nat. Rev. Mater. 2017, 2, 16098. (7) Manzeli, S.; Ovchinnikov, D.; Pasquier, D.; Yazyev, O. V.; Kis, A. 2D Transition Metal Dichalcogenides. Nat. Rev. Mater. 2017, 2, 17033. (8) 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. (9) Chhowalla, M.; Shin, H. S.; Eda, G.; Li, L.-J.; Loh, K. P.; Zhang, H. The Chemistry of Two-Dimensional Layered Transition Metal Dichalcogenide Nanosheets. Nat. Chem. 2013, 5, 263. (10) Ilatikhameneh, H.; Tan, Y.; Novakovic, B.; Klimeck, G.; Rahman, R.; Appenzeller, J. Tunnel Field-Effect Transistors in 2D Transition Metal Dichalcogenide Materials. IEEE J. Explor. Solid-State Computat. Devices Circuits 2015, 1, 12−18. (11) Yue, Q.; Shao, Z.; Chang, S.; Li, J. Adsorption of Gas Molecules on Monolayer MoS2 and Effect of Applied Electric Field. Nanoscale Res. Lett. 2013, 8, 425. (12) Naylor, C. H.; Kybert, N. J.; Schneier, C.; Xi, J.; Romero, G.; Saven, J. G.; Liu, R.; Johnson, A. T. C. Scalable Production of Molybdenum Disulfide Based Biosensors. ACS Nano 2016, 10, 6173− 6179. (13) Lee, J.; Wang, Z.; He, K.; Shan, J.; Feng, P. X. L. High Frequency MoS2 Nanomechanical Resonators. ACS Nano 2013, 7, 6086−6091. (14) Splendiani, A.; Sun, L.; Zhang, Y.; Li, T.; Kim, J.; Chim, C.-Y.; Galli, G.; Wang, F. Emerging Photoluminescence in Monolayer MoS2. Nano Lett. 2010, 10, 1271−1275. (15) Lopez-Sanchez, O.; Lembke, D.; Kayci, M.; Radenovic, A.; Kis, A. Ultrasensitive Photodetectors Based on Monolayer MoS2. Nat. Nanotechnol. 2013, 8, 497. (16) Duerloo, K.-A. N.; Reed, E. J. Structural Phase Transitions by Design in Monolayer Alloys. ACS Nano 2016, 10, 289−297. (17) Augustin, J.; Eyert, V.; Böker, T.; Frentrup, W.; Dwelk, H.; Janowitz, C.; Manzke, R. Electronic Band Structure of the Layered Compound Td-WTe2. Phys. Rev. B: Condens. Matter Mater. Phys. 2000, 62, 10812−10823. (18) Soluyanov, A. A.; Gresch, D.; Wang, Z.; Wu, Q.; Troyer, M.; Dai, X.; Bernevig, B. A. Type-II Weyl Semimetals. Nature 2015, 527, 495. (19) Sánchez-Barriga, J.; Vergniory, M. G.; Evtushinsky, D.; Aguilera, I.; Varykhalov, A.; Blügel, S.; Rader, O. Surface Fermi Arc Connectivity in the Type-II Weyl Semimetal Candidate WTe2. Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 94, 161401. (20) Wang, C.; Zhang, Y.; Huang, J.; Nie, S.; Liu, G.; Liang, A.; Zhang, Y.; Shen, B.; Liu, J.; Hu, C.; et al. Observation of Fermi Arc and Its Connection with Bulk States in the Candidate Type-II Weyl Semimetal WTe2. Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 94, 241119.

states. The edge states projected along the ky direction of the BZ are shown in Figure 6b. We can clearly observe that the metallic edge states between the valence and conduction bands linearly cross at around the X̅ point of the projected onedimensional BZ. This provides powerful evidence that the QSH state is present when ML 1T′-WTe2 is fully oxidative. Note that ML 1T′-WTe2 undergoes a topological transition from low-coverage oxidation to fully oxidative ML WTe2O. In summary, we systematically investigated the effect of surface oxidation on the topological electronic properties of ML 1T′-WTe2. The oxidation of ML 1T′-WTe2 can be promoted with the existence of H2O, which accounts for the fast surface degradation of WTe2 film. As expected, the oxidation of the ML 1T′-WTe2 under the low oxygen coverage condition is accompanied by transition from a QSH insulator to a trivial one. Strikingly, the fully oxidized 2D ML WTe2O presents a QSH insulator state with topological nontrivial edge states, suggesting the occurrence of topological phase transition from low-coverage oxidation to fully oxidative ML WTe2O. Our findings are expected to shed light on changes of the topological properties of ML 1T′-WTe2 under oxidation and improve the understanding of oxidation for WTe2.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpclett.8b01999. Dissociation structures of O2 and results of molecular dynamics at 300 K, including Figures S1 and S2, and detailed proof of the low-energy effective model (PDF)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected] (B.Z.). *E-mail: [email protected] (R.W.). *E-mail: [email protected] (H.X.). ORCID

Baobing Zheng: 0000-0003-4292-1556 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by the National Natural Science Foundation of China (NSFC, Grant Nos. 11674148, 11334004, 91634106, and 11404159), the Guangdong Natural Science Funds for Distinguished Young Scholars (No. 2017B030306008), the Basic Research Program of Science, Technology, and Innovation Commission of Shenzhen Municipality (Grant Nos. JCYJ20160531190054083 and JCYJ20170412154426330), the Fundamental Research Funds for the Central Universities of China (Nos. 106112017CDJXY300005 and cqu2018CDHB1B01), the Natural Science Basic Research plan in Shaanxi Province of China (Grant No. 2018JQ1083), and the Scientific Research Program Funded by Shaanxi Provincial Education Department (Grant No. 17JK0041).



REFERENCES

(1) Berner, S.; Corso, M.; Widmer, R.; Groening, O.; Laskowski, R.; Blaha, P.; Schwarz, K.; Goriachko, A.; Over, H.; Gsell, S.; et al. Boron 4787

DOI: 10.1021/acs.jpclett.8b01999 J. Phys. Chem. Lett. 2018, 9, 4783−4788

Letter

The Journal of Physical Chemistry Letters (21) Bruno, F. Y.; Tamai, A.; Wu, Q. S.; Cucchi, I.; Barreteau, C.; de la Torre, A.; McKeown Walker, S.; Riccò, S.; Wang, Z.; Kim, T. K.; et al. Observation of Large Topologically Trivial Fermi Arcs in the Candidate Type-II Weyl Semimetal WTe2. Phys. Rev. B: Condens. Matter Mater. Phys. 2016, 94, 121112. (22) Lv, Y.-Y.; Li, X.; Zhang, B.-B.; Deng, W. Y.; Yao, S.-H.; Chen, Y. B.; Zhou, J.; Zhang, S.-T.; Lu, M.-H.; Zhang, L.; et al. Experimental Observation of Anisotropic Adler-Bell-Jackiw Anomaly in Type-II Weyl Semimetal WTe1.98 Crystals at the Quasiclassical Regime. Phys. Rev. Lett. 2017, 118, 096603. (23) Pan, X.-C.; Chen, X.; Liu, H.; Feng, Y.; Wei, Z.; Zhou, Y.; Chi, Z.; Pi, L.; Yen, F.; Song, F.; et al. Pressure-Driven Dome-Shaped Superconductivity and Electronic Structural Evolution in Tungsten Ditelluride. Nat. Commun. 2015, 6, 7805. (24) Wang, L.; Gutiérrez-Lezama, I.; Barreteau, C.; Ubrig, N.; Giannini, E.; Morpurgo, A. F. Tuning Magnetotransport in a Compensated Semimetal at the Atomic Scale. Nat. Commun. 2015, 6, 8892. (25) Ali, M. N.; Xiong, J.; Flynn, S.; Tao, J.; Gibson, Q. D.; Schoop, L. M.; Liang, T.; Haldolaarachchige, N.; Hirschberger, M.; Ong, N. P.; et al. Large, Non-Saturating Magnetoresistance in WTe2. Nature 2014, 514, 205. (26) O’Hare, P. A. G.; Hope, G. A. Thermodynamic Properties of Tungsten Ditelluride (WTe2) II. Standard Molar Enthalpy of Formation at the Temperature 298.15 K. J. Chem. Thermodyn. 1992, 24, 639−647. (27) Zheng, F.; Cai, C.; Ge, S.; Zhang, X.; Liu, X.; Lu, H.; Zhang, Y.; Qiu, J.; et al. On the Quantum Spin Hall Gap of Monolayer 1T′WTe2. Adv. Mater. 2016, 28, 4845−4851. (28) Tang, S.; Zhang, C.; Wong, D.; Pedramrazi, Z.; Tsai, H.-Z.; Jia, C.; Moritz, B.; Claassen, M.; Ryu, H.; Kahn, S.; et al. Quantum Spin Hall State in Monolayer 1T′-WTe2. Nat. Phys. 2017, 13, 683. (29) Lee, C.-H.; Silva, E. C.; Calderin, L.; Nguyen, M. A. T.; Hollander, M. J.; Bersch, B.; Mallouk, T. E.; Robinson, J. A. Tungsten Ditelluride: A Layered Semimetal. Sci. Rep. 2015, 5, 10013. (30) Kim, Y.; Jhon, Y. I.; Park, J.; Kim, J. H.; Lee, S.; Jhon, Y. M. Anomalous Raman Scattering and Lattice Dynamics in Mono- and Few-Layer WTe2. Nanoscale 2016, 8, 2309−2316. (31) Ye, F.; Lee, J.; Hu, J.; Mao, Z.; Wei, J.; Feng, P. X. L. Environmental Instability and Degradation of Single- and Few-Layer WTe2 Nanosheets in Ambient Conditions. Small 2016, 12, 5802− 5808. (32) Naylor, C. H.; et al. Large-Area Synthesis of High-Quality Monolayer 1T′-WTe2 Flakes. 2D Mater. 2017, 4, 021008. (33) Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 47, 558−561. (34) Kresse, G.; Furthmüller, J.; Hafner, J. Ab Initio Force Constant Approach to Phonon Dispersion Relations of Diamond and Graphite. EPL (Europhysics Letters) 1995, 32, 729. (35) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (36) Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15−50. (37) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 1758−1775. (38) Dion, M.; Rydberg, H.; Schröder, E.; Langreth, D. C.; Lundqvist, B. I. Van Der Waals Density Functional for General Geometries. Phys. Rev. Lett. 2004, 92, 246401. (39) Burns, L. A.; Vazquez-Mayagoitia, Á .; Sumpter, B. G.; Sherrill, C. D. Density-Functional Approaches to Noncovalent Interactions: A Comparison of Dispersion Corrections (DFT-D), Exchange-Hole Dipole Moment (XDM) Theory, and Specialized Functionals. J. Chem. Phys. 2011, 134, 084107. (40) Henkelman, G.; Arnaldsson, A.; Jónsson, H. A Fast and Robust Algorithm for Bader Decomposition of Charge Density. Comput. Mater. Sci. 2006, 36, 354−360.

(41) Henkelman, G.; Uberuaga, B. P.; Jónsson, H. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901−9904. (42) Togo, A.; Chaput, L.; Tanaka, I. Distributions of Phonon Lifetimes in Brillouin Zones. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 91, 094306. (43) Zhu, L.; Wang, S. S.; Guan, S.; Liu, Y.; Zhang, T.; Chen, G.; Yang, S. A. Blue Phosphorene Oxide: Strain-Tunable Quantum Phase Transitions and Novel 2D Emergent Fermions. Nano Lett. 2016, 16, 6548−6554. (44) Wang, S. S.; Liu, Y.; Yu, Z. M.; Sheng, X. L.; Zhu, L.; Guan, S.; Yang, S. A. arXiv:1806.08063, 2018 (45) Mostofi, A. A.; Yates, J. R.; Lee, Y.-S.; Souza, I.; Vanderbilt, D.; Marzari, N. Wannier90: A Tool for Obtaining Maximally-Localised Wannier Functions. Comput. Phys. Commun. 2008, 178, 685−699. (46) Marzari, N.; Mostofi, A. A.; Yates, J. R.; Souza, I.; Vanderbilt, D. Maximally Localized Wannier Functions: Theory and Applications. Rev. Mod. Phys. 2012, 84, 1419−1475. (47) Wu, Q.; Zhang, S.; Song, H.-F.; Troyer, M.; Soluyanov, A. A. Wanniertools: An Open-Source Software Package for Novel Topological Materials. Comput. Phys. Commun. 2018, 224, 405−416.

4788

DOI: 10.1021/acs.jpclett.8b01999 J. Phys. Chem. Lett. 2018, 9, 4783−4788