Macroscopic properties of high-harmonic generation from molecular ions

A. Jaron; T Joyce

REVIEW 3 major objections 5 minor 19 references

Macroscopic properties of high-harmonic generation from molecular ions

T0 review · 3 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Mollow sidebands in high-harmonic spectra from molecular ions survive macroscopic intensity averaging and appear at larger angles than the main harmonics.

desk verdict 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. read the letter →

arxiv 2508.08626 v1 pith:XMEIVRGW submitted 2025-08-12 physics.chem-ph physics.app-phphysics.atm-clusquant-ph

T Joyce , A. Jaron This is my paper · ORCID

classification physics.chem-phphysics.app-phphysics.atm-clusquant-ph

keywords high-harmonicgenerationMollowsidebandsmacroscopicresponsemolecularionsnonadiabaticdynamicsTDDFTRabioscillationsangle-resolvedspectra

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the 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.

What'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.

The 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.

The 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.

In 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.

Referee Report

3 major / 5 minor

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.

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 (3)

[II, Eq. (3); IV B, Fig. 6] 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.
[IV A, Fig. 4] 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.
[Abstract; IV B, Fig. 6] 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.

minor comments (5)

[IV B] 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.
[III A; III C] The citations for the OEP-KLI approximation and for exceptional points/avoided crossings are missing, appearing as '[ ? ]' placeholders.
[Abstract; III A] 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'.
[Fig. 1] The caption mixes 'TDSE evaluations' and 'TDDFT evaluations' when describing the same interpolation; please use consistent terminology.
[References] Reference [11] has inconsistent spelling and capitalization ('Loopez Carrenno') and omits volume/page information; please standardize the bibliography format.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the macroscopic spectrum is a deterministic transform of independently computed single-molecule dipoles; no fitted parameter or self-citation chain is load-bearing.

full rationale

The paper's derivation chain is self-contained. The far-field spectrum U(omega,n-hat) follows from a standard electrodynamic formula (Eq. 1) and the thin-medium source model (Eq. 3), where the single-molecule response D(E,omega) is computed directly from TDSE for the 1D double well (Eq. 15) and from TDDFT for N2+ (Eq. 23). The macroscopic spectrum (Eq. 4) is a coherent, deterministic superposition of these independently computed dipoles weighted by the beam geometry; no macroscopic data are fitted and no parameter is adjusted to force the sidebands to appear. The factorization in Eqs. (7)-(9) is an exact rearrangement under the stated slab-geometry assumption, and the survival of the sidebands is a genuine result rather than an input: the integral over E involves the complex, intensity-dependent dipole phase (e.g., Fig. 5), so cancellations are possible, and the sideband positions shift with E via omega = N omega0 +/- mu E, making their persistence in the averaged spectrum nontrivial. Refs. [1,2] are self-citations, but they are contextual: the paper recomputes the microscopic Mollow sidebands ab initio and does not rely on those references as the proof of the main claim, nor does it import any uniqueness theorem from the authors' prior work. The thin-medium assumption (z0 << zR, with z0 = 0.5 mm vs zR ~ 5 mm) is a physical approximation that affects the accuracy of the angular prediction, but it is not circular: it does not define the output in terms of itself. No renaming of a known result, no fitted input called a prediction, and no ansatz smuggled in via citation are present.

Assumptions & free parameters 6 free parameters · 6 assumptions · 0 invented entities

The central calculation is a forward simulation with no external data fitting. The main input choices are physical parameters of the laser and target, plus the TDDFT functional and the 1D model potential. The only potential free parameter is the transition dipole mu used in the Floquet sideband formula, since the paper does not state whether it is computed or adjusted. The key hidden assumption is the vacuum-phase, no-dispersion propagation model in Eq. (3).

free parameters (6)

Transition dipole moment mu for Floquet sideband slopes = not stated
The analytic sideband lines use omega = N omega0 +/- mu E. The paper does not state whether mu is computed from the ab initio wavefunctions or adjusted to match the numeric spectra. If adjusted, the analytic 'prediction' is a fit.
Beam waist b = 30 micrometers
Chosen as an experimentally realistic input; sets the angular scale of K(omega,kappa). The paper says results are insensitive to macroscopic parameters within the assumed range.
Target thickness z0 = 0.5 mm
Chosen to satisfy z0 << zR for the thin-medium factorization. Affects the geometric factor G and the validity of the slab approximation.
Pulse duration tau = 400 a.u. (1D), 200 a.u. (N2+)
Chosen laser envelope parameter. The macroscopic calculation assumes the same temporal envelope at every point in the target.
Sweep range for peak field E0 = 0.03 a.u. (1D), 0.05 a.u. (N2+)
The intensity range covered in the macroscopic average. Not fitted to data; chosen to cover the strong-field regime.
1D double-well potential parameters = wells at x = +/-2, Gaussian depth 1 a.u.
A model potential chosen to represent a generic aligned diatomic molecule. Not fitted to an experimental target.

