{"id":"10cf43e9-4834-41d3-bc16-7793bb6fa57a","arxiv_id":"2501.08106","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A time-dependent Schrödinger equation model for H- photodetachment reproduces the standard cross-section result at moderate laser power but diverges at high power.","lead":"This paper computes how efficiently a laser can knock one electron off a negative hydrogen ion beam, and compares a full time-dependent quantum calculation with a simpler cross-section formula. It matters for designing laser-based stripping systems for particle accelerators, where replacing a magnet with a laser could simplify injection.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'full' TDSE reference is only a single-channel model with an empirical dipole; its agreement with Eq. (19) does not establish physical photodetachment at E0=0.01 a.u.","rationale":"The reader's weakest assumption is essentially correct, but I would sharpen it. The problem is not only that other channels are missing; the one continuum channel used is not an eigenstate of the full Hamiltonian. Equation (13) ignores electron-electron repulsion in the final state. Therefore Eq. (6) is not the full TDSE, and the phrase 'precise method' overstates the reference. The agreement of Eq. (19) with Eq. (6) is partly trivial because Eq. (19) is derived from the same equations with a(t)=1; the only missing effect in Eq. (19) is depletion. Without reporting the efficiency window and grid convergence, the 10% figure is not a physical benchmark. I do not see an internal algebraic error in Eq. (19): with the rotating-wave approximation it follows from Eq. (6), and the prefactor checks against Eq. (8). The conditional verdict is appropriate; the authors should add the extended-channel and convergence checks. Hence the verdict is unchanged.","tokens_in":6800,"tokens_out":13037,"duration_ms":137917,"concrete_test":"Perform a converged two-electron TDSE calculation in a basis that includes the single-ionization channels H(nl)+e−, the double-ionization continuum, and, at the Feshbach photon energy, the doubly-excited resonance, for the same laser parameters (E0=0.01 a.u., ω=2Eb). Compute the cumulative norm in all channels outside the single-channel ansatz (3) at the final time of Fig. 6. If that omitted probability exceeds 10% of the total, the one-channel efficiencies and the claimed 10% agreement do not describe physical photodetachment; if it is below 1%, the truncation is validated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (3) restricts the wavefunction to the H− ground state plus one continuum channel H(1s)+e−, and Eqs. (13)–(15) represent that continuum as a symmetrized product of the hydrogen 1s orbital and a free P-wave. This ψ_E is not an eigenfunction of the two-electron Hamiltonian (2): the 1/|r1−r2| term is dropped in the continuum channel entirely. Consequently, the projection leading to Eq. (6) does not solve Eq. (1); it defines a single-active-electron model. The dipole element (18) is not computed from these wavefunctions but fitted to the zero-field cross section (17), so the energy dependence of the coupling is taken from low-intensity perturbation theory and inserted into a nonlinear equation at E0=0.01 a.u. (3.5×10^12 W/cm2), where the outer electron binding energy is only 0.0274 a.u. Equation (19) is the a(t)=1 limit of the same truncated equations with the same fitted μ_E, so the reported <10% agreement is internal consistency of the model, not validation of the predicted efficiency. The paper's own Feshbach section concedes that a relevant bound channel is omitted, and no estimate is given for double-ionization or continuum-continuum amplitude.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a time-dependent Schrödinger-equation model for photodetachment of the negative hydrogen ion H^- + γ → H^0 + e^-, motivated by the laser-assisted charge exchange injection (LACE) program at the Spallation Neutron Source. The authors expand the two-electron wavefunction as the H^- ground state plus a single continuum channel H(1s) + e^-, derive coupled integro-differential equations (Eq. 6), and use a dipole matrix element that is empirically fitted to the measured zero-field photodetachment cross section (Eqs. 17-18). They compare three approaches: a linear exponential attenuation model (Eq. 10), a first-order analytic approximation with a(t)=1 (Eq. 19), and a numerical solution of the coupled equations (Eq. 6). They report that the linear model agrees with the wave-equation result at low power (E0=0.001 a.u.) but disagrees at E0=0.01 a.u., while the first-order approximation agrees with the full numerical solution to within about 10%. The Feshbach-resonance detachment route is discussed qualitatively but is not calculated.","tokens_in":7030,"tokens_out":4215,"duration_ms":45571,"significance":"If substantiated, the proposed time-dependent wave-equation method could be a useful design tool for estimating laser power requirements for H^- photodetachment in accelerator applications, going beyond the standard linear cross-section model. The manuscript has the strength of explicitly comparing linear and nonlinear formulations and of providing an analytic first-order expression for the energy-dependent dipole case. However, the central numerical claim is not yet validated: the 'full' calculation is a single-channel model with an empirically fitted dipole, and the agreement between Eq. (6) and Eq. (19) is, to a large extent, built into the model construction. The paper would be significantly strengthened by a convergence study, an ab initio dipole calculation, and a quantitative statement of the model's validity domain.","major_comments":[{"comment":"The wavefunction ansatz (3) and the continuum channel form (13)-(15) constitute a single-active-electron model: the continuum function is a symmetrized product of the hydrogen 1s orbital and a free P-wave, and the e^2/|r1-r2| interaction is not included in the continuum channel. This is not an eigenfunction of the full two-electron Hamiltonian (2). At the photon energy used (ω=2Eb=0.0548 a.u.), the Feshbach state and excited hydrogen bound states are not energetically open, so their omission may be acceptable at low intensity, but at E0=0.01 a.u. multiphoton and continuum-continuum processes could populate high-energy continuum states; no estimate of their amplitude is given. The claim that Eq. (6) is a 'precise' solution of Eq. (1) is therefore overstated; it is a truncated model whose truncation error is unquantified.","section":"MODEL OF H− PHOTODETACHMENT, Eqs. (2)-(3) and (13)-(15)"},{"comment":"The dipole matrix element μ_E is not computed from the wavefunctions (11) and (13) but is inverted from the empirical cross-section fit (17) to the low-intensity experimental data of Ref. [8]. Consequently, Eq. (6) is a nonlinear version of the same empirical model that produces the linear cross-section (10). The agreement between Eq. (6) and Eq. (19) shown in Figure 6 is therefore an internal consistency check of the assumed μ_E, not a validation of the predicted photodetachment efficiency at high intensity. The manuscript should either compute μ_E from the stated wavefunctions or compare the model's intensity-dependent predictions (e.g., detachment probability versus pulse energy) against independent experimental data.","section":"PHOTODETACHMENT CROSSECTION, Eq. (18)"},{"comment":"The numerical solution of the coupled equations is described only as discretization of the continuum and integration from E=0 to E_max=1, with no information on the number of grid points, time step, integration method, or convergence criteria. Without a convergence study, the reported 'agreement to within 10%' in Figure 6 is not quantitatively supported. The choice E_max=1 a.u. is also arbitrary; the manuscript should show that results are insensitive to this cutoff and to the discretization density.","section":"PHOTODETACHMENT CALCULATIONS, numerical solution of Eq. (6)"},{"comment":"Equation (19) is derived under the assumption a(t)=1, which is valid only when the ground-state depletion is small. At E0=0.01 a.u., the detachment probability is large enough (as seen in Figure 6) that a(t) must deviate significantly from unity, so the claimed 10% agreement between Eq. (19) and the full solution of Eq. (6) is surprising and unexplained. The authors should either provide an analytic argument for why the a(t)=1 approximation remains accurate at high intensity, or clarify the time window and efficiency range over which this agreement holds.","section":"PHOTODETACHMENT CALCULATIONS, Eq. (19) and Figure 6"}],"minor_comments":[{"comment":"The manuscript contains numerous typographical errors and inconsistent spacing (e.g., 'Schrodinger' for 'Schrödinger', 'photodetachmnet', 'crossection', 'thehydrogenion'), which should be corrected before submission.","section":"Throughout"},{"comment":"The continuum integral is written from 0 to ∞ in Eq. (3) and in the text, but the numerical calculation truncates at E_max=1; this inconsistency should be explained and the truncation justified.","section":"MODEL OF H− PHOTODETACHMENT, Eq. (3) and following text"},{"comment":"Figure captions for Figures 5 and 6 give laser power density in MW/mm^2, which is an unusual unit; the authors should verify and, if possible, state the intensity in W/cm^2 consistently with the main text.","section":"PHOTODETACHMENT CALCULATIONS"},{"comment":"The Feshbach-resonance discussion is only qualitative and no calculation is presented; if this mechanism is central to the LACE goal, the manuscript should at least provide an estimate of the required laser power or state explicitly that such a calculation is future work.","section":"FESHBACH RESONANCE, Eq. (20)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript from ORNL has an accelerator-physics motivation and would be of interest to a specialized audience. However, in its current form it is closer to a workshop/proceedings contribution than to a fully developed archival paper: the numerical details are sparse, the dipole is empirical, and the central comparison is an internal consistency check. The authors should be encouraged to add a convergence study, an ab initio or independent validation of the dipole, and a clear statement of the model's limitations. The paper is not suitable for acceptance without these additions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know: this paper is not a breakthrough in atomic physics and does not claim to be. It is an interim report from the SNS group on laser-based stripping for charge exchange injection. The new content is applying a standard single-channel time-dependent Schrödinger equation to LACE-relevant laser powers and showing that the widely used linear cross-section formula (10) significantly overestimates photodetachment efficiency at E0 = 0.01 a.u. (35,000 MW/mm2). That warning is real and useful for accelerator design.\n\nThe paper does several things well. The coupled equations (6) are derived cleanly, the first-order approximation (19) is a sensible improvement over the constant-dipole linear model, and the authors are explicit that the full two-electron problem, including the Feshbach resonance that would actually make the scheme practical, is not yet solved. The comparison between the numerical TDSE and the analytic approximation is a nice internal consistency check.\n\nThe soft spots are mostly the ones the stress-test note flags, and they are genuine. Equation (3) restricts the wavefunction to one continuum channel; the electron-electron repulsion is dropped in that channel, so this is effectively a single-active-electron model. The dipole matrix element (18) is not computed from the wavefunctions; it is inverted from the experimental cross section, so the energy dependence is imported from low-intensity perturbation theory and inserted into a nonlinear equation at high intensity. The agreement between (19) and the numerical solution is therefore internal consistency of the model, not validation against experiment. Add to that the missing numerical details: no convergence tests, no discretization parameters, no code or data supplied. These are fixable in revision.\n\nThat said, the paper is honest about its own scope. The single-channel assumption is stated up front, the Feshbach section explicitly lists what is missing, and the summary says the main LACE goal still requires more work. For an accelerator physicist designing a laser stripping scheme, the conclusion that the linear model is unsafe at high intensity is worth having, even if the quantitative efficiencies are model-dependent.\n\nI would send this to peer review. It is not ready as is, but with convergence tests, a clearer description of the numerical method, and perhaps a benchmark against known photodetachment yields, it becomes a legitimate contribution to the LACE literature. It does not deserve a desk reject.","headline":"A useful engineering warning that the linear cross-section model fails at high laser power, wrapped in a single-channel TDSE calculation with an empirical dipole; the promised Feshbach mechanism is deferred.","tokens_in":7580,"tokens_out":1603,"would_cite":true,"duration_ms":18621,"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":"The paper argues that for high-power laser photodetachment of $\\mathrm{H}^-$, the standard linear cross-section model fails, while a first-order analytic approximation from the time-dependent wave equation stays within 10% of the full…","keywords":["H- photodetachment","laser-assisted charge exchange injection","time-dependent Schrodinger equation","photodetachment cross section","Feshbach resonance","dipole matrix element","strong-field laser stripping"],"falsifier":"A pump-probe measurement of the neutral hydrogen yield from a pulsed laser with $E_0=0.01$ a.u. at $\\omega=2E_b$ would settle the matter: if the measured detachment fraction follows the linear formula (10) instead of the wave-equation curves (6) and (19), or if it deviates from (6) by more than 10%, the central claim is wrong. A numerical check with a double-ionization continuum channel added to (6) would show whether the omitted channels break the single-channel truncation at that field strength.","tokens_in":1834,"feed_emoji":"⚡","tokens_out":2349,"duration_ms":82339,"temperature":0.7,"pith_summary":"This paper tackles a practical beam-physics problem: stripping an electron off a negative hydrogen ion with a laser alone, instead of a stripping magnet, for laser-assisted charge exchange injection into an accelerator. The authors compute photodetachment efficiency with the full time-dependent Schrödinger equation for the two-electron system, truncated to the ground state of $\\mathrm{H}^-$ plus the $\\mathrm{H}^0+e^-$ continuum. They find that the standard linear cross-section model matches the full calculation at moderate laser power but disagrees substantially at high power ($E_0=0.01$ a.u., about 35 000 MW/mm$^2$). A first-order analytic approximation that keeps the energy dependence of the dipole matrix element agrees with the full wave-equation result to within 10% of photodetachment efficiency. The paper also lays out the Feshbach-resonance route through a doubly excited $\\mathrm{H}^-$ state as the promising path to the >99% efficiency that practical injection would need.","feed_headline":"Standard H- stripping model fails at high laser power","feed_subtitle":"A first-order analytic wave-equation approximation stays within 10% of the full calculation","key_machinery":"The machinery is the two-electron time-dependent Schrödinger equation with the ansatz $\\Psi(\\mathbf{r}_1,\\mathbf{r}_2,t) = a(t)\\psi_{\\mathrm{H}^-} e^{-iE_0t} + \\int_0^\\infty c_E(t)\\psi_E e^{-i(E_{1s}+E)t} dE$, which keeps only the ground state of $\\mathrm{H}^-$ and a single $\\mathrm{H}^0+e^-$ continuum channel. Substitution gives the integro-differential system (6) for $a(t)$ and $c_E(t)$. The load-bearing object is the dipole matrix element $\\mu_E = \\langle\\psi_E|z_1+z_2|\\psi_{\\mathrm{H}^-}\\rangle$, supplied by inverting the empirical cross-section fit (17) of measured photodetachment data, so the nonlinear calculation uses measured physics rather than a constant $\\mu$. The first-order analytic formula (19) keeps the energy integral and the finite lower bound of the continuum, which is why it outperforms the constant-$\\mu$ linear approximation (8)--(10).","core_discovery":"The central claim is that photodetachment probability of $\\mathrm{H}^-$ under a strong laser is accurately described by solving the coupled integro-differential equations (6) for the bound-state amplitude $a(t)$ and the continuum amplitudes $c_E(t)$, with the dipole matrix element $\\mu_E$ taken from the empirical cross-section fit, and that the standard exponential formula $1-\\exp(-N_\\gamma \\sigma t)$ is not valid in the high-power regime. At $E_0=0.01$ a.u. the linear model disagrees significantly with the wave-equation result, while the first-order analytic expression (19)---which integrates the sinc-squared resonance factor over the continuum with the full energy-dependent $\\mu_E$---reproduces the full calculation to within 10% of photodetachment efficiency. The paper therefore presents (19) as a cheap, accurate replacement for full numerics in parameter scans, and identifies the Feshbach resonance, not direct detachment, as the mechanism that would need to be engineered for practical >99% laser stripping.","pith_inferences":["If the 10% agreement between (19) and (6) holds over a wider range of frequencies and pulse shapes, one could invert (19) to design laser pulses that hit a target detachment fraction without any full wave-equation runs.","The same single-channel machinery could be applied to other negative ions, such as $\\mathrm{D}^-$, by replacing the $\\mathrm{H}^-$ wave function and dipole matrix element, giving quick estimates of whether laser stripping is feasible at their beam energies.","A natural next step is to add the Feshbach bound state as a second discrete channel in (20) and solve the coupled equations; the paper sets up that formalism but leaves the calculation to future work.","Because the single-channel truncation omits double ionization and excited hydrogen channels, an experiment at 35 000 MW/mm$^2$ might see additional $\\mathrm{H}^+$ or proton yield that the model cannot represent, which would set the boundary of the approximation."],"forward_implications":["At moderate laser power ($E_0=0.001$ a.u., about 350 MW/mm$^2$), the standard linear cross-section formula remains adequate for estimating $\\mathrm{H}^-$ photodetachment efficiency.","At high power ($E_0=0.01$ a.u., about 35 000 MW/mm$^2$), design calculations that rely on the linear model will misestimate efficiency, so the full wave equation or the analytic approximation (19) must be used.","Expression (19) gives a fast surrogate for the full numerical solution, making laser-power requirement scans practical without solving the integro-differential system each time.","Direct photodetachment still requires extremely high power density, so the Feshbach-resonance path through the doubly excited $\\mathrm{H}^-$ state is the route that could plausibly reach the >99% stripping efficiency needed for charge-exchange injection.","If the single-channel model is confirmed experimentally, laser-only stripping becomes a viable replacement for the first stripping magnet in injection systems, simplifying the accelerator layout."],"supporting_citations":[{"why":"Supplies the measured photodetachment cross section that the empirical fit (17) is normalized to, fixing the dipole matrix element $\\mu_E$.","marker":"[8]"},{"why":"Provides the empirical cross-section formula (17) used to define $\\mu_E$ for the wave-equation calculations.","marker":"[10]"},{"why":"Supplies the high-accuracy 20-parameter $\\mathrm{H}^-$ ground-state wave function used as the initial bound state.","marker":"[14]"},{"why":"Gives the energy-dependent photodetachment cross section and the resonance peaks that motivate the Feshbach path.","marker":"[11]"},{"why":"Identifies the Feshbach resonance of $\\mathrm{H}^-$ and its doubly excited state as the high-efficiency detachment mechanism.","marker":"[12]"},{"why":"Provides the experimental observation of the Feshbach resonance near $h\\nu \\approx 10.93$ eV used for the resonance parameters.","marker":"[15]"},{"why":"Supplies the time-dependent perturbation theory baseline that the linear model (8)--(10) is built on.","marker":"[9]"}],"fun_headline_variants":["Standard H- stripping model fails at high laser power","New analytic method beats standard H- photodetachment model","H- photodetachment: first-order wave equation within 10%","Laser stripping of H-: beyond the exponential model","Strong-laser H- detachment: simple formula rivals full numerics"],"cache_read_input_tokens":9728,"weakest_assumption_plain":"The calculation assumes the evolving wavefunction contains only the $\\mathrm{H}^-$ ground state plus a single continuum channel $\\mathrm{H}^0+e^-$; if double ionization, excited hydrogen channels, or the Feshbach state carry significant amplitude at the laser powers studied, the computed efficiencies and the claimed 10% agreement would not describe the real process.","fun_headline_variants_meta":{"raw":{"variants":["Standard H- stripping model fails at high laser power","New analytic method beats standard H- photodetachment model","H- photodetachment: first-order wave equation within 10%","Laser stripping of H-: beyond the exponential model","Strong-laser H- detachment: simple formula rivals full numerics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000124,"raw_usage":{"total_tokens":1060,"prompt_tokens":855,"completion_tokens":205,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":471,"completion_tokens_details":{"reasoning_tokens":120}},"tokens_in":471,"tokens_out":205,"duration_ms":2931,"temperature":1.0,"reasoning_tokens":120,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:29:55.930348+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A pump-probe measurement of the neutral hydrogen yield from a pulsed laser with $E_0=0.01$ a.u. at $\\omega=2E_b$ would settle the matter: if the measured detachment fraction follows the linear formula (10) instead of the wave-equation curves (6) and (19), or if it deviates from (6) by more than 10%, the central claim is wrong. A numerical check with a double-ionization continuum channel added to (6) would show whether the omitted channels break the single-channel truncation at that field strength.","supporting_citations":[{"cited_title":"Smith and D","cited_arxiv_id":null,"evidence_quote":"Supplies the measured photodetachment cross section that the empirical fit (17) is normalized to, fixing the dipole matrix element $\\mu_E$."},{"cited_title":"Armstrong, Empirical analysis of the H− photodetach- ment cross section, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the empirical cross-section formula (17) used to define $\\mu_E$ for the wave-equation calculations."},{"cited_title":"Hart and G","cited_arxiv_id":null,"evidence_quote":"Supplies the high-accuracy 20-parameter $\\mathrm{H}^-$ ground-state wave function used as the initial bound state."},{"cited_title":"Broad and W.P","cited_arxiv_id":null,"evidence_quote":"Gives the energy-dependent photodetachment cross section and the resonance peaks that motivate the Feshbach path."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the Feshbach resonance of $\\mathrm{H}^-$ and its doubly excited state as the high-efficiency detachment mechanism."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the experimental observation of the Feshbach resonance near $h\\nu \\approx 10.93$ eV used for the resonance parameters."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the time-dependent perturbation theory baseline that the linear model (8)--(10) is built on."}],"review_version":1}