Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2103.11783 · Journal of Mathematical Physics

Various variational approximations of quantum dynamics

Caroline Lasser, Chunmei Su

Published Jul 1, 2022• 12 citations
arXivHTMLPDF

Abstract

We investigate variational principles for the approximation of quantum dynamics that apply for approximation manifolds that do not have complex linear tangent spaces. The first one, dating back to McLachlan [Mol. Phys. 8, 39–44 (1964)], minimizes the residuum of the time-dependent Schrödinger equation, while the second one, originating from the lecture notes of Kramer and Saraceno [Geometry of the Time-Dependent Variational Principle in Quantum Mechanics, Lecture Notes in Physics Vol. 140 (Springer, Berlin, 1981)], imposes the stationarity of an action functional. We characterize both principles in terms of metric and symplectic orthogonality conditions, consider their conservation properties, and derive an elementary a posteriori error estimate. As an application, we revisit the time-dependent Hartree approximation and frozen Gaussian wave packets.

DOI
10.1063/5.0088265

Related papers

On variational approximations in quantum molecular dynamics

Christian Lubich

The Dirac-Frenkel-McLachlan variational principle is the basic tool for obtaining computationally accessible approximations in quantum molecular dynamics. It determines equations of motion for an approximate time-dependent wave function on an approximation manifold of reduced dimension. This paper gives a near-optimality result for variational approximations. It bounds the error in terms of the distance of the exact wave function to the approximation manifold and identifies the parameters that control the deviation of the variational approximation from the best approximation on the manifold.

Mathematics of Computation•May 25, 2004•54 citations

Local-in-Time Error in Variational Quantum Dynamics

Rocco Martinazzo, Irène Burghardt

The McLachlan "minimum-distance" principle for optimizing approximate solutions of the time-dependent Schrödinger equation is revisited, with a focus on the local-in-time error accompanying the variational solutions. Simple, exact expressions are provided for this error, which are then evaluated in illustrative cases, notably the widely used mean-field approach and the adiabatic quantum molecular dynamics. Based on these findings, we demonstrate the rigorous formulation of an adaptive scheme that resizes on the fly the underlying variational manifold and, hence, optimizes the overall computational cost of a quantum dynamical simulation. Such adaptive schemes are a crucial requirement for devising and applying direct quantum dynamical methods to molecular and condensed-phase problems.

arXiv:1907.00841•Physical Review Letters•Apr 13, 2020•34 citations

Chaos in time-dependent variational approximations to quantum dynamics

Fred Cooper, John F. Dawson, Salman Habib +1 more

Dynamical chaos has recently been shown to exist in the Gaussian approximation in quantum mechanics and in the self-consistent mean field approach to studying the dynamics of quantum fields. In this study, we first note that any variational approximation to the dynamics of a quantum system based on the Dirac action principle leads to a classical Hamiltonian dynamics for the variational parameters. Since this Hamiltonian is generically nonlinear and nonintegrable, the dynamics thus generated can be chaotic, in distinction to the exact quantum evolution. We then restrict our attention to a system of two biquadratically coupled quantum oscillators and study two variational schemes, the leading order large-$N$ (four canonical variables) and Hartree (six canonical variables) approximations. The chaos seen in the approximate dynamics is an artifact of the approximations: this is demonstrated by the fact that its onset occurs on the same characteristic time scale as the breakdown of the approximations when compared to numerical solutions of the time-dependent Schr\"odinger equation.

arXiv:quant-ph/9610013•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Feb 1, 1998•46 citations

Time-dependent variational principle for quantum dynamics in the Schrodinger-picture

Fred Cooper, So-Young Pi

Starting with the exact variational formulation of quantum mechanics (or quantum field theory) in the Schroedinger picture, and using Gaussian trail wave functions, an exact functional Schroedinger equation is obtained. The time-dependent Hartree-Fock approximation is obtained and extended to many oscillators and to field theory. The problem of quantum roll, in which an initial Gaussian wave packet starting at an unstable classical equilibrium point evolves in time, is treated using this formalism. Quantum tunnelling is also studied. (LEW)

•Jan 1, 1986•1 citation

Geometry and symmetry of quantum and classical-quantum dynamics

Esther Bonet-Luz

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.

•Feb 28, 2017•0 citations

A strategy for time dependent quantum mechanical calculations using a Gaussian wave packet representation of the wave function

Shin‐ichi Sawada, Robert Heather, Bret Jackson +1 more

We develop methodology for performing time dependent quantum mechanical calculations by representing the wave function as a sum of Gaussian wave packets (GWP) each characterized by a set of parameters such as width, position, momentum, and phase. The problem of computing the time evolution of the wave function is thus reduced to that of finding the time evolution of the parameters in the Gaussians. This parameter motion is determined by minimizing the error made by replacing the exact wave function in the time dependent Schrödinger equation with its Gaussian representation approximant. This leads to first order differential equations for the time dependence of the parameters, and those describing the packet position and the momentum of each packet have some resemblance with the classical equations of motion. The paper studies numerically the strategy needed to achieve the best GWP representation of time dependent processes. The issues discussed are the representation of the initial wave function, the numerical stability and the solution of the differential equations giving the evolution of the parameters, and the analysis of the final wave function. Extensive comparisons are made with an approximate method which assumes that the Gaussians are independent and their width is smaller than the length scale over which the potential changes. This approximation greatly simplifies the calculations and has the advantage of a greater resemblance to classical mechanics, thus being more intuitive. We find, however, that its range of applications is limited to problems involving localized degrees of freedom that participate in the dynamic process for a very short time. Finally we give particular attention to the notion that the GWP representation of the wave function reduces the dynamics of one quantum degree of freedom to that of a set of pseudoparticles (each represented by one packet) moving according to a pseudoclassical (i.e., classical-like) mechanics whose ‘‘phase space’’ is described by a position and momentum as well as a complex phase and width.