assumptions (6)

standard math Far-field radiation formula and dipole approximation for the source polarization (Eqs. 1 and 3)
The macroscopic spectrum is computed from the induced polarization via classical electrodynamics. This is standard but carries the dipole-approximation and local-response assumptions.
domain assumption Thin medium (z0 << zR) and slab geometry allow factorization into G and K (Eqs. 7-9)
Stated in Sec. II. Requires the medium thickness to be much smaller than the Rayleigh length and the density to depend only on z. The paper chooses z0 = 0.5 mm to satisfy this.
domain assumption The source polarization phase is e^{-i omega z / c}, i.e., no gas dispersion, plasma, or phase mismatch
Eq. (3) assumes both driving field and emitted field propagate at vacuum speed c. This neglects Delta-k effects and can change the angular pattern of harmonics in a real medium.
domain assumption D(E,omega) is an odd function of E, enabling symmetry-reduced Chebyshev interpolation (Eq. 32)
Requires inversion symmetry along the laser polarization axis. True for the symmetric 1D double well and for N2+ oriented along x.
domain assumption Floquet picture: resonant driving creates an exceptional point at Ec = 0, and two Floquet states dominate (Eqs. 34-39)
Used to identify sidebands at omega = N omega0 +/- mu E. The paper acknowledges higher Rydberg states and 3-level dynamics complicate this simple picture.
domain assumption TDDFT with PZ-SIC and OEP-KLI gives reliable single-molecule HHG dipoles for N2+
The macroscopic result inherits all errors from the microscopic input. TDDFT accuracy for strong-field ionization, especially near Rydberg states, is an approximation.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Macroscopic properties of high-harmonic generation from molecular ions." pith.science (2026). https://pith.science/paper/XMEIVRGW

@misc{pith2026250808626,
  author       = {Pith},
  title        = {Pith review of: Macroscopic properties of high-harmonic generation from molecular ions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XMEIVRGW}},
  note         = {Machine review of arXiv:2508.08626}
}

read the original abstract

We extend the existing framework of macroscopic HHG to combine it with high-accuracy ab initio calculations for molecules as microscopic input. This approach is applied to HHG spectra exhibiting Mollow sidebands, for open shell molecules undergoing nonadiabatic dynamics. We demonstrate the details of the method and analyze how the predicted features in the microscopic HHG response unambiguously survive macroscopic response calculations, and furthermore they exhibit a interesting angular pattern in the far-field. We calculate the macroscopic harmonic spectrum by combining many single-molecule calculations at different intensities, obtained in one case from time-dependent density functional theory calculations for N+ 2 , in second case for one electron time dependent Schr\"odinger equation for a 1D double well model potential. For both cases one can observe that the resulting macroscopic spectra exhibit Mollow sidebands of approximately the same intensity as the main harmonics, while being radiated at wider angles, meaning they could be isolated more easily in an experiment.

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

19 extracted references · 19 canonical work pages

[1]

Xia and A

Y. Xia and A. Jaron-Becker, Opt. Express 24, 4689-4697 (2016)

work page 2016
[2]

Joyce, A

T. Joyce, A. Jaron-Becker, Optics Letters 45, 1954 (2020)

work page 2020
[3]

B. R. Mollow, Phys. Rev.188, 1969–1975 (1969)

work page 1969
[4]

Stalgies, I

Y. Stalgies, I. Siemers, B. Appasamy, T. Altevogt,and P. E. Toschek, EPL (Europhysics Letters) 35, 259 (1996)

work page 1996
[5]

Wrigge, I

G. Wrigge, I. Gerhardt, J. Hwang, G. Zumofen, andV. Sandoghdar, Nature Physics 4, 60 (2007)

work page 2007
[6]

Ulhaq, S

