planes ( 31 ). Another piece of experimental evi-
dence is the critical exponentb= 0.321(3) of
the (1/3, 1/3, 0) peak intensity nearT 2 (fig. S4A),
which indicates the 3D Ising nature of the
magnetic order ( 42 ).
Discussion
Although the Monte Carlo simulations of the
classical spin model described above are in
partial agreement with our experiments, they
do not explain the experimental value of the
magnetic entropySm=10.38J·mol−^1 K−^1 ≈
1.248RatT 2 , which is very different from the
0.231Rgiven by the model. Qualitatively, the
discrepancy should be a result of the thermal
population of multiple low-lying CEF levels
of Ho3+; however, this leads to the question of
why the classical Ising Hamiltonian is appli-
cable. In pyrochlore systems such as Dy 2 Ti 2 O 7 ,
the classical Ising behavior is a consequence of
the dominance of CEF splitting over exchange
and dipolar interactions. This leads to an ef-
fective pseudospin-1/2 Hamiltonian with only
the pseudospin components along the local easy
axes present ( 43 ). These pseudospin components
are thus good quantum numbers, justifying the
use of the classical Ising Hamiltonian.
In HoAgGe, metallicity simultaneously sup-
presses the CEF splitting of Ho3+ions and en-
hances the exchange coupling between them,
making the two energy scales comparable at
least for the low-lying CEF levels. Thus, the
large (J=8)Ho3+moments in HoAgGe at
moderately low temperatures, when multiple
CEF levels are occupied, are closer to semi-
classical spins with strongsingle-ion anisotropy.
Such a semiclassical model can still be mapped
to an Ising model at the expense of introducing
further-neighbor exchange interactions ( 44 ),
which serves as an explanation for the apparent
validity of the classical Ising Hamiltonian for
HoAgGe. A complete understanding of the
entropy data awaits a full quantum mechanical
description of the system. It is worth noting
that the deviation from an ideal spin-1/2 system
can also lead to stronger quantum fluctuations, as
in the cases of Tb 2 Ti 2 O 7 ( 45 )andTb 2 Sn 2 O 7 ( 46 ).
The metallic nature of HoAgGe not only
makes it a high-temperature (compared with
pyrochlore spin ices) kagome ice, but may also
lead to exotic phenomena such as the interac-
tion between electric currents and the magnetic
monopoles or the toroidal moments, the rela-
tionship between the noncollinear ordering
and the anomalous Hall effect ( 47 – 50 ), and
metallic magnetoelectric effects caused by
broken inversion symmetry ( 51 ). Our results
suggest that ZrNiAl-type intermetallic com-
pounds are a prototypical family of kagome
spin systems that may host other exotic phases
beyond the classical spin liquid ( 52 , 53 )and
deserve further investigation.
REFERENCES AND NOTES
- L. Balents,Nature 464 , 199–208 (2010).
- M. J. Harris, S. T. Bramwell, D. F. McMorrow, T. Zeiske,
K. W. Godfrey,Phys. Rev. Lett. 79 , 2554–2557 (1997). - A. P. Ramirez, A. Hayashi, R. J. Cava, R. Siddharthan,
B. S. Shastry,Nature 399 , 333–335 (1999). - S. T. Bramwell, M. J. P. Gingras,Science 294 , 1495– 1501
(2001). - C. Castelnovo, R. Moessner, S. L. Sondhi,Nature 451 ,42– 45
(2008). - L. Pauling,J. Am. Chem. Soc. 57 , 2680–2684 (1935).
- D. J. P. Morriset al.,Science 326 , 411–414 (2009).
- T. Fennellet al.,Science 326 , 415–417 (2009).
- A. S. Wills, R. Ballou, C. Lacroix,Phys. Rev. B 66 , 144407
(2002). - G. Möller, R. Moessner,Phys. Rev. B 80 , 140409 (2009).
- G.-W. Chern, P. Mellado, O. Tchernyshyov,Phys. Rev. Lett. 106 ,
207202 (2011). - R. F. Wanget al.,Nature 439 , 303–306 (2006).
- Y. Qi, T. Brintlinger, J. Cumings,Phys. Rev. B 77 , 094418
(2008). - S. Ladak, D. Read, G. K. Perkins, L. F. Cohen, W. Branford,
Nat. Phys. 6 , 359–363 (2010). - E. Mengottiet al.,Nat. Phys. 7 ,68–74 (2011).
- S. Zhanget al.,Nature 500 , 553– 557 (2013).
17.C. Nisoli, R. Moessner, P. Schiffer,Rev. Mod. Phys. 85 ,
1473 – 1490 (2013). - J. C. Gartsideet al.,Nat. Nanotechnol. 13 ,53–58 (2018).
- K. Matsuhira, Z. Hiroi, T. Tayama, S. Takagi, T. Sakakibara,
J. Phys. Condens. Matter 14 , L559–L565 (2002). - Y. Tabataet al.,Phys. Rev. Lett. 97 , 257205 (2006).
- T. Fennell, S. T. Bramwell, D. F. McMorrow, P. Manuel,
A. R. Wildes,Nat. Phys. 3 , 566–572 (2007). - Z. L. Dunet al.,Phys. Rev. Lett. 116 , 157201 (2016).
- J. A. M. Paddisonet al.,Nat. Commun. 7 , 13842 (2016).
- E. Lhotelet al.,Phys. Rev. Lett. 115 , 197202 (2015).
- E. Lhotelet al.,Nat. Commun. 9 , 3786 (2018).
- R. Gibson, R. Pöttgen, R. K. Kremer, A. Simon, K. R. A. Ziebeck,
J. Alloys Compd. 239 ,34–40 (1996). - E. Morosan, S. L. Bud’ko, P. C. Canfield, M. S. Torikachvili,
A. H. Lacerda,J. Magn. Magn. Mater. 277 , 298–321 (2004). - H. T. Stokes, D. M. Hatch,J. Appl. Cryst. 38 , 237– 238
(2005). - S. Baranet al.,J. Alloys Compd. 281 ,92–98 (1998).
- E. Morosan, S. L. Bud’ko, P. C. Canfield,Phys. Rev. B Condens.
Matter Mater. Phys. 71 , 014445 (2005). - Materials and methods are available as supplementary
materials. - V. M. Dubovik, V. V. Tugushev,Phys. Rep. 187 ,1 45 – 202
(1990). - A. Dönniet al.,J. Phys. Condens. Matter 8 , 11213–11229 (1998).
- S. A. M. Mentinket al.,Phys. Rev. Lett. 73 , 1031–1034 (1994).
- T. Takagi, M. Mekata,J. Phys. Soc. Jpn. 62 , 3943– 3953
(1993). - M. Wolf, K. D. Schotte,J. Phys. Math. Gen. 21 , 2195– 2209
(1988). - G.-W. Chern, O. Tchernyshyov,Phil. Trans. R. Soc. A 370 ,
5718 – 5737 (2012). - E. D. Mun, S. L. Bud’ko, H. Ko, G. J. Miller, P. C. Canfield,
J. Magn. Magn. Mater. 322 , 3527–3543 (2010). - D. I. Gorbunovet al.,Phys. Rev. B 97 , 184412 (2018).
40. M. Russinaet al.,Physica B 551 , 506–511 (2018).
41. P. Javorský, M. Diviš, H. Sugawara, H. Sato, H. Mutka,
Phys. Rev. B 65 , 014404 (2001).
42. A. Scheieet al.,Phys. Rev. B 93 , 180407 (2016).
43. J. G. Rau, M. J. P. Gingras,Phys. Rev. B 92 ,144417(2015).
44. Y.-Z. Chou, Y.-J. Kao,Phys. Rev. B 82 , 132403 (2010).
45. I. Mirebeauet al.,Phys. Rev. Lett. 94 , 246402 (2005).
46. M. J. P. Gingraset al.,Phys. Rev. B 62 ,6496–6511 (2000).
47. R. Shindou, N. Nagaosa,Phys.Rev.Lett. 87 , 116 801 (2001).
48. Y. Taguchi, Y. Oohara, H. Yoshizawa, N. Nagaosa, Y. Tokura,
Science 291 , 2573–2576 (2001).
49. H. Chen, Q. Niu, A. H. MacDonald,Phys. Rev. Lett. 112 , 017205
(2014).
50. S. Nakatsuji, N. Kiyohara, T. Higo,Nature 527 ,212– 215
(2015).
51. H. Saitoet al.,J. Phys. Soc. Jpn. 87 , 033702 (2018).
52. J. Carrasquilla, Z. Hao, R. G. Melko,Nat. Commun. 6 , 7421
(2015).
53. Y. Huang, M. Hermele,Phys. Rev. B 95 , 075130 (2017).
54. K. Zhao, H. Deng, H. Chen, K. A. Ross, V. Petříček, G. Günther,
M. Russina, V. Hutanu, P. Gegenwart, Data for: Realization of
the kagome spin ice state in a frustrated intermetallic
compound, Zenodo (2020); https://doi.org/10.5281/zenodo.
3623256.
ACKNOWLEDGEMENTS
We thank Y. Tokiwa, S. Nakatsuji, O. Tchernyshyov, J. Xu,
J. Song, V. Peçanha-Antonio, and M. Meven for helpful
discussions; A. Jesche for providing access to the MPMS;
and T. Delazzer for help with the CEF calculations. K.A.R.
acknowledges J. Rau and M. Gingras for CEF discussions. We
also thank K. Kiefer and colleagues from Helmholz Zentrum
Berlin for providing access and assistance using the VM4
magnet at MLZ. The technical assistance of W. Luberstetter
in setting up the Oxford Instruments magnet on POLI is
appreciated.Funding:The work in Augsburg was supported
by the German Science Foundation through SPP1666 (project
no. 220179758) and TRR80 (project no. 107745057). The
instrument POLI at Heinz Maier-Leibnitz Zentrum (MLZ), Garching,
Germany, was operated by RWTH Aachen University in cooperation
with FZ Jülich (Jülich Aachen Research Alliance JARA). Inelastic
neutron experiments were conducted at the time-of-flight
spectrometer NEAT operated by Helmholtz Zentrum Berlin. This
work utilized the RMACC Summit supercomputer, which is
supported by the National Science Foundation (awards
ACI-1532235 and ACI-1532236), the University of Colorado–Boulder,
and Colorado State University. The Summit supercomputer is
a joint effort of the University of Colorado–Boulder and Colorado
State University. The work in Prague was supported by the Czech
Science Foundation through project no. 18-10504S.Author
contributions:K.Z. and P.G. proposed the experiments; K.Z.
synthesized single crystals and measured magnetic properties
and specific heat; H.D. and V.H. conducted single-crystal elastic
neutron scattering; V.P., H.D., and K.Z. refined magnetic structures;
G.G. and M.R. conducted inelastic neutron scattering; H.C.
provided theoretical analysis and MC simulation; K.A.R. performed
CEF calculations; K.Z., H.C., and P.G. wrote the manuscript with
input from all authors.Competing interests:The authors declare
no competing interests.Data and materials availability:The
data presented in this paper can be found on Zenodo ( 54 ).
SUPPLEMENTARY MATERIALS
science.sciencemag.org/content/367/6483/1218/suppl/DC1
Materials and Methods
Supplementary Text
Figs. S1 to S22
Tables S1 to S4
References ( 55 – 65 )
26 November 2018; accepted 14 February 2020
10.1126/science.aaw1666
Zhaoet al.,Science 367 , 1218–1223 (2020) 13 March 2020 6of6
RESEARCH | RESEARCH ARTICLE