{"id":"35764aaa-a783-4f70-830b-d17dac6acdeb","arxiv_id":"2502.04808","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A time-resolved criterion distinguishes true electron capture from strong collisions in molecular dynamics simulations, enabling recombination broadening to be included in spectral line calculations for Z>1 ions.","lead":"A new algorithm for classical plasma simulations tracks when an electron gets captured by an ion, by checking how long the electron lingers nearby and whether its total energy is negative. This lets researchers cut the electric field history at the right moment, improving spectral line shape calculations for ions with charge greater than one.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The validation in Sec. III B checks time-averaged mean charge, not the capture-event timing that the algorithm must get right for field-sequence truncation.","rationale":"The reader's weakest-assumption identification is the heuristic tau_bound = 3 tau_T and the possibility of long-lived transient encounters being misclassified as bound. That is a real concern and is part of what I flag. My load-bearing concern goes one step further: even the validation the paper does provide does not test the temporal correctness of capture events. Fig. 5 compares a time-averaged population, not the location or duration of capture events. Since the paper's stated purpose is to select correct time-histories of the electric microfield for subsequent Stark line-shape calculations, a validation that only matches mean charge leaves the central use case under-determined. The algorithm could mis-time every capture and still satisfy Fig. 5. This is not an accusation of error; it is a gap in the evidence. I do not think the paper should be rejected: the approach is plausible, the method description is concrete, and the potential-energy-distribution comparison is a useful internal consistency check. However, the claim of 'precisely identifying' capture events is stronger than the validation supports. The recommended remedy is a sensitivity analysis of tau_bound and a timing-level comparison against an independent orbital diagnostic or an atomic-kinetics reference. The abstract/full-text discrepancy about atomic kinetics is a secondary issue, but it reinforces the need for an external validation of capture rates and times. I therefore keep the reader's CONDITIONAL verdict unchanged.","tokens_in":9908,"tokens_out":7006,"duration_ms":84555,"concrete_test":"Re-analyze the same FMD trajectories with an independent instantaneous bound-state diagnostic: for each electron-ion pair, compute the specific orbital energy and angular momentum from relative position and velocity using the regularized potential of Eq. 1, and declare a capture only when the pair has negative energy and an elliptical character for at least one full orbit. Compare the resulting capture start/end times and the distribution of valid field-sequence lengths against the paper's algorithm. If the sequence-length distributions or truncation times differ materially despite matching <Q>, the ionization-balance validation is insufficient. As a cheaper robustness check, recompute Fig. 5 and the field-sequence length statistics for tau_bound = 1, 2, 3, and 6 tau_T; if either quantity changes by more than a few percent, the 3x heuristic is load-bearing for the line-shape application.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central product of the paper is a set of correctly truncated electric-field time-histories for Stark line-shape calculations. The validation in Sec. III B, however, compares only the mean ion charge <Q> from the new criterion with <Q> obtained by integrating the ion potential-energy distribution (Fig. 5). Mean charge is a time-averaged species population; it is insensitive to when a capture is declared and to how long a bound interval lasts, provided the total bound-time fraction is roughly correct. A systematic timing error can therefore leave <Q> essentially unchanged while materially altering the field-sequence lengths and cut points on which recombination broadening depends. One concrete source of such an error is the algorithm's rule that, once an interval longer than tau_bound is found and the average pair energy is negative, the electron is 'considered to be bound for the whole interval' (Section II). This backdates the capture to the start of the interval, even if the actual binding occurs later. Likewise, a transient looping encounter with negative time-averaged pair energy and duration greater than 3 tau_T is classified as a capture and truncates a sequence that should have continued. The heuristic factor 3 in tau_bound = 3 tau_T controls exactly this timing, but no sensitivity analysis is presented. Thus the agreement in Fig. 5 does not establish that the event times and durations are 'precisely identified', which is the assumption on which the spectral-line application rests. The abstract also promises a comparison with atomic kinetic simulations, but the full text does not provide one; this omission removes a second, more physically independent check on the timing and rates of captures.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a history-based criterion for identifying when an electron becomes bound to an emitting ion in full molecular dynamics (FMD) simulations, extending an earlier Z=1 criterion to Z>=1. For each ion, the algorithm records the potential energy and the labels of the Z+1 closest electrons; when the potential energy drops below Eth, it checks whether the closest electron has remained among the Z+1 neighbors for longer than tau_bound=3*tau_T and, if the average electron-ion pair energy over that interval is negative, declares the electron bound for the entire interval. The algorithm is applied to a Z=2 helium plasma, and the resulting mean ion charge is compared with values obtained by integrating the ion potential-energy distribution. The paper argues that this validates the algorithm and enables correct truncation of electric-field time histories for Stark line-shape calculations.","tokens_in":10187,"tokens_out":11929,"duration_ms":131612,"significance":"The work addresses a genuine gap: instantaneous distance/energy criteria misclassify strong collisions and fail to count multiple captured electrons. The derivation of tau_T is clean, the algorithm is sufficiently specified for reproduction, and it is not fitted to the validation target, which is a strength. If the classification and timing are correct, the algorithm would let FMD simulations include recombination broadening for Z>1 emitters, which is currently missing. The main weakness is that the validation tests a time-averaged population statistic, not the event timing that the algorithm must deliver for its stated purpose. The paper would be considerably strengthened by direct validation of capture/release times or a clear statement that the algorithm provides only interval-level information.","major_comments":[{"comment":"The validation is an internal consistency check rather than an external test, and it does not exercise the timing information that the algorithm is designed to provide. Both the new criterion and the potential-energy-distribution method operate on the same regularized pair potential (Eq. 1) and the same simulated trajectories, so agreement between them shows that the two post-processing procedures agree with each other, not that the captured-electron intervals match a reference reality. The comparison in Fig. 5 is made on the time-averaged mean charge <Q>, which is insensitive to when within a bound interval the capture is declared and to the exact interval boundaries, provided the total bound-time fraction is approximately correct. Since the stated purpose of the algorithm is to produce correctly truncated electric-field time-histories, I request a direct test of event timing, for example by running synthetic trajectories with known injection and release times, or by comparing spectra computed with and without the truncation. In addition, the abstract's claim that results are compared with atomic kinetic simulations is not supported by Section III B, which contains no such comparison.","section":"Sec. III B, Fig. 5"},{"comment":"The threshold tau_bound=3*tau_T is a hand-set heuristic, and the paper gives no evidence that it separates genuine captures from long transient encounters. The derivation of tau_T assumes a circular orbit at the thermal velocity in the hyperbolic potential, but the simulation uses the regularized potential of Eq. (1); the paper only checks the consistency condition r_T>a (Eq. 9). Most importantly, no sensitivity scan is reported, so the reader cannot tell whether the ionization balance and, more crucially, the number and timing of bound intervals depend strongly on the factor 3. I recommend reporting results for at least a few values of tau_bound/tau_T (e.g., 1, 2, 4, 6) and discussing the behavior at intermediate temperatures, where the statement that recombination is negligible at T about 3E0 does not apply.","section":"Sec. II, determination of tau_bound"},{"comment":"The algorithm backdates captures to the beginning of the interval in which the electron remained among the Z+1 closest neighbors. The text states that if the interval is longer than tau_bound and the average pair energy is negative, the electron is 'considered to be bound for the whole interval.' This is a causal assumption, not a demonstrated property: the actual binding may start at any point within that interval, and in the Fig. 1 example electron 298 is labeled bound from t about 0 even though the potential-energy threshold is first crossed at t about 2.7 t0. For the intended spectral-line application, cutting the field sequence at the interval start removes emission that may have occurred while the electron was still free. The manuscript should define the capture onset time explicitly (e.g., the first time at which the energy and neighbor conditions are jointly satisfied) and quantify how much the interval-start truncation changes the generated field sequences.","section":"Sec. II, bound-interval declaration"},{"comment":"The potential-energy-distribution method used for validation assumes that the ion potential-energy lobes are sufficiently separated to be identified with charge states. The examples in Figs. 2 and 3 use a large ionization potential (Vi=30E0) and one value of the coupling; under those conditions the lobes are reasonably distinct. The paper does not show that the method, or the agreement in Fig. 5, survives for lower Vi, stronger coupling, or Z>2, where the lobes overlap and the integral boundaries become ambiguous. The claim that the algorithm is validated for emitters with Z>=1 therefore needs either additional scan cases or a precise statement of the conditions under which the validation method is reliable.","section":"Sec. III B"}],"minor_comments":[{"comment":"The sign convention for V_i is inconsistent: the text speaks of negative binding energies, while Eq. (1) introduces V_i as a positive quantity. Please define the convention explicitly.","section":"Sec. II, Eq. (1)"},{"comment":"The text says the criterion is valid if the ionization potential is larger than twice the thermal energy, but the derived condition is V_i > 3 k_B T_e; please correct the wording.","section":"Sec. II, after Eq. (9)"},{"comment":"The phrase 'the potential energy jumps below 6 Eth' is ambiguous and appears inconsistent with the plotted range (Eth=-10E0, so 6Eth would be -60E0). Please rephrase, e.g., 'drops to about 6E0 below the threshold.'","section":"Sec. III A"},{"comment":"The 'average total energy of the electron-ion pair' is not defined. State explicitly which kinetic terms and which potential terms enter the average and in which reference frame.","section":"Sec. II, trapped-electron criterion"},{"comment":"Each point is said to correspond to a group of simulations, but no error bars or run-to-run scatter are shown. Please quantify the statistical uncertainty.","section":"Figs. 4 and 5"},{"comment":"The algorithm is described in prose; a pseudocode listing or flowchart would make the interval-merging and backdating behavior unambiguous.","section":"Sec. II"}],"recommendation":"major_revision","confidential_remarks":"The manuscript fits the journal's scope and the method is a useful contribution, but the central claim of 'precisely identifies' outruns the evidence because the validation is an internal time-averaged check. The authors should be asked for either a timing-level validation or a more limited statement of what the algorithm guarantees. No issues with novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this paper does something genuinely useful and mostly checks out, but its validation is weaker than its language suggests. The new piece is a time-history criterion for deciding when an electron is captured by an ion in full molecular dynamics simulations for Z>1 emitters: an electron is counted as bound if it stays among the Z+1 closest neighbors for longer than 3 tau_T and the pair's average total energy over that interval is negative. That's a real extension of the earlier Z=1 distance/energy rule and it targets a known gap in Stark-broadening simulations, namely recombination broadening from electron capture on multi-charged emitters.\n\nThe algorithm is clearly specified, the derivation of tau_T is clean classical mechanics, and the example in Fig. 1 is helpful. I also appreciate the explicit limitations paragraph. The check against the potential-energy distribution method does show good agreement for the mean ion charge over a wide temperature range, which is a reasonable internal consistency test.\n\nThe soft spots are real, though not fatal. The validation compares only time-averaged mean charge, which is insensitive to when captures are declared and how long bound intervals last. The algorithm backdates the capture to the start of the interval once it finds a qualifying interval, and the factor 3 in tau_bound = 3 tau_T is a hand-set heuristic with no sensitivity analysis. A systematic timing error could leave mean charge almost unchanged while corrupting the very field-sequence truncation points the algorithm is supposed to produce. The abstract promises a comparison with atomic kinetic simulations, but the full text does not include one; that's a mismatch that should be fixed. There are also no error bars, and no code or data release, which makes independent checking harder than it should be.\n\nOn balance, the central argument holds up as a plausible and useful algorithm; it just isn't demonstrated as precisely as claimed. The paper deserves a serious referee, but I'd send it back for revisions: add a sensitivity study on tau_bound and Eth, report uncertainties, and either add the atomic-kinetics comparison or drop it from the abstract.","headline":"A useful time-history algorithm for detecting electron capture in FMD simulations that merits revision, not rejection, but its validation is too coarse to support the 'precisely identifies' claim.","tokens_in":10740,"tokens_out":1786,"would_cite":false,"duration_ms":17582,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A new time-history criterion identifies when a plasma electron is captured by an ion in full molecular dynamics simulations, separating true recombination from strong collisions.","keywords":["electron capture","full molecular dynamics","Stark broadening","recombination broadening","ionization balance","plasma spectroscopy","electric microfield","bound-free classification"],"falsifier":"Run the same simulated plasma conditions with $\\tau_{\\mathrm{bound}}$ set to $2\\tau_T$, $3\\tau_T$, and $5\\tau_T$; if the inferred ionization balance shifts beyond statistical scatter, the heuristic threshold is doing the labeling work rather than physical capture. A more direct check is to follow a pair that the algorithm labels bound and verify that its relative trajectory closes on itself as a stable orbit for many periods instead of slowly separating.","tokens_in":9732,"feed_emoji":"⚛️","tokens_out":7927,"duration_ms":80167,"temperature":0.7,"pith_summary":"This paper introduces an algorithm for full molecular dynamics (FMD) simulations of plasmas that decides when an electron has actually been captured by an ion of charge $Z\\ge 1$. The criterion uses the history of each electron-ion pair rather than a single snapshot of distance or energy, so it can distinguish a briefly close strong collision from an electron that settles into a stable orbit. That distinction is what lets simulations cut an emitter's electric-field time-history at the moment of true recombination, adding a broadening contribution that line-shape calculations have so far omitted. The authors validate the algorithm by comparing the ionization balance it produces with one obtained from ion potential-energy distributions and find agreement over a range of temperatures. If correct, the method makes recombination broadening accessible in Stark line shape calculations.","feed_headline":"Time-history test separates electron capture from strong collisions","feed_subtitle":"Lets line-shape simulations add recombination broadening instead of discarding every close encounter.","key_machinery":"The load-bearing object is a retrospective time-history classifier built on three stored quantities per ion: its potential energy, the electric field at its position, and the labels of the $Z+1$ closest electrons at every timestep. When the potential energy crosses the threshold $E_{\\mathrm{th}}=-2V_i/3$, the algorithm looks back over the residence time of the first neighbor among that window; if the residence exceeds $\\tau_{\\mathrm{bound}}=3\\tau_T$ and the pair's mean total energy is negative, the electron is labeled bound for the whole interval. The $Z+1$ window matters because a fast free electron can briefly pass closer to the core than an already bound electron without unbinding it. This converts a point-in-time bound/free decision into a history-based decision, which is what lets the simulation tell recombination from strong collisions and choose where to cut the electric-field sequence.","core_discovery":"Within a classical molecular dynamics simulation, an electron is declared bound to an emitter of charge $Z$ when, at the moment the emitter's potential energy dips below $E_{\\mathrm{th}}=-2V_i/3$, the closest electron has already remained among the $Z+1$ closest electrons for longer than $\\tau_{\\mathrm{bound}}=3\\tau_T$, and the time-averaged total energy of that electron-ion pair over the interval is negative. Here $\\tau_T$ is the period of a circular orbit at the distance where the Coulomb attraction balances the thermal kinetic energy of an electron. A strong collision fails the test because the interloper does not stay among the nearest neighbors long enough; a genuinely captured electron passes it. Once an electron is declared bound, the field sequence for the original ion is cut, the ion is reclassified to charge $Z-1$ with its emission coherently lost, and the sequence resumes only when the ion is reionized.","pith_inferences":["Beyond the paper, the same neighbor-residence logic could be adapted to other snapshot-ambiguous events in classical plasma simulations, such as Debye-sphere cluster membership or transient molecule formation.","Beyond the paper, a direct scan of the factor 3 in $\\tau_{\\mathrm{bound}}$ would test whether the ionization-balance curve is robust or partly controlled by that heuristic choice.","Beyond the paper, the capture time itself is a spectral parameter: the abrupt loss of coherence at the cut point should produce a measurable recombination-broadening component in strongly coupled plasmas that could be compared with experiment."],"forward_implications":["Valid electric-field time-histories from FMD simulations can now be produced without discarding either recombination or strong collisions, so Stark line shape calculations can include both effects.","The algorithm doubles as an ionization-balance diagnostic, since the number of electrons declared bound to each emitter gives the fractional populations of charge states.","The earlier $Z=1$ treatment is extended to $Z\\ge 1$, where simple distance-based tests cannot even count how many electrons are trapped.","The agreement between the algorithm's ionization balance and potential-energy-distribution counting is evidence that the captured-electron labels are physically meaningful rather than arbitrary.","Recombination broadening becomes a calculable contribution in strongly coupled plasmas, where it was previously omitted."],"supporting_citations":[{"why":"Supplies the regularized Coulomb potential, the simulation energetics, and the earlier Z=1 bound/free criterion that this algorithm extends to Z>1.","marker":"[15]"},{"why":"Identifies recombination broadening in emitter field sequences as the effect the new algorithm is meant to handle.","marker":"[13]"},{"why":"Defines the distance-based bound/free criterion whose failures motivate the history-based method.","marker":"[16]"},{"why":"Documents that for Z>1 the simple criterion cannot count how many electrons are trapped, the gap this algorithm fills.","marker":"[17]"},{"why":"Catalogues regularization procedures for the Coulomb potential, against which the chosen parabolic-hyperbolic regularization is situated.","marker":"[18]"},{"why":"Uses FMD simulations for line broadening from charged emitters and motivates why full particle interaction matters for such radiators.","marker":"[12]"}],"fun_headline_variants":["New algorithm flags true electron capture in plasma MD","Bound vs flyby: MD algorithm separates capture from collisions","Time-and-energy test spots captured electrons for line shapes","Capture criterion keeps recombination in spectral line models"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The algorithm assumes that an electron which remains among the $Z+1$ closest neighbors for more than $3\\tau_T$ with negative mean total energy is genuinely bound; a long-lived three-body encounter or resonant looping trajectory could satisfy the same conditions without being a stable capture.","fun_headline_variants_meta":{"raw":{"variants":["New algorithm flags true electron capture in plasma MD","Bound vs flyby: MD algorithm separates capture from collisions","Time-and-energy test spots captured electrons for line shapes","Capture criterion keeps recombination in spectral line models"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000219,"raw_usage":{"total_tokens":1447,"prompt_tokens":953,"completion_tokens":494,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":433}},"tokens_in":569,"tokens_out":494,"duration_ms":6330,"temperature":1.0,"reasoning_tokens":433,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T21:22:35.778598+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same simulated plasma conditions with $\\tau_{\\mathrm{bound}}$ set to $2\\tau_T$, $3\\tau_T$, and $5\\tau_T$; if the inferred ionization balance shifts beyond statistical scatter, the heuristic threshold is doing the labeling work rather than physical capture. A more direct check is to follow a pair that the algorithm labels bound and verify that its relative trajectory closes on itself as a stable orbit for many periods instead of slowly separating.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the regularized Coulomb potential, the simulation energetics, and the earlier Z=1 bound/free criterion that this algorithm extends to Z>1."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies recombination broadening in emitter field sequences as the effect the new algorithm is meant to handle."},{"cited_title":"Gigosos, D","cited_arxiv_id":null,"evidence_quote":"Defines the distance-based bound/free criterion whose failures motivate the history-based method."},{"cited_title":"Lara, Calculations of Stark spectra of strongly coupled plasmas by molecular dynamics simulation , Ph.D","cited_arxiv_id":null,"evidence_quote":"Documents that for Z>1 the simple criterion cannot count how many electrons are trapped, the gap this algorithm fills."},{"cited_title":"Gonz´ alez-Herrero,Study of the Stark broadening of the 3s-3p spectral lines in Be-like ions –Molecular Dynam- ics simulations– , Ph.D","cited_arxiv_id":null,"evidence_quote":"Catalogues regularization procedures for the Coulomb potential, against which the chosen parabolic-hyperbolic regularization is situated."},{"cited_title":"Stambulchik, D","cited_arxiv_id":null,"evidence_quote":"Uses FMD simulations for line broadening from charged emitters and motivates why full particle interaction matters for such radiators."}],"review_version":1}