Visualizing Type-II Weyl Points in Tungsten Ditelluride by

Oct 23, 2017 - We measured energy-dependent quasiparticle interference patterns with a cryogenic scanning tunneling microscope, revealing the position...
1 downloads 11 Views 1MB Size
Subscriber access provided by READING UNIV

Article

Visualizing Type-II Weyl Points in Tungsten Ditelluride by Quasiparticle Interference Chun-Liang Lin, Ryuichi Arafune, Ro-Ya Liu, Masato Yoshimura, Baojie Feng, Kazuaki Kawahara, Zeyuan Ni, Emi Minamitani, Satoshi Watanabe, Youguo Shi, Maki Kawai, Tai-chang Chiang, Iwao Matsuda, and Noriaki Takagi ACS Nano, Just Accepted Manuscript • DOI: 10.1021/acsnano.7b06179 • Publication Date (Web): 23 Oct 2017 Downloaded from http://pubs.acs.org on October 25, 2017

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 free 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 accessible to all readers and 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.

ACS Nano 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 28

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

ACS Nano

Visualizing Type-II Weyl Points in Tungsten Ditelluride by Quasiparticle Interference Chun-Liang Lin1, Ryuichi Arafune2, Ro-Ya Liu3#, Masato Yoshimura1, Baojie Feng3 ‡, Kazuaki Kawahara1§, Zeyuan Ni4, Emi Minamitani4, Satoshi Watanabe4, Youguo Shi5, Maki Kawai1†,TaiChang Chiang6, Iwao Matsuda3 and Noriaki Takagi1*

1

Department of Advanced Materials Science, The University of Tokyo, 5-1-5 Kashiwanoha,

Kashiwa, Chiba 277-8561, Japan 2

International Center for Materials Nanoarchitectonics (MANA), National Institute for Materials

Science, 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan 3

Institute for Solid State Physics, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba

277-8561, Japan 4

Department of Materials Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku,

Tokyo 113-8656, Japan 5

Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China

6

Department of Physics, University of Illinois, Urbana, Illinois 61801, USA

*

n-takagi@ k.u-tokyo.ac.jp

ACS Paragon Plus Environment

1

ACS Nano

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 28

ABSTRACT

Weyl semimetals (WSMs) are classified into two types, type I and II, according to the topology of the Weyl point, where the electron and hole pockets touch each other. Tungsten ditelluride (WTe2) has garnered a great deal of attention as a strong candidate to be a type-II WSM. However, the Weyl points for WTe2 are located above the Fermi level, which has prevented us from identifying the locations and the connection to the Fermi arc surface states by using angleresolved photoemission spectroscopy. Here, we present experimental proof that WTe2 is a type-II WSM. We measured energy-dependent quasiparticle interference patterns with a cryogenic scanning tunneling microscope, revealing the position of the Weyl point and its connection with the Fermi arc surface states, in agreement with prior theoretical predictions. Our results provide an answer to this crucial question and stimulate further exploration of the characteristics of WSMs.

Keywords Weyl Semimetals, topological matter, transition metal dichalcogenides, Fermi arc, quasiparticle interference, scanning tunneling microscopy, WTe2.

ACS Paragon Plus Environment

2

Page 3 of 28

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

ACS Nano

Topology in abstract mathematics has revolutionized our conventional understanding of condensed matter physics, resulting in the emergence of exotic quantum phases such as topological insulators (TIs),1,

2

topological superconductors,3,

4

Dirac semimetals,5 Weyl

semimetals (WSMs),6-10 etc. WSMs have been the subject of a deal of great attention because the quasiparticles in WSMs behave as Weyl fermions, massless chiral fermions, long sought after in high-energy particle physics. Two Weyl fermions with opposite chirality can be combined to generate a massless Dirac fermion. Such particles to realize the prediction in quantum field theory have not yet been discovered in nature but firstly observed as the quasiparticle excitation in WSMs. In WSMs, the crossing of two linear bands constitutes a Weyl point (WP). At their surfaces, unclosed surface states (SSs) emerge, called Fermi arcs. These Fermi arcs connect a pair of WPs associated with the opposite chiral charges of ±1. Similarly to the nontrivial topological SSs of the TIs, the Fermi arc SSs are protected topologically by bulk symmetry and robust from external perturbations in contrast to conventional surface states generated only by the discontinuity of the wave function at surfaces. WSMs are classified into type-I and type-II. In type-I WSMs, the WPs are located at the Fermi level (EF) to form a point-like Fermi surface (see Figure 1a). Such intriguing characteristics have been confirmed very recently for transition metal monopnictides such as TaAs, NbAs and NbP, by measuring the Fermi arc SSs through angle-resolved photoemission spectroscopy (ARPES)2-4 as well as scanning tunneling microscopy quasiparticle interference (STM-QPI).11, 12 Soon after the discovery of type-I WSMs, type-II WSMs were proposed theoretically.10 In type-II WSM, electron and hole pockets (EPs and HPs) form a tilted Weyl cone with a finite density of states at the WP, as shown in Figure 1b. In contrast to the type-I WSMs where the Weyl fermions exhibit the character of Lorentz invariance13 (i.e. the characteristic same as the invariants during Lorentz

ACS Paragon Plus Environment

3

ACS Nano

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 4 of 28

transformation, such as the speed of light and the Plank’s constant), Lorentz invariance is violated in type-II WSMs due to the symmetry breaking by its tilted-cone band structure (see Figure 1b). It generates a different type of Weyl fermions which has just recently been defined in the condense matter, type-II WSM.10 Therefore the type-II WSM provides a playground for the investigation of Lorentz-violating theories beyond the Standard Model,14, 15 giving rise to several peculiar properties such as non-saturating magnetoresistance.16 WTe2 and MoTe2 are strong candidates for type-II WSMs. Unclosed Fermi arc SSs for MoTe217, 18

