Quantum tunneling dynamics in multidimensional systems: A matching-pursuit description

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.

Grid-based methods for chemistry simulations on a quantum computer

First-quantized, grid-based methods for chemistry modeling are a natural and elegant fit for quantum computers. However, it is infeasible to use today's quantum prototypes to explore the power of this approach because it requires a substantial number of near-perfect qubits. Here, we use exactly emulated quantum computers with up to 36 qubits to execute deep yet resource-frugal algorithms that model 2D and 3D atoms with single and paired particles. A range of tasks is explored, from ground state preparation and energy estimation to the dynamics of scattering and ionization; we evaluate various methods within the split-operator QFT (SO-QFT) Hamiltonian simulation paradigm, including protocols previously described in theoretical papers and our own techniques. While we identify certain restrictions and caveats, generally, the grid-based method is found to perform very well; our results are consistent with the view that first-quantized paradigms will be dominant from the early fault-tolerant quantum computing era onward.

Functional Tensor-Train Chebyshev Method for Multidimensional Quantum Dynamics Simulations

Methods for efficient simulations of multidimensional quantum dynamics are essential for theoretical studies of chemical systems where quantum effects are important, such as those involving rearrangements of protons or electronic configurations. Here, we introduce the functional tensor-train Chebyshev (FTTC) method for rigorous nuclear quantum dynamics simulations. FTTC is essentially the Chebyshev propagation scheme applied to the initial state represented in a continuous analogue tensor-train format. We demonstrate the capabilities of FTTC as applied to simulations of proton quantum dynamics in a 50-dimensional model of hydrogen-bonded DNA base pairs.

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.

Two-Dimensional Vibrational Spectroscopy of a Dissipative System with the Optimized Mean-Trajectory Approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical method for computing vibrational response functions from action-quantized classical trajectories connected by discrete transitions representing radiation-matter interactions. Here we apply this method to an anharmonic chromophore coupled to a harmonic bath. A forward-backward trajectory implementation of the OMT method is described that addresses the numerical challenges of applying the OMT to large systems with disparate frequency scales. The OMT is shown to well reproduce line shapes and waiting time dynamics in the pure dephasing limit of weak coupling to an off-resonant bath. The OMT is also shown to describe a case where energy transfer is the predominant source of line broadening.

Steered quantum dynamics for energy minimization

We introduce a quantum optimal control algorithm for energy minimization that combines the diffeomorphic modulation under observable response preserving homotopy (D-MORPH) gradient and the Broyden Fletcher Goldfarb Shanno (BFGS) iterative scheme for nonlinear optimization. An extended set of controls defining the time-dependent mass, dipole moment, and external perturbational field are optimized to find an effective Hamiltonian that steers the dynamics of the system into the global minimum without getting trapped into local minima. The algorithm is illustrated as applied to energy minimization on rugged surfaces and golf potentials comparable to those previously explored for testing quantum annealing methodologies.

Two-dimensional spectroscopy of coupled vibrations with the optimized mean-trajectory approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical representation of the nonlinear vibrational response function used to compute multidimensional infrared spectra. In this method, response functions are calculated from a sequence of classical trajectories linked by discontinuities representing the effects of radiation-matter interactions, thus providing an approximation to quantum dynamics using classical inputs. This approach was previously formulated and assessed numerically for a single anharmonic degree of freedom. Our previous work is generalized here in two respects. First, the derivation of the OMT is extended to any number of coupled anharmonic vibrations by determining semiclassical approximations for pairs of double-sided Feynman diagrams. Second, an efficient numerical procedure is developed for calculating two-dimensional infrared spectra of coupled anharmonic vibrations in the OMT approximation. The OMT approximation is shown to reproduce the fundamental features of the quantum response function including both coherence and population dynamics.

Broadband infrared and Raman probes of excited-state vibrational molecular dynamics: simulation protocols based on loop diagrams

