{"id":"22e415ec-702e-4bab-8f44-7af5a5514f69","arxiv_id":"2608.00469","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"TbFeO3 shows a tiny volume increase on cooling, phonon anomalies beyond anharmonic decay, and a new Raman mode emerging below about 175 K, attributed to spin-lattice coupling without a structural transition.","lead":"This paper studies how the crystal, magnetism, and atomic vibrations of TbFeO3 change with temperature from 5 to 300 K. It finds a tiny expansion on cooling and a new vibration appearing below about 175 K, pointing to hidden links between magnetism and the crystal lattice.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"NTE claim rests on a 0.017% volume change that is at the mercy of uncorrected temperature-dependent XRD systematics; no calibration standard or error budget is reported.","rationale":"The reader's weakest assumption correctly identifies the principal vulnerability: the NTE signal is 0.017% in volume, which is smaller than the systematic uncertainties typically present in laboratory temperature-dependent powder XRD unless explicitly calibrated. This is not an internal mathematical inconsistency, but a correctness risk from an incomplete error budget. The paper gives no internal standard, no displacement parameter refinement, no comparison of cooling versus warming, and no replicate measurement; the quoted error bars are Rietveld statistical uncertainties. The NTE is the headline finding and the basis for the broader claim of coupled lattice-magnetic anomalies, so if this concern lands, the central claim is weakened substantially. Other issues, such as the 'first time' wording conflicting with Section III.A's 'established experimentally' and the speculative origin of the 175 K mode, are secondary. I would keep the reader's CONDITIONAL verdict, with the condition being a quantitative XRD error analysis and/or internal-standard calibration.","tokens_in":19206,"tokens_out":2866,"duration_ms":35475,"concrete_test":"Mix the TbFeO3 powder with NIST SRM 640e silicon internal standard, measure temperature-dependent XRD on the same instrument from 5 to 300 K, and Rietveld-refine both phases while freely varying the sample-displacement parameter. If the TbFeO3 volume still increases by ≥0.017% on cooling while the Si lattice parameter follows its known contraction, the NTE claim is supported. Also repeat the measurement on a second sample mount and in both cooling and warming directions to rule out holder drift and thermal hysteresis.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central novelty is the negative thermal expansion claim: a +0.017% unit-cell volume change between 303 K and 5 K (Section III.A, Fig. 1(e)). This is an average linear strain of roughly +0.006%, i.e. about 0.0003 Å on the a-axis. The error bars shown in Fig. 1(e) are described only as 'standard uncertainties obtained from the Rietveld refinement' (figure caption), which capture statistical fitting error, not systematic errors of the temperature-dependent powder XRD measurement. Low-temperature XRD setups are susceptible to temperature mis-calibration, thermal contraction of the sample holder, slow drift in sample surface height, zero-point offset, and detector alignment drift; in parallel-beam geometry some height errors are suppressed, but the paper provides no calibration standard, no repeated reproducibility check, and no refinement or reporting of a sample-displacement parameter as a function of temperature. The room-temperature agreement with published lattice parameters (Table I) validates the absolute scale at one temperature, not the temperature derivative. Because the claimed effect is only 0.017%, an uncontrolled systematic error at the 0.01% level would erase or reverse the conclusion. The paper's own Section IV calls the NTE report 'for the first time,' which raises the evidentiary bar further. This is the load-bearing concern: the 'pronounced interplay' conclusion depends on the NTE being real, and the current XRD evidence does not yet establish it at the required precision.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a temperature-dependent (5–300 K) study of polycrystalline TbFeO3 combining powder XRD, DC magnetization, XPS, and Raman scattering. The main claims are: (i) a small negative thermal expansion of the unit-cell volume (~0.017% increase on cooling from 303 K to 5 K) with no structural phase transition; (ii) mixed Tb3+/Tb4+ and Fe3+/Fe2+ surface states suggesting oxygen non-stoichiometry; (iii) Raman-active phonon modes whose temperature dependence deviates from the Klemens anharmonic-decay model, interpreted as spin-phonon coupling; (iv) two-magnon excitations near 960 and 1108 cm−1, assigned with the help of linear spin-wave theory using exchange parameters from Ref. [27]; and (v) a broad Raman mode near 206 cm−1 emerging below ~175 K with an order-parameter-like intensity increase. The paper concludes that TbFeO3 shows pronounced interplay among lattice, spin, and local magnetic degrees of freedom.","tokens_in":19641,"tokens_out":5900,"duration_ms":71840,"significance":"If the NTE claim is correct, it would be a new, small negative-thermal-expansion effect in an orthoferrite over a wide temperature range, and the combined structural/magnetic/spectroscopic dataset would be a useful reference. The paper is commendably explicit about several limitations: the spin-phonon coupling analysis is described as qualitative, the power-law exponent for the 206 cm−1 mode is called an effective fitting parameter, and the two-magnon comparison is acknowledged in the SI as qualitative. These statements are in the paper itself and reduce the risk of overclaiming. However, the central NTE conclusion currently rests on an XRD analysis with no calibration or systematic error budget, and several secondary conclusions are built on qualitative fits. The paper is a good candidate for major revision: the multi-technique framework is appropriate, but the evidence for the headline result needs to be strengthened before the conclusions can be accepted.","major_comments":[{"comment":"The NTE claim rests on a 0.017% volume increase between 303 K and 5 K. The error bars shown are described only as standard uncertainties from Rietveld refinement, i.e., statistical fitting errors. No calibration standard (e.g., NIST SRM 640/660), no sample-displacement refinement versus temperature, and no reproducibility check are reported. An uncorrected systematic effect at the 0.01% level (sample-holder contraction, height drift, zero-point offset, alignment drift) would be comparable to or larger than the claimed effect. The authors should add a temperature-dependent calibration measurement, report the refined displacement parameter as a function of T, or otherwise quantify the systematic floor. Without this, the statement in Section IV that NTE is established 'for the first time' is not supported.","section":"Section III.A, Fig. 1(e)"},{"comment":"The claim of spin-phonon coupling is based on deviations of phonon frequencies from Klemens-model fits, but the deviations are not quantified. The red dotted curves are extrapolations of fits to high-temperature data, yet the paper does not report fit residuals, uncertainties in the fitted parameters (ω0, C, Γ0, Γ), or a comparison against an alternative model. Since some deviations are only a few cm−1, a quantitative Δω(T) = ω_exp − ω_Klemens plot with error bars is needed to substantiate the 'clear deviations' stated in the text. The authors explicitly describe the discussion as qualitative, but the conclusion of spin-phonon coupling is load-bearing for the paper's central message.","section":"Section III.E, Fig. 8"},{"comment":"The two-magnon assignment is supported only by a qualitative overlap between the calculated 2×MDOS and two broad experimental features at ~960 and ~1108 cm−1. The calculation does not include two-magnon Raman matrix elements, exchange-striction vertices, or orientational averaging for a powder, and the exchange parameters are imported from Ref. [27] rather than determined here. The SI already disclaims this comparison as qualitative, but the main text presents the two-magnon assignment as established and later uses W2M as a spin-correlation proxy. Please either include a more direct two-magnon scattering calculation or clearly label the assignment as tentative in the main text.","section":"Section III.D, Fig. 5 and inset"},{"comment":"The emergent 206 cm−1 mode is characterized through difference spectra and a power-law fit A(T) ∝ (T*−T)^β with T*≈175 K and β≈0.36. No uncertainties are given for T* and β, and no alternative functional forms (activated, BCS-like, or Gaussian onset) are tested. Because the mode is weak and broad, it is important to show that the difference-spectrum procedure and the background model do not create or distort the feature; raw Bose-corrected spectra and fit residuals should be displayed. The text properly says β is an effective fitting parameter, but the conclusion of an 'order-parameter-like' evolution goes beyond what a single effective power-law fit can establish.","section":"Section III.E, Fig. 9"}],"minor_comments":[{"comment":"The phrase 'for the first time' should be backed by a comparison with previous thermal-expansion or lattice-parameter studies of TbFeO3; as written it is an unsupported novelty claim.","section":"Section IV"},{"comment":"The lattice parameters for this work are quoted without uncertainties, even though the text refers to standard uncertainties from Rietveld refinement. Add the uncertainties and the χ² value in the table.","section":"Table I"},{"comment":"The onset temperatures for the Gaussian-to-Lorentzian crossover are inconsistent: 'above ~40 K' for the Ag mode in one paragraph and 'near ~50 K' in the next. Please make these values consistent and explain any discrepancy.","section":"Section III.E"},{"comment":"The difference spectrum is plotted as |χ''(ω,T)−χ''(ω,300 K)|. Using an absolute value is unconventional for a difference spectrum and should be justified, since it can create artificial cusps at crossings.","section":"Fig. 9(a)"},{"comment":"Several modes have very large fitting uncertainties (e.g., 257.15±6.18 cm−1 with FWHM 13.60±3.09 cm−1; 403.17±1.32 cm−1 with FWHM 9.20±4.40 cm−1). Please comment on whether these modes are reliably resolved; otherwise the comparison with literature is difficult to evaluate.","section":"Table IV"},{"comment":"The sentence 'the B1g mode, involving Tb atomic displacements along the a direction, gaussian nature aligns with...' is grammatically incomplete and should be revised.","section":"Section III.E"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid multi-technique characterization, but the headline NTE result is not yet established because of missing XRD calibration and systematic error analysis. If the authors can provide the requested calibration/reproducibility data, or reframe the NTE claim as 'apparent NTE' pending such data, the paper could be acceptable. The 'first time' claim should also be checked editorially against prior thermal-expansion work on TbFeO3. The Raman and magnetization parts are more robust and could support a revised manuscript even if the NTE claim is softened."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid, honest multi-probe study of polycrystalline TbFeO3 that gives you a useful temperature map of lattice, magnetic, and phonon behavior. The genuinely new pieces are the reported negative thermal expansion (0.017% volume increase on cooling from 303 to 5 K) and a broad Raman mode at ~206 cm-1 that grows in below about 175 K. The paper does a lot right: the Raman analysis is careful, the Klemens-model deviations are shown mode-by-mode, the two-magnon assignment uses literature exchange parameters in SpinW as a check rather than a fit, and the authors are frank about what is qualitative—β is an effective parameter, T0 is unresolved, the spin-phonon coupling constants are not extracted, and the 175 K feature's origin is left open. That candor is real and makes the paper easier to trust than most.\n\nThe soft spots are in proportion: the NTE is the load-bearing novelty and it rests on a 0.017% volume change with no calibration standard, no stated systematic error budget, and no temperature-dependent sample-displacement correction. Rietveld \"standard uncertainties\" only capture statistical error; the room-temperature lattice parameters matching literature validates the absolute scale, not the temperature derivative. Uncorrected systematics at the 0.01% level would erase the effect. The paper also undermines itself by writing in Section III.A that NTE \"is established experimentally\" and in Section IV that it is reported \"for the first time\"—those need to be reconciled, and the prior literature checked. Two smaller things: the 206 cm-1 mode is tracked via difference spectra, which can create apparent peaks from shifts of nearby phonons, so the raw spectra should be shown more prominently; and the calculated two-magnon DOS peaks near 1056 cm-1 while the observed features are at 960 and 1108 cm-1—calling that \"close proximity\" is charitable, though the SI does flag the comparison as qualitative.\n\nOverall: no load-bearing math error, no obvious fabrication, and the central claim of coupled lattice-spin anomalies is plausible but not yet proven because the NTE evidence is too thin. The paper deserves a serious referee; a good referee should ask for an error budget and a reproducibility check on the XRD, a resolution of the novelty inconsistency, and a robustness check on the 206 cm-1 mode. If the NTE survives that, this is a useful reference for TbFeO3 and related orthoferrites. I'd send it to review.","headline":"A useful, honest multi-probe temperature map of TbFeO3, whose most novel claim—negative thermal expansion—rests on a 0.017% volume change that needs a real error budget before it can be believed.","tokens_in":20121,"tokens_out":2883,"would_cite":true,"duration_ms":33202,"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":"TbFeO3 expands slightly as it cools from 300 K to 5 K, and its Raman spectrum shows a new spin-coupled mode near 206 cm^-1 below ~175 K.","keywords":["TbFeO3","orthoferrite","negative thermal expansion","spin-phonon coupling","two-magnon Raman scattering","antiferromagnetism","mixed valence","Klemens anharmonic model"],"falsifier":"A high-resolution diffraction measurement over 5-300 K on the same sample using an internal standard, with explicit correction for sample displacement and sample-holder thermal expansion, would settle the NTE claim: if the volume change disappears or reverses sign, the negative thermal expansion is not supported. A complementary check is capacitive dilatometry on a dense pellet.","tokens_in":19108,"feed_emoji":"🧲","tokens_out":5973,"duration_ms":69232,"temperature":0.7,"pith_summary":"This paper argues that TbFeO3, a rare-earth orthoferrite, shows a small but genuine negative thermal expansion: its crystal lattice expands as it is cooled from 300 K to 5 K, even though the crystal structure never changes. The same study ties this lattice behavior to magnetism: several Raman phonon modes deviate from the standard anharmonic decay model, two modes switch from Gaussian to mixed line shapes, and a new broad Raman mode appears below about 175 K with order-parameter-like growth. The authors also identify two-magnon excitations from the iron sublattice, matching linear spin-wave calculations, and find that these magnetic excitations stay nearly constant while the lattice anomalies develop. If correct, the work makes TbFeO3 a concrete example of coupled spin-lattice physics in which thermal expansion, local disorder, and magnetism respond together.","feed_headline":"Cooling makes TbFeO3 expand — a new Raman mode appears at 175 K","feed_subtitle":"XRD, magnetization, and Raman data tie the lattice expansion to spin-phonon coupling in an orthoferrite.","key_machinery":"The argument is carried by four interlocked measurements: Rietveld-refined temperature-dependent XRD supplies the volume anomaly; Voigt line-shape analysis with the Thomas-Cox-Hastings pseudo-Voigt parameter separates Gaussian (inhomogeneous) from Lorentzian (lifetime) broadening; the Klemens decay model provides the anharmonic baseline from which phonon frequency deviations are measured; and linear spin-wave theory converts the Fe-sublattice magnon density of states into the predicted two-magnon Raman response. The coupling constant lambda, defined by the dependence of exchange on atomic displacement, is the conceptual link between the phonon anomalies and magnetic correlations.","core_discovery":"On its own terms, the paper establishes that polycrystalline TbFeO3 exhibits negative thermal expansion over 5-300 K—a total volume increase of about 0.017% on cooling—with no structural phase transition. Raman scattering reveals that several phonon modes deviate from the Klemens anharmonic decay model, which the authors interpret as spin-phonon coupling, and that two modes of Ag and B1g symmetry cross over from Gaussian-dominated low-temperature line shapes to mixed Gaussian-Lorentzian profiles at higher temperature, indicating a change from inhomogeneous broadening to lifetime-driven dynamics. High-energy Raman features near 960 and 1108 cm^-1 are assigned to two-magnon scattering from the","pith_inferences":["A direct test of exchange striction as the NTE driver would be to measure thermal expansion across TN ~ 650 K: if the volume anomaly strengthens near the magnetic ordering temperature, spin-lattice coupling is the likely source.","The surface mixed valence seen by XPS makes oxygen stoichiometry a plausible control knob; comparing as-made and oxygen-annealed samples would show whether the NTE and the 175 K mode are intrinsic or vacancy-enhanced.","Because the ~206 cm^-1 mode grows without a two-magnon anomaly, it may signal short-range or glassy spin reconfiguration rather than a thermodynamic transition; local probes such as muon spin rotation or neutron pair-distribution analysis could distinguish these possibilities."],"forward_implications":["If confirmed, TbFeO3 becomes a candidate component for compensating positive thermal expansion in composite materials.","The ~175 K onset of the new Raman mode becomes a benchmark temperature for spin-lattice coupling in orthoferrites, even though no long-range structural or magnetic transition occurs there.","The two-magnon assignment provides a basis for future neutron or optical studies of magnon-phonon hybridization in the Pbnm orthoferrite family.","The Gaussian-to-Lorentzian crossover shows that line-shape analysis, not only peak position, can expose local magnetic or structural disorder in polycrystalline correlated oxides."],"supporting_citations":[{"why":"Supplies the single-crystal Raman mode frequencies and symmetry assignments used to identify the polycrystalline modes.","marker":"[19]"},{"why":"Supplies first-principles phonon assignments and the ~650 cm^-1 upper cutoff used to rule out a first-order origin for high-energy Raman features.","marker":"[51]"},{"why":"Supplies the Fe-sublattice exchange, anisotropy, and DM-interaction parameters used in the linear spin-wave calculation of the two-magnon density of states.","marker":"[27]"},{"why":"Provides the Rietveld refinement method used to extract lattice parameters and unit-cell volume from the XRD patterns.","marker":"[23]"},{"why":"Provides the Klemens anharmonic decay model used as the baseline for phonon frequency and linewidth fits.","marker":"[60]"},{"why":"Supplies the Thomas-Cox-Hastings line-shape formalism used to quantify the Gaussian versus Lorentzian broadening weights.","marker":"[57]"},{"why":"Supplies the low-temperature Fe and Tb magnetic structures used to interpret the magnetization and spin-correlation behavior.","marker":"[43]"},{"why":"Reports earlier deviations from the Klemens model in TbFeO3 that the present work extends with combined structural and magnetic data.","marker":"[18]"}],"fun_headline_variants":["TbFeO3 expands when cooled — new Raman mode at 175 K","Spin-phonon coupling drives negative thermal expansion in TbFeO3","Cooling expands TbFeO3, reveals hidden magnetic-lattice mode","New Raman mode signals magnetic-lattice interplay in TbFeO3","TbFeO3: negative thermal expansion tied to spin-lattice coupling"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The claimed 0.017% volume increase on cooling is larger than systematic errors in the powder XRD experiment, such as sample displacement, holder expansion, and instrument drift.","fun_headline_variants_meta":{"raw":{"variants":["TbFeO3 expands when cooled — new Raman mode at 175 K","Spin-phonon coupling drives negative thermal expansion in TbFeO3","Cooling expands TbFeO3, reveals hidden magnetic-lattice mode","New Raman mode signals magnetic-lattice interplay in TbFeO3","TbFeO3: negative thermal expansion tied to spin-lattice coupling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000173,"raw_usage":{"total_tokens":1143,"prompt_tokens":797,"completion_tokens":346,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":541,"completion_tokens_details":{"reasoning_tokens":245}},"tokens_in":541,"tokens_out":346,"duration_ms":4165,"temperature":1.0,"reasoning_tokens":245,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T00:54:51.997796+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-resolution diffraction measurement over 5-300 K on the same sample using an internal standard, with explicit correction for sample displacement and sample-holder thermal expansion, would settle the NTE claim: if the volume change disappears or reverses sign, the negative thermal expansion is not supported. A complementary check is capacitive dilatometry on a dense pellet.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the single-crystal Raman mode frequencies and symmetry assignments used to identify the polycrystalline modes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies first-principles phonon assignments and the ~650 cm^-1 upper cutoff used to rule out a first-order origin for high-energy Raman features."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Fe-sublattice exchange, anisotropy, and DM-interaction parameters used in the linear spin-wave calculation of the two-magnon density of states."},{"cited_title":"Rodríguez-Carvajal, Satellite Meeting on Powder Diffraction of the XV Congress of the IUCr127, 127 (1990)","cited_arxiv_id":null,"evidence_quote":"Provides the Rietveld refinement method used to extract lattice parameters and unit-cell volume from the XRD patterns."},{"cited_title":"Granado, A","cited_arxiv_id":null,"evidence_quote":"Provides the Klemens anharmonic decay model used as the baseline for phonon frequency and linewidth fits."},{"cited_title":"Thompson, D","cited_arxiv_id":null,"evidence_quote":"Supplies the Thomas-Cox-Hastings line-shape formalism used to quantify the Gaussian versus Lorentzian broadening weights."},{"cited_title":"Artyukhinet al., Nat","cited_arxiv_id":null,"evidence_quote":"Supplies the low-temperature Fe and Tb magnetic structures used to interpret the magnetization and spin-correlation behavior."},{"cited_title":"Vilarinhoet al., Sci","cited_arxiv_id":null,"evidence_quote":"Reports earlier deviations from the Klemens model in TbFeO3 that the present work extends with combined structural and magnetic data."}],"review_version":1}