Large anomalous Hall effect in the chiral-lattice antiferromagnet CoNb3S6.

An ordinary Hall effect in a conductor arises due to the Lorentz force acting on the charge carriers. In ferromagnets, an additional contribution to the Hall effect, the anomalous Hall effect (AHE), appears proportional to the magnetization. While the AHE is not seen in a collinear antiferromagnet, with zero net magnetization, recently it has been shown that an intrinsic AHE can be non-zero in non-collinear antiferromagnets as well as in topological materials hosting Weyl nodes near the Fermi energy. Here we report a large anomalous Hall effect with Hall conductivity of 27 Ω-1 cm-1 in a chiral-lattice antiferromagnet, CoNb3S6 consisting of a small intrinsic ferromagnetic component (≈0.0013 μB per Co) along c-axis. This small moment alone cannot explain the observed size of the AHE. We attribute the AHE to either formation of a complex magnetic texture or the combined effect of the small intrinsic moment on the electronic band structure.

Introduction

The anomalous Hall effect (AHE) is a signature of emergent electromagnetic fields in solids that affect the motion of the electrons, and hence it has been a recent topic of intense study in the context of the topological materials[1,2]. The Hall effect in general is an intrinsic property of a conductor due to the Lorentz force experienced by the charge carriers. In systems with spontaneously broken time-reversal symmetry, an additional contribution, independent of the Lorentz force, is observed, the AHE[1]. AHE was first observed in ferromagnets where its origin lies in the interplay between spin–orbit coupling (SOC) and magnetization. Reformulation of the SOC-induced intrinsic mechanism of AHE in ferromagnets to the Berry phase curvature in momentum space has been fruitful in predicting and describing the AHE in several other systems, including Weyl (semi)metals[3], non-collinear antiferromagnets[4], non-coplanar magnets[5-7], and other nontrivial spin textures[8-11]. Recent observations of the large anomalous Hall effect in metals with possible Weyl[12-14] and massive Dirac fermions[15,16] and/or complex spin textures, e.g., skyrmion bubbles[17], have generated interest in such materials, especially for the role of correlated topological states in the emergent electronic properties. Here we present a large AHE in CoNb3S6 that cannot be understood in terms of conventional mechanisms of the AHE.

CoNb3S6 is a member of a large class of intercalated transition metal dichalcogenides, where a 3d-transition metal sandwiches between layers of a 5d-transition metal dichalcogenide that are coupled by a weak Van der Waals force. CoNb3S6 represents the 1/3 fractional intercalation of Co atoms between the layers of NbS2[18]. It crystallizes in the hexagonal chiral space group P63 22[19,20]. It orders antiferromagnetically below ~26 K and is known to have a rather complex susceptibility for the magnetic field applied parallel to the c-axis[21,22]. At 4 K, neutron diffraction measurement has suggested a collinear antiferromagnetic state in which the spins are in the ab-plane pointing along a certain crystallographic axis[20]. By itself, however, such a spin structure cannot give rise to the anomalous Hall effect.

The key finding of this work is a large c-axis anomalous Hall effect in the antiferromagnetic CoNb3S6. Although CoNb3S6 shows a small, intrinsic ferromagnetic component (≈0.0013 μB per Co) along the c-axis, this small moment alone cannot explain the observed size of the AHE. Based on its chiral crystal structure and the calculated band structure, we attribute the AHE in CoNb3S6 either to the formation of a complex magnetic texture or to an influence of the small intrinsic ferromagnetic moment on the underlying electronic band structure.

Results

Crystal structure and physical properties of CoNb3S6