Vibrational motions in electronically excited states can be observed either by time and frequency resolved infrared absorption or by off resonant stimulated Raman techniques. Multipoint correlation function expressions are derived for both signals. Three representations which suggest different simulation protocols for the signals are developed. These are based on the forward and the backward propagation of the wavefunction, sum over state expansion using an effective vibrational Hamiltonian or a semiclassical treatment of a bath. We show that the effective temporal (Δt) and spectral (Δω) resolution of the techniques is not controlled solely by experimental knobs but also depends on the system dynamics being probed. The Fourier uncertainty ΔωΔt > 1 is never violated.

An optimized semiclassical approximation for vibrational response functions

The observables of multidimensional infrared spectroscopy may be calculated from nonlinear vibrational response functions. Fully quantum dynamical calculations of vibrational response functions are generally impractical, while completely classical calculations are qualitatively incorrect at long times. These challenges motivate the development of semiclassical approximations to quantum mechanics, which use classical mechanical information to reconstruct quantum effects. The mean-trajectory (MT) approximation is a semiclassical approach to quantum vibrational response functions employing classical trajectories linked by deterministic transitions representing the effects of the radiation-matter interaction. Previous application of the MT approximation to the third-order response function R(3)(t3, t2, t1) demonstrated that the method quantitatively describes the coherence dynamics of the t3 and t1 evolution times, but is qualitatively incorrect for the waiting-time t2 period. Here we develop an optimized version of the MT approximation by elucidating the connection between this semiclassical approach and the double-sided Feynman diagrams (2FD) that represent the quantum response. Establishing the direct connection between 2FD and semiclassical paths motivates a systematic derivation of an optimized MT approximation (OMT). The OMT uses classical mechanical inputs to accurately reproduce quantum dynamics associated with all three propagation times of the third-order vibrational response function.

Signatures of discrete breathers in coherent state quantum dynamics

In classical mechanics, discrete breathers (DBs) - a spatial time-periodic localization of energy - are predicted in a large variety of nonlinear systems. Motivated by a conceptual bridging of the DB phenomena in classical and quantum mechanical representations, we study their signatures in the dynamics of a quantum equivalent of a classical mechanical point in phase space - a coherent state. In contrast to the classical point that exhibits either delocalized or localized motion, the coherent state shows signatures of both localized and delocalized behavior. The transition from normal to local modes have different characteristics in quantum and classical perspectives. Here, we get an insight into the connection between classical and quantum perspectives by analyzing the decomposition of the coherent state into system's eigenstates, and analyzing the spacial distribution of the wave-function density within these eigenstates. We find that the delocalized and localized eigenvalue components of the coherent state are separated by a mixed region, where both kinds of behavior can be observed. Further analysis leads to the following observations. Considered as a function of coupling, energy eigenstates go through avoided crossings between tunneling and non-tunneling modes. The dominance of tunneling modes in the high nonlinearity region is compromised by the appearance of new types of modes - high order tunneling modes - that are similar to the tunneling modes but have attributes of non-tunneling modes. Certain types of excitations preferentially excite higher order tunneling modes, allowing one to study their properties. Since auto-correlation functions decrease quickly in highly nonlinear systems, short-time dynamics are sufficient for modeling quantum DBs. This work provides a foundation for implementing modern semi-classical methods to model quantum DBs, bridging classical and quantum mechanical signatures of DBs, and understanding spectroscopic experiments that involve a coherent state.

Quantized Hamilton dynamics describes quantum discrete breathers in a simple way

We study the localization of energy in a nonlinear coupled system, exhibiting so-called breather modes, using quantized Hamilton dynamics (QHD). Already at the lowest order, which is only twice as complex as classical mechanics, this simple semiclassical method incorporates quantum-mechanical effects. The transition between the localized and delocalized regimes is instantaneous in classical mechanics, while it is gradual due to tunneling in both quantum mechanics and QHD. In contrast to classical mechanics, which predicts an abrupt appearance of breathers, quantum mechanics and QHD show an alternation of localized and delocalized behavior in the transient region. QHD includes zero-point energy that is reflected in a shifted energy asymptote for the localized states, providing another improvement on the classical perspective. By detailed analysis of the distribution and transfer of energy within classical mechanics, QHD, and quantum dynamics, we conclude that QHD is an efficient approach that accounts for moderate quantum effects and can be used to identify quantum breathers in large nonlinear systems.

