{"id":"3ea57343-fa9a-4afd-a822-bf7c094ece85","arxiv_id":"2509.02026","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Backpropagated classical trajectories separate tunneling from over-barrier ionization in a model helium atom, showing that the over-barrier path sets in at smaller initial transverse electron momentum and at a Stark-shift-corrected intensity threshold.","lead":"This paper analyzes how electrons escape atoms in very strong laser fields, separating the classic tunnel-out path from the over-the-barrier path using a backpropagation technique. It finds that which path an electron takes depends not only on laser strength but also on the electron's sideways motion, a distinction experiments may be able to read directly from momentum images.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The OBI tag (Fig. 2b: 'orbits the ion' + local speed minimum) is never validated against the paper's own energy-over-saddle criterion, and since orbiting kinematically requires small k_perp, the central k_perp-dependence finding may be pre-encoded by the tagging rule.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern I found: the unvalidated, qualitative orbit criterion behind the OBI tag. I considered and rejected two alternative candidates. (a) The mechanical fit problem (P_OBI = −0.288 I^−1.758 + 0.809 is negative over part of the plotted range; I_OBI = 0.556 a.u. is a fit zero compared with I_Stark_th = 0.537 a.u., a 3.5% intensity / 1.7% field difference). This is real but secondary: it attacks only the threshold sub-claim; if the tag were validated, a physical fit form would likely preserve the qualitative Stark-shift conclusion. (b) The Fig. 3(c) sum check: as the tag is exhaustive, the sum reproducing the full PMD is essentially tautological, so this cannot serve as the validation the paper implies; this strengthens, rather than replaces, the tag concern. My concern sharpens the reader's: the topology tag is never compared with the energy-over-saddle condition that the paper itself states as the definition of OBI, and the 'orbiting' requirement kinematically favors small k_perp, so the central k_perp-dependence finding is at risk of being pre-encoded. The proposed test uses only quantities already in the paper (the effective-barrier energy expression and the saddle picture of Fig. 1b), so it is directly executable by the authors on existing data. Credit where due: the TDSE backpropagation machinery is established for TI (Ni et al., PRL 117, 023002), the two-cycle circular pulse cleanly suppresses intercycle interference, the depletion correction (Eq. 2) is the standard form, and Eq. (5) is a self-consistent Stark-shift estimate requiring only α; if the confusion-matrix test passes, the CONDITIONAL verdict should be upgraded to ACCEPT. If it fails, the central deliverable loses support. Either way, my recommendation does not change the reader's verdict, so verdict_should_be is UNCHANGED.","tokens_in":9583,"tokens_out":17622,"duration_ms":167446,"concrete_test":"For every backpropagated trajectory, record the tag point (k∥=0 for TI; local speed minimum before orbiting for OBI), the instantaneous field F(t*), position r*, and transverse momentum k_perp*. Evaluate E* = ½k_perp*² − 1/r* + r*·F(t*) and compare against the saddle energy E_saddle(k_perp*, F(t*)) of the effective potential in Fig. 1(b): tag OBI if E* exceeds E_saddle, else TI. Build the 2×2 confusion matrix against the topology tag, and compare the PMD boundary radius with the solution of E_ground = E_saddle(k_perp, F) at the emission times. If agreement is ≥98% and the boundary matches within envelope uncertainties (~10%), the topology tag is validated; if agreement is substantially worse, the orbit criterion is not a faithful OBI diagnostic and the k_perp-dependence, probability competition, and I_OBI threshold claims are unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the OBI tagging rule in Fig. 2(b): a backpropagated trajectory is OBI if it 'returns to the proximity of the parent ion and orbits around it,' with the barrier top at a local speed minimum (Fig. 2b2). The rule is described qualitatively—no tolerances, no convergence tests—and it is never checked against the paper's own independent OBI condition, the effective-barrier energy criterion stated in the text: 'OBI occurs for specific values of transverse momentum k⊥ when the ground-state energy exceeds the corresponding saddle surface level' (E = ½k_perp² − 1/r + r·F_c, Fig. 1b). Three consequences. (1) Circularity: at the tag point, 'orbiting the ion' selects trajectories with small angular momentum about the ion, which is kinematically near-synonymous with small k_perp; the paper's finding that OBI electrons occupy the central (inner-arc) PMD region (Fig. 3b) is therefore partially pre-encoded, and the 'sharp border' is guaranteed by the binary rule rather than demonstrated. (2) The Fig. 3(c) sum check merely shows the tag is exhaustive and backpropagation preserves final momenta; it does not validate the mechanism attribution. (3) All downstream claims—the k_perp-dependent TI/OBI boundary, the competitive probability curves (Fig. 5, Eqs. 3–4), and the Stark-shift threshold comparison (I_OBI = 0.556 a.u. vs I_Stark_th = 0.537 a.u., Eq. 5)—inherit any error in the tag. Secondary but real: the fit P_OBI = −0.288 I^−1.758 + 0.809 is negative over the lower part of the plotted intensity range, so I_OBI is the zero of an unphysical fit; that is the number compared with the analytic threshold. If the orbit criterion accidentally tags low-angular-momentum tunneling trajectories as OBI, every headline finding loses support.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript extends the classical backpropagation method to over-barrier ionization (OBI) in a two-dimensional model helium atom driven by a short, circularly polarized laser pulse. Starting from a TDSE solution, the ionized wave packet is converted into classical trajectories that are propagated backward in time. Trajectories are classified as tunneling ionization (TI) if they rebound at a tunnel exit identified by k_parallel = 0, or as OBI if they return to the vicinity of the parent ion, orbit it, and have the barrier top at a local speed minimum. The authors use this classification to separate the photoelectron momentum distribution (PMD) and ionization-time distribution into TI and OBI contributions, to argue that the TI/OBI boundary depends on the initial transverse momentum k_perp as well as field strength, to extract intensity-dependent probabilities with empirical fits P_OBI = -0.288 I^(-1.758) + 0.809 and P_TI+OBI = -0.313 I^(-1.377) + 1.355, and to compare the numerical threshold I_OBI = 0.556 a.u. with a Stark-shift-corrected saddle threshold I_Stark = 0.537 a.u. The central claim is that this provides the first trajectory-level dynamical picture of OBI in a full-Coulomb Hamiltonian.","tokens_in":9954,"tokens_out":5844,"duration_ms":68997,"significance":"If the classification rule is physically faithful, the work would fill a genuine gap: most strong-field theories either assume a zero-range potential and therefore exclude OBI, or treat OBI only through rates. The manuscript offers a concrete, falsifiable pipeline: full TDSE with the Coulomb potential, classical backpropagation, an exact-by-construction decomposition of the PMD, and an explicit Stark-shifted threshold that is compared with a numerical value. These are strengths. The paper also makes a specific and testable prediction that the TI/OBI border depends on k_perp, and it supplies explicit fitting functions. However, the central classification rule is described qualitatively and is not validated against an independent OBI diagnostic. Because nearly all quantitative conclusions—the PMD separation, the k_perp dependence, the competitive probability curves, and the threshold comparison—inherit the tagging rule, the significance of the paper is currently conditional on this validation.","major_comments":[{"comment":"The OBI tagging rule is not specified as an algorithm. 'Returns to the proximity of the parent ion and orbits around it' requires a distance threshold and a definition of 'orbits' (minimum winding angle, number of turns, allowed time window), and 'local speed minimum' requires a tolerance or a bracketing criterion. No convergence tests are reported. Since every PMD pixel, ionization-time bin, and probability in Fig. 5 is built from this tag, the load-bearing step must be reproducible and its sensitivity quantified. Please provide an explicit, implementable rule and tests showing the results do not depend on the chosen tolerances.","section":"Backpropagation OBI identification, Fig. 2(b)"},{"comment":"The sum check in Fig. 3(c) does not validate the mechanism attribution. By construction, each backpropagated trajectory is assigned to exactly one of two classes, so the sum of the two distributions equals the total distribution up to numerical discretization. It confirms exhaustiveness and that backpropagation preserves final momenta, but it says nothing about whether the 'OBI' class actually corresponds to over-barrier ionization. This statement should be softened, and an independent check should be supplied.","section":"Fig. 3(c), 'This confirms the accuracy...'"},{"comment":"The manuscript states that 'OBI occurs for specific values of transverse momentum k_perp when the ground-state energy exceeds the corresponding saddle surface level,' and uses this to explain the central OBI arc. This is an independent, quantitative criterion, but it is never compared with the numerical TI/OBI boundary extracted from the trajectory tag. A direct test would be to plot the observed OBI region in the (field, k_perp) plane against the saddle energy condition E = 1/2 k_perp^2 - 1/r + r F_c. Without this comparison, the k_perp-dependence claim rests entirely on the qualitative tagging rule, and the 'sharp border' statement in the following paragraph is circular: a binary classification always produces a sharp boundary, regardless of whether the underlying physics is sharp.","section":"Effective potential barrier discussion, Fig. 1(b) and Eq. for E"},{"comment":"The sentence 'This is because the categorization of TI and OBI is based on trajectory topology, which is a binary condition' explains a property of the classification, not a property of the ionization dynamics. The physical question is whether the binary mechanism boundary, after appropriate binning and finite statistics, appears as a sharp contrast in physical observables. The current wording may mislead readers into thinking that sharpness of the extracted PMD boundary is a discovery rather than an artifact of the tagging procedure. Please replace this explanation with a quantitative analysis of the boundary width and its relation to the saddle criterion.","section":"'Remarkably, this border is sharp'"},{"comment":"The fit P_OBI = -0.288 I^(-1.758) + 0.809 is negative for intensities below the fitted zero crossing (for example, at I = 0.3 a.u. it evaluates to approximately -0.13). As written, it cannot represent an ionization probability over the lower range of Fig. 5. If the fit is intended only for the OBI-active range above threshold, the text should state the fitting domain and the number of data points used; otherwise a bounded functional form should be adopted.","section":"Eq. (3), fitting function"}],"minor_comments":[{"comment":"The caption refers to 'yellow solid lines' while the body text says 'orange solid lines'. Please unify.","section":"Fig. 4 caption vs. text"},{"comment":"The solving of Eq. (5) is stated without the numerical method or the value of I_p used for the model helium atom. Given that I_p is needed to reproduce F_Stark_th = 0.222, please state it explicitly.","section":"Eq. (5) and threshold values"},{"comment":"The distinction between the tunnel exit in Fig. 2(a1) and the barrier top in Fig. 2(b1) would be clearer if the coordinates of the ion and the field direction were marked on both panels.","section":"Notation for barrier top in Fig. 2(b)"},{"comment":"The sentence 'the PMD for TI encircles that for OBI' is qualitative; adding an angular or radial histogram with error bars would strengthen the claim.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is promising and addresses a real gap, but the central trajectory tag is not yet a validated observable. The referee report should ask for an explicit algorithm, a convergence test, and a quantitative comparison with the effective-saddle criterion. If the authors can provide that, the paper may be acceptable as a Letter; without it, the central claims are not yet supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is adapting backpropagation to the over-barrier regime and classifying trajectories by topology: TI trajectories rebound at the tunnel exit, OBI trajectories return to the ion and orbit it. That split, plus the PMD partitioning in Fig. 3, is a real step beyond the tunneling-only backpropagation work in refs 40-43. The claim that the TI/OBI border moves with initial transverse momentum is also not in the cited OBI literature. The internal pipeline - TDSE forward, classical backward - is coherent, and the exact sum of TI and OBI PMDs is a nice check. The soft spots are exactly where the reader's stress test lands, and I think it holds up. The OBI tagging rule in Fig. 2 is qualitative: 'orbits around the ion' with the barrier top at a local speed minimum. No tolerances, no convergence tests, and no comparison against the paper's own energy-over-saddle criterion. The worry is not hypothetical - orbiting kinematically favors small angular momentum, which is near-synonymous with small k_perp. So the finding that OBI electrons sit in the central arc may be partly pre-encoded by the tagging rule. The paper also says the border is sharp because the classification is binary, which is close to circular. The sum check only proves exhaustiveness, not mechanism attribution. The other mechanical issue is the threshold. The fit P_OBI = -0.288 I^-1.758 + 0.809 goes negative over part of the plotted range, so the zero at I_OBI = 0.556 a.u. is not a robust number, and the agreement with the Stark-shifted threshold (0.537 a.u.) is weaker than the text implies. That said, the Stark-shift correction is a sensible idea and the direction of the effect is probably right. None of this is fatal. The underlying TDSE numerics seem sound, and the classification method is plausible. What's missing is specification and validation. The authors should release the exact algorithm, show the tag is stable under reasonable tolerance choices, and test it against an independent diagnostic. The fit should use a form that is nonnegative over the fitting range. Who gets value: anyone working on strong-field ionization, attoclock measurements, or trajectory-based analysis of ionization dynamics. It would be a useful contribution once the method is nailed down. My recommendation: send it to peer review, but the reviewers should push for the classification details and a physical fit. Right now the paper is a good idea with under-supported quantitative claims.","headline":"A promising but underspecified extension of backpropagation to over-barrier ionization; the central TI/OBI separation needs a precise algorithm and an independent validation before its headline claims are fully supported.","tokens_in":793,"tokens_out":715,"would_cite":false,"duration_ms":22232,"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":"Backpropagated electron trajectories separate over-barrier from tunneling ionization, showing the boundary depends on transverse momentum and the threshold needs the Stark-shifted energy.","keywords":["over-barrier ionization","tunneling ionization","backpropagation","strong-field ionization","photoelectron momentum distribution","ionization time","Stark shift","trajectory topology"],"falsifier":"A reader could settle the central claim by taking the computed backpropagated trajectories and checking, trajectory by trajectory, whether the topological OBI label (orbit around the ion, local speed minimum) coincides with the independent saddle condition E_ground >= max_r[(1/2)k_perp^2 - 1/r + r·F_c] for the same instantaneous field. If the two diagnostics disagree systematically near the boundary, or if the recovered k_perp values jump discontinuously at the OBI/TI border, the topological tagging rule is not a faithful mechanism label.","tokens_in":9381,"feed_emoji":"⚛️","tokens_out":7639,"duration_ms":79522,"temperature":0.7,"pith_summary":"This paper tries to establish a trajectory-level dynamical picture of over-barrier ionization (OBI) in strong-field ionization of helium, a regime that standard tunnel-ionization theories cannot describe because they use a zero-range potential. Using backpropagation of classical trajectories launched from a numerical time-dependent Schrödinger wavefunction under the full Coulomb Hamiltonian, it claims that ionizing electrons split into two topological classes: tunneling-ionization (TI) trajectories, which rebound at the tunnel exit where the velocity along the field direction vanishes, and OBI trajectories, which return to the parent ion and orbit it, with the barrier top at a local speed minimum. From this split, the paper derives a clean decomposition of the photoelectron momentum distribution and ionization-time distribution, shows that the OBI/TI boundary depends on initial transverse momentum as well as field strength, reveals a competitive intensity dependence of the two probabilities, and argues that the threshold field must use the Stark-shifted binding energy. A sympathetic reader would care because this is the first dynamic, not merely rate-based, account of OBI in a full-Coulomb model, and it produces concrete, testable features in momentum and timing observables.","feed_headline":"Trajectory topology splits tunneling from over-barrier ionization","feed_subtitle":"OBI electrons orbit the ion and occupy a central momentum arc; the threshold intensity requires the Stark shift.","key_machinery":"The central object is the backpropagated classical trajectory obtained from the ionized wavefunction under the full Hamiltonian, classified by trajectory topology. For TI, the stopping criterion is k_parallel = 0, the vanishing of velocity along the instantaneous field direction at the tunnel exit. For OBI, the stopping criterion is that the backpropagated electron returns close to the parent ion and orbits it, with the barrier top set at a local minimum of speed before the orbital motion. This binary topological classification is what converts a quantum ionization calculation into separately labeled photoelectron momentum distributions and ionization-time distributions, and it is what makes","core_discovery":"On its own terms, the central claim is that over-barrier ionization is not just a rate threshold but a distinct topology of electron motion in the full Coulomb-plus-laser Hamiltonian. The method works by propagating the ionized part of the wavefunction backward in time as classical trajectories until a stopping criterion is reached. TI trajectories are identified by the velocity criterion k_parallel = 0 at the tunnel exit, while OBI trajectories are identified by the electron returning near the parent ion and orbiting it, with the barrier top located at the local minimum of speed just before the orbit begins. Because the classification is binary, the TI/OBI boundary is sharp. The resulting m","pith_inferences":["Editorial inference: the same topological tagging could be transferred to three-dimensional atoms, where 'orbiting the parent ion' corresponds to low angular momentum about the ion; OBI electrons should then carry small magnetic quantum number, testable in magnetic-sublevel-resolved momentum tomographies.","Editorial inference: the stranded-trajectory singularity at the OBI/TI boundary implies a sharp, subcycle discontinuity in the recovered ionization time as the intensity crosses the threshold; an attoclock scan across the threshold could look for this feature.","Editorial inference: the quantitative agreement between the numerical threshold and the Stark-shifted condition depends on the model polarizability alpha = 1.57; repeating the extraction for atoms or molecules with known polarizabilities would show whether the Stark-shift explanation is generic.","Editorial inference: the backpropagation labels could be applied to trajectories that rescatter or drive high-harmonic generation, potentially separating OBI and TI contributions in those observables rather than only in direct ionization."],"forward_implications":["OBI electrons occupy a central arc in the photoelectron momentum distribution while TI electrons form an encircling outer arc, so channel-resolved momentum spectra become separable.","The OBI/TI boundary depends on initial transverse momentum as well as field strength: even above the barrier-suppression intensity, large-transverse-momentum electrons still ionize by tunneling.","After depletion correction, TI and OBI ionization-time distributions are symmetric about the pulse peak with no significant relative delay, except for a peak/dip feature at the boundary caused by stranded trajectories.","As intensity rises, OBI probability increases while TI probability decreases, so the two mechanisms compete; their contributions become comparable near the threshold.","Accurate threshold determination for OBI requires the Stark-shifted binding energy: F_Stark_th = 0.222 a.u. matches the numerically inferred intensity threshold, while the unshifted value F_th = 0.204 a.u. does not.","The sum of the separately labeled TI and OBI momentum distributions reproduces the total photoelectron momentum distribution, confirming that the trajectory classification exhaustively partitions the ionized population."],"supporting_citations":[{"why":"Supplies the original backpropagation method that converts the ionized wave packet into classical trajectories integrated backward in time.","marker":"[40]"},{"why":"Extends backpropagation to extract tunnel-exit characteristics and notes that the velocity criterion fails when the nontunneled fraction is high, motivating new OBI-specific criteria.","marker":"[41]"},{"why":"Supplies the k_parallel = 0 velocity criterion and the nonadiabatic analysis used for TI trajectories, and the numerically determined polarizability alpha = 1.57 used in the Stark-shift threshold equation.","marker":"[42]"},{"why":"ADK theory and its anatomy; grounds the initial transverse momentum k_perp at the tunnel exit and therefore the effective potential barrier model E = (1/2)k_perp^2 - 1/r + r·F_c.","marker":"[13, 14]"},{"why":"Gives the simple Coulomb-barrier threshold F_th = I_p^2/(4Z) that the paper extends by including the Stark shift.","marker":"[19]"},{"why":"Show that including the Stark shift in the binding energy improves threshold estimates, which the paper turns into a quantitative threshold condition.","marker":"[22, 23]"},{"why":"Supplies the single-active-electron model used for the helium target in the numerical simulation.","marker":"[30]"},{"why":"Provides the depletion-correction formula used to obtain the symmetric ionization-time distributions.","marker":"[51]"},{"why":"Describe stranded trajectories at the effective barrier top, used to explain the peak/dip structure at the TI/OBI boundary in the ionization-time distribution.","marker":"[52, 53]"}],"fun_headline_variants":["Orbiting electrons mark over-barrier ionization","Stark shift sets over-barrier ionization threshold","Transverse momentum steers ionization mechanism","Backpropagation reveals over-barrier's electron orbits","Tunneling and over-barrier ionization: topology splits them"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The load-bearing premise is that an OBI trajectory is reliably identified as one which, when backpropagated, returns close to the parent ion and orbits it, with the barrier top at a local minimum of speed; this rule is stated without tolerances and already steers OBI electrons toward small transverse momentum, so if it mislabels trajectories, the momentum decomposition, the transverse-momentum dependence claim, the probability curves, and the Stark-shift threshold conclusion","fun_headline_variants_meta":{"raw":{"variants":["Orbiting electrons mark over-barrier ionization","Stark shift sets over-barrier ionization threshold","Transverse momentum steers ionization mechanism","Backpropagation reveals over-barrier's electron orbits","Tunneling and over-barrier ionization: topology splits them"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000369,"raw_usage":{"total_tokens":1795,"prompt_tokens":703,"completion_tokens":1092,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":447,"completion_tokens_details":{"reasoning_tokens":1019}},"tokens_in":447,"tokens_out":1092,"duration_ms":11183,"temperature":1.0,"reasoning_tokens":1019,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T12:00:13.072643+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A reader could settle the central claim by taking the computed backpropagated trajectories and checking, trajectory by trajectory, whether the topological OBI label (orbit around the ion, local speed minimum) coincides with the independent saddle condition E_ground >= max_r[(1/2)k_perp^2 - 1/r + r·F_c] for the same instantaneous field. If the two diagnostics disagree systematically near the boundary, or if the recovered k_perp values jump discontinuously at the OBI/TI border, the topological tagging rule is not a faithful mechanism label.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the original backpropagation method that converts the ionized wave packet into classical trajectories integrated backward in time."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Extends backpropagation to extract tunnel-exit characteristics and notes that the velocity criterion fails when the nontunneled fraction is high, motivating new OBI-specific criteria."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the k_parallel = 0 velocity criterion and the nonadiabatic analysis used for TI trajectories, and the numerically determined polarizability alpha = 1.57 used in the Stark-shift threshold equation."},{"cited_title":"Augst, D","cited_arxiv_id":null,"evidence_quote":"Gives the simple Coulomb-barrier threshold F_th = I_p^2/(4Z) that the paper extends by including the Stark shift."},{"cited_title":"hole” in the trajectory at the ion, which ef- fectively makes the trajectory “topological","cited_arxiv_id":null,"evidence_quote":"Supplies the single-active-electron model used for the helium target in the numerical simulation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the depletion-correction formula used to obtain the symmetric ionization-time distributions."}],"review_version":1}