Pith. sign in

REVIEW 3 major objections 5 minor 28 references

Versatile, open-source program for simulating high-harmonic generation

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

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

Pith's one-line read A single open-source program simulates high-harmonic generation across the full range of modern table-top sources, from few-cycle 800 nm drivers to multi-bar soft-X-ray generation, by coupling the microscopic single-atom response with macro

desk verdict Useful open-source HHG simulator, but Eq. (7) omits the Lewenstein prefactor and validation is only qualitative; the accuracy claim doesn't hold yet. read the letter →

arxiv 2509.02867 v1 pith:2HRROTIW submitted 2025-09-02 physics.optics physics.atom-ph

Christian A. Schröder , Reinhard Kienberger , Stephen R. Leone This is my paper · ORCID
classification physics.opticsphysics.atom-ph PACS 42.65.Ky
keywords high-harmonicgenerationmacroscopicpropagationLewensteinmodelphasematchingisolatedattosecondpulsestwo-colordrivingfieldssoftX-raysourcesopen-sourcesimulation
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 presents a single, open-source program that simulates high-harmonic generation over the parameter range used by modern table-top XUV and soft-X-ray sources: few-cycle 800 nm drivers and multi-cycle 1300 nm pulses, gas pressures from 100 mbar to 4 bar, and single- or multi-color fields. The program works by propagating the fundamental and harmonic fields together while recomputing the per-atom nonlinear response from the instantaneous driving field, so plasma generation, Kerr self-focusing, ionization losses, dispersion, and absorption all affect the same run as the harmonic generation itself. It is validated on four established cases: phase-matching maxima as a short gas cell moves through the focus, isolated attosecond pulse synthesis, even-order harmonics from a two-color drive, and soft-X-ray spectra from high-pressure helium. A reader designing or improving such a source would care because the tool offers one documented, BSD-licensed calculation path instead of a chain of separate codes.

What carries the argument

The central mechanism is the coupled solution of the paraxial propagation equation in a frame traveling with the laser pulse, applied separately to the fundamental and harmonic fields and advanced in lockstep. At each z-plane the microscopic response is computed from the instantaneous fundamental field: a Lewenstein-model integral gives the single-atom dipole, and multiplying by the neutral fraction and gas density gives the macroscopic harmonic polarization. The fundamental field, in turn, receives Kerr self-focusing, plasma defocusing, and ionization-loss currents derived from ADK ionization rates. A Crank-Nicholson integrator advances the fields, while Sellmeier formulas and tabulated opt

What would settle it

Run the program on the exact parameters of a well-characterized HHG source—pulse energy, duration, spectral phase, beam waist, gas type, cell length, pressure, and cell position—and compare the calculated harmonic spectrum, including the relative plateau-to-cutoff yield and the pressure or position of maximum yield, against measured spectra using a defined error metric. If the yields differ by orders of magnitude or the phase-matching maxima shift significantly, the paper's accuracy claim would fail.

Watch

Extended reading notes

Core claim

The paper's central claim is that a monolithic HHG simulator can accurately take both macroscopic and microscopic aspects into account and is applicable across the broad parameter range of current HHG-based light sources. On the authors' terms, that claim is established by the program's design—solving the propagation equation in a frame moving with the pulse, with the fundamental and harmonic fields advanced in lockstep—and by the four validation examples, which reproduce the expected phase-matching positions and beam-shape changes, an isolated sub-femtosecond XUV pulse from amplitude gating, even-order harmonics in a two-color driving field, and soft-X-ray spectra in high-pressure helium. T

Load-bearing premise

The central accuracy claim stands on the premise that matching previously published simulations and expected physics broadly—without a quantitative comparison to a measured spectrum or an independent reference code—is enough to prove the program is correct across the entire stated parameter range.

Editorial extensions

