Surface to bulk fermi arcs via weyl nodes as topological defects

ABSTRACT A hallmark of Weyl semimetal is the existence of surface Fermi arcs. An intriguing question is what determines the connectivity of surface Fermi arcs, when multiple pairs of Weyl

nodes are present. To answer this question, we show that the locations of surface Fermi arcs are predominantly determined by the condition that the Zak phase integrated along the

normal-to-surface direction is . The Zak phase can reveal the peculiar topological structure of Weyl semimetal directly in the bulk. Here, we show that the winding of the Zak phase around

each projected Weyl node manifests itself as a topological defect of the Wannier–Stark ladder, energy eigenstates under an electric field. Remarkably, this leads to bulk Fermi arcs,

open-line segments in the bulk spectra. Bulk Fermi arcs should exist in conjunction with surface counterparts to conserve the Weyl fermion number under an electric field, which is supported

by explicit numerical evidence. SIMILAR CONTENT BEING VIEWED BY OTHERS TOPOLOGICALLY DISTINCT WEYL FERMION PAIRS Article Open access 11 January 2021 EMERGENCE OF FERMI ARCS DUE TO MAGNETIC

SPLITTING IN AN ANTIFERROMAGNET Article 23 March 2022 INTERTWINED WEYL PHASES EMERGENT FROM HIGHER-ORDER TOPOLOGY AND UNCONVENTIONAL WEYL FERMIONS VIA CRYSTALLINE SYMMETRY Article Open

access 27 January 2022 INTRODUCTION Weyl semimetal is a gapless, topological phase of matter, which can be generated quite generally by breaking either time-reversal or inversion symmetry

near the phase boundary between topological and trivial insulators. One of the most dramatic properties of Weyl semimetal is that surface states have a Fermi surface consisting of open-line

segments called surface Fermi arcs1,2,3,4,5,6,7,8,9,10,11,12. Each surface Fermi arc connects two surface-projected Weyl nodes with opposite chiralities, playing an important role in

resolving the chiral anomaly of Weyl fermion13,14,15,16,17,18,19,20,21. An intriguing question is what determines the connectivity of surface Fermi arcs, when multiple pairs of Weyl nodes

are present. In this work, we answer this question by showing that the locations of surface Fermi arcs are predominantly determined by the condition that the Zak phase integrated along the

normal-to-surface direction is . The Zak phase is the Berry phase integrated along a straight, but closed path in the momentum space traversing the entire one-dimensional Brillouin zone22.

More importantly, the Zak phase can reveal the peculiar topological structure of Weyl semimetal directly in the bulk. It has been shown in a previous work23 that the non-trivial topological

order of topological insulator can be directly manifested in the winding number of the Wannier–Stark ladder (WSL), which is ultimately governed by the topological structure of the Zak phase.

The WSL is the energy eigenstates of electrons confined in the lattice under an electric field. Physically speaking, the WSL can be thought as the quantized modes of the Bloch oscillation

in a similar way that phonons are those of the lattice vibration. Here, we show that, in Weyl semimetal, the Zak phase winds by around each projected Weyl node, creating a screw dislocation

in the energy spectrum of WSL eigenstates. Eventually, these topological defects induce open-line segments in the momentum spectra of WSL, which we call bulk Fermi arcs. We provide an

argument that the existence of bulk Fermi arcs is actually required to conserve the Weyl fermion number under an electric field, which is supported by explicit numerical evidence. RESULTS

CONNECTIVITY OF SURFACE FERMI ARCS We ask if a certain bulk property of the system can determine the connectivity of surface Fermi arcs. An answer to this question would provide valuable

information to characterize Weyl semimetal without solving complicated eigenvalue equations of the microscopic Hamiltonian with an open boundary condition. To this end, let us begin by

considering graphene, which is a two-dimensional Dirac/Weyl semimetal. In graphene, edge states appear depending on the edge orientation. Delplace _et al_.24 proposed an idea that the

existence of edge states is related with the condition that the Zak phase integrated along the normal-to-edge direction is . This can be proved rigorously for certain edge orientations,

while numerically confirmed in general. Meanwhile, Mong and Shivamoggi25 provided a related, but somewhat more general proof for the existence condition of edge/surface states in

two-/three-dimensional (3D) topological insulators. Specifically, they considered the Dirac Hamiltonian, which can be written as where K⊥ and K|| (_k_||=|K|||) are the momenta perpendicular

and parallel to the normal direction to the edge/surface, respectively. Γ is a vector composed of the gamma matrices satisfying the Clifford algebra. _a_ is the lattice constant along K||.

