Abstract
Magnetic Weyl semimetals are a newly discovered class of topological materials that may serve as a platform for exotic phenomena, such as axion insulators or the quantum anomalous Hall effect. Here, we use angle-resolved photoelectron spectroscopy and ab initio calculations to discover Weyl cones in CoS2, a ferromagnet with pyrite structure that has been long studied as a candidate for half-metallicity, which makes it an attractive material for spintronic devices. We directly observe the topological Fermi arc surface states that link the Weyl nodes, which will influence the performance of CoS2 as a spin injector by modifying its spin polarization at interfaces. In addition, we directly observe a minority-spin bulk electron pocket in the corner of the Brillouin zone, which proves that CoS2 cannot be a true half-metal.
INTRODUCTION
Since the experimental discovery of topological insulators and topological semimetals, there has been a substantial effort to functionalize topological materials for spintronic applications, most prominently in the family of Bismuth chalcogenides (1–4). However, because the field of spintronics predates the first prediction of topological insulators (5, 6) and topological semimetals (7), topological phenomena in many well-studied spintronic materials may have been overlooked during initial investigations of their electronic properties. Here, we reveal that CoS2, a material that has been long studied because of its itinerant ferromagnetism and potential for half-metallicity, actually hosts Weyl fermions and Fermi arc surface states in its band structure close to the Fermi level EF, as well as topological nodal lines below the Fermi level. CoS2 is, therefore, a rare example of the recently discovered class of experimentally verified magnetic topological metals (8–10), which have been proposed to realize new mechanisms of spin-to-charge conversion (11), and are of broader interest for fundamental science, for instance, as a platform to realize axion insulators (7), the intrinsic anomalous Hall effect (12), or the anomalous fractional quantum Hall effect (13).
Beyond the discovery of its topological properties, we also clarify the question of whether CoS2 is a true half-metal. Half-metallic ferromagnets, i.e., materials that are metallic in one spin channel but gapped in the other, have been considered as promising materials for spintronic devices, where they could act as sources of highly spin-polarized charge carriers. There has been a longstanding debate on whether CoS2 or the related alloys Co1−xFexS2 are true half-metals, which has important implications for materials and device engineering. Calculations within the local spin density approximation (LSDA) (14, 15) [see Fig. 1 (A and B)] conclude that CoS2 hosts a minority-spin electron pocket at the R point that leads to a peak in the minority-spin density of states (DOS), which would make CoS2 a minority-spin conductor. This scenario is supported by magnetotransport experiments that suggest a sign flip in the spin polarization upon hole doping with iron that may turn CoS2 into a half-metal (i.e., vanishing minority-spin DOS at EF) due to the depopulation of a minority-spin electron pocket (16–20). On the other hand, calculations based on the generalized gradient approximation (GGA) (15, 21) [see Fig. 1 (B and C)] conclude that half-metallicity is already obtained in the undoped compound, and quantum oscillation experiments were unable to detect the putative minority-spin electron pocket (21). To resolve the debate about the half-metallicity of CoS2, a direct spectroscopic band structure investigation would be clearly desirable, also because the spin polarization measured with transport probes can be reduced by surface defects (22) and can therefore not entirely rule out intrinsic half-metallicity. It should be noted that finite temperature effects can also reduce the spin polarization of nominally half-metallic ferromagnets (23).
(A) Band structure obtained from LSDA. (B) Comparison of DOS from both spin channels for LSDA and GGA. The red solid area shows the majority spin, and the blue solid area shows the minority spins. (C) Band structure obtained from GGA. (D) Illustration of the Weyl points in the bulk Brillouin zone, which are located on the high symmetry planes kx = 0 and ky = 0 parallel to the magnetization direction
Besides the bulk band structure, electronic surface states also influence the spin polarization at the Fermi level, which becomes particularly important for heterostructure interfaces in spintronic devices where half-metals could function as spin injectors. One important class of surface states are topological Fermi arcs in Weyl semimetals, which are considered to be particularly robust against passivation due to the protection by topological bulk invariants. Our LSDA calculations of CoS2 shown in Fig. 1 (A and B) predict Weyl nodes close to the Fermi level on the kx = 0 and ky = 0 planes parallel to the magnetization direction (see Fig. 1D), which should give rise to topological Fermi arc surface states. In addition, we find topological nodal lines a few hundreds of milli–electron volts below the Fermi level, which are protected by a {mz|½ 0 ½} glide operation (also see the Supplementary Materials). Since the concept of topological semimetals and Fermi arc surface states was not yet established at the time, previous theoretical studies of CoS2 overlooked these important features in its band structure, which can have a decisive influence on spin transport properties.
Previous angle-resolved photoelectron spectroscopy (ARPES) experiments (24–26) on CoS2 were performed with photon energies between hv = 20 to 120 eV where the inelastic mean free path of the photoelectrons can be expected to be close to its minimum ~0.5 to 1 nm, resulting in very surface sensitive probes. This implies that the signal from bulk bands is substantially broadened along the momentum direction kz perpendicular to the sample surface (27), which makes the identification of any bulk band dispersion in these experiments challenging, and no minority-spin bulk electron pocket in the vicinity of the R point could be resolved. While previous studies did report the presence of surface states (26), it remained unclear whether these surface states are crossing the Fermi level and whether they are of topological origin.
Here, we overcome these limitations by using complementary bulk-sensitive soft x-ray ARPES (SX-ARPES; hv > 400 eV) and surface-sensitive vacuum ultraviolet ARPES (VUV-ARPES; hv < 120 eV) to disentangle the bulk and surface electronic structure of CoS2. As a result, we are able to observe its elusive minority-spin bulk electron pocket at the R point of the Brillouin zone directly, settling the longstanding debate of half-metallicity in CoS2. We also detect a topological line node located at ~150 meV below the Fermi level, as well as a Weyl cone that suggest a Weyl point slightly above the Fermi level, which gives rise to topological Fermi arc surface states that cross the Fermi level. By comparison with ab initio calculations, we find that these Fermi arcs are of majority spin character and therefore affect the spin polarization at a heterostructure interface.
RESULTS
Our samples were synthesized using chemical vapor transport. We studied their elemental composition with core-level spectroscopy and energy-dispersive x-ray spectroscopy, which confirmed the expected stoichiometry (see Fig. 1E). Powder x-ray diffraction confirmed the previously reported cubic space group 205 and a lattice constant of a = 5.5287 (5) (see Materials and Methods and the Supplementary Materials for more details about synthesis and characterization). We measured a Curie temperature of Tc = 124 K and a saturation magnetization of 0.91 μB/Co at 9 T (which becomes linear at 0.2 T and 0.89 μB), in good agreement with the literature values (28) (see Fig. 1F for magnetization curves). The coercivity in our data is small, about 10 to 15 Oe, which agrees also with the literature values (16).
When cleaving CoS2 for ARPES experiments, we obtained two distinct cleavage planes with the surface normal pointing along the (111) and (100) directions. We performed photon energy–dependent ARPES measurements with soft x-ray photons (hv = 350 to 800 eV) to locate the high-symmetry planes along the kz direction normal to the sample surface (see the Supplementary Materials). The band structure in the kz = π plane containing the R point is for the (111) surface displayed in Fig. 2 (see fig. S9 for a definition of the high-symmetry planes). We can clearly identify circular Fermi surface pockets at the R point in the corner of the Brillouin zone in Fig. 2A. Our calculated Fermi surface (Fig. 2B; using the LSDA) is in good qualitative agreement with the experimental data, confirming the existence of Fermi surface pockets at the R point. When inspecting the experimental band dispersion along the R-X-R direction (Fig. 2C), we see that the circular pockets at the R point are electron-like and are related to another parabolic band with a minimum at around ~0.65 eV by the exchange splitting. The magnitude of the exchange splitting extracted from the energy distribution curve at the R point (Fig. 2D) is ΔE = 0.60 eV (3). Our LSDA calculations of the band dispersion shown in Fig. 2E are in good qualitative agreement with the experimental data and indicate that the observed electron pocket is of minority-spin character, which implies that CoS2 is not a true half-metal. However, the experimentally observed exchange splitting is ~250 meV smaller than in the LSDA calculations, such that the majority spin bands are located closer to the Fermi level in the experiment than expected from the calculations. The data measured on the (100) surface also show electron pockets at the R point and are displayed in the Supplementary Materials.
(A) Experimental Fermi surface was measured on the (111) cleavage plane with photon energy hv = 602 eV and linear-vertical polarization, integrated over 50 meV below the Fermi energy. The red arrow indicates an electron pocket located at the R point. The black arrow indicates the position of the line cut shown in (C). (B) Calculated Fermi surface spectral function A(k,EF) for the same plane as shown in (A) obtained with the LSDA. (C) Line cut along the R-X-R direction as shown in (A); red arrows indicate the electron pockets at the R point. (D) Energy distribution curves for the two R points shown in (E). The black arrow indicates the magnitude of the exchange splitting of ΔE = 0.60(3) eV. (E) Calculated band structure obtained with LSDA; red arrows indicate minority-spin electron pocket.
To search for the topological nodal line and Weyl nodes in CoS2, we also probed the bulk band structure in the kz = 0 plane, containing the Γ point, as illustrated in Fig. 3. Figure 3 (A to D) displays the experimental and calculated Fermi surfaces for the (111) and (100) cleavage planes, which are in good qualitative agreement. Figure 3E shows the band dispersion along the M-Γ direction (black arrow in Fig. 3A), measured on the (111) surface (black arrow in Fig. 3A), and the Γ-M direction, measured on the (100) surface (black arrow in Fig. 3C). Note that the observable bands along these two directions are very different, possibly due to matrix element effects. The line cut obtained from the (111) surface shows a V-shaped feature centered at the M point and a quasi-parabolic band centered at the Γ point. In contrast, the dispersion obtained from the (100) surface shows a single band dispersing in the opposite direction from the quasi-parabolic band. To enhance the contrast of our data, we also show the corresponding second derivative spectrum in Fig. 3F. The combined band dispersion from both surfaces is illustrated in Fig. 3G, which displays the peak positions from a fit of the momentum distribution curves. By comparison with the calculated band dispersion shown in Fig. 3H, we can see that the band crossing between the blue and red bands [from the (111) surface and (100) surface, respectively] is part of a topological line node, while the blue bands form a Weyl cone. The Weyl point that corresponds to the Weyl cone is shown in Fig. 3I, which displays the calculated band dispersion along the M*-Γ direction, where M*
(A) Experimental Fermi surface measured on the (111) cleavage plane with photon energy hv = 512 eV and linear-vertical polarization, integrated over 50 meV below the Fermi energy. The black arrow indicates the position of the line cut shown in (E). (B) Calculated Fermi surface spectral function for the same plane as shown in (A) obtained with LSDA. (C) Experimental Fermi surface spectral function A(k,EF) measured on the (100) cleavage plane with photon energy hv = 475 eV and linear-vertical polarization, integrated over 50 meV below the Fermi energy. The black arrow indicates the position of the line cut shown in (E). (D) Calculated Fermi surface for the same plane as shown in (C) obtained with LSDA. (E) Line cuts along the M-Γ direction from the (111) surface, as shown in (A), and the Γ-M direction from the (100) surface, as shown in (C). (F) Second-derivative spectrum of (E). (G) Result of the momentum distribution curve (MDC) fitting of the bands along the M-Γ-M direction, where blue circles originate from data of the (111) plane, and red circles originate from data of the (100) plane. (H) Calculated band dispersion along the M-Γ direction. (I) Band dispersion along the M*-Γ direction where M*
Fermi arc surface states are a hallmark of Weyl points in topological semimetals. Therefore, the Weyl points in CoS2 must be accompanied by Fermi arc surface states that are connecting the projections of the Weyl points in the surface Brillouin zone. We used surface-sensitive VUV-ARPES to investigate the surface electronic structure of the (100) surface in CoS2, the results of which are displayed in Fig. 4. Figure 4A shows the experimentally obtained Fermi surface, which was measured on a strongly tilted crystal plane. The photon energy dependence of the Fermi surface maps (see the Supplementary Materials) indicates that all Fermi surface pockets measured by VUV-ARPES are surface states and that the signal from bulk states is mostly suppressed in this photon energy range. The positions of the projections of the Weyl point in the (100) surface Brillouin zone are shown in Fig. 1D. From our ab initio calculations, we expect that there should be two Weyl cone projections with total Chern numbers C = ±2 in the vicinity of the
(A) Experimental Fermi surface measured on the (100) cleavage plane with photon energy hv = 100 eV and linear-horizontal polarization. The dashed black line indicates the boundary of the surface Brillouin zone, and the black dashed arrow indicates the momentum direction of the line cuts shown in (C) to (F). (B) Illustration of the Fermi arcs and Weyl point projections in the vicinity of the surface Brillouin zone boundary (black dashed line) at the
DISCUSSION
One may wonder whether the appearance of the conical bulk-band pockets connected by surface states could just be a coincidence, and whether the circular surface state pocket observed here may actually be just a trivial surface state. However, such a coincidence would require an unrealistic degree of fine-tuning, because there are no symmetry constraints that should pin a trivial surface state to the projections of the conical bulk bands. The only plausible explanation is, therefore, that the surface states mentioned above are Fermi arcs, and the observed linear bulk bands are part of a Weyl cone. Although we cannot directly observe the Weyl point since it is located in the unoccupied part of the DOS, the overall agreement between the experimental data and LSDA calculations in conjunction with experimentally observed momentum space alignment between surface states and bulk bands provides sufficient evidence to conclude that CoS2 is a magnetic Weyl semimetal.
Because of the relatively low effective mass of the minority-spin electron pocket at the R point, the minority-spin DOS at the Fermi level is large, and some calculations suggested that CoS2 is a minority-spin conductor (18). Our discovery of the majority-spin Fermi arc surface states implies that the total interface spin polarization of a heterostructure involving CoS2 will be reduced compared to the bulk value. Compared to trivial dangling bond surface states, Fermi arc surface states are more robust against attempts of passivation, because they are protected by the topological invariants (the Chern numbers) of the bulk Weyl points; hence, engineering of the interface potential cannot easily remove the Fermi arcs. On the other hand, hole doping with iron has been suggested to transform CoS2 from a minority-spin to a majority-spin conductor and ultimately to a full half-metal by depopulation of the minority-spin bulk electron pocket (18). Such a transition will be facilitated by the majority-spin Fermi arcs at interfaces in heterostructures since they compensate the DOS of the minority-spin electron pocket. Therefore, CoS2 provides a prime example of a spintronic material whose performance is affected by its topologically nontrivial band structure. Note that nonmagnetic Weyl semimetals have already been successfully implemented for applications in spin-to-charge conversion and magnetic switching (29, 30), and theoretically, proposals have highlighted the importance of topological Fermi arcs for interfacial spin accumulation (31). We therefore expect that our results will inform the design of future spintronic prototypes that use magnetic Weyl-semimetals, for instance, for spin-to-charge conversion (11). There may be magnetic domains on the sample surface that are smaller than the size of our beamspot (diameter of 50 to 70 μm), but on the basis of our ab initio calculations, we believe that the direction of the magnetization vector will have a negligible effect on the interpretation of the Fermi arc structure. This is because spin-orbit coupling is relatively small in CoS2, and its inclusion has only a negligible effect on the band dispersion in our ab initio calculations (see the Supplementary Materials for further discussion). CoS2 is also known as a good catalyst, e.g., for the hydrogen evolution reaction (32), and it has recently been speculated that Fermi arcs in Pt- and Pd-based chiral topological semimetals (33, 34) may play a role in catalysis due to their d-electron character and their robustness against hydrogen passivation (35, 36). Since the Fermi arcs in CoS2 are also derived from bulk bands of d orbital character (see the Supplementary Materials), they may contribute to the catalytic performance of CoS2.
MATERIALS AND METHODS
Sample growth
Single crystals of CoS2 were synthesized using chemical vapor transport. The elements cobalt (200 mg, 1 eq, 99.5%; Sigma-Aldrich) and sulfur (217 mg, 2 eq, 99.5%; Alfa Aesar) were mixed with 50 mg of iodine (99.999%; Sigma Aldrich) and sealed in an evacuated quartz glass ampoule. The ampoule was slowly (within 12 hours) heated to 1000°C and kept there for 140 hours. The growth was along a natural temperature gradient of a tube furnace. An increase in crystal size can be accomplished using an additional slow cooling step from 1000° to 800°C in 100 hours.
ARPES measurements
SX-ARPES measurements were performed at the SX-ARPES endstation (37) of the ADRESS beamline (38) at the Swiss Light Source, Switzerland, with a SPECS analyzer with an angular resolution of 0.1°. The photon energy varied from 350 to 1000 eV, and the combined energy resolution ranged between 50 and 150 meV. The temperature during sample cleaving and measurements was below 20 K, and the pressure was lower than 1 × 10−10 mbars. The increase of the photoelectron mean free path in the soft x-ray energy range results, by the Heisenberg uncertainty principle, in a higher kz resolution of the ARPES experiment compared to measurements at lower photon energies (39), which was critical to measure the bulk band structure of CoS2.
VUV-ARPES measurements were performed at the high-resolution ARPES branch line of the beamline I05 at the Diamond Light Source, UK (40). Measurements at the high-resolution branch were performed with a Scienta R4000 analyzer, and a photon energy range between 90 and 130 eV, at a temperature below 20 K. Measurements in the VUV-ARPES regime are more surface sensitive than SX-ARPES and therefore most suitable to image the Fermi arcs in CoS2.
Ab initio calculations
We used density functional theory as implemented in the Vienna Ab Initio Simulation Package. For the GGA calculations, the exchange correlation term is described according to the Perdew-Burke-Ernzerhof prescription together with projected augmented wave (PAW) pseudopotentials, while for the LSDA, the exchange correlation term is described according to the Dudarev simplified, rotationally invariant approach, together with PAW pseudopotentials. The kinetic energy cutoff was set to 400 eV. For the self-consistent calculation, a grid of 7 × 7 × 7 k points was used. For DOS calculation, a grid of 11 × 11 × 11 k points was used, with 1000 energy points. The spectral function of the Fermi surfaces shown in Figs. 2 and 3 was computed with the programs Wannier90 (41) and WannierTools (42). Using ab initio calculation of bands, we constructed maximally localized Wannier functions considering d orbitals for cobalt sites and p orbitals for sulfur sites using the Wannier90 software (41). We computed an effective tight-binding model that accurately reproduces the band structure of the bulk system. We then computed the surface spectrum by the method of iterative Green’s function as implemented in WannierTools (42) software for a slab of 50 unit cells.
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/51/eabd5000/DC1
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).