and MoxW1-xTe219, 20 have been captured by ARPES and STM-QPI, confirming that these

dichalcogenides are indeed type-II WSM. In contrast, it remains an open question whether pure WTe2 is a type-II WSM or not, despite of several ARPES measurements performed so far.21-26 The ARPES studies have shown an occupied SS (denoted as SS1) consistent with the detailed theoretical analysis by the combination of density functional theory (DFT) calculations and ab initio tight binding model obtained by maximally localized Wannier function. However, the WPs calculated for WTe2 are located above the EF, which prevents both experimental identification of WP locations and determining their connectivity with the Fermi arc SSs. These crucial issues still remain to be resolved so that it is elusive whether WTe2 is a type-II WSM or not. Experimental approaches other than ARPES are therefore required. The STM-QPI technique allows us to probe both occupied and unoccupied electronic states simultaneously. Its application to WSMs has successfully elucidated detailed electronic structures.27-30 Thus STM-QPI is a suitable tool to clarify the topological properties of WTe2. Here, we identify the location of the WPs and reveal its connectivity with both occupied and unoccupied SSs through detailed analysis of energy-dependent QPI patterns.

ACS Paragon Plus Environment

4

Page 5 of 28

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

ACS Nano

RESULTS We measured spectroscopic mapping images as a function of sample voltage with lowtemperature STM for the [001] surface of an in-situ cleaved WTe2 crystal. Figure 2a shows a topographic STM image acquired soon after cleaving. As shown in the inset of Figure 2a, the diagonal lines of bright spots appear together with less-brighter zigzag stripes in between the lines. Figure 2b shows the structure model of WTe2. Three atomic layers constitute one unit layer, which stacks to crystalize into Td configuration; these atomic layers are the top-Te, centerW and bottom-Te layers in WTe2, respectively. The top-Te layer is buckled, in which the Te atoms take higher and lower vertical positions. In the high-resolution STM image, the brighter lines are derived from the higher Te atoms and the zigzag stripes are the lower Te atoms together with W atoms. Figure 2c shows the spectroscopic mapping image acquired at 40 mV. Wavy structures appear around the bright round protrusions. These protrusions are associated with defects formed mainly due to the missing of surface Te atoms. The wavy structures arise from the interference of quasiparticles scattered by the defects. We acquired the QPI patterns through the fast Fourier transformation (FFT) of the mapping image. Figure 2d shows a QPI pattern obtained from the mapping image shown in Figure 2c. We determined the ГX and ΓY directions in the QPI pattern from the spots marked by the red circles in Figure 2d. These spots originate from the surface periodic lattice shown in the inset of Figure 2d. Figure 3a shows the energy-dependent evolution of QPI patterns. Each QPI pattern contains several fine structures caused by the interferences between specified pairs of incident and

ACS Paragon Plus Environment

5

ACS Nano

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 6 of 28

scattered waves with the wave vectors of k1 and k2, respectively. To reveal the relation of these fine structures with the electronic band structure, we cut each QPI pattern along the direction specified by q = k1 - k2, and plot the cross-sectional intensity curves as functions of energy (E) and q. An example of the energy-dependent evolution of the cross-sectional intensity curves (denoted as EDE-CIC hereafter) is shown in Figure 3b, which is acquired by cutting the QPI patterns along the ГX direction, ࢗ = (‫ݍ‬௫ , 0). In the analysis of EDE-CIC, we at first consider the scattering of k1 = -k2 as a dominant channel. When the scattering does not explain the EDE-CIC, we consider the other scattering channels satisfying q = k1 - k2. Several recent researches also suggest that the scattering on the WSM can origin from the channels in which k1 ≠ -k2.27, 28 Recent spin-resolved ARPES studies WTe2 demonstrate that even though the spin functions of +k and –k mostly take antiparallel directions, they are not orthogonal to each other.22,

