Arxiv:2008.10628V3 [Cond-Mat.Str-El] 14 Jan 2021
Total Page:16
File Type:pdf, Size:1020Kb
Prediction of Spin Polarized Fermi Arcs in Quasiparticle Interference of CeBi Zhao Huang,1 Christopher Lane,1, 2 Chao Cao,3 Guo-Xiang Zhi,4 Yu Liu,5 Christian E. Matt,5 Brinda Kuthanazhi,6, 7 Paul C. Canfield,6, 7 Dmitry Yarotski,2 A. J. Taylor,2 and Jian-Xin Zhu1, 2, * 1Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 2Center for Integrated Nanotechnology, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 3Department of Physics, Hangzhou Normal University, Hangzhou 310036, China 4Department of Physics, Zhejiang University, Hangzhou 310013, China 5Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA 6Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA 7Department of Physics and Astronomy, Iowa State University, Ames, Iowa 50011, USA (Dated: January 15, 2021) We predict that CeBi in the ferromagnetic state is a Weyl semimetal. Our calculations within density func- tional theory show the existence of two pairs of Weyl nodes on the momentum path (0,0,kz) at 15 meV above and 100 meV below the Fermi level. Two corresponding Fermi arcs are obtained on surfaces of mirror- symmetric (010)-oriented slabs at E = 15 meV and both arcs are interrupted into three segments due to hy- bridization with a set of trivial surface bands. By studying the spin texture of surface states, we find the two Fermi arcs are strongly spin-polarized but in opposite directions, which can be detected by spin-polarized ARPES measurements. Our theoretical study of quasiparticle interference (QPI) for a nonmagnetic impurity at the Bi site also reveals several features related to the Fermi arcs. Specifically, we predict that the spin polariza- tion of the Fermi arcs leads to a bifurcation-shaped feature only in the spin-dependent QPI spectrum, serving as a fingerprint of the Weyl nodes. I. INTRODUCTION cols. Only very recently have Co3Sn2S2 and Co2MnGa been reported as possible magnetic Weyl semimetals with six Weyl 4–6 Despite being predicted more than 90 years ago as an el- nodes and nodal lines, respectively . This is far from enough egant special solution of the relativistic wave equation, Weyl to explore the wide range of possible exotic states. Inter- fermions have never been observed in nature. However, just estingly, recent experimental reports of a pronounced nega- as the sun was setting on the Weyl fermion, they were ob- tive magnetoresistance and observed band inversion in cerium monopnictides35–37 and the geometrically frustrated Shastry- served as a quasiparticle in the electronic structure of TaAs, 38 thereby pushing Weyl physics back to the vanguard and invig- Sutherland lattice GdB4 could provide clues to a possible orating the search for other exotic relativistic particles within new material family harboring Weyl physics. solids. Weyl fermions appear in condensed matter systems when a single Dirac cone is split into two by the breaking of In this paper, we argue that CeBi in the ferromagnetic either time-reversal1–6 or inversion symmetry7–13. This pro- (FM) state is a magnetic Weyl semimetal using first-principles cess generates a pair of cone-like features with opposite chi- LDA+U calculations and we predict robust quasiparticle inter- rality defined by their Berry curvature14. As a direct con- ference (QPI) signatures of spin polarized Fermi arcs directly sequence, anomalous surface states appear connecting Weyl accessible by spin resolved scanning tunneling spectroscopy nodes of opposite topological charge, providing an experi- (STS). By aligning the Ce magnetic moments along the c- mentally accessible fingerprint of the underlying non-trivial axis, we find two pairs of Weyl nodes on the zone diagonal 4,15–21 band topology . along the kz-axis at 15 meV above and 100 meV below the Recently, magnetic Weyl semi-metals have been gaining Fermi level. Then by examining a (010)-oriented slab which interest for providing a new platform to study the inter- is mirror symmetric, two Fermi arcs of opposite spin polar- play between chirality, magnetism, correlation, and topolog- ization are observed on both the top and bottom of surfaces, ical order, thus opening up new routes to novel quantum stretching across the Brillouin zone. These Fermi arcs then states, spin polarized chiral transport22,23, and exotic opti- hybridize with a set of trivial surface states producing a spin cal phenomena24–27. Compared with inversion breaking Weyl vortex encircling the M¯ point in the Brillouin zone. Finally, arXiv:2008.10628v3 [cond-mat.str-el] 14 Jan 2021 materials, magnetic Weyl materials have the advantage of re- by considering a weak potential scatterer at the Bi site on the ducing the minimum number of allowed Weyl nodes from surface, we predict several features directly related to the scat- four to two. This allows simpler band structures thereby tering among the Fermi arcs in the QPI spectra. Specifically, facilitating clearer comparisons with theoretical predictions, there is a bifurcation-shaped feature originating from the scat- and more robust applications in spintronics28–30 and quantum tering between the two symmetric outer segments of the spin- computation31–34. Another advantage is the magnetic field up Fermi arc. We propose this feature as a direct indicator of tunability which provides us a strong tool on controlling the the Weyl nodes for future experiments. Moreover, due to the band structure and related electromagnetic functionality. isolation of the Weyl nodes near the Fermi level, we find CeBi To date very few magnetic Weyl materials have been syn- to be a robust platform for theoretical and experimental anal- thesized, demonstrating an urgent and important need for ysis of Weyl physics with minimal interference from trivial more theoretical material predictions and robust growth proto- bulk states. 2 II. METHOD First-principles band structure calculations were carried out within density functional theory framework using the generalized-gradient approximation (GGA) as implemented in the all-electron code WIEN2K39, which is based on the augmented-plane-wave + local-orbitals(APW+lo) basis set. Spin-orbit coupling was included in the self-consistency cy- cles. Furthermore, the Hubbard Coulomb interaction on Ce- 4 f electrons was incorporated to ensure a fully localized Ce-4 f state, which also gives a total magnetic moment of 2µB/Ce consistent with earlier studies oncerium monopnic- tides35,40,41. Although these compounds exhibit several field- FIG. 1. (a) Crystal structure of CeBi with the spin polarization along induced magnetic phases, a ferromagnetic-like state can be the c-axis direction. Here we display a body-centered tetragonal cell stabilized using state-of-the-art cryogenic STM in 3D vector for an otherwise FCC crystal structure of CeBi in bulk. (b) DFT magnetic fields above 4 Tesla35,41,42. The topological analysis electronic band structure of CeBi in the ferromagnetic state with a was performed by employing a real-space tight-binding model red circle highlighting a pair of Weyl nodes near the Fermi level. (c) Same as (b) except a zoom-in around Fermi level along Γ Z with the Hamiltonian, which was obtained by using the wien2wannier − 43 character of various bands indicated. (d) FCC Brillouin zone of FM interface . Bi-6p and Ce-5d states were included in gen- CeBi with the various high-symmetry points marked. The locations erating Wannier functions. It is worth mentioning that, the of the Weyl nodes is given by teal (orange) dots denoting the Weyl DFT+U method for a non-magnetic phase will still leave node as a source (sink) of Berry curvature. f -electron pinned around the Fermi level, which is com- pletely different from the spin-polarized case44. Therefore, the DFT+U approach is not appropriate to describe strongly correlated materials in a paramagnetic state, for which one function composed of 6 Bi-6p and 10 Ce-5d orbitals, and needs to resort to more powerful methods like DFT+DMFT.45 T is the scattering T-matrix containing all multiple scat- To identify the presence of a pair of Weyl nodes in a given tering effects off a single-site impurity. Since the infi- band structure we must check if the band crossings are topo- nite sum of multiple scattering terms form a geometric se- logically trivial or non-trivial. That is, we must calculate ries, the T-matrix can be compactly written as T(ω) = 2 2 1 the Berry curvature at and surrounding the nodal points to I Vimp(2π)− d kG(k,ω) − Vimp, where Vimp is the im- see if the field is divergent or not. To do so, we calculate purity− potential. Here we considered a Born approximation of R the Berry curvature Ωn,αβ (k) as implemented in the Wannier- the impurity scattering such that T = Vimp and further chose 46–48 Tools package , for a given band n and momenta k as Vimp to be a constant diagonal matrix with non-zero elements only for p (d) orbitals for a Bi (Ce) impurity site. Our gen- vnm,α (k)vmn,β (k) eral methodology for QPI prediction in CeBi is justified by Ωn,αβ (k) = 2Im ∑ 2 (1) − m=n [εmk εnk] the good agreement between calculated and experimental STS 6 − data on the (001)-oriented slab, shown in detail in the Supple- 50 where εnk is the eigenvalue of Hamiltonian H for band n and mental Material (SM) . momenta k and the matrix elements of the Cartesian velocity operators are given by III. BAND STRUCTURE AND WEYL NODES ∂Hˆ (k) v , (k) = φ φ . (2) nm α h nk ∂k mki α Figure1 (left panel) shows the crystal structure of CeBi in the FM phase along with its corresponding face-centered- Finally the Berry curvature vector is given by cubic (FCC) Brillouin zone, where various high symmetry lines are marked.