Linear dependence and energy conservation in Gaussian wavepacket basis sets

Assessing the Accuracy of Quantum Dynamics Performed in the Time-Dependent Basis Representation

A full quantum-mechanical (QM) description of large amplitude nuclear motion, associated with chemical reactions or isomerization of high-dimensional molecular systems, is inherently challenging due to the exponential scaling of the QM complexity with system size. To ameliorate the scaling bottleneck in studies of realistic systems, typically modeled in the configuration space, the nuclear wave functions are represented in terms of time-dependent basis functions. Such bases are expected to give an accurate description with a modest number of basis functions employed, by adapting them to the wave function solving the time-dependent Schrödinger equation. It is not, however, straightforward to estimate the accuracy of the resulting solution: in QM the energy conservation, a convenient such measure for a classical trajectory evolving in a time-independent potential, is not a sufficient criterion of the dynamics' accuracy. In this work, we argue that the expectation value of the Hamiltonian's "variance", quantifying the basis completeness, is a suitable practical measure of the quantum dynamics' accuracy. Illustrations are given for several chemistry-relevant test systems, modeled employing time-independent as well as time-dependent bases, including the coupled and variational coherent states methods and the quantum-trajectory guided adaptable Gaussians (QTAG) as the latter basis type. A novel semilocal definition of the QTAG basis time-evolution for placing the basis functions "in the right place at the right time" is also presented.

Quantum Dynamics from Classical Trajectories

Nonadiabatic molecular dynamics plays an essential role in exploring the time evolution of molecular systems. Various methods have been developed for this study, with varying accuracy and computational cost. One very successful among them is trajectory surface hopping, which propagates nuclei as classical trajectories using forces from a quantum description of the electrons and incorporates nonadiabatic effects through stochastic state changes during each trajectory propagation. A statistical analysis of an ensemble of the independent trajectories recovers the simulated system's behavior. This approach can give good results, but it is known to overlook nuclear quantum effects, leading to inaccurate predictions. Here, we present quantum dynamics from classical trajectories (QDCT), a new protocol to recover the quantum wavepacket from the classical trajectories generated by surface hopping. In this first QDCT implementation, we apply it to recover results at the multiple spawning level from postprocessing surface hopping precomputed trajectories. With a series of examples, we demonstrate QDCT's potential to improve the accuracy of the dynamics, correct decoherence effects, and diagnose problems or increase confidence in surface hopping results. All that comes at virtually no computational cost since no new electronic calculation is required.

Probing the sensitivity of ab initio multiple spawning to its parameters

Full multiple spawning (FMS) offers a strategy to simulate the nonadiabatic dynamics of molecular systems by describing their nuclear wavefunctions by a linear combination of coupled trajectory basis functions (TBFs). Applying a series of controlled approximations to the full multiple spawning (FMS) equations leads to the ab initio multiple spawning (AIMS), which is compatible with an on-the-fly propagation of the TBFs and an accurate description of nonadiabatic processes. The AIMS strategy and its numerical implementations, however, rely on a series of user-defined parameters. Herein, we investigate the influence of these parameters on the electronic-state population of two molecular systems- trans-azomethane and a two-dimensional model of the butatriene cation. This work highlights the stability of AIMS with respect to most of its parameters, underlines the specific parameters that require particular attention from the user of the method, and offers prescriptions for an informed selection of their value.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00214-023-03004-w.

Variational Approach for Linearly Dependent Moving Bases in Quantum Dynamics: Application to Gaussian Functions

In this paper, we present a variational treatment of the linear dependence for a non-orthogonal time-dependent basis set in solving the Schrödinger equation. The method is based on (i) the definition of a linearly independent working space and (ii) a variational construction of the propagator over finite time steps. The second point allows the method to properly account for changes in the dimensionality of the working space along the time evolution. In particular, the time evolution is represented by a semi-unitary transformation. Tests are carried out on a quartic double-well potential with Gaussian basis functions whose centers evolve according to classical equations of motion. We show that the resulting dynamics converges to the exact one and is unitary by construction.