26

This

indicates that the restriction of tying the QPI to the spin texture is not strictly enforced. Thus, we neglect the contribution from the spin texture to the QPI. In prior theoretical calculations of the QPI patterns observed for several WSMs,27-30 similar analysis reasonably explains the QPI patterns. Taking a look at the EDE-CIC shown in Figure 3b, one can see six main features (denoted as AF). These features are observed clearly and symmetrically in both sides with respect to ‫ݍ‬௫ = 0. Each peak arises from the scattering specified by a pair of waves with (݇ଵ୶ , ݇ଵ୷ ) and (݇ଶ୶ , ݇ଶ୷ ) where ݇ଵ୷ = ݇ଶ୷ and ‫݇ = ୶ݍ‬ଵ୶ − ݇ଶ୶ .

DISCUSSIONS

ACS Paragon Plus Environment

6

Page 7 of 28

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

ACS Nano

Visualization of the WPs. In these six features, first of all we notice the feature C. As shown in the contour map in Figure 4a, the feature C appears as a high intensity spot marked by the circle at ‫୶ݍ‬۱ = 0.46 ± 0.02 Å-1 and E = 50 meV. The prior theoretical study indicates that the WP1 and WP2 are located at (kx, ky, E) = (±0.219 Å-1, ±0.038 Å-1, 52 meV) and (±0.219 Å-1, ±0.045 Å-1, 58 meV), respectively.10 The interferences relevant to the WPs with different chiralities in the opposite sides of the Brillouin zone simply should provide high intensity spots around (0.438 Å1

, 0.076 Å-1, 52 meV) and (0.438 Å-1, 0.090 Å-1, 58 meV) in the QPI patterns. The spot observed

at ‫୶ݍ‬۱ = 0.46 ± 0.02 Å-1 reasonably coincides with 0.438 Å-1, leading to the assignment that the feature C originates from the scattering between the pair of WPs. Note that we could not reveal a pair of separated spots for the WPs because two WPs are too close to each other (only 6 meV difference which is under the energy resolution in our experiment). Figure 4b and 4c shows the CICs cut along the ГY direction around 50 meV and their derivatives, respectively. The CICs reveal some shoulders within the huge peak at ࢗ = 0 and the derivatives further help us to uncover the evolution which is worthy to be noticed. Figure 4c shows the several peaks around ‫ = ୷ݍ‬0.10 ± 0.02 Å-1. The evolution of the major peaks as marked by the dash line indicates that a change happens at 50 meV. From the band structures parallel to ΓY direction as shown by Figure 4d, we notice the intervalley scattering along ΓY direction for the HP in lower energy is shorter than that for the EP in higher energy. This change from ky = 0.04 to 0.06 Å-1 well matches the shift from qy = 0.08 to 0.11 Å-1 in the QPI observation. Therefore the shift found in Figure 4c can be rationalized as the variation in intervalley scattering from HP to EP with the increase of energy, which should happen near the WPs because the WPs are created by the touch of EP and HP. This result assures the assignment of the feature C and gives the ݇୷ௐ௉ = 0.05 ±

ACS Paragon Plus Environment

7

ACS Nano

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 8 of 28

0.01 Å-1. The present QPI experiments identify the location of WPs to be (0.23 ± 0.01 Å-1, 0.05 ± 0.01 Å-1, 50 meV), providing a nice support to the calculated WPs. Connection of Fermi arc surface states to the WPs. The evolution of QPI patterns also captures the SS1 predicted theoretically and the connectivity with the WPs identified by the present QPI measurements. Figure 5a is an enlarged EDE-CIC along the ΓX direction, focusing on the feature D. As marked by green crosses in Figure 5a, the feature D starts from -70 meV and crosses the EF, connecting to the WPs identified in this study. The dispersion deduced by ‫୶ݍ‬۲ = 2݇୶۲ from the feature D nicely matches with that of SS1 determined by the previous collaborative studies of ARPES and theoretical calculations.23-25 As shown in Figure 5b to 5f, the superposition with the bulk band structures cut at several values of ݇୷ from 0 to 0.08 Å-1 demonstrates that the feature D extends in the band gap region, and that the feature D cannot be rationalized by the bulk electronic states. Besides, the red dash curve in Figure 5b indicates the Fermi arc SS obtained by previous APRES.23 It also well matches the current QPI observation. These results give direct evidence to conclude that the SS1 is the Fermi arc SS connecting to the WPs. Furthermore, we analyzed the EDE-CIC by cutting the QPI patterns along the directions making angle θ with respect to the ΓX direction as indicated by the scheme in Figure 6a. This analysis enables us to discuss the topology of the Fermi arc SS1 in the ݇୶ -݇୷ plane. The variation of feature D as a function of θ reflects the scattering processes as shown by the arrows in Figure 6b. Figure 6c shows the distribution of feature D at three different energies. The feature D distributes in a limited region of ‫ = ୷ݍ‬± 0.2 Å-1 near EF. It well reproduces the previous ARPES studies

