Pith. sign in

REVIEW 4 major objections 4 minor 2 cited by

Ce2Hf2O7 magnetization and magnetostriction data are reproduced by a spin Hamiltonian with dominant octupole exchange, indicating the emergence of a dipole-octupole quantum spin ice.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Magnetization and magnetostriction data on Ce2Hf2O7 show a kink anomaly best reproduced by a Hamiltonian with dominant octupole interactions, supporting an octupole-based quantum spin ice ground state.

T0 review reviewed 2026-08-04 challenge →

load-bearing objection Solid new low-temperature magnetization and magnetostriction data on Ce2Hf2O7; the QSI conclusion, however, hangs on a 16-site ED cluster that disagrees with classical MC and is never checked for AIAO compatibility. the 4 major comments →

arxiv 2509.09189 v1 pith:66MGQFG2 submitted 2025-09-11 cond-mat.str-el

Magnetization and magnetostriction measurements of the dipole-octupole quantum spin ice candidate Ce2Hf2O7

classification cond-mat.str-el
keywords quantum spin icedipole-octupoleoctupole interactionpyrochloremagnetizationmagnetostrictionCe2Hf2O7exact diagonalization
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper reports magnetization and magnetostriction measurements on the pyrochlore compound Ce2Hf2O7 down to 50 mK. It finds magnetic hysteresis below around 300 mK and a kink-like anomaly in the magnetization under B || [111], with a matching feature in the magnetostriction. The authors argue that these observations are quantitatively reproduced by a nearest-neighbor spin Hamiltonian whose dominant interaction couples octupole moments, with a U(1)_pi quantum spin ice ground state. If true, this would establish Ce2Hf2O7 as a concrete realization of an octupole-based quantum spin liquid, where fractionalized excitations emerge from octupolar degrees of freedom rather than the usual dipolar spins.

Core claim

The central claim is that the low-temperature magnetization and magnetostriction of Ce2Hf2O7 are consistent with a dipole-octupole pseudospin-1/2 Hamiltonian in which the dominant exchange is octupolar and the ground state is a U(1)_pi quantum spin ice. Classical Monte Carlo and exact diagonalization calculations, using exchange parameters inherited from prior heat-capacity fits, reproduce the kink anomaly under B || [111] only when the leading interaction is octupolar; dipolar-dominated parameters instead yield a conventional plateau. The magnetostriction along [111] shows a convex deviation at the same field, fitted with the same Hamiltonian after adjusting magnetoelastic couplings. The ma

What carries the argument

The key object is the nearest-neighbor dipole-octupole exchange Hamiltonian for the pyrochlore lattice, with coupling parameters (Jx, Jy, Jz, Jxz) acting on Ce3+ pseudospins. A global rotation reduces this Hamiltonian to an XYZ model; the rotation angle theta tunes the balance between dipolar and octupolar exchange. The paper combines classical Monte Carlo on a 6912-spin lattice and exact diagonalization on a 16-site cluster to compute magnetization curves, scanning theta and exchange permutations to fit the experimental data. Magnetostriction is modeled by coupling average spin components to elastic strain, with magnetoelastic constants obtained by fitting the measured response.

Load-bearing premise

The analysis assumes that Ce2Hf2O7 is described by the nearest-neighbor dipole-octupole Hamiltonian within the two exchange-parameter sets taken from heat-capacity fits; if further-neighbor terms or entirely different exchange regimes apply, the conclusion of an octupole-dominated U(1)_pi quantum spin ice ground state would not follow.

What would settle it

