Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2601.17164

Highly accurate semiclassical strong-field Herman-Kluk propagator method for high-harmonic generation

Phi-Hung Tran, Hao Quan Truong, R. Esteban Goetz, Anh-Thu Le

physics.atom-phquant-phSubmitted Jan 23, 2026• Announced Jan 27, 2026
arXivHTMLPDF

Abstract

We extend our recently developed semiclassical strong-field Herman-Kluk propagator (SFHK) method to calculate high-order harmonic generation (HHG) for atoms in intense lasers. We show that our method, based on a combination of the Herman-Kluk propagator and the strong-field approximation, can provide highly accurate results for both HHG yield and phase, nearly identical to those from the exact numerical solutions of the time-dependent Schrödinger equation. We provide detailed analyses of our method and its applications to the HHG process, particularly the recombination time. The main computational task in this method is to solve the classical Newton equations for the active electron in the combined atomic potential and laser-electron interaction. The motion of the centers of the electron wave packets, modeled by coherent states, is governed by independent classical trajectories so that the computation can therefore be parallelized very efficiently.

Comments
8 figures

Related papers

Strong-field Herman-Kluk propagator method for high-harmonic generation in molecules

Phi-Hung Tran, Hao Quan Truong, R. Esteban Goetz +1 more

We extend our recently developed semiclassical strong-field Herman-Kluk (SFHK) propagator method to calculate high-order harmonic generation (HHG) in diatomic molecules driven by few-cycle intense laser fields. On the example of applications to H2 and N2, we show that our method, based on a combination of the Herman-Kluk propagator and the strong-field approximation, can provide very accurate results for both HHG yield and phase, nearly identical to those from the exact numerical solutions of the time-dependent Schrodinger equation. To compare with experimental measurements, averaging over molecular orientations must be performed. Here we demonstrate a distinct and powerful advantage of the SFHK, as its Monte Carlo sampling for the integration over the alignment distribution can be efficiently combined with the integration over the initial momentum distributions of electron wave-packet right after the tunnel exit. Therefore, the total number of trajectories used for the alignment-averaged HHG spectrum does not increase much compared to that for a single fixed alignment. Similar to atomic targets, the main computational task in the SFHK is to solve the classical Hamiltonian equations for the active electron in the combined electron-target ion potential and electron-laser interaction. The motion of the center of each electron wave packet in the continuum, represented by a coherent state, is governed by an independent classical trajectory so that the computation can be parallelized very efficiently.

physics.atom-phCross-listedarXiv:2608.13686•Aug 17, 2026

Efficient Algorithms for Semiclassical Quantum Dynamics

David Sattlegger

We study the well-known Herman-Kluk propagator in order to calculate approximate solutions to the time--dependent semiclassical Schroedinger equation in high dimensions. We derive a discretisation scheme and provide its approximation properties. This scheme, as well as other semiclassical methods, are of an intrinsically parallel nature. This allows for the design of highly efficient algorithms employing state of the art parallelization and vectorization techniques.

•Jan 1, 2016•0 citations

Comment on Semiclassical approximations in phase space with coherent states

Frank Großmann, Michael F. Herman

The semiclassical Herman–Kluk initial value approximation of the quantum-mechanical propagator, which is heavily used in dynamical studies of atomic and molecular systems, takes into account nonlinearity of the underlying classical dynamics by using multiple classical initial value solutions. In applications to the propagation of Gaussian wavepackets non-Gaussian distortions in the course of time and therefore a realistic description of the quantum dynamics are possible. This is in contrast to the single-trajectory approximations investigated in a recent paper by Baranger et al , one of which was erroneously termed Herman–Kluk approximation.

Journal of Physics A Mathematical and General•Oct 21, 2002•38 citations

Quantum optical signatures in strong-field laser physics: Infrared photon counting in high-order-harmonic generation

Ivan Gonoskov, Nikolaos Tsatrafyllis, Iannis K. Kominis +1 more

We analytically describe the strong-field light-electron interaction using a quantized coherent laser state with arbitrary photon number. We obtain a light-electron wave function which is a closed-form solution of the time-dependent Schrödinger equation (TDSE). This wave function provides information about the quantum optical features of the interaction not accessible by semi-classical theories. With this approach we can reveal the quantum optical properties of high harmonic generation (HHG) process in gases by measuring the photon statistics of the transmitted infrared (IR) laser radiation. This work can lead to novel experiments in high-resolution spectroscopy in extreme-ultraviolet (XUV) and attosecond science without the need to measure the XUV light, while it can pave the way for the development of intense non-classical light sources.

arXiv:1601.02764•Scientific Reports•Sep 7, 2016•51 citations

High-Order Harmonic Generation with Beyond-Semiclassical Emitter Dynamics: A Strong-Field Quantum Optical Heisenberg Picture Approach