An important assumption above is that hopping occurs only between nearest neighbours along the normal direction to the edge/surface. Under this assumption, the curve traced by H as a

function of K|| forms an ellipse, whose semi-major and semi-minor axes are 2Re[B(K⊥)] and 2Im[B(K⊥)] with its centre located at B0(K⊥). It is proved in ref. 25 that an edge/surface state

exists at K⊥ if and only if the projection of the H curve onto the Re[B(K⊥)]−Im[B(K⊥)] plane encloses the origin of H=0. Moreover, the energy of such an edge/surface state is equal to the

distance between the origin and the plane containing the H curve. This means that zero-energy edge/surface states occur when the origin is enclosed by the H curve, lying within the same

plane. For two-band models, this existence condition for zero-energy edge/surface states can be nicely rephrased in terms of the Zak phase. In two-band models, where Γ is replaced by _σ_,

there is a Dirac monopole with monopole strength _q_=±1/2 at the origin of H=0, generating the radial Berry curvature. Then, the above existence condition is precisely equivalent to the

condition that the Berry phase integrated along the H curve is , which is half the solid angle of an equator. In turn, this particular Berry phase is nothing but the Zak phase integrated

along the normal direction to the edge/surface, that is, with the Berry connection , where is the periodic part of the Bloch wave function in the -th band, which can be either valence or

conduction band in two-band models. Strictly, the applicability of the above existence condition is limited to two-band models with nearest neighbour hopping. However, this limitation can be

somewhat relaxed considering that, by its intrinsic nature, the microscopic Hamiltonian of every Weyl semimetal can be accurately approximated as a two-band, low-energy effective

Hamiltonian, which is obtained by expanding the microscopic Hamiltonian up to second order of momenta near Weyl nodes. By performing and , one can then construct a minimally

lattice-regularized two-band Hamiltonian with nearest neighbour hopping. Provided that the connectivity of surface Fermi arcs is well captured by such a minimally lattice-regularized

Hamiltonian, we predict that the locations of surface Fermi arcs (which are the zero-energy surface states) are predominantly determined by the condition that the Zak phase integrated along

the normal-to-surface direction is . This prediction is confirmed to be accurate in various theoretical models. TOPOLOGICAL DEFECTS OF THE WANNIER–STARK LADDER The Zak phase can reveal the

peculiar topological structure of Weyl semimetal directly in the bulk through the WSL emerging under an electric field. To get physical intuitions for what this means and how this is

possible, it is helpful to first understand that the WSL is the quantized modes of the Bloch oscillation, which is the semiclassical motion of electrons confined in the lattice under an

electric field. In this situation, electrons are accelerated, travelling through the momentum space until they hit the Brillouin zone boundary. Then, due to the periodic boundary condition,

they reappear at the opposite end of the Brillouin zone, eventually forming closed orbits. Such closed orbits can be quantized upon switching from the semiclassical to quantum description.

The so-obtained quantized modes are the WSL eigenstates. The quantization procedure becomes particularly simple and intuitive under the adiabatic condition that the electric field is not too

strong to cause mixing between different bands, that is, there is no Zener tunnelling. Under this condition, the Bloch oscillation originating from each band can be individually quantized

via the Bohr–Sommerfeld quantization rule. Specifically, the energy of WSL eigenstates is given as follows23: where is the one dimensionally averaged energy of the -th band, _E_ is the

electric-field strength and _n_ is the WSL index. is same as the above Zak phase except that, here, K⊥ and K|| are the momenta perpendicular and parallel to the electric field, respectively.

Physically, the Zak phase indicates the shift in the mean position of the WSL eigenstate wave packets, identified as the polarization22,26. To demonstrate concretely how the Zak phase

reveals the peculiar topological structure of Weyl semimetal, let us consider the model Hamiltonian proposed by Yang _et al_.2 which describes a time-reversal symmetry-broken Weyl semimetal:

which has two Weyl nodes at K=(±_k_0, 0, 0). From this forward, all momenta are denoted in units of 1/_a_ unless stated otherwise. Figure 1a shows the zero-energy momentum spectrum of

surface states residing in a _y_ axis-cut surface, where _k__x_ and _k__z_ are good quantum numbers. As one can see, there exists a surface Fermi arc connecting two surface-projected Weyl

nodes at (_k__x_, _k__z_)=(±_k_0, 0). It was predicted in the previous section that the locations of surface Fermi arcs are predominantly determined by the condition that the Zak phase