A measurement that would settle the claim: low-temperature inelastic neutron scattering on Ce2Hf2O7, looking for the gapless photon-like continuum expected for a U(1)_pi quantum spin ice; alternatively, observation of a static ordered moment via muon spin rotation or neutron diffraction below 50 mK would directly contradict the quantum spin liquid ground state.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • If the interpretation is correct, Ce2Hf2O7 becomes a concrete example of a U(1)_pi octupolar quantum spin ice, adding to the short list of materials where octupole moments carry the fractionalized degrees of freedom.
  • The kink anomaly under B || [111] serves as an experimental fingerprint of dominant octupolar exchange, distinguishing this state from conventional dipolar spin ice in field-dependent measurements.
  • The magnetostriction convexity provides an independent thermodynamic corroboration of the same field-induced transition, reinforcing the magnetization data.
  • The slow relaxation and hysteresis below 300 mK imply that the low-field spin-liquid state has a strongly suppressed population of mobile excitations, which should be reflected in ac susceptibility and thermal conductivity measurements.
  • The fitted parameter sets give concrete predictions for future neutron scattering, resonant X-ray, or NMR experiments on Ce2Hf2O7, guiding searches for the predicted photon-like and spinon excitations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same fitting strategy could be applied to other dipole-octupole pyrochlores such as Ce2Zr2O7 and Ce2Sn2O7 to test whether a kink-like anomaly is a universal signature of octupole-dominant exchange or is specific to this compound.
  • The study leaves open the possibility that further-neighbor exchange, not included in the Hamiltonian, could shift the optimum parameters; a systematic extension to second-neighbor couplings would sharpen the assignment of the ground state.
  • The MC- and ED-based optimal parameters differ slightly, indicating combined quantum-classical and finite-size systematics; direct measurement of the octupolar structure factor would provide a cleaner test than magnetization alone.
  • If the hysteresis indeed reflects an excitation gap, one expects the kink to shift or vanish with changing field-sweep rate; a systematic sweep-rate study could separate equilibrium thermodynamics from glassy dynamics.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper reports field-dependent magnetization measurements down to 50 mK for B along [100], [110], and [111], plus magnetostriction along [111] under B || [111], on the pyrochlore Ce2Hf2O7. It finds magnetic hysteresis below roughly 300 mK and a kink-like feature in M(B) under B || [111] that is also reflected in the magnetostriction. Using a nearest-neighbor dipole-octupole Hamiltonian with exchange parameters taken from prior heat-capacity fits (Sets A and B), the authors perform classical Monte Carlo and 16-site exact diagonalization and claim that the magnetization curves are well reproduced by parameters with dominant octupolar interactions and a U(1)_pi quantum spin ice ground state, indicating a dipole-octupole QSI.

Significance. The experimental data are careful and useful: they extend to 50 mK, are reproducible on two samples, and the coincidence of the magnetization kink with a convex magnetostriction feature is a valuable phenomenological constraint. A confirmed dipole-octupole QSI in Ce2Hf2O7 would be significant for frustrated magnetism. The authors also make a good faith effort to tie their calculations to previously constrained parameter sets and to use two complementary numerical methods. However, the central QSI assignment is not as secure as the abstract suggests, because the decisive calculations rely on a small ED cluster with admitted finite-size bias, and the level of agreement is largely an in-sample fit.

major comments (4)
  1. [SM S5, Table 1] The distinction between QSI (Set A) and AIAO (Set B) is not settled by the calculations as presented. The classical MC minimum deviation is for Set B (a,c,b) with 0.137, whereas the ED minimum is for Set A (a,b,c) with 0.848 (best Set B: 1.167). The manuscript argues that ED is more reliable for quantum spin liquid physics, but SM S5 states that ED's 'principal systematic error is the finite-size effect and boundary condition bias, which can imprint discrete-level features and shift optimal parameters on small clusters.' The paper does not quantify this bias or show it is smaller than the ED misfit gap of ~0.32 between sets. The specific concern that the 16-site cluster cannot accommodate q=0 AIAO is not literally correct, since a 16-site cubic cell can represent q=0 order; nevertheless, the cluster may still bias the relative stability of QSI vs AIAO. A cluster-size scaling, DMRG check,
  2. [Fig. 4 and SM S5] The 'good agreement' is an in-sample fit rather than an ab initio reproduction. The theta minima are selected by minimizing the deviation from the measured magnetization (SM S5: 'We then investigated the theta dependence of the sum of the deviation...'), and the magnetostriction curve is obtained after six magnetoelastic coupling combinations are fit to the same data ('determined by fitting the experimental data'). The kink under B || [111] is therefore not a falsifiable prediction. The paper should either fix all parameters from independent data and then compare with M(B) and magnetostriction, or explicitly reframe the analysis as parameter extraction with covariance/uncertainty estimates and out-of-sample checks (e.g., leave out one field direction or temperature and predict it).
  3. [Table 1] The labels 'QSI' and 'AIAO' are inherited from Sets A and B of Ref. [16], but the fitted parameters in Table 1 are not the original triples; they are obtained after a theta rotation and a permutation. The manuscript never shows where the final (Jx,Jy,Jz,Jxz) values lie in the phase diagram of Ref. [6]. For example, Set A (a,b,c) at theta=0.14 pi gives Jx=0.042, Jy=0.021, Jz=0.012, Jxz=0.018 meV; it is not automatic that this point remains in the U(1)_pi QSI region. If any of the Table 1 points fall outside the claimed phases, the correlation between 'best fit' and 'QSI/AIAO' is broken. Please plot the fitted points on the phase diagram or provide the phase-boundary values.
  4. [Main text Discussion] The conclusion is stated more strongly than the paper's own caveats. The abstract says the calculations 'demonstrate' a QSI ground state, and the Discussion calls it 'direct evidence,' while the main text also concedes that 'more decisive distinction between the candidate parameters requires further experimental investigation' and SM S5 says 'identifying the interaction parameters ... requires further experiments, which remains a future issue.' The central claim should be softened to 'consistent with' or 'suggestive of' unless the robustness checks requested above are provided.