If this is right

  • A single calculation can track driving-pulse reshaping and harmonic generation together, so source design no longer requires stitching separate modules for the fundamental and harmonic fields.
  • The program reproduces the two phase-matching maxima when a short gas cell is scanned across focus, along with the associated near-field beam-profile changes, supporting gas-cell placement studies.
  • It reproduces isolated-attosecond-pulse conditions for a few-cycle pulse with selected carrier-envelope phase and spectral filtering, supporting pulse-synthesis design.
  • It predicts even-order harmonics for a fundamental-plus-second-harmonic drive and supports an arbitrary number of driving colors with independent spectral phases.
  • It can simulate soft-X-ray generation at multi-bar pressures and produce coherence-length maps that identify where to place the exit face of a semi-infinite gas cell for optimal outcoupling.

Reading between the lines

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

  • A natural next step beyond the paper is to benchmark the code's absolute or relative harmonic yields against calibrated experimental photon spectra; the qualitative validations shown would then become a quantitative test of the stated accuracy claim.
  • Because the harmonic field does not feed back into the fundamental, the model's boundary is most likely to appear in regimes where frequency mixing of the strong driving field is non-negligible; high-intensity multi-color driving would be the sharpest test.
  • The one-text-file configuration and monolithic structure suggest the program could be used for automated parameter scans and optimization of phase-matching conditions, an application the authors imply but do not demonstrate.
  • Coherence-length maps generated by the code could be used as a fast design layer for semi-infinite gas cells, decoupling the phase-matching analysis from full harmonic-spectrum calculations.
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.

Referee Report

3 major / 5 minor

Summary. The manuscript presents an open-source C++ program (gitlab.com/caschroeder/hhg) for simulating macroscopic high-harmonic generation. The model follows the standard Lewenstein single-atom response (Eqs. 7–9) coupled to a paraxial propagation equation (Eq. 1) with plasma, Kerr, and linear dispersion effects. Four example calculations are presented: phase-matching scans (Sec. 4.1), isolated attosecond pulse synthesis (Sec. 4.2), two-color driving (Sec. 4.3), and soft-X-ray HHG (Sec. 4.4), plus a coherence-length analysis (Sec. 4.5). The paper claims the program accurately accounts for both microscopic and macroscopic aspects across a broad parameter range, and that the examples demonstrate validity.

Significance. If the accuracy claim is supported, this would be a valuable community resource: a single monolithic, open-source simulator that couples single-atom strong-field response with macroscopic propagation, without interfacing external modules, applicable to both low- and high-pressure regimes. The open-source BSD-licensed code with a git repository and documentation is a concrete strength, as is the demonstrated capability to model multi-color driving and ionization-induced reshaping. However, the validation in Sections 4.1–4.4 is entirely qualitative, and the microscopic model as written contains a concrete omission from the Lewenstein expression. The significance of the contribution therefore depends on whether these issues are corrected; the tool could be widely used, but its current presentation does not establish the quantitative accuracy asserted in the abstract.

major comments (3)
  1. The single-atom dipole in Eq. (7) is written without the Lewenstein prefactor (π/(ε+iτ/2))^{3/2} (see Refs. [3,17]). This prefactor is not a normalization detail: it arises from the Gaussian integration over the two transverse momentum components, and it strongly affects the τ-weighting of the integral. Without it, long/short trajectory contributions are weighted incorrectly and absolute yields are wrong; it also provides the regularization that makes the oscillatory τ integral well behaved on the real axis. If the code actually includes the prefactor, Eq. (7) is a misprint that makes the printed model non-reproducible; if the code follows Eq. (7) literally, the quantitative accuracy claim is unsupported. Either way, this is a load-bearing inconsistency in the central model equation.
  2. The validation is qualitative. Section 4.1 states spectra “compare well” to Balcou et al., and Section 4.3 states even harmonics are “the expected result, as validated.” No error metric, no comparison to an experimental measurement, and no independent reference code are provided. The qualitative features (positions of phase-matching maxima, presence of even harmonics, smooth cutoff) are insensitive to the kind of quantitative errors that the omitted prefactor in Eq. (7) would introduce. To support the abstract's claim that the program “accurately” takes both macroscopic and microscopic aspects into account, the authors should add a quantitative benchmark—for example, comparing absolute or relative harmonic yields to experimental data or to a published independent simulation code—or temper the accuracy claim to match the evidence.
  3. The paper does not report convergence tests or numerical parameter sensitivity. Since the code introduces a Crank-Nicholson integrator, adaptive Runge-Kutta steps, and user-defined grids, it is important to show that the presented spectra are converged with respect to the radial grid, axial step, and time/frequency sampling, particularly for the high-pressure (2–4 bar) case in Sec. 4.4. Without such tests, the “validity” demonstrated by the examples is not distinguishable from numerical artifact. This is a load-bearing point for a program whose main purpose is predictive simulation.
