{"id":"9f24e5fb-5542-467e-9efd-31843d7573b0","arxiv_id":"2507.04096","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Field-induced orbital distortion and higher parabolic channels, not excited states, can reverse the angular pattern of molecular tunneling ionization, turning a minimum into a maximum.","lead":"This paper shows that as the ionizing laser field gets stronger, two effects, distortion of the orbital shape and extra tunnelling channels, can change how molecular ionization depends on orientation, even turning a minimum into a maximum. The authors demonstrate the effect in CH3Br and introduce a faster way to compute orientation-dependent ionization rates.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"CH3Br min-to-max result rests on Gaussian-basis asymptotic tails: the basis set is unspecified, no basis/Lmax convergence is shown, and the predicted field range is ~4x below TDCI; an inaccurate tail could change the channel ratio at beta=0 and make the transformation a model artifact.","rationale":"The reader's conditional verdict and weakest assumption align with my read. The paper's own limitation statement in Sec. III B admits the Gaussian basis tail is imperfect and WFAT magnitudes are unreliable at a given field, and the CH3Br crossover appears at F~0.011 rather than the TDCI value F~0.05. Because the min-to-max crossover is controlled by the relative weights of the (0,0) and (0,+-1) channels at beta=0, a tail error can shift or erase the effect. The independent TDCI panel confirms the phenomenon at higher fields, so I do not regard the claim as false; rather, the mechanism and the quantitative field scale are not yet pinned down. The manuscript also does not specify the CH3Br basis or show Lmax/basis convergence for this molecule, which weakens the reproducibility of the headline result. A focused recomputation with augmented diffuse functions and higher Lmax would resolve whether the concern lands. Since the reader already conditioned the verdict on this type of issue, I recommend no verdict change.","tokens_in":30441,"tokens_out":20831,"duration_ms":233133,"concrete_test":"Recompute the CH3Br OE-WFAT(1) rates at F=0.0105, 0.0110, 0.0115, 0.0125 with: (i) an augmented diffuse basis (e.g., aug-cc-pvtz, and aug-cc-pvqz if feasible), and (ii) Lmax=20 and 25 plus the explicit-integration method for beta=0,10,20,30 (gamma=90 deg). Determine whether the normalized rate at beta=0 changes from a local minimum to a local maximum over the same field interval, and whether the crossover field shifts by less than a factor of ~2. If the transformation disappears, or if Lmax=15 differs materially from Lmax=20/25, the central mechanism claim fails; if it persists, the Gaussian-tail concern is mitigated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most fragile link in the central CH3Br claim is the Gaussian-basis asymptotic tail, and the manuscript does not provide the convergence evidence needed to rule out a tail-driven artifact. End of Sec. III B concedes that 'the limited accuracy of the ionization potential and the imperfect asymptotic tail modeled by Gaussian basis' makes WFAT rates 'less reliable for predicting the magnitude at any particular field strength.' For CH3Br the basis set is not even stated; for CO/OCS cc-pvtz (no diffuse augmentation) is used, and if the same is used for CH3Br the tail beyond ~3-5 bohr is dominated by Gaussians rather than the physical exp(-kappa r), with direction-dependent errors from the C and Br centers. The min-to-max transformation at beta=0 occurs when the (0,+-1) parabolic channels overtake the near-zero (0,0) channel; that crossover field is set by the ratio |g_{0,+-1}|^2/|g_{00}|^2, which is precisely a tail-sensitive quantity. The factor-of-~4 shift between the OE-WFAT(1) crossover (F~0.0105-0.0125) and the TDCI crossover (F~0.045-0.055) shows the tail error is large. Independent TDCI validates the phenomenon, but not the orbital-distortion/parabolic-channel decomposition, so the mechanism claim is not independently confirmed.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a partial-wave reformulation of the one-electron weak-field asymptotic theory including first-order corrections (OE-WFAT(1)) in the integral representation. It validates the reformulation on noble gases and H2+, then applies it to field-dependent orientation-resolved ionization rates of CO, OCS, and CH3Br. The central physical claim is that two effects independent of excited states—first-order orbital distortion and contributions from parabolic channels (0,±1)—can change the orientation dependence of tunneling ionization as the field increases. In CO the global maximum shifts from the O-to-C to the C-to-O orientation; in CH3Br a local minimum at beta=0 (field along C-Br) turns into a local maximum at higher fields. The authors also demonstrate a large speed-up of the orientation scan via the partial-wave method (up to about 1.5 x 10^4 for CH3Br).","tokens_in":30754,"tokens_out":5329,"duration_ms":54248,"significance":"The claim, if established, challenges the standard orbital-image intuition that the orientation dependence of tunnel ionization maps the field-free HOMO shape at all but the lowest fields, and it provides a mechanism that does not require excited-state participation. The paper's strengths include a detailed analytical derivation of the partial-wave decomposition of the first-order parabolic-channel function, validation of atomic parameters against the tail-representation values, qualitative agreement with RT-TDDFT for CO and OCS, a public dataset, and an NWChem implementation. The main risk is that the CH3Br min-to-max transformation is governed by ratios of asymptotic Gaussian-orbital tails, whose accuracy is explicitly conceded to be limited, and the H2+ validation leaves an unresolved discrepancy with the tail-representation reference.","major_comments":[{"comment":"The CH3Br calculations do not state the basis set used to generate the HOMO; for CO and OCS the text specifies cc-pvtz without diffuse augmentation. The central min-to-max transformation at beta=0 is decided by the field at which the (0,±1) parabolic channels overtake the (0,0) channel, i.e. by the ratio |g_{0,±1}|^2 / |g_{00}|^2, which is an asymptotic-tail quantity. The manuscript itself states (end of Sec. III B) that Gaussian-basis WFAT rates are 'less reliable for predicting the magnitude at any particular field strength,' and the OE-WFAT(1) crossover field in Fig. 7(a) (F about 0.0105-0.0125) is about a factor of four below the TDCI crossover (F about 0.045-0.055). Please specify the CH3Br basis, provide basis-set and Lmax convergence tests, and show that the crossover and the channel-ratio mechanism are stable under augmentation of the diffuse tail.","section":"III C, Fig. 7"},{"comment":"The H2+ results for both the 1s-sigma and 2p-pi+ states show noticeable deviations from the tail-representation results of Ref. [34], and the authors attribute this to the radial basis (FEDVR versus Laguerre) and to IR versus TR differences. Since H2+ is an exact one-electron system, both representations should converge to the same asymptotic rates in the complete-basis limit; an unresolved discrepancy here weakens the validation of exactly the first-order correction that drives the molecular conclusions. Please quantify the size of the discrepancy, test convergence with box size and radial-element parameters, and either resolve the difference or state its expected impact on the CO, OCS, and CH3Br rates.","section":"III A, Fig. 4"},{"comment":"The statement that Gaussian-basis WFAT rates cannot reliably predict field magnitudes is in tension with the paper's use of a specific field-dependent crossover as the central result. The paper needs an explicit sensitivity analysis—for example, rerunning CH3Br with a diffuse-augmented basis, with a range of tuned functionals, or with an exactly-tailed model HOMO—to show that the min-to-max transformation and its mechanism are not artifacts of the finite Gaussian tail. Without such a test, the claim that the mechanism is independent of excited-state effects is not fully separated from numerical artifacts.","section":"III B, last paragraph; III C"}],"minor_comments":[{"comment":"There is a typo in the Introduction: 'excited the state contribution' should be 'excited-state contribution'.","section":"I, Abstract"},{"comment":"The mixed-order approximation in Eq. (18), with Gamma_00 taken to first order and Gamma_{0,±1} taken to zeroth order, should be justified explicitly in the text rather than only by referencing the order of F in W_nu.","section":"II B, Eq. (18)"},{"comment":"The left and right halves of Fig. 7(a) use different field ranges (TDCI: 0.045-0.055; OE-WFAT(1): 0.0105-0.0125); the caption should state this explicitly so readers do not interpret the two panels as directly comparable field strengths.","section":"Fig. 7(a)"},{"comment":"The dataset DOI in Ref. [71] is welcome; please state in the main text that the data and code underlying the figures are available, rather than only in a reference.","section":"IV, Data availability"},{"comment":"The comparison with Ref. [36] reports a00 and B00, but A00 is also listed in the table; please state in the caption or text whether A00 was compared and with what agreement.","section":"Table I"}],"recommendation":"major_revision","confidential_remarks":"The unresolved H2+ discrepancy and the Gaussian-tail sensitivity of the CH3Br crossover are the key risks. The manuscript is otherwise well within scope and the partial-wave reformulation is a useful methodological contribution. I would ask for the basis-set and convergence analysis for CH3Br and a sensitivity study before considering acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Candid take: the paper's headliner is the partial-wave expansion of the first-order correction h_nu (Eqs. 25-28), which is a real technical contribution, and the CH3Br minimum-to-maximum transformation is a new physical claim, not just a fit to data. Validation on noble gases matches the TR reference well; CO and OCS compare qualitatively well with RT-TDDFT; the internal decomposition via Eq. (18) isolates parabolic channel contributions cleanly. No parameter fitting to the target yields—the screening parameters are standard inputs.\n\nSoft spots, in order: (1) The H2+ comparison with Ref. [34] leaves an unresolved discrepancy in normalized rates, and the suggested causes (basis vs. representation) are not tested. (2) For CH3Br, the basis set is not stated, no basis-set convergence is shown, and Lmax=15 is asserted without a convergence check. Given the paper's own concession that Gaussian asymptotic tails make rate magnitudes 'less reliable,' the exact crossover field in Fig. 7 is model-dependent. (3) The CH3Br field range is about a factor of four lower than the TDCI comparison, which the paper notes but does not analyze. These caveats weaken the quantitative field scale, but not the qualitative mechanism: Fig. 7(b) directly demonstrates the channel crossover within the model, and TDCI independently shows the same min-to-max trend.\n\nThe stress-test worry about asymptotic tails is real but lands on the field scale, not on the mechanism. I would want the authors to state the CH3Br basis, add a convergence check, and discuss the field-ratio difference explicitly.\n\nThis paper is for strong-field and molecular-imaging people who rely on the simple orbital-shape mapping at higher intensities. They should read it. I'd send it to peer review; a referee should ask for the missing basis/convergence details and a more careful comparison of field ranges.","headline":"Real technical contribution in the partial-wave reformulation, a credible min-to-max mechanism in CH3Br, but the quantitative field scale is model-dependent and needs convergence evidence.","tokens_in":31324,"tokens_out":2582,"would_cite":true,"duration_ms":29073,"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":"This paper claims that field-induced orbital distortion and higher parabolic channels, not excited-state effects, can reverse the orientation dependence of molecular tunneling ionization, turning a minimum in CH3Br into a maximum.","keywords":["tunneling ionization","weak-field asymptotic theory","orbital distortion","parabolic channels","orientation dependence","CH3Br","partial-wave expansion","strong-field ionization"],"falsifier":"Recompute the CH3Br orientation-dependent rates at F=0.0105 to 0.0125 a.u. with an orbital whose large-distance tail reproduces the true Coulomb decay, for example a grid or exponential basis fitted to the same ionization potential; if the beta=0 rate remains below its neighbors at every field instead of rising into a maximum, the minimum-to-maximum transformation is an artifact of the Gaussian asymptotic tail.","tokens_in":1701,"feed_emoji":"⚛️","tokens_out":1617,"duration_ms":93515,"temperature":0.7,"pith_summary":"This paper argues that two field-induced effects—distortion of the ionized orbital and the contribution of higher parabolic channels—can reshape the angular dependence of molecular tunneling ionization as the field grows, even where excited-state effects are negligible. The central demonstration is CH3Br: as the ionizing field increases, a local minimum in the orientation-dependent rate at the angle where the field points along the C–Br axis turns into a local maximum. In CO, orbital distortion alone moves the global maximum to the opposite orientation, bringing the theory into qualitative agreement with experiments that zeroth-order theory misses. The claim is carried by the one-electron weak-field asymptotic theory with first-order corrections (OE-WFAT(1)) and by agreement with independent time-dependent calculations. If the mechanism is right, angular ionization yields at high field cannot be read as a simple image of the field-free HOMO.","feed_headline":"CH3Br ionization minimum flips to a maximum as the field grows","feed_subtitle":"Field-induced orbital distortion plus extra parabolic channels, not excited states, drive the flip.","key_machinery":"The central object is OE-WFAT(1) in the integral representation, a tunneling theory in which the outer electron's wave function is expanded as psi approximately psi^(0) + F psi^(1) and the partial rate to each final parabolic channel nu=(n_xi,m) is kept to first order in the field. The paper's new mathematical result is a partial-wave reformulation: the distorted final-state functions $\\Omega$^(0)_nu and $\\Omega$^(1)_nu are expanded in spherical harmonics so that the expensive orientation-independent integrals are computed once and stored, and only cheap angular sums (Eqs. 25–28) are done per orientation. This makes full scans over Euler angles ($\\beta$,gamma) affordable and isolates the two physical effects being studied.","core_discovery":"At zeroth order, OE-WFAT(0), the field and orientation dependences of the rate factorize, so the angular shape of the ionization yield is field-independent. The first-order theory couples them through the field-induced orbital distortion and through the parabolic channels (0,±1). The paper's central claim is that in CH3Br the (0,0) channel's rate near beta=0 is suppressed by orbital distortion as the field grows, and the (0,±1) channels then take over and turn what was a local minimum into a local maximum (Eq. 18 and Fig. 7). In CO and OCS the same machinery identifies orbital distortion as the dominant field-dependent effect, with parabolic-channel contributions minor; in all three molecules the changes are shown not to come from HOMO-1 (excited-state) ionization, whose angular pattern peaks at the opposite orientation.","pith_inferences":["Generalizing from CH3Br, any molecule whose nodal structure produces a rate minimum aligned with the field could show a similar minimum-to-maximum reversal once its (0,±1) channels become competitive; the CO and OCS cases suggest the ingredients are generic even though the reversal is not observed in every molecule.","The range of field strengths over which the reversal occurs is basis-sensitive in the paper's own account, so an experimental test at the TDCI field window (around F=0.05 a.u.) and at the OE-WFAT(1) window (around F=0.011 a.u.) would distinguish the physical mechanism from asymptotic-tail artifacts.","The computational speed-up makes it reasonable to use distortion-corrected rates in molecular orbital tomography and laser-induced electron diffraction inversions; those analyses currently assume a field-free HOMO shape, and correcting for distortion could change reconstructed orbitals at higher intensities."],"forward_implications":["In CO, first-order orbital distortion is necessary to reproduce the experimentally observed shift of the global rate maximum to the C-to-O orientation as the field increases; the paper reports qualitative agreement with experiment and with RT-TDDFT.","In OCS, orbital distortion suppresses the rate around beta=0 degrees and leaves the peak positions unchanged, while contributions from the nu=(0,±1) parabolic channels remain negligible.","In CH3Br, omitting the (0,±1) channels leaves the beta=0 feature a minimum at all studied fields, so the minimum-to-maximum conversion requires both orbital distortion and parabolic-channel contributions.","OE-WFAT(1) in partial-wave form reproduces TR reference rates for noble gases and H2+ and gives speed-ups that grow linearly with the number of orientation angles, reaching about 1.5e4 for 18,029 orientations.","The first-order correction couples field and orientation dependences, so OE-WFAT(0)-based orbit-to-yield mapping fails as intensity increases; techniques that assume the HOMO-shape mapping become unreliable at higher fields."],"supporting_citations":[{"why":"Supplies the OE-WFAT(1) integral representation and the central rate formula that the paper builds on.","marker":"[21]"},{"why":"Defines the first-order correction terms in WFAT and provides the noble-gas reference values used for validation.","marker":"[36]"},{"why":"Applies first-order-corrected OE-WFAT to molecules and provides the H2+ reference rates used as a benchmark.","marker":"[34]"},{"why":"Gives the Becke-cell grid integration methodology underlying the explicit method that the partial-wave reformulation speeds up.","marker":"[38]"},{"why":"Provides independent haCC field-scan results for CO that the paper's OE-WFAT(1) rates reproduce.","marker":"[57]"},{"why":"Supplies the TDCI angular yields of CH3Br whose minimum-to-maximum behavior is reproduced and explained.","marker":"[60]"},{"why":"Provides the experimental CO orientation-dependent ionization data used as the qualitative target.","marker":"[55]"},{"why":"Established field-induced orbital distortion in aligned molecules, which this paper extends with a field scan.","marker":"[29]"}],"fun_headline_variants":["CH3Br ionization dip flips to peak at higher field","Field flips CH3Br ionization minimum into maximum","Orbital distortion turns CH3Br ionization minimum into peak","Orbital distortion, not excited states, flips CH3Br ionization minimum"],"cache_read_input_tokens":33408,"weakest_assumption_plain":"The load-bearing premise is that the field-free HOMO obtained from a tuned range-separated DFT calculation with a finite Gaussian basis has an asymptotic tail accurate enough for the orientation-dependent rates; the paper itself cautions that this tail limits the reliability of rate magnitudes at a particular field strength, so the field-dependent switch could in principle be an artifact of the model.","fun_headline_variants_meta":{"raw":{"variants":["CH3Br ionization dip flips to peak at higher field","Field flips CH3Br ionization minimum into maximum","Orbital distortion turns CH3Br ionization minimum into peak","Orbital distortion, not excited states, flips CH3Br ionization minimum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001326,"raw_usage":{"total_tokens":5395,"prompt_tokens":940,"completion_tokens":4455,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":556,"completion_tokens_details":{"reasoning_tokens":4382}},"tokens_in":556,"tokens_out":4455,"duration_ms":36317,"temperature":1.0,"reasoning_tokens":4382,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T19:56:28.718594+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the CH3Br orientation-dependent rates at F=0.0105 to 0.0125 a.u. with an orbital whose large-distance tail reproduces the true Coulomb decay, for example a grid or exponential basis fitted to the same ionization potential; if the beta=0 rate remains below its neighbors at every field instead of rising into a maximum, the minimum-to-maximum transformation is an artifact of the Gaussian asymptotic tail.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the OE-WFAT(1) integral representation and the central rate formula that the paper builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the first-order correction terms in WFAT and provides the noble-gas reference values used for validation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Applies first-order-corrected OE-WFAT to molecules and provides the H2+ reference rates used as a benchmark."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Becke-cell grid integration methodology underlying the explicit method that the partial-wave reformulation speeds up."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides independent haCC field-scan results for CO that the paper's OE-WFAT(1) rates reproduce."},{"cited_title":"Hoerner and H","cited_arxiv_id":null,"evidence_quote":"Supplies the TDCI angular yields of CH3Br whose minimum-to-maximum behavior is reproduced and explained."},{"cited_title":"Zhang, J","cited_arxiv_id":null,"evidence_quote":"Provides the experimental CO orientation-dependent ionization data used as the qualitative target."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Established field-induced orbital distortion in aligned molecules, which this paper extends with a field scan."}],"review_version":1}