minor comments (4)
  1. [SM S5] Typo: 'magnetostraiction' should be 'magnetostriction'.
  2. [SM S5] SM S5 refers to 'Fig. 5 in the main text,' but the relevant figure with the theta dependence is Fig. 4(b) in the main text.
  3. [Fig. 3] The figure caption does not explicitly state that the ED calculations are equilibrium calculations at 50 mK, whereas the experimental demagnetization curves may be influenced by the long relaxation and hysteresis described in the text. Please clarify what the comparison means in a non-equilibrium setting.
  4. [References] Ref. [23] is an arXiv preprint; if a published version is available, it should be cited instead.

Circularity Check

2 steps flagged

The magnetization and magnetostriction 'reproductions' are in-sample fits (θ and g_i optimized against the same data); partial independent support from the prior kink prediction in Ref. [23] keeps this from being fully forced.

specific steps
  1. fitted input called prediction [Discussion (main text, after Fig. 3); Table 1; SM S5]
    "These parameters arise as the best-fit when both our MC and ED calculations (discussed below) are considered. As shown in Fig. 3(b), the calculated magnetization curves well reproduce our experimental results, in particular the kink-like anomaly observed under B || [111]."

    The 'best-fit' parameters are obtained by varying θ to minimize Σ|M_MC − M_exp| for B ≤ 2 T, using the same 50 mK demagnetization data (B || [100] and B || [111]) that contains the kink (SM S5; Table 1 lists parameters 'corresponding to the minimum of Σ|M_MC − M_exp|'). The calculated curves at that minimum are then presented as reproducing the same data, including the kink. This is an in-sample fit reported as a reproduction, not an independent prediction; the kink is not an out-of-sample confirmation within this paper.

  2. fitted input called prediction [SM S5, magnetostriction subsection (after Eq. for ΔL/L)]
    "What remains unknown are the magnetoelastic coupling constants g_i, which are determined by fitting the experimental data as follows. ... As shown in Fig. S7, our fitting gives a good agreement with the experimental data at 100 mK and 200 mK, supporting the validity of our calculations done for the magnetization data."

    The magnetoelastic constants are fitted to the very magnetostriction curve that the calculation is then said to 'reproduce'; the agreement is therefore guaranteed in-sample. Calling the fit 'supporting the validity' of the magnetization calculations uses a fit of the same data as confirmation, rather than an independent check. This does not affect the magnetization data themselves, but it inflates the evidentiary weight of the magnetostriction match.

full rationale

The core magnetization comparison is also partly a fitting exercise: the allowed parameter space (Sets A/B) comes from the prior heat-capacity study [16], and θ is then scanned for each permutation to minimize the deviation from the new magnetization data (SM S5, Fig. 4b, Table 1). Claiming that the resulting curves 'well reproduce' the data, including the kink, is a fit-quality statement rather than a prediction. This is the main circularity and justifies score 6. It is not higher because (i) the kink under B || [111] was noted to 'appear in the theoretical calculations [23]' before measurement, providing genuine prior support; (ii) Sets A and B are constrained by an independent heat-capacity observable in [16]; and (iii) the model could in principle have failed to reproduce the data even at the best θ. The magnetostriction match is more clearly in-sample, with g_i explicitly fitted to the same data. Self-citations to Refs [16,23] exist, but they are not load-bearing in a circular way: [16] is a different observable and [23] is a prior prediction. The paper itself flags the main non-circular concern: ED 'principal systematic error is the finite-size effect and boundary condition bias, which can imprint discrete-level features and shift optimal parameters on small clusters' (SM S5), and it concedes 'more decisive distinction between the candidate parameters requires further experimental investigation.' Because classical MC actually favors a Set B (AIAO) permutation (deviation 0.137 vs 0.318 for the chosen Set A), the QSI conclusion rests heavily on the finite-size ED method; this is a correctness/systematic risk, not a reduction of the result to its inputs, and is weighed here without raising the circularity score further.

