A quantum spin-liquid is a spin disordered state of matter in which spins are strongly correlated and highly entangled with low-energy excitations. It has been often found in two-dimensional S = ½, highly frustrated spin networks but rarely observed in three-dimensional (3D) frustrated quantum magnets. Here, KSrFe2(PO4)3, forming a complicated 3D frustrated lattice with a spin moment S = 5/2, is investigated by thermodynamic, neutron diffraction measurements and electronic structure calculations. Despite the relatively sizable Curie–Weiss temperature θCW = −70 K, a conventional magnetic long-range order is confirmed to be absent down to 0.19 K. The magnetic heat capacity data follow the power-law behavior at the lowest temperature region, supporting gapless excitations in a 3D spin-liquid state. Strong geometrical spin frustration responsible for the spin-liquid feature is understood as originating from the almost comparable five competing nearest-neighbor antiferromagnetic exchange interactions, which form the complicated 3D frustrated spin network. All these results suggest that the compound KSrFe2(PO4)3, representing a unique 3D spin frustrated network, could be a rare example of forming a gapless spin-liquid state even with a large spin moment of S = 5/2.
I. INTRODUCTION
Geometrically frustrated antiferromagnetic (GFAFM) systems have attracted considerable attention in condensed matter physics due to their exotic magnetic properties, such as the spin-liquid state. In those GFAFM systems, spins are put into a bizarre situation that does not allow a single conventional antiferromagnetic (AFM) ground state because they fail to satisfy all the AFM interactions simultaneously. Moreover, it has been known that the materials with a small spin value of S = ½ bear strong zero-point quantum fluctuation even at T = 0 K.1,2 Therefore, the spin lattices with geometrically frustrated interaction, with low dimensions and low spin values, have been actively investigated in the search for exotic spin liquid states. A few ideal examples for the two-dimensional (2D) edge-shared triangular lattice and 2D corner-shared triangular (kagome) lattice are κ-(BEDT-TTF)2Cu2(CN)3 and ZnCu3(OH)6Cl2, respectively.3–7 These materials are promising examples of the quantum spin-liquid (QSL) state with S = 1/2 moments and exhibit strong quantum fluctuations. As quantum fluctuations are known to be dominant in low-dimensional materials and diminish in higher dimensions, it has been challenging to discover a QSL candidate with no static spin ordering in three-dimensional (3D) systems, especially with large spins, such as S = 5/2.
Recently, emergent magnetic behavior has been observed in 3D geometrically frustrated magnets due to the interplay between competing spin interactions and quantum fluctuations. A few examples of 3D GFAFM systems can be found in corned-shared triangle (hyperkagome and trillium) lattices and corned-shared tetrahedral (pyrochlore) lattices; the 3D corner-shared triangular networks include Na4Ir3O8,8–10 PbCuTe2O6,11–14 Yb3Ga5O12,15 and K2Ni2(SiO4)3.16 Na4Ir3O8 has a 3D hyperkagome structure with Ir4+ ions forming a 3D corner-shared triangular network with z = 4, where z is the number of nearest neighbor (NN) atoms.8 Ir4+ (5d5) bears the magnetic moment of Jeff = ½ due to the strong spin–orbit coupling. Initially, magnetic data analysis shows no long-range order (LRO), and the heat capacity follows T2 behavior, indicating a formation of a spin-liquid phase with gapless excitations.8,9 The recent nuclear magnetic resonance (NMR) and neutron scattering studies reveal slowing down of spin fluctuations and a frozen spin state in the system.17,18 On the other hand, the system PbCuTe2O6 has an S = ½ (Cu2+), 3D corner-shared triangular lattice.11 Bulk magnetic susceptibility suggests that it is a paramagnetic system with a Curie–Weiss temperature (θCW) of −22 K. Neither magnetic ordering nor spin freezing is seen down to 20 mK. Electronic structure calculations suggest that it has a 3D hyperkagome network with additional frustrated couplings. Low-temperature NMR, muon-spin rotation, and inelastic neutron scattering studies have revealed that it is a spin-liquid system with gapless, fractionalized spinon excitations.12,13 Recently, a ferroelectric ordering with quantum critical behavior has been observed at 1 K due to the strong lattice distortion.14
In search of 3D QSLs, recently, interest has been raised in another type of 3D frustrated system called the trillium lattice, which constitutes a 3D network of corner-shared triangles with z = 6. The trillium lattice is usually found in B20 structures. A few example materials for B20 structures are MnSi, FeGe, and FeSi.19–23 Due to the lack of inversion symmetry, these trillium lattice systems exhibit skyrmion-type ordering.24 A few more frustrated trillium lattice network systems include CeIrSi,25 demonstrating a spin-ice ground state with Ising interactions, and EuPtX (Ge, Si),26 showing strong magnetic frustration with a finite LRO at 4.1 K. Recently, the field-induced quantum spin liquid behavior has been found in an S = 1 3D frustrated system K2Ni2(SO4)3. Magnetic susceptibility data follow the Curie–Weiss behavior with θCW = −18 K, suggesting AFM interactions. Heat capacity data (Cp) show the second order transition at 1.14 K and the first order transition at 0.74 K. The magnetic field H suppresses these transitions, and Cp data, following T2 behavior, suggest the field-induced QSL ground state with gapless spinon excitations. It was claimed that K2Ni2(SO4)3 is an inter-connected trillium lattice of two Ni2+ ions; however, the non-trillium coupling is, indeed, more dominant than trillium interactions, according to the reported electronic structure calculations.16 The presence of dominant non-trillium couplings could render the system exhibit LRO at finite temperatures.
To explore exotic properties in the class of 3D complicated network systems, we have investigated a compound KSrFe2(PO4)3 that has a similar structure to K2Ni2(SO4)3 but has a spin S = 5/2. Quite interestingly, no long-range order (LRO) is found down to 0.19 K although it has relatively strong AFM correlation with θCW = −70 K. Magnetic heat capacity follows the power law T α (α ∼ 2.3) behavior, suggesting that our titled compound is a possible candidate of forming a spin-liquid state with gapless excitations in the family of 3D frustrated spin network with S = 5/2.
II. METHODS
Polycrystalline samples of KSrFe2(PO4)3 were synthesized using the conventional solid-state heating method at 1000 °C.27 The non-magnetic compound KSrSc2(PO4)328 with the same structure as KSrFe2(PO4)3 was prepared at 1100 °C. Polycrystalline powder samples were used to measure x-ray diffraction (XRD), magnetization, and heat capacity measurements. The magnetic and heat capacity measurements on polycrystalline samples were carried out using the Physical Properties Measurement System (PPMS, Quantum Design). The low-temperature heat capacity was measured on a tiny pellet using a CF-150 dilution fridge (Leiden cryogenics). The Rietveld refinement of the XRD profile has been done using FULLPROF SUITE software [see Fig. 1(a)]. The extracted lattice parameters are a = 9.780 Å and χ2 = 2.30, Rp = 20.2%, Rwp = 19.2%, and Rexp = 12.63%. As reported previously,27 there is nearly 20% of site exchange between K and Sr atoms.
(a) The refinement of x-ray diffraction of the titled compound KSrFe2(PO4)3. (b) Crystal structure of KSrFe2(PO4)3 in a unit cell with Fe1O6, Fe2O6 octahedral, and PO4 tetrahedral units.
(a) The refinement of x-ray diffraction of the titled compound KSrFe2(PO4)3. (b) Crystal structure of KSrFe2(PO4)3 in a unit cell with Fe1O6, Fe2O6 octahedral, and PO4 tetrahedral units.
III. STRUCTURAL DESCRIPTION
The compound KSrFe2(PO4)3 holds a cubic crystal system with space group P213 (space group number 198). The 3D framework is built by Fe1O6 and Fe2O6 octahedral units, PO4 tetrahedral units, K, and Sr atoms, as shown in Fig. 1(b). The Fe1O6 octahedral unit is slightly distorted than the Fe2O6 octahedral unit. The magnetic interactions are exchanged among Fe1 and Fe2 atoms through PO4 units. The first, second, and third NN couplings (J1, J2, and J3) are between Fe1 and Fe2 atoms. The combination of first NN and second NN couplings forms a J1–J2 anisotropic diamond lattice, known to be a non-frustrated system. After adding the third NN coupling (J3), the system still remains non-frustrated, as shown in Fig. 2(a). On the other hand, the fourth and fifth NN couplings (J4 and J5) are through Fe1–Fe1 and Fe2–Fe2 bonds, respectively (see Table I). The couplings J4 and J5 alone built up two independent trillium lattices, respectively [see Fig. 2(b)], which are highly frustrated 3D triangular networks. Overall, the compound KSrFe2(PO4)3 has the combination of all exchange couplings J1, J2, J3, J4, and J5 and stands for a complicated 3D frustrated network system, as shown in Fig. 2(c). The super-super-exchange interaction between Fe3+ ions is through the O–P–O path. The details of various magnetic couplings are listed in Table I. The magnetic exchange energies were extracted from DFT calculations, which will be discussed in Sec. IV.
(a) The non-frustrated 3D network formed by J1, J2, and J3 couplings between Fe1 (magenta) and Fe2 (blue) atoms. (b) The formation of two trillium networks by J4 and J5, respectively. (c) The complicated 3D network is formed by J1, J2, J3, J4, and J5 couplings.
(a) The non-frustrated 3D network formed by J1, J2, and J3 couplings between Fe1 (magenta) and Fe2 (blue) atoms. (b) The formation of two trillium networks by J4 and J5, respectively. (c) The complicated 3D network is formed by J1, J2, J3, J4, and J5 couplings.
Details of magnetic exchange couplings in KSrFe2(PO4)3.
Magnetic couplings . | Magnetic path . | Bond distance (Å) . | Exchange path . | Exchange energy (K) . | Relative energy (Ji/J5) . |
---|---|---|---|---|---|
J1 | Fe1–Fe2 | 4.55 | Fe1–O4–P–O1–Fe2 | −4.87 | 0.91 |
J2 | Fe1–Fe2 | 4.92 | Fe1–O2–P–O3–Fe2 | −2.67 | 0.5 |
J3 | Fe1–Fe2 | 5.97 | Fe1–O2–P–O4–Fe2 | −4.99 | 0.93 |
J4 | Fe1–Fe1 | 6.03 | Fe1–O2–P–O1–Fe1 | −3.94 | 0.73 |
J5 | Fe2–Fe2 | 6.10 | Fe2–O4–P–O3–Fe2 | −5.34 | 1 |
Magnetic couplings . | Magnetic path . | Bond distance (Å) . | Exchange path . | Exchange energy (K) . | Relative energy (Ji/J5) . |
---|---|---|---|---|---|
J1 | Fe1–Fe2 | 4.55 | Fe1–O4–P–O1–Fe2 | −4.87 | 0.91 |
J2 | Fe1–Fe2 | 4.92 | Fe1–O2–P–O3–Fe2 | −2.67 | 0.5 |
J3 | Fe1–Fe2 | 5.97 | Fe1–O2–P–O4–Fe2 | −4.99 | 0.93 |
J4 | Fe1–Fe1 | 6.03 | Fe1–O2–P–O1–Fe1 | −3.94 | 0.73 |
J5 | Fe2–Fe2 | 6.10 | Fe2–O4–P–O3–Fe2 | −5.34 | 1 |
IV. RESULTS
A. Magnetic measurements
Temperature-dependent magnetic susceptibility χ(T) measurements were carried out on polycrystalline samples of KSrFe2(PO4)3 down to 2 K at H = 10 kOe, as shown in Fig. 3(a). No conventional magnetic LRO is observed in the χ(T) data down to 2 K. A minor splitting between zero-field-cooled (ZFC) and field-cooled (FC) susceptibility data is noted at 3.5 K in the low fields, as shown in the inset of Fig. 3(a). However, this splitting disappears at a small H = 5 kOe. It could be due to the site exchange between K and Sr atoms, as also noted in many geometrically frustrated magnets, such as Na4Ir3O88 and LiGaCr4S8,29 due to the disorder. The inverse χ(T) data follow the Curie–Weiss behavior at high temperatures, as shown in Fig. 3(b). The data fit to the expression 1/χ = 1/[χdia + χvv + C/(T − θCW)]. The extracted θCW is about −70 K, suggesting dominant AFM couplings. The calculated diamagnetic susceptibility (χdia) from the individual ions of K+, Sr2+, Fe3+, P5+, and O2− is −9.75 ⨯ 10−5 cm3/mol.30 The yielded Van Vleck susceptibility χvv from the fit is 1 ⨯ 10−5 cm3/mol-Fe. The calculated effective magnetic moment from the Curie constant (C) = 4.43 cm3 K/mol-Fe is about μeff = 5.95 μB, which is consistent with the expected value for S = 5/2 of the Fe3+. The χ(T) data show field-independent behavior down to 2 K in the magnetic fields of up to 50 kOe [see Fig. 3(c)]. Interestingly, the log–log scale of χ(T) data shows nearly temperature-independent behavior at low temperatures, as seen in the gapless spin-liquid materials, such as PbCuTe2O6, EtMe3Sb[Pd(dmit)2]2, and κ-H3(Cat-EDT-TTF)2.12,31,32 The H-dependence of magnetization M(H) was also measured at 1.9 K by varying H up to 160 kOe, as shown in Fig. 3(d). The magnetization did not saturate and reached only 1.42 μB even at 160 kOe, indicating the presence of AFM correlations in this system. The absence of hysteresis indicates the antiferromagnetic nature of the samples.
(a) Temperature dependence of magnetic susceptibility χ(T) at 10 kOe of KSrFe2(PO4)3. The inset in (a) shows the ZFC and FC χ(T) data at various magnetic fields H = 500 Oe, 5 kOe, and 10 kOe. (b) Inverse magnetic susceptibility χ−1(T) with the Curie–Weiss fit. (c) χ(T) data measured in different fields 10, 20, 30, and 50 kOe. (d) Isotherm magnetization as a function of H of up to 160 kOe at T = 1.9 K.
(a) Temperature dependence of magnetic susceptibility χ(T) at 10 kOe of KSrFe2(PO4)3. The inset in (a) shows the ZFC and FC χ(T) data at various magnetic fields H = 500 Oe, 5 kOe, and 10 kOe. (b) Inverse magnetic susceptibility χ−1(T) with the Curie–Weiss fit. (c) χ(T) data measured in different fields 10, 20, 30, and 50 kOe. (d) Isotherm magnetization as a function of H of up to 160 kOe at T = 1.9 K.
B. Heat capacity measurements
Temperature-dependent heat capacity Cp(T) measurements were done down to 0.19 K at H = 0 and 30 kOe on KSrFe2(PO4)3 samples. Cp(T) for the nonmagnetic analog KSrSc2(PO4)3 is also measured down to 2 K at 0 kOe, as shown in Fig. 4(a). The Cp(T) data of KSrSc2(PO4)3 are well-matched with its Cp(T) data of KSrFe2(PO4)3 at high temperatures. It deviates at low temperatures due to the presence of magnetic contribution of KSrFe2(PO4)3. Following the analysis on KSrSc2(PO4)3,28 the Cp(T) data fit the polynomial AT3 + BT5 in the range from 1.8 to 7 K. It is extrapolated down to 0.19 K. The magnetic heat capacity Cm(T) is obtained after subtracting Cp(T) of KSrSc2(PO4)3 from that of KSrFe2(PO4)3. The Cm(T) data do not show any sharp peak down to 0.19 K, suggesting the absence of magnetic ordering. In particular, Cm(T) does not have an anomaly at 3.5 K, where ZFC-FC magnetic susceptibility shows the bifurcation, which rules out the possibility of magnetic LRO. However, the possibility of spin-glass behavior and its percentage of the frozen spin state are to be further investigated through the detailed Muon spin rotation measurements. The frustration parameter (f), the ratio of Curie–Weiss temperature (θCW) to Néel temperature (TN), is larger than 370. Cm(T) data show a broad maximum at 9 K, indicating the presence of short-range spin correlations [see Fig. 4(b)] as observed in most of the geometrically frustrated magnets. The Cm data are compared with the T2 and T3 power-law fits. The Cm(T) data fall between T2 and T3 behavior, meaning that the power is between 2 and 3. To understand further, we have fitted the Cm data with the power-law, Tα, over the temperature range from 0.19 to 1.7 K, resulting in power-law behavior with a coefficient of α ≈ 2.3 [see Fig. 4(c)]. In general, the Cm(T) data follow the T3 behavior for conventional AFM LRO systems due to 3D magnon excitations, which, thus, rule out the presence of AFM LRO in our system. Conversely, the sublinear (T2/3) dependent behavior is expected for the 2D QSL with spinon Fermi surface. Table II shows the experimental 3D spin-liquid systems following the nearly quadratic behavior. The peculiar power-law behavior of Cm(T) in our system could be due to the spin-liquid nature of the system. Furthermore, detailed theoretical and experimental studies are needed to understand the precise nature of excitations. Despite the application of H = 110 kOe, no field dependence was observed in this compound, suggesting that the low-energy excitations are robust in this spin-liquid candidate.8
(a) Variation of heat capacity Cp of KSrFe2(PO4)3 and KSrSc2(PO4)3 measured at various magnetic fields. (b) Temperature dependence of magnetic heat capacity Cm and the normalized magnetic entropy Sm. The magnetic entropy reaches nearly a maximum R ln(6) (R: universal gas constant) expected value for the S = 5/2 system. (c) The log–log scale Cm vs T data with different power-law fits.
(a) Variation of heat capacity Cp of KSrFe2(PO4)3 and KSrSc2(PO4)3 measured at various magnetic fields. (b) Temperature dependence of magnetic heat capacity Cm and the normalized magnetic entropy Sm. The magnetic entropy reaches nearly a maximum R ln(6) (R: universal gas constant) expected value for the S = 5/2 system. (c) The log–log scale Cm vs T data with different power-law fits.
Details of 3D GFAFM spin-liquid materials.
Compound . | Spin . | Lattice . | θCW (K) . | Tmax in Cm (K) . | α (in Cm ∼ Tα) . | Ground state . | References . |
---|---|---|---|---|---|---|---|
Na4Ir3O8 | Jeff = ½ | 3D hyperkagome | −650 | 30 | 2 | Gapless QSL | 8–10 |
PbCuTe2O6 | S = ½ | 3D hyperkagome network | −22 | 1.15 | 1.9 | Gapless QSL | 11–14 |
with additional interactions | |||||||
K2Ni2(SO4)3 | S = 1 | 3D coupled trillium | −18 | 5 | 2 | Field-induced QSL | 16 |
NaCaNi2F7 | S = 1 | 3D pyrochlore | −129 | 18 | 2.2 | SL | 33,43 |
LiGa0.2In0.8Cr4O8 | S = 3/2 | 3D breathing pyrochlore | −386 | 70 | 1.9 | QSL | 44 |
KSrFe2(PO4)3 | S = 5/2 | 3D frustrated network | −70 | 9 | 2.3 | SL | This work |
Compound . | Spin . | Lattice . | θCW (K) . | Tmax in Cm (K) . | α (in Cm ∼ Tα) . | Ground state . | References . |
---|---|---|---|---|---|---|---|
Na4Ir3O8 | Jeff = ½ | 3D hyperkagome | −650 | 30 | 2 | Gapless QSL | 8–10 |
PbCuTe2O6 | S = ½ | 3D hyperkagome network | −22 | 1.15 | 1.9 | Gapless QSL | 11–14 |
with additional interactions | |||||||
K2Ni2(SO4)3 | S = 1 | 3D coupled trillium | −18 | 5 | 2 | Field-induced QSL | 16 |
NaCaNi2F7 | S = 1 | 3D pyrochlore | −129 | 18 | 2.2 | SL | 33,43 |
LiGa0.2In0.8Cr4O8 | S = 3/2 | 3D breathing pyrochlore | −386 | 70 | 1.9 | QSL | 44 |
KSrFe2(PO4)3 | S = 5/2 | 3D frustrated network | −70 | 9 | 2.3 | SL | This work |
C. Neutron diffraction measurements
To confirm the absence of magnetic LRO, neutron diffraction (ND) experiments were performed on the polycrystalline samples of KSrFe2(PO4)3 down to 2 K. The measurements were carried out using the neutron powder diffractometer PD-I (λ = 1.094 Å) at various temperatures from 6 to 300 K. The ND data 2 and 20 K were obtained at another neutron powder diffractometer PD-II (λ = 1.244 300 Å) with three linear position-sensitive detectors at the Dhruva reactor, Bhabha Atomic Research Center, India. Rietveld refinement on ND data was performed using the FULLPROF software package using the inputs from the structural parameters obtained from the powder XRD data. Neither the additional magnetic Bragg peaks nor enhancement in the intensity of the fundamental nuclear Bragg peaks from 300 to 2 K was found (see Fig. 5). The measured ND pattern could be fitted by considering only the nuclear phase. In addition, no significant difference was found upon high-temperature ND data at 300 K being subtracted from that at 6 K. This suggests that the broad maximum in the heat capacity may stem from short-range correlations. Similarly, no significant difference was seen in the subtracted ND data at 20 K from that at 2 K, ruling out the presence of magnetic LRO down to 2 K.
(a) and (b) Neutron diffraction data recorded with the wavelength of λ = 1.094 Å at temperatures 300 and 6 K. (c) The ND data of 6 K after subtracting the ND data at 300 K. (d) and (e) ND data recorded at the wavelength of λ = 1.2443 Å at temperatures 20 and 2 K. (f) The ND data of 2 K after subtracted the ND data at 20 K. The red closed circles indicate the measured data (Iobs), the black lines indicate the calculated pattern (Ical), the blue lines represent the difference (Iobs − Ical), and the green vertical lines represent the Bragg positions.
(a) and (b) Neutron diffraction data recorded with the wavelength of λ = 1.094 Å at temperatures 300 and 6 K. (c) The ND data of 6 K after subtracting the ND data at 300 K. (d) and (e) ND data recorded at the wavelength of λ = 1.2443 Å at temperatures 20 and 2 K. (f) The ND data of 2 K after subtracted the ND data at 20 K. The red closed circles indicate the measured data (Iobs), the black lines indicate the calculated pattern (Ical), the blue lines represent the difference (Iobs − Ical), and the green vertical lines represent the Bragg positions.
D. Electronic structure calculations
The electronic structure calculations were performed using a full-potential linearized muffin-tin orbital method as implemented in the RSPt code.40–36 For the exchange-correlation functional, we have employed a generalized gradient approximation + Hubbard U (GGA+U) approach. The Brillouin zone integration is done using the thermal smearing method with an 8 × 8 × 8 k-mesh. For charge density and potential angular decomposition inside the muffin tin sphere, maximum angular momentum is considered to be lmax = 8. To account for electronic correlation within the rotationally invariant GGA+U approach,37 we have considered coulomb interactions U = 4.0 eV38 and Hund coupling JH = 0.8 eV for Fe-3d states.
To begin with, we have checked the relative stability of a few possible magnetic states within the GGA+U approach and found that an AFM state where the NN Fe spins are oppositely aligned comes out to be the lowest in energy. The computed total and partial density of states for Fe-d and O-p states are displayed in Fig. 6. The total density of states [Fig. 6(a)] clearly shows that the ground state is insulating with a gap of 2.2 eV. From the projected density of states [Figs. 6(b) and 6(c)], we find that the majority Fe-t2g and eg, states are completely occupied, while the corresponding minority states are completely empty. According to Hund’s rule, such distribution of Fe-states is expected for Fe3+ (d5). The gap primarily arises between the majority and minority d-states, which are exchange split due to the formation of local moments and the correlation effect of the Fe-3d subshell. The magnetic moment on the Fe-site is estimated to be 4.3 µB, which is smaller than the moment expected from the nominal d5 configuration. This also implies strong Fe-d and O-p hybridization in this system, as evident from the spectral distribution of the O-p states [Fig. 6(c)], which appears in a broader energy range between −4 eV and Fermi energy. We will see below that such strong d-p hybridization will facilitate a super-super-exchange mechanism for the magnetic interactions between the Fe-spins, which are even quite far from each other.
(a)–(d) Total and partial density of states in the magnetic ground state using GGA+U calculations (U = 4 eV).
(a)–(d) Total and partial density of states in the magnetic ground state using GGA+U calculations (U = 4 eV).
Furthermore, we have attempted to understand the nature of magnetism and employed the magnetic force theorem to extract magnetic exchange energies.45–41 The computed inter-site exchange interactions are displayed in Table I. Surprisingly, the fifth NN coupling (J5) is dominant among all five exchange couplings although the bond distance of J5 is larger than those of other couplings, which forms the trillium lattice by the Fe2 atoms [see Fig. 2(b)]. As per the structure, the fourth NN coupling J4 also forms another trillium lattice by Fe1 atoms. The magnitude of J4 is slightly smaller than that of J5, which could be due to the distortion in Fe1O6 being slightly larger than that of Fe2O6 environments. In addition to these J4 and J5 couplings, there are strong and compatible AFM couplings J1, J2, and J3. All these spin interactions make a complicated 3D spin system.
V. DISCUSSION
To understand the origin of the highly frustrated magnetism and spin liquid features in the complicated 3D spin network, one needs to look into the details of structural features and magnetic coupling strengths. It is clear from Fig. 2(a) that the system is non-frustrated upon the magnetic couplings of J1, J2, and J3 only being considered. Hence, we need to look at other couplings, J4 and J5, to trace the origin of magnetic frustration. Figure 2(b) clearly shows that J4 and J5 interactions built two uniform equilateral triangular networks (3D perfect trillium lattices). From the DFT electronic structure calculations, it is noted that J5 is a dominant AFM coupling, which is the starting point to explain the geometrical frustration in this compound as the J1–J2 anisotropic diamond lattice is not frustrated even with the inclusion of the J3 coupling. However, the isolated trillium lattice model is insufficient to describe the geometrical frustration of this compound. There exist two trillium lattices formed by J4 and J5 couplings, respectively. The connections between two trillium lattices are made by the three interactions of J1, J2, and J3. These are AFM of whose strengths are almost comparable to those of J4 and J5. Importantly, these AFM inter-trillium connections again produce the triangular-type of interactions between J4 and J5 triangles [see Fig. 2(c)]. These connections make the J4–J2 and J5–J3 anisotropic tetrahedral environments, which are another type of frustrating motifs [see Fig. 2(c)]. Overall, the compound forms a more complicated 3D frustrated spin network than the trillium lattice. The formation of a complex 3D frustrated network with competing and compatible magnetic interactions (J1, J2, J3, J4, and J5) is likely the primary source of high frustration and the related gapless spin liquid feature in the compound.
Since the crystal structure and space group of KSrFe2(PO4)3 are the same as those of K2Ni2(SO4)3, it will be worthwhile to compare the structural and magnetic features of the two compounds. There are a few differences noticed, leading to different magnetic ground states. While K2Ni2(SO4)3 shows magnetic LRO at zero magnetic fields, the titled compound KSrFe2(PO4)3 shows the spin-liquid feature. In the case of K2Ni2(SO4)3, one of the interactions (J2) is, indeed, ferromagnetic; the combination of AFM and FM couplings likely reduces the levels of magnetic frustration. In addition, one of the inter-trillium couplings (J3) in the compound K2Ni2(SO4)3 is nearly three times more dominant than the intra-trillium couplings.16 Unlike K2Ni2(SO4)3, all the magnetic couplings in the KSrFe2(PO4)3 are AFM and comparable with each other, which may enhance the levels of magnetic frustration. The elevated levels of magnetic frustration in the compatible magnetic interactions could lead to the robust spin-disorder in KSrFe2(PO4)3. Moreover, the system has a site disorder between K and Sr atoms and thus, induces, some randomness in magnetic exchange couplings. The randomness could destabilize the LRO and help the system stay in a spin-liquid ground state.42 Our experimental and theoretical results, thus, reveal that the titled compound is a suitable candidate for forming the spin-liquid state with gapless excitations.
VI. CONCLUSION
We have investigated thermodynamic properties and spin interactions by electronic structure calculations in the compound KSrFe2(PO4)3, which forms a 3D frustrated lattice with S = 5/2 magnetic moments. The absence of a conventional magnetic LRO is confirmed down to 0.19 K, despite the relatively sizable Curie–Weiss temperature θCW = −70 K. The power-law T2.3 behavior of the magnetic heat capacity Cm(T) at lower temperatures, followed by a broad maximum at 9 K, supports the formation of a gapless ground state. The appearance of geometrical spin frustration is understood as due to the complicated frustrated AFM exchange couplings, rendering the system remains in an unusual spin-liquid state, at least down to 0.19 K. Our results suggest that the compound KSrFe2(PO4)3 is an appealing candidate for forming a gapless spin-liquid state among the family of S = 5/2, 3D spin networks.
ACKNOWLEDGMENTS
B.K. acknowledges the support from DST, India. The work at SNU was supported by the NRF (Grant No. 2019R1A2C2090648) and the Ministry of Education (Grant No. 2021R1A6C101B418). E.K. acknowledges the financial support from the labex PALM for the QuantumPyroMan project (Project No. ANR10-LABX-0039-PALM).
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
K. Boya and K. Nam contributed equally to this work.
K. Boya: Conceptualization (equal); Investigation (equal); Methodology (equal). K. Nam: Data curation (equal); Formal analysis (equal); Investigation (equal); Validation (equal). K. Kargeti: Investigation (equal); Software (equal); Validation (equal). A. Jain: Data curation (equal); Formal analysis (equal). R. Kumar: Data curation (equal); Formal analysis (equal). S. K. Panda: Software (equal); Supervision (equal). S. M. Yusuf: Data curation (equal); Formal analysis (equal). P. L. Paulose: Data curation (equal); Formal analysis (equal). U. K. Voma: Data curation (equal). E. Kermarrec: Data curation (equal); Formal analysis (equal); Validation (equal). Kee Hoon Kim: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Supervision (equal). B. Koteswararao: Conceptualization (equal); Investigation (equal); Methodology (equal); Project administration (equal); Supervision (equal); Validation (equal); Visualization (equal).
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding authors upon reasonable request.