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:2608.19985

Resonance Raman spectroscopy from ab initio Hagedorn wavepacket dynamics

Davide Barbiero, Léa Zupan, Jiří J. L. Vaníček

physics.chem-phphysics.comp-phquant-phSubmitted Aug 20, 2026• Announced Aug 21, 2026
arXivHTMLPDF

Abstract

We present a practical, ab initio time-dependent method using Hagedorn wavepackets to simulate resonance Raman (RR) spectra of polyatomic molecules. Hagedorn functions---Gaussians multiplied by specific polynomials---are used to represent RR initial and final states because these functions are exact solutions to the time-dependent Schrödinger equation for at-most-quadratic potentials and can be propagated at zero cost beyond that of propagating the guiding Gaussian. Using efficient recursive formulae to compute overlaps between Hagedorn wavepackets, we can evaluate RR excitation profiles for arbitrary spectral signals, such as fundamental, overtone, combination, and hot bands. We then construct the Stokes and anti-Stokes RR spectra from these profiles. We first validate the method in a two-dimensional displaced, distorted, and Duschinsky-rotated harmonic model against numerically exact split-operator calculations. Then, we apply the method to compute RR spectra of anthracene by performing dynamics on a 66-dimensional harmonic potential energy surface constructed from density functional theory calculations.

Related papers

On Hagedorn wavepackets associated with different Gaussians

Jiří J. L. Vaníček, Zhan Tong Zhang

Abstract Hagedorn functions are carefully constructed generalizations of Hermite functions to the setting of many-dimensional squeezed and coupled harmonic systems. Wavepackets formed by superpositions of Hagedorn functions have been successfully used to solve the time-dependent Schrödinger equation exactly in harmonic systems and variationally in anharmonic systems. To evaluate typical observables, such as position or kinetic energy, it is sufficient to consider orthonormal Hagedorn functions with a single Gaussian center. Instead, we derive various relations between Hagedorn bases associated with different Gaussians, including their overlaps, which are necessary for evaluating quantities nonlocal in time, such as the time correlation functions needed for computing spectra. First, we use the Bogoliubov transformation to obtain the commutation relations between the ladder operators associated with different Gaussians. Then, instead of using numerical quadrature, we employ these commutation relations to derive exact recurrence relations for the overlap integrals between Hagedorn functions with different Gaussian centers. Finally, we present numerical experiments that demonstrate the accuracy and efficiency of our algebraic method as well as its suitability for treating problems in spectroscopy and chemical dynamics.

arXiv:2405.07880•Journal of Physics A Mathematical and Theoretical•Jan 30, 2025•6 citations

Nonzero-temperature vibronic spectra of polyatomic molecules from a zero-temperature classical trajectory

Davide Barbiero, Jiří J. L. Vaníček

By combining coherence thermofield dynamics with the single-Hessian approximation, we enable simulations of low- to medium-resolution vibronic spectra of weakly anharmonic systems at nonzero temperatures, at negligible additional cost relative to zero-temperature calculations. Single-Hessian coherence thermofield Gaussian wavepacket dynamics is exact in any harmonic potential, provided that the reference Hessian is that of the final surface. When applied to Morse systems of increasing anharmonicity and varying temperature, this method successfully captures excited-state anharmonicity and key temperature-dependent spectral features, including hot bands and broadening. By combining the method with on-the-fly ab initio dynamics, we demonstrate its utility by computing the absorption spectra of naphthalene, aminocoumarin C450, and phenyl radical, and the photoelectron spectrum of SeO$_{2}^{-}$ . Within the ab initio single-Hessian approximation, after the zero-temperature spectrum is obtained at the cost of classical molecular dynamics (on the order of hours), all nonzero-temperature spectra are computed in seconds.

physics.chem-phCross-listedarXiv:2608.20075•Aug 21, 2026

Nonequilibrium DMFT approach to time-resolved Raman spectroscopy

Philipp Werner, Martin Eckstein, Naoto Tsuji

Raman spectroscopy uses light scattering to extract information on low-energy excitations of solids. The Raman process is described by diagrams which are fourth order in the light-matter interaction and in particular the resonant contribution, which involves four different space-time arguments, is difficult to evaluate. If one instead simulates explicitly the incoming (classical) light pulse, the Raman signal is given by the outgoing photon flux and can be determined from a two-point correlation function. Such a formalism can be used to compute the time-resolved Raman spectrum of nonequilibrium systems, as well as nonlinear signals which are higher order in the incoming field, such as hyper-Raman scattering. Here we explain how to implement this time-dependent formalism within the dynamical mean field theory framework. The method is illustrated with applications to the Holstein-Hubbard model in the strong electron-phonon coupling regime. We demonstrate hyper-Raman scattering in measurements with strong probe fields and frequency mixing signals in the presence of a pump field and simulate the evolution of Stokes and anti-Stokes features after photoexcitations of metallic and Mott insulating systems.

arXiv:2309.02866•Physical review. B./Physical review. B•Dec 27, 2023•2 citations