minor comments (5)
  1. The caption of Fig. 1 refers to panels (c-f) while the text mentions (c-f) but the spatial profiles for the cut-off region are labeled (e) and (f); check panel order and cross-references.
  2. The phrase “as validated” in Sec. 4.3 seems to refer to prior literature; it would be clearer to state that even harmonics are a well-known consequence of a broken inversion symmetry in a two-color field, rather than treating the calculation as its own validation.
  3. Typos: “uniqe” (Abstract), “developements” (Sec. 1), “non-pertutbative” (Sec. 1), “righ” (Sec. 2.1.1), “c.f.” (Fig. 1 caption). These are presentation issues only.
  4. Eq. (1) is stated in a moving frame, but the relation between the retarded time coordinate and the lab frame is not defined; a reader cannot determine how the group-velocity term is derived. Adding one sentence of context would help.
  5. In Eq. (12), τ=5.08 and τ=3.16 are from Balcou et al. for specific conditions; it would be useful to state that these coefficients are empirical and cite the relevant limits of validity.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the program implements published HHG models with externally fixed parameters and validates against independent qualitative benchmarks; no prediction reduces to a fitted input or to a self-citation chain.

full rationale

The paper's derivation chain is a direct implementation of established models, not a circular argument. Eq. (1) is the standard traveling-frame propagation equation; Eqs. (3)-(10) are the standard Kerr, plasma, and Lewenstein source terms taken from the cited literature (Refs. [3,5,15]). The fixed coefficients τ = 5.08 and 3.16 in Eq. (12) are taken explicitly from Balcou et al. [4], an external source, so no parameter is fitted to the paper's own outputs. The validation examples in Secs. 4.1-4.4 compare qualitative features (phase-matching maxima, presence of even harmonics, cutoff smoothness, isolated attosecond pulse formation) against prior published calculations and expected physics; they do not involve fitting a model to data and then predicting that same data. The self-citations (Refs. [16,21,25]) are background references or optional corrections and are not load-bearing for the central claim. The skeptic's point about the omitted (π/(ε+iτ/2))^{3/2} prefactor in Eq. (7) is a possible correctness or misprint issue in the printed equation, not a circularity: it does not make any output equal to an input by construction. Under the review rules, a correctness concern of that kind belongs outside the circularity score, which is therefore 0.

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

The central claim rests on models imported from the cited literature: the Lewenstein single-atom response (ref. 3), ADK ionization rates (ref. 15), the Gaarde macroscopic propagation framework (ref. 5), and Balcou phase-matching coefficients (ref. 4). The only fixed numerical coefficients introduced in the paper's own analysis, tau = 5.08 and 3.16 in Eq. 12, are copied from ref. 4 rather than fitted. The axioms below are the domain assumptions that these standard models remain quantitatively valid across the extended parameter range (high pressure, high intensity) the program claims to cover.

