{"id":"719fd59e-a5ae-4df0-8e08-6c05274685f4","arxiv_id":"2608.03850","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Ho2TiO5 retains local <111> Ising anisotropy and ice-rule-compatible spin configurations on a magnetically incomplete tetrahedral network, making it a vacancy-disrupted spin ice.","lead":"Neutron scattering experiments on the stuffed pyrochlore Ho2TiO5 show that holmium moments keep their spin-ice-like local Ising character even though many magnetic sites are missing. The material behaves as a vacancy-disrupted spin ice: local ice-rule tendency survives on an incomplete tetrahedral network.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"RMC config fractions in Ho2TiO5 are close to random-dilution expectation; no null-model comparison is given, so the 'local ice-rule tendency' claim is not established.","rationale":"The reader's baseline concern is valid for the Ho2Ti2O7 reference but is less load-bearing for the Ho2TiO5 central claim, because the Ho2TiO5 40 K mPDF moment is near zero. The more fundamental issue is that the key quantitative evidence for 'local ice-rule tendency' in Ho2TiO5 is the tetrahedral configuration distribution from RMC, and that distribution has not been compared to a random-dilution null model. The paper's own structural model already forces 80% of tetrahedra to be incomplete, and 2-in/1-out/1-in/2-out configurations are statistically dominant among incomplete tetrahedra even with zero spin-spin correlations. The complete-tetrahedron fractions, the only place where ice-rule preference could be cleanly identified, are nearly random. Without a null-model control, the headline characterization 'vacancy-disrupted spin ice' overstates what the data constrain. The paper has independent support from pNPD for <111> single-ion anisotropy and is transparent about the Pnma overfitting problem, but those do not establish short-range ice-rule correlations. The verdict should remain conditional, with the explicit requirement of a null-model comparison; if the null model reproduces the data, the central claim would need to be rejected.","tokens_in":13006,"tokens_out":9404,"duration_ms":100533,"concrete_test":"Generate at least 100 random spin configurations on the same 5x5x5 supercell with 2/3 Ho occupancy on A and B sites (random vacancy placement) and each Ho spin independently chosen to point in or out along its local <111> axis. Compute (i) the tetrahedral configuration fractions and (ii) the magnetic pair distribution function using the same Qmax and Fourier filter as the paper. Compare to the observed Ho2TiO5 0.3 K fractions and mPDF. If the null ensemble reproduces the observed fractions and mPDF within statistical uncertainty, the data provide no evidence for a local ice-rule tendency; if the 2-in/2-out fraction and nearest-neighbor mPDF peak differ significantly from the null model, the claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that Ho2TiO5 retains a local ice-rule tendency rests on tetrahedral configuration fractions from RMC: incomplete tetrahedra are mostly 2-in/1-out and 1-in/2-out, and complete tetrahedra are 10% 2-in/2-out, 9% 3-in/1-out, 1% all-in/all-out. These numbers are close to a trivial null model. For random 2/3 Ho occupancy and independent random in/out spins, incomplete tetrahedra are 2-in/1-out or 1-in/2-out with probability 3/8+3/8 = 75%, and complete tetrahedra are 2-in/2-out with probability 6/16 = 37.5%, 3-in/1-out with 50%, and all-in/all-out with 12.5%. The observed complete-tetrahedron fractions (~50%, 45%, 5%) are within statistical reach of random, yet no error bars or null-model comparison are reported for Ho2TiO5. The RMC was run with spins constrained along <111>; the unconstrained fit was dismissed as overfitting, which implies the diffuse scattering may be too weak to determine spin correlations independently. Thus the data as presented do not demonstrate a residual ice-rule tendency beyond random dilution plus single-ion anisotropy.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a multi-technique neutron scattering study of the stuffed pyrochlore Ho2TiO5 and the parent spin-ice compound Ho2Ti2O7. Using atomic and magnetic pair distribution function (PDF/mPDF) analysis, reciprocal-space reverse Monte Carlo (RMC) modeling of diffuse magnetic scattering, and half-polarized neutron powder diffraction (pNPD), the authors determine the average and local atomic structure and the local magnetic correlations. They find that Ho2TiO5 retains the average Fd-3m pyrochlore structure with Ho/Ti mixing on both A and B sites at approximately 2/3 occupancy, and that the Ho moments on both sublattices retain uniaxial <111> Ising anisotropy, with a reduced anisotropy ratio relative to Ho2Ti2O7. For the parent compound, RMC and mPDF yield a nearly ideal 2-in/2-out spin-ice state at 0.3 K (about 95% ice-rule compliance). For Ho2TiO5, the RMC analysis indicates that most tetrahedra are magnetically incomplete, with dominant 2-in/1-out and 1-in/2-out configurations on incomplete tetrahedra and a small excess of 2-in/2-out over 3-in/1-out on complete tetrahedra. The authors conclude that Ho2TiO5 is a vacancy-disrupted spin ice in which a local ice-rule tendency survives on an incomplete tetrahedral network.","tokens_in":13423,"tokens_out":5630,"duration_ms":55440,"significance":"If fully established, the conclusion that a local ice-rule tendency survives strong magnetic dilution in a stuffed pyrochlore would be of genuine interest to the frustrated-magnetism community, as it bears on the robustness of local spin-ice correlations under disorder. The paper has notable strengths: the pNPD analysis is performed with an unconstrained susceptibility tensor and independently recovers uniaxial <111> anisotropy for both A- and B-site Ho, providing direct microscopic evidence for the persistence of single-ion Ising character. The combination of real-space mPDF and reciprocal-space RMC is appropriate and, for the reference compound Ho2Ti2O7, yields a clean benchmark (95% 2-in/2-out at 0.3 K). The authors also explicitly flag the Pnma overfitting issue and choose the more conservative Fd-3m model. However, the central claim regarding Ho2TiO5 is undermined by the lack of a statistical null-model comparison for the RMC-derived tetrahedral configuration fractions; the observed fractions appear close to those expected for random in/out spins on a randomly diluted lattice, so the paper's central assertion is not yet quantitatively supported.","major_comments":[{"comment":"The central claim of a residual 'local ice-rule tendency' in Ho2TiO5 rests on the RMC-derived tetrahedral configuration fractions. No null-model comparison or error bars are reported for these fractions. For random 2/3 Ho occupancy and independent random in/out spins, complete tetrahedra have P(2-in/2-out)=6/16=37.5%, P(3-in/1-out)=8/16=50%, and P(4-in/0-out)=2/16=12.5%; incomplete 3-spin tetrahedra have P(2-in/1-out or 1-in/2-out)=6/8=75%. The observed complete-tetrahedron fractions, approximately 50% 2-in/2-out, 45% 3-in/1-out, and 5% all-in/all-out (derived from the text's 10%, 9%, and 1% of all tetrahedra), are only modestly different from this null distribution, and the dominance of 2-in/1-out or 1-in/2-out among incomplete tetrahedra is exactly the null expectation. Furthermore, the statement that these configurations are 'locally compatible with the ice rule if the missing Ho spin","section":"§IV.2, Fig. 9(c-i)"},{"comment":"The 50 K dataset is used as a magnetic-free reference for both the mPDF and the temperature-subtracted RMC analyses. The manuscript states that magnetic correlations at 50 K are 'negligible' but does not demonstrate this. For the reference compound Ho2Ti2O7, Fig. 6(c) shows a finite effective correlated moment of 1.47 μB at 40 K, indicating that residual magnetic correlations can persist well above the 'negligible' range in related materials. If a similar residual exists at 50 K in Ho2TiO5 or Ho2Ti2O7, the subtraction removes some real magnetic signal and biases the extracted low-temperature amplitudes and hence the configuration fractions. Since this baseline assumption underlies the quantitative separation of magnetic from nuclear scattering throughout Sections III.C and IV, the authors should either justify the 50 K reference with data (e.g., show that the 50 K mPDF/RMC signal is feat","section":"§III.C and §IV.2"},{"comment":"The mPDF refinement of Ho2TiO5 uses a pre-assumed 2-in/2-out spin configuration extended to include both Ho sites. Because the model explicitly imposes spin-in/spin-out ordering on the tetrahedra, the good fit to the negative nearest-neighbor peak at ~3.6 Å does not independently confirm a spin-ice-like correlation; it only shows that such a model is consistent with the data. The mPDF analysis is therefore best viewed as a forward-model estimate of the correlated moment amplitude, not as evidence for the ice-rule tendency. The manuscript should acknowledge this circularity more explicitly, or ideally use a more agnostic starting model (e.g., a configuration with random in/out spins) to demonstrate that the data require the ice-rule-like correlations.","section":"§III.C.2, Fig. 6(b)"}],"minor_comments":[{"comment":"Typo: 'neligble' should be 'negligible'.","section":"Fig. 6 caption"},{"comment":"The text is ambiguous about whether the quoted fractions (10% 2-in/2-out, 9% 3-in/1-out, 1% all-in/all-out) are relative to all tetrahedra or to the 20% of complete tetrahedra. Please clarify and, if possible, report fractions conditional on the number of occupied vertices. Also, error bars are shown in Fig. 8 for Ho2Ti2O7 but not in Fig. 9(c-i); adding them would directly support the statistical comparison requested in the major comments.","section":"§IV.2, Fig. 9"},{"comment":"The text says the mPDF model for Ho2TiO5 'used Ho2Ti2O7 refinement as the starting point' and extends to both sites. It would be helpful to describe how the B-site (16d) Ho spin directions are defined on their local <111> axes, since the B-site tetrahedral network is different from the A-site network. A brief description of the spin construction would improve reproducibility.","section":"§III.C.2"},{"comment":"The statement that the unconstrained RMC for Ho2TiO5 produced 'randomly oriented spins with no physically meaningful correlations, indicative of overfitting' is important, but it is not quantified. For a fair assessment, please report the fit quality (e.g., R-factor or chi-square) for both the constrained and unconstrained models so readers can judge whether the diffuse scattering carries enough information to determine spin orientations independently of the <111> constraint.","section":"§IV.1"}],"recommendation":"major_revision","confidential_remarks":"The null-model concern is the central issue: the paper's headline claim is not yet supported by the reported RMC fractions, which appear close to random-dilution expectations. The pNPD result is strong and likely publishable on its own, but the 'vacancy-disrupted spin ice' narrative needs either a quantitative null-model comparison or a softened conclusion. I would encourage the editors to send the manuscript back with a request for that analysis; it is within scope and should be feasible given the existing RMC ensembles."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Name],\n\nThe one thing to know: this paper contains one genuinely solid, new result — the half-polarized neutron diffraction shows uniaxial <111> susceptibility tensors on both the A- and B-site Ho in Ho2TiO5, and that's an independent determination, not an input constraint. The other main claim, that a local ice-rule tendency survives in the vacancy-disrupted network, is plausible but not yet demonstrated. The RMC configuration fractions are presented without a null-model comparison or error bars, and for incomplete tetrahedra the \"dominant\" 2-in/1-out / 1-in/2-out pattern is exactly what random dilution plus random in/out spins already gives.\n\nWhat the paper does well: the pNPD analysis is careful. The unconstrained refinement yields nearly identical anisotropic tensors for both sites, so the single-ion Ising character survives stuffing — a real and useful result. The structural story is also convincing: 2/3 Ho on both sites in Fd-3m, with low-r PDF misfits attributed to local bond disorder rather than symmetry lowering. The authors are transparent about the Pnma overfitting problem, which is good practice. The mPDF benchmark on Ho2Ti2O7 is a nice first — no one had directly fit mPDF to data for the canonical spin ice before.\n\nSoft spots, in proportion. The biggest is the missing null model for the RMC fractions. For the complete tetrahedra, the observed 50% 2-in/2-out, 45% 3-in/1-out, 5% all-in/all-out actually leans toward the ice rule relative to the random 37.5/50/12.5 — so the stress-test note that says \"within statistical reach\" is not obviously right, and the direction of the deviation supports the paper. But without error bars (the 100-refinement averaging only gives a mean, no spread) and without the authors telling us how many complete tetrahedra are in the supercell, we cannot assess significance. For the incomplete tetrahedra, the argument is genuinely weak: random spins already make 2-in/1-out the dominant configuration. If the authors want to claim a residual ice-rule tendency, they need to compare their full configuration histogram to a random-occupancy/random-spin null, and show that the ice-rule-compatible fraction is enriched.\n\nSecond, the 50 K subtraction baseline. Figure 6(c) shows 1.47 uB at 40 K for Ho2Ti2O7, so \"magnetic correlations are negligible\" at 50 K is an assumption, not a measurement. If there is residual signal at 50 K, the mPDF and temperature-subtracted RMC lose some real magnetic scattering. This is a minor-to-moderate issue and fixable with a quick check.\n\nThird, the mPDF fit to Ho2TiO5 uses a pre-assumed 2-in/2-out model, so the good fit is partly circular for the central question.\n\nWho this is for: frustrated-magnetism experimentalists working on stuffed pyrochlores will want the anisotropy result and the structural characterization. The vacancy-disrupted spin-ice question is worth a serious referee, but the manuscript as written doesn't close it. I'd send it to review with explicit requests for a null-model comparison, error bars on the config fractions, a 50 K baseline check, and deposition of the RMC spin configurations and pNPD refinement parameters — \"data available upon request\" is not enough for this kind of analysis.\n\nRecommendation: send to peer review.","headline":"Solid site-resolved <111> anisotropy result from half-polarized diffraction, but the 'vacancy-disrupted spin ice' conclusion needs a null-model test and error bars before it lands.","tokens_in":13847,"tokens_out":9195,"would_cite":true,"duration_ms":95457,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The stuffed pyrochlore Ho2TiO5 preserves local ice-rule-compatible order even though most tetrahedra are missing a magnetic spin.","keywords":["spin ice","stuffed pyrochlore","Ho2TiO5","magnetic pair distribution function","reverse Monte Carlo","Ising anisotropy","magnetic dilution","frustrated magnetism"],"falsifier":"Measure the same Ho2TiO5 powder with 100 K (or higher) as the magnetic-free reference instead of 50 K and repeat the mPDF and RMC analyses; if the extracted low-temperature correlation amplitudes and tetrahedral configuration fractions change, the baseline assumption is wrong.","tokens_in":12963,"feed_emoji":"🧊","tokens_out":12470,"duration_ms":122185,"temperature":0.7,"pith_summary":"The paper asks what happens to spin-ice physics when the magnetic lattice is punched with vacancies. Starting from the canonical spin ice Ho2Ti2O7, the authors study Ho2TiO5, a stuffed pyrochlore where holmium and titanium share both cation sites so that roughly 80% of tetrahedra lack one or more magnetic spins. Combining real-space magnetic pair distribution function analysis, reverse Monte Carlo modeling of diffuse neutron scattering, and half-polarized neutron powder diffraction, they find that the holmium moments keep their local <111> Ising easy axes on both sublattices. Most incomplete tetrahedra adopt 2-in/1-out or 1-in/2-out configurations, which become ice-rule-valid two-in/two-out tetrahedra if the missing spin is restored. The conclusion is that Ho2TiO5 is a vacancy-disrupted spin ice: the single-ion prerequisite survives, the local ice-rule tendency persists, but collective two-in/two-out ordering is destroyed by quenched disorder rather than by thermal fluctuations.","feed_headline":"Ice-rule tendency survives when most tetrahedra lack a spin","feed_subtitle":"Neutron data show Ho2TiO5 keeps Ising spins and ice-rule-compatible local patterns despite dilution.","key_machinery":"The object that carries the argument is the tetrahedron of the pyrochlore network, classified by how many of its four vertices carry a Ho moment pointing into or out of the center. The paper's central named result is 'vacancy-disrupted spin ice': the local ice-rule constraint survives on a topologically incomplete tetrahedral network. The mechanism is tested with three coordinated probes—real-space magnetic pair distribution function for nearest-neighbor spin-in/spin-out correlations, reciprocal-space reverse Monte Carlo for counting tetrahedral configurations in a supercell with random 2/3 Ho occupancy, and half-polarized neutron powder diffraction for the site-resolved susceptibility tenso","core_discovery":"Both compounds keep the average pyrochlore Fd-3m structure; Ho2TiO5 is best modeled with two-thirds Ho on both cation sites and local bond disorder. In Ho2Ti2O7, magnetic pair distribution function and reverse Monte Carlo analysis find ~95% two-in/two-out tetrahedra at 0.3 K, with progressive thermal disordering on warming. In Ho2TiO5, ~80% of tetrahedra are magnetically incomplete; the dominant configurations are 2-in/1-out and 1-in/2-out, which become ice-rule-valid if the missing Ho spin is restored. Half-polarized neutron powder diffraction shows <111> Ising anisotropy is retained on both sublattices, with a reduced anisotropy ratio. The paper's claim: Ho2TiO5 is a vacancy-disrupted spin","pith_inferences":["A testable extension the paper does not perform: vary the stuffing level x in Ho2(Ti1−xHox)2O7−δ and map tetrahedral configuration fractions versus x; the same analysis should show a crossover from near-ideal ice-rule compliance to vacancy-disrupted behavior, with 2-in/1-out dominating once dilution is strong.","If the local ice-rule tendency is genuine, heat-capacity or magnetocaloric measurements on Ho2TiO5 may reveal a reduced residual entropy plateau, although vacancy disorder should smear any sharp Pauling-type value.","The similar <111> anisotropy found for A-site and B-site Ho suggests other pyrochlores with magnetic ions on both sublattices may retain Ising character under stuffing; half-polarized neutron diffraction offers a direct test without single crystals."],"forward_implications":["In Ho2TiO5, the ice-rule-compatible fraction is set by quenched Ho/Ti disorder rather than by thermal population of defects: warming changes configuration fractions only weakly, so the magnetic response is disorder-dominated rather than monopole-dominated.","Local spin-in/spin-out correlations extend only about 6–7 Å in Ho2TiO5, so collective ice-rule order is genuinely short-range despite the retained single-ion anisotropy.","The ~20% magnetically complete tetrahedra match the (2/3)^4 expectation for random site occupancy, validating the structural model against correlated-occupancy alternatives.","The reduced anisotropy ratio (χ∥/χ⊥ ≈ 2.1) implies that the local crystal-field environments in the stuffed lattice are broadened, so future crystal-field calculations must include site disorder to reproduce the susceptibility."],"supporting_citations":[{"why":"Establishes Ho2Ti2O7 as a geometrically frustrated pyrochlore with local <111> Ising moments, the parent spin-ice benchmark.","marker":"[2]"},{"why":"Documents dipolar spin-ice correlations in Ho2Ti2O7, the reference low-temperature state the paper benchmarks against.","marker":"[4]"},{"why":"Characterizes the stuffed pyrochlore Ho2TiO5 structural disorder and magnetic properties, supplying the compound and the linear M/H regime for pNPD.","marker":"[7]"},{"why":"Shows long- and short-range order in stuffed titanate pyrochlores, supporting the Fd-3m average structure with mixed A/B occupancy.","marker":"[9]"},{"why":"Provides the magnetic pair distribution function formalism used to extract real-space spin-ice-like correlations.","marker":"[10, 11]"},{"why":"Supplies the reverse Monte Carlo program used to count tetrahedral in/out configuration fractions from diffuse magnetic scattering.","marker":"[12]"},{"why":"Gives the polarized neutron powder diffraction method for extracting the local susceptibility tensor, establishing the <111> easy axis.","marker":"[18]"},{"why":"Documents the mPDF and half-polarized pNPD setup at HB-2A and provides the Ho2Ti2O7 susceptibility comparison data.","marker":"[13]"},{"why":"Supplies the real-space temperature-subtraction procedure used to isolate the magnetic PDF signal from atomic contributions.","marker":"[19]"}],"fun_headline_variants":["Vacancy-disrupted spin ice still favors ice rule","Missing Ho spins don't erase ice-rule tendency","Ice-rule logic survives in partial spin ice","Ho2TiO5: diluted spin ice keeps local Ising order"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The analysis treats the 50 K data as magnetically empty and subtracts them from all lower-temperature data; if short-range magnetic correlations still exist at 50 K, the extracted magnetic signal is biased.","fun_headline_variants_meta":{"raw":{"variants":["Vacancy-disrupted spin ice still favors ice rule","Missing Ho spins don't erase ice-rule tendency","Ice-rule logic survives in partial spin ice","Ho2TiO5: diluted spin ice keeps local Ising order"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000197,"raw_usage":{"total_tokens":1261,"prompt_tokens":865,"completion_tokens":396,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":609,"completion_tokens_details":{"reasoning_tokens":331}},"tokens_in":609,"tokens_out":396,"duration_ms":4825,"temperature":1.0,"reasoning_tokens":331,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T11:05:50.161042+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the same Ho2TiO5 powder with 100 K (or higher) as the magnetic-free reference instead of 50 K and repeat the mPDF and RMC analyses; if the extracted low-temperature correlation amplitudes and tetrahedral configuration fractions change, the baseline assumption is wrong.","supporting_citations":[{"cited_title":"The isolation of the magnetic PDF signal was done in a similar manner to Ho 2Ti2O7 usingF d¯3m space group with Ho occupancy 2/3 on both A and B sites","cited_arxiv_id":null,"evidence_quote":"Establishes Ho2Ti2O7 as a geometrically frustrated pyrochlore with local <111> Ising moments, the parent spin-ice benchmark."},{"cited_title":"To isolate the diffuse magnetic signal a temperature subtraction ap- proach was employed","cited_arxiv_id":null,"evidence_quote":"Documents dipolar spin-ice correlations in Ho2Ti2O7, the reference low-temperature state the paper benchmarks against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Characterizes the stuffed pyrochlore Ho2TiO5 structural disorder and magnetic properties, supplying the compound and the linear M/H regime for pNPD."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the reverse Monte Carlo program used to count tetrahedral in/out configuration fractions from diffuse magnetic scattering."},{"cited_title":"Juh´ as, C","cited_arxiv_id":null,"evidence_quote":"Gives the polarized neutron powder diffraction method for extracting the local susceptibility tensor, establishing the <111> easy axis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the mPDF and half-polarized pNPD setup at HB-2A and provides the Ho2Ti2O7 susceptibility comparison data."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the real-space temperature-subtraction procedure used to isolate the magnetic PDF signal from atomic contributions."}],"review_version":1}