{"id":"bb6b4098-0a93-4b79-bfc3-0eaba512eb18","arxiv_id":"2607.11253","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"In two-temperature MD of Si cascades, the ion–electron coupling function controls defect yield, clustering, and recombination, with effects that reverse between Stillinger–Weber and Tersoff/ZBL potentials.","lead":"Simulations show that the mathematical form of atom–electron coupling strongly changes how many defects form, cluster, and recombine in irradiated silicon. That choice, together with the atomic potential, can flip recombination trends and must be fixed carefully for reliable radiation-damage models.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"Fixed low-T Ce and κe mute bidirectional feedback by construction, so the claim that coupling form (not electronic feedback) dominates rests on an untested electronic-transport regime.","rationale":"The Reader correctly isolates the fixed low-T electronic coefficients as the weakest modeling assumption and assigns CONDITIONAL with high confidence. My stress-test confirms that this assumption is load-bearing for the central claim: the paper’s own argument that “coupling form matters more than feedback” is licensed only by the muted fixed/variable-Te comparison performed under those coefficients. No internal contradiction or data-fabrication issue appears; the comparative trends across energies and the SRIM/TDDFT consistency checks remain supportive. The concrete test (temperature-dependent Ce) is the minimal, decisive check that would either validate or qualify the claim. Because the Reader already flagged the same soft spot and recommended CONDITIONAL pending further scrutiny, no change of verdict is required.","tokens_in":12167,"tokens_out":716,"duration_ms":7019,"concrete_test":"Re-run the 20 keV ensemble (both SW and T3/ZBL, both coupling functions) with a temperature-dependent electronic heat capacity that rises toward the free-electron value once Te exceeds ~1000 K (or with Ce increased by 1–2 orders of magnitude while keeping κe fixed). If the fixed-Te versus variable-Te difference in surviving FPs or recombination efficiency R becomes comparable to the quadratic-versus-four-density difference, the claim that coupling form dominates over electronic feedback is weakened.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The strongest claim is that the functional form of the ion–electron coupling (quadratic vs four-density) critically controls defect production, clustering and recombination, and that its effect on cooling-phase recombination further depends on the interatomic potential (SW vs T3/ZBL). The paper supports this by showing that fixed-Te and variable-Te runs give nearly identical FP counts (Fig. 2, §3.1) and therefore concludes that “the functional form of the coupling and the resulting energy losses … have a stronger influence … than the subsequent feedback from a heated electronic system.” That comparison, however, is performed with Ce and κe fixed at the low-temperature values 5\times10⁻⁶ eV K⁻¹ Å⁻³ and 5\times10⁻³ eV K⁻¹ Å⁻¹ ps⁻¹ for every voxel throughout the cascade (Simulation details). Under these parameters the electronic subsystem cannot store or transport heat efficiently once Te rises, so bidirectional energy return is suppressed by construction. Consequently the fixed-vs-variable-Te null result does not demonstrate that feedback is physically unimportant; it only shows that feedback is unimportant when the electronic transport coefficients are held at their cold values. The ranking of the two coupling functions and the reported potential-dependent recombination trends therefore rest on an untested electronic-transport regime.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript studies primary radiation damage in diamond-structure silicon with two-temperature MD (LAMMPS + USER-EPH), comparing a quadratic local-density coupling function and a four-density coupling function fitted to rt-TDDFT trajectories. Cascades are run for PKA energies 0.1–20 keV with Stillinger–Weber and T3/ZBL potentials, both fixed and variable electronic temperature, and against friction-based electronic-stopping (ESP) baselines with several kinetic-energy cut-offs. Residual defects are quantified by Wigner–Seitz and coordination analysis; clustering and recombination efficiency R = 1 − Nsurv/Nmax are reported. The central claim is that the functional form of the ion–electron coupling controls defect production, clustering and recombination more strongly than electronic-temperature feedback under the chosen parameters, and that the effect of coupling on cooling-phase recombination further depends on the interatomic potential (compact amorphous-like SW cascades versus more extended T3/ZBL clusters).","tokens_in":12530,"tokens_out":1200,"duration_ms":11454,"significance":"If the trends hold under more realistic electronic transport, the work supplies a concrete, transferable demonstration that threshold-free density-dependent couplings are not interchangeable for cascade predictions in Si, and that electronic and atomic models must be co-validated. Strengths include transparent methods (cell sizes, adaptive timesteps, explicit cut-offs, PCA cluster lengths), systematic cross-checks against SRIM/ESP friction baselines and against TDDFT-motivated energy-loss curves, and the explicit potential-dependent recombination analysis (Figs. 7–11). The four-density coupling’s closer agreement with SRIM in the ballistic regime is a useful external benchmark. The results are of direct interest for radiation-damage modelling of Si detectors and for the broader community using UTTM-style electronic stopping.","major_comments":[{"comment":"Simulation details and §3.1 (Figs. 2–3, 6): Ce and κe are fixed at the low-temperature values 5×10⁻⁶ eV K⁻¹ Å⁻³ and 5×10⁻³ eV K⁻¹ Å⁻¹ ps⁻¹ for every voxel throughout the cascade. Under these coefficients the electronic subsystem cannot store or transport heat efficiently once Te rises, so bidirectional energy return is suppressed by construction. The near-identity of fixed-Te and variable-Te FP counts therefore does not demonstrate that feedback is physically unimportant; it only shows that feedback is unimportant when transport coefficients remain cold. The ranking of the two coupling functions and the potential-dependent recombination trends rest on this untested electronic-transport regime. At least a limited sensitivity study (temperature-dependent Ce/κe, or literature high-Te values) is needed before the claim that “coupling form … [has] a stronger influence … than the subsequent fe","section":null},{"comment":"§3.1–3.2 and Conclusions: the four-density coupling is presented as more physically consistent because its ballistic energy losses track SRIM and because it was fitted including close-collision trajectories. SRIM, however, is an empirical average over straight-path, intact-lattice conditions and is not a first-principles ground truth for disordered cascade cores. The paper should state more carefully that SRIM agreement is a useful consistency check rather than decisive validation, and should quantify how much of the defect-number difference survives when both couplings are constrained to the same integrated Se (or when only the high-density core region of the four-density function is varied).","section":null}],"minor_comments":[{"comment":"Section numbering is inconsistent: the text states “Section 2 describes … Section 2 presents and discusses the results. Section 4 summarizes…” (Introduction, final paragraph). Renumber Results as §3 and Conclusions as §4 consistently.","section":null},{"comment":"Fig. 1 caption and body: “Xin ESP-Xstands” and similar missing spaces/typos appear in several figure captions and in the abstract (“threshold-free approach”). A careful proof-read is needed.","section":null},{"comment":"Table 1 and Simulation details: electronic voxel side length is given as “≈25 Å”; the precise value used and the resulting number of electronic voxels relative to the atomic cell should be stated for reproducibility.","section":null},{"comment":"Fig. 4 and Fig. 5: the gray ESP bands are helpful, but the number of independent cascades averaged for each PKA energy (especially 100 keV) is not stated; adding N and error bars (or standard error) would strengthen the comparison.","section":null},{"comment":"Recombination efficiency R is defined with coordination defects (Eq. 1) while surviving damage is also reported as Frenkel pairs (WS). A short note on how the two defect definitions correlate would avoid ambiguity when comparing R across potentials.","section":null},{"comment":"References [38] and [39] are listed as “Research data” and “Submitted”; if still unpublished at acceptance, a brief statement of availability (or deposition) would help readers reproduce the four-density function.","section":null}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid, well-executed parametric study that belongs in a materials-modelling or radiation-effects journal. The Ce/κe limitation is the only load-bearing methodological gap; once addressed (even with a short sensitivity appendix) the paper should be publishable. No concerns about novelty disclosure or citation pattern."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The paper's real payload is a clean head-to-head of two UTTM couplings (Jarrin's quadratic vs the authors' four-density form) across 0.1–20 keV with both SW and T3/ZBL. Surviving FPs, peak coordination defects, cluster size/length, recombination efficiency R = 1 − Nsurv/Nmax, and opposite- vs same-type separation PDFs all move with the coupling, and the recombination response flips sign between the two potentials. That joint dependence is the new, usable result.\n\nThey do the work carefully. Cell sizes, adaptive timesteps, WS and coordination cutoffs, PCA lengths, and the ESP/SRIM friction baselines are all stated. Ballistic energy loss for the four-density form tracks SRIM (and the TDDFT-motivated trajectories it was fitted to) better than the quadratic form; that external check is real. The cascade-morphology discussion (compact amorphous-like SW pockets vs more elongated T3/ZBL clusters with preferred separations) explains why recombination efficiency moves in opposite directions without hand-waving.\n\nThe soft spot the stress-test flags is real but limited. Ce and κe are locked at the low-T values from Jarrin et al. for every voxel, so the electronic system cannot store or carry heat once Te rises. The near-null fixed-vs-variable-Te comparison therefore only shows that feedback is unimportant under cold transport coefficients; it does not prove feedback is physically secondary. The ranking of the two couplings and the potential-dependent recombination trends still stand on the energy-loss differences during the ballistic phase, which are less sensitive to that choice. Free parameters (cutoffs, voxel size, coupling coefficients) are inherited from prior work rather than re-tuned to the defect numbers, so circularity is mild. No code or ensembles are shipped, and error bars on the averages are thin.\n\nThis is for people who run or consume primary-damage statistics in Si (or who care about how electronic stopping is parameterized in MD). It is not a methods breakthrough, but it is a clear warning that the coupling functional form is not a secondary detail. I would send it to referees; the central claim is supported and the modeling limitation is easy to state and discuss. Worth citing when you next need numbers or a caution on electronic-dissipation choices.","headline":"Solid comparative cascade study: coupling form and potential jointly set Si primary-damage metrics; the fixed low-T Ce/κe choice mutes feedback by construction but does not sink the main result.","tokens_in":13151,"tokens_out":599,"would_cite":true,"duration_ms":6670,"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":"The functional form of ion–electron coupling decides how many defects, and how they cluster, in silicon radiation cascades—and that effect itself depends on which atomic potential is used.","keywords":["radiation damage","silicon","two-temperature model","electronic stopping","electron-phonon coupling","molecular dynamics","defect clustering","cascade recombination"],"falsifier":"Repeat the same 20 keV cascade suite with electronic heat capacity and conductivity allowed to rise with electronic temperature; if the ranking of defect yields and recombination efficiencies between the two coupling functions then reverses or collapses, the present conclusions are parameter-dependent rather than general.","tokens_in":13054,"feed_emoji":"⚛️","tokens_out":945,"duration_ms":12087,"temperature":0.7,"pith_summary":"Radiation damage in silicon begins when a knock-on atom rips through the lattice, but part of its energy is also lost to electrons. Classical molecular dynamics usually ignores that channel or approximates it with a crude cut-off. This paper shows that a threshold-free, density-based two-temperature model still produces very different numbers of surviving defects, different cluster sizes, and different recombination efficiencies once the precise functional form of the ion–electron coupling is changed. Two couplings fitted to the same first-principles stopping data—one quadratic, one a four-density form that better captures close collisions—yield systematically different damage when the same cascade is run with either the Stillinger–Weber or the Tersoff/ZBL potential. The four-density form dissipates more energy during the ballistic phase, suppresses peak disorder, and reduces clustering; yet whether it also increases or decreases final recombination depends on which interatomic potential is chosen. The practical message is that both the electronic coupling and the atom–atom potential must be physically grounded if primary-damage statistics are to be trusted for device-lifetime modelling.","feed_headline":"Coupling form, not just cut-off, sets silicon defect yields","feed_subtitle":"Two density-based ion–electron models give different Frenkel pairs and recombination—and the ranking flips with the atomic potential","key_machinery":"The threshold-free unified two-temperature model (UTTM) that couples atoms to a local electronic density via either a quadratic or a four-density coupling function, both fitted to real-time TDDFT stopping powers.","core_discovery":"Within a unified two-temperature molecular-dynamics framework, the functional form of the local-density ion–electron coupling controls defect production efficiency, clustering, and recombination in silicon cascades; the same coupling can raise or lower recombination efficiency depending on whether the Stillinger–Weber or Tersoff/ZBL potential is used, because the two potentials generate different cascade morphologies and different spatial correlations between vacancies and interstitials.","pith_inferences":["If the same density-based couplings were applied to germanium or silicon carbide, the potential-dependent recombination crossover observed here would likely reappear wherever cascade compactness differs strongly between empirical potentials.","Device-level rate-theory models that currently ingest MD defect yields as fixed input may need an additional “coupling-form” uncertainty band once electronic effects are treated at this fidelity.","A natural next experiment is to re-run the identical PKA set with a machine-learned potential known to reproduce DFT defect formation energies; any residual coupling-form effect would then isolate the electronic contribution more cleanly."],"forward_implications":["Primary-damage databases used for silicon detector lifetime models must quote both the coupling function and the interatomic potential, not merely the PKA energy.","Friction-only electronic-stopping models with arbitrary kinetic-energy cut-offs will continue to scatter widely in predicted Frenkel-pair numbers.","The four-density coupling’s closer match to SRIM in the ballistic phase suggests it should be preferred for high-energy cascade statistics until a better TDDFT-fitted form appears.","Cascade morphology (compact amorphous pockets versus elongated crystalline defect strings) is co-determined by electronic dissipation and the short-range repulsion of the chosen potential."],"fun_headline_variants":["Ion–electron coupling form sets silicon defect yields and recombination","Local-density coupling shape controls Frenkel pairs by potential","Coupling function, not cutoff, governs cascade defects in Si","Same coupling raises or lowers recombination by SW vs Tersoff","Functional form of density-based coupling drives Si damage outcomes"],"cache_read_input_tokens":128,"weakest_assumption_plain":"Electronic heat capacity and conductivity are held fixed at their low-temperature values even in voxels that become strongly heated, so energy feedback from electrons to the lattice is largely muted by construction.","fun_headline_variants_meta":{"raw":{"variants":["Ion–electron coupling form sets silicon defect yields and recombination","Local-density coupling shape controls Frenkel pairs by potential","Coupling function, not cutoff, governs cascade defects in Si","Same coupling raises or lowers recombination by SW vs Tersoff","Functional form of density-based coupling drives Si damage outcomes"]},"model":"grok-4.5","effort":"low","cost_usd":0.003054,"raw_usage":{"total_tokens":1055,"prompt_tokens":723,"num_sources_used":0,"completion_tokens":85,"cost_in_usd_ticks":30540000,"prompt_tokens_details":{"text_tokens":723,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":247,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":723,"tokens_out":85,"duration_ms":3239,"temperature":1.0,"reasoning_tokens":247,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T05:58:50.484165+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Repeat the same 20 keV cascade suite with electronic heat capacity and conductivity allowed to rise with electronic temperature; if the ranking of defect yields and recombination efficiencies between the two coupling functions then reverses or collapses, the present conclusions are parameter-dependent rather than general.","supporting_citations":[],"review_version":1}