We verified the room temperature crystal structure of CoNb3S6 in the chiral space group P6322 by means of single crystal X-ray diffraction (see Supplementary Note 1 and Supplementary Fig. 1). A sketch of the crystal structure adopted by CoNb3S6 is shown in Fig. 1a, where the Co atoms occupy the octahedral position between the triangular prismatic layers of the parent compound 2H-NbS2. Figure 1b shows the magnetic susceptibility of CoNb3S6 as a function of temperature measured with a magnetic field of 0.1 T applied along a- and c-axis (χ and χ, respectively). Over the entire temperature range, χ exceeds χ, consistent with behavior reported in the literature[21]. χ shows a kink at TN = 27.5 K corresponding to the antiferromagnetic transition and decreases on further cooling down to 1.8 K. The field cooled (FC) and zero field cooled (ZFC) measurements show identical behavior. Along the c-axis, the susceptibility below TN shows irreversibility between FC and ZFC measurements. The ZFC χ shows a behavior similar to that of χ. However, the FC χ increases on cooling below TN, becomes maximum at 25 K and decreases on further cooling. This increase in FC susceptibility along the c-axis near the transition temperature has been reported in previous studies[20,21] and implies the presence of a small ferromagnetic component along the c-axis. With increasing magnetic field, the magnitude of this jump decreases, and at 3 T, FC susceptibility shows a continuous decrease below TN, as shown in the inset of Fig. 1b, suggesting that the moment along the c-axis is suppressed by a sufficiently large magnetic field. At higher temperatures, the susceptibility shows Curie–Weiss behavior. A Curie–Weiss fit to the data between 200 and 300 K (Supplementary Note 2 and Supplementary Fig. 2) yields the powder averaged effective moment of 3.0 μB per Co. This value is smaller than the spin-only moment expected for Co2+ (3.87 μB) established by neutron diffraction and optical studies[20,23]. Resistivity shows metallic behavior over the measured temperature range of 1.8–300 K, with a sudden drop below 27.5 K (TN), presumably due to the reduction of electron scattering with the onset of the spin ordering (Fig. 1c). CoNb3S6 has a very small magnetoresistance of 0.2% at 2 K in the field of 9 T applied along c-axis, as shown in the inset of Fig. 1c.

Fig. 1

Crystal structure and characteristics of CoNb3S6. a Sketch of crystal structure of CoNb3S6. Each Co atom is octahedrally coordinated with S atoms. b Magnetic susceptibility measured with magnetic field along a- and c-axis. c Electrical resistivity as a function of temperature measured with current along a-axis. Inset shows the magnetoresistance defined by MR = (ρ − ρ0)/ρ0 × 100%, where ρ is the resistivity measured at the magnetic field H and ρ0 the resistivity measured at H = 0. The MR was measured with current along a-axis and magnetic field along c-axis

Magnetization

M vs H measured between 29 and 22 K, with the magnetic field along c-axis, is shown in Fig. 2a, and the corresponding first-order derivatives are shown in Fig. 2b. Despite the nearly linear M–H curves, a hysteresis with very small remanent magnetization appears below TN, which becomes maximum at 25 K (0.0013 μB per formula unit) and decreases at lower temperatures, as shown in Fig. 2b, c. The hysteresis becomes more clear in dM/dH vs H plots shown in Fig. 2b. At 29 K, which is above TN, dM/dH is featureless. At 27 K, when decreasing the field from 6 T to −6 T, a peak in dM/dH appears at −1 T. When increasing the field from −6 T to 6 T, no peak is seen at −1 T. Instead, the peak is observed at 1 T. These peaks appear at the coercive field of the hysteresis in M vs H. With the decrease in temperature, the coercive field increases and becomes 3 T at 24 K. No such peaks are seen below 24 K. These data suggest that there is an intrinsic, hard ferromagnetic component along the c-axis in the ordered state that is flipped by a temperature-dependent critical field. Switching of this ferromagnetic component gives rise to the hysteresis in M vs H. Below 24 K, a magnetic field of 6 T cannot switch the ferromagnetic component. As a result, no hysteresis is observed. Application of a larger magnetic field is expected to reveal the hysteresis even below 24 K. In fact, this hysteresis is evident in the Hall effect measurement presented below. In Fig. 2c we show the ferromagnetic component as a function of the magnetic field obtained by subtracting a linear antiferromagnetic background from the M vs H curves presented in Fig. 2a. On the other hand, when magnetic field is applied along the a-axis, M vs H shows linear behavior at all temperatures measured, as depicted in Fig. 2d–f.

Fig. 2

Magnetization of CoNb3S6. a M vs H measured with magnetic field applied along the c-axis. b dM/dH vs H for magnetic field applied along the c-axis. The black (red) curves represent dM/dH for decreasing (increasing) H. c Ferromagnetic moment (ΔM) along the c-axis obtained by subtracting the antiferromagnetic background from the measured M vs H presented in a. d M vs H for magnetic field applied along the a-axis. e dM/dH vs H for magnetic field applied along the a-axis. The black (red) curves represent dM/dH for decreasing (increasing) H. f Low field M vs H measured with magnetic field along a-axis. dM/dH plots in b and e are shifted by a constant factor for clarity

Hall effect