Axiom & Free-Parameter Ledger

4 free parameters · 5 axioms · 0 invented entities

The central claim depends on a small number of fitted parameters: theta, g_z, and several magnetoelastic constants. The spin exchange parameters themselves originate from prior heat-capacity fits (Ref [16]) and are not independently derived here. No new physical entities are introduced; the octupole moments and U(1)π QSI are established concepts in the cited literature. The main model assumptions are the dipole-octupole doublet description, nearest-neighbor exchange truncation, and the validity of the numerical methods used.

free parameters (4)
  • theta (pseudospin rotation angle) = e.g., 0.14π for (Ja,Jb,Jc) permutation of Set A; values differ per permutation in Table 1
    Determined by minimizing sum |M_MC - M_exp| for B <= 2 T (SM S5, Fig. 4b).
  • g_z = 2.36
    Obtained from the saturated moment M_parallel = 1.18 mu_B fitted to high-field magnetization (SM S5).
  • Magnetoelastic coupling combinations (8√2g1 - 4g2, 8√2g7 - 4g8, g3, g4, g9, g10) = -23.93, 7.978, 11.54, -20.37, -0.3009, 0.6017
    Fitted to the measured magnetostriction at 100 mK and 200 mK (SM S5).
  • Spin exchange parameters (Jx, Jy, Jz, Jxz) = e.g., (0.042, 0.021, 0.012, 0.018) meV for best Set A permutation
    Derived from prior heat-capacity Sets A/B plus fitted theta; not independently fitted here but treated as adjustable via theta.
axioms (5)
  • domain assumption The CEF ground state of Ce3+ in Ce2Hf2O7 is a dipole-octupole Kramers doublet described by an effective spin-1/2.
    Assumed from Refs [12-16]; if the doublet assignments were wrong, the entire model Hamiltonian would not apply.
  • domain assumption The nearest-neighbor exchange Hamiltonian (Eq. 1) is sufficient to describe the low-field magnetization and magnetostriction.
    Further-neighbor terms and long-range dipolar interactions are neglected; the paper notes weak further-neighbor terms as an uncertainty (SM S5).
  • standard math The global pseudospin rotation eliminates the Jxz term without loss of generality.
    Standard transformation shown in SM S5; this is a mathematical identity for the two-site exchange.
  • domain assumption Classical Monte Carlo at 25 mK and exact diagonalization on a 16-site cluster adequately capture the experimental magnetization at 50 mK.
    The paper acknowledges both have limitations: MC neglects quantum effects, ED has finite-size bias (SM S5).
  • domain assumption The phase diagram of Benton (Ref [6]) correctly identifies Set A as U(1)π QSI and Set B as all-in-all-out ordered.
    Used to connect the fitted parameters to the ground state label; the authors rely on the prior phase diagram.

reviewed 2026-08-04 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Magnetization and magnetostriction measurements of the dipole-octupole quantum spin ice candidate Ce2Hf2O7." pith.science (2026). https://pith.science/paper/66MGQFG2

@misc{pith2026250909189,
  author       = {Pith},
  title        = {Pith review of: Magnetization and magnetostriction measurements of the dipole-octupole quantum spin ice candidate Ce2Hf2O7},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/66MGQFG2}},
  note         = {Machine review of arXiv:2509.09189}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

We investigate the magnetization and the magnetostriction of the dipole-octupole quantum spin ice candidate Ce2Hf2O7 down to 50 mK. We find that the magnetization curves observed with the magnetic field applied along all the principal axes ([100], [110], and [111]) exhibit a magnetic hysteresis below around 300 mK. In addition, a kink-like feature is observed in the magnetization under B || [111], at which the magnetostriction also shows a convex field dependence. Our classical Monte-Carlo and quantum exact diagonalization calculations demonstrate that these features in the magnetization are well reproduced by the spin Hamiltonian with a dominant interaction between the octupole moments and with a QSI ground state, indicating the emergence of a dipole-octupole QSI in this compound.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Observation of Dipolar Spin-ice--like Correlations in the Quantum Spin Ice Candidate Ce$_2$Sn$_2$O$_7$

    cond-mat.str-el 2026-01 conditional novelty 7.0

    In Ce2Sn2O7, low-temperature diffuse neutron scattering matches Dy2Ti2O7 dipolar spin ice and shows no all-in-all-out order, falsifying the nearest-neighbor XYZ model.

  2. Quantum Fisher Information as a Thermal Probe in Frustrated Magnets through Insights from Quantum Spin Ice

    cond-mat.str-el 2025-10 conditional novelty 7.0

    Quantum Fisher information computed for the pyrochlore quantum-spin-ice model maps its phase diagram and thermal crossover scales, establishing it as a thermal/dynamical probe accessible to neutron scattering.

