# Multi-terminal Conductance at the Surface of a Weyl Semimetal

###### Abstract

Weyl semimetals are a new paradigmatic topological phase of matter featuring a gapless spectrum. One of its most distinctive features is the presence of Fermi arc surface states. Here, we report on atomistic simulations of the dc conductance and quantum Hall response of a minimal Weyl semimetal. By using scattering theory we show that a quantized Hall conductance with a non-vanishing longitudinal conductance emerges associated to the Fermi arc surface states with a remarkable robustness to high concentrations of defects in the system. Additionally, we predict that a slab of a Weyl semimetal with broken time-reversal symmetry bears persistent currents fully determined by the system size and the lattice parameters.

## I Introduction

Topological but gapless, Weyl semimetals Hosur and Qi (2013); Vafek and Vishwanath (2014) challenge our physical intuition with striking properties, thereby constituting a new paradigm in the physics of topological quantum matter. One of their distinctive features is the presence of Weyl nodes, discrete points in the three-dimensional (3D) Brillouin zone where the valence and conduction bands meet forming a 3D analog of graphene. Unlike the Dirac points in graphene, Weyl nodes are non-degenerate and a perturbation does not gap the system, it only shifts the Weyl nodes in k-space Yan and Felser (2017). Interestingly, the existence of Weyl nodes in the sample’s bulk dispersion, which appears in pairs of opposite chirality acting as Berry phase monopoles, can be tied to the existence of Fermi arc surface states at the boundary of a finite system. Control over the exotic electronic properties of a Weyl semimetal is expected to lead to the development of novel devices. In particular, high temperature quantum devices are envisioned to take advantage of the topological protection of their states. Superconducting Weyl semimetals, expected to host novel forms of Majorana fermions, are likely to aid in the development of quantum information devices. Jia et al. (2016); Parameswaran et al. (2014); Chan et al. (2016); Wang et al. (2013); Chan et al. (2017)

Fueled by the experimental observation of Weyl semimetals in recent years Xu et al. (2015); Lu et al. (2015); Lv et al. (2015), a plethora of different flavors is now known to exist Bradlyn et al. (2016) and there is a fierce competition to unveil them in the laboratory. While most experiments rely heavily on probing the electronic dispersion through angle resolved photoemission spectroscopy (ARPES), transport experiments are also expected to provide for valuable information Reis et al. (2016); Belopolski et al. (2016). In the bulk, the conductivity tensor is expected to be anisotropic and, depending on the crystal symmetry, a giant anomalous Hall effect may arise Yang et al. (2011) (, with being the momentum separation between the Weyl nodes).

In spite of the intense research, with few exceptions Baireuther et al. (2016); Igarashi and Koshino (2017); Wang et al. (2017), most studies have focused on the bulk rather than the surface transport properties of Weyl semimetals Yang et al. (2011); Hosur et al. (2012); Nandy et al. (2017). Besides being an easier to access property than its bulk counterpart, the transport response measured at the surface may guard new hallmarks of the Weyl semimetal phase. Theoretical simulations may help to unveil them while clarifying many interesting aspects that are unique to this topological phase. For example, in contrast to topological insulators, the bulk states in Weyl Semimetals are expected to play a role in transport Gorbar et al. (2016). But the intricate electronic structure of many Weyl semimetals featuring sometimes dozens of Weyl nodes poses an obstacle in this young field. Studying a simple ”hydrogen-like” model of a Weyl semimetal is, therefore, a first necessary step.

Here, we report on the electronic and transport properties of a minimal Weyl semimetal with only two Weyl nodes. This minimal model, which requires broken time-reversal symmetry, is used as a playground for the study of the multiterminal conductance at the surface of a finite slab within a coherent approach. To such end, we combine a scattering formulation of transport together with a tight-binding model. Our simulations show that associated with the Fermi arc surface states, there is a nearly quantized Hall conductance which remains robust to imperfections and defects. Furthermore, and unlike to Chern insulators, this Hall response is accompanied by a non-vanishing longitudinal conductance. These two observations, the deviations of the Hall conductance from perfect quantization and the non-vanishing longitudinal conductance, are a consequence of the unavoidable leak of the surface states into the bulk. We also test the topological nature of our quantized conductance results by adding random lattice defects and find that they remain robust at concentrations of the order of 25% which is remarkably higher than what has been previously reported on 2D topological systems Castro et al. (2015).

