{"id":"39224122-2bf8-4f37-8522-72191de62caa","arxiv_id":"2608.00797","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A channel-mixed RPA stopping model with average-atom orbitals reproduces proton stopping data above the Bragg peak and predicts bound-bound and mixing effects in warm dense plasmas.","lead":"This paper builds a new quantum model for how fast ions lose energy in hot, partially ionized plasmas, combining atomic bound-electron transitions with the usual free-electron response. It matches experiments for protons near and above the Bragg peak, and suggests that bound electrons matter in fusion-relevant tungsten but do not explain a known warm dense carbon puzzle.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Free-free Lindhard assumption is the load-bearing weakness: Eq. (4)'s nonlinear mixing amplifies any error in the f-f channel, and Appendix A shows strong sensitivity to AA continuum states, yet the full AA f-f response is never computed.","rationale":"The reader's weakest assumption—the inconsistent Lindhard free-free response—is also the most load-bearing concern I find. The paper validates against ambient experiments and one WDM carbon point, but the claim of a predictive all-electron tool across wide parameter ranges depends on the f-f channel being trustworthy where free electrons are not plane-wave-like. Appendix A shows that a partial correction to the f-f response causes an unphysical low-velocity peak, demonstrating that the channel is sensitive to the treatment of AA continuum states. Because Eq. (4) mixes channels nonlinearly, the channel-mixing effects in Section III C and the tungsten predictions in Section III D inherit this uncertainty. The concrete test—computing the full AA f-f response for a known resonance case—would directly settle whether Lindhard is adequate. I do not see a more fundamental flaw: the derivations are standard, the code appears to implement the stated equations, and the experimental comparisons are appropriate. Therefore the reader's CONDITIONAL verdict should stand unchanged. My agreement is with the same weakest assumption, though I would frame the issue as spectral distribution rather than the 1% f-sum violation, since Appendix A shows sum-rule preservation is not sufficient.","tokens_in":18034,"tokens_out":13460,"duration_ms":183111,"concrete_test":"Compute the full AA f-f contribution using Eq. (6) with both initial and final tartarus continuum states (not Lindhard) for aluminum at T=100 eV, ρ=2.7 g/cm^3, the case with the quasibound 3d resonance, using a reduced continuum grid (e.g., N_E≈200, l_max≈4) and checking convergence of the f-sum rule. Evaluate stopping from Eq. (1) for proton velocities 1–10 a.u. and compare with the Lindhard-based cmRPA result. If the stopping differs by more than ~5% anywhere above the Bragg peak, the f-f approximation is load-bearing for the claimed regimes; if it differs by less, the central claim is supported. If the full calculation is too costly at full resolution, the same convergence check on a coarse grid would still reveal whether the spectral redistribution is consequential.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of predictive all-electron stopping rests on treating the finite-temperature Lindhard function as the free-free channel inside the mixed RPA denominator of Eq. (4). This denominator is nonlinear in χ_ff, so an error in f-f does not merely add a small correction; it changes the interference between the plasmon and the bound-free/bound-bound channels that Section III C highlights. The paper's own Appendix A is the clearest evidence of sensitivity: replacing the Lindhard DOS with the AA DOS while keeping plane-wave matrix elements (Eq. A1) satisfies the f-sum rule yet produces an 'unphysical peak at low velocity' (Fig. A, panel b). That shows the f-f spectral distribution is not robustly controlled by sum rules alone, and the fully consistent AA f-f polarizability—with both initial and final continuum states in Eq. (6)—is never computed. The 1% f-sum violation in Fig. 1(b,d) is a symptom, not the core issue: even an exact sum rule does not fix the distribution. The ambient validations in Fig. 2 involve small Z_bar and nearly plane-wave-like free electrons, so they do not constrain regimes with quasibound resonances, such as Al at T=100 eV or the HED tungsten cases in Section III D. If the true AA continuum response differs from Lindhard in those regimes, the channel-mixing corrections and the WDM/HED predictions could be miscalibrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a channel-mixed RPA (cmRPA) dielectric function for electronic stopping in partially ionized plasmas, combining a finite-temperature Lindhard free-free response with average-atom (AA) bound-bound and bound-free polarizabilities in a single RPA denominator (Eq. 4). The model is validated against ambient IAEA data and TD-DFT for Al, Fe, Ne, and W, and applied to temperature-dependent Al stopping, HED tungsten, and the Malko et al. warm dense carbon experiment. The authors report good agreement above the Bragg peak, identify non-trivial channel-mixing and bound-bound contributions, and conclude that improved bound-state modeling is unlikely to explain the Malko discrepancy.","tokens_in":18438,"tokens_out":6895,"duration_ms":76943,"significance":"If the central claim holds, cmRPA is an efficient, physically transparent all-electron linear-response stopping model for partially ionized plasmas, with explicit bound-channel transitions and a natural mechanism for bound-free/free-free interference. The paper's strengths are the explicit orbital matrix elements, the broad ambient validation set, the direct comparison with TD-DFT, and the application to ICF-relevant conditions. The main limitation is the inconsistent treatment of the free-free response: the Lindhard assumption is not validated in regimes where AA continuum states deviate strongly from plane waves, and Appendix A shows that a sum-rule-preserving modification of the f-f kernel produces large, unphysical stopping changes. This makes the predictive claim for such regimes currently unsupported.","major_comments":[{"comment":"The free-free Lindhard assumption is load-bearing. Because Eq. (4) mixes chi_ff nonlinearly with chi_bb and chi_bf in the denominator, an error in the f-f spectral distribution is not a small additive correction; it modifies the channel interference that is a central result (Sec. III.C). Appendix A shows that replacing the Lindhard DOS with the AA DOS in a sum-rule-preserving manner creates an unphysical low-velocity peak in stopping (Fig. 12b). This demonstrates that f-sum-rule conservation is insufficient to control the stopping-relevant spectrum. The full AA f-f polarizability, with both initial and final continuum states in Eq. (6), is not computed in this work. The ambient validations in Fig. 2 involve small Z_bar and nearly free-electron-like free states, so they do not constrain the problematic regime, e.g., Al at T=100 eV (Fig. 1c-d) or HED tungsten (Sec. III.D). Please provide a","section":"Sec. II.A, Eqs. (4),(6); Appendix A, Fig. 12"},{"comment":"The width parameter c_Gamma is set to 1 with the statement that varying it changes only computational expense and 'does not change the result,' but no convergence study is shown. Given the continuum sensitivity demonstrated in Appendix A, this insensitivity must be supported numerically. Please provide the stopping power as a function of c_Gamma (or, equivalently, continuum mesh spacing) for at least one representative case, and state the criterion used to choose c_Gamma in the production runs.","section":"Sec. II.A, Eq. (12)"},{"comment":"The paper uses the ion-sphere-restricted radial integrals for all production results, explicitly sacrificing orthogonality between bound and free states. The authors show that the effect on stopping is small for carbon at T=10 eV, but this is a single, relatively low-Z case. For the high-Z, tightly compressed tungsten conditions of Sec. III.D, bound-state radii may be comparable to the ion-sphere radius, and the j0 -> j0-1 correction described in Appendix B is an ad hoc fix. Please provide a targeted check of the ion-sphere truncation error for at least one of the tungsten conditions, or state a criterion for when the truncation is valid.","section":"Appendix B, Fig. 13"}],"minor_comments":[{"comment":"The angular frequency appears as w in Eqs. (6), (10)-(12), and (A1), but as omega elsewhere. Use \\omega consistently.","section":"Notation throughout"},{"comment":"The phrase 'the f-f part only (orange) versus ⟨Z⟩' is unclear; please say 'the free-free contribution normalized to ⟨Z⟩'.","section":"Fig. 1 caption, panel b"},{"comment":"The initial proton beam energy distribution is shown 'in green,' but the RPA+CBC curve is also green; use different colors/linestyles to avoid ambiguity.","section":"Fig. 8"},{"comment":"Define g_AA_DOS(E) and g_ideal_DOS(E) explicitly; currently only their ratio is used, which makes the construction hard to reproduce.","section":"Appendix A"},{"comment":"The journal name 'Npj Comput. Mater.' should be 'npj Computational Materials'.","section":"Reference [36]"},{"comment":"Typographical issue: 'we have a)aluminum...' lacks spaces after the panel labels in several places.","section":"Sec. III.A"}],"recommendation":"major_revision","confidential_remarks":"This is a well-executed study with a clear physical framework and honest acknowledgment of its main approximation. The f-f inconsistency is the weakest point: Appendix A shows that the stopping result can be highly sensitive to the continuum treatment even when f-sum rules are preserved, so the current validation set does not justify the predictive claim in regimes with quasibound continuum resonances. A major revision that benchmarks the f-f channel against a full AA continuum response or TD-DFT at such a condition would substantially strengthen the paper. The c_Gamma insensitivity claim should also be backed by a convergence plot. No concerns about novelty or citation practice; the paper fits the journal's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this is a genuine step forward for ion stopping in partially ionized plasmas. Mixing Lindhard free response with explicit average-atom bound-bound and bound-free transitions in one RPA denominator (Eq. 4) is new for stopping applications, and the paper shows it matters: the difference from Chihara-style unmixed decomposition is non-trivial, and bound-bound contributions appear at low velocity. They validate against IAEA data and TD-DFT across four ambient materials, including the neutral neon case with only bound response, and the agreement above the Bragg peak is credible. The Malko analysis is careful and appropriately hedged.\n\nThe paper is honest about its main soft spot: the free-free channel uses finite-T Lindhard rather than AA continuum states. That inconsistency shows up as ~1% f-sum rule violation and is potentially worse where continuum resonances are strong—their own Appendix A shows a modified DOS plus plane-wave matrix elements produces an unphysical low-velocity peak. Since Eq. (4) is nonlinear in the f-f response, an error there does not just add; it can shift the interference between channels they emphasize in Sec. III C. So the tungsten predictions in the resonant regime are the least secure part of the paper. I'd want to see either a full AA f-f calculation for one resonant case (even if expensive) or a sensitivity test that varies the f-f model and shows the channel-mixing results are stable. That would address the concern more directly than the sum-rule statement.\n\nOther, smaller issues: no code or data released, and no quantitative error metrics; the c_Gamma=1 width insensitivity claim is not backed by a convergence plot; the ion-sphere truncation sacrifices orthogonality but Appendix B shows the effect is minor. These are fixable in revision.\n\nWho is this for: anyone doing stopping modeling in WDM/HED plasmas, and groups working on average-atom response. It deserves a serious referee—the central idea is sound and the paper is clearly argued. I'd recommend conditional acceptance; ask for the f-f consistency check and reproducibility details before publication.\n\nBest.","headline":"A solid, usable all-electron linear-response stopping model with a real but contained free-free weak spot; worth refereeing, conditional on reproducibility and f-f consistency checks.","tokens_in":18835,"tokens_out":1818,"would_cite":true,"duration_ms":21283,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A channel-mixed RPA dielectric reproduces proton stopping near and above the Bragg peak in partially ionized plasmas.","keywords":["ion stopping","dielectric response","random phase approximation","average atom","bound electrons","warm dense matter","inertial confinement fusion","linear response"],"falsifier":"Measure proton stopping in warm dense carbon near 0.5 g/cm3 and 10 eV across the Bragg peak with enough precision to distinguish cmRPA's hybridized curve from the unmixed additive model; if the data trace the additive curve, the paper's channel-mixing claim is falsified.","tokens_in":17939,"feed_emoji":"⚛️","tokens_out":11551,"duration_ms":125900,"temperature":0.7,"pith_summary":"This paper builds a single dielectric response function for ion stopping in partly ionized plasmas by putting bound-electron transitions, computed from average-atom orbitals, on the same footing as the standard Lindhard free-electron response, and mixing them at the level of the random phase approximation rather than adding their stopping contributions separately. The aim is a predictive, all-electron, linear-response stopping tool cheap enough to sweep across the densities and temperatures relevant to inertial confinement fusion. Against ambient experiments and TD-DFT simulations, the resulting channel-mixed RPA stopping power agrees well for protons at and above the Bragg peak. The paper's main new claims are that bound and free channels interfere through the dielectric denominator, that bound-bound transitions contribute at low projectile velocity when bound states are partially occupied, and that bound electrons still matter for alpha stopping in highly compressed tungsten. Applied to a recent warm dense carbon stopping experiment, the model indicates that inadequate bound-state modeling is unlikely to explain the measured deficit within linear response.","feed_headline":"Mixing bound and free electron response predicts plasma stopping","feed_subtitle":"A fast average-atom model matches proton stopping from ambient to fusion conditions.","key_machinery":"Equation (4), the channel-mixed energy loss function Im[ε⁻¹]_{cmRPA} = V_k Im[χ_bb + χ_bf + χ_ff] / |1 − V_k(χ_bb + χ_bf + χ_ff)|², is the central object. It is what distinguishes cmRPA from the Chihara-style unmixed form, Eq. (5), because taking the imaginary part of the inverse dielectric puts every channel's real and imaginary response into one denominator; that mixing is the mechanism by which bound-free transitions and the free-electron plasmon interfere. The bound channel polarizabilities come from radial matrix elements of average-atom orbitals, with a small width that smooths the continuum, and the free-free piece is the finite-temperature Lindhard function with exactly ⟨Z⟩ electrons","core_discovery":"The paper claims that the electronic stopping of an ion in a partially ionized plasma can be captured by a dielectric function in which bound-bound, bound-free, and free-free transitions all enter a single random-phase-approximation denominator, with the bound transitions built from average-atom orbitals and the free-free piece taken as finite-temperature Lindhard response. The resulting energy-loss function, Eq. (4), is the channel-mixed RPA (cmRPA) form. It is distinguished from a Chihara-style unmixed decomposition, Eq. (5), because the real and imaginary parts of all channels screen each other in one denominator, so bound and free responses interfere. Against ambient-condition stopping d","pith_inferences":["Editorial inference: because the paper's local-field-correction test raises low-velocity stopping toward experiment, adding a consistent local-field correction to the mixed denominator is a natural next step that could push cmRPA's validity down across the Bragg peak.","Editorial inference: the same dielectric function governs inelastic x-ray scattering, so the predicted bound-free/plasmon hybridization should be visible as a redistribution of spectral weight in warm dense matter x-ray Thomson scattering or electron-energy-loss measurements; a direct spectral measurement would independently test the channel-mixing claim.","Editorial inference: the f-sum violation and the Appendix A low-velocity artifact both trace to using plane-wave Lindhard for continuum states that the average-atom model sees as quasibound; replacing the Lindhard free-free response with an average-atom-consistent continuum appears to be the highest-leverage fix and would likely recalibrate the tungsten predictions at extreme compression.","Editorial inference: at projectile charges beyond protons, nonlinear effects enter precisely in the low-velocity region where bound-bound features appear, so comparing cmRPA against all-electron TD-DFT for alpha or heavier ions in partially ionized targets would map where linear response breaks and where bound contributions must be folded into binary-collision models."],"forward_implications":["Stopping curves for protons and alphas in warm dense and inertial-confinement fusion conditions can be generated orders of magnitude faster than TD-DFT, enabling parameter sweeps and diagnostic design.","At temperatures where bound states are partially occupied, bound-bound transitions give a low-velocity stopping contribution that a free-electron-only model misses.","A proper RPA treatment of bound-free transitions can either enhance or suppress stopping relative to additive channel models, depending on where the bound edge sits relative to the plasmon, so the bound correction is condition-dependent.","In double-shell inertial-confusion tungsten at 10–100 times compression, omitting bound electrons under-predicts stopping, and cold versus hot tungsten give noticeably different alpha ranges.","For the recent warm dense carbon measurement, the model implies that bound-state treatment cannot fix the measured deficit under linear response, directing attention to other physics or experimental effects."],"supporting_citations":[{"why":"supplies the ambient-condition experimental stopping curves used to validate the model across Al, Fe, Ne, and W.","marker":"[14]"},{"why":"provides TD-DFT stopping data for warm dense carbon and the all-electron comparison used to assess bound-state treatment.","marker":"[4]"},{"why":"reports the warm dense carbon stopping measurement whose theory-experiment discrepancy the model is used to interpret.","marker":"[19]"},{"why":"defines the CBC bound-state model whose contribution is replaced by the paper's 1s contribution in the carbon comparison.","marker":"[22]"},{"why":"supplies the average-atom density-of-states free-free alternative used as a corrected Lindhard limit and examined in Appendix A.","marker":"[3]"},{"why":"gives the unmixed Chihara-style decomposition, Eq. (5), against which channel mixing is defined.","marker":"[26]"},{"why":"the average-atom code that generates the bound and continuum orbitals entering the matrix elements.","marker":"[28]"},{"why":"source of the f-sum rules used to verify the dielectric and expose the free-bound inconsistency.","marker":"[35]"}],"fun_headline_variants":["Channel-mixed RPA captures ion stopping in plasmas","Bound and free electrons jointly predict plasma ion stopping","New dielectric model matches proton stopping across plasma conditions","Fast plasma stopping model validated up to fusion conditions","Channel mixing alters stopping predictions in warm dense matter"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"Everything rests on the assumption that the small inconsistency between how free electrons (plane waves) and bound electrons (average-atom orbitals) are treated does not materially change the stopping power the model predicts.","fun_headline_variants_meta":{"raw":{"variants":["Channel-mixed RPA captures ion stopping in plasmas","Bound and free electrons jointly predict plasma ion stopping","New dielectric model matches proton stopping across plasma conditions","Fast plasma stopping model validated up to fusion conditions","Channel mixing alters stopping predictions in warm dense matter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000726,"raw_usage":{"total_tokens":3102,"prompt_tokens":769,"completion_tokens":2333,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":513,"completion_tokens_details":{"reasoning_tokens":2260}},"tokens_in":513,"tokens_out":2333,"duration_ms":16073,"temperature":1.0,"reasoning_tokens":2260,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T00:11:04.905558+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure proton stopping in warm dense carbon near 0.5 g/cm3 and 10 eV across the Bragg peak with enough precision to distinguish cmRPA's hybridized curve from the unmixed additive model; if the data trace the additive curve, the paper's channel-mixing claim is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the ambient-condition experimental stopping curves used to validate the model across Al, Fe, Ne, and W."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides TD-DFT stopping data for warm dense carbon and the all-electron comparison used to assess bound-state treatment."},{"cited_title":"Malko, W","cited_arxiv_id":null,"evidence_quote":"reports the warm dense carbon stopping measurement whose theory-experiment discrepancy the model is used to interpret."},{"cited_title":"The CBC model has also been extended to use AA input, however it was found to typically overestimate stopping [4]","cited_arxiv_id":null,"evidence_quote":"defines the CBC bound-state model whose contribution is replaced by the paper's 1s contribution in the carbon comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the average-atom density-of-states free-free alternative used as a corrected Lindhard limit and examined in Appendix A."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the unmixed Chihara-style decomposition, Eq. (5), against which channel mixing is defined."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"the average-atom code that generates the bound and continuum orbitals entering the matrix elements."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"source of the f-sum rules used to verify the dielectric and expose the free-bound inconsistency."}],"review_version":1}