This work describes a geometric framework on molecular reaction dynamics based on the variational principle, where the Schr{ö}dinger equation must be solved to ``see'' how a reaction occurs. First, the mathematical preliminaries are given by discussing the principle of least action and the mountain pass theorem. Second, we discuss the physical preliminaries, including the principle of equivalence for deriving the kinetic energy operator (KEO) and artificial intelligence (AI) techniques to build the potential energy surface (PES) in general spacetime. Moreover, we simplified electromagnetic interactions in curved spacetime within the molecular system and consequently, we are able to construct the nuclear Hamiltonian in nonzero curvature spacetime. This indicates possibility to introduce gauge fields through the curvature, such as additional term in the nuclear KEO near a conical intersection. Third, the single-particle approximation provides a powful ansatz to solve the Schr{ö}dinger equation by variational principle. Thus, one can formulate the variational approaches for either electronic structure or quantum dynamics. In this work, based on previous discussions ({\it Phys. Chem. Chem. Phys.} {\bf 27} (2025), 20397) we unified them by a geometric description, where the geometric phase is naturally introduced. Finally, due to optimization characteristic of the present theory, further discussions on the present theory from optimization insight are also given, including two postulates, generative AI techniques, role of perturbation, and Markov process in optimization.
This work describes a geometric framework on molecular reaction dynamics based on the variational principle, where the Schr{ö}dinger equation must be solved to ``see'' how a reaction occurs. First, the mathematical preliminaries are given by discussing the principle of least action and the mountain pass theorem. Second, we discuss the physical preliminaries, including the principle of equivalence for deriving the kinetic energy operator (KEO) and artificial intelligence (AI) techniques to build the potential energy surface (PES) in general spacetime. Moreover, we simplified electromagnetic interactions in curved spacetime within the molecular system and consequently, we are able to construct the nuclear Hamiltonian in nonzero curvature spacetime. This indicates possibility to introduce gauge fields through the curvature, such as additional term in the nuclear KEO near a conical intersection. Third, the single-particle approximation provides a powful ansatz to solve the Schr{ö}dinger equation by variational principle. Thus, one can formulate the variational approaches for either electronic structure or quantum dynamics. In this work, based on previous discussions ({\it Phys. Chem. Chem. Phys.} {\bf 27} (2025), 20397) we unified them by a geometric description, where the geometric phase is naturally introduced. Finally, due to optimization characteristic of the present theory, further discussions on the present theory from optimization insight are also given, including two postulates, generative AI techniques, role of perturbation, and Markov process in optimization.
ConspectusSimulating molecular dynamics (MD) within a comprehensive quantum framework has been a long-standing challenge in computational chemistry. An exponential scaling of computational cost renders solving the time dependent Schrödinger equation (TDSE) of a molecular Hamiltonian, including both electronic and nuclear degrees of freedom (DOFs), as well as their couplings, infeasible for more than a few DOFs. In the Born-Oppenheimer (BO), or adiabatic, picture, electronic and nuclear parts of the wave function are decoupled and treated separately. Within this framework, the nuclear wave function evolves along potential energy surfaces (PESs) computed as solutions to the electronic Schrödinger equation parametrized in the nuclear DOFs. This approximation, together with increasingly elaborate numerical approaches to solve the nuclear time dependent Schrödinger equation (TDSE), enabled the treatment of up to a few dozens of degrees of freedom (DOFs). However, for particular applications, such as photochemistry, the BO approximation breaks down. In this regime of non-adiabatic dynamics, solving the full molecular problem including electron-nuclear couplings becomes essential, further increasing the complexity of the numerical solution. Although valuable methods such as multiconfigurational time-dependent Hartree (MCTDH) have been proposed for the solution of the coupled electron-nuclear dynamics, they remain hampered by an exponential scaling in the number of nuclear DOFs and by the difficulty of finding universal variational forms.In this Account, we present a perspective on novel quantum computational algorithms, aiming to alleviate the exponential scaling inherent to the simulation of many-body quantum dynamics. In particular, we focus on the derivation and application of quantum algorithms for adiabatic and non-adiabatic quantum dynamics, which include efficient approaches for the calculation of the BO potential energy surfaces (PESs). Thereafter, we study the time-evolution of a model system consisting of two coupled PESs in first and second quantization. In a first application, we discuss a recently introduced quantum algorithm for the evolution of a wavepacket in first quantization and exploit the potential quantum advantage of mapping its spatial grid representation to logarithmically many qubits. For the second demonstration, we move to the second quantization framework and review the scaling properties of two alternative time-evolution algorithms, namely, a variational quantum algorithm (VQA) (based on the McLachlan variational principle) and conventional Trotter-type evolution (based on a Lie-Trotter-Suzuki formula). Both methods clearly demonstrate the potential of quantum algorithms and their favorable scaling compared to the available classical approaches. However, a clear demonstration of quantum advantage in the context of molecular quantum dynamics may require the implementation of these algorithms in fault-tolerant quantum computers, while their application in near-term, noisy quantum devices is still unclear and deserves further investigation.
In this work, we explore the implementation possibility of quantum simulation for quantum molecular dynamics, in particular for reaction dynamics, though several implementations have already reported through quantum-classical mixed simulations ({\it Acc. Chem. Res.} {\bf 54} (2021), 4229 and {\it J. Phys. Chem. Lett.} {\bf xx} (2026), XXXX). To analyze this aspect, we examine (1) the conjugacy relation between quantum simulator and the target molecular system, (2) the wave function correspondence in quantum algorithm and classical algorithm for multi-dimensional dynamics, (3) problems arisen from real-valued classical algorithms, and finally (4) geometric phase arisen from the separation among the degrees of freedom (DOFs). As is well known, the aforementioned first and second points play fundamental roles in quantum simulation of quantum many-body systems, and the third and fourth points are theoretical issues that might introduce problems in classical and quantum computing. In this work, we mainly focus on the third and fourth points by analysis of the first two points by reviewing previously reported quantum-classical mixed implementations of quantum simulation. We also consider gauge freedom in high-dimensional quantum molecular dynamics that has been introduced recently, and then discuss possibility of advantages and disadvantages of quantum simulation for molecular reaction dynamics.
The symmetry properties of quantum variational principles are considered. Euler-Poincaré reduction theory is applied to the Dirac-Frenkel variational principle for Schrödinger's dynamics producing new variational principles for the different pictures of quantum mechanics (Schrödinger's, Heisenberg's and Dirac's) as well as for the Wigner-Moyal formulation. In addition, new variational principles for mixed states dynamics have been formulated. The already known geometric characterization of quantum mechanics on the complex projective space is shown to emerge naturally from the Euler-Poincaré variational principle. Semidirect-product structures are seen to produce new variational principles for Dirac's picture. \n \nOn the other hand, the variational and Hamiltonian approach to Ehrenfest expectation values dynamics is proposed. In the Lagrangian framework, the quantum variational principle for Schrödinger's dynamics is extended to account for both classical and quantum degrees of freedom. First, it is shown that the mean field model of any quantum mechanical system can be derived from a classical-quantum Euler-Poincaré Lagrangian on the direct sum Lie algebra of the Heisenberg and unitary groups. Then, the semidirect-product structure (named Ehrenfest group), is constructed using the displacement operator from the theory of coherent quantum states ( the unitary action of the Heisenberg group on the space of wavefunctions). New classical-quantum equations for Ehrenfest's expectation values dynamics are derived redefining the mean-field model Euler-Poincaré Lagrangian on the Lie algebra of the Ehrenfest group. In the Hamiltonian framework, first expectation values of the canonical observables are shown to be equivariant momentum maps for the unitary action of the Heisenberg group on quantum states. Then, the Hamiltonian structure for Ehrenfest's dynamics is shown to be Lie-Poisson for the Ehrenfest group. The variational formulation is then given a corresponding Hamiltonian structure. The classical-quantum Ehrenfest dynamics equations produce classical and quantum dynamics as special limit cases. \n \nIn the particular case of Gaussian states, expectation values couple to second order moments, so that GS are completely characterized by first and second moments. When the total energy is computed with respect to a Gaussian state, higher moments can be expressed in terms of the first two, so that the moment hierarchy closes for Gaussian states. Second moments are shown to be equivariant momentum maps for the action of the symplectic group on the space of Wigner functions. Eventually, Gaussian states are shown to possess a Lie-Poisson structure on the Jacobi group. This structure produces an energy-conserving variant of a class of Gaussian moment models that have appeared in the chemical physics literature.
The applicability of the time-dependent Schrödinger equation to molecular reaction dynamics is discussed from several perspectives. First, with reference to Schrödinger's original paper, we point out that he himself stated that what is now called the time-dependent Schrödinger equation is insufficient to solve the problems of non-conservative systems like chemical reactions. Second, the detonation velocity of some explosives has been demonstrated as a new parameter that directly illustrates the timescales of elementary chemical reactions. The calculated values of ~30 fs are in good agreement with the timescale of 10–100 fs for chemical reactions, as experimentally demonstrated by the Nobel laureate Zewail. Third, we consider the state of an electron in a box when the length of the box suddenly increases from L to 2L, which is often found in textbooks as the simplest example of the application of the time-dependent Schrödinger equation to chemical reactions, and highlight the inconsistency between this model and an actual chemical reaction. Finally, we re-iterate that only quantum steady states can be accurately described by the Schrödinger equation and that it is the most natural to connect changes between adjacent quantum steady states by the variational principle.
Quantum computers offer the potential to simulate chemical systems beyond the capabilities of classical computers. However, current quantum computing approaches to ab initio molecular dynamics (AIMD) and vibrational analysis largely rely on the Variational Quantum Eigensolver (VQE), whose circuit depth, iterative optimization, and measurement overhead hinder the practical simulation of molecular systems. Here, we extend Sample-based Quantum Diagonalization (SQD), a quantum sampling-based electronic-structure method, to AIMD and vibrational analysis on noisy quantum hardware. The resulting framework enables nuclear forces and Hessians to be derived from quantum-computed electronic energies, facilitating molecular dynamics and characterization of potential energy surfaces. We demonstrate the approach using diverse molecular processes, including H2 bond vibration, the proton-exchange reaction (H+ + H2), the bimolecular nucleophilic substitution (SN2) reaction (Cl- + CH3Cl), and proton-transfer dynamics between two water molecules involving transient formation of the Zundel cation (H2O···H+···OH2). Across these systems, SQD-AIMD trajectories and vibrational frequencies obtained with noisy quantum hardware/noise-models show excellent agreement with classical reference calculations and substantially improved accuracy relative to previously reported VQE-based simulations. Importantly, the SQD-based approach eliminates costly variational optimization and significantly reduces measurement requirements relative to VQE, enabling AIMD simulations spanning approximately 500 timesteps (∼200 fs), substantially longer than previously reported AIMD simulations on quantum hardware. These results establish SQD-AIMD and SQD-vibrational analysis as practical and scalable frameworks for simulating chemical reaction dynamics and elucidating reaction mechanisms on near-term quantum hardware, extending quantum computing applications beyond static electronic structure calculations.
Abstract Time-dependent variational principle (TDVP) provides powerful methods in solving the time-dependent Schröinger equation. As such Kan developed a TDVP (Kan 1981 Phys. Rev. A 24, 2831) and found that there is no Legendre transformation in quantum variational principle, suggesting that there is no place for the Maupertuis reduced action to appear in quantum dynamics. This claim is puzzling for the study of quantum–classical correspondence, since the Maupertuis least action principle practically sets the very basic foundation of classical mechanics. Zambrini showed within the theory of stochastic calculus of variations that the Maupertuis least action principle can lead to the Nelson stochastic quantization theory (Zambrini 1984 J. Math. Phys. 25, 1314). We here revisit the basic aspect of TDVP and reveal the hidden roles of Maupertuis-Hamilton least action in the Schrödinger wavepacket dynamics. On this basis we propose a dual least (stationary) action principle, which is composed of two variational functionals; one responsible for ‘energy related dynamics’ and the other for ‘dynamics of wave-flow’. The former is mainly a manifestation of particle nature in wave-particle duality, while the latter represents that of matter wave. It is also shown that by representing the TDVP in terms of these inseparably linked variational functionals the problem of singularity, which is inherent to the standard TDVPs, is resolved. The structure and properties of this TDVP are also discussed.