ACS Paragon Plus Environment

8

Page 9 of 28

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

ACS Nano

show that each SS1 disperses in the range of ݇୷ = ±0.09 Å-1 in the isoenergetic surfaces near EF. Theoretically, another surface state (denoted as SS2) has been demonstrated near the electron bands above the EF.10, 23 This state has not been identified yet by using ARPES because it runs above the EF. The present QPI measurements capture the SS2. In contrast to SS1, SS2 is not very clearly seen in Figure 3b as it overlaps with the feature B. As shown by an enlarged EDE-CIC in Figure 7a, each peak of the feature B is broad, indicating that the feature B arises from the superposition of several components (marked by the green bars). Overlapping a component at smaller ‫ݍ‬௫ (green crosses in Figure 7a) by the use of ‫ = ୶ݍ‬2݇୶ with the bulk band structures provides a band branch running very close to the bulk electron bands and as the ݇௬ increases as shown by Figures 7b-f, SS2 clearly distributes in the bulk band gap. This is reasonably matched with the dispersion of SS2 predicted theoretically.10, 23 We are supposed to give some comments on both theoretical predictions from the current experiment results since Soluyanov et al. suggest that two surface states should be topologically protected10 while Bruno et al. claim that they could be trivial except the short section connecting with the WPs.23 It is difficult to judge which one is correct by the present QPI experiments since two calculations show almost the same distributions of SS1 and SS2. From the two calculations, we can only conclude that the current observation is important because the connection to WPs of two Fermi arc surface states (SS1 and SS2) are clearly revealed. The comparison of both SS1 and SS2 between the QPI measurements and the theoretical prediction23 is shown in Figure 8. As a result, because both SS1 and SS2 fit the calculation reasonably, the Fermi arc SSs observed by QPI are reliable and almost reproduce the theoretical prediction.

ACS Paragon Plus Environment

9

ACS Nano

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 28

The origins of other features. Finally, we discuss the origins of remaining four features (A, B, E and F). These features arise from the bulk electronic states. Figure 9a shows the calculated bulk band structure. We cut the bulk band structure at characteristic energies and discuss the origins specified by ‫ ۯ୶ݍ‬, ‫୶ݍ‬۰ , ‫୶ݍ‬۳ and ‫୶ݍ‬۴ based on the topologies of the isoenergetic cuts. We assign the origins to be the scatterings specified by the vectors shown in Figures 9b-d. The features A originates from the scattering between the electron pockets as shown by the arrows in Figure 9b. The feature B is the intervalley scattering between two electron pockets on the opposite sides of the Brillouin zone as shown in Figure 9c. The features E and F are caused by the scatterings between hole bands as shown in Figure 9d. Several electron and hole bands distribute around the EF, which provides various fine structures in the QPI patterns.

CONCLUSIONS The electronic structure of pure WTe2 is investigated through a detailed analysis of the energydependent evolution of QPI patterns. We identified the positions of WPs to be (kx, ky, E) = (0.23±0.01 Å-1, 0.05±0.01 Å-1, 50 meV) which are in a good agreement with the theoretical prediction10. Furthermore, we also demonstrated the connection of topological Fermi arc SSs to the WPs; SS1 disperses in the occupied state and crosses the EF to connect to the WPs, and SS2 disperses all above the EF. These features nicely agree with the theoretical prediction. The appearance of WPs and the connection with topological SSs are the two decisive factors of Weyl semimetals. The present results provide positive proof that WTe2 is a type-II WSM.

ACS Paragon Plus Environment

10

Page 11 of 28

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

ACS Nano

METHODS Sample preparation and STM measurement. We used high quality WTe2 crystals (HQ Graphene Company, The Netherlands). A piece of WTe2 sample of ~5 mm2 area and ~0.1 mm thick was stuck on a Cu plate and then mounted on the sample holder. A carbon tape was used for in situ exfoliating WTe2 sample to get a clean surface in an ultrahigh vacuum (UHV) chamber. The base pressure of the UHV chamber was kept under 2×10-10 torr. The STM images, spectroscopic images and tunneling spectra were acquired at 5 K. The modulation voltages used for the lock-in detection were set 10 mV and 4 mV of 616 Hz for measuring spectroscopic images and tunneling spectra, respectively. Calculation of Band structure. All the DFT calculations were done by using Vienna ab initio Simulation Package (VASP)31,

32

with the projected augmented wave (PAW) method.33

Exchange and correlation functional are described at the level of generalized gradient approximation (GGA) parametrized by Perdew-Burke-Ernzerhof.34 The cut off energy of the plane wave expansion used was 400 eV. The lattice constants and atom positions in the unit cell are taken from the literature.35 In the self-consistent calculations, the Brillouin zone was sampled with (24 × 12 × 8) k-points.