integrated along the normal-to-surface direction, which is the _y_ direction here, is . Figure 1b shows that this prediction is accurate. More importantly, each projected Weyl node creates a

screw dislocation in the Zak phase. See Fig. 1c for the 3D plot of the Zak phase. Such a screw dislocation in the Zak phase manifests itself as a topological defect of the WSL. Figure 1d

shows the zero-energy momentum spectrum of WSL eigenstates generated from the valence band, which is obtained via the adiabatic formula in equation (2). Specifically, in Fig. 1d, we plot the

following spectral function at : which exhibits various spectral peaks following the trajectory of . Note that becomes multi-valued if there is a screw dislocation in the Zak phase.

Consequently, the zero-energy momentum spectrum of WSL eigenstates can show, in addition to many closed loops, an open-line segment connecting two projected Weyl nodes with opposite

chiralities similar to the surface Fermi arc. We call this open-line segment the bulk Fermi arc. In fact, the conduction band generates a similar bulk Fermi arc (as well as other closed-loop

WSL eigenstates), which, incidentally, is exactly overlapped with the valence-band counterpart at zero energy in the above model Hamiltonian. Fortunately, it turns out that the bulk Fermi

arc remains robust despite mixing between WSL eigenstates generated from both valence and conduction bands. In other words, the bulk Fermi arc can persist even beyond the strictly valid

regime of the adiabatic condition, that is, _eaE_/_t_<<1. To confirm this, we compute the momentum spectrum of WSL eigenstates by directly diagonalizing the microscopic model

Hamiltonian under an electric field. Specifically, we compute the following spectral function, which is constructed in terms of the exact eigenstates of the microscopic Hamiltonian in the

presence of electrostatic potential: where is obtained by Fourier-transforming the microscopic model Hamiltonian with respect to _k__y_; , where _n__y_ is the layer index along the _y_

direction. Note that the trace Tr is taken over both _n__y_ and pseudospin index. The electrostatic potential term is given by , where the electrostatic potential is set to be zero at the

middle of the system. Figure 2 shows various zero-energy cuts of the above spectral function as a function of electric-field strength. In particular, Fig. 2a is computed at the same

electric-field strength as Fig. 1d. As one can see, the two figures are essentially identical, showing that the adiabatic formula provides an excellent approximation to the exact results, at

least at this range of electric-field strengths. Figure 2b–d show that the bulk Fermi arc persists up to reasonably strong electric fields. The model Hamiltonian in equation (3) provides a

convenient platform to study various topological properties of Weyl semimetal. The applicability of this model, however, is somewhat limited since it requires a breaking of the time-reversal

symmetry. Another pathway to generate Weyl semimetal is to break the inversion symmetry while preserving the time-reversal symmetry, which may be more relevant in view of recent

experimental confirmations of Weyl semimetal in TaAs6,7,8,9,10,11,12. In this work, we focus on the tight-binding model Hamiltonian proposed by Ojanen5, which describes a time-reversal

invariant Weyl semimetal: where _d_1(K)=_t _(1+cos K·A1+cos K·A2+cos K·A3), _d_2(K)=_t_ (sin K·A1+sin K·A2+sin K·A3) and with , , and . (Here, we reintroduce the lattice constant _a_ for

clarity.) (_σ__x_, _σ__y_, _σ__z_) and (_s__x_, _s__y_, _s__z_) are the Pauli matrices acting on the sublattice and spin basis, respectively. The other components, _D__y_(K) and _D__z_(K),

are obtained by permuting A_i_ (_i_=1, 2, 3) cyclically from the expression of _D__x_(K). The above Hamiltonian has four bands composed of two conduction and two valence bands, among which

the middle two bands, that is, the top valence and bottom conduction bands constitute a Weyl semimetal. Concretely, the Hamiltonian can be decomposed conveniently into two block-diagonalized

sublattice-basis Hamiltonians, , by first diagonalizing the spin-basis part of the Hamiltonian, : where . For _ϵ_>0 (_ϵ_<0), describes a Weyl semimetal composed of the middle two

bands, provided that , in which case there are 12 inequivalent Weyl nodes in the first Brillouin zone. These Weyl nodes are located at , , , , and , where . Meanwhile, describes a

topologically trivial insulator composed of the outer two bands. We check if such a time-reversal invariant Weyl semimetal can be also characterized by the existence of bulk Fermi arcs. To

this end, it is important to realize that the band structure of time-reversal invariant Weyl semimetal is generally more complicated than that of the time-reversal symmetry-broken

