Dirac crystals have attracted great attention in the past years [1] . Topological semimetals with linear energy dispersion in momentum space host variety of quantum transport properties, such as Quantum Hall Effect (QHE), extremely large magnetoresistance (MR) and high mobility [2] [3] [4] . In particular, manipulation of topologically protected Dirac states by antiferromagnetic (AFM) spins in crystals not only gives rise to rich quantum effects but also enables stable devices that dissipate less energy when compared to traditional FM heterostructures [5] [6] [7] .
In contrast to local-moment compounds with felectron levels deep below Fermi energy or heavy fermions that feature magnetic spin exchange with conduction electron bands, mixed i.e. non-integer chemical valence (MV) on 4f or 5f electronic orbitals facilitates charge transfer and absence of long-range magnetic order [8, 9] . Charge fluctuations exert considerable influence on the conduction electron bands in general and in particular on the formation of Dirac states at Topological Kondo Insulator surfaces [10] [11] [12] [13] . On the other hand, magnetic MV compounds with Dirac topological states are hitherto unknown.
Here we report discovery of a new material, EuZnSb 2 , with 4f magnetic order below T N = 20 K. Magnetic moment of Eu is estimated to be 3.5 µ B , smaller to that of 7.9 µ B in EuMnBi 2 due to the presence of MV Eu 2+ and Eu 3+ atoms, similar to SmB 6 [13] . Angular-and temperature-dependent magnetoresistance (MR) are consistent with quasi-two-dimensional Fermi surface. First-principles calculations reveal that square lattice Sb 5p x /5p y -derived bands have a Dirac point between Γ and M points in the Brillouen zone. This in-dicates both the topological electronic bands and magnetic texture. Such arrangement in a tunable ABX 2 (A = alkaline earth and/or rare earth, B = transition metal) crystallographic structure where graphenelike two-dimensional (2D) quantum charge transport and high mobility were found [14] [15] [16] [17] [18] [19] [20] can lead to the emergence of rich quantum phases [5, 21] .
Single crystals of EuZnSb 2 were grown by the Bridgeman method. First, Eu chunks, Zn particles, and Sb lumps were mixed in a stoichiometric ratio, placed into alumina crucible and then sealed in a fused silica tube. The quartz tube was slowly heated to 1030 • C. After 2 hours, the sample was fast-cooled in four hours to 850 • C and then slow-cooled to 400 • C at a rate of 3 • C/h. Finally, the quartz tube was cooled to room temperature with the furnace power shut off. Single crystals with size up to 4 × 2 × 0.5 mm 3 can be cleaved from the melted ingot.
Magnetotransport and heat capacity measurement up to 9 T were measured in a Quantum Design PPMS-9. Resistivity ρ xx was measured by a standard four-probe method. Hall resistivity ρ xy was measured by four terminal technique by switching the polarity of the magnetic filed to eliminate the contribution of ρ xx .
Powder X-ray diffraction (XRD) data were taken with Cu K α (λ = 0.15418 nm) radiation of Rigaku Miniflex powder diffractometer. The element analysis was performed using an energy-dispersive X-ray spectroscopy (EDX) in a JEOL LSM-6500 scanning electron microscope, confirming stoichiometric EuZnSb 2 . The X-ray absorption spectroscopy measurement was performed at 8-ID beamline of NSLS II (BNL). The X-ray absorption near edge structure (XANES) and the extended x-ray absorption fine structure (EXAFS) spectra were processed using the Athena software package. The extracted EX-AFS signal, χ(k), was weighed by k 2 to emphasize the high-energy oscillation and then Fourier-transformed in a k range from 2 to 11Å −1 to analyze the data in the R space.
X-ray photoemission spectroscopy (XPS) experiments were carried out in an ultrahigh vacuum (UHV) system with base pressures less than 2·10 −9 Torr equipped a hemispherical electron energy analyzer (SPECS, PHOI-BOS 100) and twin anode X-ray source (SPECS, XR50). Al K α (1486.7 eV) radiation was used at 10 kV and 30 mA. The angle between the analyzer and X-ray source is 45 • and photoelectrons were collected along the sample surface normal.
For accurate crystal structure determination, single crystal intensity data sets were collected at 200 K in Bruker SMART APEX II CCD diffractometer using graphite monochromatized Mo K α radiation (λ = 0.71073Å). Approximately a quarter of sphere of reciprocal space data was collected in 2 batch runs at different ω and φ angles with an exposure time of 6 sec/frame. A total of 1415 reflections (2θ max ∼56 • ) were collected, 191 of which were unique (T min /T max = 0.170/0.327, R int = 0.066). The data collection, data reduction and integration, as well as refinement of the cell parameters were carried out using the Bruker-provided programs with applied semi-empirical absorbtion correction [22, 23] . The structure was subsequently solved by direct methods and refined on F 2 (12 parameters) with the aid of the SHELXL package [24] . All atoms were refined with anisotropic displacement parameters with scattering factors (neutral atoms) and absorption coefficients [25] . The final Fourier map is featureless with the highest residual density and deepest hole of about 3.5 e − /Å 3 , situated 0.9 A and 1.3Å away from Sb2 and Sb1, respectively. Full occupancies of all sites were confirmed by refinement of the site occupation factors for all sites, none of which deviated from unity within more than 3 standard deviations. Note that Zn (Z = 30) and Sb (Z = 51) have sufficiently different X-ray scattering factors [25] , which makes them clearly distinguishable and ascertains the established, devoid of positional disorder structural model.
For first-principles band structure calculations, we applied the WIEN2K [26] implementation of the full potential linearized augmented plane wave method in the generalized gradient approximation (GGA) [27] + U eff = 6 eV on the Eu 4f orbitals [28] [29] [30] with the spin-orbit (SO) coupling treated in the second variation method. The basis size was determined by R mt K max = 7 and the Brillouin zone was sampled with a regular 18 × 18 × 3 mesh containing 162 irreducible k points to achieve energy convergence of 1 meV. A 10000 k-point mesh was used for the Fermi surface calculations.
The EuZnSb 2 refined structure [ Fig. 1(a,b) , Tables  I, II] is ZrCuSiAs type, Pearson index tP 8. EuZnSb 2 is a new ternary compound in the Eu-Zn-Sb phase diagram. The other known ternary compounds are EuZn 2 Sb 2 , Eu 2 ZnSb 2 , Eu 21 Zn 4 Sb 18 and Eu 11 Zn 6 Sb 12 [31] [32] [33] [34] . EuZnSb 2 is isotypic to ZrCuSiAs [space group P 4/nmm (No. 129), Z=2]. There are close to 150 compounds of this type, whose structures are confirmed and deposited in structural databases. Of those, only one other pnictide containing an alkaline-earth metal forms with the same structure and it is the bismuthide CaMnBi 2 [15, 35] . No other Eu-based pnictides are known to exist with this structure. This is also true for the EuMnBi 2 [36] which is structurally similar, but not identical to EuZnSb 2 . Strontium, whose solid state chemistry often mirrors that of Eu does not form pnictide phase isotypic to EuZnSb 2 . A Sr-analog with the 1:1:2 stoichiometric ratio, SrZnSb 2 is known, but forms with the orthorombic CaMnSb 2 structure type (P nma, Pearson index oP 16 [37] . Importantly, the vast majority of pnictides with the tetragonal ZrCuSiAs structure type are with trivalent rare-earth metals and their structures show a discernable tendency to be transition metaldeficient. A good and relevant example of this notion is the series of RZn 1−x Sb 2 (R = La, Ce, Pr, Nd, Sm, Gd, Tb), where the Zn-site can be vacant with amount of defects being as high as 40 % [38] . In this regard, EuZnSb 2 represents a rare compound with the tetragonal ZrCuSiAs/HfCuSi 2 structure, not only considering its chemical makeup, but also the fact that it is a zincpnictide from this family that is apparently as close to being a stoichiometric phase as possible.
The structure prototype is known and will not be described at length. Results of structural refinement are shown in Table I . There are four independent sites in the asymmetric unit (Table II) and based on the refinements, the structure appears to be devoid of disorder on any of them. The nearly spherical shape of the anisotropic displacement parameters is also a testament to this conjecture. The structure is based on PbO-type layers made up of fused ZnSb 4 -tetrahedra and square-nets of Sb atoms [ Fig. 1(a-b) ]. All Zn-Sb interatomic distances are within the common range (Table II) and within the range for the sum of the single-bonded covalent radii. The Sb-Zn-Sb angles within the ZnSb 4 -tetrahedra deviate from the ideal 109.5 • value and range from 104.4(2) • to 112.1(1) • . There are no indications from the structure refinements that the Sb square lattice is a subject to a Peierls distortion. The Sb-Sb interatomic distances within the square nets are longer than what is expected for a single covalent Sb-Sb bond, but the atomic interactions in such topology are hypervalent, consistent with the longer distances We note that ABX 2 (A = alkaline earth and/or rare earth, B = transition metal) materials feature rather tunable crystal structure. The square Bi or Sb slabs host graphene-like Dirac states of X=Bi/Sb orbital character with charge-or spin-density-wave order coupled to topological states [39] [40] [41] [42] [43] . Weyl states were proposed in Sr 1−y Mn 1−z Sb 2 and YbMnBi 2 [44, 45] . Topological Dirac states were also found in Zn-based 112 materials such as BaZnBi 2 and SrZnSb 2 [46, 47] ; EuMnBi 2 hosts QHE with field tunable Eu 4f magnetic order and magnetopiezoelectric effect [48, 49] .
Figures 1(c) and 1(d) show the Zn K-edge oscillations and corresponding Fourier transform magnitudes of Extended X-ray absorption fine structure (EXAFS) spectra of EuZnSb 2 , respectively. In the single-scattering approximation, the EXAFS could be described by the fol-lowing equation [50] :
where N i is the number of neighbouring atoms at a distance R i from the photoabsorbing atom. S 2 0 is the passive electrons reduction factor, f i (k, R i ) is the backscattering amplitude, λ is the photoelectron mean free path, δ i is the phase shift of the photoelectrons, and σ 2 i is the correlated Debye-Waller factor measuring the mean square relative displacement of the photoabsorber-backscatter pairs. The corrected main peak around R ∼ 2.77Å in Fig. 1(d) corresponds to the Zn-Sb1 bond distances in ZnSb 4 tetrahedra and is in good agreement with single crystal refinement result (Table II) . The Eu atoms are coordinated by eight nearest neighbor antimony atoms in a square-antiprismatic fashion [ Figure 1 (e)] with Eu-Sb1 and Eu-Sb2 distances longer than 3.3Å (Table II) . Unit cell from the powder X-ray diffraction experiment on pulverized crystals [ Fig. 1(f) ] can be fitted well with the P 4/nmm structural model (Table I) , confirming single crystal refinement result and phase purity.
Temperature dependence of magnetic susceptibility measured with zero field cooling and field cooling mode is shown in Fig. 1(g) for two primary crystallographic directions. A sharp peak at 20 K was observed for the magnetic susceptibility with H//a, while a saturation behavior is observed below 20 K with H//c, indicating that the magnetic easy axis is along the a-axis. The peak at 20 K can be explained by the antiferromagnetic transition of Eu 2+ . Magnetic susceptibility above 20 K exhibits a typical Curie-Weiss behavior, which can be well fitted to χ = χ 0 + C T −θ , where χ 0 is the temperature independent Pauli contribution and C is related to the effective moments. The effective moment of Eu 2+ is 3.5 µ B and θ = -9.4 K for H//ab while 3.7 µ B and θ = -11.2 K for H//c. Negative Curie temperatures indicate dominant antiferromagnetic interactions of Eu 2+ . In order to obtain more information about the magnetic structure of Eu 2+ , we perform the field dependent measurement of magnetization loop, as shown in the inset in Fig. 1(g) . The magnetization along the c-axis exhibits perfect linear behavior at 2 K and 30 K, i.e. at temperatures below and above T N . When field is applied along the ab plane, the magnetization shows linear behavior in 30 K whereas a very small deviation from linearity is observed at 1.1 T in M(H) taken at 2 K. The temperature and field dependence magnetization of EuZnSb 2 [ Fig. 1(g) ] resembles that of EuZnBi 2 whose magnetic structure is analogous to EuMnBi 2 spin-flop AFM phase [36, 49] . However, these Bi-based materials belong to I4/mmm space group containing different Eu-sublattice structure and fixed Eu valence when compared to EuZnSb 2 . The weak in-plane anisotropy of Eu 2+ gives rise to the kinks around 1. no spin-flop is observed, in contrast to that in EuZn 2 Sb 2 [51] . Thus, the magnetic field that might induce the spinflop of Eu 2+ can be higher than 20 T [49] . X-ray photoelectron 3d core level spectra [ Fig. 1(h) ] shows spin-orbit split states 3d 2+ 5/2 , 3d 2+ 3/2 , 3d 3+ 5/2 and 3d 3+ 3/2 , at binding energies 1125, 1155, 1135 and 1165 eV, respectively, confirming the presence of both Eu 2+ and Eu 3+ . We also observe a weak satellite peak at somewhat higher binding energies from 3d 3+ 5/2 , originating from the multielectronic excitations in the photoelectron emission. The XPS spectra is rather similar to Eu 2 SrBi 2 S 4 F 4 [52] and provide explanation for the observed paramagnetic Curie-Weiss moment reduced from 7.9 µ B . Interestingly, both Eu 2+ and Eu 3+ share the same 2c atomic site in the P 4/nmm unit cell; this might reduce Zn vacancies when compared to La 1−x ZnSb 2 [53] . The absence of Zn defects and mixed valent Eu might suggest that the squareplanar nets in EuZnSb 2 can not be treated as hypervalent Sb since for trivalent rare earth defects on the Zn site are necessary to lower the valence electron count [53, 54] .
Temperature dependence of the specific heat of EuZnSb 2 is shown in Fig. 2(a) . The anomaly corresponding to the AFM transition of Eu 2+ is observed below 20 K [ Fig. 2(a) inset]. The anomaly is suppressed to low temperature by magnetic field, consistent with resistivity and magnetization measurement. Room temperature heat capacity is close to 3N R, where N is the atomic number per chemical formula and R is the universal gas constant 8.314 J mol −1 K −1 . Specific heat above 25 K is well fitted by the Debye-Einstein model [55] :
Where θ D and θ E are the Debye and Einstein temperatures, respectively, and α denotes the relative contribution of Debye and Einstein terms to phonon heat capacity. The θ D = 233 K and θ E = 73 K can be obtained from the fitting. We note that inclusion of Einstein term was necessary to fit C(T ), suggesting the presence of optical phonon modes [56, 57] . The entropy change induced by the AFM transition of Eu 2+ is obtained by ∆S = (C p − C ph+el )/T dT , and ∆S is estimated to be 11 J mol −1 K −1 , smaller than the theoretical value ∆S = R ln(2J + 1) = 17.3 J mol −1 K −1 for Eu 2+ with J = 7 2 . However, we can infer the J = 1.38 from ∆S = 11 J mol −1 K −1 . Then, the effective moment of Eu 2+ is estimated to be 3.6 in agreement with the magnetization measurement and also with XPS investigation that reveals the presence of both Eu 2+ and Eu 3+ .
Hall resistivity ρ xy exhibits positive slope and small temperature dependence [ Fig. 2(b) ], suggesting dominant hole carriers. Clear nonlinear behavior at high field indicates multiband transport, similar to YbMnSb 2 [58] . We fit ρ xy with semiclassical two-band model [59] :
where n e (n h ) and µ e (µ h ) denote the carrier concentrations and mobilities of electrons and holes, respectively. The obtained carrier concentration and mobility are shown in Fig. 2 (c) and 2(d), respectively. The hole concentration is nearly one order of magnitude larger than that of electrons whereas at low temperatures electron-type carriers have much larger mobility, in contrast to values at 200 K. Resistivity of EuZnSb 2 decreases on cooling; anomaly corresponding to the AFM transition of Eu 2+ is observed at 20 K [ Fig. 3(a) ]. There is an upward shift of the temperature-dependent resistivity and the T N is suppressed to 12 K in 9 T applied along the c-axis [θ = 0 in Fig. 3(a) inset]. The resistivity is similar to that of AZnSb 2 (A = Ba and Sr) and EuZnBi 2 , but is in contrast Angle θ dependence of resistivity at 2 and 30 K, respectively measured in 9 T. (c,d) Magnetoresistance of EuZnSb2 for θ = 0. Green lines show fits using two-band orbital magnetoresistance obtained using carrier concentrations and mobilities from Fig. 2(c,d) ; we note that such fits are possible for carrier concentration and mobility variations for up to ±10% and ±20% from values reported in Fig. 2 . (e) Magnetoresistance at different tilt angles θ between B and c-axis. Data were measured at 2 K and the current was always perpendicular to the magnetic field.
to semiconducting EuMnSb 2 [60] [61] [62] . Since MR is determined by the extremal cross section of the Fermi surface perpendicular to the magnetic field direction, information about geometric structure of Fermi surface can be obtained from the angular dependence of MR. As shown in Fig. 3(b) , MR angular dependence is anisotropic. It shows a typical 2D behavior, B ∝ |cosθ| [14, 15, 60] . Similar angular dependence at 2 K and at 30 K suggests that the magnetic order of Eu 2+ is only a small perturbation on ρ(θ).
Magnetic field dependence of resistivity at different temperatures for θ = 0 is shown in Fig. 3(c,d) . The maximum MR ratio, MR = [ρ(B) − ρ(0)] /ρ(0)×100%, is 10 % at 50 K and 9 T, which is relatively small. MR(B) hints towards a low field B 2 dependence to high field linear MR crossover, as is sometimes observed in Dirac materials [14, 15, 60, 61] . Semiclassical MR in metals features B 2 dependence in the low field and a saturating MR in high fields [63] .
Strong external magnetic field might lead to a complete quantization of the perpendicular orbital motion of carriers with linear energy dispersion and quantized Landau Levels E n = sgn(n)v F 2e B|n| where n = 0, ±1, ±2, · · · is the LL index and v F is the Fermi velocity [64, 65] . In such quantum limit, the energy difference from the lowest to the 1 st Landau level is △ LL = ±v F √ 2e B; △ LL is larger than Fermi energy E F and the thermal fluctuations k B T so that all carriers are at the lowest Landau level where linear MR is expected [63, 66] . Crossover magnetic field B * above which the quantum limit is satisfied at specific temperature T B * = 1 [67] could then be used to estimate Fermi velocity. Such analysis, assuming quantum limit, gives the Fermi velocity v F ∼ 5.13 × 10 5 ms −1 , i.e. ∆ LL = 5 meV . We estimate Fermi energy from E F = 2 /m 0 (3π 2 n) 2/3 where n is the carrier density. By taking measured carrier density [ Fig. 2(c) ] we obtain E F = 16.5 meV for electron pocket and E F = 250 meV for the hole pocket, suggesting that MR arises from the first few Landau levels. Indeed, the two-band orbital magnetoresistance MR: [68] 
is satisfactory only at high temperatures [ Fig. 3(c) ].
Other possible reasons for linear MR in high fields include mobility fluctuations in an inhomogeneous crystal [69] or open orbits. The former can be excluded since EuZnSb 2 is stoichiometric crystal. The latter usually arises for electronic motion on orbits associated with magnetic field oriented along elongated necks of the Fermi surface, for example in Cu [70] , in topological materials with extremely large magnetoresistance with compensated charge carriers [71] or in two-dimensional conductors where magnetic field is applied parallel to the conducting layers [72] . Whereas we note that EuZnSb 2 is not a compensated metal with equal electron and hole concentration [ Fig. 2(c) ], and that in Fig. 3 (c,d) magnetic field is oriented orthogonal to quasi-2D conduction direction, we can not exclude the possibility for open Fermi surface pockets. It also should be noted that small negative MR appears when magnetic field B is tilted away from c-axis at 2 K [ Fig. 2(d) ]. This likely arises due to in-plane spin reorientation which could indicate shift of the magnetic easy axis with external field rotation.
First-principles band structure calculations reveal the mechanism of Dirac point formation and point to possible magnetic structure after comparing the total energy for five different magnetic structures [ Fig. 4(a-e) , Table III and (c) for the definition of the exchange paths]. These numbers indicate that the AAF3 pattern is the ground state. Fig.4(e) shows the low-energy magnetic structure with the in-plane easy axis obtained by GGA+U +SO calculations with U = 6 eV for the Eu 4f orbitals. However, the state is almost degenerate with the one with the out-of-plane easy axis, meaning that the magnetic anisotropy is very weak. The moments of Eu 2+ feature ferromagnetic arrangement in plane, and are stacked antiferromagnetic by every two layers in the out plane direction. In its corresponding electronic structure [ Fig.4 Fig.4(g) ] shows a coronavirus-like electron pocket at Γ point surrounded by four hole pockets. The Dirac states around (π/2, π/2) are flat along the c-axis, i.e. they are restricted to the Sb2 square lattice. Overall, the volume of the hole pockets are larger than the electron pocket, in agreement with the experiments. The structure at the Γ point features elongated necks which could contribute to open orbits and linear MR, as discussed above. The antiferromagnetic J Sb2 (i.e., the exchange interaction between two nearest Eu ions above and below the square-lattice Sb2 layer) means that the time reversal symmetry is not broken for the itinerant electrons on the Sb2 layers at the level of the mean-field-like calculations, leading to the four-fold degenerated Dirac states. Yet, it would have been broken by the random distribution of Eu 2+ and Eu 3+ ions on the same crystal 2c site-leading to the formation of Weyl states-which cannot be described by standard DFT calculations with the 1 × 1 × 2 minimal structure model. Further stud-ies of topological electronic structure and putative Weyl states in EuZnSb 2 using angular resolved photoemission (ARPES) and large-supercell DFT calculations are of high interest.
In conclusion, we report discovery of EuZnSb 2 , a new magnetic semimetal. EuZnSb 2 exhibits MV Eu, i.e. both Eu 2+ and Eu 3+ are present in the unit cell. The Eu 2+ moments order antiferomagnetically below 20 K. Firstprinciple calculations are consistent with the presence of Dirac states in the band structure whereas magnetotransport suggests dominant quasi-2D Fermi surface sheets. Further ARPES and neutron experiments are of interest to shed more light on Dirac dispersion in momentum space and magnetic space group in the ordered state. Moreover, spectroscopic studies are of interest to investigate possible coupling of dynamical charge fluctuations and the long-range magnetic order with Dirac bands.
Work 