Hall resistivity (ρ) vs H measured between 28 and 23 K with current along a-axis and magnetic field along c-axis is depicted in Fig. 3a. The Hall effect at each temperature was measured by cooling the sample in a magnetic field of 9 T and then by changing the magnetic field from 9 to −9 T and subsequently from −9 to 9 T. At 28 K, which is just above TN, the Hall resistivity is linear, as expected for a normal conductor. The sign of the Hall resistivity indicates that holes are the majority charge carriers in CoNb3S6. Within the single band model, the estimated carrier concentration (n = 1/|eR0|) is 2.49 × 1021 cm−3, where e is the charge of an electron and R0 is the ordinary Hall coefficient. When the temperature is lowered below TN, an additional contribution appears in the Hall resistivity. At 27 K, this anomalous Hall signal is small. It increases both in magnitude and in field with decreasing temperature, becoming maximum at 23 K. Below 23 K, no switching behavior of the Hall resistivity is observed when measured with the maximum field of 9 T. However, the Hall resistivity increases down to 20 K (Supplementary Note 3 and Supplementary Fig. 3). At 15 K, a very small AHE is observed, and at 2 K there is no anomalous Hall effect. These observations reveal two important points. (1) Hysteresis in the AHE is observed when the magnetic field switches the small FM component (Fig. 2a–c). (2) The AHE is observed due to the stabilization of the FM component, which is observed between TN and 23 K. At 20 K, a magnetic field of 9 T cannot switch the FM component, and hence the AHE does not switch sign, but gives the same large value. At 2 K, a magnetic field of 9 T cannot induce the FM component at all. As a result, there is no AHE observed for |H| ≤ 9 T.

Fig. 3

Hall effect in CoNb3S6 with magnetic field along c-axis. a Hall resistivity as a function of magnetic field. b Anomalous Hall resistivity as a function of magnetic field obtained by subtracting the ordinary Hall resistivity ρ as discussed in the text. c Scaling of the Hall resistivity and the ferromagnetic moment. d Anomalous Hall resistivity as a function of magnetization. The straight line represents the slope of the high field vs M that gives the anomalous Hall part due to the ferromagnetic component from the field-induced spin tilt. e Anomalous Hall resistivity obtained by subtracting the estimated Hall contribution form the ferromagnetic component of magnetic field-induced spin tilt as a function of magnetic field. f Hall conductivity as a function of magnetic field

Conventionally, the Hall resistivity in a ferromagnet is given by the relation , where R0 is the ordinary Hall coefficient, is the anomalous Hall resistivity[1]. In Fig. 3a, we see that the Hall resistivity measured in the paramagnetic state, i.e., at 28 K and in the antiferromagnetic state, except near the switching fields has the same slope. This indicates that R0B is constant between 28 and 23 K, and allows us to subtract the ordinary Hall resistivity to obtain the anomalous part. Figure 3b shows vs H, where is obtained by subtracting the data measured at 28 K from ρ measured below TN. As there is a finite ferromagnetic component along the c-axis, we first examine the effect of this small ferromagnetic component on the anomalous Hall effect. For a ferromagnet, , where RS is the anomalous Hall coefficient, μ0 is the permeability, and M is the magnetization. RS scales the M–H curve to the anomalous part of the Hall resistivity. To check for this scaling in CoNb3S6, we obtained the magnetization of the ferromagnetic component by subtracting a straight line from the M vs H curve measured at a representative temperature of 25 K, which shows the maximum remanent magnetization. The obtained magnetization multiplied by −1, (−ΔM) vs H is plotted in Fig. 3c. The anomalous Hall resistivity measured at 25 K is also plotted in the same figure. These plots show that the anomalous Hall effect does scale with magnetization. RS is negative and has an exceedingly large value compared to R0. RS estimated from this scaling at 25 K is −1.41 × 10−4 m3 C−1, which is four orders of magnitude larger than R0 = 2.4 × 10−9 m3 C−1. In a canonical ferromagnet-like Fe[24] or ferromagnets with closely related and even the same crystal structure as CoNb3S6, such as Fe1/4 TaS2[25] and Fe1/3TaS2[26], RS/R0 ≈ 100, two orders of magnitude smaller than that in CoNb3S6. This significant boost of the coefficient of AHE in CoNb3S6 is unlikely to be accounted for by the effect of the weak c-axis ferromagnetic component alone and therefore must have an additional contribution.

