{"id":"d1fae7b0-f2e0-4906-ae81-96d015285bbe","arxiv_id":"2412.17996","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"For a model atom in a flat-top laser pulse, the nondipole shift of multiphoton photoelectron peaks depends on emission angle, dominated by electron recoil, with only a small redshift from retardation.","lead":"This paper calculates how a laser pulse's forward push changes the directions and energies of electrons knocked out of a model atom. It separates two nondipole effects, electron recoil and field retardation, and shows recoil causes the largest directional shifts.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The analytic path to Eq. (36) is internally inconsistent: Eq. (33)'s recoil term uses n·p = p0 − p2, which duplicates the ponderomotive phase and reverses the sign of the 1/c shift, so Eq. (36) does not follow from Eqs. (33)–(34) as printed.","rationale":"The reader's conditional verdict is appropriate, but my main concern differs from the reader's weakest assumption about finite-pulse validity of the quasi-energy picture. The flat-top pulse with 18 nearly constant cycles makes that assumption reasonable. The more concrete problem is that Eq. (33) as printed cannot generate Eq. (36): the recoil term n·p = p0 − p2 contains a rest-mass contribution that duplicates the explicit Φ^(0) term and gives the O(1/c) correction with the wrong sign. An independent expansion of Eq. (32) does yield Eq. (36) and matches the sign of the TDSE stripes, so the physics conclusion is probably right despite the flawed printed derivation. This strengthens the case for CONDITIONAL rather than ACCEPT or REJECT: the numerics and final formula are plausible, but the analytic derivation must be corrected and the sign/notation issue resolved before the attribution can be taken as rigorously established. The numerical method (Suzuki-Trotter, with reported convergence parameters) is a positive feature; no evidence of a fatal numerical flaw was found. This is a technical inconsistency in the manuscript text, not a criticism of the authors' intent or integrity.","tokens_in":14172,"tokens_out":19208,"duration_ms":186736,"concrete_test":"Independently re-derive Eq. (36) directly from the exact Volkov phase Eq. (32) by expanding c/(p0 − p2) to first order in 1/c for a constant-amplitude field, without using Eq. (33). If the coefficient of p2 in the time-linear phase is −U_p/(m_e c), Eq. (36) survives and only Eq. (33) needs correction; if the sign is +U_p/(m_e c), the central recoil attribution would need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central attribution 'recoil dominates the directional shift' rests on Eq. (36), whose stated derivation uses Eqs. (33) and (34). As printed, Eq. (33) has the term (i/m_e c) Φ_p^(0)(t) n·p with n·p = p0 − p·n (defined before Eq. (32)). For nonrelativistic electrons p0 ≈ m_e c + E_p/c, so n·p/(m_e c) ≈ 1 + E_p/(m_e c^2) − p2/(m_e c). Substituting this into Eq. (33) adds an extra full Φ^(0) phase beyond the explicit Φ^(0) already present, and gives an O(1/c) term with the opposite sign to the exact expansion of c/(n·p) in Eq. (32). The correct expansion of Eq. (32) for a constant-amplitude field yields the phase correction (Φ^(0)/(m_e c))(p2 − E_p/c), which leads to E_n(φ_p) = E_n − Ω_n sin φ_p, i.e., Eq. (36). Thus Eq. (36) is believable, but not because of Eqs. (33)–(34) as written; the printed derivation is invalid. If the typo were taken literally, the predicted angular shift would have the opposite sign, contradicting both Eq. (36) and the TDSE stripes. The load-bearing step is therefore the unstated corrected expansion, and the paper should show it.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies nondipole effects in multiphoton ionization of a two-dimensional hydrogen-like model atom driven by a linearly polarized flat-top laser pulse propagating along x2. The authors solve the time-dependent Schrödinger equation with a Suzuki-Trotter split-step method, treating the propagating vector potential A(t - x2/c) exactly, and compare the results with those obtained in the dipole approximation. For moderate and strong intensities, they observe that the multiphoton stripes in the energy-angle photoelectron distributions become tilted, and individual peaks exhibit angle-dependent red and blue shifts. They attribute the tilt to the electron recoil (Nordsieck) correction and the small polarization-direction shift to retardation, quantified by Eq. (36), E_n(phi_p) = E_n - U_p sqrt(2 m_e E_n)/(m_e c) sin(phi_p). They also report double-hump structures, angular distribution asymmetries, and the gradual disappearance of multiphoton peaks for strong propagating pulses.","tokens_in":14516,"tokens_out":28699,"duration_ms":253427,"significance":"The numerical treatment of a genuinely propagating laser pulse beyond the dipole approximation is a useful contribution, and the central directional effect is directly visible in the comparison between propagating-pulse and dipole distributions. The paper's main strength is that it produces a simple, falsifiable analytic prediction, Eq. (36), for the angle-dependent peak shift. If the derivation of that formula is repaired and the post hoc fitting in the strong-field section is addressed, the conclusion that the recoil term governs the directional dependence of photoelectron energy spectra would be of significant interest to the strong-field ionization community. The explicit numerical parameters and the machine-checkable split-step scheme are also positive features.","major_comments":[{"comment":"The derivation of Eq. (36) is not valid as printed. With n·p = p0 − p·n defined before Eq. (32), the recoil term in Eq. (33) reads (i/m_e c) Phi_p^(0)(t) (p0 − p·n). Inserting p0 = m_e c + E_p/c + ... produces a zeroth-order contribution i Phi_p^(0)(t), which doubles the nonrelativistic Volkov phase, and a first-order contribution −(i/m_e c) Phi_p^(0)(t)(p·n − E_p/c), which has the opposite sign to the first-order term obtained by directly expanding Eq. (32). Thus Eq. (36) does not follow from Eqs. (33)-(34) as written. The authors should either correct the definition of n·p in Eq. (33), present the actual expansion of Eq. (32), or derive Eq. (36) directly from the nonrelativistic TDSE, and they should state the sign convention explicitly.","section":"Sec. IV, Eqs. (32)-(36)"},{"comment":"The strong-pulse comparison is not a parameter-free test of Eq. (36). The dipole reference energies E^(d)_1 = 95.8E0 and E^(d)_2 = 99.8E0 are first used to draw the white lines in Fig. 7(b); when these do not follow the stripes, the authors change to E^(p)_1 = 95.5E0 and E^(p)_2 = 99.5E0 for the red lines. The zero-angle offset is therefore adjusted to the propagating-pulse data. This tests only the sinusoidal angular shape and the tilt amplitude, not the absolute peak positions. The authors should either compute the expected zero-angle redshift from the retardation term or explicitly state that E^(p)_1 and E^(p)_2 are fit parameters.","section":"Sec. V, Fig. 7"},{"comment":"No quantitative convergence or uncertainty information is reported. The text says convergence is checked by varying parameters, but no convergence data, error bars, or resolution estimates are given. Since the central quantitative claims rely on small energy shifts (for example, the tiny redshift in Figs. 4(b) and 8(b) and the red/blue shifts in Figs. 4(c,d) and 8(c,d)), the authors should provide a convergence test for the peak positions and estimate the numerical uncertainty of the reported shifts.","section":"Sec. II C and Figs. 4 and 8"}],"minor_comments":[{"comment":"There is a spurious space in 'ap proximation' in the title; the phrase 'moderate and a high intensity' is also awkward and could be rephrased.","section":"Title and abstract"},{"comment":"The differences between panels (a) and (b) of Fig. 1 are said to be hardly visible; a difference plot or a color-scale annotation would help the reader identify the nondipole changes before moving to the polar-coordinate figures.","section":"Sec. I, Fig. 1"},{"comment":"The solid and dashed curves are identified only in the main text; adding a legend or explicit line labels in the figure panels would improve readability.","section":"Sec. III, Fig. 4"},{"comment":"The symbol n·p is overloaded: it is first defined as p0 − p·n, then used inside a small correction in Eq. (33), where it introduces a large zeroth-order term. A separate symbol, such as q = p·n, would clarify the expansion.","section":"Sec. IV, Eq. (33)"},{"comment":"The paper does not report the total norm after propagation; it gives only the ground-state population and the escape probability. A norm-conservation check would be a useful validation of the numerical scheme.","section":"Sec. II C"}],"recommendation":"major_revision","confidential_remarks":"The numerical part of the paper appears credible and the central effect is clearly demonstrated in the computed distributions. The main blocker is the derivation of Eq. (36), which is internally inconsistent as printed. This is repairable, so I recommend major revision rather than rejection. The editor may also wish the authors to clarify the relationship between the nonrelativistic TDSE used numerically and the relativistic Volkov expansion used analytically, since the sign conventions in that expansion are not fully pinned down."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is worth a look: it solves the 2D TDSE for a hydrogen-like atom in a propagating laser pulse, compares with the dipole approximation, and clearly demonstrates that the directional tilt of multiphoton peaks is dominated by the recoil (Nordsieck) term, with retardation giving only a tiny redshift. The classification of 'moderate' vs 'strong' pulses via Ω_n vs ω/2 is simple and useful. The numerical work appears honest: the method is established, parameters are explicit, and the comparison to dipole results is directly visible.\n\nThe soft spot is the analytic derivation of Eq. (36), the formula that carries the main quantitative claim. As printed, Eq. (33) contains the recoil term (i/(m_e c)) Φ_p^(0)(t) n·p, with n·p = p0 - p·n. For nonrelativistic electrons, n·p ≈ m_e c + E_p/c - p2. Substituting that into Eq. (33) gives an extra full Φ^(0) phase and, more importantly, a p2 term with the opposite sign to what a correct leading-order expansion of Eq. (32) yields. The correct expansion gives a phase correction (Φ^(0)/(m_e c))(p2 - E_p/c), which indeed leads to Eq. (36) and the observed redshift for forward emission. But Eq. (36) does not follow from Eqs. (33)-(34) as printed; the derivation has a sign error that is load-bearing. This should be fixed and checked carefully.\n\nA second issue: in Sec. V, the strong-pulse comparison uses reference energies E^(p)_1,2 that are shifted by -0.3E0 from the dipole values, attributed to retardation, but the shift looks post hoc. A quantitative estimate of that redshift would strengthen the interpretation. Also, the paper gives no error bars or convergence measures beyond a statement that convergence was checked; that's minor for a numerics paper of this type, but a few numbers would help.\n\nThe 2D model limits direct transfer to experiment, but the qualitative picture is convincing and the classification criterion is a nice takeaway. The citations to the relevant nondipole literature are appropriate and the comparison is fair.\n\nRecommendation: send it to peer review. The numerical results are solid and the central physics claim is likely correct, but the derivation of Eq. (36) must be corrected and the strong-pulse overlay justified. After that, I'd be happy to cite it.","headline":"Solid numerical study of nondipole multiphoton ionization, but the analytic derivation of the key peak-shift formula has a sign error that needs fixing.","tokens_in":15026,"tokens_out":11073,"would_cite":true,"duration_ms":92676,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A propagating laser pulse's recoil effect, not field retardation, controls the directional shifts of multiphoton photoelectron peaks.","keywords":["nondipole effects","multiphoton ionization","Volkov state","recoil correction","retardation correction","photoelectron momentum distribution","strong-field ionization","propagating laser pulse"],"falsifier":"For a chosen multiphoton peak, compute the peak energy as a function of emission angle in the full 2D TDSE and test whether it follows $E_n(\\phi_p)=E_n-U_p\\sqrt{2E_n/(m_e c^2)}\\sin\\phi_p$; a deviation larger than the next-order $1/c^2$ contribution, or a persisting tilt when the Nordsieck term is artificially removed from the strong-field-approximation phase, would overturn the recoil-dominated reading.","tokens_in":13931,"feed_emoji":"⚛️","tokens_out":7793,"duration_ms":67939,"temperature":0.7,"pith_summary":"The paper asks what breaks the dipole picture of multiphoton ionization when the driving laser pulse propagates in space instead of being treated as a spatially uniform oscillating field. By solving the two-dimensional Schrödinger equation for a hydrogen-like atom driven by a flat-top pulse, the authors show that multiphoton energy peaks become angle-dependent: they shift red for some emission directions and blue for others. Using a leading-order relativistic expansion of the Volkov electron state, they trace this directionality to the electron recoiling off the laser pulse (the Nordsieck correction), while the field-retardation correction only produces a tiny redshift. If this is right, the distinct roles of the two nondipole effects are cleanly separated, and the size of the angle-dependent shift is fixed by a simple formula involving the ponderomotive energy and the photoelectron energy.","feed_headline":"Electron recoil, not retardation, tilts photoelectron spectra","feed_subtitle":"The emission angle decides whether multiphoton peaks shift red or blue, following a simple recoil formula.","key_machinery":"The carrying object is the leading-order relativistic expansion of the electron Volkov state, Eq. (33), in which the nonrelativistic Volkov phase is augmented by two 1/c terms. The first, pi_p(t) n·x, is the retardation correction coming from the space- and time-dependence of the field; the second, $Phi_p^{{(0)}}$(t) n·p/(m_e c), is called the recoil or Nordsieck correction. For photoelectrons emitted along the polarization axis the recoil term vanishes exactly, which isolates the retardation redshift; for other angles the recoil term produces the sinusoidal quasi-energy shift of Eq. (36). The numerical stripe tilt is matched by applying that formula with a constant-amplitude plane-wave assumption inside the flat-top pulse.","core_discovery":"In the dipole approximation the multiphoton peaks of the photoelectron energy-angular distribution appear as vertical stripes in polar coordinates: for a given peak, the energy does not depend on the electron emission angle. When the same pulse is treated as a wave propagating along the x2-direction, the stripes tilt, so the peak energy acquires a sinusoidal dependence on the emission angle. The paper's central result is an analytic formula for this dependence, E_n(phi_p) = E_n - U_p sqrt(2 m_e E_n)/(m_e c) sin phi_p, obtained from the leading 1/c expansion of the Volkov state (the exact electron state in a plane-wave field, ignoring the binding potential). Through this formula the tilt is attributed almost entirely to the electron recoil (Nordsieck) term in the Volkov phase, with the retardation term contributing only a small energy redshift that is most clearly seen for emission along the polarization axis. For stronger pulses the sinusoidal shift grows until it exceeds half the photon frequency, and the angle-integrated multiphoton peaks wash out even though they remain visible in the dipole approximation; the rescattering plateau also shrinks because the propagating pulse pushes the electron wave packet away from the ion.","pith_inferences":["If the paper's interpretation is right, angle-resolved photoelectron spectra could be fit directly to $E_n(\\phi_p)=E_n-\\Omega_n\\sin\\phi_p$; a residual angle dependence beyond this sinusoidal form would signal corrections beyond the leading $1/c$ order.","Because the same Volkov expansion is dimension-independent, one would expect analogous tilted stripes and the same recoil-versus-retardation split in three dimensions, with the emission angle measured relative to the propagation axis.","The moderate/strong boundary $\\Omega_n=\\omega/2$ offers a parameter-free diagnostic: once the predicted shift for a peak of interest exceeds half the photon energy, propagation corrections can no longer be treated as a small perturbation of the dipole spectrum.","The paper's statement that higher-order relativistic corrections are beyond its scope suggests a natural extension: repeating the analysis at x-ray photon energies, where the leading $1/c$ approximation is known to be insufficient, would show whether the recoil term continues to dominate the directionality or whether higher-order terms alter the sinusoidal shift."],"forward_implications":["For a given multiphoton line, emission directions with a positive component along the pulse propagation direction should show a redshift, and opposite directions a blueshift, with amplitude $U_p\\sqrt{2E_n}/(m_e c)$.","When the shift amplitude $\\Omega_n$ exceeds $\\omega/2$, the angle-integrated energy distribution loses resolved multiphoton structure; this serves as the paper's quantitative boundary between moderate and strong pulses.","The retardation correction is most visible for electrons emitted along the polarization axis and shows up as a small redshift of the peak, consistent with earlier x-ray pulse studies.","The double-hump shape seen on intermediate-energy peaks combines the dipole-approximation sidelobes with the tilted-stripe geometry of the propagating pulse.","At higher intensity the rescattering plateau shrinks and eventually fades, because the propagating pulse separates the electron wave packet from the parent ion before in-pulse rescattering can occur."],"supporting_citations":[{"why":"Supplies the Suzuki-Trotter split-step numerical scheme and the 2D propagating-pulse treatment on which the TDSE results rest.","marker":"[8]"},{"why":"Provides the mask-and-Fourier method used to extract photoelectron momentum distributions from the final wave packet.","marker":"[22, 23]"},{"why":"Defines the Volkov state used as the zeroth-order strong-field-approximation electron state.","marker":"[28]"},{"why":"Supplies the relativistic Volkov solution of the Klein-Gordon equation that is the starting point for the 1/c expansion.","marker":"[29]"},{"why":"Gives the leading-order 1/c form of the Volkov phase with the retardation and recoil terms and its prior use in radiative recombination.","marker":"[30]"},{"why":"Names and identifies the Nordsieck (recoil) correction that the paper finds dominates the directionality dependence.","marker":"[31]"},{"why":"Earlier x-ray-pulse studies whose redshift conclusion the paper confirms and classifies as a retardation effect.","marker":"[34, 35]"}],"fun_headline_variants":["Recoil, not retardation, tilts photoelectron energy peaks","Angle-dependent photoelectron energies traced to electron recoil","Propagating pulse reveals recoil-driven tilt in multiphoton peaks","Recoil term sets sinusoidal angular shift in photoionization spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The interpretation assumes that a finite flat-top laser pulse behaves locally like a constant-amplitude plane wave, so that the angle-dependent peak shift computed for plane waves applies to the real pulse; if pulse fronts, edges, or the two-dimensional geometry break that picture, the separation of recoil and retardation effects would have to be redone.","fun_headline_variants_meta":{"raw":{"variants":["Recoil, not retardation, tilts photoelectron energy peaks","Angle-dependent photoelectron energies traced to electron recoil","Propagating pulse reveals recoil-driven tilt in multiphoton peaks","Recoil term sets sinusoidal angular shift in photoionization spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000305,"raw_usage":{"total_tokens":1755,"prompt_tokens":957,"completion_tokens":798,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":573,"completion_tokens_details":{"reasoning_tokens":727}},"tokens_in":573,"tokens_out":798,"duration_ms":7355,"temperature":1.0,"reasoning_tokens":727,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:09:34.843407+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a chosen multiphoton peak, compute the peak energy as a function of emission angle in the full 2D TDSE and test whether it follows $E_n(\\phi_p)=E_n-U_p\\sqrt{2E_n/(m_e c^2)}\\sin\\phi_p$; a deviation larger than the next-order $1/c^2$ contribution, or a persisting tilt when the Nordsieck term is artificially removed from the strong-field-approximation phase, would overturn the recoil-dominated reading.","supporting_citations":[{"cited_title":"Ehlotzky, A","cited_arxiv_id":null,"evidence_quote":"Defines the Volkov state used as the zeroth-order strong-field-approximation electron state."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the relativistic Volkov solution of the Klein-Gordon equation that is the starting point for the 1/c expansion."},{"cited_title":"Kanti, M","cited_arxiv_id":null,"evidence_quote":"Gives the leading-order 1/c form of the Volkov phase with the retardation and recoil terms and its prior use in radiative recombination."}],"review_version":1}