Gaussian wavepacket dynamics has proven to be a useful semiclassical approximation for quantum simulations of high-dimensional systems with low anharmonicity. Compared to Heller's original local harmonic method, the variational Gaussian wavepacket dynamics is more accurate, but much more difficult to apply in practice because it requires evaluating the expectation values of the potential energy, gradient, and Hessian. If the variational approach is applied to the local cubic approximation of the potential, these expectation values can be evaluated analytically, but they still require the costly third derivative of the potential. To reduce the cost of the resulting local cubic variational Gaussian wavepacket dynamics, we describe efficient high-order geometric integrators, which are symplectic, time-reversible, and norm-conserving. For small time steps, they also conserve the effective energy. We demonstrate the efficiency and geometric properties of these integrators numerically on a multidimensional, nonseparable coupled Morse potential.
Download full-text PDF |
Source |
|---|---|
| http://dx.doi.org/10.1063/5.0180070 | DOI Listing |
Publication Analysis
Top Keywords
Similar Publications
Thawed Gaussian Wavepacket Dynamics with Δ-Machine-Learned Potentials.
J Phys Chem A
October 2024
OTI Lumionics Inc., 3415 American Drive Unit 1, Mississauga, Ontario L4V 1T4, Canada.
A method for performing variable-width (thawed) Gaussian wavepacket (GWP) variational dynamics on machine-learned potentials is presented. Instead of fitting the potential energy surface (PES), the anharmonic correction to the global harmonic approximation (GHA) is fitted using kernel ridge regression─this is a Δ-machine learning approach. The training set consists of energy differences between ab initio electronic energies and values given by the GHA.
View Article and Find Full Text PDFSingle vibronic level fluorescence spectra from Hagedorn wavepacket dynamics.
J Chem Phys
September 2024
Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland.
In single vibronic level (SVL) fluorescence experiments, the electronically excited initial state is also excited in one or several vibrational modes. Because computing such spectra by evaluating all contributing Franck-Condon factors becomes impractical (and unnecessary) in large systems, here we propose a time-dependent approach based on Hagedorn wavepacket dynamics. We use Hagedorn functions-products of a Gaussian and carefully generated polynomials-to represent SVL initial states because in systems whose potential is at most quadratic, Hagedorn functions are exact solutions to the time-dependent Schrödinger equation and can be propagated with the same equations of motion as a simple Gaussian wavepacket.
View Article and Find Full Text PDFDynamic protected states in the non-Hermitian system.
Sci Rep
September 2024
School of Science, Jiangnan University, Wuxi, 214122, China.
The non-Hermitian skin effect and nonreciprocal behavior are sensitive to the boundary conditions, which are unique features of non-Hermitian systems. In such systems, eigenenergies can become complex, and all eigenstates tend to localize at the boundary, a phenomenon that contrasts with Hermitian topologies. In this work, we theoretically study the dynamic behavior of the propagation of Gaussian wavepackets inside a non-Hermitian lattice and analyze the self-acceleration process of bulk state or Gaussian wavepackets toward the system's boundary.
View Article and Find Full Text PDFSpatiotemporal optical vortices with controllable radial and azimuthal quantum numbers.
Nat Commun
June 2024
School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, Shanghai, China.
Article Synopsis
- This research investigates optical spatiotemporal vortices with transverse photon orbital angular momentum, focusing on spatiotemporal Laguerre-Gaussian (STLG) wavepackets that display unique topological features.
- The study highlights the differences between conventional optical vortices and STLG wavepackets, which carry a composite transverse OAM with two opposing components.
- Additionally, the research demonstrates a method for converting STLG wavepackets to ST Hermite-Gaussian (STHG) wavepackets, paving the way for advancements in quantum information, photonic topology, and various wave applications.
Controlling Electronic Coherences and the Curvature Induced by the Derivative Coupling at a Conical Intersection: A Quantum Ehrenfest (QuEh) Protocol for Reaction Path Following Application to "Channel 3" Benzene Photochemistry.
J Phys Chem A
July 2024
Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, W12 0BZ London, U.K.
We report a protocol for the implementation of "reaction path following" from a transition state through a conical intersection, including both the path curvature induced by the derivative coupling and the corresponding induced electronic coherences. This protocol focuses on the "central" Gaussian wavepacket (initially unexcited) in the quantum Ehrenfest (QuEh) method. Like the reaction path following, the normal mode corresponding to the imaginary frequency at the transition state is given an initial momentum.
View Article and Find Full Text PDFWant AI Summaries of new PubMed Abstracts delivered to your In-box?
Enter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!