The contribution of the ferromagnetic component to the anomalous Hall resistivity, i.e., RS can be estimated from the high magnetic field part of the anomalous Hall resistivity vs magnetization plot. As , the slope of vs M gives μ0RS. In Fig. 3d, we see that vs M resembles vs H such that changes sign at the coercive fields. After the abrupt sign change, shows a linear variation with M. This linear behavior of vs M at high field arises from the continuous spin tilting by the magnetic field toward the c-axis and represents the contribution from the field-induced ferromagnetic component. The slope in this linear regime gives μ0RS for this compound, which is represented by dashed lines in Fig. 3d. RS estimated from this slope is −1.12 × 10−7 m3 C−1, which is comparable to that of the ferromagnet Fe1/3TaS2[26]. The extra component in the anomalous Hall resistivity, which we label as can be obtained by subtracting this linear part from the anomalous Hall effect: and is shown in Fig. 3e, which shows that contribution from the FM component from spin tilting is negligibly small compared to . In Fig. 3f, we plot the anomalous Hall conductivity as a function of magnetic field between 27 and 23 K. The anomalous Hall conductivity at 23 K takes the value 27 Ω−1 cm−1. A similar behavior, in which a very small ferromagnetic component and its anomalously large scaling with the anomalous Hall resistivity, has been observed in Mn3Sn and Mn3Ge. Here, the anomalous Hall effect has been attributed to the combined effect of non-collinear antiferromagnetic spin texture and electronic band structure[12,13]. A large anomalous Hall coefficient, reported for Fe3Sn2[15,16] cannot be described by the ferromagnetic component alone and is believed to be due to the massive quasi-2D Dirac cones near the Fermi energy. Likewise, the large AHE observed in frustrated systems such as Pr2Ir2O7 has been attributed to the non-coplanar spin texture[27].

Discussion

We now discuss possible origins of the anomalous Hall effect in CoNb3S6. First, we can rule out impurity-related extrinsic mechanisms of the AHE from the fact that the AHE vanishes at low temperatures, where the laboratory magnetic fields are insufficient to alter spins state, e.g., at 2 K as shown in Fig. 3a and Supplementary Fig. 3a. Second, a collinear antiferromagnetic state cannot give rise to the anomalous Hall effect. Finally, we showed above that an FM component due to simple magnetic field-induced spin tilt alone cannot account for the observed AHE. We thus consider other scenarios for the origin of the AHE in CoNb3S6: (1) non-collinear, non-coplanar or other complex magnetic textures, or (2) interplay of electronic and magnetic degrees of freedom.

CoNb3S6 has a chiral crystal structure. In a magnet with such a chiral lattice, competition among the exchange interaction, Dzyaloshinskii–Moriya interaction (DMI), magnetocrystalline anisotropy, and Zeeman energy may result in complex magnetic textures[28-30], including non-collinear or non-coplanar spins that can give rise to large AHE[10-13,27,31]. The neutron diffraction experiment reported by Parkin et al.[20] revealed no such complex spin structure. However, this is not unexpected given (1) the measurement is carried out in zero field at 4 K, where our data does not show AHE when measured up to the magnetic field of 9 T, and (2) a diffraction experiment such as that carried out by Parkin et al. would be insensitive to the large-scale spin structures. Thus, the nature of the spin structure in CoNb3S6 near TN remains an open question.

Figure 4a–c shows the paramagnetic band structure of CoNb3S6 without inclusion of SOC. The metallic character of CoNb3S6 is evident by the presence of the hole pockets along Γ−A, which is consistent with the measured hole-like character of the charge carriers. There are two linear electron and hole band crossings along Γ−M and K−Γ. Inclusion of SOC (Fig. 4b) opens a small gap at the band crossings along these lines and splits the bands due to the lack of inversion symmetry[32,33]. The band structure along Γ−A is enlightening in this regard. In the calculation with SOC (Fig. 4a) both at Γ and A, two doubly degenerate bands touch at a point ≈40 meV above the Fermi energy (EF). The splitting of these bands due to SOC along Γ−A (Fig. 4b) is remarkably large and results in the formation of several linearly crossing bands and avoided crossings. In Fig. 4d, we show the band dispersion along A–Γ–A near EF (see Supplementary Note 4 and Supplementary Fig. 4 for symmetry analysis). Bands cross linearly at all of these points. Two of these points lie within ≈15 meV of EF. It has been established rigorously that in non-magnetic crystals and in directions parallel to a 63 screw axis (along Γ−A in CoNb3S6), band degeneracies at the Γ points are Weyl nodes[34]. The AHE in CoNb3S6 is observed only when the external magnetic field induces a small moment along c-axis, which suggests that the large enhancement in AHE may result due to the combined effect of such a field in the presence of the near EF Weyl nodes.

Fig. 4