Interestingly, our results also show that a slab of a Weyl semimetal can bear persistent currents even in the absence of a magnetic field. These currents are governed by the system size and the intrinsic time-reversal symmetry breaking lattice parameters and can explain why the surface transport is impeded in the direction perpendicular to the Fermi arcs.

## Ii Hamiltonian model

For our calculations we use a model of spinful electrons on a cubic lattice (one orbital per site with lattice constant ) Yang et al. (2011); Vafek and Vishwanath (2014),

(1) | ||||

As noticed in Ref. Yang et al., 2011, this model has two, four, six or eight Weyl nodes depending on the interplay between the hoppings , , and the mass term determined by . Each pair of nodes is related by inversion symmetry which is preserved as there is a single orbital per site and only one site in the lattice basis.

The breaking of the time reversal symmetry is governed by the parameter in the Hamiltonian. Let us analyze this in more detail. The time-reversal operator can be written as , where is a unitary operator and is the complex conjugation operator ( (see for example Haake (2001), chapter 2), one can also show that , and ). For spin particles, the unitary operator must have the form: , with an arbitrary phase. This satisfies all the required conditions for a time-reversal operator: it is anti-unitary, changes the sign of spin and is a matrix. The second condition can be expressed as: .

## Iii Electronic Structure

For the electronic structure and transport calculations that follow we set and so that our model contains only two Weyl nodes located at as shown in the bulk dispersion calculated along high symmetry points in Fig. 1(a). This phase is characterized by a Chern number in the 2D Brillouin planes located between the Weyl nodes, and a Chern number in the planes located outside of the nodes (see Fig. 1(b)). The topological nature of this phase motivates us to ask about the transport response originated by the surface states in this minimal case, but before turning to this question, we need to review its electronic structure in more detail.

The numerical calculations were implemented in Python using the modules: Kwant Groth et al. (2014) and Python tight-binding
^{1}^{1}1The Python Tight-Binding package is available at this link.

As mentioned earlier, the Fermi arcs are gapless bands that connect the projections of the Weyl nodes in a given plane of the Brillouin zone. Since we are in a phase with a single pair of nodes, only one Fermi arc is to be expected on this surface. To evidence the Fermi arcs we need to terminate the system at the correct surface. Figure 2(a) and (b) show the band structure for a system with lattice sites in the z-direction (and infinite on the other two). In Fig. 2(a), the two gapless bands at zero energy form the Fermi arc. Nevertheless, these dispersionless bands associated to the Fermi arc in the direction do have a non-vanishing group velocity along as shown in Fig. 2(b).

Next, we focus on the electronic structure of a Weyl semimetal wire with a square cross-section of lattice sites in the , and planes respectively. The band structures for the wires are shown in Fig. 2(c) and (d). Figure 2(c) corresponds to a wire extending along (i.e. with the cross-section in the plane), the results for a wire extending along are identical. The projection of the two Weyl Nodes is a single Dirac cone that is doubly degenerate and hosts edge states around its center. The wire with the cross-section in the yz plane has a different band structure as shown in Fig. 2(d). Flat bands appear around that do not carry associated edge states. These flat bands are the manifestation of the Fermi arcs: due to the to finite size effects of the wire geometry they present a splitting which scales with .

## Iv Multiterminal transport response

Let us now turn to the transport response of this ‘hydrogen-model’ of a Weyl semimetal. We are interested, in particular, in a Hall setup at one of the surfaces of a three-dimensional sample. Specifically, we consider six two-dimensional electrodes connected (in an H-configuration) to the surface of a rectangular cuboid as represented in Fig. 3(g). Note that this differs from other works where the sample is accessed with the electrodes on the full volume Takane (2016); Kobayashi et al. (2015) and not just the surface. In the following we present simulations for a rectangular cuboid of lattice sites in the , and -directions, respectively. Each lead is described by a two-dimensional square lattice with lattice constant matching that of the 3D cuboid and lattice sites in width. The separation between the voltage probes is set to lattice sites. The 2D lattice has nearest neighbors and spin independent hoppings equal to . This value also characterizes the hopping from sites in the leads to the matching points in our structure.