A. Ulhaq, S. Weiler, S. M. Ulrich, R. Roßbach, M. Jet- ter, and P. Michler, Nature Photonics 6, 238 (2012)

work page 2012
[7]

E. B. Flagg, A. Muller, J. W. Robertson, S. Founta, D. G.Deppe, M. Xiao, W. Ma, G. J. Salamo, and C. K. Shih, Nature Physics 5, 203 (2009)

work page 2009
[8]

Peiris, B

M. Peiris, B. Petrak, K. Konthasinghe, Y. Yu, Z. C.Niu, and A. Muller, Phys. Rev. B91, 195125 (2015)

work page 2015

Show all 19 references

[9]

K. G. Lagoudakis, K. A. Fischer, T. Sarmiento,P. L. 9 McMahon, M. Radulaski, J. L. Zhang, Y. Kelaita, C. Dory, K. M¨ uller, and J. Vuckovic, Phys. Rev. Lett. 118, 013602 (2017)

2017
[10]

Aspect, G

A. Aspect, G. Roger, S. Reynaud, J. Dalibard, and C. Cohen-Tannoudji,Phys. Rev. Lett. 45, 617–620 (1980)

1980
[11]

J. C. Loopez Carrenno, E. del Valle, and F. P. Laussy, Laser & Photonics Reviews 11, 1700090
[12]

Y. Xia, A. Jaron-Becker, Optics Letters 39, 1461 (2014)

2014
[13]

Cheng Jin, Anh-Thu Le, and C. D. Lin, Phys. Rev. A 83, 023411 (2011)

2011
[14]

Marques, A

M.A.L. Marques, A. Castro, G.F. Bertsch, and A. Rubio, Comput. Phys. Commun. 151, 60 (2003)

2003
[15]

Castro, H

A. Castro, H. Appel, M. Oliveira, C.A. Rozzi, X. An- drade, F. Lorenzen, M.A.L. Marques, E.K.U. Gross, and A. Rubio, Phys. Stat. Sol. B 243, 2465 (2006)

2006
[16]

Andrade, D.A

X. Andrade, D.A. Strubbe, U. De Giovannini, A.H. Larsen, M.J.T. Oliveira, J. Alberdi-Rodriguez, A. Varas, I. Theophilou, N. Helbig, M. Verstraete, L. Stella, F. Nogueira, A. Aspuru-Guzik, A. Castro, M.A.L. Mar- ques, and A. Rubio, Phys. Chem. Chem. Phys. 17 31371 (2015)

2015
[17]

J. B. Krieger, Yan Li, and G. J. Iafrate, Phys. Rev. A, 46, 5453, (1992)

1992
[18]

Ogata, Publ

H. Ogata, Publ. Res. Inst. Math. Sci., 41, 949, (2005)

2005
[19]

See. e.g.M. M¨ uller, I. Rotter, Journal of Physics A: Math- ematical and Theoretical 41, 244018 (2008)

2008

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.

[1] [1]

Xia and A

Y. Xia and A. Jaron-Becker, Opt. Express 24, 4689-4697 (2016)

work page 2016

[2] [2]

Joyce, A

T. Joyce, A. Jaron-Becker, Optics Letters 45, 1954 (2020)

work page 2020

[3] [3]

B. R. Mollow, Phys. Rev.188, 1969–1975 (1969)

work page 1969

[4] [4]

Stalgies, I

Y. Stalgies, I. Siemers, B. Appasamy, T. Altevogt,and P. E. Toschek, EPL (Europhysics Letters) 35, 259 (1996)

work page 1996

[5] [5]

Wrigge, I

G. Wrigge, I. Gerhardt, J. Hwang, G. Zumofen, andV. Sandoghdar, Nature Physics 4, 60 (2007)

work page 2007

[6] [6]

Ulhaq, S

A. Ulhaq, S. Weiler, S. M. Ulrich, R. Roßbach, M. Jet- ter, and P. Michler, Nature Photonics 6, 238 (2012)

work page 2012

[7] [7]

E. B. Flagg, A. Muller, J. W. Robertson, S. Founta, D. G.Deppe, M. Xiao, W. Ma, G. J. Salamo, and C. K. Shih, Nature Physics 5, 203 (2009)

work page 2009

[8] [8]

Peiris, B

M. Peiris, B. Petrak, K. Konthasinghe, Y. Yu, Z. C.Niu, and A. Muller, Phys. Rev. B91, 195125 (2015)

work page 2015