{"id":"7b325b7a-fb06-497d-bb8e-e756bd80ee98","arxiv_id":"2502.08811","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A parity-based separation of two-photon ionization amplitudes lets one continuous RABBITT simulation resolve whole series of resonant states in He and H2, and tests the logarithmic Hilbert transform linking magnitude and phase.","lead":"This paper introduces a computational trick, crRABBITT, that extracts attosecond interferometry parameters continuously across the photoelectron spectrum in a single simulation run. It uses this trick to map autoionizing states in helium and molecular hydrogen and to test a mathematical relation between signal magnitude and phase.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Parity separation extracts even-photon, not strictly two-photon, amplitudes; uncontrolled higher-order contamination could bias continuous RABBITT parameters.","rationale":"The reader's weakest_assumption correctly identifies the parity-based separation as the load-bearing premise of crRABBITT. My analysis agrees that this is the most central uncontrolled approximation, but refines the reader's phrasing: the symmetrization does not isolate the pure two-photon amplitude; it isolates the even-parity component, which also contains four-photon and potentially two-XUV-photon contributions. The two-XUV-photon term is independent of the IR delay and therefore only affects the constant A in the sideband oscillation, not B and C. However, the four-photon XUV+3IR term can interfere with the two-photon term at the 2ωτ frequency and thus contaminate B and C. The paper's intensities (IR 1e10 W/cm^2) suggest this contamination is very small, but no numerical bound or comparison with the direct subtraction method is provided. This is a concrete, checkable gap rather than a demonstrated failure. The rest of the central claim - continuous extraction across resonances and the LHT application - appears internally consistent, and the comparisons with literature support the method's accuracy. The absence of code and data, already noted by the reader, reinforces the need for a reproducibility check. I therefore recommend keeping the CONDITIONAL verdict, with an added condition that the parity-separated amplitude be validated against direct subtraction or an IR-intensity scan.","tokens_in":9078,"tokens_out":16593,"duration_ms":155468,"concrete_test":"Run the He sp2+ case with the present intensities and extract the sideband amplitude by (i) parity symmetrization and (ii) direct subtraction An0(k)-An_XUV(k) using a separate XUV-only TDSE run. If the two amplitudes agree to within the numerical accuracy of the reported B and C, the even-photon contamination is negligible. Additionally, repeat at IR intensities of 1e11 and 1e12 W/cm^2; the B and C parameters should remain unchanged within the stated uncertainty if the 2ωτ contamination is negligible.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"In Sec. II, the paper states that 'the pure two-photon ionization amplitude can be extracted by a simple symmetrization An(k)=[An0(k)+Pn An0(-k)]/2'. Strictly, this operation selects the even-parity component of the single-pulse amplitude, i.e., all even-photon orders: the two-photon XUV+IR term plus the four-photon XUV+3IR term (and any two-XUV-photon term). The one-photon XUV contribution (odd parity) is removed, but even-order terms are retained. The RABBITT B and C parameters are extracted from the 2ωτ Fourier component of the sideband signal; the four-photon term can interfere with the two-photon term to produce a 2ωτ contribution that is not distinguishable from the desired two-photon signal. The paper asserts the XUV and IR fields are weak, but gives no quantitative estimate of the relative magnitude of the four-photon amplitude. For the IR intensity 1e10 W/cm^2 used here, the four-photon amplitude scales as E_IR^3 and is likely small, but this is not demonstrated. The absence of a comparison with the direct subtraction method An0(k)-An_XUV(k), which the same section presents as the general alternative, leaves the purity of the 'two-photon' amplitude unverified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a computational variant of rainbow RABBITT, called continuous rainbow RABBITT (crRABBITT), in which the two-photon XUV+IR ionization amplitude is isolated from a single TDSE run by exploiting parity: for a symmetric target, symmetrizing the ionization amplitude over photoelectron momentum removes the odd-parity one-photon contribution and leaves an even-parity amplitude that the authors identify with the two-photon process. From this amplitude the RABBITT magnitude B and phase C are obtained continuously across the photoelectron spectrum, allowing the authors to resolve whole series of autoionizing states in He and H2 as well as below-threshold 1s→np excitations in He without scanning the laser frequency. The extracted amplitude and phase are then used to test the logarithmic Hilbert transform (LHT) as a relation between the two-photon RABBITT magnitude and phase, which the authors claim as the first successful application of LHT to two-photon ionization.","tokens_in":9364,"tokens_out":5745,"duration_ms":57325,"significance":"If the parity-based separation were fully validated, crRABBITT would be a useful computational shortcut: one TDSE run would map out resonant structure over a broad spectral range that currently requires multiple frequency-tuned RABBITT calculations. The paper's TDSE results are internally consistent, and the Fano parameters in Table I agree reasonably with the cited literature for both He and H2, which gives some confidence in the underlying numerical method. The LHT test for two-photon parameters, if established, would be of methodological interest for attosecond interferometry. However, the central methodological step — the identification of the symmetrized amplitude with the pure two-photon amplitude — is not quantitatively verified, and the numerical convergence of the TDSE results is not documented, so the strength of the paper currently rests on an unexamined assumption.","major_comments":[{"comment":"The parity symmetrization An(k) = [An0(k) + Pn An0(-k)]/2 removes the odd-parity one-photon amplitude but retains all even-photon orders, not just the XUV+IR two-photon amplitude. In particular, an XUV+3IR (four-photon) amplitude has the same even parity and can interfere with the two-photon term to produce a contribution at the same 2ωτ Fourier component of the sideband signal. The paper states that the XUV and IR fields are weak and quotes the IR intensity as 1×10^10 W/cm^2, but it provides no quantitative estimate of the relative magnitude of the four-photon amplitude, nor does it compare the symmetrization result with the direct subtraction An0 − An_XUV that is described in the same paragraph. Without such a comparison, the extracted B and C parameters are not demonstrated to be pure two-photon RABBITT parameters.","section":"Sec. II, after Eq. (8)"},{"comment":"No convergence tests or numerical uncertainties are reported for the TDSE calculations. The quantitative outputs include Fano resonance widths as narrow as 8 meV in He and continuous phase C across resonances, so the resolution of the calculation matters directly for the central claims. The paper should specify the numerical parameters (radial grid and finite-element sizes, angular momentum truncation, number of configurations Ns, time step, ECS parameters, and APT pulse durations) and provide tests showing that the extracted B and C, and hence the Fano parameters in Table I, are converged with respect to these choices. Without this information the claimed 'very fine energy resolution' cannot be assessed.","section":"Sec. III A, Figs. 1–3 and Table I"},{"comment":"The manuscript never states explicitly how the RABBITT parameters B and C of Eq. (1) are obtained from the extracted two-photon amplitude An(k). The reader is left to infer whether B = |An(k)| and C = arg An(k) for some fixed delay, or whether a delay scan is performed, and how the continuous spectra in Figs. 2–4 were constructed. The extraction formulas and the number of delays (if any) need to be stated so that the results are reproducible.","section":"Sec. II, Eqs. (7)–(8) and Figs. 2–4"},{"comment":"The evidence for the claimed 'first successful application of the LHT to two-photon ionization' is limited. In the above-threshold case, the LHT comparison is shown only for the He sp2+ resonance shifted by ±ω, where a single resonance is embedded in one RABBITT arm; the paper itself acknowledges that the relation would be difficult to establish when two or more resonances intertwine. In the below-threshold case, Fig. 4 shows only 'qualitatively similar' agreement with experiment. To make the LHT claim convincing, the authors should provide a quantitative measure of the difference between the LHT-derived phase and the direct TDSE phase in every shown case, and should specify the conditions under which the relation is expected to hold.","section":"Sec. III B, Eq. (11) and Figs. 2–4"}],"minor_comments":[{"comment":"The word 'Menawhile' should be 'Meanwhile'.","section":"Introduction, paragraph 2"},{"comment":"The word 'exciations' should be 'excitations'.","section":"Sec. III B, title"},{"comment":"The caption reads 'B ans C'; it should read 'B and C'.","section":"Fig. 2 caption"},{"comment":"The phrase 'resolve tje whole series' contains a typo; it should be 'resolve the whole series'.","section":"Conclusion"},{"comment":"The notation ⌊NAPT/2⌋ is used without specifying that NAPT is odd; for NAPT = 41 the floor is unambiguous but a brief statement would help the reader.","section":"Sec. II, Eq. (4)"}],"recommendation":"major_revision","confidential_remarks":"The main methodological concern — that parity separation retains all even-photon orders and may contaminate the two-photon RABBITT parameters — is central to the paper and should be addressed with a quantitative estimate or a direct-subtraction benchmark before publication. In addition, the absence of convergence tests is a significant omission for a numerical methods paper. The authors should be asked to include a data-availability statement if possible, as the community would benefit from the TDSE outputs underlying Figs. 2–4."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper does something genuinely new: it shows that for atomic or symmetric-molecular targets, a single TDSE run with XUV+IR can give continuous RABBITT B and C parameters across the whole photoelectron spectrum, by symmetrizing the amplitude to isolate the even-parity part. That replaces the usual frequency-scanning rRABBITT workflow, and the demonstration on He and H2 is convincing. The Fano parameters extracted from their XUV-only spectra match the literature values, and the phase from the logarithmic Hilbert transform tracks the direct TDSE phase, with a qualitative match to the Neoricic experiment. The LHT application to two-photon parameters is a useful consistency test, even if it's on their own simulation output rather than new experimental data.\n\nThe main soft spot is the absence of numerical hygiene: no convergence tests with respect to grid parameters, basis size, or pulse durations, and no uncertainty estimates on the quoted Fano parameters. For a method paper this matters. There's also no code or data deposit, so the reader can't reproduce the numbers without reimplementing the whole TDSE machinery.\n\nOn the stress-test concern about parity separation: the symmetrization does select even-photon orders, not strictly two-photon. But at their IR intensity of 1e10 W/cm^2 (E ~ 5e-4 a.u.), the XUV+3IR amplitude is smaller than the XUV+IR amplitude by a factor on the order of E^2 ~ 3e-7. That is not a plausible source of error in these results. The paper could have said this in one line, but the omission is minor, not a load-bearing flaw. The direct-subtraction comparison would be nice, but the agreement with literature Fano parameters and the LHT check already provide independent validation.\n\nThe paper is short, clearly written, and the central claim holds up. I'd send it to review. The main request to the authors would be to add convergence tests and scatter some error bars, and to state explicitly why higher-order IR terms are negligible. If they do that, it's a solid contribution to the attosecond metrology toolbox.","headline":"Parity-separated crRABBITT is a real computational shortcut that works; the paper's only real gap is missing convergence detail, and the higher-order contamination worry is negligible at their intensities.","tokens_in":9879,"tokens_out":2971,"would_cite":true,"duration_ms":27782,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["32.80.Rm","32.80.Fb","42.50.Hz"],"model":"deepseek-v4-flash","headline":"A parity-based separation of the ionization amplitude recovers the pure two-photon part and maps entire resonant-series spectra in one TDSE run, with the logarithmic Hilbert transform linking phase and magnitude.","keywords":["continuous rainbow RABBITT","parity-based amplitude separation","attosecond interferometry","autoionizing states","under-threshold excitations","logarithmic Hilbert transform","time-dependent Schrödinger equation","Fano parameters"],"falsifier":"Take a target with no definite parity, or raise the infrared intensity until higher-order sidebands appear, and compare the parity-separated B and C against a full frequency-tuned rRABBITT scan in the same energy window; a mismatch would show the symmetrized amplitude is not purely two-photon.","tokens_in":8894,"feed_emoji":"⚛️","tokens_out":6156,"duration_ms":53573,"temperature":0.7,"pith_summary":"This paper establishes a computational scheme, continuous rainbow RABBITT (crRABBITT), that turns the usual RABBITT sideband analysis into a continuous function of photoelectron energy. Because the target has definite parity, the one-photon ionization amplitude is odd and the two-photon amplitude is even, so a single symmetrization of the amplitude separates the two-photon component from the background without tuning the dressing-laser frequency. That separation yields RABBITT magnitude and phase across the whole spectrum at once, resolving series of autoionizing states above threshold in He and H2 and discrete 1s→np excitations below threshold in He, all from one time-dependent Schrödinger equation run. The same extracted parameters are then used to check a logarithmic Hilbert transform connecting RABBITT phase and magnitude, which the authors report as the first successful application of this transform to two-photon ionization.","feed_headline":"One parity trick maps whole resonance series in one TDSE run","feed_subtitle":"Continuous rainbow RABBITT resolves autoionizing and below-threshold states in a single shot and ties phase to magnitude.","key_machinery":"The load-bearing object is the parity-symmetrized ionization amplitude An(k)=[An0(k)+Pn An0(−k)]/2, where An0(k) is the ionization amplitude of the ground-state target in the combined XUV+IR field and Pn is the parity of the residual ionic state n. Since one-photon absorption changes parity while two-photon absorption preserves it, this formula isolates the pure XUV+IR two-photon amplitude and is what allows B and C to be read off continuously. The paper combines it with a multiconfiguration two-electron TDSE solver on a discrete-variable-representation grid with exterior complex scaling and t-SURFFc amplitude extraction, and with Fano-profile fits for the resonant lineshapes, to produce the spectra and to test the logarithmic Hilbert transform C(E)=−(1/π)P∫ dx lnB(x)/(x−E).","core_discovery":"In the paper's own terms, the central discovery is that the parity of the target can be exploited to isolate the two-photon XUV+IR ionization amplitude exactly: with a positive-parity initial state, the one-photon amplitude is negative parity and the two-photon amplitude positive parity, so the combination An(k)=[An0(k)+Pn An0(−k)]/2 removes the one-photon background. This isolation makes the RABBITT parameters B and C continuous functions of photoelectron energy, resolving the entire series of resonant features in a single TDSE simulation and eliminating the need to repeat calculations at different IR photon frequencies. Applied to He and H2, it reproduces Fano parameters of autoionizing states, maps under-threshold np excitations, and provides the data needed to show that the logarithmic Hilbert transform relates the resonant RABBITT phase to the magnitude, extending the Kramers-Kronig program from single-photon to two-photon ionization.","pith_inferences":["A direct experimental analogue would need to separate XUV-only and XUV+IR wavepacket components by their distinct angular distributions, as the paper notes; a testable design would verify whether the resulting continuous B and C match the parity-separated predictions.","If the logarithmic Hilbert transform holds for this two-photon case, it may also apply to circularly polarized RABBITT observables, whose under-threshold phase-magnitude relation has not yet been examined.","The requirement that only one ionic orbital be continuous suggests the method could scale to heavier atoms and molecules by enlarging the multiconfiguration basis, though that scalability remains an extrapolation rather than a demonstrated result.","Comparing crRABBITT results against conventional frequency-tuned rRABBITT in the same energy window would provide a quantitative cross-check of the parity-separation assumption."],"forward_implications":["All sidebands in a photoelectron spectrum can be analyzed at once rather than one by one, so resonant structure is mapped across a wide energy window in a single run.","Below-threshold discrete-state series appear as sharp magnitude peaks with damped π phase jumps, making target electronic structure visible without scanning the infrared frequency.","The logarithmic Hilbert transform check means RABBITT phase can be reconstructed from magnitude measurements in resonant two-photon ionization, potentially easing experimental phase retrieval.","The same parity separation should extend to any atom or symmetric molecule whose residual ionic state parity Pn is known.","The technique yields Fano parameters for two-photon resonances that have no synchrotron analogue, since those states decay into two non-resonant continua."],"supporting_citations":[{"why":"Supplies the APT summation rule that builds the full ionization amplitude from single-pulse amplitudes, the basis of crRABBITT.","marker":"[19]"},{"why":"Provides the multiconfiguration time-dependent Schrödinger equation method used for the two-electron targets.","marker":"[23]"},{"why":"Earlier treatment of resonant two-photon ionization that supplies the Fano-parameter extraction for XUV+IR resonances.","marker":"[24]"},{"why":"Supplies the discrete-variable-representation finite-element discretization used in the TDSE solver.","marker":"[25]"},{"why":"Supplies the implicit fourth-order time-propagation scheme used to advance the wavefunction.","marker":"[26]"},{"why":"Supplies the t-SURFFc surface-flux method that extracts ionization amplitudes from the wavefunction.","marker":"[27]"},{"why":"Supplies the logarithmic Hilbert transform / Kramers-Kronig relation that connects RABBITT phase and magnitude, which the paper tests.","marker":"[21]"},{"why":"Supplies the high-resolution experimental Fano parameters for He double-excitation states used as a benchmark.","marker":"[28]"},{"why":"Supplies the earlier theoretical Fano parameters for H2 double-excited states used for comparison.","marker":"[30]"},{"why":"Supplies the experimental near-threshold RABBITT phase data used for qualitative comparison in the below-threshold case.","marker":"[6]"}],"fun_headline_variants":["Parity trick isolates two-photon amplitude for continuous resonance mapping","Single TDSE run resolves full resonance series via parity selection","Parity-based amplitude separation maps He and H2 resonances continuously","Continuous rainbow RABBITT: parity trick ties phase to magnitude via Hilbert transform","Parity trick enables single-shot continuous resonance spectra in He and H2"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole construction assumes a target of definite parity probed by weak fields, so that symmetrizing the amplitude isolates a pure two-photon term from everything else; if parity is mixed or the infrared intensity is not weak, the extracted RABBITT parameters would be contaminated.","fun_headline_variants_meta":{"raw":{"variants":["Parity trick isolates two-photon amplitude for continuous resonance mapping","Single TDSE run resolves full resonance series via parity selection","Parity-based amplitude separation maps He and H2 resonances continuously","Continuous rainbow RABBITT: parity trick ties phase to magnitude via Hilbert transform","Parity trick enables single-shot continuous resonance spectra in He and H2"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000596,"raw_usage":{"total_tokens":2745,"prompt_tokens":856,"completion_tokens":1889,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":472,"completion_tokens_details":{"reasoning_tokens":1798}},"tokens_in":472,"tokens_out":1889,"duration_ms":13339,"temperature":1.0,"reasoning_tokens":1798,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T23:36:01.151775+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a target with no definite parity, or raise the infrared intensity until higher-order sidebands appear, and compare the parity-separated B and C against a full frequency-tuned rRABBITT scan in the same energy window; a mismatch would show the symmetrized amplitude is not purely two-photon.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the APT summation rule that builds the full ionization amplitude from single-pulse amplitudes, the basis of crRABBITT."},{"cited_title":"Time-dependent convergent close coupling method for molecular ionization in laser fields","cited_arxiv_id":"2405.12455","evidence_quote":"Provides the multiconfiguration time-dependent Schrödinger equation method used for the two-electron targets."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier treatment of resonant two-photon ionization that supplies the Fano-parameter extraction for XUV+IR resonances."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the discrete-variable-representation finite-element discretization used in the TDSE solver."},{"cited_title":"Puzynin, A","cited_arxiv_id":null,"evidence_quote":"Supplies the implicit fourth-order time-propagation scheme used to advance the wavefunction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the t-SURFFc surface-flux method that extracts ionization amplitudes from the wavefunction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the logarithmic Hilbert transform / Kramers-Kronig relation that connects RABBITT phase and magnitude, which the paper tests."},{"cited_title":"Domke, K","cited_arxiv_id":null,"evidence_quote":"Supplies the high-resolution experimental Fano parameters for He double-excitation states used as a benchmark."},{"cited_title":"S´ anchez and F","cited_arxiv_id":null,"evidence_quote":"Supplies the earlier theoretical Fano parameters for H2 double-excited states used for comparison."},{"cited_title":"Neorici´ c, D","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental near-threshold RABBITT phase data used for qualitative comparison in the below-threshold case."}],"review_version":1}