Christian Saugbjerg Lange, Ella Elisabeth Lassen, Rasmus Vesterager Gothelf +1 more

Quantum-optical descriptions of strong-field processes have attracted significant attention in recent years. Typically, the theoretical modeling has been conducted in the Schrödinger picture, where results are only obtainable under certain approximations, while, in contrast, the Heisenberg picture has remained relatively unexplored. In this work, we develop an accurately controlled perturbative expansion of the time-evolution operator in the Heisenberg picture and derive beyond-semiclassical corrections to the emitter dynamics due to the coupling to the quantized electromagnetic field, capturing effects of the quantum fluctuations present in the latter. We focus on high-order harmonic generation (HHG), where the approach is accurate in parameter regimes of current interest and it gives closed-form expressions for key observables. This formulation not only simplifies numerical calculations compared to the Schrödinger-picture approach but also provides a clear correspondence between nonclassical features of the emitted light and the underlying induced dynamics of the generating medium including quantum fluctuations. Moreover, the Heisenberg framework naturally yields scaling relations with the number of independent emitters, enabling us to assess whether nonclassical behavior should persist under typical experimental conditions involving large emitter ensembles. Interestingly, we find that the degree of squeezing increases with the number of emitters, whereas the photon statistics approaches a classical Poissonian distribution in the many-emitter limit. We also find that the beyond-semiclassical emitter dynamics significantly enhances the degree of squeezing of the emitted light. Our work advances the theoretical understanding of quantum-optical HHG and introduces an accessible and well-controlled framework to describe realistic experiments.

quant-pharXiv:2512.14174•Physical Review A•Dec 17, 2025•1 citation

Nonclassical light generation and control from laser-driven semiconductor intraband excitations

Ivan Gonoskov, René Sondenheimer, Christian Hünecke +3 more

The generation of higher-order harmonic radiation originating from the interaction of intense laser pulses with matter is typically described semiclassically: While the electronic structure and dynamics of matter is described quantum mechanically, the intense light field is described classically---and accordingly the generated harmonic radiation. However, pioneering experiments on quantum optical properties of high-order harmonic generation (HHG) in atomic gases from the group of P. Tzallas [I. Gonoskov et al., Sci. Rep. 6, 32821 (2016) and N. Tsatrafyllis et al., Nat. Commun. 8, 15170 (2017)], and theoretical investigations from the group of M. Lewenstein [M. Lewenstein et al., Nat. Phys. 17, 1104 (2021) and P. Stammer et al., Phys. Rev. Lett. 128, 123603 (2022)] have impressively demonstrated that light's quantum properties are observable in the strong-field realm. In this paper, we develop a quantum optical description of HHG in a bulk semiconductor originating from the nonlinear Bloch current. This mechanism of HHG, known as the intraband current, constitutes the major contribution to the emission spectrum for harmonic with energies of quanta roughly below the bandgap. Under certain approximations, employing a quantum description of both light and matter, we obtain analytical solutions, which allow us to analyze the classical and quantum optical characteristics of the fundamental mode of light and the harmonic modes. We find intricate but sufficiently mild modifications of the photon statistics of the fundamental mode and coherent displacements depending on the parameters of the driving laser field. Similar to high-harmonic generation in atoms, the fundamental and emitted harmonic field modes are entangled. Moreover, our analytical model predicts parameter ranges where these quantum optical properties will be most pronounced, allowing protocols for quantum information processing with high photon numbers over a large range of frequencies.

