Introduction {#Sec1}
============

The ordinary Hall effect that arises from the Lorentz force deflecting the moving charge carriers has been well understood^[@CR1]^. In contrast, the anomalous Hall effect (AHE)^[@CR2],[@CR3]^ has attracted tremendous interests because of the fundamental physics and great potential in technical application^[@CR4],[@CR5]^, but the microscopic origin of AHE is still not fully solved. The key issue is whether the effect is intrinsic or extrinsic. It is now recognized that there are three mechanisms to account for the AHE^[@CR6]^. One is the extrinsic mechanism related to the scattering affected by the spin--orbit interaction, i.e., the skew scattering and side jump mechanisms^[@CR7]--[@CR9]^. The skew scattering model predicts that the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{{A}}$$\end{document}$ is linearly proportional to *ρ*~*xx*~, whereas, the side jump models give $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}} \propto \rho _{{{xx}}}^{\mathrm{2}}$$\end{document}$. The other is intrinsic Kaplus--Luttinger (KL) mechanism related to spin--orbit interaction of Bloch electronic bands, originally proposed by Karplus and Luttinger, which also gives $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}} \propto \rho _{{{xx}}}^{\mathrm{2}}$$\end{document}$^[@CR10]^. Importantly, recent studies indicate that there is an intimate relation between the AHE and the Berry curvature of occupied electronic Bloch states^[@CR11]--[@CR14]^.

Recently-discovered topological semimetals (TSMs) are characterized by the topologically robust or symmetry-protected bulk band crossings near the Fermi energy (*E*~F~)^[@CR15]--[@CR18]^. Classified by the degeneracy of nodes, Dirac and Weyl semimetals (WSMs) with the 4- and 2-fold degenerate Dirac or Weyl points, respectively, have been theoretically predicted and experimentally confirmed^[@CR18]--[@CR22]^. Especially, in a magnetic WSM with broken time-reversal symmetry (TRS), Weyl node can be seen as a magnetic monopole in momentum space^[@CR23]^, around which there are significant contributions to Berry curvature^[@CR24]^. Thus, if the position of Weyl nodes are proper, especially close to *E*~F~, and the individual Fermi surface (FS) sheets have nonzero Chern numbers, there should be large intrinsic AHE in magnetic WSMs nearly proportional to the distance between a pair of Weyl nodes with opposite chirality^[@CR24]^. On the other hand, the half-metallic ferromagnets (HMFMs) can also exhibit AHE, such as ferromagnetic (FM) Heusler alloy Co~2~MnGa, NiMnSb and related intermetallic compounds^[@CR25],[@CR26]^. Different from traditional ferromagnets, the HMFMs have attracted extensive interests because of the (nearly) 100% spin polarization of conduction electrons at *E*~F~^[@CR27]^. They can be perfectly applied to spintronic devices.

In this work, we study the AHE and electronic structure of HMFM Co~3~Sn~2~S~2~ with kagome layers of Co. As previous experimental results reported on polycrystalline samples, the compound orders ferromagnetically at *T*~C~ = 177 K and the spontaneous magnetic moment is about 0.3 μ~B~/Co^[@CR28]--[@CR30]^. By employing systematic first-principles calculations, transport, and angle-resolved photoemission spectroscopy (ARPES) measurements on Co~3~Sn~2~S~2~ single crystals, we reveal the quadratically scaling relationship between anomalous Hall resistivity $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}}$$\end{document}$ and *ρ*~*xx*~, and the linear behavior between intrinsic anomalous Hall conductivity (AHC) and the magnetization. Furthermore, the agreement between the experimental band structures and theoretical calculations support the existence of magnetic Weyl fermions in Co~3~Sn~2~S~2~. Our findings suggest the main contribution of observed large AHE originates from the intrinsic mechanism, which is intimately related to the Weyl nodes near *E*~F~.