The setup studied here (see scheme in Fig. 3(g)) is typically used in Hall measurements and allow for the determination of the longitudinal and Hall resistances associated with the edge states. Furthermore, a planar architecture for the contacts may offer advantages, such as easier gating and tuning, over three-dimensional counterparts. Such contacts problems together with difficulty in tuning the Weyl semimetals themselves Belopolski et al. (2016). In our case, we would like to determine the fingerprints of the topology of this 3D system in the transport response measured at a surface. The answer is a priori not obvious since the states associated to the Fermi arcs are surface states, and not edge states as in a Chern insulator. Interestingly, one could also study the crossover between those two cases by varying the dimension of the cuboid so that it becomes 2D (we will come back to this when discussing Fig. 4), a transition that has elicited much interest Burkov and Balents (2011); Yokomizo and Murakami (2017).

Our numerical simulations are restricted to the coherent regime, where scattering theory is valid. Specifically, we apply the Landauer-Büttiker formalism Datta (1995) to this three-dimensional system (see schemes in Fig. 3(g)) with the aid of the Kwant module Groth et al. (2014). Fig. 3 shows the two-terminal (panel a), longitudinal (panel b) and Hall conductances (panel c) as a function of the Fermi energy in the zero-temperature limit. These results correspond to a setup where the leads are attached to the surface. We emphasize that our results (as well as those of Ref. Takane (2016)) correspond to conductances rather than conductivities as calculated by other authors Chen et al. (2015); Shapourian and Hughes (2016). In the latter case, since they are bulk quantities, the effect of Fermi arcs are less evident.

In the first place, as shown in Fig. 3(c), we find a nearly-quantized Hall conductance . This value of the Hall conductance is connected with the Chern number on the plane in the Brillouin zone that is outside the Weyl nodes because the planes with non-zero Chern number can be regarded as the stacking of 2D Chern insulators in a quantum Hall effect phase Yang et al. (2011). The differences between the simulated values and perfect quantization are in the order of , this is , except at a few discrete points that coincide with the onset of new subbands either in the semimetal or in the leads.

The non-vanishing value of the Hall conductance even in the absence of an external magnetic field is a consequence of the intrinsic breaking of time reversal symmetry in our model Hamiltonian, i.e. by the role played by the parameter , much as in the (anomalous) quantum Hall effect Haldane (1988); Nagaosa et al. (2010). But besides this similarity, there are also differences that we point out below. They concern the longitudinal conductivity and the robustness to disorder.

In addition, we observe a non-vanishing longitudinal conductance as shown in Fig. 3(b). This is connected with the observed small deviations from perfect quantization of the Hall conductance and can be interpreted in terms of the finite and non-vanishing penetration length of the surface states, which decay exponentially into the bulk.

To get further insight, we test the robustness of these results to defects. To such end, we add vacancies at random within to the sample and computed how the two-terminal, longitudinal and Hall conductances changed in the same setup. The results are shown in the panels (d), (e) and (f) of Fig. 3, averaged realizations of disorder. For a given defects concentration one chooses a set of sites and transform them into vacancies. The results turn out to be remarkably robust to disorder, with the conductance being substantially reduced only at about 25% of defects. This is much higher than expected, for example, in the two-dimensional Haldane model Castro et al. (2015).

Indeed, an interesting question is how does the transport response evolve as the cuboid becomes thinner and, eventually, two-dimensional. In the present setup this can be done without changing the measurement probes which are on the surface. Fig. 4 shows the results for the two-terminal conductance , the longitudinal conductance , and the hall conductance (computed at ) as a function of the dimension of the cuboid along , . spans from the two-dimensional limit () to as used in Fig. 3.