Physical review. B./Physical review. B•Mar 7, 2024•23 citations
We extend our recently developed semiclassical strong-field Herman-Kluk propagator (SFHK) method to calculate high-order harmonic generation (HHG) for atoms in intense lasers. We show that our method, based on a combination of the Herman-Kluk propagator and the strong-field approximation, can provide highly accurate results for both HHG yield and phase, nearly identical to those from the exact numerical solutions of the time-dependent Schrödinger equation. We provide detailed analyses of our method and its applications to the HHG process, particularly the recombination time. The main computational task in this method is to solve the classical Newton equations for the active electron in the combined atomic potential and laser-electron interaction. The motion of the centers of the electron wave packets, modeled by coherent states, is governed by independent classical trajectories so that the computation can therefore be parallelized very efficiently.
We extend our recently developed semiclassical strong-field Herman-Kluk (SFHK) propagator method to calculate high-order harmonic generation (HHG) in diatomic molecules driven by few-cycle intense laser fields. On the example of applications to H2 and N2, we show that our method, based on a combination of the Herman-Kluk propagator and the strong-field approximation, can provide very accurate results for both HHG yield and phase, nearly identical to those from the exact numerical solutions of the time-dependent Schrodinger equation. To compare with experimental measurements, averaging over molecular orientations must be performed. Here we demonstrate a distinct and powerful advantage of the SFHK, as its Monte Carlo sampling for the integration over the alignment distribution can be efficiently combined with the integration over the initial momentum distributions of electron wave-packet right after the tunnel exit. Therefore, the total number of trajectories used for the alignment-averaged HHG spectrum does not increase much compared to that for a single fixed alignment. Similar to atomic targets, the main computational task in the SFHK is to solve the classical Hamiltonian equations for the active electron in the combined electron-target ion potential and electron-laser interaction. The motion of the center of each electron wave packet in the continuum, represented by a coherent state, is governed by an independent classical trajectory so that the computation can be parallelized very efficiently.
The semiclassical Herman–Kluk initial value approximation of the quantum-mechanical propagator, which is heavily used in dynamical studies of atomic and molecular systems, takes into account nonlinearity of the underlying classical dynamics by using multiple classical initial value solutions. In applications to the propagation of Gaussian wavepackets non-Gaussian distortions in the course of time and therefore a realistic description of the quantum dynamics are possible. This is in contrast to the single-trajectory approximations investigated in a recent paper by Baranger et al , one of which was erroneously termed Herman–Kluk approximation.
We analytically describe the strong-field light-electron interaction using a quantized coherent laser state with arbitrary photon number. We obtain a light-electron wave function which is a closed-form solution of the time-dependent Schrödinger equation (TDSE). This wave function provides information about the quantum optical features of the interaction not accessible by semi-classical theories. With this approach we can reveal the quantum optical properties of high harmonic generation (HHG) process in gases by measuring the photon statistics of the transmitted infrared (IR) laser radiation. This work can lead to novel experiments in high-resolution spectroscopy in extreme-ultraviolet (XUV) and attosecond science without the need to measure the XUV light, while it can pave the way for the development of intense non-classical light sources.
Quantum-optical descriptions of strong-field processes have attracted significant attention in recent years. Typically, the theoretical modeling has been conducted in the Schrödinger picture, where results are only obtainable under certain approximations, while, in contrast, the Heisenberg picture has remained relatively unexplored. In this work, we develop an accurately controlled perturbative expansion of the time-evolution operator in the Heisenberg picture and derive beyond-semiclassical corrections to the emitter dynamics due to the coupling to the quantized electromagnetic field, capturing effects of the quantum fluctuations present in the latter. We focus on high-order harmonic generation (HHG), where the approach is accurate in parameter regimes of current interest and it gives closed-form expressions for key observables. This formulation not only simplifies numerical calculations compared to the Schrödinger-picture approach but also provides a clear correspondence between nonclassical features of the emitted light and the underlying induced dynamics of the generating medium including quantum fluctuations. Moreover, the Heisenberg framework naturally yields scaling relations with the number of independent emitters, enabling us to assess whether nonclassical behavior should persist under typical experimental conditions involving large emitter ensembles. Interestingly, we find that the degree of squeezing increases with the number of emitters, whereas the photon statistics approaches a classical Poissonian distribution in the many-emitter limit. We also find that the beyond-semiclassical emitter dynamics significantly enhances the degree of squeezing of the emitted light. Our work advances the theoretical understanding of quantum-optical HHG and introduces an accessible and well-controlled framework to describe realistic experiments.
The generation of higher-order harmonic radiation originating from the interaction of intense laser pulses with matter is typically described semiclassically: While the electronic structure and dynamics of matter is described quantum mechanically, the intense light field is described classically---and accordingly the generated harmonic radiation. However, pioneering experiments on quantum optical properties of high-order harmonic generation (HHG) in atomic gases from the group of P. Tzallas [I. Gonoskov et al., Sci. Rep. 6, 32821 (2016) and N. Tsatrafyllis et al., Nat. Commun. 8, 15170 (2017)], and theoretical investigations from the group of M. Lewenstein [M. Lewenstein et al., Nat. Phys. 17, 1104 (2021) and P. Stammer et al., Phys. Rev. Lett. 128, 123603 (2022)] have impressively demonstrated that light's quantum properties are observable in the strong-field realm. In this paper, we develop a quantum optical description of HHG in a bulk semiconductor originating from the nonlinear Bloch current. This mechanism of HHG, known as the intraband current, constitutes the major contribution to the emission spectrum for harmonic with energies of quanta roughly below the bandgap. Under certain approximations, employing a quantum description of both light and matter, we obtain analytical solutions, which allow us to analyze the classical and quantum optical characteristics of the fundamental mode of light and the harmonic modes. We find intricate but sufficiently mild modifications of the photon statistics of the fundamental mode and coherent displacements depending on the parameters of the driving laser field. Similar to high-harmonic generation in atoms, the fundamental and emitted harmonic field modes are entangled. Moreover, our analytical model predicts parameter ranges where these quantum optical properties will be most pronounced, allowing protocols for quantum information processing with high photon numbers over a large range of frequencies.