Modeling vibrational resonance in linear hydrocarbon chain with a mixed quantum-classical method

The quantum dynamics of a vibrational excitation in a linear hydrocarbon model system is studied with a new mixed quantum-classical method. The method is suited to treat many-body systems consisting of a low dimensional quantum primary part coupled to a classical bath. The dynamics of the primary part is governed by the quantum corrected propagator, with the corrections defined in terms of matrix elements of zeroth order propagators. The corrections are taken to the classical limit by introducing the frozen Gaussian approximation for the bath degrees of freedom. The ability of the method to describe dynamics of multidimensional systems has been tested. The results obtained by the method have been compared to previous quantum simulations performed with the quasiadiabatic path integral method.

Approximating the ground state of fermion systems by multiple determinant states: matching pursuit approach

We present a simple and stable numerical method to approximate the ground state of a quantum many-body system by multiple determinant states. This method searches these determinant states one by one according to the matching pursuit algorithm. The first determinant state is identical to that of the Hartree-Fock theory. Calculations for a two-dimensional Hubbard model serve as a demonstration.

On the properties of a primitive semiclassical surface hopping propagator for nonadiabatic quantum dynamics

A previously developed nonadiabatic semiclassical surface hopping propagator [M. F. Herman J. Chem. Phys. 103, 8081 (1995)] is further studied. The propagator has been shown to satisfy the time-dependent Schrodinger equation (TDSE) through order h, and the O(h2) terms are treated as small errors, consistent with standard semiclassical analysis. Energy is conserved at each hopping point and the change in momentum accompanying each hop is parallel to the direction of the nonadiabatic coupling vector resulting in both transmission and reflection types of hops. Quantum mechanical analysis and numerical calculations presented in this paper show that the h2 terms involving the interstate coupling functions have significant effects on the quantum transition probabilities. Motivated by these data, the h2 terms are analyzed for the nonadiabatic semiclassical propagator. It is shown that the propagator can satisfy the TDSE for multidimensional systems by including another type of nonclassical trajectories that reflect on the same surfaces. This h2 analysis gives three conditions for these three types of trajectories so that their coefficients are uniquely determined. Besides the nonadiabatic semiclassical propagator, a numerically useful quantum propagator in the adiabatic representation is developed to describe nonadiabatic transitions.

Including quantum effects in the dynamics of complex (i.e., large) molecular systems

The development in the 1950s and 1960s of crossed molecular beam methods for studying chemical reactions at the single-collision molecular level stimulated the need and desire for theoretical methods to describe these and other dynamical processes in molecular systems. Chemical dynamics theory has made great strides in the ensuing decades, so that methods are now available for treating the quantum dynamics of small molecular systems essentially completely. For the large molecular systems that are of so much interest nowadays (e.g., chemical reactions in solution, in clusters, in nanostructures, in biological systems, etc.), however, the only generally available theoretical approach is classical molecular dynamics (MD) simulations. Much effort is currently being devoted to the development of approaches for describing the quantum dynamics of these complex systems. This paper reviews some of these approaches, especially the use of semiclassical approximations for adding quantum effects to classical MD simulations, also showing some new versions that should make these semiclassical approaches even more practical and accurate.

Matching-pursuit/split-operator-Fourier-transform simulations of excited-state nonadiabatic quantum dynamics in pyrazine

A simple approach for numerically exact simulations of nonadiabatic quantum dynamics in multidimensional systems is introduced and applied to the description of the photoabsorption spectroscopy of pyrazine. The propagation scheme generalizes the recently developed matching-pursuit/split-operator-Fourier-transform (MP/SOFT) method [Y. Wu and V. S. Batista, J. Chem. Phys. 121, 1676 (2004)] to simulations of nonadiabatic quantum dynamics. The time-evolution operator is applied, as defined by the Trotter expansion to second order accuracy, in dynamically adaptive coherent-state expansions. These representations are obtained by combining the matching-pursuit algorithm with a gradient-based optimization method. The accuracy and efficiency of the resulting computational approach are demonstrated in calculations of time-dependent survival amplitudes and photoabsorption cross sections, using a model Hamiltonian that allows for direct comparisons with benchmark calculations. Simulations in full-dimensional potential energy surfaces involve the propagation of a 24-dimensional wave packet to describe the S(1)S(2) interconversion of pyrazine after S(0)-->S(2) photoexcitation. The reported results show that the generalized MP/SOFT method is a practical and accurate approach to model nonadiabatic reaction dynamics in polyatomic systems.

