{"id":"87827305-6898-44e9-90f3-9581934e9e25","arxiv_id":"2508.08626","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Mollow sidebands in high-harmonic spectra of molecular ions survive macroscopic intensity averaging and are emitted at wider angles than the main harmonics, making them easier to isolate.","lead":"This paper simulates how the light emitted by many molecules in a gas target adds together, and shows that predicted extra spectral peaks survive and spread out at wider angles. That matters because these sidebands carry a fingerprint of the molecule's quantum motion, and the wider angles give experiments a cleaner way to detect them.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The off-axis sideband prediction depends on the vacuum-phase thin-medium assumption in Eq. (3); without a propagation check including dispersion, plasma, and Gouy phase, the central observability claim is not yet secure.","rationale":"The paper's central claim is that Mollow sidebands survive macroscopic intensity averaging and appear at larger angles. The derivation of that angular pattern is Eqs. (3)-(11), and the least secure step is the phase of the source polarization: Eq. (3) sets every spectral component to the vacuum phase e^{-i omega z/c}, and Sec. II justifies this only by z0 << z_R. That condition controls Gaussian beam parameter variation, not gas dispersion, plasma, or resonant absorption. The factorized form G(kappa) K(omega, kappa) in Eq. (7) makes the angular prediction directly dependent on this phase: the z-integral G(kappa) is a Fourier transform of rho(z) with kappa = (omega/c) sin^2(theta/2), and any additional phase Delta k(omega) z shifts the argument of G. For a density of 10^18 cm^-3 in N2, n-1 ~ 10^-5 gives Delta k z0 ~ 0.07 rad; near resonance the effect can be larger, and the intensity-dependent dipole phase, already present in D, can combine with propagation to move angular nodes. The reader is right that no full Maxwell propagation step is included. The concern is not internal inconsistency: Eqs. (1)-(11) are standard, the numerical convergence in Fig. 1 is credible, and the single-molecule physics is standard TDDFT/TDSE. The issue is an untested physical assumption about the experimental environment. Because the abstract's 'approximately the same intensity' and 'isolated more easily' claims depend on the angular pattern, this is the load-bearing point. A quantitative propagation-sensitivity check would settle whether the concern lands; if the pattern does not change, the CONDITIONAL verdict could be upgraded toward ACCEPT. The reader's weak assumption and my concern coincide, so no verdict change is needed.","tokens_in":10131,"tokens_out":11250,"duration_ms":126468,"concrete_test":"Run a propagation-sensitivity check: recompute the macroscopic angular spectrum with Eq. (3)'s phase replaced by e^{-i k(omega) z}, where k(omega) = omega n(omega)/c, n(omega) = 1 + chi_N2(omega)/2 - omega_p^2/(2 omega^2), and with the driving-field Gouy phase zeta(z) = arctan(z/z_R) added to the phase of D. Use the paper's parameters (b = 30 micron, z0 = 0.5 mm, wavelengths 446 nm and 550 nm) and gas densities between 10^17 and 10^19 cm^-3. Compare the sideband/main harmonic intensity ratio at theta = 0 and the angle of the off-axis maximum in Fig. 6. If either shifts by more than about 20%, or if the sidebands drop below 'approximately the same intensity,' the central claim needs revision; if the pattern is essentially unchanged, the concern does not land.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central observable is the angular pattern of Mollow sidebands in the macroscopic spectrum. That pattern is computed from Eq. (4) with the source polarization P(r,omega) = xhat rho(r) e^{-i omega z/c} D(E0 e^{-(x^2+y^2)/b^2}, omega), i.e., Eq. (3). This fixes every spectral component of the source to the vacuum phase e^{-i omega z/c}. The factorized angular spectrum in Eqs. (7)-(9) and the conclusion that sidebands radiate at larger angles than the main harmonics rest on that phase. In a real gas target, the fundamental and harmonic fields propagate with a frequency-dependent refractive index: neutral dispersion adds Delta k = (omega/c)(n-1), free-electron plasma contributes a negative correction, and the focused Gaussian beam adds a Gouy phase approximately z/z_R and wavefront curvature. The condition z0 << z_R stated in Sec. II controls only the z-dependence of the beam width and envelope, not these phase terms. For a nitrogen density of order 10^18 cm^-3, n-1 ~ 10^-5 gives Delta k z0 ~ 0.07 rad, and near the resonant transition the effect can be larger. Because the z-integral G(kappa) in Eq. (8) is a Fourier transform of rho(z), an additional phase Delta k(omega) z directly shifts the argument of G and can move or suppress the predicted off-axis maxima. The paper does not include a full Maxwell propagation step or a phase-matching sensitivity analysis. The same assumption also weakens the abstract's claim that sidebands have 'approximately the same intensity' as the main harmonics and could be 'isolated more easily': that claim depends on the angular separation surviving phase-matching effects, which is exactly what is untested. This is a load-bearing concern, not an internal inconsistency: the derivation from Eq. (1) to Eq. (11) is standard, and the numerical method is clearly and credibly described. The missing piece is a physical propagation check.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript extends the standard macroscopic high-harmonic generation (HHG) formulation, in which the far-field spectrum is obtained from a coherent integral over the induced polarization, to accept single-molecule dipole spectra from ab initio time-dependent calculations. The authors derive a factorized angular far-field expression for a thin slab under a Gaussian beam [Eqs. (7)-(9)], introduce an interpolation/quadrature scheme to evaluate the rapidly oscillating Bessel integral using only about 40 intensity samples, and apply the scheme to two systems: a 1D double-well model and 3D TDDFT for N2+ driven near resonance. The microscopic spectra show Mollow sidebands whose positions scale as omega = N omega0 +/- mu E, and the macroscopic spectra retain these sidebands with angular patterns that differ from the main harmonics. The paper concludes that these nonadiabatic signatures should be observable in experiments and that the framework is general.","tokens_in":10425,"tokens_out":6318,"duration_ms":70017,"significance":"If the central observability claim holds, the paper provides a useful and falsifiable prediction: Mollow sidebands in molecular-ion HHG survive macroscopic intensity averaging and appear at distinct far-field angles, offering a route to isolate them experimentally. The manuscript also contributes a practical numerical strategy for coupling expensive ab initio single-molecule calculations to the macroscopic HHG propagator, with explicit convergence evidence in Fig. 1. At the same time, the strength of the claim is limited by the model's idealized propagation assumptions and by unquantified statements about sideband intensity; the underlying single-molecule physics and the factorization derivation are standard and clearly presented.","major_comments":[{"comment":"The source polarization in Eq. (3) fixes every frequency component to the vacuum phase e^{-i omega z/c}, and the factorization in Eqs. (7)-(9) together with the predicted off-axis sideband brightening inherits this assumption. The condition z0 << zR stated in Sec. II controls only the z-dependence of the beam width and envelope; it does not control phase mismatch from neutral dispersion, free-electron plasma, or the Gouy phase. Since G(kappa) in Eq. (8) is a Fourier transform of rho(z), an additional phase Delta k(omega) z directly shifts the argument of G and can move or suppress the predicted off-axis maxima. For a target density around 10^18 cm^-3, (n-1) ~ 10^-5 gives Delta k z0 ~ 0.07 rad, and near resonances the effect can be larger. The manuscript includes no full Maxwell propagation step and no phase-matching sensitivity analysis. Please add such an analysis or explicitly restrict the observability claim to the thin-medium, vacuum-phase regime.","section":"II, Eq. (3); IV B, Fig. 6"},{"comment":"The analytical sideband positions shown as dashed red lines use omega = N omega0 +/- mu E, but the value of mu used for N2+ is not stated, and the text does not explain whether mu is computed from field-free orbitals or adjusted to match the numerical spectra. If mu is adjusted, the agreement in Fig. 4 is by construction and the Floquet interpretation loses predictive weight. Please report mu for both the 1D and TDDFT calculations, specify how it is obtained, and state explicitly that it is fixed before comparing with the computed spectra.","section":"IV A, Fig. 4"},{"comment":"The abstract's claims that sidebands have 'approximately the same intensity as the main harmonics' and are 'radiated at wider angles' are not quantified. Fig. 6 is on a logarithmic scale truncated at 0.1% of the maximum, and no lineout, integrated-power ratio, or angle-resolved contrast measure is provided. Without a quantitative definition, for example the peak sideband amplitude versus the adjacent harmonic at fixed theta or the angle-integrated power in a sideband window, the central experimental promise is not yet supported. Please add explicit quantitative comparisons and state their dependence on laser intensity, wavelength, and target thickness.","section":"Abstract; IV B, Fig. 6"}],"minor_comments":[{"comment":"The paragraph introducing Fig. 3 appears to be copied from the single-molecule discussion in Sec. IV A; Fig. 3 is the macroscopic angular spectrum for the 1D model, and the text should describe it accordingly.","section":"IV B"},{"comment":"The citations for the OEP-KLI approximation and for exceptional points/avoided crossings are missing, appearing as '[ ? ]' placeholders.","section":"III A; III C"},{"comment":"There are several typographical errors: 'a interesting angular pattern' should read 'an interesting angular pattern', and 'lest of the laser pulse' should read 'rest of the laser pulse'.","section":"Abstract; III A"},{"comment":"The caption mixes 'TDSE evaluations' and 'TDDFT evaluations' when describing the same interpolation; please use consistent terminology.","section":"Fig. 1"},{"comment":"Reference [11] has inconsistent spelling and capitalization ('Loopez Carrenno') and omits volume/page information; please standardize the bibliography format.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for physics.chem-ph and the numerical method is useful, but the observability claims in the abstract and conclusions currently outrun the propagation model used in the calculations. The missing phase-matching analysis is the main technical gap and should be addressed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this is a serious computational paper that does something new — it feeds ab initio single-molecule dipoles (TDDFT for N2+, TDSE for a 1D double well) into a macroscopic coherent sum, and it shows that Mollow sidebands can survive intensity averaging and radiate at wider angles than the main harmonics. The central physics is probably right within the stated model, and the paper deserves a proper referee.\n\nWhat's genuinely good: the interpolation scheme outlined in Sec. III B is clever and the convergence shown in Fig. 1 is convincing. The factorization into a geometric factor G and a universal factor K in Eqs. (7)–(9) is tidy, and the Floquet interpretation of the sideband positions is helpful. The computations are described in enough detail that someone could reproduce them, though no code or data are provided.\n\nThe soft spots are real but not fatal. The biggest is Eq. (3): every spectral component of the source polarization is assigned the vacuum phase e^{-iωz/c}, and the predicted angular pattern—especially the off-axis sideband brightening—follows directly from that. The paper never checks what happens with neutral dispersion, free-electron plasma, or the Gouy phase. At a gas density around 1e18 cm^-3, the phase mismatch is modest (order 0.07 rad over 0.5 mm) but not negligible, and near a resonance it can be larger. That's exactly the kind of effect that can reshape angular distributions. I'd want a sensitivity analysis or a quick full Maxwell propagation step before trusting the “isolate more easily in an experiment” claim.\n\nThe abstract says the sidebands have “approximately the same intensity” as the main harmonics; the body shows they are visible, but I didn't find a quantitative intensity comparison. Also there are mechanical defects: two “[?]” citation placeholders and a paragraph in Section IV.B that duplicates the one in Section IV.A (the Figure 3 text repeats Figure 2's text). These are easy fixes.\n\nIn short: the method is a useful step forward, and the qualitative result—sidebands survive macroscopic averaging and show angular structure—is likely robust. But the paper currently overclaims the experimental observability without testing the propagation physics. Send it to review, but ask for a phase-matching check and a cleanup.","headline":"Serious computational study showing Mollow sidebands survive macroscopic intensity averaging in molecular-ion HHG, with an off-axis angular signature, but the observability claim rests on an untested vacuum-phase propagation model and the abstract overstates the intensity comparison.","tokens_in":11097,"tokens_out":2268,"would_cite":true,"duration_ms":24039,"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":"Mollow sidebands in high-harmonic spectra from molecular ions survive macroscopic intensity averaging and appear at larger angles than the main harmonics.","keywords":["high-harmonic generation","Mollow sidebands","macroscopic response","molecular ions","nonadiabatic dynamics","TDDFT","Rabi oscillations","angle-resolved spectra"],"falsifier":"An angle-resolved high-harmonic experiment on aligned N2+ at 446 nm, or a full Maxwell propagation of the same single-molecule response including neutral dispersion and free electrons, should show whether Mollow sidebands appear at roughly main-harmonic intensity at angles of order a milliradian; their absence there would refute the paper's central claim.","tokens_in":9852,"feed_emoji":"⚛️","tokens_out":7411,"duration_ms":72432,"temperature":0.7,"pith_summary":"The paper argues that Mollow sidebands in high-harmonic spectra from molecular ions survive the averaging over laser intensity that any real macroscopic sample imposes. In both a one-dimensional double-well model and a three-dimensional TDDFT description of aligned N2+, the sidebands remain at roughly the strength of the main harmonics and are radiated at wider angles, which would make them experimentally isolable. The authors build a bridge between ab initio single-molecule spectra and far-field macroscopic spectra, rather than relying on approximate semiclassical models that miss below-threshold and excited-state physics. Their broader claim is that signatures of nonadiabatic dynamics in open-shell molecules should be observable under realistic experimental conditions.","feed_headline":"Mollow sidebands survive realistic molecular intensity averaging","feed_subtitle":"In angle-resolved N2+ spectra, sidebands match main harmonics in strength and radiate wider, easing detection.","key_machinery":"The central object is the factorization $$U(\\omega,\\hat n)= \\frac{\\$omega^{4}$ $b^{4}$}{$2c^{3}$}|\\hat n \\times \\hat x|^2\\, G\\!\\left(\\frac{\\omega}{c}\\$sin^{2}$\\frac{\\$\\theta$}{2}\\right) K\\!\\left(\\omega,\\frac{\\omega b}{c}\\sin\\$\\theta$\\right)$$, which separates the target geometry $G$ from a universal molecular-response factor $K$. The factor $K$ is a Bessel-weighted integral of the single-molecule polarization over the Gaussian intensity profile of the driving beam, and the authors evaluate it efficiently by interpolating $D(E,\\omega)$ in the electric field at Chebyshev nodes and integrating with Bessel-root quadrature. This factorization is what lets expensive ab initio single-molecule calculations be converted into far-field angular spectra.","core_discovery":"The central discovery is that the Mollow sidebands predicted in the single-molecule high-harmonic response are not washed out by macroscopic intensity averaging. For N2+ driven at 446 nm and for a resonant one-dimensional double-well model, the macroscopic angle-resolved spectrum retains the sidebands around the main harmonics at comparable intensity, while radiating at larger angles and bending slightly inward. The sidebands also carry different group delays from the main harmonics, and several intensity-dependent features such as the bifurcation near $6\\omega_0$ survive with intricate angular patterns. The paper concludes that nonadiabatic molecular dynamics leaves measurable fingerprints in macroscopic high-harmonic spectra.","pith_inferences":["The thin-medium factorization assumes vacuum phase velocity for every frequency component; a full Maxwell propagation including neutral dispersion and free-electron plasma could shift the off-axis brightening, so a realistic 3D propagation test would show how robust the wider-angle sidebands are.","The same wider-angle fingerprint could be used to spatially gate high-harmonic spectra of other open-shell molecules with resonantly coupled states, not just N2+, provided their Rabi-frequency sidebands are strong enough.","The angular separation between sidebands and main harmonics likely scales with the Rabi frequency, suggesting that measuring the emission angle could extract transition dipole moments, an application the paper does not discuss."],"forward_implications":["Macroscopic spectra of aligned N2+ retain Mollow sidebands around the 5th and 11th harmonics at intensities comparable to the main harmonics, so the sidebands should be detectable without isolating single molecules.","Because the sidebands radiate at wider angles than the main harmonics, a spatial aperture or imaging setup could separate the nonadiabatic-dynamics signatures from ordinary harmonics.","The same method can promote any ab initio single-molecule spectrum, including TDDFT, to a macroscopic prediction, replacing semiclassical models where below-threshold, excited-state, or Rydberg physics matters.","Harmonic phase and group-delay maps retain angular structure, so far-field angular measurements can carry time-delay information about the underlying nonadiabatic dynamics."],"supporting_citations":[{"why":"Predicted Mollow sidebands in molecular-ion high-harmonic generation, the microscopic effect whose macroscopic fate is tested.","marker":"[1]"},{"why":"Provides the analytic sideband position as $\\omega = N\\omega_0 \\pm \\mu E$, used to label the spectral features.","marker":"[2]"},{"why":"Defines the Mollow triplet in quantum optics, the namesake of the sidebands.","marker":"[3]"},{"why":"Documents the failure of semiclassical strong-field models for below-threshold and excited-state features, motivating full TDSE and TDDFT input.","marker":"[12]"},{"why":"Supplies the TDDFT implementation used to compute the N2+ microscopic spectra.","marker":"[14-16]"},{"why":"Supplies Bessel-root quadrature, making the macroscopic integral converge with few microscopic evaluations.","marker":"[18]"}],"fun_headline_variants":["Mollow sidebands survive averaging, radiate wider in N2+ HHG","Sidebands match main harmonics, spread wider in molecular HHG","Nonadiabatic fingerprints in HHG: sidebands resist averaging","Wide-angle Mollow sidebands ease detection in N2+ HHG"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes each frequency component of the source travels through the gas at the vacuum speed of light, with no dispersion, plasma, or frequency-dependent refractive index; if phase matching in a real gas jet breaks this assumption, the predicted off-axis sideband pattern could change.","fun_headline_variants_meta":{"raw":{"variants":["Mollow sidebands survive averaging, radiate wider in N2+ HHG","Sidebands match main harmonics, spread wider in molecular HHG","Nonadiabatic fingerprints in HHG: sidebands resist averaging","Wide-angle Mollow sidebands ease detection in N2+ HHG"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000577,"raw_usage":{"total_tokens":2667,"prompt_tokens":834,"completion_tokens":1833,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":450,"completion_tokens_details":{"reasoning_tokens":1763}},"tokens_in":450,"tokens_out":1833,"duration_ms":14947,"temperature":1.0,"reasoning_tokens":1763,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:35:21.767601+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An angle-resolved high-harmonic experiment on aligned N2+ at 446 nm, or a full Maxwell propagation of the same single-molecule response including neutral dispersion and free electrons, should show whether Mollow sidebands appear at roughly main-harmonic intensity at angles of order a milliradian; their absence there would refute the paper's central claim.","supporting_citations":[{"cited_title":"Xia and A","cited_arxiv_id":null,"evidence_quote":"Predicted Mollow sidebands in molecular-ion high-harmonic generation, the microscopic effect whose macroscopic fate is tested."},{"cited_title":"Joyce, A","cited_arxiv_id":null,"evidence_quote":"Provides the analytic sideband position as $\\omega = N\\omega_0 \\pm \\mu E$, used to label the spectral features."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Mollow triplet in quantum optics, the namesake of the sidebands."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the failure of semiclassical strong-field models for below-threshold and excited-state features, motivating full TDSE and TDDFT input."},{"cited_title":"Ogata, Publ","cited_arxiv_id":null,"evidence_quote":"Supplies Bessel-root quadrature, making the macroscopic integral converge with few microscopic evaluations."}],"review_version":2}