Controlling energy conservation in quantum dynamics with independently moving basis functions: Application to multi-configuration Ehrenfest

Application of the time-dependent variational principle to a linear combination of frozen-width Gaussians describing the nuclear wavefunction provides a formalism where the total energy is conserved. The computational downside of this formalism is that trajectories of individual Gaussians are solutions of a coupled system of differential equations, limiting implementation to serial propagation algorithms. To allow for parallelization and acceleration of the computation, independent trajectories based on simplified equations of motion were suggested. Unfortunately, within practical realizations involving finite Gaussian bases, this simplification leads to breaking the energy conservation. We offer a solution for this problem by using Lagrange multipliers to ensure the energy and norm conservation regardless of basis function trajectories or basis completeness. We illustrate our approach within the multi-configurational Ehrenfest method considering a linear vibronic coupling model.

Predicting Molecular Photochemistry Using Machine-Learning-Enhanced Quantum Dynamics Simulations

The processes which occur after molecules absorb light underpin an enormous range of fundamental technologies and applications, including photocatalysis to enable new chemical transformations, sunscreens to protect against the harmful effects of UV overexposure, efficient photovoltaics for energy generation from sunlight, and fluorescent probes to image the intricate details of complex biomolecular structures. Reflecting this broad range of applications, an enormously versatile set of experiments are now regularly used to interrogate light-driven chemical dynamics, ranging from the typical ultrafast transient absorption spectroscopy used in many university laboratories to the inspiring central facilities around the world, such as the next-generation of X-ray free-electron lasers.Computer simulations of light-driven molecular and material dynamics are an essential route to analyzing the enormous amount of transient electronic and structural data produced by these experimental sources. However, to date, the direct simulation of molecular photochemistry remains a frontier challenge in computational chemical science, simultaneously demanding the accurate treatment of molecular electronic structure, nuclear dynamics, and the impact of nonadiabatic couplings.To address these important challenges and to enable new computational methods which can be integrated with state-of-the-art experimental capabilities, the past few years have seen a burst of activity in the development of "direct" quantum dynamics methods, merging the machine learning of potential energy surfaces (PESs) and nonadiabatic couplings with accurate quantum propagation schemes such as the multiconfiguration time-dependent Hartree (MCTDH) method. The result of this approach is a new generation of direct quantum dynamics tools in which PESs are generated in tandem with wave function propagation, enabling accurate "on-the-fly" simulations of molecular photochemistry. These simulations offer an alternative route toward gaining quantum dynamics insights, circumventing the challenge of generating ab initio electronic structure data for PES fitting by instead only demanding expensive energy evaluations as and when they are needed.In this Account, we describe the chronological evolution of our own contributions to this field, focusing on describing the algorithmic developments that enable direct MCTDH simulations for complex molecular systems moving on multiple coupled electronic states. Specifically, we highlight active learning strategies for generating PESs during grid-based quantum chemical dynamics simulations, and we discuss the development and impact of novel diabatization schemes to enable direct grid-based simulations of photochemical dynamics; these developments are highlighted in a series of benchmark molecular simulations of systems containing multiple nuclear degrees of freedom moving on multiple coupled electronic states. We hope that the ongoing developments reported here represent a major step forward in tools for modeling excited-state chemistry such as photodissociation, proton and electron transfer, and ultrafast energy dissipation in complex molecular systems.

Computational approaches to efficient generation of the stationary state for incoherent light excitation

Light harvesting processes are often computationally studied from a time-dependent viewpoint, in line with ultrafast coherent spectroscopy experiments. Yet, natural processes take place in the presence of incoherent light, which induces a stationary state. Such stationary states can be described using the eigenbasis of the molecular Hamiltonian, but for realistic systems, a full diagonalization is prohibitively expensive. We propose three efficient computational approaches to obtain the stationary state that circumvents system Hamiltonian diagonalization. The connection between the incoherent perturbations, decoherence, and Kraus operators is established.

Multidimensional Tunneling Dynamics Employing Quantum-Trajectory Guided Adaptable Gaussian Bases

