{"id":"3e3704d3-25d5-4a5b-82ca-5219884d015c","arxiv_id":"2607.14029","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In Gd2Ru2O7, terahertz spectroscopy directly measures the Ru-induced exchange field on Gd (≈5–6 T, θ≈70–79°) and, via phonon anomalies, reveals a hidden Gd-ordering transition near 10 K.","lead":"Terahertz spectra of Gd2Ru2O7 resolve an exchange-split Gd3+ magnetic mode and an optical phonon in the same measurement, yielding an internal exchange field of about 5–6 T with orientation near 78° and evidence for a Gd-ordering transition that bulk probes missed.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Powder sample: extracted θ is not a well-defined single-crystal angle; the claimed agreement with θ_Gd≈80° is unsupported without powder-averaging analysis or single-crystal measurement.","rationale":"The central claim hinges on the quantitative agreement between the extracted θ and the theoretical value. The data were taken on a powder, but the model assumes a unique angle. Because the distribution width is an order of magnitude larger than the observed linewidth, the single-θ interpretation is geometrically impossible for a random powder. This is a stronger and more concrete flaw than the reader's mode-identity concern, because it does not require speculation about alternative excitations; it follows directly from the sample geometry. The paper does not mention powder averaging, texturing, or grain orientation, and the reader's verdict did not flag this. I therefore recommend keeping the conditional verdict but adding this specific requirement: the authors must either provide single-crystal data or a proper powder-averaged analysis. If neither is provided, the headline agreement should be removed. Thus verdict_should_be = CONDITIONAL, with the additional condition.","tokens_in":10894,"tokens_out":11693,"duration_ms":114556,"concrete_test":"Perform the same THz-TDS measurement on a single crystal of Gd2Ru2O7 with known crystallographic orientation relative to B_ext, and repeat the Lorentzian fits and the Δ² analysis. If the single-crystal spectra show a comparable narrow line and yield θ≈70–80°, the powder concern is resolved. Alternatively, compute the powder-averaged absorption spectrum using the reported Δ_ext, Δ_int and the theoretical θ_Gd≈80° with a random distribution of grain orientations, and compare it with Fig. 3(a); a predicted broad continuum versus an observed narrow peak would demonstrate that the single-θ extraction is invalid.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The most load-bearing concern is that the analysis treats the sample as a single crystal. The measurements are on a pressed powder pellet (Supplemental Material, Sample Synthesis). The field-dependent splitting is fit with Δ² = Δ_ext² + Δ_int² + 2Δ_extΔ_int cosθ, where θ is a single angle between B_ext and B_int. In a random powder, θ varies among crystallites; the absorption would be an inhomogeneous sum spanning |Δ_ext−Δ_int| to Δ_ext+Δ_int, with spectral weight proportional to Δ. With Δ_ext≈0.202 THz and Δ_int≈0.155 THz, this spans ≈0.05–0.36 THz, far broader than the fitted Lorentzian (τ_m≈10–20 ps, FWHM≈0.02 THz). A single narrow resonance is incompatible with random powder averaging, so the extracted θ≈70–79° is not a well-defined crystallographic angle and cannot be compared with the predicted θ_Gd≈80°. The zero-field Δ_int may still give B_int≈5.4 T, but the orientation and the field-fit B_int=6.3(5) T are not meaningful without a powder-averaging model or a single-crystal measurement. This is an internal inconsistency between the stated sample form and the analysis, independent of the mode-assignment question.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports magneto-THz time-domain spectroscopy of a pressed powder pellet of Gd2Ru2O7. It identifies a low-frequency mode that appears below ~30 K, assigns it to a transition within the exchange-split J=7/2 Gd3+ multiplet, and uses its frequency under applied fields up to 7 T to extract an internal exchange field B_int ≈ 5–6 T and an angle θ ≈ 70–79° between the internal and external fields, which the authors compare with the first-principles prediction θ_Gd ≈ 80° for cluster-multipolar order. It additionally tracks a ~1.4 THz optical phonon whose frequency, linewidth, and oscillator strength show anomalies near the Ru Néel temperature (≈110 K) and near 10 K (0 T) / 30 K (7 T); the latter anomaly is interpreted as evidence for a genuine Gd-sublattice ordering transition that bulk thermodynamic probes miss. The paper claims that a single THz spectrum can simultaneously quantify the internal field at the rare-earth site and reveal hidden rare-earth order through spin-phonon coupling.","tokens_in":11308,"tokens_out":5913,"duration_ms":65208,"significance":"If the central claims held, the work would be significant: it would provide a spectroscopic route to internal exchange fields and hidden multipolar order in rare-earth pyrochlores, and it would connect magnetic excitations and phonons in a single measurement. The authors include useful internal consistency checks, notably the integrated spectral weight analysis in Fig. S3, which shows that the main spectral trends are not purely artifacts of the Lorentzian fitting. The g-factor extracted from the field-dependent fit (g = 1.86 ± 0.17) is also close to the expected Gd3+ value. However, the central quantitative claims — the orientation θ and the identification of a genuine Gd ordering transition — rest on assumptions that are not met for the measured sample and on an unproven mode assignment. The paper is therefore not yet a reliable basis for the advertised unified framework.","major_comments":[{"comment":"The sample is a pressed powder pellet, but the central analysis treats the measured resonance as arising from a single relative angle θ between B_ext and B_int. In a random powder, θ is distributed over crystallites with uniform cosθ, so the resonance frequency at fixed B_ext should be a broad inhomogeneous distribution spanning from |Δ_ext − Δ_int| to Δ_ext + Δ_int. With Δ_ext ≈ 0.202 THz and Δ_int ≈ 0.155 THz at 1.6 K and 7 T, this spans ≈ 0.05–0.36 THz, more than an order of magnitude wider than the fitted Lorentzian (τ_m ≈ 10–20 ps corresponds to FWHM ≈ 0.01–0.02 THz). The observation of a single narrow resonance is incompatible with random powder averaging. Consequently, the extracted θ ≈ 70–79° is not a well-defined crystallographic angle and cannot be compared with the predicted θ_Gd ≈ 80°. The zero-field B_int ≈ 5.4 T is less affected, but the field-dependent B_int = 6.3(5) T and","section":"Supplemental Material, Sample Synthesis; Figures 2(c) and 3(c); law-of-cosines formula"},{"comment":"The assignment of the low-frequency mode to a transition between exchange-split Gd3+ levels is inferred from its low-temperature emergence and field-induced hardening, but no definitive evidence is given. The mode could in principle be an impurity excitation, a two-magnon process, a phonon, or a Ru-related spin-gap excitation; the authors do not provide polarization analysis, high-field saturation data, or a quantitative CEF/exchange model that predicts the mode's oscillator strength, temperature dependence, or selection rules. Since all extracted values of B_int and θ depend on this assignment, the mode identification is load-bearing and needs direct support.","section":"Results and Discussion, identification of the low-frequency mode; Figs. 2(a)–(d) and 3(a)–(d)"},{"comment":"The claim that the phonon anomalies near 10 K (0 T) and 30 K (7 T) identify a 'genuine ordering transition' rather than a gradual crossover is not supported quantitatively. The anomalies are identified by eye in Lorentzian fit parameters, with no error bars, no statistical comparison to a smooth crossover model, and no independent confirmation (e.g., specific heat, neutron diffraction, or Mössbauer). Spin-phonon coupling can produce smooth or broad anomalies across a crossover, and the absence of a corresponding feature in the magnetic mode is explained by the broad Schottky contribution rather than by a sharp transition. Without a quantitative criterion or a corroborating probe, the interpretation of the low-temperature phonon anomaly as an elusive ordering transition is overclaimed.","section":"Results and Discussion, phonon anomalies and Gd-ordering claim; Figs. 2(f)–(h)"},{"comment":"The field-dependent fit uses g, B_int, and θ as free parameters to describe a single resonance frequency curve. Even setting aside the powder-averaging problem, three free parameters for one smooth curve can accommodate a wide range of B_int–θ combinations, and the reported θ = 70(8)° has a large uncertainty that is not propagated into the later comparison with the ab initio value. The agreement claimed in the abstract and conclusions is therefore weaker than the presentation suggests.","section":"Figure 3(c) and the field-dependent fit"}],"minor_comments":[{"comment":"The term 'pellet' is used without emphasizing that the sample is polycrystalline; this is central to the analysis and should be stated prominently in the main text, not only in the Supplemental Material.","section":"Throughout"},{"comment":"No error bars are shown for the fitted Lorentzian parameters. At minimum, representative uncertainties should be given, especially for the phonon parameters used to claim anomalies.","section":"Figures 2 and 3"},{"comment":"The equation Δ² = Δ_ext² + Δ_int² + 2Δ_extΔ_int cosθ is central but is not numbered or defined carefully; the physical meaning of Δ_ext, Δ_int, and θ should be stated explicitly.","section":"Main text, law-of-cosines formula"},{"comment":"The vertical dashed lines are listed as marking 10, 30, 75, and 110 K, but the text refers to both 75 K and 110 K as Ru-related scales; the assignment of each line should be clarified.","section":"Figure 2 caption"},{"comment":"Reference [32] contains '[URL will be inserted by publisher]'; this placeholder should be resolved before publication.","section":"References"},{"comment":"The data are not publicly available; the central claims rely on fitting procedures, so making the raw conductivity spectra and fit parameters available would materially strengthen reproducibility.","section":"Data availability"}],"recommendation":"reject","confidential_remarks":"The powder-averaging problem is not a presentation issue: it invalidates the paper's headline orientation result and the comparison with the predicted cluster-multipolar angle. The mode assignment and the Gd-ordering interpretation are also under-supported. A revision could potentially salvage the zero-field B_int estimate and the observation of a field-dependent low-energy mode, but the central claims as stated would require single-crystal measurements and substantially more evidence. In my view, these are load-bearing issues that cannot be fixed within the current manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: if this were a single crystal, the story would be genuinely exciting. But the authors state in the Supplemental Material that all measurements were done on a pressed powder pellet, and the analysis that yields θ≈70–79° uses a single-angle Zeeman formula that only applies to one orientation. That is a load-bearing inconsistency, and it undercuts the paper's headline agreement with θ_Gd≈80°.\n\nThe genuinely new and useful part is the measurement itself. No prior THz study of Gd2Ru2O7 exists, and the low-frequency mode that appears below ~30 K and hardens with field is a plausible candidate for a transition within the exchange-split Gd3+ J=7/2 multiplet. At zero external field, the resonance frequency gives an internal exchange field B_int ≈ 5.4 T. That number is nearly independent of crystallite orientation and is a real spectroscopic result worth taking seriously. The temperature evolution around the Ru Néel temperature is coherent, and the integrated spectral weight in the supplement shows the main trends are not artifacts of Lorentzian fitting.\n\nThe soft spots are serious. First, the powder issue. In a random polycrystal, the resonance frequency for each crystallite depends on the angle between B_ext and the local B_int; the observed spectrum should be an inhomogeneous sum spanning roughly 0.05–0.36 THz at 7 T. That is far broader than the narrow Lorentzian they fit. Extracting a single θ from that is only valid if every crystallite is aligned, something they do not demonstrate. The agreement with the predicted 80° is therefore not established. Second, the mode assignment is asserted rather than proven. A g-factor near 2 is consistent with Gd3+, but a magnon or an impurity mode could behave similarly. Third, the phonon anomaly near 10 K is called a genuine ordering transition with no error bars and no order-parameter measurement; it may be a real spin–phonon effect, but the claim goes beyond the evidence. Data not being public makes independent checking harder.\n\nNotice the asymmetry: the zero-field B_int survives the powder critique, but the field-orientation result and the cluster-multipolar claim do not. The paper deserves a serious referee because it reports new data and a falsifiable claim, but the current version needs major revision—ideally single-crystal data or a proper powder-averaging model—before the orientation result can stand.","headline":"New THz data on Gd2Ru2O7, but the multipolar-orientation claim is undermined by treating a powder sample as a single crystal.","tokens_in":11717,"tokens_out":3043,"would_cite":false,"duration_ms":30705,"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":"A single terahertz spectrum of Gd2Ru2O7 exposes the internal magnetic field at the rare-earth site and reveals an ordering transition of the Gd sublattice that bulk thermodynamic probes miss.","keywords":["terahertz time-domain spectroscopy","Gd2Ru2O7","pyrochlore ruthenate","exchange-split rare-earth levels","internal exchange field","spin-phonon coupling","multipolar order","hidden magnetic ordering"],"falsifier":"A zero-field neutron-diffraction or 155Gd Mössbauer study that resolves the Gd sublattice below 10 K and shows no onset of long-range order, or no change in the Gd local environment at the temperature of the sharp phonon anomaly, would falsify the claim that the anomaly is a genuine ordering transition. Alternatively, measuring the low-frequency mode on a single crystal with the field along a known crystallographic axis and finding that the splitting does not follow ν² = νext² + νint² + 2νextνint cosθ with g≈2 would refute the exchange-split-level assignment.","tokens_in":10850,"feed_emoji":"🧲","tokens_out":5531,"duration_ms":49131,"temperature":0.7,"pith_summary":"The paper claims that a single magneto-terahertz spectrum of the pyrochlore Gd2Ru2O7 can simultaneously measure the internal exchange field acting on the Gd3+ rare-earth moments and detect a magnetic ordering transition of the Gd sublattice that bulk specific-heat and susceptibility measurements have missed. By tracking a low-frequency resonance attributed to exchange-split Gd3+ levels, the authors extract an internal field of about 5–6 T whose orientation (70–79°) matches the value predicted for a cluster-multipolar ground state. In parallel, an infrared-active optical phonon shows sharp anomalies in frequency, linewidth, and spectral weight at the established Ru-ordering temperature and at a lower temperature scale; the low-temperature anomaly is interpreted as a genuine Gd-ordering transition, revealed because the phonon is sensitive to spin correlations through spin–phonon coupling. The larger point is a route: terahertz spectroscopy can expose hidden ordering transitions in frustrated magnets where thermodynamic probes are blind.","feed_headline":"Terahertz light uncovers a hidden transition in Gd2Ru2O7","feed_subtitle":"A single spectrum measures the internal field at Gd sites and exposes an ordering that heat-capacity and susceptibility missed.","key_machinery":"The central object is the exchange-split Gd3+ magnetic excitation, whose resonance frequency νm is proportional to the magnitude of the total effective field through the Zeeman-like relation Δ = gμBBeff/h with Beff² = Bext² + Bint² + 2BextBint cosθ. This quadratic (cosine) law is the load-bearing identity: it converts two measured splittings into two independent estimates of Bint and one estimate of θ, the angle between internal and applied fields. The second probe is the 1.4 THz optical phonon, treated as a Lorentz oscillator whose scattering time τp acts as a thermometer for spin-disorder scattering; sharp maxima in τp mark the abrupt suppression of spin fluctuations at the two ordering te","core_discovery":"At low temperature and zero field, Gd2Ru2O7 displays a weak resonance near 0.15 THz that hardens and gains weight under an applied field. The paper identifies this mode as a transition between exchange-split levels of the Gd3+ J=7/2 multiplet, split by an internal molecular field from the antiferromagnetically ordered Ru sublattice. Calibrating the splitting against the 7 T response above the Ru Néel temperature gives g≈2 and an internal field Bint≈5.4 T; combining the zero-field and 7 T splittings through the law Δ² = Δext² + Δint² + 2ΔextΔint cosθ yields an orientation θ≈70–79°, consistent with the θGd≈80° predicted for a cluster-multipolar Gd ground state. Independently, the 1.4 THz optic","pith_inferences":["If the framework generalizes, the same two-channel observation could be used to search for hidden order in other R2Ru2O7 (e.g., Tb, Dy) where rare-earth transitions may likewise be concealed by Schottky backgrounds; the predicted signature is a sharp phonon-lifetime maximum at a temperature where no thermodynamic anomaly appears.","A testable extension is to measure the phonon anomaly's field dependence in finer steps and look for hysteresis or a kink in the anomaly position versus field; a first-order or continuous phase boundary would strengthen the 'genuine transition' interpretation beyond what the current data alone establish.","The single-crystal prediction of the model — that the extracted θ rotates with crystal orientation — could be checked directly, and would independently confirm both the multipolar ground state and the validity of the effective-field picture.","If Bint were found (by, say, high-field measurements beyond 7 T) to be field-dependent, the paper's extraction method would need revision; the clean g≈2 and θ≈80° consistency is the main evidence that it does not."],"forward_implications":["If the identification holds, THz conductivity becomes a quantitative probe of the internal exchange field at rare-earth sites, giving Bint and its orientation from one measurement.","The sharp phonon anomaly provides a thermodynamic-independent marker of rare-earth ordering, applicable to other pyrochlores and frustrated magnets where Schottky anomalies mask transitions.","The θ≈70–79° result independently supports the first-principles prediction of a cluster-multipolar ground state for the Gd sublattice, a state that bulk thermodynamic probes do not clearly detect.","The field shifts of both anomalies (110→75 K for Ru order, 10→30 K for Gd order) map a field–temperature phase diagram of the two sublattices in a single sample.","Because both magnetic and lattice responses appear in the same spectrum, the framework ties spin–phonon coupling directly to the magnetic energy scales, enabling tests of magnetoelastic coupling models."],"fun_headline_variants":["Terahertz spectrum catches hidden Gd ordering in pyrochlore","One THz scan measures internal field, finds missed transition","Gd2Ru2O7's hidden order exposed by exchange-split mode","Single THz spectrum spots both Gd mode and phonon to reveal transition"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire extraction rests on the assumption that the low-frequency mode is a simple two-level transition between exchange-split Gd3+ states of a single J=7/2 multiplet with a fixed g-factor and a field- and temperature-independent internal field Bint; if the mode is a different excitation, or Bint changes with field or temperature, the quoted Bint and θ values, and the conclusion that the phonon anomaly marks a Gd ordering transition, would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Terahertz spectrum catches hidden Gd ordering in pyrochlore","One THz scan measures internal field, finds missed transition","Gd2Ru2O7's hidden order exposed by exchange-split mode","Single THz spectrum spots both Gd mode and phonon to reveal transition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000851,"raw_usage":{"total_tokens":3505,"prompt_tokens":677,"completion_tokens":2828,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":421,"completion_tokens_details":{"reasoning_tokens":2752}},"tokens_in":421,"tokens_out":2828,"duration_ms":18411,"temperature":1.0,"reasoning_tokens":2752,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T02:56:50.463500+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A zero-field neutron-diffraction or 155Gd Mössbauer study that resolves the Gd sublattice below 10 K and shows no onset of long-range order, or no change in the Gd local environment at the temperature of the sharp phonon anomaly, would falsify the claim that the anomaly is a genuine ordering transition. Alternatively, measuring the low-frequency mode on a single crystal with the field along a known crystallographic axis and finding that the splitting does not follow ν² = νext² + νint² + 2νextνint cosθ with g≈2 would refute the exchange-split-level assignment.","supporting_citations":[],"review_version":1}