ACS Paragon Plus Environment

11

ACS Nano

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 12 of 28

Single-column Figure 1. Band topologies of Type-I and Type-II WSMs near the WPs. (a) (b) Band structures of Type-I and Type-II WSMs, respectively. The upper electron bands touch the lower hole bands at the WP. In Type-II WSM, Planes A and B indicate band dispersion along selected k paths on the surface, which do not cross in the case of Plane A, and cross at the WP in the case of Plane B. The band dispersions in Planes A and B are shown in the right panels, respectively. In Plane A, the electron pocket (EP) and the hole pocket (HP) are still far away from each other.

ACS Paragon Plus Environment

12

Page 13 of 28

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

ACS Nano

Single-column Figure 2. STM-QPI observations. (a) Topographic image of an in situ cleaved WTe2 surface (Vs = 40 mV, It = 2 nA, 50 ×50 nm2). The inset is a high-resolution topographic image (Vs = 50 mV, It = 0.5 nA, 5.6 ×5.6 nm2). (b) Top view and side view of the structure model of Td phase WTe2. (c) dI/dV mapping image taken at the sample voltage of Vs = 40 mV (It = 2 nA, 50 ×50 nm2) for the same area in (a). (d) FFT image (QPI pattern) of (c). The FFT image of topographic image is shown in the inset. Red circles mark the same spots observed in both FFT images of the topographic image and QPI pattern.

ACS Paragon Plus Environment

13

ACS Nano

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 14 of 28

Double-column Figure 3. Energy-dependent evolution (EDE) of QPI pattern. (a) QPI patterns from -160 meV to 160 meV. (b) The EDE of cross-sectional intensity curves (CIC) cut along the X-Г-X direction (qx = 0 corresponds to the Г point). The labels (unit: meV) on the right show the energy of curves. Some of curves are multiplied by the factors (×0.3 for 50 meV, ×0.5 for 30, 40, 60, 70, 80, 90 meV, ×0.6 for 100, 110 meV, ×0.7 for 130 meV, ×0.8 for 140, 150 meV, ×0.9 for 160, 170 meV, ×3.5 for -10 meV) to avoid overlapping with other curves. Six main features observed symmetrically on the opposite side with respect to Г point are marked from A to F.

ACS Paragon Plus Environment

14

Page 15 of 28

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

ACS Nano

Single-column Figure 4. Manifestation of WPs in QPI observation. (a) Contour map of EDE-CIC near feature C. Feature C is marked by the blue circle. (b and c) EDE-CIC along ГY direction from 30 meV to 90 meV (bottom to top) and its derivatives, respectively. The evolution of major peaks is indicated by the dash line. (d) Calculated band structures parallel to the ГY direction at kx = 0.20Å-1 (left) and 0.25Å-1 (right). The positions of the apex for HP and EP are marked by red and green arrows, respectively.

ACS Paragon Plus Environment

15

ACS Nano

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 16 of 28

Single-column Figure 5. Manifestation of SS1 in QPI observation. (a) The enlarged EDE-CIC near feature D. The green crosses mark the peaks corresponding to feature D. They start from -60 meV and increase up to 40 meV and the blue circle is the feature C. (b)-(f) Superposition of the feature D with the calculated band structure parallel to the ГX direction at ky = 0, 0.02, 0.04, 0.06, 0.08 Å-1. The red dash line in (b) is from the ARPES observation19.

ACS Paragon Plus Environment

16

Page 17 of 28

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

ACS Nano

Single-column Figure 6. Distribution of SS1 in isoenergetic ݇୶ -݇୷ planes. (a) Illustration of angle (θ) at which the EDE-CIC is cut. (b) Schematic illustration of scatterings between the Fermi arc SS1 that can be observed at different θ. (c) Distributions of the feature D at three different energy near EF.

ACS Paragon Plus Environment

17

ACS Nano

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 18 of 28

Single-column Figure 7. Manifestation of SS2 in QPI observation. (a) Enlarged EDE-CIC near observed SS2. The green crosses mark the peak corresponding to SS2 and the green bars mark the distribution of feature B. (b)-(f) Superposition of SS2 with the calculated band structures parallel to the ГX direction at ky = 0, 0.02, 0.04, 0.06, 0.08 Å-1.

ACS Paragon Plus Environment

18

Page 19 of 28

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

ACS Nano

Single-column Figure 8. Comparison of the Fermi arc SSs between our observation and the previous theoretical calculation.23 The crosses are the results from the QPI measurements (green and pink represent SS1 and SS2, respectively. The blue circle indicates the WPs) and the gray dash curve is the prediction from Ref. 23.

ACS Paragon Plus Environment

19

ACS Nano

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 20 of 28