Reference graph

Works this paper leans on

40 extracted references · 1 linked inside Pith · cited by 2 Pith papers

  1. [1]

    M. J. P. Gingras and P. A. McClarty, Quantum spin ice: a search for gapless quantum spin liquids in pyrochlore magnets, Rep. Prog. Phys. 77, 056501 (2014)

  2. [2]

    Castelnovo, R

    C. Castelnovo, R. Moessner, and S. L. Sondhi, Magnetic monopoles in spin ice, Nature 451, 42 (2008)

  3. [3]

    Hermele, M

    M. Hermele, M. P. A. Fisher, and L. Balents, Pyrochlore photons: The U(1) spin liquid in a S=1/2 three-dimensional frustrated magnet, Phys. Rev. B 69, 064404 (2004)

  4. [4]

    J. G. Rau and M. J. P. Gingras, Frustrated Quantum Rare-Earth Pyrochlores, Annu. Rev. Condens. Matter Phys. 10, 357 (2019)

  5. [6]

    Benton, Ground-state phase diagram of dipolar-octupolar pyrochlores, Phys

    O. Benton, Ground-state phase diagram of dipolar-octupolar pyrochlores, Phys. Rev. B 102, 104408 (2020)

  6. [8]

    Desrochers, L

    F. Desrochers, L. E. Chern, and Y. B. Kim, Competing U (1) and Z 2 dipolar-octupolar quantum spin liquids on the pyrochlore lattice: Application to Ce2Zr2O7, Phys. Rev. B 105, 035149 (2022)

  7. [9]

    Hosoi, E

    M. Hosoi, E. Z. Zhang, A. S. Patri, and Y. B. Kim, Uncovering Footprints of Dipolar- Octupolar Quantum Spin Ice from Neutron Scattering Signatures, Phys. Rev. Lett. 129, 097202 (2022)

  8. [10]

    Sanders, H

    A. Sanders, H. Yan, C. Castelnovo, and A. H. Nevidomskyy, Experimentally Tunable QED in Dipolar-Octupolar Quantum Spin Ice, arXiv:2312.11641

  9. [11]

    E. M. Smith, E. Lhotel, S. Petit, and B. D. Gaulin, Experimental Insights into Quantum Spin Ice Physics in Dipole– Octupole Pyrochlore Magnets, Annu. Rev. Condens. Matter Phys. 16, 387 (2025)

  10. [12]

    Sibille, E

    R. Sibille, E. Lhotel, V. Pomjakushin, C. Baines, T. Fennell, and M. Kenzelmann, Candidate Quantum Spin Liquid in the Ce 3+ Pyrochlore Stannate Ce2Sn2O7, Phys. Rev. Lett. 115, 097202 (2015)

  11. [13]

    Gaudet et al., Quantum Spin Ice Dynamics in the Dipole-Octupole Pyrochlore Magnet Ce2Zr2O7, Phys

    J. Gaudet et al., Quantum Spin Ice Dynamics in the Dipole-Octupole Pyrochlore Magnet Ce2Zr2O7, Phys. Rev. Lett. 122, 187201 (2019)

  12. [14]

    Gao et al., Experimental signatures of a three- dimensional quantum spin liquid in 9 effective spin-1/2 Ce2Zr2O7 pyrochlore, Nat

    B. Gao et al., Experimental signatures of a three- dimensional quantum spin liquid in 9 effective spin-1/2 Ce2Zr2O7 pyrochlore, Nat. Phys. 15, 10 (2019)

  13. [15]

    Porée et al., Crystal -field states and defect levels in candidate quantum spin ice Ce2Hf2O7, Phys

    V. Porée et al., Crystal -field states and defect levels in candidate quantum spin ice Ce2Hf2O7, Phys. Rev. Mater. 6, 044406 (2022)

  14. [17]

    Sibille et al., A quantum liquid of magnetic octupoles on the pyrochlore lattice, Nat

    R. Sibille et al., A quantum liquid of magnetic octupoles on the pyrochlore lattice, Nat. Phys. 16, 5 (2020)

  15. [18]

    Porée, H

    V. Porée, H. Yan, F. Desrochers, S. Petit, E. Lhotel, M. Appel, J. Ollivier, Y. B. Kim, A. H. Nevidomskyy, and R. Sibille, Evidence for fractional matter coupled to an emergent gauge field in a quantum spin ice, Nat. Phys. 21, 83 (2024)

  16. [19]

    D. R. Yahne et al., Dipolar Spin Ice Regime Proximate to an All -In-All-Out Néel Ground State in the Dipolar-Octupolar Pyrochlore Ce 2Sn2O7, Phys. Rev. X 14, 011005 (2024)

  17. [20]

    E. M. Smith et al., Case for a U(1) π Quantum Spin Liquid Ground State in the Dipole - Octupole Pyrochlore Ce2Zr2O7, Phys. Rev. X 12, 021015 (2022)

  18. [21]

    Bhardwaj, S

    A. Bhardwaj, S. Zhang, H. Yan, R. Moessner, A. H. Nevidomskyy, and H. J. Changlani, Sleuthing out exotic quantum spin liquidity in the pyrochlore magnet Ce 2Zr2O7, Npj Quantum Mater. 7, 1 (2022)

  19. [22]

    E. M. Smith et al., Quantum spin ice response to a magnetic field in the dipole- octupole pyrochlore Ce2Zr2O7, Phys. Rev. B 108, 054438 (2023)

  20. [23]

    Porée et al., Dipolar-Octupolar Correlations and Hierarchy of Exchange Interactions in Ce2Hf2O7, arXiv:2305.08261

    V. Porée et al., Dipolar-Octupolar Correlations and Hierarchy of Exchange Interactions in Ce2Hf2O7, arXiv:2305.08261

  21. [24]

    Bhardwaj, V

    A. Bhardwaj, V. Porée, H. Yan, N. Gauthier, E. Lhotel, S. Petit, J. A. Quilliam, A. H. Nevidomskyy, R. Sibille, and H. J. Changlani, Thermodynamics of the dipole- octupole pyrochlore magnet Ce2Hf2O7 in applied magnetic fields, Phys. Rev. B 111, 155137 (2025)

  22. [25]

    Fukazawa, R

    H. Fukazawa, R. G. Melko, R. Higashinaka, Y. Maeno, and M. J. P. Gingras, Magnetic anisotropy of the spin-ice compound Dy2Ti2O7, Phys. Rev. B 65, 054410 (2002)

  23. [26]

    Sakakibara, T

    T. Sakakibara, T. Tayama, Z. Hiroi, K. Matsuhira, and S. Takagi, Observation of a Liquid- Gas-Type Transition in the Pyrochlore Spin Ice Compound Dy 2Ti2O7 in a Magnetic Field, Phys. Rev. Lett. 90, 207205 (2003)

  24. [27]

    Placke, R

    B. Placke, R. Moessner, and O. Benton, Hierarchy of energy scales and field- tunable order by disorder in dipolar-octupolar pyrochlores, Phys. Rev. B 102, 245102 (2020)

  25. [28]

    Snyder, B

    J. Snyder, B. G. Ueland, J. S. Slusky, H. Karunadasa, R. J. Cava, and P. Schiffer, Low- temperature spin freezing in the Dy 2Ti2O7 spin ice, Phys. Rev. B 69, 064414 (2004)

  26. [29]

    Matsuhira, C

    K. Matsuhira, C. Paulsen, E. Lhotel, C. Sekine, Z. Hiroi, and S. Takagi, Spin Dynamics at Very Low Temperature in Spin Ice Dy 2Ti2O7, J. Phys. Soc. Jpn. 80, 123711 (2011). 10

  27. [30]

    S. R. Giblin, S. T. Bramwell, P. C. W. Holdsworth, D. Prabhakaran, and I. Terry, Creation and measurement of long- lived magnetic monopole currents in spin ice, Nat. Phys. 7 , 252 (2011). [31]D. Pomaranski, L. R. Yaraskavitch, S. Meng, K. A. Ross, H. M. L. Noad, H. A. Dabkowska, B. D. Gaulin, and J. B. Kycia, Absence of Pauling’s residual entropy in therm...

  28. [32]

    Paulsen, M

    C. Paulsen, M. J. Jackson, E. Lhotel, B. Canals, D. Prabhakaran, K. Matsuhira, S. R. Giblin, and S. T. Bramwell, Far- from-equilibrium monopole dynamics in spin ice, Nat. Phys. 10, 135 (2014)

  29. [33]

    E. R. Kassner, A. B. Eyvazov, B. Pichler, T. J. S. Munsie, H. A. Dabkowska, G. M. Luke, and J. C. S. Davis, Supercooled spin liquid state in the frustrated pyrochlore Dy 2Ti2O7, Proc. Natl. Acad. Sci. 112, 8549 (2015). Acknowledgements: We thank Cristian D. Batista and Y. B. Kim for fruitful discussions. This work was supported by Grants-in-Aid for Scient...

  30. [34]

    E. M. Smith et al., Two -Peak Heat Capacity Accounts for Rln(2) Entropy and Ground State Access in the Dipole-Octupole Pyrochlore Ce2Hf2O7, Phys. Rev. Lett. 135, 086702 (2025)

  31. [35]

    Shimizu, Y

    Y. Shimizu, Y. Kono, T. Sugiyama, S. Kittaka, Y. Shimura, A. Miyake, D. Aoki, and T. Sakakibara, Development of high-resolution capacitive Faraday magnetometers for sub- Kelvin region, Rev. Sci. Instrum. 92, 123908 (2021)

  32. [36]

    Küchler, T

    R. Küchler, T. Bauer, M. Brando, and F. Steglich, A compact and miniaturized high resolution capacitance dilatometer for measuring thermal expansion and magnetostriction, Rev. Sci. Instrum. 83, 095102 (2012)

  33. [37]

    Doerr, T

    M. Doerr, T. Stöter, M. Rotter, and A. A. Zvyagin, Magnetostriction of the spin-ice system Yb2Ti2O7, J. Magn. Magn. Mater. 449, 378 (2018)

  34. [38]

    Stöter, M

    T. Stöter, M. Doerr, S. Granovsky, M. Rotter, S. T. B. Goennenwein, S. Zherlitsyn, O. A. Petrenko, G. Balakrishnan, H. D. Zhou, and J. Wosnitza, Extremely slow nonequilibrium monopole dynamics in classical spin ice, Phys. Rev. B 101, 224416 (2020)

  35. [39]

    Tang et al., Crystal field magnetostriction of spin ice under ultrahigh magnetic fields, Phys

    N. Tang et al., Crystal field magnetostriction of spin ice under ultrahigh magnetic fields, Phys. Rev. B 110, 214414 (2024)

  36. [40]

    A. S. Patri, M. Hosoi, and Y. B. Kim, Distinguishing dipolar and octupolar quantum spin ices using contrasting magnetostriction signatures, Phys. Rev. Res. 2, 023253 (2020)

  37. [41]

    Huang, G

    Y.-P. Huang, G. Chen, and M. Hermele, Quantum Spin Ices and Topological Phases from Dipolar-Octupolar Doublets on the Pyrochlore Lattice, Phys. Rev. Lett. 112, 167203 (2014)

  38. [42]

    Bhardwaj, V

    A. Bhardwaj, V. Porée, H. Yan, N. Gauthier, E. Lhotel, S. Petit, J. A. Quilliam, A. H. Nevidomskyy, R. Sibille, and H. J. Changlani, Thermodynamics of the dipole- octupole pyrochlore magnet Ce 2Hf2O7 in applied magnetic fields, Phys. Rev. B 111, 155137 (2025)

  39. [43]

    A. S. Patri, M. Hosoi, S. Lee, and Y. B. Kim, Theory of magnetostriction for multipolar quantum spin ice in pyrochlore materials, Phys. Rev. Res. 2, 033015 (2020)

  40. [100]

    The data obtained in the magnetization and demagnetization processes are shown by filled and open circles, respectively

    (c). The data obtained in the magnetization and demagnetization processes are shown by filled and open circles, respectively. The data is vertically shifted for clarity. 13 Fig. 3 The magnetic field dependence of the magnetization (a,b) and its field derivative (c,d) at 50 mK observed in the magnetization (a,c) and the demagnetization (b ,d) process. The ...

This paper was first reviewed by deepseek-v4-flash on August 4, 2026.