Electronic band structure of CoNb3S6. a Without inclusion of spin–orbit coupling and b with inclusion of spin–orbit coupling. Spin–orbit coupling causes splitting of the bands, which is most pronounced along Γ−A. c Brillouin zone of CoNb3S6. d Electronic band structure (with SOC) near the Fermi energy along the A–Γ–A high symmetry line

In summary, CoNb3S6 is a member of the class of antiferromagnets that show a large anomalous Hall effect. The large AHE is not explained by the reported collinear antiferromagnetic structure, and we suggest the formation of complex, noncollinear magnetic textures or an interplay between the magnetic texture and electronic band structure as two possible mechanisms for the large AHE. CoNb3S6 is a member of a large family of the 1/3 intercalated TX2 compounds, where T is a 4- or 5-d transition metal element, and X is a chalcogen (S, Se). As the change in the intercalated 3d element changes the nature of both magnetic and electronic structure, our realization of the AHE in CoNb3S6 opens a platform to explore and perhaps manipulate the interplay among spin texture, electronic band structure, and the associated emergent phenomena in a large class of poorly explored materials.

Methods

Crystal growth and characterization

Single crystals of CoNb3S6 were grown by chemical vapor transport using iodine as the transport agent. First, a polycrystalline sample was prepared by heating stoichiometric amounts of cobalt powder (Alfa Aesar 99.998%), niobium powder (Johnson Matthey Electronics 99.8%), and sulfur pieces (Alfa Aesar 99.9995%) in an evacuated silica ampoule at 900 °C for 5 days. Subsequently, 2 g of the powder was loaded together with 0.5 g of iodine in a fused silica tube of 14 mm inner diameter. The tube was evacuated and sealed under vacuum. The ampoule of 11 cm length was loaded in a horizontal tube furnace in which the temperature of the hot zone was kept at 950 °C and that of the cold zone was ≈850 °C for 7 days. Several CoNb3S6 crystals formed with a distinct, well-faceted flat plate-like morphology. The crystals of CoNb3S6 were examined by single crystal X-ray diffraction at beamline 15-ID-D at the APS, Argonne National Laboratory (ANL), where the data were collected with an APEX2 Area Detector using synchrotron radiation (λ = 0.41328 Å) at 293 K. Compositional analysis was done using an energy dispersive X-ray spectroscopy (EDS) at the Electron Microscopy Center, ANL.

Magnetic and transport property measurements

Magnetization measurements were made using a Quantum Design VSM SQUID. Both in FC and ZFC mode, susceptibility data were measured by sweeping temperature up from 1.8 K. In the derivative of M vs H, we observed a small peak around H = 0 both above and below TN and both along a- and c-axis, possibly due to an unknown paramagnetic impurity. For data at all temperatures presented in Fig. 2, a background measured at 30 K was subtracted. A small asymmetry of the peak position in the dM/dH vs H was observed when only one loop of M vs H was measured. We did not see such an asymmetry when a second loop of the M vs H was measured, which is presented in Fig. 2. Transport measurements were performed on a quantum design PPMS following a conventional 4-probe method. Au wires of 25 μm diameter were attached to the sample with Epotek H20E silver epoxy. An electric current of 1 mA was used for the transport measurements. In magnetoresistance measurements, the contact misalignment was corrected by field symmetrizing the measured data. The following method was adopted for the contact misalignment correction in Hall effect measurements. The Hall resistance was measured at H = 0 by decreasing the field from the positive magnetic field (R), where H represents the external magnetic field. Again the Hall resistance was measured at H = 0 by increasing the field from negative magnetic field (R). Average of the absolute value of (R) and (R) was then subtracted from the measured Hall resistance. The conventional antisymmetrization method was also used for the Hall resistance measured at 28 K (above TN) and at 2 K (where no anomalous Hall effect was observed), which gave same result as obtained from the former method.

Electronic structure calculations

The electronic structure calculations were carried out within density functional theory (DFT) using the all-electron, full potential code WIEN2K[35] based on the augmented plane wave plus local orbital (APW + lo) basis set[36]. The Perdew–Burke–Ernzerhof (PBE) version of the generalized gradient approximation (GGA)[37] was chosen as the exchange correlation potential. Spin–orbit coupling (SOC) was introduced in a second variational procedure[38]. A dense k-mesh of 28 × 28 × 12 was used for the Brillouin zone (BZ) sampling. A RMTKmax of 7.0 was chosen for all calculations. Muffin tin radii were 2.5 a.u. for Nb, 2.45 a.u. for Co, and 2.01 a.u. for S.

Data availability

The authors declare that the main data supporting the findings of this study are available within the article and its Supplementary Information files. Extra data are available from the corresponding author on request.