Oxidation-Induced Topological Phase Transition in Monolayer 1T

Monolayer (ML) tungsten ditelluride (WTe2) is a well-known quantum spin Hall (QSH) ... phase transition induced by oxidation provides an exotic insigh...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of South Dakota

Surfaces, Interfaces, and Catalysis; Physical Properties of Nanomaterials and Materials 2

Oxidation-Induced Topological Phase Transition in Monolayer 1T’-WTe Jiali Yang, Yuanjun Jin, Wangping Xu, Baobing Zheng, Rui Wang, and Hu Xu

J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.8b01999 • Publication Date (Web): 06 Aug 2018 Downloaded from http://pubs.acs.org on August 7, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Oxidation-induced topological phase transition in monolayer 1T’-WTe2 Jiali Yang1, 2, Yuanjun Jin2, Wangping Xu2, Baobing Zheng2, 3, 4,*, Rui Wang1, 2, *

, Hu Xu2, *

1

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

Corresponding Authors *E-mail: [email protected] *E-mail: [email protected] *E-mail: [email protected]

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 22

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 state. The topological phase transition induced by oxidation provides an exotic insight into understanding the topological features of layered dichalcogenide materials.

Table of Contents (TOC)

2

ACS Paragon Plus Environment

transition-metal

Page 3 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Since 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 crystalize 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 effect,24, 25 and 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 ML WTe2 is a 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, few-layer 1T’-WTe2 films are relatively unstable under ambient conditions because of their large exposure area.29 The surface of 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 WTe2 film from further atmospheric 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

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 oxidized ML 1T’-WTe2 during the surface reaction, including the band structure, density of states, and differential charger 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 the optimized lattice constant 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 4

ACS Paragon Plus Environment

Page 4 of 22

Page 5 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

vacuum of 18 Å. All 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 an H2O molecule. The top and side views of the ML 1T’-WTe2 with a p(1×1) cell are shown in Figs. 1(a) and 1(b), respectively. W atoms are sandwiched between two Te layers along the c axis, whose Te atoms are asymmetric due to the structural distortion from 1T phase to 1T’ phase. Each W atom is surrounded by six Te atoms, and the nearest-neighboring three 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 − EWTe2 − 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. 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 Fig. S1 of the Supporting Information (SI)) of O2 on ML 1T’-WTe2 are considered, and the corresponding adsorption energies in the dissociation 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

form are also listed. We find that the energetically favorable configuration [Fig. S1(f)], 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 the fast degradation on the 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 Fig. 2. During the 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 the decomposition of the O2 molecule on the ML 1T’-WTe2 surface may take place at 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 eV and 0.50 eV, which are marked in the potential energy profile (see Fig. 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 ML WTe2 film in the coexistence of H2O and O2 molecules is much closer to actual condition. Therefore, we will study the oxidation of ML 1T’-WTe2 combined with H2O. 6

ACS Paragon Plus Environment

Page 6 of 22

Page 7 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

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 Fig. 3(a)]. More importantly, once H2O reaches the desirable adsorption site, O2 immediately dissociates into two isolated O atoms by overcoming a lower reaction barrier of only 0.3 eV [as shown in Fig. 3(b)]. 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 Fig. S2 in 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 the fast degradation of the WTe2 flake in air. 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 non-magnetic. The calculated band structure, density of states, and differential charge density are illustrated in Fig. 4. The valence band maximum (VBM) is located along the A-X direction, while the conduction band minimum (CBM) is located at 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 the crystalline symmetries are destroyed after partially oxidation. The space group of the low-coverage 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

oxidized 1T’-WTe2 is C1; that is, there is no any symmetry which can protect the topological features. Therefore, the oxidized ML 1T’-WTe2, whose band structure is no longer sensitive to spin-orbit 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 orbits of W, p orbits of Te, and p orbits of O, as shown in Fig. 4(a). An inspection of 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 Fig. 4(b), 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, the charge depletion for the oxidized WTe2 layer is found near the O adsorption sites. The 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. During the above discussions, we only considered one 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 8

ACS Paragon Plus Environment

Page 8 of 22