An efficient basis representation of time-dependent wavefunctions is essential for theoretical studies of high-dimensional molecular systems exhibiting large-amplitude motion. For fully coupled anharmonic systems, the complexity of a general wavefunction scales exponentially with the system size; therefore, for practical reasons, it is desirable to adapt the basis to the time-dependent wavefunction at hand. Often times on this quest for a minimal basis representation, time-dependent Gaussians are employed, in part because of their localization in both configuration and momentum spaces and also because of their direct connection to classical and semiclassical dynamics, guiding the evolution of the basis function parameters. In this work, the quantum-trajectory guided adaptable Gaussian (QTAG) bases method [ J. Chem. Theory Comput. 2020, 16, 18-34] is generalized to include correlated, i.e., non-factorizable, basis functions, and the performance of the QTAG dynamics is assessed on benchmark system/bath tunneling models of up to 20 dimensions. For the popular choice of initial conditions describing tunneling between the reactant/product wells, the minimal "semiclassical" description of the bath modes using essentially a single multidimensional basis function combined with the multi-Gaussian representation of the tunneling mode is shown to capture the dominant features of dynamics in a highly efficient manner.

Computation of the S1 ← S0 Vibronic Absorption Spectrum of Formaldehyde by Variational Gaussian Wavepacket and Semiclassical IVR Methods

The vibronic absorption spectrum of the electric dipole forbidden and vibronically allowed S₁(¹ A₂) ← S₀(¹ A₁) transition of formaldehyde is calculated by Gaussian wavepacket and semiclassical methods, along with numerically exact reference calculations, using the potential energy surface of Fu, Shepler, and Bowman ( J. Am. Chem. Soc. 2011, 133, 7957). Specifically, the variational multiconfigurational Gaussian (vMCG) approach and the Herman-Kluk semiclassical initial value representation (HK-SCIVR) are compared to assess the accuracy and convergence of these methods, benchmarked against numerically exact time-dependent wavepacket propagation (TDWP) on the reference potential energy surface. The vMCG calculation is shown to converge quite well with about 100 variationally evolving Gaussian functions and using a local cubic expansion instead of the conventional local harmonic approximation. By contrast, the HK-SCIVR approach with ∼10⁵ trajectories reproduces the vibrationally structured spectral envelope correctly but yields a strongly broadened spectrum. The comparison of the computed absorption spectrum with experiment shows that the relevant vibronic progressions are reasonably reproduced by all computations, but deviations of the order of 10-100 cm^-1 occur, underscoring that both electronic structure calculations and dynamical approaches remain challenging in the calculation of typical small-molecule excited-state spectra by trajectory-based methods.

Variational nonadiabatic dynamics in the moving crude adiabatic representation: Further merging of nuclear dynamics and electronic structure

A new methodology of simulating nonadiabatic dynamics using frozen-width Gaussian wavepackets within the moving crude adiabatic representation with the on-the-fly evaluation of electronic structure is presented. The main feature of the new approach is the elimination of any global or local model representation of electronic potential energy surfaces; instead, the electron-nuclear interaction is treated explicitly using the Gaussian integration. As a result, the new scheme does not introduce any uncontrolled approximations. The employed variational principle ensures the energy conservation and leaves the number of electronic and nuclear basis functions as the only parameter determining the accuracy. To assess performance of the approach, a model with two electronic and two nuclear spacial degrees of freedom containing conical intersections between potential energy surfaces has been considered. Dynamical features associated with nonadiabatic transitions and nontrivial geometric (or Berry) phases were successfully reproduced within a limited basis expansion.

Multiconfigurational quantum propagation with trajectory-guided generalized coherent states

A generalized version of the coupled coherent states method for coherent states of arbitrary Lie groups is developed. In contrast to the original formulation, which is restricted to frozen-Gaussian basis sets, the extended method is suitable for propagating quantum states of systems featuring diversified physical properties, such as spin degrees of freedom or particle indistinguishability. The approach is illustrated with simple models for interacting bosons trapped in double- and triple-well potentials, most adequately described in terms of SU(2) and SU(3) bosonic coherent states, respectively.