free parameters (2)
  • tau coefficient for cutoff harmonics = 5.08
    Constant in Eq. 12 mapping the ponderomotive potential gradient to the harmonic wave vector, taken from the Balcou model (ref. 4) and used for the coherence-length map in Section 4.5. Not fitted in this paper.
  • tau coefficient for plateau harmonics = 3.16
    Second fixed constant from ref. 4 used in Eq. 12 for plateau-region harmonics in Section 4.5.
assumptions (6)
  • domain assumption The Lewenstein strong-field approximation yields quantitatively accurate single-atom HHG responses in all simulated regimes
    Eq. 7 is adopted from refs. 3 and 17 without re-derivation; its validity at 2 to 4 bar and 1.3e15 W/cm2 (Section 4.4) is assumed.
  • domain assumption ADK tunnel ionization rates are accurate across the parameter range
    Ionization rates in Eqs. 5 and 6 come from ref. 15; known deviations in transitional regimes are not discussed.
  • domain assumption The hydrogenic dipole of Eq. 9 is adequate for the target gases
    One hydrogenic transition dipole is used for all species (helium in every example), without species-specific corrections.
  • domain assumption Harmonic back-action on the fundamental is negligible
    Stated in Section 2: frequency mixing of the weak harmonic field with the driving field is 'expected to be negligible', so the two fields are propagated one-way in lockstep.
  • domain assumption First-order Neumann boundary conditions at r_max do not introduce significant artifacts
    Section 2 states Neumann conditions are enforced at r = r_max with first-order accuracy; no absorbing boundary layer or reflection test is reported.
  • standard math Ideal gas behavior relates pressure and density
    Section 3 computes gas density from pressure and temperature via the ideal gas law.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Versatile, open-source program for simulating high-harmonic generation." pith.science (2026). https://pith.science/paper/2HRROTIW

@misc{pith2026250902867,
  author       = {Pith},
  title        = {Pith review of: Versatile, open-source program for simulating high-harmonic generation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2HRROTIW}},
  note         = {Machine review of arXiv:2509.02867}
}
read the original abstract

Light sources based on high-harmonic generation (HHG) underpin ultrafast spectroscopy experiments across a large range of photon energies, spanning from the extreme-ultraviolet to the soft x-ray. To this day their design, implementation and improvement presents uniqe challenges, but can be aided by numerical tools. Here we present a new simulation program designed for this purpose, which accurately takes both macroscopic and microscopic aspects of high-harmonic generation into account and is therefore applicable across the broad range of parameters that HHG based light sources are today utilized in. The program is validated by calculating harmonic emission in four common experimental configurations.

Figures

Figures reproduced from arXiv: 2509.02867 by the authors.

