{"id":"a2225fe8-2c61-4fca-886d-2dd9bd780625","arxiv_id":"1908.03530","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A reaction-field model with no fitted Hamiltonian parameters reproduces the magnetic excitation spectrum of Ho3Ga5O12 and attributes its magnetic ordering to nuclear hyperfine coupling.","lead":"This paper models published neutron-scattering data on the garnet Ho3Ga5O12 with a reaction-field model and concludes that the nuclear hyperfine interaction, not geometrical frustration, drives magnetic ordering. It matters because it provides a framework for interpreting two-singlet crystal-field magnets, including candidate quantum spin liquids such as Pr-based pyrochlores.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The A=0 counterfactual rests on an uncontrolled reaction-field approximation; Rc=0.71 should be checked by an unbiased QMC simulation of the equivalent transverse-field Ising model.","rationale":"The reader identified the reaction-field approximation as an uncontrolled mean-field correction, and I agree that this is the most load-bearing point. The A=0 model maps exactly to a transverse-field Ising model, so the reaction-field prediction Rc=0.71 can be checked by a sign-problem-free QMC calculation. If QMC finds no order, the paper's central mechanism is strongly supported; if QMC finds order, the central claim loses its basis. This does not change the reader's conditional verdict, because the condition is now concrete: perform the unbiased ground-state calculation before accepting the hyperfine-driven-ordering conclusion. I did not elevate exchange or quadrupole omissions to the primary concern because the experimental constraints on those terms are less directly tied to the counterfactual, and the reaction-field issue alone is sufficient to motivate the same conditional stance.","tokens_in":12656,"tokens_out":18275,"duration_ms":205017,"concrete_test":"Run a sign-problem-free stochastic series expansion quantum Monte Carlo simulation of the effective spin-1/2 transverse-field Ising model on the Ho3Ga5O12 lattice with the same dipolar interaction tensor K_dip (Ewald summed), Delta=7.4 K, and A=0. Extrapolate to T=0 and to the thermodynamic limit using finite-size cells of increasing linear size (e.g., L=2,4,6 conventional cells), and compute the squared order parameter at k=(0,0) corresponding to the experimental magnetic structure. If the ground state orders, Rc=0.71 is an artifact of the reaction-field approximation and the conclusion that hyperfine drives ordering fails; if the ground state remains a quantum paramagnet, the counterfactual is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that nuclear hyperfine coupling drives magnetic ordering depends on the counterfactual that with A=0 the electronic system remains a correlated paramagnet to T=0. That counterfactual is expressed as Rc=(2m^2/Delta)max[lambda_mu(Q)-lambda]=0.71<1 (Section III C). This number is not a direct ground-state calculation of the spin model; it is the instability criterion of the reaction-field susceptibility (Eqs. 8 and 21), with lambda fixed by the sum rule (Eqs. 13-14). For this N=12, long-range dipolar, two-singlet system the reaction-field approximation is uncontrolled. The sensitivity of the method is visible in the factor-of-two difference between the standard RPA ordering temperature (0.66 K) and the reaction-field value (0.31 K). Because Rc=0.71 is only about 29% below the ordering threshold, a modest approximation error could move Rc above 1, meaning the electronic system would order by itself and hyperfine coupling would not be required. The A=0 model is exactly a spin-1/2 transverse-field Ising model with dipolar couplings, so this crucial counterfactual can be tested without the reaction-field approximation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper models published inelastic neutron-scattering data on Ho3Ga5O12 using a three-parameter Hamiltonian: the two-singlet crystal-field splitting Δ = 7.4 K, the long-range dipolar coupling D = 0.0183 K, and the nuclear hyperfine coupling A = 0.039 K, all fixed from previous experiments or first-principles values. The magnetic susceptibility is computed in a reaction-field extension of RPA for the N = 12 garnet structure, with the reaction field λ determined by a total-angular-momentum sum rule. The calculations reproduce the main features of the published spectra, give TN = 0.31 K compared with the experimental 0.3 K, and yield a magnetic ground state matching the reported propagation vector. The paper's central claim is that without the hyperfine coupling the system would remain a correlated paramagnet to zero temperature (Rc = 0.71 < 1), so that hyperfine interactions drive the ordering; it also argues that geometric frustration is unimportant for Ho3Ga5O12 and that reaction-field theory improves on standard RPA.","tokens_in":12865,"tokens_out":11334,"duration_ms":112195,"significance":"If correct, the paper provides a concrete example of hyperfine-assisted magnetic ordering in a non-Kramers two-singlet system, and it demonstrates a practical reaction-field framework for two-singlet magnets. The Hamiltonian parameters are not fitted to the modeled spectra, which is a genuine strength; the ordering temperature and the soft-mode dispersion are nontrivial predictions that match experiment to good visual accuracy. The A = 0 counterfactual is falsifiable, making the central ordering claim testable. The paper is a useful benchmark for modeling singlet ground-state magnets and a caution that hyperfine couplings may matter in non-Kramers spin-liquid candidates.","major_comments":[{"comment":"The matrix-element formulas are inconsistent with the standard diagonalization. For H = hx σx + hz σz with E = sqrt(hx^2 + hz^2), the magnitudes are |<0|σz|0>| = |hz|/E, |<1|σz|1>| = |hz|/E, and |<0|σz|1>| = |hx|/E. The displayed values (hz)^2/E^2 and (hx)^2/E^2 are the squares of these magnitudes, not the magnitudes themselves, and the factor m^2 connecting σz matrix elements to the |m_ij|^2 used in Eq. (9) is absent. As written, these equations cannot be substituted into Eq. (9); if the numerical implementation used them literally, the single-ion susceptibility, the spectra in Figs. 2 and 3, and the derived TN would all be incorrect. Please state whether these are typographical errors and specify the exact expressions used in the calculation.","section":"Appendix A, Eqs. (28)-(31)"},{"comment":"The A = 0 counterfactual is computed within the reaction-field approximation, whose λ is fixed by the sum rule in Eq. (14). This approximation is uncontrolled for the present N = 12 long-range dipolar system, and its sensitivity is evidenced by the factor-of-two difference between the reaction-field TN (0.31 K) and the standard RPA TN (0.66 K) reported in Section III.D. Since Rc = 0.71 lies only about 29% below the ordering threshold, a moderate error in λ could plausibly move the system into the ordered regime, which would invalidate the central claim that hyperfine coupling is necessary for ordering. Because the A = 0 model is exactly a spin-1/2 transverse-field Ising model with dipolar couplings, this counterfactual can be tested with quantum Monte Carlo or exact diagonalization; I recommend adding such a check or otherwise quantifying the robustness of Rc to the reaction-field approximation.","section":"Section III.C, Eq. (21)"},{"comment":"The model-data comparison is made by visual inspection of Fig. 2, and the manuscript explicitly concedes that the persistent elastic diffuse scattering below TN is not reproduced. Given that the model's validation rests on this comparison, the conclusion that the Hamiltonian is 'unambiguously determined' (Section IV) is stronger than the evidence warrants. I ask for a quantitative residual analysis or, failing that, a more cautious statement of what the data establish.","section":"Sections III.B and IV"}],"minor_comments":[{"comment":"The displayed Hamiltonian contains a duplicated equals sign ('H = =').","section":"Section II.A, Eq. (1)"},{"comment":"The phrase 'reported in several other several other rare-earth garnets' contains a duplicated phrase and should be corrected.","section":"Section III.B"},{"comment":"The value (or the limiting value) of the single-ion relaxation rate Γ is not stated anywhere, although the spectral intensities in Fig. 2 depend on it; please report the value used and indicate whether it was fixed or adjusted.","section":"Section II.B, Eq. (9)"},{"comment":"The instrumental resolution convolution is stated as a Gaussian with FWHM 1.0 K, but no details of the convolution grid, background subtraction, or treatment of elastic incoherent scattering are provided; adding these details would improve reproducibility.","section":"Section III.A"},{"comment":"The abstract and introduction describe a 'three-parameter model', but the calculation also involves an overall intensity scale factor and the relaxation rate Γ; the text should clarify that these are not Hamiltonian parameters and that Δ, D, and A are not fitted.","section":"Abstract and Section II"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope and is likely to be of interest. The main technical concern is the Appendix A matrix-element typo, which must be resolved; the authors should also provide an unbiased check of the A=0 counterfactual or temper the corresponding claim. The admitted diffuse-scattering discrepancy and the unstated value of Γ should be addressed in the revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this paper is worth engaging with. It makes a real methodological contribution—extending the reaction-field approximation to N>1 magnetic ions per unit cell—and uses it to argue, with parameters all fixed from independent experiments, that nuclear hyperfine coupling drives magnetic ordering in Ho3Ga5O12. That's a concrete, testable physics claim, and the model reproduces the published neutron data and TN well enough that I'd bet the qualitative conclusion is right.\n\nWhat's genuinely new: the full excitation spectrum calculation for a two-singlet garnet with dipolar interactions, including the hyperfine term, without fitting any Hamiltonian parameter to the data being modeled. The treatment of the sum rule, the Ewald summation, and the exact diagonalization of the single-ion space including nuclear spins are all clean. The paper also gives a useful counterfactual: with A=0, the system would remain a correlated paramagnet (Rc=0.71), and with neither hyperfine nor CEF splitting it would order at 1.74 K. That pins the suppression of ordering on the two-singlet physics, not frustration, which is a nice clarification of earlier speculation.\n\nWhere it gets soft. The comparison to experiment is visual. No residual analysis, no error bars, and a free overall intensity scale. The single-ion relaxation rate Gamma is never quantified, just called small. The claim in the conclusion that the Hamiltonian is 'unambiguously determined' is stronger than the evidence supports. More importantly, the A=0 counterfactual rests entirely on the reaction-field instability criterion. The reaction-field method is an uncontrolled mean-field correction, and the factor-of-two discrepancy between its TN (0.31 K) and standard RPA (0.66 K) shows the method's sensitivity. Since Rc=0.71 is only 29% below the threshold, a modest systematic error would flip the conclusion. The paper's central claim would be much more solid if the A=0 model were checked by an unbiased simulation—say, quantum Monte Carlo for the equivalent transverse-field Ising model. That is a doable calculation.\n\nBut don't let that stop you from reading it. The paper is honest about its unexplained diffuse scattering below TN, gives a clear physical picture, and is likely to be useful for people working on Pr-based pyrochlores and other two-singlet magnets. I'd send it to peer review and ask for quantitative residuals and, if possible, the QMC check. If the authors can't do the QMC, they should at least soften 'unambiguously'.","headline":"A solid, parameter-free model of Ho3Ga5O12 that makes a plausible case for hyperfine-driven ordering; the main open question is whether the A=0 counterfactual survives an independent calculation.","tokens_in":13438,"tokens_out":3904,"would_cite":true,"duration_ms":38715,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["75.10.Jm","78.70.Nx","75.40.Gb","75.47.Lx","75.10.Dg"],"model":"deepseek-v4-flash","headline":"In Ho3Ga5O12, the nuclear hyperfine interaction—not electronic interactions alone—drives magnetic ordering; a three-parameter model matches TN=0.31 K.","keywords":["Ho3Ga5O12","two-singlet crystal-field ground state","nuclear hyperfine interaction","reaction-field approximation","transverse-field Ising model","magnetic excitation spectrum","soft-mode magnetism","garnet magnet"],"falsifier":"Compute the zero-temperature soft-mode ratio $R_c = (2m^2/\\Delta)\\max[\\lambda_\\mu - \\lambda]$ in a model that includes nearest-neighbor exchange up to the paper's own bound $|K_{nn}| = 0.1D$; if $R_c$ reaches or exceeds 1, magnetic ordering would occur without hyperfine coupling and the claim that hyperfine coupling drives the ordering would be refuted.","tokens_in":12461,"feed_emoji":"🧲","tokens_out":13290,"duration_ms":120891,"temperature":0.7,"pith_summary":"Ho$^{3+}$Ga$_5$O$_{12}$ is a garnet whose Ho$^{3+}$ ions have a two-singlet crystal-field ground state, which maps exactly onto Ising spins in a transverse magnetic field. The paper shows that a three-parameter model---crystal-field splitting $\\Delta = 7.4$ K, long-range dipolar coupling $D = 0.0183$ K, and nuclear hyperfine coupling $A = 0.039$ K---reproduces published inelastic neutron-scattering spectra at 4 K, 0.6 K, and 0.05 K. The central result is that the nuclear hyperfine interaction drives the magnetic ordering: with $A$ included the model gives $T_N = 0.31$ K, matching experiment, while with $A = 0$ the soft-mode ratio $R_c = (2m^2/\\Delta)\\max[\\lambda_\\mu-\\lambda] = 0.71 < 1$ at $T = 0$, so the system would remain a correlated paramagnet to zero temperature. The paper also finds that the suppression of $T_N$ is caused by the two-singlet crystal field rather than by geometrical frustration, and that a reaction-field correction to RPA is required to account for the spectra.","feed_headline":"Hyperfine coupling is what orders the garnet Ho3Ga5O12","feed_subtitle":"Fixed parameters reproduce the 0.31 K transition; without hyperfine coupling, the material never orders.","key_machinery":"The load-bearing object is the projected single-ion Hamiltonian $H = (\\Delta/2)\\, \\mathbb{1}\\otimes\\sigma^x + A\\, m\\, I^z\\otimes\\sigma^z$, acting in the 16-dimensional space of two electronic singlets and the eight nuclear-spin states of $^{165}$Ho, together with the long-range Ising dipolar interaction $K_{ij}$ between local $z$-components. Because the Hamiltonian commutes with $I^z$, each nuclear-spin sector is a two-level system with transverse field $h_x = \\Delta/2$ and longitudinal field $h_z = A m I^z$; the hyperfine coupling splits the two electronic singlets into eight doublets and generates the elastic susceptibility that drives ordering. The reaction-field approximation fixes the local-field correction $\\lambda$ by enforcing the sum rule $(1/N N_q)\\sum_{i,q} \\langle J^z_i(q,0) J^z_i(-q,0)\\rangle = m^2$ through the fluctuation-dissipation theorem, preventing standard RPA from over-softening the mode. Without hyperfine coupling, ordering requires $R_c = (2m^2/\\Delta)\\max[\\lambda_\\mu - \\lambda] \\ge 1$ at $T=0$; Ho$^{3+}$Ga$_5$O$_{12}$ sits at $R_c = 0.71$.","core_discovery":"The paper's central claim is that Ho$^{3+}$Ga$_5$O$_{12}$ realizes hyperfine-assisted magnetic ordering starting from a two-singlet crystal-field ground state. Projecting each Ho$^{3+}$ ion onto its two lowest singlets produces an effective transverse-field Ising Hamiltonian in which the intrinsic transverse field $\\Delta$ opposes the Ising dipolar interaction; these electronic terms alone never order because the maximum soft-mode ratio is $R_c = 0.71 < 1$. The hyperfine term $A J^z I^z$ adds an elastic, zero-frequency contribution to the susceptibility, and this elastic channel diverges first to give long-range order at $T_N = 0.31$ K, in close agreement with the observed 0.3 K. The same model predicts the zero-wavevector magnetic structure found by powder neutron diffraction, and it assigns the low ordering temperature to single-ion quantum fluctuations from the two-singlet crystal field rather than to frustration: with $\\Delta$ set to zero the model orders at 1.74 K, close to the dipolar energy scale.","pith_inferences":["Inference: a decisive numerical test is to include nearest-neighbor exchange at the paper's own upper bound $|K_{nn}| = 0.1D$ and recompute $R_c$; if the ratio crosses 1, hyperfine coupling would no longer be necessary for ordering, which would sharpen or overturn the central claim.","Inference: because the mechanism relies on the $A I^z J^z$ term, substituting ions or isotopes that change the nuclear spin or hyperfine constant should shift $T_N$ predictably; the paper does not make this comparison.","Inference: the authors' analogy with an externally applied transverse field in a doublet system suggests that a transverse magnetic field applied to Ho$^{3+}$Ga$_5$O$_{12}$ should arrest or modify the hyperfine-assisted ordering, an experiment they do not report.","Inference: for other non-Kramers singlet-ground-state magnets, the paper implies that hyperfine coupling should not be neglected in spin-liquid models; it notes Pr$^{3+}$ systems as candidates but does not model them."],"forward_implications":["The magnetic Hamiltonian of Ho$^{3+}$Ga$_5$O$_{12}$ is fixed without fitted exchange: $\\Delta = 7.4$ K, $D = 0.0183$ K, $A = 0.039$ K, and adding nearest-neighbor exchange $|K_{nn}| > 0.1D$ makes the calculated spectra qualitatively worse.","Without nuclear hyperfine coupling the material would remain a correlated paramagnet to zero temperature, since $R_c = 0.71 < 1$; the observed transition requires $A$.","The low $T_N$ is set by single-ion quantum fluctuations from the two-singlet crystal field, not by geometrical frustration; setting $\\Delta = 0$ while keeping $A = 0$ raises the ordering temperature to 1.74 K.","Standard RPA is quantitatively inadequate here: it gives $T_N = 0.66$ K, too high by a factor of two, and overestimates the bandwidth and softening; the reaction-field sum-rule correction is needed.","The lowest mode of the model has the same symmetry as the experimentally observed zero-wavevector magnetic structure, so the long-range dipolar interaction plus easy-axis anisotropy selects that arrangement."],"supporting_citations":[{"why":"Supplies the single-crystal inelastic neutron-scattering data at 4 K, 0.6 K, and about 0.05 K that the three-parameter model must reproduce.","marker":"[25]"},{"why":"Provides the reaction-field formalism for the Ising model in a transverse field that is extended here from one ion per cell to the twelve-ion primitive cell.","marker":"[39]"},{"why":"Introduces the reaction-field correction to mean-field theory that underlies the sum-rule self-consistency calculation.","marker":"[40]"},{"why":"Establishes the two-singlet soft-mode criterion (critical interaction-to-splitting ratio) from which the no-ordering result $R_c = 0.71$ follows.","marker":"[8]"},{"why":"Fixes the two lowest crystal-field singlets and their splitting $\\Delta = 7.4$ K used in the projection to the two-state model.","marker":"[24]"},{"why":"Supplies the Ho nuclear hyperfine coupling constant $A = 0.039$ K used as a fixed parameter.","marker":"[48]"},{"why":"Gives the earlier theoretical prediction that hyperfine coupling can induce magnetic ordering in rare-earth garnets, cited as qualitative support for the mechanism.","marker":"[35]"}],"fun_headline_variants":["Hyperfine coupling orders Ho3Ga5O12 at 0.31 K","Hyperfine, not frustration, drives garnet ordering","Two-singlet garnet orders via hyperfine channel","Ho3Ga5O12's 0.31 K order is hyperfine-driven"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that each Ho$^{3+}$ ion is fully captured by its two lowest crystal-field singlets interacting only through Ising-like dipolar and longitudinal hyperfine couplings, with exchange, nuclear quadrupole, transverse hyperfine, and higher crystal-field levels all negligible; if any of these omitted terms matters, the predicted ordering mechanism and $T_N$ could change.","fun_headline_variants_meta":{"raw":{"variants":["Hyperfine coupling orders Ho3Ga5O12 at 0.31 K","Hyperfine, not frustration, drives garnet ordering","Two-singlet garnet orders via hyperfine channel","Ho3Ga5O12's 0.31 K order is hyperfine-driven"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000193,"raw_usage":{"total_tokens":1346,"prompt_tokens":940,"completion_tokens":406,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":329}},"tokens_in":556,"tokens_out":406,"duration_ms":4153,"temperature":1.0,"reasoning_tokens":329,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:10:16.412448+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the zero-temperature soft-mode ratio $R_c = (2m^2/\\Delta)\\max[\\lambda_\\mu - \\lambda]$ in a model that includes nearest-neighbor exchange up to the paper's own bound $|K_{nn}| = 0.1D$; if $R_c$ reaches or exceeds 1, magnetic ordering would occur without hyperfine coupling and the claim that hyperfine coupling drives the ordering would be refuted.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the single-crystal inelastic neutron-scattering data at 4 K, 0.6 K, and about 0.05 K that the three-parameter model must reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the reaction-field formalism for the Ising model in a transverse field that is extended here from one ion per cell to the twelve-ion primitive cell."},{"cited_title":"Brout, H","cited_arxiv_id":null,"evidence_quote":"Introduces the reaction-field correction to mean-field theory that underlies the sum-rule self-consistency calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the two-singlet soft-mode criterion (critical interaction-to-splitting ratio) from which the no-ordering result $R_c = 0.71$ follows."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Fixes the two lowest crystal-field singlets and their splitting $\\Delta = 7.4$ K used in the projection to the two-state model."},{"cited_title":"Krusius, A","cited_arxiv_id":null,"evidence_quote":"Supplies the Ho nuclear hyperfine coupling constant $A = 0.039$ K used as a fixed parameter."},{"cited_title":"Hammann, P","cited_arxiv_id":null,"evidence_quote":"Gives the earlier theoretical prediction that hyperfine coupling can induce magnetic ordering in rare-earth garnets, cited as qualitative support for the mechanism."}],"review_version":1}