Page 9 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

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 generates 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 Fig. 5(a). 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 code42; the corresponding result is shown in Fig. 5(b). There is no imaginary frequency in the Brillouin zone, suggesting that this structure is dynamically stable. As shown in Fig. 6(a), the electronic band structures 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, two-fold rotational symmetry C2 along the z axis, and mirror-reflection symmetries mx and my. Without SOC, all the points along X-S share the unitary mirror symmetry Cs, and thus the 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 inverted center X point ( k X ) in the ML WTe2O can be described by v y  D q y2  H (q) = vx q x σ x +  q y − D  σ y , q y  2 

(2)

where σ i are Pauli matrixes, q = k − k X are two components of the 

momentum vectors, and the two bands linearly cross at ( , ± ) indicating 

the positions of Dirac cone in the X-S. The parameters vx and v y are 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 22

anisotropic Fermi velocities along x and y axes, respectively.43,

44

The

detailed proof of Eq. (2) is provided in 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 eigenvalues 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 states. The edge states projected along the ky direction of the Brillouin zone (BZ) are shown in Fig. 6(b). We can clearly observe that the metallic edge states between the valence and conduction bands linearly cross around X point of the projected one-dimensional 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 topologically 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 the transition from a QSH insulator to a trivial one. Strikingly, the fully oxidized 2D ML WTe2O presents the QSH insulator 10

ACS Paragon Plus Environment

Page 11 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

state with topologically 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 the changes of the topological properties of ML 1T’-WTe2 under oxidation and improve the understanding of oxidation for WTe2.

AUTHOR INFORMATION Corresponding Authors *E-mail: [email protected] *E-mail: [email protected] *E-mail: [email protected] Notes 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,

JCYJ20170412154426330), the Fundamental Research Funds for the Central Universities of China (Nos. 106112017CDJXY300005 and 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

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).

ASSOCIATED CONTENT Supporting Information Dissociation structures of O2 and results of molecular dynamics at 300 K, including Figures S1-S2, and the detailed proof of low-energy effective model.

12

ACS Paragon Plus Environment

Page 12 of 22

Page 13 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

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 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 Rashba-Like Spin Splitting in Monolayers of the 1T’-Type Transition-Metal Dichalcogenides MX2 (M=Mo, W; X=S, Se, Te). Phy. Rev. B 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. 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

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. Nature Nanotech. 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. Phy. Rev. B 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. Phy. Rev. B 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. Phy. Rev. B 2016, 94, 241119. 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. Phy. Rev. B 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 14

ACS Paragon Plus Environment

Page 14 of 22

Page 15 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

24.

25.

26.

27.

28.

29.

30.

31.

32.

33. 34.

35.

Superconductivity and Electronic Structural Evolution in Tungsten Ditelluride. Nature Commun. 2015, 6, 7805. 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. Nature Commun. 2015, 6, 8892. 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. 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. Feipeng, Z.; Chaoyi, C.; Shaofeng, G.; Xuefeng, Z.; Xin, L.; Hong, L.; Yudao, Z.; Jun, Q.; Takashi, T.; Kenji, W.; et al. On the Quantum Spin Hall Gap of Monolayer 1T’-WTe2. Adv. Mater. 2016, 28, 4845-4851. 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. Nature Phys. 2017, 13, 683. 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. 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. 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. Carl, H. N.; William, M. P.; Zhaoli, G.; Hojin, K.; Mehmet, N.; Robert, B. W.; Liang, Z. T.; Youngkuk, K.; Christopher, E. K.; Frank, S.; et al. Large-Area Synthesis of High-Quality Monolayer 1T’-WTe2 Flakes. 2D Materials 2017, 4, 021008. Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phy. Rev. B 1993, 47, 558-561. 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. Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

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. Phy. Rev. B 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.; Mayagoitia, Á. V.-.; 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. Phy. Rev. B 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. 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. Modern 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.

16

ACS Paragon Plus Environment

Page 16 of 22

Page 17 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

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

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

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

18

ACS Paragon Plus Environment

Page 18 of 22

Page 19 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

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.

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. (a) Band structure and density of states 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.

20

ACS Paragon Plus Environment

Page 20 of 22

Page 21 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

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

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. (a) Band structure of ML WTe2O without and with SOC. (b) The edge states of ML WTe2O linearly cross at X point.

22

ACS Paragon Plus Environment

Page 22 of 22