Results {#Sec2}
=======

Structure and Ferromagnetic state in Co~3~Sn~2~S~2~ {#Sec3}
---------------------------------------------------

Co~3~Sn~2~S~2~ is crystalized in a hexagonal lattice with space group of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R\bar 3m$$\end{document}$ (No. 166). As shown in Fig. [1a](#Fig1){ref-type="fig"}, the crystal structure of Co~3~Sn~2~S~2~ is composed of the slabs of CoSn~4~S~2~ octahedra stacking in a hexagonal packing (A-B-C fashion) along the *c* axis. Each Co atom is surrounded by four Sn and two S atoms, forming a distorted octahedron. The Co and S atoms are located at Wyckoff position 9*e* (1/2, 0, 0) and 6*c* (0, 0, 0.219), respectively. The CoSn~4~S~2~ octahedra connect each other along *ab* plane by face-sharing and along *c* axis by corner-sharing. On the other hand, the Co atoms form perfect kagome layer with corner-sharing triangles of Co atoms in Co--Sn layer. There are two kinds of Sn sites. A half of Sn atoms lie in the centers of the kagome hexagons (Sn2 sites, Wyckoff position 3*a* (0, 0, 0)) (Fig. [1b](#Fig1){ref-type="fig"}) and another half of Sn atoms located between the Co--Sn bilayers (Sn1 sites, Wyckoff position 3*b* (0, 0, 1/2)), connecting adjacent Co--Sn2 layers. S atoms are located below and above the Co--Sn2 layers. For the plate-like Co~3~Sn~2~S~2~ single crystal, the *c* axis is perpendicular to the crystal surface (Supplementary Fig. [1](#MOESM1){ref-type="media"}).Fig. 1Structure, magnetization and longitudinal resistivity of Co~3~Sn~2~S~2~. **a**, **b** Crystal structure of Co~3~Sn~2~S~2~ and kagome layer made of Co atoms. The small blue and yellow balls represent Co and S atoms, respectively, and the big green and red balls represent Sn atoms at Sn~1~ and Sn~2~ sites, respectively. **c** Temperature dependence of magnetic susceptibility *χ*(*T*) with ZFC and FC modes at *μ*~0~*H* = 1 T for *H*\|\|*c*. Inset: field dependence of magnetization *M*(*μ*~0~*H*) at 5 and 300 K for *H*\|\|*c*. **d** Longitudinal resistivity *ρ*~*xx*~(*T*, *μ*~0~*H*) as a function of temperature *T* at *μ*~0~*H* = 0 and 9 T along the *c* axis

The *χ*(*T*) curves with zero-field-cooling (ZFC) and field-cooling (FC) modes at *μ*~0~*H* = 1 T for *H*\|\|*c* increase rapidly when *T* is below the Curie temperature *T*~C~ (\~174 K) (Fig. [1c](#Fig1){ref-type="fig"}), well agree with previous results in the literatures^[@CR28]--[@CR30]^. Moreover, the ZFC and FC *χ*(*T*) curves overlap each other very well, suggesting that magnetic domains have aligned with the direction of external field when *μ*~0~*H* = 1 T. The field dependence of magnetization *M*(*μ*~0~*H*) at *T* = 5 K for *H*\|\|*c* further confirms the ferromagnetism at low temperature (inset of Fig. [1c](#Fig1){ref-type="fig"}). The *M*(*μ*~0~*H*) curve shows a pronounced hysteresis. The hysteresis loop is square and changes its direction at very little coercive field *μ*~0~*H*~*c*~ close to 0.08--0.1 T. The saturated magnetization *M*~s~ is about 0.3 μ~B~/Co, consistent with previous results^[@CR28],[@CR30]^. In contrast, when $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$T \gg T_{\mathrm{C}}$$\end{document}$, the hysteresis behavior vanishes and the *M*(*μ*~0~*H*) curve only shows a paramagnetic behavior (inset of Fig. [1c](#Fig1){ref-type="fig"}).

As shown in Fig. [1d](#Fig1){ref-type="fig"}, the zero-field longitudinal resistivity *ρ*~*xx*~(*T*, 0) of Co~3~Sn~2~S~2~ shows a metallic behavior in the whole measuring temperature range. A kink can be clearly observed in the *ρ*~*xx*~(*T*, 0) curve around *T*~C~. The resistivity decreases rapidly with decreasing temperature below *T*~C~ due to the decrease of spin disorder scattering. When the external magnetic field is applied (*μ*~0~*H* = 9 T), the metallic behavior of *ρ*~*xx*~(*T*) has almost no change, but the kink of *ρ*~*xx*~(*T*) near *T*~C~ becomes smooth and the staring temperature where the *ρ*~*xx*~(*T*) curve changes the slope shifts to higher temperature.

Calculated electronic structure, Weyl nodes and AHC {#Sec4}
---------------------------------------------------

The results of the first principle calculation indicate the half-metallic feature of Co~3~Sn~2~S~2~ at FM state (Fig. [2a](#Fig2){ref-type="fig"}) and reveal that the bands dispersion near the *E*~F~ are contributed mainly by the 3*d* orbitals of Co with the polarized magnetic momentum about 0.33 μ~B~/Co irrespective of the presence of spin--orbit coupling (SOC). Co~3~Sn~2~S~2~ belongs to the type *I*~*A*~ half metallic ferromagnets and the spins of all electrons are fully polarized along the up direction (defined in *c* axis)^[@CR31]^. When the SOC is absent, the valence and conductance bands inverse near the *L* point and form a nodal ring (signed by the dashed circles in Fig. [2a](#Fig2){ref-type="fig"}), which is protected by the mirror symmetry of the plane M(010) shown in Fig. [2f](#Fig2){ref-type="fig"}. When the SOC is further included, it will decay into a pair of Weyl points (WPs) with opposite chirality off the high-symmetry lines (Fig. [2b](#Fig2){ref-type="fig"}). The precise positions of the WPs have been obtained by calculating the Wilson-loop evolution^[@CR32]^ on a much denser *k*-grid in the Brillouin zone (BZ). The evolution of the Wannier charge center (WCC) on a sphere that enclosing a WP in Brillion Zone is demonstrated in Supplementary Fig. [2](#MOESM1){ref-type="media"}. There are three pairs of WPs in the BZ in total related by *C*~3~ rotation symmetry and inversion symmetry, as illustrated in Fig. [2c](#Fig2){ref-type="fig"} and Supplementary Table [1](#MOESM1){ref-type="media"}. From the energy dispersion along the direction connecting the WPs of W1 and W2, the splitting distance of them is about 0.3 Å^−1^ (Fig. [2d](#Fig2){ref-type="fig"}), which is far enough to guarantee the robustness of the WPs. We have also calculated the Fermi surface (FS) of (001) surface at the energy of WPs with Green function method^[@CR33],[@CR34]^. As shown in Fig. [2c](#Fig2){ref-type="fig"}, the pairs of the WPs with opposite chirality projected on (001) surface are indicated by green and white dots, respectively. Since all of the WPs are mixed with bulk states, it is not easy to clarify the details of how the Fermi arcs connecting the WPs although there are arc-like FSs with large weight of local density of states on surface and preserving *C*~3~ rotation symmetry. In order to distinguish the topological surface states in Fig. [2c](#Fig2){ref-type="fig"}, we have redefined a *k*-dependent Fermi energy *E*~F~(*k*) = *α*\[*E*~*N*~(*k*) + *E*~*N* + 1~(*k*)\]/2, where 0 ≤ *α* ≤ 1. *E*~*N*~(*k*) is the low branch of the two bands crossing at the Weyl nodes and *E*~*N* + 1~(*k*) is the up branch. Thus, if *α* = 1, the *k*-dependent *E*~F~(*k*) will lead to an ideal Weyl semimetal state with and only with Weyl nodes appearing at the Fermi level. When *α* ranges from 0.1 to 0.9 (as shown in Supplementary Fig. [3](#MOESM1){ref-type="media"}), we can gradually remove the projected bulk state and identify the Fermi-arc patterns from the other topologically trivial surface states. In Supplementary Fig. [4](#MOESM1){ref-type="media"}, the calculations using unrelaxed experimental crystal structure are also shown. The Weyl nodes and magnetic properties are the same as those calculated from relaxed structure. But the WPs and Fermi arcs altered slightly and separated from the projected bulk states. Both the Fermi arcs and the corresponding topological surface states can be identified.Fig. 2Band structure obtained from first-principle calculations. **a** The spin-resolved band structure without SOC. **b** The band structure with SOC included. The definition of high-symmetry points are shown in Fig. [3a](#Fig3){ref-type="fig"}. **c** Fermi surface on the (001) surface cutting on the energy of WPs. The white and green dots represent the projected WPs with positive and negative chirality respectively. **d** The energy dispersion along the direction connecting the WPs W1 and W2 in **c**. **e** The intrinsic anomalous Hall conductivity with the magnetic momentum of 0.09, 0.15, 0.19, and 0.33 μ~B~/Co at 0 K. **f** The evolution of WPs in the BZ with the magnetic momentum decreased from 0.33 to 0.09 μ~B~/Co. The directions of the evolution are indicated by the red arrows. The *E*~F~ is set to 0

The existence of WPs implies that there would be a large intrinsic AHE in Co~3~Sn~2~S~2~. The AHC *σ*~*xy*~ obtained from the integral of Berry curvature along *k*~*z*~ in the BZ has a very large value of about −1310 Ω cm^−1^ at 0 K (Fig. [2e](#Fig2){ref-type="fig"}). Moreover, because the time-reversal symmetry is broken in the FM state, the AHE has a direct relation with the intrinsic magnetic momentum. We constrained the local magnetic momentum as 0.09, 0.15, 0.19, and 0.33 μ~B~/Co and as illustrated in Fig. [2e](#Fig2){ref-type="fig"}, the AHC decreases almost linearly with the decreasing of magnetic momentum. Meanwhile, the WPs evolve monotonically with the magnetic momentum decreasing from 0.33 to 0.09 μ~B~/Co as indicated by the red arrows in Fig. [2f](#Fig2){ref-type="fig"}. The WPs with opposite chirality approach each other along *k*~*z*~ direction when the spin polarization is reduced. In addition, the AHC of magnetic Weyl semimetals can also be related to the distance of WPs. We can estimate the intrinsic AHC by the formula of *σ*~*xy*~ = *Ke*^2^/4*π*^2[@CR24]^, where *K* is the distance of WPs with opposite chirality along the direction of magnetization. The result also confirms the decrease of intrinsic AHC with decreasing magnetization strength (Supplementary Fig. [6](#MOESM1){ref-type="media"}). However, before the pair of Weyl nodes meet each other to annihilate, the two bands composing them become spin degenerate with their time-reversal symmetric partners due to Kramer degeneracy in paramagnetic state.

Experimental electronic structure observed by ARPES measurements {#Sec5}
----------------------------------------------------------------

In order to further uncover the electronic origin of AHE in Co~3~Sn~2~S~2~, we performed ARPES experiments on the (001) surface of the samples. The overall band structures are summarized in Fig. [3](#Fig3){ref-type="fig"}. The experimental FSs in Fig. [3b](#Fig3){ref-type="fig"} are in good agreement with the projections of calculated bulk FSs on the (001) surface in Fig. [3c](#Fig3){ref-type="fig"}. The spindle-shaped FSs centered at $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline M$$\end{document}$ points, the triangle-shaped FSs centered at $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline K$$\end{document}$ point, and the ring-like ones connecting the spindle-shaped FSs are all captured by ARPES. Specifically, the agreement of the spindle-shaped FSs is significant since the projection of bands on the (001) surface that creates the WPs are located along the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline {\mathrm{\Gamma }} - \overline M$$\end{document}$ direction, i.e., the (001)-surface projections of WPs W1 and W2 are located along this direction as illustrated in Fig. [2c](#Fig2){ref-type="fig"}, although the WPs are buried into bulk continuum states at the chemical potential.Fig. 3Near-*E*~F~ ARPES spectra compared with first-principle calculations of Co~3~Sn~2~S~2~. **a** Schematic primitive BZ and 2D projected BZ of the (001) surface. Blue lines parallel to the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline {\mathrm{\Gamma }} - \overline M$$\end{document}$ direction indicate the momentum locations of calculated band structures in **g**, with the *k*~*z*~ momenta of 0.96, 0.92, 0.88, and 0.84 *π* from top to bottom, respectively. Cuts 3 and 4 illustrate the momentum locations of the predicted WPs, which are parallel to the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline {\mathrm{\Gamma }} - \overline M$$\end{document}$ direction in the *k*~z~ = ± 0.1455 *π* planes. **b** Constant energy ARPES image obtained by integrating the spectral weight within *E*~F~ ± 10 meV recorded with 65 eV photons. Cuts 1 and 2 indicate the momentum locations of the experimental band structures in **d**--**i**. Red lines represent the (001) surface BZs. Red arrows indicate the triangle-shaped FSs centered at $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline K$$\end{document}$ points in **c**. The intensity within the cyan dashed rectangle has been multiplied by 10. **c** Projections of calculated bulk FSs on the (001) surface. White lines represent the (001) surface BZs. **d**, **e** ARPES intensity plot and corresponding momentum distribution curves along the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline K - \overline {\mathrm{\Gamma }} - \overline K$$\end{document}$ direction \[cut 1 in **b**\], respectively. Red and black dashed curves illuminate the electron bands forming the triangle-shaped FSs centered at $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline K$$\end{document}$ points. **f** Projections of calculated bulk bands along the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline K - \overline {\mathrm{\Gamma }} - \overline K$$\end{document}$ direction. **g** ARPES intensity plot along the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline M - \overline {\mathrm{\Gamma }} - \overline M$$\end{document}$ direction \[cut 2 in **b**\]. Black curves represent the calculated bulk bands along the four directions indicated by blue lines in **a**, where the bands from bottom to top correspond to the cuts from top to bottom, respectively. **h**, **i** Second derivative intensity plots along cuts 3 and 4 recorded at *hν* = 135 and 128 eV (*k*~z~ = ± 0.1455 *π*), respectively, with corresponding bulk calculations (superimposed red curves). The band crossings for WPs are highlighted by black curves. The crossing-*E*~F~ features of the *α*~1~ bands are marked out by white dashed ellipses

The detailed near-*E*~F~ band structures along the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline {\mathrm{\Gamma }} - \overline K$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline {\mathrm{\Gamma }} - \overline M$$\end{document}$ directions are presented in Fig. [3d--g](#Fig3){ref-type="fig"}. We first focus on the results along the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline K - \overline {\mathrm{\Gamma }} - \overline K$$\end{document}$ direction (Fig. [3d--f](#Fig3){ref-type="fig"}), the projections of calculated bulk bands along this direction in Fig. [3f](#Fig3){ref-type="fig"} well reproduce the experimental band features in Fig. [3d, e](#Fig3){ref-type="fig"}, i.e., the hole bands at near half of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline {{\mathrm{\Gamma }}K}$$\end{document}$ corresponding to the ring-like FSs and the electron bands (dashed curves in Fig. [3d, e](#Fig3){ref-type="fig"}) at $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline K$$\end{document}$ points for the triangle-like FSs, whose intensities are relatively weak due to the matrix element effect. The band structures along the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline M - \overline {\mathrm{\Gamma }} - \overline M$$\end{document}$ direction are displayed in Fig. [3g](#Fig3){ref-type="fig"}. One can obtain that the observed band dispersions are consistent with calculated bulk bands close to the *k*~z~ \~ *π* plane, as illustrated by black curves in Fig. [3g](#Fig3){ref-type="fig"}. Indicated by the four blue cuts in Fig. [3a](#Fig3){ref-type="fig"}, these bulk bands are located in the *k*~*x*~−*k*~*z*~ plane with *k*~*y*~ = 0 and *k*~*z*~ = 0.96, 0.92, 0.88, and 0.84 *π* from bottom to top, respectively, which are parallel to the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline {\mathrm{\Gamma }} - \overline M$$\end{document}$ direction. This consistency demonstrates the electronic states at high-symmetry *k*~*z*~ planes (*k*~z~ = *π* here) have main contributions to the ARPES spectra in the frame of *k*~*z*~ broadening effect^[@CR35]--[@CR37]^.

Based on the first-principles calculations, the WPs are not located on the high-symmetry *k*~*z*~ planes, they are embedded in the *k*~*z*~ = ±0.1455 *π* planes instead. By performing detailed photon-energy-dependent measurements in Supplementary Fig. [7](#MOESM1){ref-type="media"}, we can determine the specific *k*~*z*~ positions of predicted WPs. As shown in Supplementary Fig. [7a](#MOESM1){ref-type="media"}, the band dispersions along the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline {\mathrm{\Gamma }} - \overline K$$\end{document}$ direction show periodic modulation when varying the photon energy. The inner potential (*V*~0~) is thus estimated to be 10 eV. We illustrate the *k*~*z*~ = ±0.1455 *π* planes, where the Weyl fermions with opposite chirality are located, respectively, as blue curves in Supplementary Fig. [7b](#MOESM1){ref-type="media"}. To demonstrate the existence of WPs without any interference from possible trivial surface states, we utilize the photons with higher energies to record the band structure. The energies of 138 and 125 eV correspond to the *k*~*z*~ = ±0.1455 *π* planes in the neighbor BZ to that in Supplementary Fig. [7b](#MOESM1){ref-type="media"}, respectively. The measured second derivative intensity plots along cuts 3 and 4 (Fig. [3a](#Fig3){ref-type="fig"}), which are parallel to the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline {\mathrm{\Gamma }} - \overline M$$\end{document}$ direction in the *k*~z~ = ±0.1455 *π* planes, are displayed in Fig. [3h, i](#Fig3){ref-type="fig"}, respectively. The overlaid bulk band calculations well reproduce the experimental results, especially the *α*~1~ and *α*~2~ bands forming the WPs. The intensity asymmetry of the *α*~1~ band in opposite *k*~*x*~ regions derives from the matrix element effect. Since the cone nature of the Weyl dispersion is an important character proving the realization of massless Weyl fermions, we perform a FS mapping with the photon energy of 135 eV, and study the evolution of band structures gradually off the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline {\mathrm{\Gamma }} - \overline M$$\end{document}$ direction. As shown in Supplementary Fig. [8b, c](#MOESM1){ref-type="media"}, when measuring off the *k*~*y*~ = 0 line, the slopes of the near-*E*~F~ band dispersions (*α*~1~ and *α*~2~) forming the Weyl nodes decrease monotonically. This evolution is in accord with a cone nature within energy--momentum space. We would like to stress that, although the predicted WPs are located at \~50 meV above *E*~F~, the good agreement between theory and experiment, including the two specific *k*~z~ positions of WPs and the band structures forming the WPs below *E*~F~, combined with the valid Weyl--cone structure in experiment, could plausibly demonstrate the existence of Weyl fermions in Co~3~Sn~2~S~2~.

AHE of Co~3~Sn~2~S~2~ {#Sec6}
---------------------

Next, we move to the AHE of Co~3~Sn~2~S~2~. As shown in Fig. [4a](#Fig4){ref-type="fig"}, the Hall resistivity *ρ*~*xy*~ of Co~3~Sn~2~S~2~ at high temperature shows the linear field dependence. When *T* is close to *T*~C~, the *ρ*~*xy*~(*μ*~0~*H*) curve starts to bend significantly at low field region. For $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$T \ll T_{\mathrm{C}}$$\end{document}$, the *ρ*~*xy*~(*μ*~0~*H*) curve increases steeply at low field and then becomes weakly field-dependent at high field (above \~0.1 T). The field dependence of *ρ*~*xy*~(*μ*~0~*H*) resembles those of *M*(*μ*~0~*H*) (Fig. [4b](#Fig4){ref-type="fig"}), typical ones for ferromagnets. But the saturation values follow opposed temperature dependencies at $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$T \ll T_{\mathrm{C}}$$\end{document}$. When the saturation magnetization increases with decreasing temperature, the saturation value of the *ρ*~*xy*~(*μ*~0~*H*) decreases below *T*~C~. The maximum saturated value of *ρ*~*xy*~(*μ*~0~*H*) is \~21 μΩ cm at *T* = 140 K. We also measured the *ρ*~*xy*~(*μ*~0~*H*) of several crystals at the same temperature. They exhibit similar values, spanning between 19.59 and 24.22 μΩ cm at *μ*~0~*H* = 5 T (Supplementary Fig. [9](#MOESM1){ref-type="media"}). The standard deviation of determined sample thickness is also below 5% (Supplementary Fig. [10](#MOESM1){ref-type="media"}). Moreover, the *ρ*~*xy*~(*μ*~0~*H*) at various field directions are also measured. When the magnetic field is titled away from the *c* axis of Co~3~Sn~2~S~2~ crystal up to ±30°, the change of *ρ*~*xy*~(*μ*~0~*H*) is small and the \[*ρ*~*xy*~(*μ*~0~*H*)(30°) − *ρ*~*xy*~(0°)\]/*ρ*~*xy*~(0°) at 5 T is only about 1% (Supplementary Fig. [11](#MOESM1){ref-type="media"}). All of these results clearly indicate that the experimental *ρ*~*xy*~(*μ*~0~*H*) is reproducible and reliable.Fig. 4AHE of Co~3~Sn~2~S~2~. **a**, **b** Hall resistivity *ρ*~*xy*~(*μ*~0~*H*) and *M*(*μ*~0~*H*) as a function of magnetic field at various temperatures for *H*\|\|*c*, respectively. **c** Plot of log $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{{A}}(T)$$\end{document}$ vs. log *ρ*~*xx*~(*T*). The red solid line is the fit using the relation $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}} \propto \rho _{{{xx}}}(T)^\alpha$$\end{document}$. **d** Plot of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}}{\mathrm{/}}\left( {M\rho _{{{xx}}}} \right)$$\end{document}$ vs. *ρ*~*xx*~. The green solid line is the linear fit when *ρ*~*xx*~ ≤ 250 μΩ cm. Insets of **c** and **d** show the temperature dependence of *R*~0~(*T*) and *R*~s~(*T*) at *T* ≤ *T*~C~, respectively. **e** Intrinsic AHC $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{xy,\mathrm{in}}}}^{\mathrm{A}}(T)$$\end{document}$ as a function of *M*. **f** AHC $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{xy}}}^{\mathrm{A}}(T)$$\end{document}$ as a function of temperature. Inset: temperature dependence of scale factor *S*~H~(*T*). Error bars in **c**--**f** are the standard deviations of the fits when determining the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}}$$\end{document}$ (Supplementary Note [1](#MOESM1){ref-type="media"})

It is known that the Hall resistivity *ρ*~*xy*~ in ferromagnets arises from two parts^[@CR4]^, *ρ*~*xy*~ = $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{O}} + \rho _{{{xy}}}^{\mathrm{A}}$$\end{document}$ = *R*~0~*B* + 4*πR*~s~*M*, where $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{O}}$$\end{document}$ is the normal Hall resistivity due to the Lorentz force, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}}$$\end{document}$ is the anomalous Hall resistivity, *R*~0~ is the ordinary Hall coefficient, and *R*~s~ is the anomalous Hall coefficient. From the *ρ*~*xy*~(*μ*~0~*H*) and *M*(*μ*~0~*H*) curves (Supplementary Fig. [12](#MOESM1){ref-type="media"} and Supplementary Note [1](#MOESM1){ref-type="media"}), the *R*~0~ and *R*~s~ can be determined. The *R*~0~ is positive at entire  temperature region (inset of Fig. [4c](#Fig4){ref-type="fig"}), indicating that the dominant carrier is hole-type. The apparent charge carrier density *n*~a~ can be deduced using the relation of *n*~a~ \~ −1/\|*e*\|*R*~0~ (Supplementary Fig. [13](#MOESM1){ref-type="media"}), and it reaches 2.1 × 10^22^ cm^−3^ at 5 K, corresponding to about 0.7 carriers per formula unit of Co~3~Sn~2~S~2~. The *R*~s~(*T*) is also positive but the absolute value is much larger than *R*~0~ (inset of Fig. [4d](#Fig4){ref-type="fig"}). It increases monotonically with increasing temperature. The obtained *R*~s~ value is 2.7 × 10^−9^ Ω cm G^−1^ at 170 K, which exhibits an enhancement of two orders of magnitude when compared with the conventional itinerant ferromagnets, such as pure Fe and Ni^[@CR38],[@CR39]^.

Discussion {#Sec7}
==========

The relation between log $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}}$$\end{document}$ and log *ρ*~*xx*~ at low temperature region (Fig. [4c](#Fig4){ref-type="fig"}) can be fitted using the formula $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}} \propto \rho _{{{xx}}}^\alpha$$\end{document}$, which gives the scaling exponent *α* = 1.87(4). It suggests that the intrinsic KL or extrinsic side-jump mechanism are mainly responsible for the AHE in Co~3~Sn~2~S~2~. On the other hand, the dominant mechanism of AHE can also be decided by fitting the relationship between $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}}$$\end{document}$ and *ρ*~*xx*~ using the formula $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}}$$\end{document}$ = $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a(M)\rho _{{{xx}}} + b(M)\rho _{{{xx}}}^{\mathrm{2}}$$\end{document}$. The first term corresponds to the skew scattering contribution, while the second term represents the intrinsic or side-jump contribution^[@CR40]^. For the skew scattering contribution, the *a*(*M*) is usually proportional to *M* linearly^[@CR41]^, and for the intrinsic contribution, the *b*(*M*) = $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}}{\mathrm{/}}\rho _{{{xx}}}^{\mathrm{2}}$$\end{document}$ is directly related to the intrinsic AHC $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{xy,\mathrm{in}}}}^{\mathrm{A}} = - \rho _{{{xy}}}^{\mathrm{A}}{\mathrm{/}}\rho _{{{xx}}}^{\mathrm{2}}$$\end{document}$. Previous study suggests that the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{xy,\mathrm{in}}}}^{\mathrm{A}}$$\end{document}$ is also proportional to *M* linearly^[@CR40]^, thus the linear fit of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}}{\mathrm{/}}\left( {M\rho _{{{xx}}}} \right)$$\end{document}$ vs. *ρ*~*xx*~ can separate the intrinsic and extrinsic contributions. The relation between $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}}{\mathrm{/}}\left( {M\rho _{{{xx}}}} \right)$$\end{document}$ and *ρ*~*xx*~ (Fig. [4d](#Fig4){ref-type="fig"}) shows the linear dependence when *ρ*~*xx*~ ≤ 250 μΩ cm. After subtracting the skew-scattering contribution, the obtained value of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$| {\sigma _{{{xy,\mathrm{in}}}}^{\mathrm{A}}} |$$\end{document}$ at 5 K is about 505 Ω^−1^ cm^−1^, which is comparable with the predicted value by using the formula of *σ*~*xy*~ = *Ke*^2^/4*π*^2^ (\~525 Ω^−1^ cm^−1^) (Supplementary Fig. [6](#MOESM1){ref-type="media"}) but smaller than that calculated from the integral of Berry curvature along *k*~*z*~ in the BZ (\~1310 Ω^−1^ cm^−1^) when setting the moment of Co as 0.33 μ~B~ (Fig. [2e](#Fig2){ref-type="fig"}). The discrepancy between experimental results and the theoretical results from the latter method is still not clear. But such kind of difference has also been observed in other systems, such as Fe \[751 Ω^−1^ cm^−1^ (experiment) vs. 1032 Ω^−1^ cm^−1^ (theory)\], and Ni \[−2073 Ω^−1^ cm^−1^ (experiment) vs. −646 Ω^−1^ cm^−1^ (theory)\]^[@CR14],[@CR42]^. On the other hand, the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{xy,\mathrm{in}}}}^{\mathrm{A}}$$\end{document}$ exhibits the linear dependence of *M* when *T* \< *T*~C~ (Fig. [4e](#Fig4){ref-type="fig"}), well consistent with the KL theory^[@CR10]^. Importantly, this result is also well in agreement with the theoretical predictions that the smaller moment of Co leads to the nearly linear decrease of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{xy}}}^{\mathrm{A}}$$\end{document}$ (Fig. [2e](#Fig2){ref-type="fig"} and Supplementary Fig. [6](#MOESM1){ref-type="media"}). In contrast to the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$| {\sigma _{{{xy,\mathrm{in}}}}^{\mathrm{A}}} |$$\end{document}$, the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$| {\sigma _{{{xy,\mathrm{sj}}}}^{\mathrm{A}}} |$$\end{document}$ for the side-jump contribution can be estimated using an expression (*e*^2^/(*ha*)(*ε*~SO~/*E*~F~), where *ε*~SO~ is the spin--orbit interaction energy^[@CR41],[@CR43]^. Using the lattice constant *a* \~ *V*^1/3^ = 9.96 Å and *ε*~SO~/*E*~F~ \~ 0.01 for the metallic ferromagnets, the derived $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{xy,\mathrm{sj}}}}^{\mathrm{A}}$$\end{document}$ is only about 3.9 Ω^−1^ cm^−1^, thus the extrinsic side-jump contribution to $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{xy}}}^{\mathrm{A}}$$\end{document}$ could be very small when compared to $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$| {\sigma _{{{xy,\mathrm{in}}}}^{\mathrm{A}}} |$$\end{document}$. To further investigate the mechanism of AHE in Co~3~Sn~2~S~2~, the temperature dependence of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{xy}}}^{\mathrm{A}}$$\end{document}$ is shown in Fig. [4f](#Fig4){ref-type="fig"}. The $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$| {\sigma _{{{xy}}}^{\mathrm{A}}} |$$\end{document}$ shows fairly large values at low temperature region and it is insensitive to the change of temperatures. Above 140 K, however, it decreases quickly with increasing temperature and reaches to about 150 Ω^−1^ cm^−1^ at 175 K. In contrast, the scale factor *S*~H~ ($\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$= {\mathrm{\mu}} _0R_{\mathrm{s}}{\mathrm{/}}\rho _{{{xx}}}^{\mathrm{2}} = - \sigma _{{{xy}}}^{\mathrm{A}}{\mathrm{/}}M$$\end{document}$) is almost constant when *T* \< *T*~C~ (inset of Fig. [4f](#Fig4){ref-type="fig"}). Therefore it further confirms that the temperature dependence of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{x}}y}^{\mathrm{A}}$$\end{document}$ originates from the temperature dependence of *M*, i.e., the dramatic decrease of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$| {\sigma _{{{xy}}}^{\mathrm{A}}} |$$\end{document}$ above 140 K is mainly due to the sharp decrease of *M*(*T*) as the temperature approaches *T*~C~. Above results clearly indicate that the large intrinsic AHC is dominant in HMFM Co~3~Sn~2~S~2~, and it is closely related to the nontrivial topology of band structure with Weyl nodes near *E*~F~.

In summary, we have investigated the AHE and electronic structure of HMFM Co~3~Sn~2~S~2~ single crystal. The nearly quadratic relationship between $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\rho _{{{xy}}}^{\mathrm{A}}$$\end{document}$ and *ρ*~*xx*~ indicates that the mechanism of AHE in Co~3~Sn~2~S~2~ is dominated by the intrinsic contribution. The consistency between the experimental band structures and first-principles calculations, especially the spindle-shaped FSs along the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\overline {\mathrm{\Gamma }} - \overline M$$\end{document}$ direction and the band dispersions forming the WPs, combined with the Weyl-cone structure presented in ARPES, indicate that the intrinsic AHE in Co~3~Sn~2~S~2~ originates from the existence of magnetic Weyl fermions near *E*~F~. Moreover, the steep decrease of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{xy}}}^{\mathrm{A}}$$\end{document}$ at high temperature can be mainly ascribed to the sharp decrease of magnetic moment, causing the change of topological properties in band structure. Current work will not only deepen our understanding on the exotic physical phenomena associated with nontrivial band topology, but also shed light on exploring novel electronic/spintronic devices based on AHE and/or half metallicity.

Methods {#Sec8}
=======

Single crystal growth and structural characterization {#Sec9}
-----------------------------------------------------

Co~3~Sn~2~S~2~ single crystals were grown by the Sn flux^[@CR44]^. The starting materials Co lumps (purity 99.99%), Sn grains (purity 99.99%) and S flakes (purity 99.999%) were mixed together in a molar ratio of Co:Sn:S = 8:86:6 and placed in an alumina crucible. Then, the crucible was sealed in the quartz ampoule under partial argon atmosphere. The sealed quartz ampoule was initially heated at 673 K for 2 h in order to avoid bursting of ampoule due to the high vapor pressure of S. And then it was heated up to 1323 K and kept there for 6 h to ensure the homogeneity of melt. Finally the ampoule was cooled down slowly to 973 K with 5 K/h. At this temperature, the ampoule was taken out from the furnace rapidly and decanted with a centrifuge to extract Co~3~Sn~2~S~2~ single crystals from the Sn flux. X-ray diffraction (XRD) patterns were measured using a Bruker D8 X-ray machine with Cu*K*~*α*~ radiation (*λ* = 0.15418 nm) at room temperature.

Magnetization and transport measurements {#Sec10}
----------------------------------------

Magnetization and electrical transport measurements were performed in a Quantum design MPMS3 and PPMS, respectively. Both longitudinal and Hall electrical resistivity were measured using a standard four-probe method on rectangular shape single crystals with current flowing in the *ab* plane. In order to effectively avoid the longitudinal resistivity contribution due to voltage probe misalignment, the Hall resistivity was measured by sweeping the field from −5 T to 5 T at various temperatures, and the total Hall resistivity was determined by *ρ*~*xy*~(*μ*~0~*H*) = \[*ρ*(+*μ*~0~*H*) − *ρ*(−*μ*~0~*H*)\]/2, where *ρ*~*xy*~(±*μ*~0~*H*) is the transverse resistivity under a positive or negative magnetic field.

ARPES measurements {#Sec11}
------------------

ARPES measurements were performed at the Dreamline beamline of the Shanghai Synchrotron Radiation Facility (SSRF) with a Scienta D80 analyzer and at the 1-cubed ARPES end-station of BESSY using a Scienta R4000 analyzer. The energy and angular resolutions were set to 15 meV and 0.2°, respectively. All samples were cleaved in situ along the (001) plane and measured at *T* = 20 K in a working vacuum better than 5 × 10^−11^ Torr.

Theoretical calculation {#Sec12}
-----------------------

We have simulated the electronic structure of Co~3~Sn~2~S~2~ with first-principle calculations using Vienna ab initio Simulation Package (VASP)^[@CR45]^, and the generalized gradient approximation (GGA) of Perdew--Burke--Ernzerhof (PBE) type exchange correlation potential^[@CR46]^ was employed. The spin--orbit coupling was taken into account in all of the calculations. We have relaxed the crystal structure and the lattice constants *a* = *b* = 5.3892 Å, and *c* = 13.1519 Å are adopted in our calculation. The surface states and anomalous Hall conductivity were calculated based on the tight binding Hamiltonian constructed by using Wannier90 package^[@CR47]^. The maximally localized Wannier functions for 3*d* orbitals on Co, 5*p* orbitals on Sn, and 3*p* orbitals on S have been used as the basis of the tight-binding Hamiltonian. We calculate the anomalous Hall conductivity (AHC) as the sum of Berry curvatures over all of the occupied bands^[@CR48]^,$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma _{{{xy}}} = - \frac{{2\pi e^2}}{h}{\int}_{{\mathrm{BZ}}} \frac{{\mathrm{d}^3{\bf{k}}}}{{(2\pi )^3}}\mathop {\sum}\limits_n {\kern 1pt} f_n({\bf{k}}){\mathrm{\Omega }}_n^{\mathrm{z}}({\bf{k}})$$\end{document}$$where *f*~*n*~(**k**) is the Fermi--Dirac distribution function, and *n* is the index of the occupied bands. The Berry curvature can be arisen from the Kubo-formula derivation,$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\mathrm{\Omega }}_n^{{z}}({\bf{k}}) = - 2{\mathrm{Im}}\mathop {\sum}\limits_{m \ne n} \frac{{\left\langle {{\mathrm{\Psi }}_{n{\bf{k}}}} \right|v_{{x}}\left| {{\mathrm{\Psi }}_{m{\bf{k}}}} \right\rangle \left\langle {{\mathrm{\Psi }}_{m{\bf{k}}}} \right|v_{{y}}\left| {{\mathrm{\Psi }}_{n{\bf{k}}}} \right\rangle }}{{\left( {E_m({\bf{k}}) - E_n({\bf{k}})} \right)^2}}$$\end{document}$$where *v*~*x*(*y*)~ is the velocity operator. The intrinsic AHC is calculated with 200 × 200 × 200 *k*-point grid based on the tight binding Hamiltonian.

Electronic supplementary material
=================================

 {#Sec14}

Supplemental Information

**Publisher\'s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Qi Wang, Yuanfeng Xu, Rui Lou.

**Change history**

10/8/2018

The original version of this Article incorrectly omitted an affiliation of Hongming Weng: "Beijing National Laboratory for Condensed Matter Physics and Institute of Physics, Chinese Academy of Sciences, 100190, Beijing, China"

This has been corrected in both the PDF and HTML versions of the Article.

Electronic supplementary material
=================================

**Supplementary Information** accompanies this paper at 10.1038/s41467-018-06088-2.

This work was supported by the National Key R&D Program of China (Grant nos. 2016YFA0300504, 2016YFA0300600), the National Natural Science Foundation of China (Nos. 11574394, 11774423, 11774421, and 11674369), the Fundamental Research Funds for the Central Universities, and the Research Funds of Renmin University of China (RUC) (15XNLF06, 15XNLQ07, 17XNH055, 18XNLG14, 18XNH057). Z.L. acknowledges Shanghai Sailing Program (17YF1422900). Y.X and H.W are also supported by the "Strategic Priority Research Program (B)" of the Chinese Academy of Sciences (Grant no. XDB07020100).

H.L., H.W., and S.W. provided strategy and advice for the research; Q.W. performed the crystal growth, magnetization and transport measurements and fundamental data analysis; R.L., Z.L., M.L., Y.H. and D.S. performed the ARPES measurements and analyzed the data; Y.X. and H.W. carried out the theoretical calculation; H.L. wrote the manuscript based on discussion with all the authors.

Data measured or analyzed during this study are available from the corresponding authors on reasonable request.

Competing interests {#FPar1}
===================

The authors declare no competing interests.