counterpart due to various band crossings. In the above Hamiltonian, it turns out that there are crossings between the top and bottom valence/conduction bands. In this situation, the WSL

eigenenergy cannot be simply given by the adiabatic formula in equation (2), but rather obtained as an eigenvalue solution of the so-called non-Abelian semiclassical Hamiltonian (NASH)23:

where is the energy dispersion and is the Berry connection with a non-Abelian structure27,28,29,30,31,32. Here, and denote the indices of all bands that cross each other. Note that the NASH

eigenvalue equation can be exactly solved by the adiabatic formula if all off-diagonal elements of the non-Abelian Berry connection are set equal to zero23. See Methods for details on how to

diagonalize the NASH efficiently to obtain the spectral function of WSL eigenstates. Figure 3a shows the zero-energy momentum spectrum of _y_ axis-cut surface states, which exhibits

multiple surface Fermi arcs. Considering that the block-diagonalized Hamiltonian of the middle two bands can be regarded as essentially a lattice-regularized Hamiltonian of Weyl semimetal

containing all Weyl nodes, it is natural to predict that the connectivity of the above surface Fermi arcs is determined by the Zak-phase condition, where the Zak phase is obtained by

integrating the Berry connection of the top valence (or bottom conduction) band along the _y_ axis. Figure 3b shows that this prediction is indeed true with excellent accuracy. More

importantly, Fig. 3c shows that each and every projected Weyl node creates a topological defect of the WSL. Specifically, see the magnified views of the boxed region (Fig. 3e–g), which

contains only two projected Weyl nodes. As one can see, there exists an edge dislocation exactly at each and every projected Weyl node. While appearing as three-way crossings at special

energy cuts, for example, in Fig. 3f, topological defects of the WSL are generically end points of an open-line segment, which is nothing but the bulk Fermi arc. Below, we provide a

heuristic explanation for the formation of these bulk Fermi arcs as well as the previous one in Fig. 1. Before doing so, it is important to mention that the above sharp structure of

topological defects gets softened in the presence of mixing between WSL eigenstates generated from all four bands including the top/bottom valence/conduction bands. Fortunately, even with

this mixing, the peculiar topological structure of Weyl semimetal is still clearly visible as a misalignment of WSL eigenstates near projected Weyl nodes. See the boxed region in Fig. 3d in

comparison with that in Fig. 3c. Figure 4 provides a heuristic explanation for the formation of bulk Fermi arcs. For simplicity, we first discuss the adiabatic situation described by

equation (2), assuming that the Zak phase plays a deciding role in determining the topology of WSL eigenstates. The Zak phase winds by either counter-clockwise or clockwise around each

projected Weyl node. This means that the WSL eigenstates with two different indices _n_ and _n_+1 can be smoothly fused together encircling a projected Weyl node. Such a fusion can cause

three-way crossings of WSL eigenstates, creating topological defects of the WSL. At general energy cuts, these three-way crossings get split in such a way that bulk Fermi arcs are formed. As

mentioned above, while softened, this structure of topological defects remains intact even in non-adiabatic situations, where mixings are allowed between WSL eigenstates generated from

different bands. WEYL FERMION NUMBER CONSERVATION UNDER AN ELECTRIC FIELD We argue that the existence of bulk Fermi arcs is actually required to conserve the Weyl fermion number under an

electric field. It was shown by Nielsen and Ninomiya13 that the chiral anomaly of Weyl fermion can be resolved by considering Weyl fermions in a crystal, or in a lattice-regularized theory.

Specifically, when parallel electric and magnetic fields are applied along the line connecting two Weyl nodes, the displacement of the Fermi surface (that is, the Weyl fermion

creation/annihilation) in one Weyl node is exactly compensated by that in the other since both Fermi surfaces are interconnected below through a one-dimensional bulk conduction channel

composed of filled states, therefore conserving the Weyl fermion number. Now, let us imagine what happens when parallel electric and magnetic fields are applied perpendicular to the

Weyl-node connecting line. In this situation, the Fermi surfaces of two Weyl nodes are not interconnected through a single one-dimensional bulk conduction channel. To conserve the Weyl

fermion number, additional conduction channels are necessary, which are provided by none other than surface Fermi arcs. Eventually, the whole conduction process forms a closed circuit

composed of two surface Fermi arcs in both sides of the surface and two one-dimensional bulk conduction channels, via which Weyl fermions can travel freely through the bulk between two