Single-column Figure 9. Origins of scattering for other features. (a) Calculated bulk band structure along the ГX direction. (b)-(d) are the isoenergetic cuts at the energies marked by the dash lines in (a). The green arrows indicate the scattering channels associated with the features A to F besides the features C and D.

ACS Paragon Plus Environment

20

Page 21 of 28

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

ACS Nano

AUTHOR INFORMATION Corresponding Author *Email: n-takagi@ k.u-tokyo.ac.jp. Present Addresses #Current address: Institute of Physics, Academia Sinica, 128 Academia Road, Section 2, Nankang, Taipei 11529, Taiwan. ‡Current address: Hiroshima Synchrotron Radiation Center, 2-313 Kagamiyama, HigashiHiroshima, 739-0046, Japan. §Current address: Institute of Engineering Innovation, The University of Tokyo, Tokyo 1138656, Japan †Current address: Institute for Molecular Science, 38 Nishigo-naka, Myodaiji, Okazaki 4448585, Japan. Author Contributions CLL, RA, IM and NT designed the experiment. CLL carried out the experiment and analyzed the data. RA carried out the DFT calculations with the help of ZN and EM, and coded the program for QPI analysis. CLL, RA and NT wrote the manuscript with feedback from all authors. Funding Sources MEXT KAKENHI Grant Number 25110008, 26102017 and 15K17465. JSPS KAKENHI Grant Number 15J01086.

ACS Paragon Plus Environment

21

ACS Nano

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 22 of 28

U.S. Department of Energy (DOE), Office of Science (OS), Office of Basic Energy Sciences, Division of Materials Science and Engineering, under Grant No. DE-FG02-07ER46383. ACKNOWLEDGMENT RA and NT acknowledge support by MEXT KAKENHI Grant Number 25110008. RA also acknowledges support by World 25 Premier International Research Center Initiative (WPI), MEXT, Japan. KK acknowledges support by JSPS KAKENHI Grant Number 15J01086. EM acknowledges support by MEXT KAKENHI Grant Numbers 26102017 and 15K17465. TCC acknowledges support by the U.S. Department of Energy (DOE), Office of Science (OS), Office of Basic Energy Sciences, Division of Materials Science and Engineering, under Grant No. DEFG02-07ER46383.

REFERENCES (1) Hsieh, D. D.; Qian, D.; Wray, L.; Xia, Y.; Hor, Y. S.; Cava, R. J.; Hasan, M. Z. A Topological Dirac Insulator in a Quantum Spin Hall Phase. Nature 2008, 452, 970-974. (2) Zhang, H.; Liu, C.-X.; Qi, X.-L.; Dai, X.; Fang, Z.; Zhang, S.-C. Topological Insulators in Bi2Se3, Bi2Te3 and Sb2Te3 with a Single Dirac Cone on the Surface. Nat. Phys. 2009, 5, 438-442. (3) Linder, J.; Tanaka, Y.; Yokoyama, T.; Sudbø, A.; Nagaosa, N. Unconventional Superconductivity on a Topological Insulator. Phys. Rev. Lett. 2010, 104, 067001. (4) Qi, X. L. and Zhang, S. C. Topological Insulators and Superconductors. Rev. Mod. Phys. 2011, 83, 1057.

ACS Paragon Plus Environment

22

Page 23 of 28

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

ACS Nano

(5) Xu, S.-Y.; Liu, C.; Kushwaha, S. K.; Sankar, R.; Krizan, J. W.; Belopolski, I.; Neupane, M.; Bian, G.; Alidoust, N.; Chang, T.-R.; Jeng, H.-T.; Huang, C.-Y.; Tsai, W.-F.; Lin, H.; Shibayev, P. P.; Chou, F.-C.; Cava, R. J.; Hasan, M. Z. Observation of Fermi Arc Surface States in a Topological Metal. Science 2015, 347, 294−298. (6) Wan, X.; Turner, A. M.; Vishwanath, A.; Savrasov, S. Y. Topological Semimetal and Fermi-Arc Surface States in the Electronic Structure of Pyrochlore Iridates. Phys. Rev. B 2011, 83, 205101. (7) Xu, S.-Y.; Belopolski, I.; Alidoust, N.; Neupane, M.; Bian, G.; Zhang, C.; Sankar, R.; Chang, G.; Yuan, Z.; Lee, C.-C.; Huang, S.-M.; Zheng, H.; Ma, J.; Sanchez, D. S.; Wang, B.; Bansil, A.; Chou, F.; Shibayev, P. P.; Lin, H.; Jia, S.; et al. Discovery of a Weyl Fermion Semimetal and Topological Fermi Arcs. Science 2015, 349, 613−617. (8) Xu, S.-Y.; Alidoust, N.; Belopolski, I.; Yuan, Z.; Bian, G.; Chang, T.-R.; Zheng, H.; Strocov, V. N.; Sanchez, D. S.; Chang, G.; Zhang, C.; Mou, D.; Wu, Y.; Huang, L.; Lee, C.C.; Huang, S.-M.; Wang, B.; Bansil, A.; Jeng, H.-T.; Neupert, T.; et al. Discovery of a Weyl Fermion State with Fermi Arcs in Niobium Arsenide. Nat. Phys. 2015, 11, 748−754. (9) Souma, S.; Wang, Z.; Kotaka, H.; Sato, T.; Nakayama, K.; Tanaka, Y.; Kimizuka, H.; Takahashi, T.; Yamauchi, K.; Oguchi, T.; Segawa, K.; Ando, Y. Direct Observation of Nonequivalent Fermi-Arc States of Opposite Surfaces in the Noncentrosymmetric Weyl Semimetal NbP, Phys. Rev. B 2016, 93, 161112(R). (10) Soluyanov, A. A.; Gresch, D.; Wang, Z.; Wu, Q. S.; Troyer, M.; Dai, X; Bernevig, B. A. Type-II Weyl Semimetals, Nature 2015, 527, 495-498.

