{"id":"b3d1254c-2bad-4e52-af04-c146f9bc58b4","arxiv_id":"2606.30978","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"The e-RPA-LDA model computes continuous ion stopping power from cold matter to HED plasmas via RPA dielectric response plus four corrections and nuclear terms, matching NIST/IAEA databases and limited plasma data.","lead":"This paper presents an enhanced RPA-LDA model that computes ion stopping power continuously from cold solids through warm dense matter to high-energy-density plasmas by combining dielectric response with four corrections and nuclear stopping. A smart generalist might read it to understand a physics-based alternative to empirical codes like SRIM for fusion and plasma energy deposition calculations.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Necessity/sufficiency of the four corrections (strong-collision, local-field, binding, Barkas-Bloch) for data agreement without implicit tuning remains unverified.","rationale":"The reader's weakest_assumption pinpoints the exact point where the argument is least secured; the open code makes the proposed test feasible and decisive. No other internal inconsistency (e.g., average-atom vs. full KS-DFT) rises to the same load-bearing level for the stated claim.","tokens_in":1899,"tokens_out":379,"duration_ms":30139,"concrete_test":"Using the released GitHub code, recompute proton stopping powers for Au at 0.1–10 MeV with the base RPA-LDA only, then successively enable each of the four corrections individually; compare each variant to NIST PSTAR values and report the RMS deviation. If any correction changes the deviation by <5% or if enabling all four still requires an overall scale factor to match data, the necessity/sufficiency claim is not supported.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The headline claim requires that the base RPA-LDA dielectric response (from average-atom muffin-tin densities via FAC) plus exactly these four corrections reproduces NIST PSTAR/IAEA tables and plasma benchmarks across Z and conditions with no post-hoc parameters. The abstract asserts the corrections are added to Wang et al.'s model and yield agreement, but does not demonstrate (a) quantitative degradation when any single correction is removed, (b) that the functional forms contain no adjustable coefficients fitted to the validation sets, or (c) that the same fixed corrections remain optimal from cold solids through WDM to HED plasmas. If any correction was selected or scaled to improve the match, the \"physics-based alternative to SRIM\" claim weakens even if the final numbers agree.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript presents an enhanced random-phase-approximation local-density-approximation (e-RPA-LDA) model for ion stopping power valid from cold solids through warm dense matter to high-energy-density plasmas. It starts from the RPA-LDA dielectric response of Wang et al. using average-atom muffin-tin densities from the Flexible Atomic Code, then augments it with four corrections (strong-collision, static local-field, electron-binding, and Barkas-Bloch). The resulting proton stopping powers are claimed to agree with NIST PSTAR and IAEA databases across the periodic table and compounds, reproduce limited plasma data and benchmarks, and are extended to total stopping power by adding nuclear and ionic contributions. An open-source GitHub implementation is provided.","tokens_in":2089,"tokens_out":542,"duration_ms":31549,"significance":"If the central claim holds without implicit tuning, the work would supply a computationally efficient, physics-based framework for stopping powers over a wide range of conditions and Z, serving as an alternative to semi-empirical codes such as SRIM. The average-atom treatment of all electrons, applicability to arbitrary degeneracy, and open-source release with tabulated data are positive features for inertial fusion and HED applications.","major_comments":[{"comment":"The headline claim requires that the base RPA-LDA response plus exactly these four corrections reproduces the NIST/IAEA tables and plasma benchmarks with no post-hoc parameters. The manuscript must demonstrate (a) that the functional forms of the corrections contain no adjustable coefficients fitted to the validation sets and (b) that the same fixed corrections remain optimal from cold solids through WDM to HED plasmas; this is not shown in the abstract or model description.","section":"Abstract and model section"},{"comment":"No quantitative evidence is provided for the necessity and sufficiency of each correction. An ablation study (stopping-power error with each correction removed in turn) is needed to establish that agreement across the full temperature-density range is not achieved by any subset or by compensatory adjustments.","section":"Results and validation sections"}],"minor_comments":[{"comment":"The abstract states agreement with databases but does not report quantitative error metrics (e.g., mean relative deviation, maximum deviation) or the precise range of projectile energies and target conditions used for the comparison.","section":"Abstract"},{"comment":"Notation for the four corrections should be introduced with explicit equations early in the model section rather than by name only.","section":"Model description"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the careful review and constructive comments on the manuscript. We address each major comment below and will revise the paper accordingly.","responses":[{"response":"We agree that the abstract and model section should explicitly document that the four corrections are taken from the literature without any coefficients adjusted to the NIST/IAEA or plasma validation data. The strong-collision, static local-field, electron-binding, and Barkas-Bloch corrections follow fixed functional forms derived in prior works and are applied uniformly. We will revise the abstract and add a short subsection in the model description that lists the exact expressions, confirms the absence of fitted parameters, and states that the identical fixed forms are used from cold matter through HED conditions.","revision_made":"yes","referee_comment":"[Abstract and model section] The headline claim requires that the base RPA-LDA response plus exactly these four corrections reproduces the NIST/IAEA tables and plasma benchmarks with no post-hoc parameters. The manuscript must demonstrate (a) that the functional forms of the corrections contain no adjustable coefficients fitted to the validation sets and (b) that the same fixed corrections remain optimal from cold solids through WDM to HED plasmas; this is not shown in the abstract or model description."},{"response":"We acknowledge that the current manuscript does not include a quantitative ablation study. While the results show progressive improvement from the base RPA-LDA to the full e-RPA-LDA model, we did not systematically remove each correction and recompute errors. We will add an ablation analysis to the results section, reporting mean absolute percentage errors against NIST/IAEA data (and available plasma benchmarks) for the model with each correction omitted individually, across representative cold, WDM, and HED regimes. This will demonstrate the contribution of every term.","revision_made":"yes","referee_comment":"[Results and validation sections] No quantitative evidence is provided for the necessity and sufficiency of each correction. An ablation study (stopping-power error with each correction removed in turn) is needed to establish that agreement across the full temperature-density range is not achieved by any subset or by compensatory adjustments."}],"tokens_in":1555,"tokens_out":465,"duration_ms":41999,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The paper adds four corrections (strong-collision, static local-field, electron-binding, Barkas-Bloch) to the Wang et al. RPA-LDA dielectric response, then folds in nuclear and ionic stopping to produce total stopping powers for protons and alphas. It supplies an open GitHub framework that runs from cold solids through WDM to HED plasmas using average-atom muffin-tin densities from FAC.\n\nWhat stands out is the practical unification: one code handles the full temperature-density range without switching models, and the repo includes tabulated outputs. That is useful for inertial-fusion work where you need consistent ranges across phases. The average-atom treatment also captures all electrons, which is a step beyond some Kohn-Sham DFT limits.\n\nThe soft spot is verification. The abstract states agreement with NIST PSTAR, IAEA tables, and plasma benchmarks, but supplies no quantitative tests showing that each correction is required, that none contain parameters fitted to the validation sets, or that the same fixed forms remain optimal from cold matter to hot plasma. If any correction was scaled to improve the match, the \"physics-based alternative to SRIM\" claim rests on weaker ground. The stress-test concern about necessity and sufficiency is therefore still open until the full derivations and sensitivity checks are examined.\n\nThis is for people who need stopping-power tables or code for HED and fusion target design. A reader already using dielectric-response models will see the incremental value and can test the repo directly.\n\nIt should go to peer review. The topic is relevant, the code is public, and the central claim is falsifiable once the validation details are on the table.","headline":"Incremental RPA-LDA extension with four corrections and open code; agreement claims need full-text checks for tuning.","tokens_in":2607,"tokens_out":396,"would_cite":false,"duration_ms":21742,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"An enhanced RPA-LDA model produces proton stopping powers that match NIST and IAEA databases from cold solids to high-energy-density plasmas.","keywords":["ion stopping power","RPA-LDA","warm dense matter","high-energy density plasmas","dielectric response","average-atom model","proton stopping","inertial fusion"],"falsifier":"A new measurement of proton stopping power in a warm-dense or plasma target whose density and temperature lie outside the published data sets, if it deviates from the model's prediction by more than the stated experimental uncertainty, would falsify the claim that the corrections suffice without tuning.","tokens_in":2786,"feed_emoji":"","tokens_out":808,"duration_ms":35622,"temperature":0.7,"pith_summary":"The paper develops a unified computational model for the energy loss of ions traveling through matter under conditions that range from room-temperature solids to hot dense plasmas. It starts from the random-phase-approximation dielectric response evaluated in the local-density approximation using an average-atom electron density and then adds four explicit corrections for strong collisions, local-field effects, electron binding, and higher-order Barkas-Bloch terms. The resulting electronic stopping powers for protons reproduce tabulated experimental values across the periodic table and for compounds, while also matching the sparse existing measurements and calculations in plasma regimes. The same framework is extended by adding nuclear and ionic stopping contributions to give a continuous total stopping power for protons and alpha particles. Because the method remains computationally tractable at arbitrary temperature and density and is released as open-source code, it supplies a single physics-based description usable in inertial-fusion and high-energy-density simulations.","feed_headline":"Enhanced RPA-LDA matches proton stopping to NIST data from solids to plasmas","feed_subtitle":"Four corrections to the dielectric response yield agreement with cold-matter databases and sparse plasma measurements in one continuous fram","key_machinery":"The RPA dielectric response function in the local-density approximation, evaluated on average-atom electron densities and augmented by the four corrections for strong collisions, local-field effects, binding, and Barkas-Bloch terms.","core_discovery":"The e-RPA-LDA model augments the RPA dielectric response in the local-density approximation, obtained from an average-atom muffin-tin potential solved with the Flexible Atomic Code, by four corrections (strong-collision, static local-field, electron-binding, and Barkas-Bloch). This construction yields proton stopping powers that agree with the NIST PSTAR and IAEA databases for cold matter across elements and compounds and that reproduce the available plasma benchmarks, time-dependent DFT results, and the first measurements of enhanced light-ion stopping, while the addition of nuclear and ionic terms produces a continuous total stopping power from cold solids to hot dense plasmas.","pith_inferences":["If the corrections remain valid at still higher densities or temperatures, the same code could generate stopping tables for unmeasured plasma conditions encountered in new fusion concepts.","Because the framework is open source, it could be coupled directly to radiation-hydrodynamics codes to replace look-up tables with on-the-fly stopping calculations.","The underlying dielectric-response approach might be extended to heavier ions by retaining the same four corrections and average-atom densities."],"forward_implications":["The model supplies a physics-based replacement for semi-empirical codes such as SRIM for cold-matter stopping powers.","It reproduces charged-particle transport workshop benchmarks and time-dependent DFT calculations for plasmas.","Nuclear and ionic stopping contributions can be added consistently to give total energy deposition from cold matter through hot plasmas.","The average-atom treatment remains applicable to both low- and high-Z targets at any degeneracy, supporting inertial-fusion design calculations."],"fun_headline_variants":["e-RPA-LDA with four corrections matches NIST proton stopping data","Unified RPA-LDA provides stopping powers from cold matter to dense plasmas","Four corrections to RPA-LDA reproduce cold and plasma ion stopping data","Open-source framework extends RPA-LDA to total stopping in all conditions"],"cache_read_input_tokens":64,"weakest_assumption_plain":"The four added corrections are both necessary and sufficient to match data over the full temperature-density range without any parameter values adjusted to the validation sets.","fun_headline_variants_meta":{"raw":{"variants":["e-RPA-LDA with four corrections matches NIST proton stopping data","Unified RPA-LDA provides stopping powers from cold matter to dense plasmas","Four corrections to RPA-LDA reproduce cold and plasma ion stopping data","Open-source framework extends RPA-LDA to total stopping in all conditions"]},"model":"grok-4.3","cost_usd":0.006926,"raw_usage":{"total_tokens":3292,"prompt_tokens":828,"num_sources_used":0,"completion_tokens":73,"cost_in_usd_ticks":69262000,"prompt_tokens_details":{"text_tokens":828,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2391,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":828,"tokens_out":73,"duration_ms":25120,"temperature":1.0,"reasoning_tokens":2391,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-01T00:37:55.303349+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A new measurement of proton stopping power in a warm-dense or plasma target whose density and temperature lie outside the published data sets, if it deviates from the model's prediction by more than the stated experimental uncertainty, would falsify the claim that the corrections suffice without tuning.","supporting_citations":[],"review_version":1}