surface-projected copies of each Weyl node. Note that this conduction process has been predicted to cause an intriguing quantum oscillation in Weyl semimetal33,34,35. There is, however, a

hidden problem when this argument is applied to the situation with finite electric fields. Under any finite electric fields, the surface Fermi arc in one side is energetically far separated

from that in the opposite (provided that the system is macroscopically large). This means that the whole conduction process cannot form a closed circuit at the same energy level. A

resolution of this problem is that there exist many bulk Fermi arcs in conjunction with surface counterparts, which form a chain of many closed circuits, eventually connecting both sides of

the surface. Below, we provide explicit numerical evidence supporting this argument. Figure 5 shows layer-by-layer constant-energy momentum spectra of WSL eigenstates in the model

Hamiltonian in equation (3) under an electric field with _y_ axis-cut surfaces. Specifically, we compute the following spectral function: where _n__y_ is the layer index and the trace is

taken over only the pseudospin index. is the same as in equation (5) except that, here, a _y_ axis-cut surface is located at _n__y_=1. It is important to note that the surface Fermi arc is

joined with a partner Fermi arc at two projected Weyl nodes, forming a closed circuit together. This partner Fermi arc is the first in a series of many bulk Fermi arcs forming the periodic

structure of the WSL. One may ask how the connectivity of surface Fermi arcs evolves into that of bulk Fermi arcs. In the above example, the two connectivities happen to be the same, but in

general can be very different, as seen in Fig. 4. As explained previously, the connectivity of surface Fermi arcs is predominantly determined by the Zak-phase condition, while that of bulk

Fermi arcs is determined by a delicate interplay between the Zak phase and the band dispersion. DISCUSSION In this work, we have shown that Weyl nodes, which are responsible for the peculiar

topological structure of Weyl semimetal, can be directly visualized as topological defects of the WSL emerging under an electric field. This opens up the possibility of a novel

spectroscopic method to characterize Weyl semimetal. Below, we discuss briefly how this method can be realized in experiments. So far, the WSL has been observed only in artificial structures

such as semiconductor superlattices36,37 and optical lattices38 due to the fact that the lattice spacing in a natural crystal is usually too small that a strong electric field is necessary

to generate sufficiently well-developed WSL spectral lines; for typical experimental situations, the necessary electric-field strength is estimated to be around the order of 100 kV cm−1

(ref. 23). To overcome this obstacle, there may be two possible strategies: (i) constructing a Weyl semimetal with a large lattice spacing; or (ii) applying a strong electric field without

damaging the sample. For the first strategy, it has been proposed39,40 that a Weyl semimetal can be constructed in a superlattice system composed of alternating layers of 3D topological

insulators and ordinary insulators. Meanwhile, there has been a recent outburst of various proposals for constructing Weyl semimetals in optical lattice systems with cold

atoms41,42,43,44,45,46. Our method can be particularly useful for such cold-atom Weyl semimetals in optical lattice systems, which are known to suffer from various detection issues: (i)

edges/surfaces are not well defined47; and (ii) transport measurements are limited, or have different characteristics from those in condensed matter systems48. Our method, which detects a

bulk property in a non-transport measurement, could be an ideal alternative. For the second strategy, various pump–probe techniques can be useful since a strong electric field can be applied

in the form of pulse or radiation without damaging the sample49,50. METHODS DIAGONALIZATION OF THE NON-ABELIAN SEMICLASSICAL HAMILTONIAN Here, we discuss how to diagonalize the NASH

efficiently to obtain the spectral function of WSL eigenstates. One method is to Fourier-transform the NASH from the momentum to the real space, which involves Fourier-transforming both

energy dispersion and non-Abelian Berry connection (ref. 23). Unfortunately, this method turns out to be inefficient in Weyl semimetal due to a slow convergence of the truncation error for

higher-order Fourier components. A more efficient alternative is to rewrite the differential operator _i_∇K in a discrete momentum representation, which is convenient for numerical

diagonalization. To this end, it is important to note that _i_∇K is in fact the position operator , which is represented as a matrix in the momentum space as follows: From this forward, let

us focus on position and momentum components parallel to the electric field, which are denoted as and _k_||, respectively. Next, we note the following representation of the delta function by

using the so-called Dirichlet kernel: where . Motivated by this equality, we replace the delta function by its discrete version: where with _j_=−_N_, −_N_+1, …, _N_−1, _N_. This leads to a

matrix representation of the position operator in the discrete momentum space: This representation may seem natural in a slightly different, but more physical perspective; what we have done

