{"id":"c475c320-42bf-42ca-ae55-76977769cb71","arxiv_id":"2607.23484","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Resistivity of elemental actinides is fit by two parallel channels—Bloch–Grüneisen electron–phonon scattering and an Arrhenius nearest-neighbor hopping term—yielding Debye temperatures consistent with heat capacity.","lead":"A two-channel fit (Bloch–Grüneisen plus Arrhenius) matches published resistivity curves for nine actinide phases and some Pu alloys, and recovers Debye temperatures close to heat-capacity values. It offers a single practical form for a long-messy dataset, but the hopping interpretation is postulated rather than independently proven.","discovery_kind":"new_application","skeptic_critique":null,"referee_report":{"model":"moonshotai/kimi-k3","summary":"The manuscript proposes a single phenomenological model (Eq. 10) for the temperature-dependent resistivity of elemental actinides: two parallel conduction channels, one Bloch–Grüneisen (electron–phonon) and one Arrhenius (identified with nearest-neighbor hopping), combined by reciprocal addition (Eq. 8) rather than Matthiessen's rule. The model is fit to publicly available ρ(T) data for nine elemental actinide phases (Th–Cm) and three δ-Pu alloys, with reported R² typically ≥ 0.999. Fitted Debye temperatures are compared to heat-capacity values; for Np an activation energy Ea = 15.9 meV is reported as the largest among elemental actinides. Global fits across samples (Th, α-Pu, Am) and a cooling/warming stability check (Cm) are presented as internal consistency tests. The author further claims (§2, §4) that Matthiessen's rule is fundamentally inapplicable to single-phase metals because transport \"always occurs in parallel.\"","tokens_in":24189,"tokens_out":3950,"duration_ms":155450,"significance":"Actinide ρ(T) curves are notoriously anomalous and resist first-principles description, so a compact, unified fitting form would be genuinely useful to the community. The paper ships several checkable, falsifiable outputs rather than a single illustrative fit: (i) cross-sample consistency of ρ_sat ≈ 136 µΩcm and Ea ≈ 5 meV for three independent α-Pu datasets (Fig. 8–9), (ii) parameter stability under cooling/warming cycling for Cm (Fig. 13–14), (iii) a concrete prediction for Cm (Θ_D ≈ 180 K) that departs from the one published heat-capacity value and can be tested, and (iv) per-sample Θ_D values benchmarked against independent calorimetry. These reproducible, parameter-explicit fits and the external Θ_D anchor are real strengths. However, the model is descriptive rather than derived, with up to five free parameters per curve, so its value currently rests on demonstrating that the fitted parameters are robust and that the functional form is discriminated from alternatives — which the manuscript does not yet do.","major_comments":[{"comment":"The load-bearing conceptual claim — that Matthiessen's rule (Eq. 1) is 'fundamentally flawed' for single-phase metals and that 'electrical transport always occurs in parallel' (Eq. 8) — is asserted without derivation, citation, or test. Standard Boltzmann transport adds scattering rates for mechanisms acting on the same carrier population; parallel (reciprocal) addition applies to distinct carrier channels (multiband, shunting phases). The empirical success of the parallel form does not by itself establish channel physics, since the parallel-resistor form is a well-known phenomenological saturation descriptor (Refs. 79–80) with no hopping content required. The interpretation of Ea as a physical NNH activation energy (the paper's headline result, e.g., Np Ea=15.9 meV) depends entirely on this assertion. The author must either provide a microscopic justification for a coexisting activated","section":"§2, Eq. 8 and §4"},{"comment":"No robustness or model-selection analysis supports the central fitting claim. With up to five free parameters (ρ0, ρD, ΘD, ρsat, Ea) fit to smooth, monotonic curves, R² ≥ 0.999 carries little discriminatory weight, and several fitted parameters betray ill-conditioning: ρD = (386±98) µΩcm for α′-U (Fig. 6a, ~25% relative error) and ρD = (2.4±1.2) mΩcm for Cm (Fig. 13a, ~50%). The paper criticizes polynomial fits (Eq. 3) precisely for parameter instability, but does not demonstrate its own parameters are better behaved. Needed: residual plots (not just confidence bands), parameter-correlation/covariance reporting, and a statistical comparison (F-test/AIC/BIC) of Eq. 10 against nested alternatives — the Ea=0 parallel-resistor limit (Eq. 11) and a series BG+saturation form — for the cases where Ea>0 is claimed (Np, α-Pu, β-Pu, Cm, δ-Pu alloys).","section":"§3, Figs. 6–15"},{"comment":"The Θ_D validation is presented more strongly than the numbers support. For α-Pu the fitted Θ_D spans 137–196 K across three samples (Figs. 8–9), and 'agreement' is claimed against literature values themselves spanning 116–207 K; with ranges this wide the benchmark is nearly non-falsifiable. For Cm the fitted Θ_D = (180±11) K disagrees with the only published heat-capacity value (121 K, Ref. 101) by ~50%, which is noted but not analyzed. For uranium, both Θ_D = 185 K (α′) and 238 K (α) are claimed to be 'in good agreement' with different literature values. The manuscript should state a quantitative criterion for agreement and present a compact table of fitted vs. calorimetric Θ_D per phase so the reader can judge the validation rather than relying on narrative comparison.","section":"§3.5, §3.8; Fig. 16"},{"comment":"Fig. 15c reports ρ0 = (155±2) µΩcm and ρsat = (155±2) µΩcm — identical values and uncertainties, which is implausible for independent fit parameters and suggests a fit degeneracy or a transcription error. Given that the δ-Pu alloys are offered as independent confirmation of the model (with Ea ≈ 16–18 meV, exceeding the elemental maximum), this panel needs to be checked and the fit reported correctly.","section":"§3.9, Fig. 15c"}],"minor_comments":[{"comment":"Abstract: 'a maximum Arrhenius activation energy (among all actinides) of Ea = 15.9 meV' is contradicted within the paper by the δ-Pu alloy values Ea = 16.1–18.0 meV (Fig. 15). Clarify that the maximum refers to elemental actinides only.","section":"Abstract"},{"comment":"Eq. 9 text: 'ρ_∞ and E_a are free-fitting parameters' — ρ_∞ is undefined; the parameter is ρ_sat. Also in the third line of Eq. 16 the subscript 'global' is missing on E_a.","section":"§2, Eqs. 9, 16"},{"comment":"Figure captions inconsistently report fit quality as 'R' (Figs. 2b, 2c, 6b) vs 'R-square (COD)'; presumably all are R². Typographical issues: 'protoactinium' (Fig. 4 caption), 'fundings' (§3.8), 'particularly successful' should read 'partially successful' (§1, first-principles paragraph).","section":"Figs. 2, 4, 6; §1, §3.8"},{"comment":"Data provenance: several datasets are digitized from published figures (e.g., King and Lee data 'digitized from Figure 17 in Reference [12]', Fig. 8a), and some are sparse (β-Pu, Fig. 10) or non-uniformly sampled (Hall Pa data, Fig. 4a). The digitized datasets and fit scripts should be deposited (e.g., Zenodo) rather than offered 'upon reasonable request'; this would materially strengthen reproducibility.","section":"§3; Data availability statement"},{"comment":"Ref. [47] (NiBi3 superconductivity) appears out of place as a citation for Matthiessen-rule-based resistivity models of actinides (Eq. 1 context); please check. Several 2026-dated references (e.g., [56], [57], [67], [81]) should be verified for accuracy at proof stage.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"Single-author manuscript with a noticeable self-citation pattern ([56, 85, 86, 87, 89, 100], several tangential to actinide transport). The empirical fitting work is competent and the internal consistency checks are above the usual standard for this genre, but the strong anti-Matthiessen framing in §2/§4 is likely to draw justified criticism if published as stated; the revision should be judged mainly on whether that claim is either substantiated or appropriately softened, and on whether the requested model-selection analysis survives scrutiny."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful thing here is simple: one explicit two-channel formula (Bloch–Grüneisen in parallel with an Arrhenius term) fits the publicly available ρ(T) curves for Th through Cm and a few δ-Pu alloys, and the extracted Debye temperatures mostly land near independent heat-capacity values. That systematic coverage is new even though every piece of the formula is old.\n\nWhat the paper does well is the bookkeeping. It stops treating Matthiessen addition as mandatory for single-phase metals, writes the reciprocal sum cleanly, shows high R² fits on the historical datasets, and checks ΘD against the literature. For Th, Pa, U, Am the Arrhenius channel collapses to ordinary saturation (Ea = 0); for Np, α-Pu, β-Pu, Cm and the δ alloys it stays finite, with Np giving the largest Ea ≈ 15.9 meV. Global fits across multiple samples of the same element keep ρsat and Ea stable while residual resistivity absorbs sample differences—exactly what you want from a descriptive model. The figures are inspectable and the citation trail to the old resistivity work is honest.\n\nSoft spots are real but proportionate. Five free parameters per curve give a lot of flexibility; several datasets are sparse or digitized from old plots; and nothing independent (pressure, doping series, Hall, optics) pins the second channel as nearest-neighbor hopping rather than another saturation or correlation form. The claim that dissipation channels “must” be parallel in elemental metals is asserted, not derived. So treat Ea as a useful fit parameter with a hopping story attached, not as a measured activation energy. First-principles failures are noted but not replaced by a microscopic calculation.\n\nThis is for people who actually fit or compare actinide and δ-Pu alloy resistivities—metallurgy, aging, alloy design—not for anyone expecting a new 5f transport theory. The math is elementary and reproducible from the published curves. I would send it to referees; they can demand error bars on the digitized points, a leave-one-out stability check, and a clearer separation of phenomenology from mechanism. Worth engaging if you work in the subfield; skip if you only care about microscopic 5f theory.","headline":"A clean phenomenological unification of actinide ρ(T) that works across the archival set and returns sensible ΘD, but the NNH interpretation is still a fit choice rather than a demonstrated mechanism.","tokens_in":24191,"tokens_out":542,"would_cite":true,"duration_ms":11781,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"A single two-channel model—Bloch–Grüneisen plus Arrhenius hopping—fits temperature-dependent resistivity across nine actinide phases and selected plutonium alloys.","keywords":["actinides","electrical resistivity","Bloch-Grüneisen","Arrhenius","nearest-neighbor hopping","parallel resistor model","Debye temperature","plutonium alloys"],"falsifier":"A high-quality ρ(T) data set for a pure elemental actinide phase that cannot be fit by Eq. 10 (or that returns a Debye temperature grossly inconsistent with heat-capacity values), or a demonstration that an alternative second-channel form (for example variable-range hopping or a pure saturation term without activation) systematically outperforms Arrhenius across the same actinide series.","tokens_in":24380,"feed_emoji":"⚡","tokens_out":1127,"duration_ms":23099,"temperature":0.7,"pith_summary":"Actinide metals have long resisted a unified description of how their electrical resistivity changes with temperature: curves look radically different from element to element, polynomial and power-law fits are unstable, and first-principles calculations often fail even qualitatively. This paper argues that the usual series addition of scattering mechanisms (Matthiessen’s rule) is the wrong starting point for single-phase metals. Instead, two conduction channels run in parallel—one classical electron–phonon channel given by the Bloch–Grüneisen formula, the other a nearest-neighbor hopping channel given by an Arrhenius term—and the measured resistivity is the reciprocal of the sum of their conductances. Applied to public data for thorium through curium (nine phases) and to δ-phase Pu–Ce and Pu–Ce–Ga alloys, the model returns high-quality fits and Debye temperatures that match independent heat-capacity values. The largest fitted hopping activation energy among the pure elements is 15.9 meV in neptunium. If correct, the same simple parallel construction supplies a practical, physically motivated way to analyze resistivity across the actinide series.","feed_headline":"One two-channel model fits actinide resistivity from Th to Cm","feed_subtitle":"Bloch–Grüneisen plus Arrhenius hopping recovers Debye temperatures and a 15.9 meV barrier in neptunium","key_machinery":"The NNH-parallel resistance model (Eq. 10): ρ(T) equals the reciprocal of the sum of the reciprocals of an Arrhenius resistivity and a Bloch–Grüneisen resistivity. It generalizes the classical parallel-resistor (saturation) model by allowing the saturated channel to carry a nonzero activation energy, and it replaces Matthiessen series addition with parallel conduction channels for single-phase metals.","core_discovery":"The temperature-dependent resistivity of elemental actinides (Th to Cm, nine phases) and of selected δ-Pu alloys is accurately described by a nearest-neighbor-hopping parallel-resistance model: one Bloch–Grüneisen channel in parallel with one Arrhenius channel. The model yields Debye temperatures in agreement with heat-capacity data and a maximum Arrhenius activation energy of 15.9 meV for neptunium.","pith_inferences":["If the parallel-channel premise is general, many historical power-law and polynomial resistivity analyses of actinides (and possibly other correlated metals) may need re-reduction with reciprocal rather than additive combination.","The element-to-element variation of fitted Ea (zero for Th, Pa, Am; finite and largest for Np) could be used as a comparative probe of 5f localization or hopping barriers once more complete data sets exist.","A natural next test is whether pressure- or self-irradiation-dependent ρ(T) in the same actinides still collapses onto Eq. 10 with smoothly varying Ea and Θ_D.","Extending the second channel to Mott or Efros–Shklovskii forms (Eq. 14) remains unneeded for present data but is a ready diagnostic if low-T curvature appears in cleaner samples."],"forward_implications":["Debye temperatures can be extracted from actinide resistivity curves with the same model that fits the full temperature dependence, and they should track heat-capacity Θ_D.","Neptunium hosts the largest nearest-neighbor hopping barrier among the pure actinides examined (Ea ≈ 15.9 meV).","δ-phase Pu–Ce and Pu–Ce–Ga alloys are described by the same two-channel form, with roughly constant Θ_D ≈ 175 K and Ea ≈ 16–18 meV at low dopant levels.","Matthiessen-style series fits and many first-principles resistivity curves for actinides are expected to remain qualitatively inadequate until parallel hopping is included.","The same parallel construction is offered as a candidate for other single-phase conductors whose resistivity has resisted unified fitting."],"fun_headline_variants":["Two-channel model unifies actinide resistivity from Th to Cm","Bloch–Grüneisen plus Arrhenius fits nine actinide phases","Parallel BG and hopping channels recover actinide Debye temperatures","Neptunium hits 15.9 meV max barrier in actinide resistivity model","Simple two-path model fits ρ(T) for actinides and δ-Pu alloys"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"In single-phase elemental metals, different dissipation mechanisms must be combined as parallel conduction channels rather than added in series, and the second channel is physically nearest-neighbor hopping with an Arrhenius form.","fun_headline_variants_meta":{"raw":{"variants":["Two-channel model unifies actinide resistivity from Th to Cm","Bloch–Grüneisen plus Arrhenius fits nine actinide phases","Parallel BG and hopping channels recover actinide Debye temperatures","Neptunium hits 15.9 meV max barrier in actinide resistivity model","Simple two-path model fits ρ(T) for actinides and δ-Pu alloys"]},"model":"grok-4.5","effort":"low","cost_usd":0.004594,"raw_usage":{"total_tokens":1380,"prompt_tokens":867,"num_sources_used":0,"completion_tokens":102,"cost_in_usd_ticks":45944000,"prompt_tokens_details":{"text_tokens":867,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":411,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":867,"tokens_out":102,"duration_ms":7603,"temperature":1.0,"reasoning_tokens":411,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-30T21:05:07.614199+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A high-quality ρ(T) data set for a pure elemental actinide phase that cannot be fit by Eq. 10 (or that returns a Debye temperature grossly inconsistent with heat-capacity values), or a demonstration that an alternative second-channel form (for example variable-range hopping or a pure saturation term without activation) systematically outperforms Arrhenius across the same actinide series.","supporting_citations":[],"review_version":1}