{"id":"1763fb34-f6b1-4c8e-875f-c5f60e71f977","arxiv_id":"2608.11305","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The temperature derivative of a transverse-field-conjugate observable in a pyrochlore magnet is proportional to the ring-exchange specific heat, peaking at the ring-exchange energy scale, with sign set by the zero-flux or pi-flux quantum spin ice state.","lead":"Quantum spin ice is a magnetic state with fractionalized excitations and an energy scale called the ring-exchange scale, which is hard to see in ordinary specific-heat measurements. This paper proposes that measuring the temperature derivative of a conjugate observable, such as anisotropic thermal expansion or dM/dT, should reveal that scale and even show which flux state the material is in.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Numerical verification tests the four-loop boundary scale, not the hexagon ring-exchange scale that enters the material predictions.","rationale":"The paper's central mechanism is a clean application of the chain rule to the effective ring-exchange Hamiltonian. The analytic derivation in SM S2 for the hexagon (B=15) is internally consistent, and the ED on the 16-site cluster confirms the proportionality and the sign law for the boundary four-loop (B=6). I find no algebraic error in Eq. (7) or the sign argument. The most load-bearing gap is that the numerical evidence does not exercise the hexagon ring-exchange sector that enters the material predictions: on the 16-site cluster the low-energy scale is the noncontractible four-loop g4, and the extracted coefficient is 6, not 15. The paper is explicit about this and the sign law is unchanged, but the quantitative claims for alpha_Eg in Pr2Zr2O7 and dM/dT in Ce2Zr2O7 rely on the hexagon mechanism. Since no larger-cluster check or experimental confirmation is provided, the reader's CONDITIONAL verdict is appropriate. The concrete test I propose would settle whether the hexagon sector obeys the same proportionality with the thermodynamic B=15, or whether the four-loop result is a finite-size artifact that changes the material-scale predictions.","tokens_in":24274,"tokens_out":21540,"duration_ms":194182,"concrete_test":"Perform exact diagonalization on a 32-site or 48-site pyrochlore cluster that contains contractible hexagons, for J±/Jzz = ±0.04 and ±0.08, with the Eg transverse field at x = lambda^2/(|J±|Jzz) = 0.014. Extract the lower-peak position T_peak and the peak-shift coefficient kappa from runs at x = 0, 0.007, 0.014, 0.021. Confirm that T_peak scales with |g| = 12|J±|^3/Jzz^2 rather than |g4| = 4J±^2/Jzz, and that kappa |J±| Jzz approaches sgn(J±) * 15 (not +/-6) as |J±| -> 0. If these hold, the material predictions are quantitatively supported; if the peak tracks g4 or kappa deviates from 15, the hexagon-versus-four-loop gap is a real caveat that should be stated in the paper.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central proportionality (7) is verified only on the periodic 16-site cluster, where the low-energy ring sector is the boundary-winding four-loop with g4 = 4J±^2/Jzz and B4(p_Eg) = 6, as stated in the Fig. 1 caption and SM S3.B. The thermodynamic-limit hexagon has g = 12J±^3/Jzz^2 and B(p_Eg) = 15, which is the value entering the material predictions for Pr2Zr2O7 and Ce2Zr2O7. The chain-rule derivation in SM S2 is algebraic and transferable, but the only nonperturbative check that the full XXZ model at finite lambda produces gamma_X proportional to C_ring with the predicted coefficient is demonstrated for the four-loop sector, not for the hexagon. Consequently, the quantitative claims that alpha_Eg (or dM/dT) peaks at the hexagon scale and shifts with kappa = B/(J±Jzz) with B = 15 are extrapolations from a finite-size artifact, even though the sign law is expected to be identical.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a thermodynamic probe of the ring-exchange energy scale in quantum spin ice. The authors consider the XXZ pyrochlore model with a weak transverse perturbation λX̂ and define the response γ_X = d⟨X̂⟩/dT. Using third-order perturbation theory in the transverse exchange, they derive Eq. (7): within the ring-exchange sector, γ_X = −2B(p)λ/(J±Jzz) C_ring(T) + O(λ³), with B(p)≥0, so the sign of γ_X is set by the flux sector and the ring-exchange peak shifts in a flux-dependent direction under the perturbation. They apply this result to anisotropic thermal expansion in non-Kramers Pr pyrochlores and to dM/dT in dipolar-octupolar Ce pyrochlores, and they support the mechanism with exact diagonalization on the 16-site cubic cluster, explicitly noting that the cluster's low-energy ring sector is a boundary four-loop rather than the thermodynamic hexagon.","tokens_in":24359,"tokens_out":14195,"duration_ms":116697,"significance":"If correct, the paper offers an experimentally accessible way to isolate the ring-exchange (photon/vison) energy scale in candidate quantum spin ice materials, addressing a long-standing challenge. The proposal is falsifiable: the proportionality, the sign law, and the peak-shift direction can be checked in measurements, and the θ dependence in the dipolar-octupolar case provides a diagnostic of the pseudospin mixing angle. Strengths include a self-contained perturbative derivation in SM S2, including explicit coefficients B(p_Eg)=15 for the hexagon and B4=6 for the 16-site cluster, and a transparent exact-diagonalization study that confirms the covariance mechanism and the sign structure. The finite-size caveat is openly disclosed. The main limitations are that the numerical check does not directly test the hexagon coefficient and that the sign convention connecting α_Eg to γ_Eg in the non-Kramers application requires an additional stated assumption about the magnetoelastic prefactor.","major_comments":[{"comment":"The numerical verification of Eq. (7) in Fig. 1(b,c) and the extracted coefficient in Fig. 1(d) are performed on the 16-site periodic cluster, where the low-energy ring sector is the boundary-winding four-loop with g4 = 4J±²/Jzz and B4(p_Eg) = 6, not the thermodynamic hexagon with g = 12J±³/Jzz² and B(p_Eg) = 15. Since the material predictions are for the hexagon scale, the main-text statement that 'the proportionality relation in Eq. (7) is verified nonperturbatively by ED' should be qualified to the four-loop variant, and the paper should state explicitly that the quantitative hexagon coefficient rests on the perturbative calculation in SM S2.B rather than on the ED check.","section":"Fig. 1 and SM S3.B / S2.D"},{"comment":"Equation (10) and SM Eq. (S89) give α_Eg = (g_Eg/c_eff^Eg) γ_Eg, so the sign of the measured anisotropic thermal expansion equals the sign of γ_Eg only up to the sign of the magnetoelastic prefactor g_Eg/c_eff^Eg. The abstract's claim that the sign of these signatures distinguishes the zero-flux and π-flux states is therefore not yet a direct statement about α_Eg unless the sign of the magnetoelastic coupling is known or specified; the peak-shift direction in Eq. (8) is robust to that prefactor and should be emphasized as the primary sign diagnostic in the non-Kramers channel.","section":"Non-Kramers pyrochlores, Eq. (10) and SM S5.B"}],"minor_comments":[{"comment":"The sentence 'which realizes Eq. (7) at B=6 on a 16-site cubic cluster' is easily misread as the material value; please add that the thermodynamic-limit value is B=15.","section":"Non-Kramers pyrochlores section"},{"comment":"The text should label the dashed lines as the cluster-specific ±6 and explicitly note that the thermodynamic-limit value would be ±15, so the reader does not confuse the extracted coefficient with the material prediction.","section":"Fig. 1(d) caption"},{"comment":"The statement that α_Eg tracks the lower anomaly of C and 'remains proportional to it through the peak' would benefit from a quantitative statement, for example a plot of α_Eg/C over the relevant temperature range, since Eq. (7) concerns the ring-sector specific heat rather than the total specific heat.","section":"Fig. 1(b,c) discussion"},{"comment":"Because κ = B(p)/(J±Jzz) and B(p) depends on the normalization of p, the normalization convention from SM Eq. (S9) should be recalled near Eq. (8) to avoid ambiguity when comparing coefficients between different probe patterns.","section":"Eq. (8)"},{"comment":"The finite-size caveat that the θ=0.1π response is only a factor 1.4 above this cluster's own g4 should be kept in the main text, as it is currently; consider adding the corresponding thermodynamic-limit factor if a clean estimate is available.","section":"DO section, Fig. 2"}],"recommendation":"minor_revision","confidential_remarks":"The finite-size concern raised in the stress test is real but not disqualifying: the hexagon coefficient B=15 is derived explicitly in SM S2.B, and the ED work transparently verifies the same structural relation on the four-loop cluster. The revision should mainly tighten the language around 'nonperturbative verification' and add the sign-of-prefactor caveat for α_Eg. The paper fits the scope of the journal and, with these clarifications, would be a solid contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the cleanest statement I've seen of a useful trick—if you hit QSI with a weak transverse perturbation, the temperature derivative of the conjugate observable is proportional to C_ring(T) at leading order, with a sign determined by the flux sector. That's new, as far as the cited literature goes, and it gives two concrete experimental channels (thermal expansion in non-Kramers, dM/dT in dipolar-octupolar). The nuclear Schottky argument is sound, and the dM/dT sign separation in Fig. 3 is a nice bonus.\n\nThe algebra in the SM is solid. The chain rule to Eq. (7), the B(p) counts, and the exact identity for the uniform Eg probe within the one-coupling ring Hamiltonian—all check out. The ED on the 16-site cluster confirms the proportionality and the sign law in the perturbative regime. The authors are also honest about the cluster's lower scale being the boundary four-loop g4, not the thermodynamic hexagon.\n\nThe soft spot is exactly that. The only nonperturbative verification of Eq. (7) is for the four-loop sector with B4=6. The material predictions for Pr2Zr2O7 and Ce2Zr2O7 use the hexagon value B=15, and the peak-shift coefficient quoted there is an extrapolation. The sign law is expected to be identical, and the derivation is transferable, but the quantitative claim that alpha_Eg peaks and shifts with the hexagon-scale coefficient is not directly tested. That's a moderate caveat, not a fatal one. The paper says so itself; it just doesn't emphasize how much of the quantitative statement rests on the algebraic extension.\n\nSecond soft spot: the DO prediction requires a purely transverse Zeeman coupling (theta=0 or small). Figure 2 shows that for theta=0.1 pi the response moves to a different scale entirely. Whether Ce2Zr2O7 has a small enough theta is model-dependent. The authors are up front about this, but it limits how sharp the dM/dT diagnostic is for any specific material right now.\n\nOverall: this is a theory paper with a real mechanism, honest numerics, and no parameter fitting. The central proportionality is likely correct under the stated assumptions. The main weakness is that the numerical check and the material prediction are at different loop orders, which should be stated more prominently. I'd send it to peer review, and I'd suggest the referee ask for a 32-site ED check or at least a more prominent caveat in the main text. The paper is worth reading for anyone working on QSI thermodynamics or pyrochlore materials.","headline":"New thermodynamic probe for ring-exchange scale in QSI; clean derivation and honest ED, but the numerical verification sits at the four-loop boundary scale, not the hexagon scale used for materials.","tokens_in":24972,"tokens_out":2445,"would_cite":true,"duration_ms":23935,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A weak transverse perturbation makes the temperature derivative of the conjugate observable proportional to the ring-exchange specific heat in quantum spin ice, providing a thermodynamic probe of the ring-exchange scale and the flux sector.","keywords":["quantum spin ice","pyrochlore lattice","ring exchange","thermal expansion","magnetization derivative","flux sector","emergent photon","covariance spectroscopy"],"falsifier":"A clean experimental test is to measure dM/dT in a dipolar-octupolar pyrochlore with $\\theta$ near zero: if no positive lobe appears at the ring-exchange scale, or if its sign disagrees with the known flux sector, Eq. (7) fails for that material. A numerical falsifier is a cluster calculation whose lowest loop scale is the hexagon $g=12J_\\pm^3/J_{zz}^2$; if the response peak there does not track $g$ with the coefficient $B(p)=15$, the thermodynamic-limit claim is unsupported.","tokens_in":23978,"feed_emoji":"🧊","tokens_out":15761,"duration_ms":123930,"temperature":0.7,"pith_summary":"Quantum spin ice is expected to possess two energy scales: a spinon gap near $J_{zz}/2$ and a lower ring-exchange scale $g=12J_\\pm^3/J_{zz}^2$ that sets the photon bandwidth and vison energy. The paper argues that this lower scale can be seen thermodynamically even when the specific heat cannot resolve it: apply a weak perturbation coupled only to the transverse pseudospin components, and take the temperature derivative of the conjugate observable. The central identity, Eq. (7), states that this derivative is proportional to the ring-exchange specific heat and therefore peaks at the ring-exchange scale. In non-Kramers pyrochlores the probe is realized as the difference of thermal expansion coefficients along $[100]$ and $[010]$; in dipolar-octupolar pyrochlores it is $dM/dT$ under a weak field. The sign of the response and the direction of its peak shift identify whether the ground state has zero flux or $\\pi$-flux.","feed_headline":"Spin ice's hidden ring-exchange scale shows in one derivative","feed_subtitle":"Strain or a weak field exposes the low-energy scale, and its sign tells the flux sector.","key_machinery":"The load-bearing object is the dressed ring-exchange Hamiltonian $H_{\\mathrm{ring}}(\\lambda)=-\\sum_\\mu g^{\\mathrm{eff}}_\\mu(\\lambda)W_\\mu$, where $W_\\mu$ flips the six spins around a hexagon of orientation $\\mu$. A weak transverse field $\\lambda\\hat X$ acts within the ice manifold only in pairs, replacing one exchange insertion with two field insertions and renormalizing each hexagon coupling by the factor $1+B_\\mu(p)\\lambda^2/(J_\\pm J_{zz})$. Because the field enters the low-energy sector only through these couplings, the chain rule converts the covariance $T^{-2}\\mathrm{Cov}_T(\\hat X,H_{\\mathrm{ring}})$ into $-2B(p)\\lambda/(J_\\pm J_{zz})\\, C_{\\mathrm{ring}}(T)$, which is Eq. (7).","core_discovery":"The paper's claim is the proportionality $\\gamma_X|_{\\mathrm{ring}} = -2B(p)\\lambda/(J_\\pm J_{zz})\\, C_{\\mathrm{ring}}(T) + O(\\lambda^3)$, where $\\gamma_X=d\\langle \\hat X\\rangle/dT$ is the thermal response of an observable conjugate to a weak transverse perturbation of strength $\\lambda$, $B(p)\\ge 0$ is a form-factor coefficient, and $C_{\\mathrm{ring}}$ is the specific heat of the effective ring-exchange Hamiltonian. Within the low-temperature ice manifold (states with zero spinon charge on every tetrahedron), the perturbation only dresses the ring-exchange coupling, so the response tracks the ring-exchange anomaly rather than the spinon anomaly; the sign of the proportionality constant is set by $J_\\pm$, positive in the $\\pi$-flux sector and negative in the zero-flux sector. For the sublattice-uniform $E_g$ probe the relation is exact within a one-coupling ring Hamiltonian. Exact diagonalization on a 16-site pyrochlore cluster verifies the proportionality and the sign law in the perturbative limit, albeit at the boundary four-loop scale $g_4=4J_\\pm^2/J_{zz}$ rather than the thermodynamic hexagon scale $g$.","pith_inferences":["If the proportionality survives in the thermodynamic limit, the same covariance spectroscopy should transfer to other frustrated magnets or lattice gauge theories whose low-energy sector is a single dominant interaction renormalized by a weak conjugate field.","A decisive next test would be exact diagonalization or quantum Monte Carlo on a cluster whose lowest loop scale is the hexagon $g=12J_\\pm^3/J_{zz}^2$ rather than the boundary four-loop $g_4$; I would expect the coefficient $B(p)$ to move from 6 toward 15 as the thermodynamic limit is approached.","The sign diagnostic suggests a practical phase-diagram tool: applying uniaxial stress or a field while tracking the sign of the thermal response could map zero-flux and pi-flux regions in material parameter space, an extension the paper does not develop.","Because the method measures a covariance between an observable and the energy rather than an energy alone, it may remain useful in disordered or partially ordered candidates where specific-heat peaks are broadened but a correlated response stays sharper."],"forward_implications":["In non-Kramers pyrochlores such as Pr2Zr2O7 and Pr2Hf2O7, the difference of thermal expansion coefficients along [100] and [010] should show a peak at the ring-exchange scale, with no direct nuclear Schottky background.","In dipolar-octupolar pyrochlores with a purely transverse Zeeman coupling (theta = 0), |dM/dT| peaks at the ring-exchange scale.","The sign of the response (positive for pi-flux, negative for zero-flux), combined with the direction of the peak shift under increasing strain or field, distinguishes the two flux sectors.","dM/dT can separate ring-exchange and spinon features that merge in the specific heat, because the two contributions enter with opposite signs.","A nonzero mixing angle theta moves the low-temperature dM/dT extremum away from the ring-exchange scale to a field-set scale, making the measurement a diagnostic of theta."],"supporting_citations":[{"why":"Derives the ring-exchange effective Hamiltonian and the emergent photon description from the XXZ pyrochlore model, fixing the energy scale g.","marker":"[1]"},{"why":"Shows that the sign of the ring coupling selects the zero-flux or pi-flux ground state, the basis of the flux-sector diagnostic.","marker":"[7, 13–16]"},{"why":"Quantum Monte Carlo results establishing the two-peak specific-heat structure, spinon and ring-exchange, that the new probe is designed to resolve.","marker":"[17–21]"},{"why":"Provides the perturbation theory in which two field insertions replace an exchange insertion, giving the dressed ring coupling.","marker":"[47, 48]"},{"why":"Supplemental material containing the full derivation of Eq. (7), the B coefficients, the thermal-expansion relation, and the absence of a direct hyperfine term.","marker":"[22]"},{"why":"Provides the magnetoelastic coupling framework through which Eg strain realizes the transverse probe and maps the response to thermal expansion.","marker":"[49–54]"},{"why":"Models the hyperfine-enhanced nuclear Schottky contribution whose obscuring of the heat-capacity peak the strain probe avoids.","marker":"[23–26]"},{"why":"Defines the mixing angle theta of the dipolar-octupolar pseudospin rotation that controls whether a magnetic field acts purely transverse.","marker":"[30, 36, 39, 56, 58]"}],"fun_headline_variants":["A single derivative exposes spin ice's ring-exchange scale","Sign of a thermal derivative tells spin ice's flux sector","Weak strain or field reads spin ice's ring-exchange peak"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim collapses if the low-energy sector of the standard nearest-neighbor XXZ pyrochlore model is not a single ring-exchange term that a weak transverse strain or field merely dresses; for dipolar-octupolar materials it also requires a purely transverse Zeeman coupling, and for the thermodynamic limit the 16-site four-loop verification must be replaced by a genuine hexagon check.","fun_headline_variants_meta":{"raw":{"variants":["A single derivative exposes spin ice's ring-exchange scale","Sign of a thermal derivative tells spin ice's flux sector","Weak strain or field reads spin ice's ring-exchange peak"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000966,"raw_usage":{"total_tokens":4189,"prompt_tokens":1103,"completion_tokens":3086,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":719,"completion_tokens_details":{"reasoning_tokens":3034}},"tokens_in":719,"tokens_out":3086,"duration_ms":21474,"temperature":1.0,"reasoning_tokens":3034,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:13:23.417894+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A clean experimental test is to measure dM/dT in a dipolar-octupolar pyrochlore with $\\theta$ near zero: if no positive lobe appears at the ring-exchange scale, or if its sign disagrees with the known flux sector, Eq. (7) fails for that material. A numerical falsifier is a cluster calculation whose lowest loop scale is the hexagon $g=12J_\\pm^3/J_{zz}^2$; if the response peak there does not track $g$ with the coefficient $B(p)=15$, the thermodynamic-limit claim is unsupported.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplemental material containing the full derivation of Eq. (7), the B coefficients, the thermal-expansion relation, and the absence of a direct hyperfine term."}],"review_version":1}