is basically equivalent to Fourier-transforming the position operator from the real lattice space with a finite-length _L_=2_N_+1 to the discrete parallel momentum space with with _j_=−_N_,

−_N_+1, …, _N_−1, _N_. It is important to note that is a simple diagonal matrix with respect to both band and discrete parallel momentum indices. On the other hand, is a 2 × 2 matrix with

generally non-zero off-diagonal elements with respect to the band index, while being an _N_ × _N_ diagonal matrix with respect to the discrete parallel momentum index. Of course, is a

diagonal matrix with respect to the band index. With the knowledge of all these operators in the above discrete momentum representation, the NASH can be diagonalized to generate WSL

eigenstates as a function of perpendicular momentum, K⊥. Specifically, we compute the following spectral function of WSL eigenstates obtained from the NASH: where the trace Tr is taken over

both band and discrete parallel momentum indices. DATA AVAILABILITY All relevant data as well as computer codes are available from the authors. ADDITIONAL INFORMATION HOW TO CITE THIS

ARTICLE: Kim, K. W. _et al_. Surface to bulk Fermi arcs via Weyl nodes as topological defects. _Nat. Commun._ 7, 13489 doi: 10.1038/ncomms13489 (2016). PUBLISHER'S NOTE: Springer Nature

remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. REFERENCES * 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_ 83, 205101 (2011). Google Scholar  * Yang, K.-Y., Lu, Y.-M. & Ran, Y. Quantum

Hall effects in a Weyl semimetal: possible application in pyrochlore iridates. _Phys. Rev. B_ 84, 075129 (2011). Google Scholar  * Xu, G., Weng, H., Wang, Z., Dai, X. & Fang, Z. Chern

semimetal and the quantized anomalous hall effect in HgCr2Se4 . _Phys. Rev. Lett._ 107, 186806 (2011). Google Scholar  * Go, A., Witczak-Krempa, W., Jeon, G. S., Park, K. & Kim, Y. B.

Correlation effects on 3D topological phases: from bulk to boundary. _Phys. Rev. Lett._ 109, 066401 (2012). Google Scholar  * Ojanen, T. Helical Fermi arcs and surface states in

time-reversal invariant Weyl semimetals. _Phys. Rev. B_ 87, 245112 (2013). Google Scholar  * Xu, S.-Y. et al. Discovery of a Weyl fermion semimetal and topological Fermi arcs. _Science_ 349,

613–617 (2015). CAS  Google Scholar  * Xu, S.-Y. et al. Discovery of a Weyl fermion state with Fermi arcs in niobium arsenide. _Nat. Phys._ 11, 748–755 (2015). CAS  Google Scholar  * Lv, B.

Q. et al. Experimental discovery of Weyl Semimetal TaAs. _Phys. Rev. X_ 5, 031013 (2015). Google Scholar  * Lv, B. Q. et al. Observation of Weyl nodes in TaAs. _Nat. Phys._ 11, 724–728

(2015). Google Scholar  * Yang, L. X. et al. Weyl semimetal phase in the non-centrosymmetric compound TaAs. _Nat. Phys._ 11, 728–732 (2015). CAS  Google Scholar  * Weng, H., Fang, C., Fang,

Z., Bernevig, B. A. & Dai, X. Weyl semimetal phase in noncentrosymmetric transition-metal monophosphides. _Phys. Rev. X_ 5, 011029 (2015). Google Scholar  * Huang, S.-M. et al. A Weyl

fermion semimetal with surface Fermi arcs in the transition metal monopnictide TaAs class. _Nat. Commun._ 6, 7373 (2015). CAS  Google Scholar  * Nielsen, H. B. & Ninomiya, M. The

Adler-Bell-Jackiw anomaly and Weyl fermions in a crystal. _Phys. Lett. B_ 130, 389–396 (1983). Google Scholar  * Zyuzin, A. A. & Burkov, A. A. Topological response in Weyl semimetals and

the chiral anomaly. _Phys. Rev. B_ 86, 115133 (2012). Google Scholar  * Son, D. T. & Spivak, B. Z. Chiral anomaly and classical negative magnetoresistance of Weyl metals. _Phys. Rev. B_

88, 104412 (2013). Google Scholar  * Kim, H.-J. et al. Dirac versus Weyl fermions in topological insulators: Adler-Bell-Jackiw anomaly in transport phenomena. _Phys. Rev. Lett._ 111, 246603

