Ab-initio calculations of laser-atom interactions reveal harmonics feedback during macroscopic propagation
Pith reviewed 2026-05-25 19:28 UTC · model grok-4.3
The pith
Generated harmonics feed back to alter ionization probability and harmonic yields during macroscopic laser propagation.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By linking the full 3D ab initio quantum evolution of the light pulse polarization in interaction with an atom with a propagation model, the calculations simulate the propagation of ultrashort laser pulses over macroscopic dimensions in the presence of self-generated harmonics up to order 11 and evidence a clear feedback of the generated harmonics on propagation, with an influence on the ionization probability as well as the yield of the harmonic generation itself.
What carries the argument
Coupled 3D ab initio quantum polarization evolution and macroscopic propagation model
If this is right
- The feedback changes the ionization probability along the propagation path.
- The feedback modifies the yield of harmonic generation itself.
- Simulations of macroscopic pulse propagation must retain the harmonic components to capture this effect.
- The influence appears for harmonics up to at least the eleventh order.
Where Pith is reading between the lines
- Models of high-harmonic generation sources that treat propagation separately from atomic response may need revision to include this loop.
- Filamentation or self-focusing experiments in gases could exhibit measurable signatures of the same feedback.
- The coupled approach could be tested by varying gas pressure or pulse energy to map how the feedback strength scales.
Load-bearing premise
The numerical coupling between the 3D ab initio atomic polarization evolution and the macroscopic propagation model faithfully represents the physics without unaccounted approximations in either component.
What would settle it
A side-by-side comparison of measured ionization rates or harmonic yields in a gas cell with versus without the generated harmonics present, showing divergence that matches the difference between the coupled simulation and an uncoupled propagation run.
Figures
read the original abstract
We couple the full 3D ab initio quantum evolution of the light pulse polarization in interaction with an atom with a propagation model to simulate the propagation of ultrashort laser pulses over macroscopic dimensions, in the presence of self-generated harmonics up to order 11. We evidence a clear feedback of the generated harmonics on propagation, with an influence on the ionization probability as well as the yield of the harmonic generation itself.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript couples a 3D ab-initio TDSE treatment of the atomic polarization response to a macroscopic propagation model for ultrashort laser pulses. It reports that self-generated harmonics up to order 11 produce a measurable feedback that alters both the ionization probability and the HHG yield during propagation.
Significance. If the numerical coupling is shown to be free of discretization or split-step artifacts and the feedback survives convergence tests, the result would strengthen the case for self-consistent inclusion of harmonic back-action in macroscopic HHG simulations, affecting quantitative predictions of ionization and conversion efficiency.
major comments (2)
- [Numerical method / propagation coupling] The description of the iterative coupling between the TDSE polarization and the propagation operator (likely in the Methods or Numerical Implementation section) must explicitly state whether the atomic response is recomputed at each spatial step using the full field (fundamental plus harmonics up to order 11) or only the fundamental; the former is required to claim genuine nonlinear feedback rather than linear superposition.
- [Results / feedback quantification] Convergence of the reported feedback with respect to propagation step size, transverse grid spacing, and harmonic bandwidth cutoff must be demonstrated; without such tests the influence on ionization probability and HHG yield could arise from inconsistent discretization or artificial dispersion in the envelope equation.
minor comments (1)
- [Abstract] The abstract states the central finding but supplies no quantitative measure (e.g., percentage change in ionization or yield) or error estimate; a brief numerical illustration should be added.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable suggestions. Below we address each major comment in detail.
read point-by-point responses
-
Referee: [Numerical method / propagation coupling] The description of the iterative coupling between the TDSE polarization and the propagation operator (likely in the Methods or Numerical Implementation section) must explicitly state whether the atomic response is recomputed at each spatial step using the full field (fundamental plus harmonics up to order 11) or only the fundamental; the former is required to claim genuine nonlinear feedback rather than linear superposition.
Authors: The coupling in our simulations recomputes the atomic polarization response using the complete electric field at each propagation step, incorporating the fundamental frequency as well as the harmonics up to the 11th order. This approach is what enables the observation of the nonlinear feedback effect. We agree that the description in the manuscript could be more explicit on this point and will revise the Numerical Implementation section to clearly state that the full field is used for the TDSE solution at every spatial step. revision: yes
-
Referee: [Results / feedback quantification] Convergence of the reported feedback with respect to propagation step size, transverse grid spacing, and harmonic bandwidth cutoff must be demonstrated; without such tests the influence on ionization probability and HHG yield could arise from inconsistent discretization or artificial dispersion in the envelope equation.
Authors: We recognize the importance of demonstrating numerical convergence to rule out artifacts. While our primary results were obtained with parameters that we believe are sufficient, we did not include explicit convergence tests for all mentioned quantities in the original submission. We will perform additional convergence studies and include them in the revised manuscript to confirm that the feedback effects are robust. revision: yes
Circularity Check
No circularity: simulation output from coupled TDSE-propagation model
full rationale
The paper reports results from numerically coupling a 3D ab-initio TDSE atomic response to a macroscopic propagation equation that includes harmonics up to order 11. The claimed feedback on ionization probability and HHG yield is an emergent output of this direct simulation rather than any algebraic reduction, parameter fit, or self-referential definition. No equations are presented that equate a derived quantity to its own input by construction, and the approach contains no load-bearing self-citations or ansatz smuggling. The derivation chain is therefore self-contained as a numerical experiment.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The 3D ab initio quantum evolution of the light pulse polarization accurately captures the atomic response under the simulated conditions.
Reference graph
Works this paper leans on
-
[1]
593 a.u. the hydrogen polarizability at 800 nm) we can write t he previous equation as: ∂x ˜Az = i ( ωc −ω v ) ˜Az + ˜/bbPz + iωǫ0ρatα0 ˜Az 2ǫ0c , (9) then ∂x ˜Az = i [ ω c ( 1 + ρatα0 2 ) −ω v ] ˜Az + 1 2ǫ0c ˜/bbPz. (10) 3 11 ω / ω0 10 8 10 10 10 12 10 14 10 16 Spectrum (arb. units) Ly α Ly β Ly γ Ly δ 1 3 5 7 9 Figure 1: Electric field spectrum of a 50 T...
work page 2048
-
[2]
In contrast, harmonic 9 reaches a peak af ter as little as 0
5 mm, and then decay. In contrast, harmonic 9 reaches a peak af ter as little as 0 . 2 mm, then decays almost fully till 3 mm, before entering in a new growth cycle. The faster rise and dec ay of the 7 th, and even more of the 9 th harmonics, can be expected to stem from their partial overlap with the Lymann resonances, that strongly affect both their abso...
-
[3]
5 mm. We attribute the complex evolution of the harmonics pul se shape to their wide band associated with their short durat ion. As a consequence, group-velocity dispersion is significant within each harmonic. The phase matching and the associated (back-) conversion efficiency are therefore di fferent among spectral slices within each harmonic. To get more in...
-
[4]
X. F. Li, A. L’Huillier, M. Ferray, L. A. Lompré, and G. Mai nfray, Physical Review A 39, 5751 (1989)
work page 1989
- [5]
-
[6]
P . Agostini and L. F. DiMauro, Reports on progress in phys ics 67, 813 (2004)
work page 2004
- [7]
-
[8]
S. L. Chin, S. A. Hosseini, W. Liu, Q. Luo, F. Théberge, N. A közbek, A. Becker, V . P . Kandidov, O. G. Kosareva, and H. Schr öder, Canadian Journal of Physics 83, 863 (2005)
work page 2005
- [9]
- [10]
-
[11]
Wolf, Reports on Progress in Physics 81, 026001 (2017)
J.-P . Wolf, Reports on Progress in Physics 81, 026001 (2017)
work page 2017
-
[12]
P . B. Corkum, Physical Review Letters 71, 1994 (1993)
work page 1994
-
[13]
K. J. Schafer, B. Yang, L. F. DiMauro, and K. C. Kulander, Physical review letters 70, 1599 (1993)
work page 1993
-
[14]
A. Rundquist, C. G. Durfee, Z. Chang, C. Herne, S. Backus , M. M. Murnane, and H. C. Kapteyn, Science 280, 1412 (1998)
work page 1998
-
[15]
T. Popmintchev, M.-C. Chen, A. Bahabad, M. Gerrity, P . S idorenko, O. Cohen, I. P . Christov, M. M. Murnane, and H. C. Ka pteyn, Proceedings of the National Academy of Sciences 106, 10516 (2009)
work page 2009
-
[16]
M. Lewenstein, P . Balcou, M. Y . Ivanov, A. L’Huillier, a nd P . B. Corkum, Physical Review A 49, 2117 (1994)
work page 1994
- [17]
-
[18]
M. B. Gaarde, J. L. Tate, and K. J. Schafer, Journal of Phy sics B: Atomic, Molecular and Optical Physics 41, 132001 (2008)
work page 2008
-
[19]
A.-T. Le, R. Lucchese, S. Tonzani, T. Morishita, and C. L in, Physical Review A 80, 013401 (2009)
work page 2009
- [20]
-
[21]
M. Kolesik and J. V . Moloney, Reports on Progress in Phys ics 77, 016401 (2013)
work page 2013
-
[22]
A. Bahl, V . P . Majety, A. Scrinzi, and M. Kolesik, Optics Letters 42, 2295 (2017), ISSN 0146-9592, 1539-4794, URL https://www.osapublishing.org/abstract.cfm?URI=ol-42-12-2295
work page 2017
-
[23]
Y . Liu, M. Durand, A. Houard, B. Forestier, A. Couairon, and A. Mysyrowicz, Optics Communications 284, 4706 (2011)
work page 2011
-
[24]
M. B. Gaarde and A. Couairon, Physical review letters 103, 043901 (2009)
work page 2009
- [25]
-
[26]
N. Aközbek, A. Iwasaki, A. Becker, M. Scalora, S. L. Chin , and C. M. Bowden, Physical review letters 89, 143901 (2002)
work page 2002
-
[27]
T. V ockerodt, D. S. Steingrube, E. Schulz, M. Kretschma r, U. Morgner, and M. Kova ˇcev, Applied Physics B 106, 529 (2012)
work page 2012
-
[28]
M. Kolesik, E. M. Wright, A. Becker, and J. V . Moloney, Ap plied Physics B 85, 531 (2006)
work page 2006
-
[29]
D. S. Steingrube, E. Schulz, T. Binhammer, M. B. Gaarde, A. Couairon, U. Morgner, and M. Kovaˇcev, New Journal of Physics 13, 043022 (2011)
work page 2011
- [30]
- [31]
-
[32]
J. Doussot, P . Béjot, and O. Faucher, Physical Review A 93, 033857 (2016)
work page 2016
- [33]
- [34]
- [35]
-
[36]
A. Paul, R. A. Bartels, R. Tobey, H. Green, S. Weiman, I. P . Christov, M. M. Murnane, H. C. Kapteyn, and S. Backus, Natur e 421, 51 (2003)
work page 2003
-
[37]
E. A. Gibson, A. Paul, N. Wagner, T. Ra’anan, D. Gaudiosi , S. Backus, I. P . Christov, A. Aquila, E. M. Gullikson, D. T. A ttwood, et al., Science 302, 95 (2003)
work page 2003
-
[38]
D. L. Weerawarne, X. Gao, A. L. Gaeta, and B. Shim, Physic al review letters 114, 093901 (2015)
work page 2015
- [39]
- [40]
- [41]
-
[42]
M. Petrarca, Y . Petit, S. Henin, R. Delagrange, P . Béjot, and J. Kasparian, Optics letters 37, 4347 (2012)
work page 2012
-
[43]
M. Richter, S. Patchkovskii, F. Morales, O. Smirnova, a nd M. Ivanov, New Journal of Physics 15, 083012 (2013)
work page 2013
-
[44]
J. H. Odhner, D. A. Romanov, E. T. McCole, J. K. Wahlstran d, H. M. Milchberg, and R. J. Levis, Physical review letters 109, 065003 (2012)
work page 2012
-
[45]
J. K. Wahlstrand, Y .-H. Cheng, and H. M. Milchberg, Phys ical Review A 85, 043820 (2012)
work page 2012
- [46]
-
[47]
O. Kosareva, J.-F. Daigle, N. Panov, T. Wang, S. Hossein i, S. Y uan, G. Roy, V . Makarov, and S. L. Chin, Optics letters36, 1035 (2011)
work page 2011
-
[48]
A. V . Mitrofanov, A. A. V oronin, D. A. Sidorov-Biryukov , A. Pugžlys, E. A. Stepanov, G. Andriukaitis, T. Flöry, S. Ališauskas, A. B. Fedotov, A. Baltuška, et al., Scientific Reports 5, 8368 (2015), ISSN 2045-2322, URL 9 https://www.nature.com/articles/srep08368
work page 2015
-
[49]
T. C. Rensink, T. M. Antonsen Jr, J. P . Palastro, and D. F. Gordon, Physical Review A 89, 033418 (2014)
work page 2014
-
[50]
E. V olkova, A. Popov, and O. Tikhonova, Journal of Exper imental and Theoretical Physics 116, 372 (2013)
work page 2013
- [51]
- [52]
-
[53]
A. D. Bandrauk and H. Shen, The Journal of chemical physi cs 99, 1185 (1993)
work page 1993
-
[54]
A. D. Bandrauk, S. Chelkowski, D. J. Diestler, J. Manz, a nd K.-J. Y uan, Physical Review A 79, 023403 (2009)
work page 2009
-
[55]
M. Nurhuda, A. Suda, and K. Midorikawa, Physical Review A 66, 041802 (2002)
work page 2002
-
[56]
M. Nurhuda, A. Suda, and K. Midorikawa, New Journal of Ph ysics 10, 053006 (2008)
work page 2008
-
[57]
P . O. Kano, M. Brio, and J. V . Moloney, Communications in Mathematical Sciences 4, 53 (2006)
work page 2006
- [58]
- [59]
-
[60]
E. Cormier, D. Garzella, P . Breger, P . Agostini, G. Chér iaux, and C. Leblanc, Journal of Physics B: Atomic, Molecula r and Optical Physics 34, L9 (2001)
work page 2001
-
[61]
M. Krüger and N. Dudovich, Nature Photonics 10, 626 (2016), ISSN 1749-4885, 1749-4893, URL http://www.nature.com/articles/nphoton.2016.189
work page 2016
-
[62]
G. Orenstein, O. Pedatzur, A. J. Uzan, B. D. Bruner, Y . Ma iresse, and N. Dudovich, Physical Review A 95 (2017), ISSN 2469-9926, 2469-9934, URL http://link.aps.org/doi/10.1103/PhysRevA.95.051401
-
[63]
C. Hofmann, A. Landsman, and U. Keller, Applied Science s 8, 341 (2018), ISSN 2076-3417, URL http://www.mdpi.com/2076-3417/8/3/341
work page 2018
- [64]
- [65]
- [66]
- [67]
- [68]
-
[69]
M. Kolesik, J. V . Moloney, and M. Mlejnek, Physical revi ew letters 89, 283902 (2002)
work page 2002
- [70]
-
[71]
C. Cohen-Tannoudji, B. Diu, and F. Laloë, Mécanique quantique. Tome I (Hermann, 1973)
work page 1973
-
[72]
G. L. Sleijpen and D. R. Fokkema, Electronic Transactio ns on Numerical Analysis 534, 86 (1993)
work page 1993
-
[73]
G. L. Sleijpen, H. A. V an der V orst, and D. R. Fokkema, Num erical Algorithms 7, 75 (1994)
work page 1994
- [74]
-
[75]
E. Cormier and P . Lambropoulos, Journal of Physics B: At omic, Molecular and Optical Physics 30, 77 (1997)
work page 1997
-
[76]
J. Weideman and B. Herbst, SIAM Journal on Numerical Ana lysis 23, 485 (1986)
work page 1986
-
[77]
M. Karplus and H. Kolker, The Journal of Chemical Physic s 39, 1493 (1963)
work page 1963
-
[78]
J. Doussot, P . Béjot, G. Karras, F. Billard, and O. Fauch er, Journal of Physics B: Atomic, Molecular and Optical Phys ics 48, 184005 (2015)
work page 2015
- [79]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.