ACS Paragon Plus Environment

23

ACS Nano

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 24 of 28

(11) Chang, G.; Xu, S.-Y.; Zheng, H.; Lee, C.-C.; Huang, S.-M.; Belopolski, I.; Sanchez, D. S.; Bian, G.; Alidoust, N.; Chang, T.-R.; Hsu, C.-H.; Jeng, H.-T.; Bansil, A.; Lin, H.; Hasan, M. Z. Signatures of Fermi Arcs in the Quasiparticle Interferences of the Weyl Semimetals TaAs and NbP. Phys. Rev. Lett. 2016, 116, 066601. (12) Batabyal, R.; Morali, N.; Avraham, N.; Sun, Y.; Schmidt, M.; Felser, C.; Stern, A.; Yan, B,; Beidenkopf, H. Visualizing Weakly Bound Surface Fermi Arcs and Their Correspondence to Bulk Weyl Fermions. Science Adv. 2016, 2, e1600709. (13) Pavlopoulos, T. G. Breakdown of Lorentz Invariance. Phys. Rev. 1967, 159, 1106-1109. (14) Colladay, D. and Kostelecký, V. A. Lorentz-Violating Extension of the Standard Model. Phys. Rev. D 1998, 58, 116002. (15) López-Fogliani, D. E. and Muñoz, C. Proposal for a Supersymmetric Standard Model. Phys. Rev. Lett. 2006, 97, 041801. (16) Ali, M. N.; Xiong, J.; Flynn, S.; Tao, J.; Gibson, Q. D.; Schoop, L. M.; Liang, T.; Haldolaarachchige, N.; Hirschberger, M.; Ong, N. P.; Cava R. J. Large, Non-Saturating Magnetoresistance in WTe2. Nature 2014, 514, 205-208. (17) Huang, L.; McCormick, T. M.; Ochi, M.; Zhao, Z.; Suzuki, M.-T.; Arita, R.; Wu, Y.; Mou, D.; Cao, H.; Yan, J.; Trivedi, N.; Kaminski, A. Spectroscopic Evidence for a Type II Weyl Semimetallic State in MoTe2. Nat. Mater. 2016, 15, 1155–1160. (18) Deng, K.; Wan, G.; Deng, P.; Zhang, K.; Ding, S.; Wang, E.; Yan, M.; Huang, H.; Zhang, H.; Xu, Z.; Denlinger, J.; Fedorov, A.; Yang, H.; Duan, W.; Yao, H.; Wu, Y.; Fan, S.;

ACS Paragon Plus Environment

24

Page 25 of 28

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

ACS Nano