(2013). Google Scholar  * Burkov, A. A. Chiral anomaly and diffusive magnetotransport in Weyl metals. _Phys. Rev. Lett._ 113, 247203 (2014). CAS  Google Scholar  * Hosur, P. & Qi, X.

Recent developments in transport phenomena in Weyl semimetals. _C. R. Phys._ 14, 857–870 (2013). CAS  Google Scholar  * Haldane, F. D. M. Attachment of surface ‘Fermi arcs’ to the bulk Fermi

surface: ‘Fermi-level plumbing’ in topological metals. Preprint at http://arxiv.org/abs/1401.0529 (2014). * Xiong, J. et al. Evidence for the chiral anomaly in the Dirac semimetal Na3Bi.

_Science_ 350, 413–416 (2015). CAS  Google Scholar  * Zhang, C.-L. et al. Signatures of the Adler-Bell-Jackiw chiral anomaly in a Weyl fermion semimetal. _Nat. Commun._ 7, 10735 (2016). CAS

  Google Scholar  * Zak, J. Berry’s phase for energy bands in solids. _Phys. Rev. Lett._ 62, 2747–2750 (1989). CAS  Google Scholar  * Lee, W.-R. & Park, K. Direct manifestation of

topological order in the winding number of the Wannier-Stark ladder. _Phys. Rev. B_ 92, 195144 (2015). Google Scholar  * Delplace, P., Ullmo, D. & Montambaux, G. Zak phase and the

existence of edge states in graphene. _Phys. Rev. B_ 84, 195452 (2011). Google Scholar  * Mong, R. S. K. & Shivamoggi, V. Edge states and the bulk-boundary correspondence in Dirac

Hamiltonians. _Phys. Rev. B_ 83, 125109 (2011). Google Scholar  * King-Smith, R. D. & Vanderbilt, D. Theory of polarization of crystalline solids. _Phys. Rev. B_ 47, 1651–1654 (1993).

CAS  Google Scholar  * Wilczek, F. & Zee, A. Appearance of gauge structure in simple dynamical systems. _Phys. Rev. Lett._ 52, 2111–2114 (1984). CAS  Google Scholar  * Shindou, R. &

Imura, K.-I. Noncommutative geometry and non-Abelian Berry phase in the wave-packet dynamics of Bloch electrons. _Nucl. Phys. B_ 720, 399–435 (2005). Google Scholar  * Culcer, D., Yao, Y.

& Niu, Q. Coherent wave-packet evolution in coupled bands. _Phys. Rev. B_ 72, 085110 (2005). Google Scholar  * Chang, M.-C. & Niu, Q. Berry curvature, orbital moment, and effective

quantum theory of electrons in electromagnetic fields. _J. Phys. Condens. Matter_ 20, 193202 (2008). Google Scholar  * Xiao, D., Chang, M.-C. & Niu, Q. Berry phase effects on electronic

properties. _Rev. Mod. Phys._ 82, 1959–2007 (2010). CAS  Google Scholar  * Kelardeh, H. K., Apalkov, V. & Stockman, M. I. Wannier-Stark states of graphene in strong electric field.

_Phys. Rev. B_ 90, 085313 (2014). Google Scholar  * Potter, A. C., Kimchi, I. & Vishwanath, A. Quantum oscillations from surface Fermi arcs in Weyl and Dirac semimetals. _Nat. Commun._

5, 5161 (2014). CAS  Google Scholar  * Moll, P. J. W. et al. Transport evidence for Fermi-arc mediated chirality transfer in the Dirac semimetal Cd3As2 . _Nature_ 535, 266–270 (2016). CAS 

Google Scholar  * Zhang, Y., Bulmash, D., Hosur, P., Potter, A. C. & Vishwanath, A. Quantum oscillations from generic surface Fermi arcs and bulk chiral modes in Weyl semimetals. _Sci.

Rep._ 6, 23741 (2016). CAS  Google Scholar  * Mendez, E. E. & Bastard, G. Wannier-Stark ladders and Bloch oscillations in superlattices. _Phys. Today_ 46, 34–42 (1993). CAS  Google

Scholar  * Wacker, A. Semiconductor superlattices: a model system for nonlinear transport. _Phys. Rep._ 357, 1–111 (2002). CAS  Google Scholar  * Raizen, M., Salomon, C. & Niu, Q. New

light on quantum transport. _Phys. Today_ 50, 30–34 (1997). CAS  Google Scholar  * Burkov, A. A. & Balents, L. Weyl Semimetal in a topological insulator multilayer. _Phys. Rev. Lett._