Matching-pursuit split-operator Fourier-transform simulations of excited-state intramolecular proton transfer in 2-(2′-hydroxyphenyl)-oxazole

The excited-state intramolecular proton-transfer dynamics associated with the keto-enolic tautomerization reaction in 2-(2(')-hydroxyphenyl)-oxazole is simulated according to a numerically exact quantum-dynamics propagation method and a full-dimensional excited-state potential energy surface, based on an ab initio reaction surface Hamiltonian. The reported simulations involve the propagation of 35-dimensional wave packets according to the recently developed matching-pursuit/split-operator-Fourier-transform (MP/SOFT) method by Wu and Batista. The underlying propagation scheme recursively applies the time-evolution operator as defined by the Trotter expansion to second order accuracy in dynamically adaptive coherent-state expansions. Computations of time-dependent survival amplitudes, photoabsorption cross sections, and time-dependent reactant(product) populations are compared to the corresponding calculations based on semiclassical approaches, including the Herman-Kluk semiclassical initial value representation method. The reported results demonstrate the capabilities of the MP/SOFT method as a valuble computational tool to study ultrafast reaction dynamics in polyatomic systems as well as to validate semiclassical simulations of complex (nonintegrable) quantum dynamics in multidimensional model systems.

A version of diffusion Monte Carlo method based on random grids of coherent states. II. Six-dimensional simulation of electronic states of H2

We report a new version of the diffusion Monte Carlo (DMC) method, based on coherent-state quantum mechanics. Randomly selected grids of coherent states in phase space are used to obtain numerical imaginary time solutions of the Schrodinger equation, with an iterative refinement technique to improve the quality of the Monte Carlo grid. Accurate results were obtained, for the appropriately symmetrized two lowest states of the hydrogen molecule, by Monte Carlo sampling and six-dimensional propagation in the full phase space.

Quantum dynamics of complex molecular systems

This Perspective presents a broad overview of the present status of theoretical capabilities for describing quantum dynamics in molecular systems with many degrees of freedom, e.g., chemical reactions in solution, clusters, solids, or biomolecular environments.

Matching-pursuit∕split-operator Fourier-transform simulations of nonadiabatic quantum dynamics

A rigorous and practical approach for simulations of nonadiabatic quantum dynamics is introduced. The algorithm involves a natural extension of the matching-pursuitsplit-operator Fourier-transform (MPSOFT) method [Y. Wu and V. S. Batista, J. Chem. Phys. 121, 1676 (2004)] recently developed for simulations of adiabatic quantum dynamics in multidimensional systems. The MPSOFT propagation scheme, extended to nonadiabatic dynamics, recursively applies the time-evolution operator as defined by the standard perturbation expansion to first-, or second-order, accuracy. The expansion is implemented in dynamically adaptive coherent-state representations, generated by an approach that combines the matching-pursuit algorithm with a gradient-based optimization method. The accuracy and efficiency of the resulting propagation method are demonstrated as applied to the canonical model systems introduced by Tully for testing simulations of dual curve-crossing nonadiabatic dynamics.

Matching-pursuit/split-operator-Fourier-transform computations of thermal correlation functions

A rigorous and practical methodology for evaluating thermal-equilibrium density matrices, finite-temperature time-dependent expectation values, and time-correlation functions is described. The method involves an extension of the matching-pursuit/split-operator-Fourier-transform method to the solution of the Bloch equation via imaginary-time propagation of the density matrix and the evaluation of Heisenberg time-evolution operators through real-time propagation in dynamically adaptive coherent-state representations.