Zhang, H.; Chen, X.; Zhou, S. Experimental Observation of Topological Fermi Arcs in Type-II Weyl Semimetal MoTe2. Nat. Phys. 2016, 12, 1105–1110. (19) Chang, T.-R.; Xu, S. Y.; Chang. G.; Lee, C.-C.; Huang, S.-M.; Wang, B. K.; Bian, G.; Zheng, H.; Sanchez, D. S.; Belopolski, I.; Alidoust, N.; Neupane, M.; Bansil, A.; Jeng, H.T.; Lin, H.; Hasan M. Z. Prediction of an Arc-Tunable Weyl Fermion Metallic State in MoxW1-xTe2. Nat. Commun. 2016, 7, 10639. (20) Belopolski, I.; Sanchez, D. S.; Ishida, Y.; Pan, X; Yu, P.; Xu, S.-Y.; Chang, G.; Chang, T.R.; Zheng, H.; Alidoust, N.; Bian, G.; Neupane, M.; Huang, S.-M.; Lee, C.-C.; Song, Y. Bu, H.; Wang, G.; Li, S.; Eda, G.; Jeng, H.-T.; et al. Discovery of a New Type of Topological Weyl Fermion Semimetal State in MoxW1−xTe2. Nat. Commun. 2016, 7, 13643. (21) Jiang, J.; Tang, F.; Pan, X.; Liu, H.; Niu, X.; Wang, Y.; Xu, D.; Yang, H.; Xie, B.; Song, F. Signature of Strong Spin-Orbital Coupling in the Large Nonsaturating Magnetoresistance Material WTe2. Phys. Rev. Lett. 2015, 115, 166601. (22) Das, P. K.; Di Sante, D.; Vobornik, I.; Fujii, J.; Okuda, T.; Bruyer, E.; Gyenis, A.; Feldman, B. E.; Tao, J.; Ciancio, R.; Rossi, G.; Ali, M. N.; Picozzi, S.; Yadzani, A.; Panaccione, G.; Cava, R. J. Layer-Dependent Quantum Cooperation of Electron and Hole States in the Anomalous Semimetal WTe2. Nat. Commun. 2016, 7, 10847. (23) Bruno, F. Y.; Tamai, A.; Wu, Q. S.; Cucchi, I.; Barreteau, C.; dela Torre, A.; McKeown Walker, S.; Riccò, S.; Wang, Z.; Kim, T. K.; Hoesch, M.; Shi, M.; Plumb, N. C.; Giannini, E.; Soluyanov, A. A.; Baumberger, F. Observation of Large Topologically Trivial Fermi Arcs in the Candidate Type-II Weyl Semimetal WTe2. Phys. Rev. B 2016, 94, 121112(R).

ACS Paragon Plus Environment

25

ACS Nano

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 26 of 28

(24) Wu, Y.; Mou, D.; Jo, N. H.; Sun, K.; Huang, L.; Bud’ko, S. L.; Canfield, P. C.; Kaminski, A. Observation of Fermi Arcs in the Type-II Weyl Semimetal Candidate WTe2. Phys. Rev. B 2016, 94, 121113(R). (25) 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 2016, 94, 161401(R). (26) Feng, B.; Chan, Y.-H.; Feng, Y.; Liu, R.-Y.; Chou, M.-Y.; Kuroda, K.; Yaji, K.; Harasawa, A.; Moras, P.; Barinov, A.; Malaeb, W.; Bareille, C.; Kondo, T.; Shin, S.; Komori, F.; Chiang, T.-C.; Shi Y.; Matsuda I. Spin Texture in Type-II Weyl Semimetal WTe2. Phys. Rev. B 2016, 94, 195134. (27) Zheng, H.; Xu, S.-Y.; Bian, G.; Guo, C.; Chang, G.; Sanchez, D. S.; Belopolski, I.; Lee, C.C.; Huang, S.-M.; Zhang, X.; Sankar, R.; Alidoust, N.; Chang, T.-R.; Wu, F.; Neupert, T.; Chou, F.; Jeng, H.-T.; Yao, N.; Bansil, A.; Jia, S.; et al. Atomic-Scale Visualization of Quantum Interference on a Weyl Semimetal Surface by Scanning Tunneling Microscopy. ACS Nano 2016, 10, 1378-1385 (2016). (28) Kourtis, S.; Li, J.; Wang, Z.; Yazdani, A.; Bernevig B. A. Universal Signatures of Fermi Arcs in Quasiparticle Interference on the Surface of Weyl Semimetals. Phys. Rev. B 2016, 93, 041109(R). (29) Mitchell, A. K. and Fritz, L. Signatures of Weyl Semimetals in Quasiparticle Interference. Phys. Rev. B 2016, 93, 035137.

ACS Paragon Plus Environment

26

Page 27 of 28

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

ACS Nano

(30) Zheng, H.; Bian, G.; Chang, G.; Lu, H.; Xu, S.-Y.; Wang, G.; Chang, T.-R.; Zhang, S.; Belopolski, I.; Alidoust, N.; Sanchez, D. S.; Song, F.; Jeng, H.-T.; Yao, N.; Bansil, A.; Jia, S.; Lin, H.; Hasan, M. Z. Atomic-Scale Visualization of Quasiparticle Interference on a Type-II Weyl Semimetal Surface. Phys. Rev. Lett. 2016, 117, 266804 (2016). (31) Kresse, G. and Furthmuller, J. Efficient Iterative Schemes for ab initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 56, 11169. (32) Kresse, G. and Furthmuller, J. Efficiency of ab-initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput Mater. Sci. 1996, 6, 16. (33) Bloch, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 24, 17953. (34) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865. (35) Mar, A.; Jobic, S.; Ibers, J. A. Metal-metal vs Tellurium-Tellurium Bonding in WTe2 and Its Ternary Variants TaIrTe4 and NbIrTe4. J. Am. Chem. Soc. 1992, 114, 8963-8971.

ACS Paragon Plus Environment

27

ACS Nano

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 28 of 28

Table of Contents

ACS Paragon Plus Environment

28