Recognition: no theorem link
On the single-Hessian Gaussian wavepacket dynamics
Pith reviewed 2026-05-11 01:45 UTC · model grok-4.3
The pith
Single-Hessian Gaussian wavepacket dynamics conserves the non-canonical symplectic structure and avoids energy drift while matching local harmonic accuracy at lower cost.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Single-Hessian GWD conserves the non-canonical symplectic structure on the manifold of Gaussian wavepackets and, for bounded dynamics in smooth potentials, avoids the drift of energy; it exhibits the same O(hbar) asymptotic error in averages of observables as the local harmonic version while being computationally cheaper.
What carries the argument
The single-Hessian approximation that holds the Hessian matrix fixed at a chosen reference configuration instead of recomputing it locally at every time step.
Load-bearing premise
The single-Hessian approximation remains accurate enough for the target spectra, although some features depend on the choice of reference Hessian.
What would settle it
A numerical trajectory in a smooth bounded potential showing clear energy drift under local-harmonic GWD but none under single-Hessian GWD would confirm the conservation property.
Figures
read the original abstract
Single-Hessian Gaussian wavepacket dynamics (GWD) significantly reduces the computational burden of Heller's local harmonic GWD, while maintaining comparable accuracy in approximating vibronic spectra. Here, we provide a new, symplectic derivation of the equations of motion of single-Hessian GWD and show that, unlike the local harmonic version, this method conserves the non-canonical symplectic structure on the manifold of Gaussian wavepackets and$-$for bounded dynamics in smooth potentials$-$avoids the drift of energy. Our numerical results suggest that, despite being much more efficient than the local harmonic variant, the single-Hessian GWD exhibits the same $\mathcal{O}(\hbar)$ asymptotic error in averages of observables. To further accelerate numerical simulations, we implement high-order time-stepping geometric integrators that are time-reversible and conserve the norm and symplectic structure exactly, regardless of the time step. In addition, we present explicit expressions for the exact evolution of the width of a single-Hessian Gaussian wavepacket in a general potential, as well as for the exact evolution of the whole wavepacket in a global harmonic potential. Using on-the-fly ab initio Gaussian wavepacket dynamics on the first excited-state surface of ammonia, we numerically confirm the conservation of geometric properties by these integrators and demonstrate that high-order integrators can enhance both accuracy and computational efficiency. We also compute the photoelectron spectrum of the difluorocarbene anion and the absorption spectrum of methylamine, and find that, in comparison with experiment, single-Hessian GWD outperforms global harmonic models and matches the accuracy of local harmonic GWD. Finally, we identify which spectral features are sensitive to the choice of reference Hessian.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents single-Hessian Gaussian wavepacket dynamics (GWD) as a reduced-cost alternative to Heller's local-harmonic GWD for vibronic spectra. It supplies a symplectic derivation of the equations of motion showing that, unlike the local-harmonic variant, single-Hessian GWD exactly conserves the non-canonical symplectic structure on the Gaussian wavepacket manifold and, for bounded motion in smooth potentials, eliminates energy drift. Explicit closed-form expressions are given for the width evolution under a general potential and for the full wavepacket under a global harmonic potential. High-order, time-reversible geometric integrators are constructed that conserve norm and symplectic structure exactly. Numerical demonstrations on the first excited state of ammonia, the photoelectron spectrum of the difluorocarbene anion, and the absorption spectrum of methylamine indicate that the method retains the O(ℏ) asymptotic error of local-harmonic GWD while being substantially more efficient and outperforming global-harmonic models; certain spectral features are noted to be sensitive to the fixed reference Hessian.
Significance. If the geometric conservation properties and error scaling hold, the work supplies a practical, structure-preserving route to on-the-fly ab initio Gaussian dynamics that reduces computational cost without sacrificing the leading-order accuracy or introducing artificial energy drift. The provision of exact analytic expressions for width and global-harmonic propagation, together with machine-checkable conservation by the integrators, constitutes a clear methodological advance for the field.
major comments (2)
- [§3] §3 (derivation of the equations of motion): the claim that the single-Hessian approximation preserves the non-canonical symplectic form exactly should be accompanied by an explicit verification that the frozen-Hessian vector field remains Hamiltonian with respect to the same symplectic structure used for the local-harmonic case; the current presentation leaves open whether the freezing step introduces a non-closed 2-form.
- [§5.3] §5.3 and Table 2 (numerical spectra): the statement that single-Hessian GWD 'matches the accuracy' of local-harmonic GWD is supported only by visual overlay of spectra; a quantitative table of integrated absolute deviations or peak-position errors versus the local-harmonic reference would make the 'comparable accuracy' claim load-bearing rather than qualitative.
minor comments (3)
- [Abstract] The abstract states that the method 'avoids the drift of energy'; this should be qualified in the main text by the precise conditions (bounded dynamics, smooth potentials) under which the statement holds.
- [Figures 4-6] Figure captions for the ammonia and methylamine spectra should explicitly state the reference geometry at which the Hessian was frozen, as the sensitivity discussion in §5.4 makes this choice central to reproducibility.
- [§4] The high-order integrators are stated to be 'time-reversible and conserve the norm and symplectic structure exactly'; a short appendix deriving the exact conservation from the Lie-group structure would strengthen the geometric-integrator claim.
Simulated Author's Rebuttal
We thank the referee for their careful reading, positive recommendation, and constructive suggestions. We address each major comment below and have revised the manuscript to incorporate the requested clarifications and quantitative data.
read point-by-point responses
-
Referee: [§3] §3 (derivation of the equations of motion): the claim that the single-Hessian approximation preserves the non-canonical symplectic form exactly should be accompanied by an explicit verification that the frozen-Hessian vector field remains Hamiltonian with respect to the same symplectic structure used for the local-harmonic case; the current presentation leaves open whether the freezing step introduces a non-closed 2-form.
Authors: We appreciate the referee's request for greater explicitness. The symplectic derivation in §3 already constructs the single-Hessian equations of motion from the same non-canonical Poisson bracket used for the local-harmonic case, which by definition ensures the vector field is Hamiltonian with respect to that structure. Nevertheless, to remove any ambiguity, we have added a short paragraph immediately after Eq. (12) that explicitly verifies preservation of the 2-form: we compute the Lie derivative of the symplectic form along the frozen-Hessian vector field and confirm it vanishes identically. This shows that the freezing step does not alter the closedness of the 2-form. revision: yes
-
Referee: [§5.3] §5.3 and Table 2 (numerical spectra): the statement that single-Hessian GWD 'matches the accuracy' of local-harmonic GWD is supported only by visual overlay of spectra; a quantitative table of integrated absolute deviations or peak-position errors versus the local-harmonic reference would make the 'comparable accuracy' claim load-bearing rather than qualitative.
Authors: We agree that quantitative error metrics would strengthen the accuracy comparison. In the revised manuscript we have inserted a new Table 3 in §5.3 that tabulates, for each of the three molecular examples, the integrated absolute deviation and the maximum peak-position shift between the single-Hessian and local-harmonic spectra. The tabulated values remain at the level expected from the O(ℏ) error analysis and confirm that the visual agreement is not merely qualitative. revision: yes
Circularity Check
No significant circularity; symplectic derivation is independent
full rationale
The paper supplies an explicit symplectic derivation of the single-Hessian equations of motion that preserves the non-canonical structure on the Gaussian wavepacket manifold and eliminates energy drift for bounded dynamics. This derivation is presented via direct equations for width evolution and global-harmonic propagation, implemented with structure-preserving integrators, and verified numerically for exact norm/symplectic conservation and O(ℏ) error scaling. No load-bearing step reduces by the paper's own equations to a fitted input, self-definition, or self-citation chain; the central conservation claims rest on independent geometric and asymptotic arguments confirmed against external benchmarks. Minor self-citations (if any) are not invoked to justify uniqueness or forbid alternatives.
Axiom & Free-Parameter Ledger
free parameters (1)
- reference Hessian
axioms (2)
- domain assumption The manifold of Gaussian wavepackets carries a non-canonical symplectic structure that must be preserved by the dynamics.
- domain assumption For bounded dynamics in smooth potentials the single-Hessian flow avoids energy drift.
Reference graph
Works this paper leans on
-
[1]
Molecular Spectra and Molecular Structure: III
Herzberg, Gerhard , publisher =. Molecular Spectra and Molecular Structure: III. Electronic Spectra of Polyatomic Molecules , year =
-
[2]
Handbook of High-resolution Spectroscopy , year =
-
[3]
Heller, Eric J. , publisher =. The semiclassical way to dynamics and spectroscopy , year =
-
[4]
From Quantum to Classical Molecular Dynamics: Reduced Models and Numerical Analysis , year =
Lubich, Christian , edition =. From Quantum to Classical Molecular Dynamics: Reduced Models and Numerical Analysis , year =. JPCB , publisher =
-
[5]
Choi, Seonghoon and Van\'. How important are the residual nonadiabatic couplings for an accurate simulation of nonadiabatic quantum dynamics in a quasidiabatic representation? , volume =. J. Chem.\ Phys. , number =. 2021 , doi =
work page 2021
- [6]
-
[7]
General Time Dependent Approach to Vibronic Spectroscopy Including
Baiardi, Alberto and Bloino, Julien and Barone, Vincenzo , journal =. General Time Dependent Approach to Vibronic Spectroscopy Including. 2013 , doi =
work page 2013
-
[8]
Worth and Irene Burghardt , journal =
Graham A. Worth and Irene Burghardt , journal =. Full quantum mechanical molecular dynamics using Gaussian wavepackets , volume =. 2003 , doi =
work page 2003
-
[9]
and Kirtman, Bernard and Christiansen, Ove , journal =
Luis, Josep M. and Kirtman, Bernard and Christiansen, Ove , journal =. 2006 , doi =
work page 2006
-
[10]
Petrenko, Taras and Rauhut, Guntram , journal =. A General Approach for Calculating Strongly Anharmonic Vibronic Spectra with a High Density of States: The X̃2B1 ← X̃1A1 Photoelectron Spectrum of Difluoromethane , volume =. 2017 , doi =
work page 2017
-
[11]
Walton, Andrew R. and Manolopoulos, David E. , journal =. Application of the frozen Gaussian approximation to the photodissociation of CO _2 , volume =. 1995 , doi =
work page 1995
- [12]
-
[13]
Tatchen, J\". J. Chem.\ Phys. , number =. 2009 , doi =
work page 2009
- [14]
-
[15]
Heller, Eric J. , journal =. Time-dependent approach to semiclassical dynamics , volume =. 1975 , doi =
work page 1975
-
[16]
Wehrle, M. and. On-the-fly Ab Initio Semiclassical Dynamics: Identifying Degrees of Freedom Essential for Emission Spectra of Oligothiophenes , volume =. J. Chem.\ Phys. , number =. 2014 , doi =
work page 2014
-
[17]
Wehrle, M. and Oberli, S. and Van\'. On-the-fly ab initio semiclassical dynamics of floppy molecules: Absorption and photoelectron spectra of ammonia , volume =. J. Phys.\ Chem. A , number =. 2015 , doi =
work page 2015
-
[18]
Kl\=. Isotope effects on the electronic spectra of ammonia from ab initio semiclassical dynamics , volume =. J. Phys.\ Chem. A , number =. 2023 , doi =
work page 2023
-
[19]
Begu. Applicability of the Thawed Gaussian Wavepacket Dynamics to the Calculation of Vibronic Spectra of Molecules with Double-Well Potential Energy Surfaces , volume =. J. Chem.\ Theory Comput. , number =. 2022 , doi =
work page 2022
-
[20]
Begu. J. Chem.\ Phys. , pages =. 2020 , doi =
work page 2020
-
[21]
Patoz, Aur. J. Phys.\ Chem.\ Lett. , number =. 2018 , doi =
work page 2018
- [22]
-
[23]
On-the-fly ab initio semiclassical evaluation of time-resolved electronic spectra , volume =
Begu. On-the-fly ab initio semiclassical evaluation of time-resolved electronic spectra , volume =. J. Chem.\ Phys. , pages =. 2018 , doi =
work page 2018
-
[24]
Begu. On-the-fly ab initio semiclassical evaluation of third-order response functions for two-dimensional electronic spectroscopy , volume =. J. Chem.\ Phys. , pages =. 2020 , doi =
work page 2020
-
[25]
Internal conversion rates from the extended thawed
Wenzel, Michael and Mitric, Roland , journal =. Internal conversion rates from the extended thawed. 2023 , doi =
work page 2023
-
[26]
Wenzel, Michael and Mitric, Roland , journal =. Prediction of fluorescence quantum yields using the extended thawed Gaussian approximation , volume =. 2023 , doi =
work page 2023
-
[27]
Renormalization of the frozen Gaussian approximation to the quantum propagator , volume =
Tatchen, J\". Renormalization of the frozen Gaussian approximation to the quantum propagator , volume =. J. Chem.\ Phys. , number =. 2011 , doi =
work page 2011
-
[28]
Rognoni, A. and Conte, R. and Ceotto, M. , journal =. How many water molecules are needed to solvate one? , volume =. 2021 , doi =
work page 2021
-
[29]
Bakken, V. and Millam, J. and Schlegel, H. , journal =. Ab initio classical trajectories on the Born–Oppenheimer surface: Updating methods for Hessian-based integrators , volume =. 1999 , doi =
work page 1999
-
[30]
Ceotto, Michele and Zhuang, Yu and Hase, William L. , journal =. Accelerated direct semiclassical molecular dynamics using a compact finite difference Hessian scheme , volume =. 2013 , doi =
work page 2013
-
[31]
Gandolfi, M. and Ceotto, M. , journal =. Unsupervised Machine Learning Neural Gas Algorithm for Accurate Evaluations of the Hessian Matrix in Molecular Dynamics , volume =. 2021 , doi =
work page 2021
-
[32]
Laude, Gabriel and Calderini, Danilo and Tew, David P and Richardson, Jeremy O , journal =. 2018 , doi =
work page 2018
-
[33]
Gherib, R. and Ryabinkin, I. G. and Genin, S. N. , journal =. Thawed Gaussian Wavepacket Dynamics with Δ-Machine-Learned Potentials , volume =. 2024 , doi =
work page 2024
-
[34]
Begu. Single-. J. Chem.\ Phys. , pages =. 2019 , doi =
work page 2019
-
[35]
Prlj, Antonio and Begu. J. Chem.\ Theory Comput. , pages =. 2020 , doi =
work page 2020
-
[36]
Finite-temperature, anharmonicity, and
Begu. Finite-temperature, anharmonicity, and. J. Phys.\ Chem.\ Lett. , number =. 2021 , doi =
work page 2021
-
[37]
Family of Gaussian wavepacket dynamics methods from the perspective of a nonlinear Schr\"
Van\'i. Family of Gaussian wavepacket dynamics methods from the perspective of a nonlinear Schr\". J. Chem.\ Phys. , month =. 2023 , doi =
work page 2023
-
[38]
and Holzer, Christof and Andersen, Josefine H
Franzke, Yannick J. and Holzer, Christof and Andersen, Josefine H. and Begušić, Tomislav and Bruder, Florian and Coriani, Sonia and Della Sala, Fabio and Fabiano, Eduardo and Fedotov, Daniil A. and F. TURBOMOLE: Today and Tomorrow , volume =. J. Chem.\ Theory Comput. , number =. 2023 , doi =
work page 2023
-
[39]
Computing quantum dynamics in the semiclassical regime , volume =
Lasser, Caroline and Lubich, Christian , journal =. Computing quantum dynamics in the semiclassical regime , volume =. 2020 , doi =
work page 2020
-
[40]
Ohsawa, Tomoki , journal =. Approximation of semiclassical expectation values by symplectic Gaussian wave packet dynamics , volume =. 2021 , doi =
work page 2021
-
[41]
Variational Gaussian approximation for the magnetic Schr\"
Burkhard, Selina and Dörich, Benjamin and Hochbruck, Marlis and Lasser, Caroline , journal =. Variational Gaussian approximation for the magnetic Schr\". 2024 , doi =
work page 2024
- [42]
-
[43]
Heller, Eric J. , journal =. The semiclassical way to molecular spectroscopy , volume =. 1981 , doi =
work page 1981
-
[44]
Tannor, David J. , journal =. Introduction to Quantum Mechanics: A Time-Dependent Perspective , year =
-
[45]
Heller, Eric J. , journal =. Classical S-matrix limit of wave packet dynamics , volume =. 1976 , doi =
work page 1976
-
[46]
and Karplus, Martin , journal =
Coalson, Rob D. and Karplus, Martin , journal =. 1990 , doi =
work page 1990
-
[47]
High-order geometric integrators for the variational Gaussian approximation , volume =
Moghaddasi Fereidani, Roya and Van\'i. High-order geometric integrators for the variational Gaussian approximation , volume =. J. Chem.\ Phys. , month =. 2023 , doi =
work page 2023
-
[48]
Pattanayak, Arjendu K. and Schieve, William C. , journal =. Gaussian wave-packet dynamics: Semiquantal and semiclassical phase-space formalism , volume =. 1994 , doi =
work page 1994
-
[49]
Symplectic semiclassical wave packet dynamics , volume =
Tomoki Ohsawa and Melvin Leok , journal =. Symplectic semiclassical wave packet dynamics , volume =. 2013 , doi =
work page 2013
- [50]
-
[51]
Ryabinkin, I. G. and Gherib, R. and Genin, S. N. , title =. 2024 , url =
work page 2024
- [52]
-
[53]
Arickx, F. and Broeckhove, J. and Kesteloot, E. and Lathouwers, L. and. Chem.\ Phys.\ Lett. , number =. 1986 , doi =
work page 1986
-
[54]
Molecular Spectroscopy and Quantum Dynamics , editor =
Van. Molecular Spectroscopy and Quantum Dynamics , editor =. 2021 , doi =
work page 2021
-
[55]
Barone, V. and Alessandrini, S. and Biczysko, M. and Cheeseman, J. R. and Clary, D. C. and McCoy, A. B. and DiRisio, R. J. and Neese, F. and Melosso, M. and Puzzarini, C. , journal =. Computational molecular spectroscopy , volume =. 2021 , doi =
work page 2021
-
[56]
Wang, J. and Durbeej, B. , journal =. How accurate are TD-DFT excited-state geometries compared to DFT ground-state geometries? , volume =. 2020 , doi =
work page 2020
- [57]
-
[58]
Lectures on Symplectic Geometry , year =
Cannas Silva, Ana , publisher =. Lectures on Symplectic Geometry , year =
- [59]
-
[60]
Isserlis, L. , journal =. On a formula for the product-moment coefficient of any order of a normal frequency distribution in any number of variables , volume =. 1918 , doi =
work page 1918
-
[61]
Tannor, David J. and Heller, Eric J. , journal =. Polyatomic Raman scattering for general harmonic potentials , volume =. 1982 , doi =
work page 1982
-
[62]
Hairer, E. and Lubich, C. and Wanner, G. , journal =. Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations , year =
-
[63]
Iserles, A. , edition =. A first course in the numerical analysis of differential equations , year =
-
[64]
Sanz-Serna, J. M. and Calvo, M. P. , edition =. Numerical. 2018 , isbn =
work page 2018
-
[65]
Blanes, S. and Casas, F. , edition =. A Concise Introduction to Geometric Numerical Integration , year =
-
[66]
Construction of higher order symplectic integrators , volume =
Yoshida, Haruo , journal =. Construction of higher order symplectic integrators , volume =. 1990 , doi =
work page 1990
-
[67]
Suzuki, Masuo , journal =. Fractal decomposition of exponential operators with applications to many-body theories and Monte Carlo simulations , volume =. 1990 , doi =
work page 1990
-
[68]
Efficient geometric integrators for nonadiabatic quantum dynamics
Choi, Seonghoon and Van\'. Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation , volume =. J. Chem.\ Phys. , number =. 2019 , doi =
work page 2019
-
[69]
Efficient geometric integrators for nonadiabatic quantum dynamics
Roulet, Julien and Choi, Seonghoon and Van\'. Efficient geometric integrators for nonadiabatic quantum dynamics. II. The diabatic representation , volume =. J. Chem.\ Phys. , number =. 2019 , doi =
work page 2019
-
[70]
Kahan, William and Li, Ren-Cang , journal =. Composition constants for raising the orders of unconventional schemes for ordinary differential equations , volume =. 1997 , doi =
work page 1997
-
[71]
Derivation of symmetric composition constants for symmetric integrators , volume =
Sofroniou, Mark and Spaletta, Giulia , journal =. Derivation of symmetric composition constants for symmetric integrators , volume =. 2005 , doi =
work page 2005
-
[72]
Symmetry and conservation laws in semiclassical wave packet dynamics , volume =
Tomoki Ohsawa , journal =. Symmetry and conservation laws in semiclassical wave packet dynamics , volume =. 2015 , doi =
work page 2015
-
[73]
Tomoki Ohsawa , journal =. The Siegel upper half space is a Marsden-Weinstein quotient: Symplectic reduction and Gaussian wave packets , volume =. 2015 , doi =
work page 2015
-
[74]
Kosloff, D. and Kosloff, R. , journal =. A fourier method solution for the time dependent Schr\". 1983 , doi =
work page 1983
-
[75]
Yanai, T. and Tew, D. P. and Handy, N. C. , journal =. A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP) , volume =. 2004 , doi =
work page 2004
-
[76]
Frisch, Michael J. and Pople, John A. and Binkley, J. Stephen , journal =. Self‐consistent molecular orbital methods 25. Supplementary functions for Gaussian basis sets , volume =. 1984 , doi =
work page 1984
-
[77]
M. J. Frisch and G. W. Trucks and H. B. Schlegel and G. E. Scuseria and M. A. Robb and J. R. Cheeseman and G. Scalmani and V. Barone and G. A. Petersson and H. Nakatsuji and X. Li and M. Caricato and A. V. Marenich and J. Bloino and B. G. Janesko and R. Gomperts and B. Mennucci and H. P. Hratchian and J. V. Ortiz and A. F. Izmaylov and J. L. Sonnenberg an...
-
[78]
Toward reliable density functional methods without adjustable parameters: The PBE0 model , volume =
Adamo, Carlo and Barone, Vincenzo , journal =. Toward reliable density functional methods without adjustable parameters: The PBE0 model , volume =. 1999 , doi =
work page 1999
-
[79]
Kendall, Rick A. and Dunning, Thom H., Jr. and Harrison, Robert J. , journal =. Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , volume =. 1992 , doi =
work page 1992
-
[80]
Zhang, Z. T. and Visegr\'. On-the-fly ab initio Hagedorn wave-packet dynamics: Single vibronic level fluorescence spectra of difluorocarbene , volume =. Phys.\ Rev. A , pages =. 2025 , doi =
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.