Spawning semiclassical wavepackets

Vasile Catrinel Grădinaru, Oliver Rietmann

The semiclassical (or Hagedorn) wavepackets depending on a fixed set of parameters are an orthonormal L 2 -basis of generalized coherent states. They have been used to solve numerically the time-dependent Schrödinger equation in its semiclassical formulation, yet their localization property makes them inefficient in case of non-local phenomena such as quantum tunneling. In order to overcome this difficulty, we use simultaneously several bases with different parameters. We propose an algorithm to expand a given wavefunction in terms of multiple families of Hagedorn wavepackets; each family can then be accurately and efficiently propagated using modern semiclassical time-splittings.

Journal of Computational Physics•Apr 18, 2024•5 citations

A theoretical study of resonance Raman scattering from molecules. II. Temperature effect

Yuichi Fujimura, S. H. Lin

The temperature effect on the resonance Raman scattering from molecules is investigated. Considering a total system consisting of a molecular system, a heat bath, and the weak radiation field, a general expression for the relaxation processes is derived utilizing the Liouville operator technique for the time evolution of the density matrix of the total system. By using the perturbation method with respect to the radiation and matter up to second order, rate constants for zero photon (radiationless processes), one-photon and two-photon processes are derived after taking the long time limit in the Markoff approximation. Assuming a weak interaction between the molecular system and heat bath, a temperature-dependent resonance Raman scattering cross section is formulated in the adiabatic basis set for the molecular systems. In the displaced harmonic oscillator model, we obtain an analytical expression for the temperature-dependent Raman scattering cross section for the multimode case, which is suitable for both weak coupling and strong coupling cases. Analytical expressions for the Stokes and anti-Stokes bands of the jth order vibrational transition in the one oscillator model can easily be obtained from the general expression for the resonance Raman scattering cross section derived. An approximate Stokes to anti-Stokes ratio is given up to the second order of the dimensionless displacement Δ. The resonance Raman scattering cross section and the excitation energy profiles are numerically calculated as a function of physical constant such as temperature, molecular vibrational frequency and Δ. It is shown that the excitation profiles of the anti-Stokes band of the jth order vibrational transition have the similar envelope of the corresponding Stokes band shifted by jh/ω. The inverse Raman effect which is closely related to the ordinary Raman effect is also discussed.

The Journal of Chemical Physics•Nov 1, 1979•38 citations

Non-Hermitian Schrödinger dynamics with Hagedorn’s wave packets

Stephanie Troppmann

Hagedorn wave packets form an orthonormal basis of the L2-functions and are known solutions of the semiclassical Schrödinger equation with quadratic Hamiltonians. In this thesis we analyse their structure, in particular the polynomial part, their connection to the Hermite functions and their representation in phase space. Based on these findings we utilise them to investigate the time evolution generated by the Schrödinger equation with a non-Hermitian, quadratic Hamiltonian. Operators of this type appear for example in the context of diffusion models and absorbing potentials and are thus of interest in physics and chemistry. We provide explicit formulas for the propagated wave packets and thereby show that the non-unitary evolution activates lower excited states. At last we apply our results to the Lindblad master equation with quadratic internal Hamiltonian and linear Lindblad terms.