Figure 4 shows that in the two-dimensional limit transport occurs as in a Chern insulator, with the Hall conductance being perfectly quantized and a vanishing longitudinal resistance. This picture changes as the sample becomes thicker until it converges to a situation where the Hall conductance aparts itself slightly from the perfect quantization with the concomitant non-vanishing longitudinal resistance. Besides clarifying the effect of the sample size on our previous results, this figure highlights the role of the dimensionality and the effect of the penetration of the surface states into the bulk.

To complete our analysis we study the transport response on the plane (the same results are obtained for the plane). We find that there is no Hall conductance response and that as shown in Fig. 5(a). These results are a consequence of the dispersion shown in Fig. 2(a) and (b). The Fermi arcs have group velocity in the and directions only, and thus do not allow for transport in the direction. The previous results motivate us to analyze how do the currents flow in a finite sample of the semimetal. Since transport is impeded in the direction, there also should not be current density flowing along .

## V Chiral States and Persistent Currents

Inspired by the previous conclusion, we calculate the bond currents Todorov (2002) at for a finite sample of the semimetal. We find that even in the absence of an external magnetic field, the system develops persistent currents around the axis. These currents flow through the boundary of the semimetal as shown schematically in Fig. 6. We attribute these results to the time-reversal breaking acting as an intrinsic magnetic field pointing parallel to the Fermi arcs. The currents in the system are only of the solenoid type; transport along the axis is impeded in the pristine system (the states at the boundary are dispersionless), thereby producing a very high longitudinal resistance in the direction compared to the other directions. This picture helps to explain our transport results.

The above picture is still not quantative. Figure 7 (a) shows the magnitude of the maximum surface current (maximum bond current) as a function of the time reversal breaking parameter (Fig. 7 (c) shows the Chern numbers for two different planes in the Brillouin zone as a function of ). Different curves correspond to different system sizes of a cube. The horizontal axis ranged at , the value which signals the beginning of the topological phase. Figure 7 (b) shows a detail of the current dependence with the system size at fixed .

## Vi Final Remarks

In summary, we have explored the multiterminal transport properties of a minimal model of a Weyl semimetal with broken time-reversal symmetry. This kind of hydrogen-atom model of a Weyl semimetal allows for a simplified analysis. Our results show a nearly quantized Hall conductance, a result that holds up to remarkably high concentrations of disorder. Interestingly, in contrast with the quantum Hall effect studied in time-reversal invariant Weyl semimetals, the Hall response in our time-reversal broken model is not associated to formation of edge states (at the surface of the Weyl semimetal) but rather to the surface states.

Being impervious to disorder, the transport properties of the surfaces of a Weyl semimetal have a great potential in applications in nanotechnological devices. In this way, our calculations provide the motivation and are the starting point of new simulations required to elucidate the role of finite geometries upon the transport properties of more complex Weyl semimetals.

Acknowledgments.– We acknowledge financial support from Universidad de Chile through Program ’Inserción Académica 2016. ASN acknowledges funding from grants Fondecyt Regular 1150072 and from Financiamiento Basal para Centros Científicos y Tecnológicos de Excelencia, under Project No. FB 0807, CEDENNA (Chile). LEFFT acknowledges the support of Fondecyt Regular 1170917.

## References

