{"id":"7c8b8c16-1240-40dc-b42a-2ad6516b4add","arxiv_id":"2607.22742","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The change of a diffusion activation barrier with temperature equals the variance of site energies, minus the variance of transition-state energies, plus a kinetic correlation term (Eq. 5).","lead":"The paper derives a formula that splits the temperature dependence of diffusion barriers in rough energy landscapes into three fluctuation terms: energy spread, transition-state energy spread, and kinetic correlations. A generalist should care because it offers a way to predict when complex alloys or oxides deviate from Arrhenius behavior, which matters for battery and reactor materials.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 5 is exact only under T-independent energies/prefactors; the abstract's 'broadly applicable' overreaches because real materials have thermal expansion and vibrational entropy that add unaccounted derivative terms.","rationale":"The paper's mathematical core is a clean, parameter-free derivation of Eq. 5. I verified the key steps in Sec. S1 C, including the variational identity that converts the ETS-weighted correlation term into the Δdβη^a positive fluctuation. The derivation is self-consistent and the Sm2O3 example provides an exact demonstration. The most load-bearing weakness is the temperature-independence assumption, which the paper states explicitly but does not elevate to a boundary condition in the abstract and conclusions. The reader's weakest_assumption identifies exactly this issue, so I agree. Since the paper already receives CONDITIONAL due to this and the validation concerns, my stress-test does not change the verdict. The proposed concrete test would quantify the severity of the omitted T-derivative terms and could move the verdict to REJECT if the additional terms turn out to be large in the target systems, but the current evidence is insufficient to do so.","tokens_in":26955,"tokens_out":12312,"duration_ms":102114,"concrete_test":"Implement a minimal two-state jump model with explicit temperature-dependent state and barrier energies, e.g. E_n(T)=E_n0 + α_n T and E_TS(T)=E_TS0 + α_TS T, and a T-dependent jump distance δx(T). Solve the master equation exactly for D(T), compute Q=−d ln D/dβ and dQ/dβ numerically at several T. Then evaluate Eq. 5 using the T-dependent energies at each T as if they were constants. If the numerical dQ/dβ differs from Eq. 5 by terms proportional to α_n, α_TS, and dδx/dT, the omitted-derivative concern is confirmed. A positive result would require the paper's claims to be limited to models satisfying the T-independence assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The derivation of Eq. 5 is mathematically sound under the stated assumptions: with T-independent E_n, E_TS, ρ's, and δx, the variational identities in Sec. S1 C correctly produce dQ^a/dβ = ΔE_n − Δκ,a E_TS + Δdβη^a. However, the abstract's claim of a 'broadly applicable statistical framework' goes beyond this restricted setting. The paper itself acknowledges (Sec. 1) that these quantities are generally T-dependent, but it does not quantify the omitted terms. In real solids, thermal expansion changes lattice parameters and jump distances δx, and vibrational free energies (quasiharmonic phonons) make E_n and E_TS T-dependent. Under such T-dependence, d(P0W)/dβ in Eq. S10 acquires extra terms proportional to dE_n/dβ, dE_TS/dβ, dρ/dβ, and dδx/dβ, and η itself becomes T-dependent through W even at fixed state space. These extra contributions do not reduce to the three fluctuation terms in Eq. 5; they constitute an additional, unquantified source of non-Arrhenius behavior. Thus the central claim — that non-Arrhenius diffusion is governed exactly by these fluctuations — is confirmed only within the idealized T-independent model, not for the 'complex solids' the abstract targets. This is a scope limitation, not an internal inconsistency, but it should be stated as a hard boundary and tested if the framework is to be applied to real materials.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives an exact closed-form expression for the temperature derivative of the per-defect diffusion activation barrier Q^a within harmonic transition state theory, starting from Trinkle's variational principle for mass transport. Defining a kinetic distribution Pκ,a over transition states, the authors obtain dQ^a/dβ = Var(E_n) - Var_κ(E_TS) + Δ_{dβη}, i.e., a decomposition into the equilibrium variance of state energies, the kinetic variance of transition-state energies, and a correlation term involving the temperature derivative of the relaxation vector. The theory is illustrated with two examples: an exact Green's-function calculation for oxygen-vacancy diffusion in Sm2O3, and a machine-learning-based (NN+SRBC) study of vacancy diffusion in Ni-20%Cu. The paper concludes that approximate Arrhenius behavior at high temperature arises from cancellation of thermodynamic and kinetic fluctuations, and that non-Arrhenius behavior appears when ordering and kinetic trapping become important.","tokens_in":27351,"tokens_out":8945,"duration_ms":90372,"significance":"If the central derivation is correct, Eq. (5) is a valuable result: it converts a formal derivative of the diffusion coefficient into three separately interpretable and computable statistical quantities. The Sm2O3 example is a genuinely exact benchmark with parameter-free inputs from published barriers, and it cleanly demonstrates how a kinetic trap can produce a nearly Arrhenius but non-trivial effective barrier. The Ni-Cu example shows how the formalism can be exported to a realistic disordered alloy using machine-learned surrogate models. The main limitations are the explicitly assumed temperature-independence of the energy landscape (energies, entropies, and jump vectors) and the restriction to per-defect diffusion; both limit the 'broadly applicable' claim made in the abstract. The alloy validation is also weaker than it appears because the extrapolations and reference data derive from the same surrogate.","major_comments":[{"comment":"The derivation of Eq. (5) is sound only under the assumption stated on p. 2 that ρ_n,n', E_TS, ρ_n, E_n, and δx are temperature-independent. When these quantities depend on T—thermal expansion, vibrational entropy, phonon renormalization—the derivative d(P0W)/dβ in Eq. S10 acquires additional terms, and Eq. (5) is not the complete dQ^a/dβ. The paper acknowledges this in Sec. 1, but the abstract and conclusions claim a 'broadly applicable statistical framework' for 'complex solids'. This is a load-bearing scope mismatch. Please either restrict the claims to the T-independent landscape model, or use the formalism of Sec. S3 (which already treats λ-dependent free energies) or an explicit quasiharmonic estimate to quantify the omitted derivative terms.","section":"Sec. 1 / Eq. (5)"},{"comment":"The Ni-Cu demonstration does not validate Eq. (5) or the 'nearly Arrhenius at high temperature' conclusion. The first- and second-order extrapolations in Fig. 3 are constructed from Q_v and dQ_v/dβ computed with the NN+SRBC surrogate, and the reference points being extrapolated are the same NN+SRBC estimates. Agreement thus only demonstrates that the truncated Taylor expansion is internally consistent with the surrogate; it is circular with respect to the physics. The subset-resampling uncertainties in Fig. 2(a) are statistical only and do not bound surrogate error. I request either a direct brute-force kMC benchmark for at least one pair of temperatures, or a clear statement that the alloy results are surrogate-model predictions rather than independent validation.","section":"Sec. 3 / Fig. 3"},{"comment":"The paper explicitly sets aside defect formation and annihilation. Consequently, all equations and examples describe per-defect mobility, not the total self-diffusion coefficient that tracer experiments measure. The comparison to the experimental tracer activation barrier of CoNiCrFeMn in Sec. 3 is indirect: an equilibrium vacancy concentration contributes its own Arrhenius temperature dependence. This restriction should appear in the abstract and conclusion; as written, 'thermal diffusion in complex solids' and 'vacancy diffusion in solid solutions' overstate the demonstrated scope.","section":"Sec. 1 (last paragraph) and Sec. 3"}],"minor_comments":[{"comment":"The statement that Pκ,a is 'independent of E_n' is too strong: η^a depends on the rates W and hence on the state energies E_n. Recommend rephrasing as 'has no explicit dependence on the equilibrium state probabilities P0_n'.","section":"Sec. 1, after Eq. (4)"},{"comment":"The symbol Δdβη^a is easy to misread as a product. Consider writing it as Δ_{d_β η} or adding a one-line definition of the subscript.","section":"Eq. (5)"},{"comment":"In the sentence 'kinetic correlations induced by fast, low-barrier mechanisms that that can hinder diffusion', the duplicated 'that' should be removed.","section":"Sec. 2"},{"comment":"Equation (46) for dΓ/dβ is stated without derivation; a one-line derivation would help the reader verify the sign.","section":"Sec. S1 D"},{"comment":"Please clarify whether the vacancy participates in the Metropolis swap moves, and whether the reported state energies include the vacancy. The fixed vacancy concentration is stated, but the MC protocol is otherwise ambiguous.","section":"Sec. S2 A"}],"recommendation":"major_revision","confidential_remarks":"The central analytical derivation appears correct and publishable, but the authors should not claim broad applicability beyond the T-independent, per-defect model without either extending the formalism or softening the claims. The alloy section needs an independent kMC check or a clear demotion to an illustrative surrogate-based prediction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: Eq. 5 and the Pκ weighting are worth reading. The derivation is straightforward in the best sense — it is a rearrangement of derivatives from Trinkle's variational formalism, but it is explicit, self-contained in the SI, and not fitted. The Sm2O3 example is genuinely clean: all quantities follow from literature barriers and an exact Green's-function method, parameter-free, and it shows the kinetic trap effect and the competing-fluctuations mechanism clearly. That is real credit.\n\nThe Ni-Cu alloy demo is the soft spot. The NN+SRBC surrogate is reused from [26] without an independent brute-force kMC cross-check, and Fig. 3 compares first- and second-order extrapolations with the surrogate values from which those extrapolations were derived. The 'nearly Arrhenius at high temperature' conclusion is partly self-confirming. Error bars in Fig. 2 are only one-tenth SD over subsets, not real uncertainties in the extrapolations. Data and code are not shipped; the neural network source is available from [26], but the exact data and workflows from this paper are not.\n\nThe stress-test concern is fair. Section 1 explicitly assumes temperature-independent ρ, E_TS, E_n, and δx, and then the abstract calls the framework 'broadly applicable.' If thermal expansion or vibrational entropy makes any of those T-dependent, extra derivative terms enter the Poisson equation for dη/dβ and the decomposition in Eq. 5 is incomplete. The paper acknowledges the assumption but does not quantify the neglected terms, so the 'non-Arrhenius diffusion is governed by these three fluctuations' claim is true within the idealized model, not proven for real complex solids. On top of that, the examples are per-defect diffusion; defect formation and annihilation kinetics are set aside, so the link to measured diffusivity is indirect. These are scope limits, not internal contradictions, and one of them (T-independence) is acknowledged in the text. Still, the abstract oversells.\n\nMy verdict: this deserves a serious referee. The core result is correct under its assumptions, the decomposition is new enough, and the Sm2O3 example is a compact parameter-free demonstration. The revision should add an independent kMC benchmark for the alloy, release data, and rewrite the abstract to state the T-independent boundary. The circularity in Fig. 3 does not sink the paper because the central identity does not depend on that figure; it limits the empirical support. If the authors add the benchmark and soften the framing, this becomes a solid methods paper for the diffusion-in-disordered-solids crowd.","headline":"A clean, parameter-free decomposition of dQ/dβ into thermodynamic and kinetic fluctuation terms, with a good exact example and a weaker alloy validation; the T-independence assumption is the real limit on 'broad applicability.'","tokens_in":27809,"tokens_out":1890,"would_cite":true,"duration_ms":18329,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["66.30.-h","05.20.-y"],"model":"deepseek-v4-flash","headline":"A new identity splits the temperature dependence of diffusion barriers into three computable fluctuations.","keywords":["diffusion","disordered materials","variational principle","machine learning","non-Arrhenius behavior","activation barrier","kinetic distribution","transition state theory"],"falsifier":"On a lattice model with temperature-independent barriers, compute D(T) by kinetic Monte Carlo, numerically differentiate ln D to obtain dQ/dβ, and compare with the right-hand side of Eq. 5 evaluated from the same jump statistics; agreement at all temperatures confirms the identity. Repeating with barriers that depend on temperature, for example through vibrational entropy, should produce a discrepancy proportional to those derivatives.","tokens_in":26831,"feed_emoji":"🌡️","tokens_out":4165,"duration_ms":37265,"temperature":0.7,"pith_summary":"The paper claims that non-Arrhenius diffusion in rough energy landscapes is governed by three microscopic fluctuation terms: the variance of state energies, the variance of transition-state energies under a kinetics-controlled distribution, and a kinetic-correlation term. It derives an exact identity for the first temperature derivative of the diffusion activation barrier and shows that approximately Arrhenius behavior emerges when these fluctuations compete or cancel. Calculations on oxygen-vacancy diffusion in a bixbyite oxide and vacancy diffusion in a Ni-Cu solid solution demonstrate that the identity can be evaluated from local energetics and machine-learned jump statistics. This matters because it replaces empirical fits with a quantitative, microscopic accounting of why diffusion plots curve.","feed_headline":"Arrhenius curvature traced to three microscopic fluctuations","feed_subtitle":"Diffusion-barrier temperature dependence separates into three computable variance terms, replacing empirical fits.","key_machinery":"The variational principle for mass transport expresses the self-diffusion coefficient as an infimum over trial functions; its minimizer, the relaxation vector η^a_n, encodes correlations between successive jumps and satisfies a discrete Poisson equation. From this, the activation barrier Q^a becomes an average of transition-state energies under the kinetic distribution Pκ, and differentiating that average yields Eq. 5. The relaxation-vector derivative dβη is shown to obey its own Poisson equation and variational principle, making the new kinetic term accessible to machine-learned approximations from single-step kinetic Monte Carlo data.","core_discovery":"The central claim is Eq. 5 of the main text: dQ^a/dβ = ΔE_n − Δ_{κ,a} E_{TS} + Δ_{dβη}. The first-order temperature dependence of the per-defect diffusion activation barrier is exactly three statistical fluctuations: the variance of state energies under the equilibrium Boltzmann distribution, the variance of transition-state energies under the kinetic distribution Pκ (which weights each jump by its squared contribution to net diffusion), and the averaged fluctuations of the temperature derivative of the relaxation vectors across the diffusion network. If this identity holds, non-Arrhenius behavior is not a free fit parameter but a consequence of calculable thermodynamic and kinetic fluctuati","pith_inferences":["If Eq. 5 survives scrutiny, the sign of an Arrhenius plot's curvature becomes a diagnostic: positive curvature would indicate thermodynamic energy fluctuations dominating, negative curvature would indicate kinetic barrier variance, and near-flat behavior would indicate cancellation.","The supplement's general treatment of derivatives with respect to macroscopic variables suggests analogous identities for activation volumes and elastodiffusion tensors, where the same fluctuation decomposition might hold under strain or external fields.","Existing kinetic Monte Carlo datasets could be post-processed to compute all three terms, offering a testable extension to high-entropy alloys where 'sluggish diffusion' would be attributed to specific fluctuations rather than assumed.","The role of configurational entropy near a miscibility gap suggests that the same identity could be used to detect the onset of chemical ordering from diffusion data alone."],"forward_implications":["Deviations from Arrhenius behavior in complex solids can be computed from local energies, barriers, and jump statistics rather than inferred from fits.","Approximately Arrhenius diffusion at high temperature in solid solutions follows from configurational entropy and cancellation between the thermodynamic and kinetic variance terms; non-Arrhenius behavior appears as short-range ordering sets in.","Kinetic traps, such as repeated back-and-forth jumps, raise the effective activation barrier at low temperature, while competing variance in transition-state energies lowers it; the competition limits the curvature of the Arrhenius plot.","The relaxation-vector derivatives satisfy their own variational principle, so the new kinetic term can be approximated with physics-conforming neural networks and linear bias corrections trained on kinetic Monte Carlo data.","The identity generalizes superbasin transition-state theory and shows that transitions with transition-state energies below the thermal average can contribute negative terms to the activation barrier."],"fun_headline_variants":["Diffusion barrier's hidden math: exactly three fluctuations","Non-Arrhenius behavior comes down to three fluctuations","Exactly three fluctuations set the diffusion barrier slope","Diffusion barrier's temperature slope: a sum of three variances","Non-Arrhenius diffusion demystified: exactly three fluctuations"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"State energies, transition-state energies, entropic prefactors, and jump displacement vectors are assumed temperature-independent, and defect formation and annihilation kinetics are set aside; if thermal expansion, vibrational entropy, or temperature-dependent barriers matter, Eq. 5 misses extra derivative terms.","fun_headline_variants_meta":{"raw":{"variants":["Diffusion barrier's hidden math: exactly three fluctuations","Non-Arrhenius behavior comes down to three fluctuations","Exactly three fluctuations set the diffusion barrier slope","Diffusion barrier's temperature slope: a sum of three variances","Non-Arrhenius diffusion demystified: exactly three fluctuations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000416,"raw_usage":{"total_tokens":1905,"prompt_tokens":585,"completion_tokens":1320,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":329,"completion_tokens_details":{"reasoning_tokens":1254}},"tokens_in":329,"tokens_out":1320,"duration_ms":8608,"temperature":1.0,"reasoning_tokens":1254,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T09:29:22.305854+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"On a lattice model with temperature-independent barriers, compute D(T) by kinetic Monte Carlo, numerically differentiate ln D to obtain dQ/dβ, and compare with the right-hand side of Eq. 5 evaluated from the same jump statistics; agreement at all temperatures confirms the identity. Repeating with barriers that depend on temperature, for example through vibrational entropy, should produce a discrepancy proportional to those derivatives.","supporting_citations":[],"review_version":1}