•Jan 1, 2017•1 citation
We present a practical, ab initio time-dependent method using Hagedorn wavepackets to simulate resonance Raman (RR) spectra of polyatomic molecules. Hagedorn functions---Gaussians multiplied by specific polynomials---are used to represent RR initial and final states because these functions are exact solutions to the time-dependent Schrödinger equation for at-most-quadratic potentials and can be propagated at zero cost beyond that of propagating the guiding Gaussian. Using efficient recursive formulae to compute overlaps between Hagedorn wavepackets, we can evaluate RR excitation profiles for arbitrary spectral signals, such as fundamental, overtone, combination, and hot bands. We then construct the Stokes and anti-Stokes RR spectra from these profiles. We first validate the method in a two-dimensional displaced, distorted, and Duschinsky-rotated harmonic model against numerically exact split-operator calculations. Then, we apply the method to compute RR spectra of anthracene by performing dynamics on a 66-dimensional harmonic potential energy surface constructed from density functional theory calculations.
Abstract Hagedorn functions are carefully constructed generalizations of Hermite functions to the setting of many-dimensional squeezed and coupled harmonic systems. Wavepackets formed by superpositions of Hagedorn functions have been successfully used to solve the time-dependent Schrödinger equation exactly in harmonic systems and variationally in anharmonic systems. To evaluate typical observables, such as position or kinetic energy, it is sufficient to consider orthonormal Hagedorn functions with a single Gaussian center. Instead, we derive various relations between Hagedorn bases associated with different Gaussians, including their overlaps, which are necessary for evaluating quantities nonlocal in time, such as the time correlation functions needed for computing spectra. First, we use the Bogoliubov transformation to obtain the commutation relations between the ladder operators associated with different Gaussians. Then, instead of using numerical quadrature, we employ these commutation relations to derive exact recurrence relations for the overlap integrals between Hagedorn functions with different Gaussian centers. Finally, we present numerical experiments that demonstrate the accuracy and efficiency of our algebraic method as well as its suitability for treating problems in spectroscopy and chemical dynamics.
By combining coherence thermofield dynamics with the single-Hessian approximation, we enable simulations of low- to medium-resolution vibronic spectra of weakly anharmonic systems at nonzero temperatures, at negligible additional cost relative to zero-temperature calculations. Single-Hessian coherence thermofield Gaussian wavepacket dynamics is exact in any harmonic potential, provided that the reference Hessian is that of the final surface. When applied to Morse systems of increasing anharmonicity and varying temperature, this method successfully captures excited-state anharmonicity and key temperature-dependent spectral features, including hot bands and broadening. By combining the method with on-the-fly ab initio dynamics, we demonstrate its utility by computing the absorption spectra of naphthalene, aminocoumarin C450, and phenyl radical, and the photoelectron spectrum of SeO2−_{2}^{-}2−​ . Within the ab initio single-Hessian approximation, after the zero-temperature spectrum is obtained at the cost of classical molecular dynamics (on the order of hours), all nonzero-temperature spectra are computed in seconds.
Raman spectroscopy uses light scattering to extract information on low-energy excitations of solids. The Raman process is described by diagrams which are fourth order in the light-matter interaction and in particular the resonant contribution, which involves four different space-time arguments, is difficult to evaluate. If one instead simulates explicitly the incoming (classical) light pulse, the Raman signal is given by the outgoing photon flux and can be determined from a two-point correlation function. Such a formalism can be used to compute the time-resolved Raman spectrum of nonequilibrium systems, as well as nonlinear signals which are higher order in the incoming field, such as hyper-Raman scattering. Here we explain how to implement this time-dependent formalism within the dynamical mean field theory framework. The method is illustrated with applications to the Holstein-Hubbard model in the strong electron-phonon coupling regime. We demonstrate hyper-Raman scattering in measurements with strong probe fields and frequency mixing signals in the presence of a pump field and simulate the evolution of Stokes and anti-Stokes features after photoexcitations of metallic and Mott insulating systems.
The semiclassical (or Hagedorn) wavepackets depending on a fixed set of parameters are an orthonormal L 2 -basis of generalized coherent states. They have been used to solve numerically the time-dependent Schrödinger equation in its semiclassical formulation, yet their localization property makes them inefficient in case of non-local phenomena such as quantum tunneling. In order to overcome this difficulty, we use simultaneously several bases with different parameters. We propose an algorithm to expand a given wavefunction in terms of multiple families of Hagedorn wavepackets; each family can then be accurately and efficiently propagated using modern semiclassical time-splittings.
The temperature effect on the resonance Raman scattering from molecules is investigated. Considering a total system consisting of a molecular system, a heat bath, and the weak radiation field, a general expression for the relaxation processes is derived utilizing the Liouville operator technique for the time evolution of the density matrix of the total system. By using the perturbation method with respect to the radiation and matter up to second order, rate constants for zero photon (radiationless processes), one-photon and two-photon processes are derived after taking the long time limit in the Markoff approximation. Assuming a weak interaction between the molecular system and heat bath, a temperature-dependent resonance Raman scattering cross section is formulated in the adiabatic basis set for the molecular systems. In the displaced harmonic oscillator model, we obtain an analytical expression for the temperature-dependent Raman scattering cross section for the multimode case, which is suitable for both weak coupling and strong coupling cases. Analytical expressions for the Stokes and anti-Stokes bands of the jth order vibrational transition in the one oscillator model can easily be obtained from the general expression for the resonance Raman scattering cross section derived. An approximate Stokes to anti-Stokes ratio is given up to the second order of the dimensionless displacement Δ. The resonance Raman scattering cross section and the excitation energy profiles are numerically calculated as a function of physical constant such as temperature, molecular vibrational frequency and Δ. It is shown that the excitation profiles of the anti-Stokes band of the jth order vibrational transition have the similar envelope of the corresponding Stokes band shifted by jh/ω. The inverse Raman effect which is closely related to the ordinary Raman effect is also discussed.
Hagedorn wave packets form an orthonormal basis of the L2-functions and are known solutions of the semiclassical Schrödinger equation with quadratic Hamiltonians. In this thesis we analyse their structure, in particular the polynomial part, their connection to the Hermite functions and their representation in phase space. Based on these findings we utilise them to investigate the time evolution generated by the Schrödinger equation with a non-Hermitian, quadratic Hamiltonian. Operators of this type appear for example in the context of diffusion models and absorbing potentials and are thus of interest in physics and chemistry. We provide explicit formulas for the propagated wave packets and thereby show that the non-unitary evolution activates lower excited states. At last we apply our results to the Lindblad master equation with quadratic internal Hamiltonian and linear Lindblad terms.