- Hosur and Qi (2013) P. Hosur and X. Qi, “Recent developments in transport phenomena in Weyl semimetals,” Comptes Rendus Physique Topological insulators / Isolants topologiquesTopological insulators / Isolants topologiques, 14, 857 (2013).
- Vafek and Vishwanath (2014) O. Vafek and A. Vishwanath, “Dirac Fermions in Solids: From High-Tc Cuprates and Graphene to Topological Insulators and Weyl Semimetals,” Annual Review of Condensed Matter Physics 5, 83 (2014).
- Yan and Felser (2017) B. Yan and C. Felser, “Topological materials: Weyl semimetals,” Annual Review of Condensed Matter Physics 8, 337 (2017).
- Jia et al. (2016) S. Jia, S.-Y. Xu, and M. Z. Hasan, “Weyl semimetals, fermi arcs and chiral anomalies,” Nat Mater 15, 1140 (2016).
- Parameswaran et al. (2014) S. A. Parameswaran, T. Grover, D. A. Abanin, D. A. Pesin, and A. Vishwanath, “Probing the chiral anomaly with nonlocal transport in three-dimensional topological semimetals,” Phys. Rev. X 4, 031035 (2014).
- Chan et al. (2016) C.-K. Chan, P. A. Lee, K. S. Burch, J. H. Han, and Y. Ran, “When chiral photons meet chiral fermions: Photoinduced anomalous hall effects in weyl semimetals,” Phys. Rev. Lett. 116, 026805 (2016).
- Wang et al. (2013) Y. H. Wang, H. Steinberg, P. Jarillo-Herrero, and N. Gedik, “Observation of floquet-bloch states on the surface of a topological insulator,” Science 342, 453 (2013).
- Chan et al. (2017) C.-K. Chan, N. H. Lindner, G. Refael, and P. A. Lee, “Photocurrents in weyl semimetals,” Phys. Rev. B 95, 041104 (2017).
- Xu et al. (2015) S.-Y. Xu, I. Belopolski, N. Alidoust, M. Neupane, G. Bian, C. Zhang, R. Sankar, G. Chang, Z. Yuan, C.-C. Lee, S.-M. Huang, H. Zheng, J. Ma, D. S. Sanchez, B. Wang, A. Bansil, F. Chou, P. P. Shibayev, H. Lin, S. Jia, and M. Z. Hasan, “Discovery of a Weyl Fermion semimetal and topological Fermi arcs,” Science , aaa9297 (2015).
- Lu et al. (2015) L. Lu, Z. Wang, D. Ye, L. Ran, L. Fu, J. D. Joannopoulos, and M. Soljačić, “Experimental observation of Weyl points,” Science 349, 622 (2015).
- Lv et al. (2015) B. Lv, H. Weng, B. Fu, X. Wang, H. Miao, J. Ma, P. Richard, X. Huang, L. Zhao, G. Chen, Z. Fang, X. Dai, T. Qian, and H. Ding, “Experimental Discovery of Weyl Semimetal TaAs,” Physical Review X 5, 031013 (2015).
- Bradlyn et al. (2016) B. Bradlyn, J. Cano, Z. Wang, M. G. Vergniory, C. Felser, R. J. Cava, and B. A. Bernevig, “Beyond Dirac and Weyl fermions: Unconventional quasiparticles in conventional crystals,” Science 353, aaf5037 (2016).
- Reis et al. (2016) R. D. d. Reis, M. O. Ajeesh, N. Kumar, F. Arnold, C. Shekhar, M. Naumann, M. Schmidt, M. Nicklas, and E. Hassinger, “On the search for the chiral anomaly in Weyl semimetals: the negative longitudinal magnetoresistance,” New Journal of Physics 18, 085006 (2016).
- Belopolski et al. (2016) I. Belopolski, D. S. Sanchez, Y. Ishida, X. Pan, P. Yu, S.-Y. Xu, G. Chang, T.-R. Chang, H. Zheng, N. Alidoust, G. Bian, M. Neupane, S.-M. Huang, C.-C. Lee, Y. Song, H. Bu, G. Wang, S. Li, G. Eda, H.-T. Jeng, T. Kondo, H. Lin, Z. Liu, F. Song, S. Shin, and M. Z. Hasan, “Discovery of a new type of topological Weyl fermion semimetal state in MoxW1mxte2,” Nature Communications 7, 13643 (2016).
- Yang et al. (2011) K.-Y. Yang, Y.-M. Lu, and Y. Ran, “Quantum Hall effects in a Weyl semimetal: Possible application in pyrochlore iridates,” Physical Review B 84, 075129 (2011).
- Baireuther et al. (2016) P. Baireuther, J. A. Hutasoit, J. Tworzydło, and C. W. J. Beenakker, “Scattering theory of the chiral magnetic effect in a Weyl semimetal: interplay of bulk Weyl cones and surface Fermi arcs,” New Journal of Physics 18, 045009 (2016).
- Igarashi and Koshino (2017) A. Igarashi and M. Koshino, “Magnetotransport in weyl semimetal nanowires,” Phys. Rev. B 95, 195306 (2017).
- Wang et al. (2017) C. M. Wang, H.-P. Sun, H.-Z. Lu, and X. C. Xie, “3d quantum hall effect of fermi arcs in topological semimetals,” Phys. Rev. Lett. 119, 136806 (2017).
- Hosur et al. (2012) P. Hosur, S. A. Parameswaran, and A. Vishwanath, “Charge Transport in Weyl Semimetals,” Physical Review Letters 108, 046602 (2012).
- Nandy et al. (2017) S. Nandy, G. Sharma, A. Taraphder, and S. Tewari, “Chiral anomaly as origin of planar hall effect in weyl semimetals,” arXiv:1705.09308 [cond-mat] (2017), arXiv: 1705.09308.
- Gorbar et al. (2016) E. V. Gorbar, V. A. Miransky, I. A. Shovkovy, and P. O. Sukhachov, “Origin of dissipative Fermi arc transport in Weyl semimetals,” Physical Review B 93, 235127 (2016).
- Castro et al. (2015) E. V. Castro, M. P. López-Sancho, and M. A. H. Vozmediano, “Anderson localization and topological transition in chern insulators,” Phys. Rev. B 92, 085410 (2015).
- Haake (2001) F. Haake, “Time Reversal and Unitary Symmetries,” in Quantum Signatures of Chaos, Springer Series in Synergetics No. 54 (Springer Berlin Heidelberg, 2001) pp. 15–36, dOI: 10.1007/978-3-662-04506-0_2.
- Groth et al. (2014) C. W. Groth, M. Wimmer, A. R. Akhmerov, and X. Waintal, “Kwant: a software package for quantum transport,” New Journal of Physics 16, 063065 (2014).
- (25) The Python Tight-Binding package is available at this link.
- Takane (2016) Y. Takane, “Disorder effect on chiral edge modes and anomalous hall conductance in weyl semimetals,” Journal of the Physical Society of Japan 85, 124711 (2016).
- Kobayashi et al. (2015) K. Kobayashi, Y. Yoshimura, K.-I. Imura, and T. Ohtsuki, “Dimensional crossover of transport characteristics in topological insulator nanofilms,” Phys. Rev. B 92, 235407 (2015).
- Burkov and Balents (2011) A. A. Burkov and L. Balents, “Weyl Semimetal in a Topological Insulator Multilayer,” Physical Review Letters 107, 127205 (2011).
- Yokomizo and Murakami (2017) K. Yokomizo and S. Murakami, “Topological phases in a weyl semimetal multilayer,” Phys. Rev. B 95, 155101 (2017).
- Datta (1995) S. Datta, Electronic Transport in Mesoscopic Systems (Cambridge University Press, 1995).
- Chen et al. (2015) C.-Z. Chen, J. Song, H. Jiang, Q.-f. Sun, Z. Wang, and X. C. Xie, “Disorder and metal-insulator transitions in weyl semimetals,” Phys. Rev. Lett. 115, 246603 (2015).
- Shapourian and Hughes (2016) H. Shapourian and T. L. Hughes, “Phase diagrams of disordered weyl semimetals,” Phys. Rev. B 93, 075108 (2016).
- Haldane (1988) F. D. M. Haldane, “Model for a Quantum Hall Effect without Landau Levels: Condensed-Matter Realization of the ”Parity Anomaly”,” Phys. Rev. Lett. 61, 2015 (1988).
- Nagaosa et al. (2010) N. Nagaosa, J. Sinova, S. Onoda, A. H. MacDonald, and N. P. Ong, “Anomalous hall effect,” Rev. Mod. Phys. 82, 1539 (2010).
- Todorov (2002) T. N. Todorov, “Tight-binding simulation of current-carrying nanostructures,” Journal of Physics: Condensed Matter 14, 3049 (2002).