The Journal of Chemical Physics•Sep 15, 1985•205 citations
We investigate variational principles for the approximation of quantum dynamics that apply for approximation manifolds that do not have complex linear tangent spaces. The first one, dating back to McLachlan [Mol. Phys. 8, 39–44 (1964)], minimizes the residuum of the time-dependent Schrödinger equation, while the second one, originating from the lecture notes of Kramer and Saraceno [Geometry of the Time-Dependent Variational Principle in Quantum Mechanics, Lecture Notes in Physics Vol. 140 (Springer, Berlin, 1981)], imposes the stationarity of an action functional. We characterize both principles in terms of metric and symplectic orthogonality conditions, consider their conservation properties, and derive an elementary a posteriori error estimate. As an application, we revisit the time-dependent Hartree approximation and frozen Gaussian wave packets.
The Dirac-Frenkel-McLachlan variational principle is the basic tool for obtaining computationally accessible approximations in quantum molecular dynamics. It determines equations of motion for an approximate time-dependent wave function on an approximation manifold of reduced dimension. This paper gives a near-optimality result for variational approximations. It bounds the error in terms of the distance of the exact wave function to the approximation manifold and identifies the parameters that control the deviation of the variational approximation from the best approximation on the manifold.
The McLachlan "minimum-distance" principle for optimizing approximate solutions of the time-dependent Schrödinger equation is revisited, with a focus on the local-in-time error accompanying the variational solutions. Simple, exact expressions are provided for this error, which are then evaluated in illustrative cases, notably the widely used mean-field approach and the adiabatic quantum molecular dynamics. Based on these findings, we demonstrate the rigorous formulation of an adaptive scheme that resizes on the fly the underlying variational manifold and, hence, optimizes the overall computational cost of a quantum dynamical simulation. Such adaptive schemes are a crucial requirement for devising and applying direct quantum dynamical methods to molecular and condensed-phase problems.
Dynamical chaos has recently been shown to exist in the Gaussian approximation in quantum mechanics and in the self-consistent mean field approach to studying the dynamics of quantum fields. In this study, we first note that any variational approximation to the dynamics of a quantum system based on the Dirac action principle leads to a classical Hamiltonian dynamics for the variational parameters. Since this Hamiltonian is generically nonlinear and nonintegrable, the dynamics thus generated can be chaotic, in distinction to the exact quantum evolution. We then restrict our attention to a system of two biquadratically coupled quantum oscillators and study two variational schemes, the leading order large-NNN (four canonical variables) and Hartree (six canonical variables) approximations. The chaos seen in the approximate dynamics is an artifact of the approximations: this is demonstrated by the fact that its onset occurs on the same characteristic time scale as the breakdown of the approximations when compared to numerical solutions of the time-dependent Schr\"odinger equation.
Starting with the exact variational formulation of quantum mechanics (or quantum field theory) in the Schroedinger picture, and using Gaussian trail wave functions, an exact functional Schroedinger equation is obtained. The time-dependent Hartree-Fock approximation is obtained and extended to many oscillators and to field theory. The problem of quantum roll, in which an initial Gaussian wave packet starting at an unstable classical equilibrium point evolves in time, is treated using this formalism. Quantum tunnelling is also studied. (LEW)
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.
We develop methodology for performing time dependent quantum mechanical calculations by representing the wave function as a sum of Gaussian wave packets (GWP) each characterized by a set of parameters such as width, position, momentum, and phase. The problem of computing the time evolution of the wave function is thus reduced to that of finding the time evolution of the parameters in the Gaussians. This parameter motion is determined by minimizing the error made by replacing the exact wave function in the time dependent Schrödinger equation with its Gaussian representation approximant. This leads to first order differential equations for the time dependence of the parameters, and those describing the packet position and the momentum of each packet have some resemblance with the classical equations of motion. The paper studies numerically the strategy needed to achieve the best GWP representation of time dependent processes. The issues discussed are the representation of the initial wave function, the numerical stability and the solution of the differential equations giving the evolution of the parameters, and the analysis of the final wave function. Extensive comparisons are made with an approximate method which assumes that the Gaussians are independent and their width is smaller than the length scale over which the potential changes. This approximation greatly simplifies the calculations and has the advantage of a greater resemblance to classical mechanics, thus being more intuitive. We find, however, that its range of applications is limited to problems involving localized degrees of freedom that participate in the dynamic process for a very short time. Finally we give particular attention to the notion that the GWP representation of the wave function reduces the dynamics of one quantum degree of freedom to that of a set of pseudoparticles (each represented by one packet) moving according to a pseudoclassical (i.e., classical-like) mechanics whose ‘‘phase space’’ is described by a position and momentum as well as a complex phase and width.