k + f− f Twisted Fermi surface of a thin-film Weyl semimetal N. Bovenzi,1 M. Breitkreiz,1 T. E. O'Brien,1 J. Tworzyd lo,2 and C. W. J. Beenakker1 1Instituut-Lorentz, Universiteit Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 2Institute of Theoretical Physics, Faculty of Physics, kz University of Warsaw, ul. Pasteura 5, 02{093 Warszawa, Poland (Dated: September 2017) + k (k) kz (k) The Fermi surface of a conventional two-dimensional electron gas is equivalent to a circle,z− up to smooth deformations that preserve the orientation of the equi-energy contour. Here we show + i− i that a Weyl semimetal confined to a thin film with an in-plane magnetization and broken spatial inversion symmetry can have a topologically distinct Fermi surface that is twisted into a figure-8 | opposite orientations are coupled at a crossing which is protected up to an exponentially small gap. The twisted spectral response to a perpendicular magnetic field B is distinct from that of a deformed Fermi circle, because the two lobes of a figure-8 cyclotron orbit give opposite contributions to the Aharonov-Bohm phase. The magnetic edge channels come in two counterpropagatingp types, 2 p a wide channel of width βlm / 1=B and a narrow channel of width lm / 1= B (with lm = ~=eB the magnetic length and β the momentum separation of the Weyl points). Only one of the two is transmitted into a metallic contact, providing unique magnetotransport signatures. I. INTRODUCTION a) ⌫ =0 b) ⌫ =1 c) ⌫ =2 ⌫1 =+1 ⌫1 =+1 The Fermi surface of degenerate electrons separates filled states inside from empty states outside, thereby ⌫2 =+1 ⌫ = 1 governing the electronic transport properties near equi- 2 − librium. In a two-dimensional electron gas (2DEG) the Fermi surface is a closed equi-energy contour in the mo- mentum plane. It is a circle for free electrons, with de- formations from the lattice potential such as the trigonal FIG. 1: Three oriented contours (black curves) with turning warping of graphene or the hexagonal warping on the sur- number ν = 0; 1; 2. The red segments show the uncrossing deformation that removes a self-intersection without changing face of a topological insulator [1]. These are all smooth the total turning number ν = P ν . deformations which do not change the orientation of the i i Fermi surface: The turning number is 1, meaning that the tangent vector makes one full rotation as we pass along the equi-energy contour. spatial inversion symmetry can have ν = 0: if the Fermi The turning number level lies in between the two Weyl points the circular Fermi surface is twisted into a figure-8 with zero total 1 I curvature [4]. ν = dl; (1) The self-intersection introduced when the Fermi level 2π Γ C passes through a Weyl point, to ensure that ν + s re- j j defined as the contour integral of the curvature in mains odd, is a crossing of Fermi arcs on the top and units of 2π, identifies topologically distinct deformationsC bottom surfaces of the thin film (width W ). These have of the circle in the plane, socalled \regular homotopy a penetration depth ξ0 into the thin film that can be classes" [2]. A theorem going back to Gauss [3] says that much less than the Fermi wavelength of the bulk states, a contour Γ with turning number ν has s ν 1 so that we can be in the 2D regime of a single occupied self-intersections and that the sum ν + s must≥ j bej − an subband [5] without appreciable overlap of the surface odd integer. Fig. 1 shows examplesj j of contours with states [6{8]. The effect of a nonzero surface state overlap W/ξ0 ν; s = 0; 1 , 1; 0 , and 2; 1 . is to open up an exponentially small gap δk e− in f g f g f g f g / The turning number is preserved by any smooth defor- the figure-8, as in Fig. 1a. arXiv:1707.01038v4 [cond-mat.mes-hall] 8 Feb 2018 mation of the contour. This includes socalled \uncross- In a perpendicular magnetic field B the signed area en- ing" deformations [2]: As illustrated in Fig. 1, uncrossing closed by the Fermi surface is quantized in units of 2π=l2 , p m breaks up a self-intersecting contour Γ into a collection of with lm = ~=eB the magnetic length. A figure-8 Fermi nearly touching oriented contours Γi, with turning num- surface of linear dimension kF has a signed area much P 2 bers νi. The total turning number ν = i νi is invariant smaller than kF, because the upper and lower loops have against uncrossing deformations, which is another result opposite orientation. We find that this twisted Fermi sur- 2 due to Gauss [3]. face produces edge states of width kFlm | much wider All familiar 2D electron gases belong to the ν = 1 than the usual narrow quantum Hall edge states of width j j universality class. Here we show that a thin-film Weyl lm. The wide and the narrow edge states are counter- semimetal with an in-plane magnetization M and broken propagating: if the wide channel moves parallel to M, 2 the narrow channel moves antiparallel. An applied volt- B. Dispersion relation age selectively populates one of the two types of edge 2 2 states, resulting in a conductance of e =h instead of 2e =h The Schr¨odingerequation H = E can be solved | even though there are two conducting edges. analytically in the low-energy regime by linearizing in The outline of the paper is as follows. In the next kx and applying the effective mass approximation [12] section we formulate the problem, on the basis of a two- kx i@=@x. Integration of the resulting first-order band model Hamiltonian [9, 10], and calculate the band differential7! − equation in x gives structure in a slab geometry. The way in which the Fermi arcs reconnect with the bulk Weyl cones is described ex- ixΞ 1 (x) = e (0); Ξ = σx[E H(0; ky; kz)]: (4) actly by a simple transcendental equation (Weiss equa- tx − tion). The Fermi surface in the thin-film regime is calcu- To ensure that an eigenstate of H satisfies the boundary lated in Sec. III, to show the topological transition from condition (3), we require that turning number 1 to turning number 0 when the Fermi level passes through a Weyl point. In Sec. IV we calcu- iW Ξ 1 late the edge states in a perpendicular magnetic field, by e = 0; = ; σy = : (5) h−| |−i |±i i |±i ±|±i semiclassical analytics and comparison with a numerical ± solution. The implications of the two types of counter- This reduces to the following dispersion relation for propagating edge channels for electrical conduction are E(ky; kz): investigated in Sec. V. We conclude with an overview 2 2 2 2 2 of possible experimental signatures of the twisted Fermi (E λ sin kz) ty sin ky mk = q ; (6) surface. − − − with transverse wave number q given by mk tan(W q=tx) + 1 = 0: (7) II. WEYL SEMIMETAL CONFINED TO A SLAB q In the mass term mk we should set kx = 0, as required A. Two-band model by the linearization in kx. For imaginary q = iκtx=W the transcendental equa- We consider the two-band model Hamiltonian of a tion (7) takes the form Weyl semimetal [9, 10], γ W m tanh κ = 1; γ = k ; (8) κ − t H(k) = txσx sin kx + tyσy sin ky + mkσz + λσ0 sin kz; x mk = tz(cos β cos kz) + t0(2 cos kx cos ky): (2) which is known as the Weiss equation in the theory of − − − ferromagnetism [13]. A unique solution with κ 0 exists ≥ The Pauli matrices are σ , α x; y; z , with σ the for γ 1, given by a generalized Lambert function [14, α 0 ≥ 2 2 unit matrix, acting on a hybrid2 f of sping and orbital 15]: degrees× of freedom. The momentum k varies over the κ = 1 (2γ; 2γ; 1): (9) Brillouin zone kα < π of a simple cubic lattice (lattice 2 W − − constant a 1,j andj we also set 1). The two Weyl 0 ~ A representative band structure is shown in Fig. 2. points are at≡ the momenta k = (0; 0;≡ K), K β, and at energies E = E , E λ sin β, displaced± along≈ the k - ± 0 0 ≈ z axis by the magnetization M = βz^ and displaced along C. Weyl cones and Fermi arcs the energy axis by the strain λ. Time-reversal symmetry and spatial inversion symmetry are broken by β and λ, respectively. In the large-W limit of a thick slab, Eq. (7) can be solved separately for the bulk Weyl cones and the sur- We take a slab geometry, unbounded in the y{z plane face Fermi arcs. We thus recover the familiar dispersion and confined in the x-direction between x = 0 and relations in the bulk and surface Brillouin zones of a Weyl x = W . The magnetization along z is therefore in the semimetal [16{19]. plane of the slab. We impose the infinite-mass boundary The bulk states have wave number q mk , quan- condition [11] on the wave function , 1 j j tized by q = (n + 2 )πtx=W , n = 0; 1; 2;:::, with disper- ( sion at x = 0; q σy = − (3) (n) 1 2 2 2 2 2 + at x = W: E = (n + ) (πtx=W ) + t sin ky + m bulk ± 2 y k + λ sin kz: (10) This boundary condition corresponds to a mass term m0(x)σz in H that vanishes inside the slab and tends The distinguishes the upper and lower halves of the to + outside.