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

Efficient Algorithms for Semiclassical Quantum Dynamics

David Sattlegger

Published Jan 1, 2016• 0 citations
Publisher

Abstract

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.

DOI
10.14288/1.0307179

Related papers

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

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

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.

physics.atom-phCross-listedarXiv:2601.17164•Jan 27, 2026

Semiclassical dynamics with quantum trajectories: Formulation and comparison with the semiclassical initial value representation propagator

Sophya V. Garashchuk, Vitaly A. Rassolov

We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.

The Journal of Chemical Physics•Jan 30, 2003•54 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

Semiclassical corrections to the Herman-Kluk propagator

Gili Hochman, Kenneth G. Kay

The approximate expression for the propagator of Herman and Kluk is systematically improved by application of semiclassical corrections for a double-well system. The results confirm the semiclassical nature of the Herman-Kluk approximation. The influence of caustics on the effectiveness of the corrections is discussed.

Physical Review A•Jun 27, 2006•14 citations

Initial-value semiclassical propagators for the Wigner phase space representation: Formulation based on the interpretation of the Moyal equation as a Schrödinger equation

Shin-ichi Koda

We formulate various semiclassical propagators for the Wigner phase space representation from a unified point of view. As is shown in several studies, the Moyal equation, which is an equation of motion for the Wigner distribution function, can be regarded as the Schrödinger equation of an extended Hamiltonian system where its "position" and "momentum" correspond to the middle point of two points of the original phase space and the difference between them, respectively. Then we show that various phase-space semiclassical propagators can be formulated just by applying existing semiclassical propagators to the extended system. As a result, a phase space version of the Van Vleck propagator, the initial-value Van Vleck propagator, the Herman-Kluk propagator, and the thawed Gaussian approximation are obtained. In addition, we numerically compare the initial-value phase-space Van Vleck propagator, the phase-space Herman-Kluk propagator, and the classical mechanical propagation as approximation methods for the time propagation of the Wigner distribution function in terms of both accuracy and convergence speed. As a result, we find that the convergence speed of the Van Vleck propagator is far slower than others as is the case of the Hilbert space, and the Herman-Kluk propagator keeps its accuracy for a long period compared with the classical mechanical propagation while the convergence speed of the latter is faster than the former.

The Journal of Chemical Physics•Dec 24, 2015•19 citations

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
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.
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 present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.
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 formulate various semiclassical propagators for the Wigner phase space representation from a unified point of view. As is shown in several studies, the Moyal equation, which is an equation of motion for the Wigner distribution function, can be regarded as the Schrödinger equation of an extended Hamiltonian system where its "position" and "momentum" correspond to the middle point of two points of the original phase space and the difference between them, respectively. Then we show that various phase-space semiclassical propagators can be formulated just by applying existing semiclassical propagators to the extended system. As a result, a phase space version of the Van Vleck propagator, the initial-value Van Vleck propagator, the Herman-Kluk propagator, and the thawed Gaussian approximation are obtained. In addition, we numerically compare the initial-value phase-space Van Vleck propagator, the phase-space Herman-Kluk propagator, and the classical mechanical propagation as approximation methods for the time propagation of the Wigner distribution function in terms of both accuracy and convergence speed. As a result, we find that the convergence speed of the Van Vleck propagator is far slower than others as is the case of the Hilbert space, and the Herman-Kluk propagator keeps its accuracy for a long period compared with the classical mechanical propagation while the convergence speed of the latter is faster than the former.
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.