Figure 1
Figure 1. Analysis of phase-matching conditions with a few-cycle pulse at a driving [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Simulation of isolated attosecond pulse production. In the left panel shows the [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Simulation of two-color high-harmonic generation with a long ( [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: Soft x-ray high harmonic generation calculated with our simulation program. [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: Intensity (left panel), phase evolution (center panel) of the driving pulses central [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages

  1. [1]

    Multiple-harmonic conversion of 1064 nm radiation in rare gases,

    M. Ferray, A. L’Huillier, X. F. Li,et al., “Multiple-harmonic conversion of 1064 nm radiation in rare gases,” J. Phys. B: At. Mol. Opt. Phys.21, L31 (1988)

    work page 1988

  2. [2]

    Multiple-harmonic generation in rare gases at high laser intensity,

    X. F. Li, A. L’Huillier, M. Ferray,et al., “Multiple-harmonic generation in rare gases at high laser intensity,” Phys. Rev. A39, 5751–5761 (1989)

    work page 1989

  3. [3]

    Theory of high-harmonic generation by low-frequency laser fields,

    M. Lewenstein, P. Balcou, M. Y. Ivanov,et al., “Theory of high-harmonic generation by low-frequency laser fields,” Phys. Rev. A49, 2117–2132 (1994)

    work page 1994

  4. [4]

    Generalizedphase-matchingconditionsforhighharmonics: The role of field-gradient forces,

    P.Balcou,P.Sali‘eres,A.L’Huillier,andM.Lewenstein,“Generalizedphase-matchingconditionsforhighharmonics: The role of field-gradient forces,” Phys. Rev. A55, 3204–3210 (1997)

    work page 1997

  5. [5]

    Macroscopic aspects of attosecond pulse generation,

    M. B. Gaarde, J. L. Tate, and K. J. Schafer, “Macroscopic aspects of attosecond pulse generation,” J. Phys. B: At. Mol. Opt. Phys.41, 132001 (2008)

    work page 2008

  6. [6]

    Scale-invariant nonlinear optics in gases,

    C. M. Heyl, H. Coudert-Alteirac, M. Miranda,et al., “Scale-invariant nonlinear optics in gases,” Optica3, 75–81 (2016)

    work page 2016

  7. [7]

    High-flux soft x-ray harmonic generation from ionization-shaped few-cycle laser pulses,

    A. S. Johnson, D. R. Austin, D. A. Wood,et al., “High-flux soft x-ray harmonic generation from ionization-shaped few-cycle laser pulses,” Sci. Adv.4, eaar3761 (2018)

    work page 2018

  8. [8]

    Multiscale approach to the description of high-harmonics generation in gases,

    J. Vábek, “Multiscale approach to the description of high-harmonics generation in gases,” Theses, Université de Bordeaux ; Université technique tchèque. Faculté de Génie nucléaire et de Sciences de l’ingénieur (Prague) (2022)

    work page 2022

Show all 28 references
  1. [9]

    Multi-scale model of HHG in gases,

    J. Vabek, S. Skupin, and F. Catoire, “Multi-scale model of HHG in gases,” in3rd Annual Workshop of the COST Action Attosecond Chemistry 2022, (Prague, Czech Republic, 2022)

    2022

  2. [10]

    High-order harmonic propagation in gases within the discrete dipole approximation,

    C. Hernández-García, J. A. Pérez-Hernández, J. Ramos,et al., “High-order harmonic propagation in gases within the discrete dipole approximation,” Phys. Rev. A82(2010)

    2010

  3. [11]

    Högner, https://mhoegner.gitlab.io/hhgmax-homepage/ (2014)

    M. Högner, https://mhoegner.gitlab.io/hhgmax-homepage/ (2014)

    2014

  4. [12]

    Dispersion measurement of inert gases and gas mixtures at 800 nm,

    A. Börzsönyi, Z. Heiner, M. P. Kalashnikov,et al., “Dispersion measurement of inert gases and gas mixtures at 800 nm,” Appl. Opt.47, 4856–4863 (2008)

    2008

  5. [13]

    X-ray interactions: Photoabsorption, scattering, transmission, and reflection at e = 50-30,000 ev, z = 1-92,

    B. Henke, E. Gullikson, and J. Davis, “X-ray interactions: Photoabsorption, scattering, transmission, and reflection at e = 50-30,000 ev, z = 1-92,” At. Data Nucl. Data Tables54, 181–342 (1993)

    1993

  6. [14]

    Nonresonant third order hyperpolarizability of rare gases and n2 determined by third harmonic generation,

    H. Lehmeier, W. Leupacher, and A. Penzkofer, “Nonresonant third order hyperpolarizability of rare gases and n2 determined by third harmonic generation,” Opt. communications56, 67–72 (1985)

    1985

  7. [15]

    Tunnel ionization of complex atoms and atomic ions in electromagnetic field,

    M. V. Ammosov, N. B. Delone, and V. P. Krainov, “Tunnel ionization of complex atoms and atomic ions in electromagnetic field,” inHigh intensity laser processes, vol. 664 (SPIE, 1986), pp. 138–141

    1986

  8. [16]

    Spatial separation of large dynamical blueshift and harmonic generation,

    M. B. Gaarde, M. Murakami, and R. Kienberger, “Spatial separation of large dynamical blueshift and harmonic generation,” Phys. Rev. A74, 053401 (2006)

    2006

  9. [17]

    Theory of high-order harmonic generation by an elliptically polarized laser field,

    P. Antoine, A. L’Huillier, M. Lewenstein,et al., “Theory of high-order harmonic generation by an elliptically polarized laser field,” Phys. Rev. A53, 1725–1745 (1996)

    1996

  10. [18]

    Practical sparse matrices in c++ with hybrid storage and template-based expression optimisation,

    C. Sanderson and R. Curtin, “Practical sparse matrices in c++ with hybrid storage and template-based expression optimisation,” Math. Comput. Appl.24, 70 (2019)

    2019

  11. [19]

    Armadillo: an efficient framework for numerical linear algebra,

    C. Sanderson and R. Curtin, “Armadillo: an efficient framework for numerical linear algebra,” in2025 17th International Conference on Computer and Automation Engineering (ICCAE), (IEEE, 2025), pp. 303–307

    2025

  12. [20]

    The design and implementation of FFTW3,

    M. Frigo and S. G. Johnson, “The design and implementation of FFTW3,” Proc. IEEE93, 216–231 (2005). Special issue on “Program Generation, Optimization, and Platform Adaptation”

    2005

  13. [21]

    Attosecond metrology,

    M. Hentschel, R. Kienberger, C. Spielmann,et al., “Attosecond metrology,” Nature414, 509–513 (2001)

    2001

  14. [22]

    Thegeneration,characterizationandapplicationsofbroadbandisolatedattosecond pulses,

    M.Chini,K.Zhao,andZ.Chang,“Thegeneration,characterizationandapplicationsofbroadbandisolatedattosecond pulses,” Nat. Photonics8, 178–186 (2014)

    2014

  15. [23]

    Femtosecond m2,3-edge spectroscopy of transition-metal oxides: Photoinduced oxidation state change in𝛼-fe2o3,

    J. Vura-Weis, C.-M. Jiang, C. Liu,et al., “Femtosecond m2,3-edge spectroscopy of transition-metal oxides: Photoinduced oxidation state change in𝛼-fe2o3,” The J. Phys. Chem. Lett.4, 3667–3671 (2013)

    2013

  16. [24]

    Femtosecond x-ray spectroscopy of an electrocyclic ring-opening reaction,

    A. R. Attar, A. Bhattacherjee, C. D. Pemmaraju,et al., “Femtosecond x-ray spectroscopy of an electrocyclic ring-opening reaction,” Science356, 54–59 (2017)

    2017

  17. [25]

    Ultrafast x-ray transient absorption spectroscopy of gas-phase photochemical reactions: A new universal probe of photoinduced molecular dynamics,

    A. Bhattacherjee and S. R. Leone, “Ultrafast x-ray transient absorption spectroscopy of gas-phase photochemical reactions: A new universal probe of photoinduced molecular dynamics,” Accounts Chem. Res.51, 3203–3211 (2018)

    2018

  18. [26]

    Femtosecond symmetry breaking and coherent relaxation of methane cations via x-ray spectroscopy,

    E. Ridente, D. Hait, E. A. Haugen,et al., “Femtosecond symmetry breaking and coherent relaxation of methane cations via x-ray spectroscopy,” Science380, 713–717 (2023)

    2023

  19. [27]

    High-order harmonic generation in a strongly overdriven regime,

    B. Major, K. Kovács, E. Svirplys,et al., “High-order harmonic generation in a strongly overdriven regime,” Phys. Rev. A107, 023514 (2023)

    2023

  20. [28]

    Experimental and theoretical analyses of a correlation between pump-pulse propagation and harmonic yield in a long-interaction medium,

    E. Takahashi, V. Tosa, Y. Nabekawa, and K. Midorikawa, “Experimental and theoretical analyses of a correlation between pump-pulse propagation and harmonic yield in a long-interaction medium,” Phys. Rev. A68(2003)

    2003

Pith tools

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