107, 127205 (2011). CAS  Google Scholar  * Halász, G. B. & Balents, L. Time-reversal invariant realization of the Weyl semimetal phase. _Phys. Rev. B_ 85, 035103 (2012). Google Scholar 

* Lan, Z., Goldman, N., Bermudez, A., Lu, W. & Öhberg, P. Dirac-Weyl fermions with arbitrary spin in two-dimensional optical superlattices. _Phys. Rev. B_ 84, 165155 (2011). Google

Scholar  * Anderson, B. M., Juzeliūnas, G., Galitski, V. M. & Spielman, I. B. Synthetic 3D spin-orbit coupling. _Phys. Rev. Lett._ 108, 235301 (2012). Google Scholar  * Jiang, J.-H.

Tunable topological Weyl semimetal from simple-cubic lattices with staggered fluxes. _Phys. Rev. A_ 85, 033640 (2012). Google Scholar  * Ganeshan, S. & Das Sarma, S. Constructing a Weyl

semimetal by stacking one-dimensional topological phases. _Phys. Rev. B_ 91, 125438 (2015). Google Scholar  * Dubček, T. et al. Weyl points in three-dimensional optical lattices: synthetic

magnetic monopoles in momentum space. _Phys. Rev. Lett._ 114, 225301 (2015). Google Scholar  * Zhang, D.-W., Zhu, S.-L. & Wang, Z. D. Simulating and exploring Weyl semimetal physics with

cold atoms in a two-dimensional optical lattice. _Phys. Rev. A_ 92, 013632 (2015). Google Scholar  * Bloch, I., Dalibard, J. & Nascimbène, S. Quantum simulations with ultracold quantum

gases. _Nat. Phys._ 8, 267–276 (2012). CAS  Google Scholar  * Chien, C.-C., Peotta, S. & Di Ventra, M. Quantum transport in ultracold atoms. _Nat. Phys._ 11, 998–1004 (2015). CAS  Google

Scholar  * Liu, M. et al. Terahertz-field-induced insulator-to-metal transition in vanadium dioxide metamaterial. _Nature_ 487, 345–348 (2012). CAS  Google Scholar  * Schultze, M. et al.

Controlling dielectrics with the electric field of light. _Nature_ 493, 75–78 (2013). Google Scholar  Download references ACKNOWLEDGEMENTS We are grateful to Jainendra K. Jain, Kai Sun,

Naoto Nagaosa, Takahiro Morimoto, Gil Young Cho, Jung Hoon Han, Hong Yao and Hyun-Woong Kwon for insightful discussions. This work is partially supported by the NSERC of Canada and Canadian

Institute for Advanced Research (Y.B.K.). AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * School of Physics, Korea Institute for Advanced Study, Seoul, 02455, Korea Kun Woo Kim, Woo-Ram Lee, 

Yong Baek Kim & Kwon Park * Quantum Universe Center, Korea Institute for Advanced Study, Seoul, 02455, Korea Woo-Ram Lee & Kwon Park * Department of Physics and Centre for Quantum

Materials, University of Toronto, Toronto, M5S 1A7, Ontario, Canada Yong Baek Kim Authors * Kun Woo Kim View author publications You can also search for this author inPubMed Google Scholar *

Woo-Ram Lee View author publications You can also search for this author inPubMed Google Scholar * Yong Baek Kim View author publications You can also search for this author inPubMed Google

Scholar * Kwon Park View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS K.W.K. performed all the relevant analytical as well as numerical

calculations under the supervision of K.P. The initial idea was conceived by K.W.K. based on the work by W.-R.L. and K.P. who provided details on the analysis techniques. Y.B.K. made various

contributions to the interpretation of the results. K.P. wrote the manuscript. CORRESPONDING AUTHOR Correspondence to Kwon Park. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare

no competing financial interests. SUPPLEMENTARY INFORMATION PEER REVIEW FILE (PDF 184 KB) RIGHTS AND PERMISSIONS This work is licensed under a Creative Commons Attribution 4.0 International

License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is

not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit

http://creativecommons.org/licenses/by/4.0/ Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Kim, K., Lee, WR., Kim, Y. _et al._ Surface to bulk Fermi arcs via Weyl nodes as

topological defects. _Nat Commun_ 7, 13489 (2016). https://doi.org/10.1038/ncomms13489 Download citation * Received: 01 August 2016 * Accepted: 03 October 2016 * Published: 15 November 2016

* DOI: https://doi.org/10.1038/ncomms13489 SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not

currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative