pith. sign in

arxiv: 2512.23156 · v2 · submitted 2025-12-29 · 🪐 quant-ph · physics.atom-ph

Emergence of nonclassical radiation in strongly laser-driven quantum systems

Ivan Gonoskov , Christian H\"unecke , Stefanie Gr\"afe This is my paper

Pith reviewed 2026-05-16 19:52 UTC · model grok-4.3

classification 🪐 quant-ph physics.atom-ph
keywords nonclassical radiationhigh-order harmonic generationPauli-Fierz HamiltonianWigner function negativitysqueezed statesstrong-field quantum opticsdipole response
0
0 comments X

The pith

Nonclassical radiation in laser-driven systems arises from nonlinear dependence of the electronic dipole on the light-mode coordinate.

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

The paper introduces an analytical framework for quantum-optical effects in strong laser-matter interactions, particularly high-order harmonic generation. It begins from the Pauli-Fierz Hamiltonian and applies a parametric factorization to separate the field-driven electronic dynamics from the quantized light mode. The central mechanism is the shape of the dipole response function: a flat response produces coherent light, a linear slope produces squeezing, and higher-order curvature produces states with negative Wigner-function values. The framework is applied to atomic and molecular models and extended to collections of emitters to show how photon number can be increased while preserving nonclassicality. This yields a direct, predictive link between classical strong-field trajectories and the quantum statistics of the emitted radiation.

Core claim

Within the parametric factorization, the emitted field inherits its quantum character directly from the functional form of the electronic dipole moment as a function of the light-mode coordinate. A constant dipole produces a coherent state, a linear dipole produces a squeezed state, and nonlinear dipole terms generate Wigner-function negativity, with the effect persisting and scaling in multi-emitter geometries.

What carries the argument

Parametric factorization of the light-matter wavefunction, which decouples the electronic state (driven by the quantized field coordinate) from the light-mode dynamics.

If this is right

  • Higher-order nonlinearities in the dipole response directly generate Wigner-function negativity in the radiated light.
  • The same mechanism produces squeezing when the dipole is linear in the light coordinate.
  • In systems with multiple emitters the nonclassical features scale to higher photon numbers while retaining the same origin.
  • The framework supplies an explicit mapping from strong-field electronic trajectories to the photon statistics of the emitted field.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Tailoring molecular structure or laser intensity could be used to select specific nonclassical features such as squeezing versus negativity.
  • The mechanism may extend to other strong-field processes such as above-threshold ionization or laser-assisted electron scattering.
  • Experimental detection could focus on homodyne measurements of the harmonic field phase noise in single-atom or low-density targets.

Load-bearing premise

The parametric factorization of the coupled light-matter wavefunction remains accurate enough to capture the quantum-optical properties of the emitted field across the intensity range of interest.

What would settle it

Full numerical solution of the Pauli-Fierz Hamiltonian for a driven atom at intensities where the factorization predicts Wigner negativity, followed by direct reconstruction of the emitted field's Wigner function, would falsify the claim if no negativity appears.

Figures

Figures reproduced from arXiv: 2512.23156 by Christian H\"unecke, Ivan Gonoskov, Stefanie Gr\"afe.

Figure 1
Figure 1. Figure 1: Visualization of the emergence of nonclassicality in high-order harmonic generation to result from the nonlinear q- dependence of the oscillating dipole moment. The upper row displays the dipole moments expanded into a series (truncated after the terms as labeled) (eq. 9) while the lower row shows the resulting output state of light in phase space, represented by its Wigner function. (a) Constant, almost q… view at source ↗
Figure 2
Figure 2. Figure 2: Gallery of Wigner functions after the interaction with a laser pulse for different systems and initial conditions. (a, d) atomic system mimicking the Ca atom; (b, e) molecular system with fixed internuclear distance; (c, f) hydrogen atom-like potential. The upper row shows the results when starting with the electronic ground state of the driven system while the lower row the corresponding results when star… view at source ↗
Figure 3
Figure 3. Figure 3: Nonclassical output states for Ω = 5th harmonic generated from the multi-emitter case with Ne spatially isolated (non-Coulomb interacting) emitters. Each emitter is a 1D model CaO model molecule). The parameters of the driving laser pulse are Fclass = 3.4 · 10−2 (a.u.), I = 4 · 1013 W/cm2 , λ = 1945nm. The solid line corresponds to the solution based on the parametric factorization. The dashed line shows t… view at source ↗
Figure 4
Figure 4. Figure 4: Nonclassical harmonic generation from various model systems discussed in the main text: 1D electronic system coupled to one mode of quantized light with harmonic frequency Ω as labeled. (A) Left: 1D soft-core model potential (blue) mimicking the Ca-atom and corresponding electron wavefunctions of the ground state (orange) and a superposition of the ground and excited electronic state (Ca atom: Ip = 6.11 eV… view at source ↗
Figure 5
Figure 5. Figure 5: As in [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Resonant part of dipole moment (eq. 24) as a function of q, corresponding to the numerical solutions in a field Fclass(t) +β · q · sin(Ωt); as described in the text, polynomial fitting this form is performed, which is used for the analytical solving the φ light(q,β,tint). The calculations were for a model hydrogen atom with the following parameters: Ω = 13, λ = 2227 nm, Tp = 4 cycles, I = 4 · 1013 W/cm2 , … view at source ↗
read the original abstract

We present an analytical framework for the emergence of nonclassical radiation in strongly laser-driven quantum systems, with a focus on high-order harmonic generation (HHG). Starting from a Pauli-Fierz description, we employ a parametric factorization of the coupled light-matter wavefunction that reduces the dynamics to coupled equations for a field-driven electronic state and a quantized light mode. Within this framework, we identify a simple and predictive mechanism for nonclassicality: it originates from the nonlinear dependence of the electronic dipole response on the light-mode coordinate. An approximately constant dipole yields coherent radiation, a linear dependence produces squeezing, and higher-order nonlinearities give rise to Wigner-function negativity. We illustrate this mechanism for atomic and molecular model systems and analyze its scaling in multi-emitter configurations, indicating routes toward high-photon-number nonclassical radiation in HHG. Our results provide a transparent connection between strong-field dynamics and quantum-optical properties of emitted light, offering a basis for engineering nonclassical states in intense laser-matter interactions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops an analytical framework starting from the Pauli-Fierz Hamiltonian and employing a parametric factorization of the light-matter wavefunction. This reduces the problem to coupled equations for a field-driven electronic state and a quantized light mode. The central claim is that nonclassical radiation emerges from the nonlinear dependence of the electronic dipole response on the light-mode coordinate: constant dipole yields coherent states, linear dependence produces squeezing, and higher-order nonlinearities generate Wigner-function negativity. The mechanism is illustrated on atomic and molecular model systems in the HHG context and extended to multi-emitter scaling for high-photon-number nonclassical states.

Significance. If the factorization ansatz remains accurate, the work supplies a transparent, parameter-free connection between strong-field electronic dynamics and the quantum statistics of emitted light. This offers a predictive route to engineering nonclassical radiation in intense laser-matter interactions without requiring full numerical solution of the entangled light-matter state.

major comments (2)
  1. [Theory development / parametric factorization ansatz] The central claim that higher-order nonlinearities in the dipole response d(x) produce Wigner negativity rests on the parametric factorization remaining valid. In the HHG-relevant regime the electronic wavefunction involves ionization, continuum propagation, and recollision, processes that generically generate photon-number entanglement. No direct benchmark is supplied that compares the approximate reduced field Wigner function obtained from the factorization against the exact reduced density matrix for the same Hamiltonian and model systems.
  2. [Multi-emitter scaling section] The scaling analysis for multi-emitter configurations assumes the single-emitter factorization carries over without additional entanglement channels. Because the paper reports negativity only within the factored ansatz, it is unclear whether the predicted high-photon-number nonclassicality survives once inter-emitter and light-matter entanglement are restored.
minor comments (2)
  1. [Introduction / Theory] Notation for the light-mode coordinate and the dipole operator should be introduced with an explicit equation number at first use to improve readability.
  2. [Results / Figures] Figure captions for the Wigner-function plots should state the intensity range and the precise model Hamiltonian used, including any truncation of the electronic basis.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and constructive feedback on our manuscript. We appreciate the recognition of the potential significance of our analytical framework. Below, we provide point-by-point responses to the major comments and outline the revisions we will make to address them.

read point-by-point responses

  1. Referee: The central claim that higher-order nonlinearities in the dipole response d(x) produce Wigner negativity rests on the parametric factorization remaining valid. In the HHG-relevant regime the electronic wavefunction involves ionization, continuum propagation, and recollision, processes that generically generate photon-number entanglement. No direct benchmark is supplied that compares the approximate reduced field Wigner function obtained from the factorization against the exact reduced density matrix for the same Hamiltonian and model systems.

    Authors: We agree that a direct benchmark against the exact reduced density matrix would strengthen the validation of the ansatz, especially given the entanglement-generating processes in HHG. Such comparisons are numerically intensive due to the large Hilbert space. Our derivation starts from the Pauli-Fierz Hamiltonian and the factorization is exact when the electronic state is independent of the field coordinate in a certain sense, but we acknowledge the approximation in strong-field regimes. In the revised manuscript, we will include a more detailed discussion of the ansatz's validity conditions and add a benchmark for a simplified two-level system where exact comparison is feasible. revision: partial

  2. Referee: The scaling analysis for multi-emitter configurations assumes the single-emitter factorization carries over without additional entanglement channels. Because the paper reports negativity only within the factored ansatz, it is unclear whether the predicted high-photon-number nonclassicality survives once inter-emitter and light-matter entanglement are restored.

    Authors: The multi-emitter section is meant to illustrate the potential for scaling up the nonclassicality within the framework of the factorization ansatz. We recognize that restoring full entanglement could modify the quantitative predictions. We will revise the text to emphasize that this is an indicative scaling analysis under the ansatz and discuss possible effects of inter-emitter entanglement as an avenue for future work. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained within parametric ansatz

full rationale

The paper introduces the parametric factorization ansatz explicitly as an approximation that reduces the Pauli-Fierz dynamics to coupled equations for the field-driven electronic state and quantized mode. Within this framework it expands the dipole response d(x) in powers of the light-mode coordinate x and maps constant/linear/higher-order terms to coherent/squeezed/negative-Wigner radiation. No quoted step defines the nonclassicality measure in terms of a fitted parameter that is then called a prediction, nor does any load-bearing claim reduce to a self-citation chain or imported uniqueness theorem. The central mapping is a direct algebraic consequence of the chosen factorization and the dipole Taylor expansion; the result is therefore independent of the specific numerical values chosen for the electronic dynamics.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on the validity of the parametric factorization ansatz and on the truncation of the dipole expansion at low orders; no free parameters are introduced in the abstract-level description, and no new entities are postulated.

axioms (2)
  • domain assumption The Pauli-Fierz Hamiltonian provides an accurate starting point for the coupled light-matter system.
    Invoked at the opening of the framework description.
  • domain assumption The parametric factorization of the wavefunction captures the essential dynamics without significant back-action errors.
    Central modeling step that reduces the problem to coupled equations.

pith-pipeline@v0.9.0 · 5471 in / 1350 out tokens · 26298 ms · 2026-05-16T19:52:09.229170+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages

  1. [1]

    D. N. Klyshko, Phys.-Usp.39, 573, (1996)

    work page 1996

  2. [2]

    P. B. Corkum and F. Krausz, Attosecond science, Nat. Phys.3, 381 (2007)

    work page 2007

  3. [3]

    Chekhova, G

    M.V . Chekhova, G. Leuchs, M. ˙Zukowski, Bright squeezed vacuum: Entanglement of macroscopic light beams, Opt Commun, 337 (2015) 27-43

    work page 2015

  4. [4]

    Zavatta et al., Quantum-to-Classical Transition with Single-Photon-Added Coherent States of Light, Science306, 660-662 (2004)

    A. Zavatta et al., Quantum-to-Classical Transition with Single-Photon-Added Coherent States of Light, Science306, 660-662 (2004)

    work page 2004

  5. [5]

    Ourjoumtsev et al., Generating Optical Schrödinger Kittens for Quantum Information Processing, Science312, 83-86 (2006)

    A. Ourjoumtsev et al., Generating Optical Schrödinger Kittens for Quantum Information Processing, Science312, 83-86 (2006)

    work page 2006

  6. [6]

    Mikheev et al., Efficient production of large-size optical Schrödinger cat states, Sci

    E.V . Mikheev et al., Efficient production of large-size optical Schrödinger cat states, Sci. Rep.9, 14301 (2019)

    work page 2019

  7. [7]

    Gonoskov, N

    I.A. Gonoskov, N. Tsatrafyllis, I.K. Kominis, P. Tzallas, Quantum optical signatures in strong-field laser physics: Infrared photon counting in high-order-harmonic generation, Sci Rep, 6 (2016) 32821

    work page 2016

  8. [8]

    Tsatrafyllis, I.K

    N. Tsatrafyllis, I.K. Kominis, I.A. Gonoskov, P. Tzallas, High-order harmonics measured by the photon statistics of the infrared driving-field exiting the atomic medium, Nat Commun, 8 (2017) 15170

    work page 2017

  9. [9]

    Lewenstein, M.F

    M. Lewenstein, M.F. Ciappina, E. Pisanty, J. Rivera-Dean, P. Stammer, T. Lamprou, P. Tzallas, Generation of optical Schrödinger cat states in intense laser–matter interactions, Nat Phys, 17 (2021) 1104-1108

    work page 2021

  10. [10]

    Theidel, V

    D. Theidel, V . Cotte, R. Sondenheimer, V . Shiriaeva, M. Froidevaux, V . Severin, A. Merdji-Larue, P. Mosel, S. Fröhlich, K.-A. Weber, U. Morgner, M. Kovacev, J. Biegert, H. Merdji, Evidence of the Quantum Optical Nature of High-Harmonic Generation, Prx Quantum, 5 (2024) 040319

    work page 2024

  11. [11]

    Rivera-Dean, P

    J. Rivera-Dean, P. Stammer, A.S. Maxwell, T. Lamprou, A.F. Ordóñez, E. Pisanty, P. Tzallas, M. Lewenstein, M.F. Ciappina, Nonclassical states of light after high-harmonic generation in semiconductors: A Bloch-based perspective, Phys Rev B, 109 (2024) 035203

    work page 2024

  12. [12]

    Lamprou, P

    Th. Lamprou, P. Stammer, J. Rivera-Dean, N. Tsatrafyllis, M. F. Ciappina, M. Lewenstein, P. Tzallas,Recent developments in the generation of non-classical and entangled light states using intense laser-matter interactions, arXiv:2410.17452, (2024)

    work page arXiv 2024

  13. [13]

    Stammer, J

    P. Stammer, J. Rivera-Dean, A.S. Maxwell, T. Lamprou, J. Argueello-Luengo, P. Tzallas, M.F. Ciappina, M. Lewenstein, Entanglement and Squeezing of the Optical Field Modes in High Harmonic Generation, Phys Rev Lett, 132 (2024) 143603

    work page 2024

  14. [14]

    Philipp Stammer, Javier Rivera-Dean, Maciej Lewenstein,Theory of quantum optics and optical coherence in high har- monic generation, arXiv:2504.13287, (2025)

    work page arXiv 2025

  15. [15]

    Stammer, J

    P. Stammer, J. Rivera-Dean, A. Maxwell, T. Lamprou, A. Ordóñez, M.F. Ciappina, P. Tzallas, M. Lewenstein, Quantum Electrodynamics of Intense Laser-Matter Interactions: A Tool for Quantum State Engineering, Prx Quantum, 4 (2023) 010201

    work page 2023

  16. [16]

    Gonoskov, S

    I. Gonoskov, S. Gräfe, Light-matter quantum dynamics of complex laser-driven systems, J Chem Phys, 154 (2021) 234106

    work page 2021

  17. [17]

    S. Yi, N.D. Klimkin, G.G. Brown, O. Smirnova, S. Patchkovskii, I. Babushkin, M. Ivanov, Generation of Massively Entangled Bright States of Light during Harmonic Generation in Resonant Media, Phys Rev X, 15 (2025) 011023

    work page 2025

  18. [18]

    Christian Saugbjerg Lange, Thomas Hansen, Lars Bojer Madsen,Electron-correlation induced nonclassicallity of light from high-harmonic generation, Phys. Rev. A 109, 033110, (2024)

    work page 2024

  19. [19]

    Gorlach, M.E

    A. Gorlach, M.E. Tzur, M. Birk, M. Krüger, N. Rivera, O. Cohen, I. Kaminer, High-harmonic generation driven by quantum light, Nat Phys, 19 (2023) 1689-1696

    work page 2023

  20. [20]

    Fang, F.X

    Y . Fang, F.X. Sun, Q. He, Y . Liu, Strong-Field Ionization of Hydrogen Atoms with Quantum Light, Phys Rev Lett, 130 (2023) 253201

    work page 2023

  21. [21]

    Even Tzur, O

    M. Even Tzur, O. Cohen, Motion of charged particles in bright squeezed vacuum, Light Sci Appl, 13 (2024) 41

    work page 2024

  22. [22]

    Even Tzur, M

    M. Even Tzur, M. Birk, A. Gorlach, M. Krüger, I. Kaminer, O. Cohen, Photon-statistics force in ultrafast electron dynam- ics, Nat Photonics, 17 (2023) 501-509

    work page 2023

  23. [23]

    Stammer, et al.High Photon Number Entangled States and Coherent State Superposition from the Extreme Ultraviolet to the Far Infrared, Phys

    P. Stammer, et al.High Photon Number Entangled States and Coherent State Superposition from the Extreme Ultraviolet to the Far Infrared, Phys. Rev. Lett. 128, 123603, (2022)

    work page 2022

  24. [24]

    I. A. Gonoskov, et al.,Nonclassical light generation and control from laser-driven semiconductor intraband excitations, Phys. Rev. B, 109, 125110 (2024)

    work page 2024

  25. [25]

    and Gross, E

    Abedi, Ali and Maitra, Neepa T. and Gross, E. K. U.,Exact Factorization of the Time-Dependent Electron-Nuclear Wave Function, Phys. Rev. Lett.,105, 123002, (2010)

    work page 2010

  26. [26]

    Härkönen, and Ivan A

    Ville J. Härkönen, and Ivan A. Gonoskov,On the diagonalization of quadratic Hamiltonians, Journal of Physics A: Mathematical and Theoretical 55 (1), 015306 (2022). 7

    work page 2022

  27. [27]

    Lachman, L

    L. Lachman, L. Slodi ˇcka, and R. Filip,Nonclassical light from a large number of independent single-photon emitters, Sci. Rep. 6, 19760 (2016)

    work page 2016

  28. [28]

    Christian Saugbjerg Lange, Lars Bojer Madsen,Hierarchy of approximations for describing quantum light from high- harmonic generation: A Fermi-Hubbard model study, Phys. Rev. A 111, 013113, (2025)

    work page 2025

  29. [29]

    Fleck, J

    J. Fleck, J. Morris, and M. Feit,Time-Dependent Propagation of High Energy Laser Beams through the Atmosphere, Applied Physics A: Materials Science and Processing, V ol. 10, No. 2, pp. 129-160, (1976)

    work page 1976

  30. [30]

    Gonoskov I.A., Vugalter G.A., Mironov V .A.,Ionization in a Quantized Electromagnetic Field, J. Exp. Theor. Phys.105, 1119, 5 (2007). METHODS/APPENDIX Derivation We start from the time-dependent Schrödinger equation for the fully quantum system in the Schödinger repre- sentation within the dipole approximation for all light modes, see also [16] and refere...

    work page 2007