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

The Journal of Chemical Physics

Quantum trajectory dynamics in imaginary time with the momentum-dependent quantum potential

Sophya V. Garashchuk

Published Jan 7, 2010• 30 citations
Publisher

Abstract

The quantum trajectory dynamics is extended to the wave function evolution in imaginary time. For a nodeless wave function a simple exponential form leads to the classical-like equations of motion of trajectories, representing the wave function, in the presence of the momentum-dependent quantum potential in addition to the external potential. For a Gaussian wave function this quantum potential is a time-dependent constant, generating zero quantum force yet contributing to the total energy. For anharmonic potentials the momentum-dependent quantum potential is cheaply estimated from the global Least-squares Fit to the trajectory momenta in the Taylor basis. Wave functions with nodes are described in the mixed coordinate space/trajectory representation at little additional computational cost. The nodeless wave function, represented by the trajectory ensemble, decays to the ground state. The mixed representation wave functions, with lower energy contributions projected out at each time step, decay to the excited energy states. The approach, illustrated by computing energy levels for anharmonic oscillators and energy level splitting for the double-well potential, can be used for the Boltzmann operator evolution.

DOI
10.1063/1.3289728

Related papers

Multidimensional Quantum Trajectory Dynamics in Imaginary Time with Approximate Quantum Potential

Sophya V. Garashchuk, Tijo Vazhappilly

The quantum trajectory dynamics in imaginary time with the momentum-dependent quantum potential [ J. Chem. Phys. 2010, 132, 014112], which can be used for calculations of low-lying energy eigenstates or for Boltzmann evolution, is extended to many dimensions. Evolution of a nodeless initial wave function represented in exponential form is described in terms of a trajectory ensemble. In the Lagrangian frame of reference the trajectories evolve according to the classical equations of motion under the influence of the inverted classical potential and the quantum potential defined by the gradient of the trajectory momenta. The initial wave function is sampled by the trajectories randomly, which makes the approach scalable to high dimensions. Another practical feature is approximate computation of the quantum potential and force from the global least squares fit of the trajectory momenta in a small basis. The importance sampling and the Eulerian frame-of-reference formulation, compensating for divergent dynamics in bound potentials, are introduced to reduce the number of trajectories in multidimensional applications. The zero-point energies are computed for model and chemical (H 2 O and SO 2 ) systems.

The Journal of Physical Chemistry C•Sep 10, 2010•16 citations

Efficient quantum trajectory representation of wavefunctions evolving in imaginary time

Sophya V. Garashchuk, James William Mazzuca, Tijo Vazhappilly

The Boltzmann evolution of a wavefunction can be recast as imaginary-time dynamics of the quantum trajectory ensemble. The quantum effects arise from the momentum-dependent quantum potential--computed approximately to be practical in high-dimensional systems--influencing the trajectories in addition to the external classical potential [S. Garashchuk, J. Chem. Phys. 132, 014112 (2010)]. For a nodeless wavefunction represented as ψ(x, t) = exp(-S(x, t)/ħ) with the trajectory momenta defined by ∇S(x, t), analysis of the Lagrangian and Eulerian evolution shows that for bound potentials the former is more accurate while the latter is more practical because the Lagrangian quantum trajectories diverge with time. Introduction of stationary and time-dependent components into the wavefunction representation generates new Lagrangian-type dynamics where the trajectory spreading is controlled improving efficiency of the trajectory description. As an illustration, different types of dynamics are used to compute zero-point energy of a strongly anharmonic well and low-lying eigenstates of a high-dimensional coupled harmonic system.

The Journal of Chemical Physics•Jul 19, 2011•18 citations

Approximate Quantum Trajectory Dynamics in Imaginary and Real Time; Calculation of Reaction Rate Constants

In this chapter we formulate quantum trajectory dynamics in imaginary time which seamlessly connects to Bohmian dynamics in real time. Approximate determination of the quantum potential and mixed representation of wavefunctions with nodes make this trajectory approach stable and scalable to high dimensions, while including the T&F Cat # K11224, Chapter 23, Page 362, 15-9-2010 leading quantum effects into wavefunction evolution. Reaction probabilities, energy eigenstates and reaction rate constants can be obtained by propagating an ensemble of quantum trajectories in real and/or imaginary time.

CRC Press eBooks•Apr 19, 2016•0 citations

Quantum Trajectory Dynamics Based on Local Approximations to the Quantum Potential and Force

Sophya V. Garashchuk, Vitaly A. Rassolov

The quantum nature of nuclei often affects molecular structure and properties, which are associated with the motion of protons and other light nuclei at low temperature. However, incorporation of the nuclear quantum effects into theoretical studies of large molecular systems is an outstanding challenge in theoretical chemistry. In this paper, the de Broglie-Bohm formulation of the time-dependent Schrödinger equation is used as a framework for the development of approximate quantum corrections to the dynamics of a trajectory ensemble, representing a time-dependent wave function. Specifically, the trajectory-centered local Least Square Fit (L-LSF) to the trajectory attributes is introduced, formally analyzed and illustrated on model potentials, typically used to benchmark approximate and semiclassical dynamics methods. The L-LSF method improves the dynamics accuracy, compared to the globally defined approximations, developed previously, and is shown to converge to the exact (generally impractical for large systems) limit of quantum dynamics.

Journal of Chemical Theory and Computation•May 24, 2019•9 citations

Quantum Dynamics with the Quantum Trajectory-Guided Adaptable Gaussian Bases

Matthew Dutra, Sachith Wickramasinghe, Sophya V. Garashchuk

The computational cost of describing a general quantum system fully coupled by anharmonic interactions scales exponentially with the system size. Thus, an efficient basis representation of wave functions is essential, and when it comes to the large-amplitude motion of high-dimensional systems, the dynamic bases of Gaussian functions are often employed. The time dependence of such bases is determined from the variational principle or from classical dynamics; the former is challenging in implementation due to singular matrices, while the latter may not cover the configuration space relevant to quantum dynamics. Here we describe a method using Quantum Trajectory-guided Adaptable Gaussian (QTAG) bases "tuned"-including the basis position, phase, and width-to the wave function evolution, thanks to the continuity of the probability density in the course of the quantum trajectory dynamics. Thus, an efficient basis in configuration space is generated, bypassing the variational equations on the parameters of the Gaussians. We also propose a time propagator with basis transformation by projections which lends efficiency and stability to the QTAG dynamics, as demonstrated on standard tests and the ammonia inversion model.

Journal of Chemical Theory and Computation•Dec 4, 2019•22 citations

Semiclassical dynamics with quantum trajectories: Formulation and comparison with the semiclassical initial value representation propagator

Sophya V. Garashchuk, Vitaly A. Rassolov

We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.

The Journal of Chemical Physics•Jan 30, 2003•54 citations
The quantum trajectory dynamics is extended to the wave function evolution in imaginary time. For a nodeless wave function a simple exponential form leads to the classical-like equations of motion of trajectories, representing the wave function, in the presence of the momentum-dependent quantum potential in addition to the external potential. For a Gaussian wave function this quantum potential is a time-dependent constant, generating zero quantum force yet contributing to the total energy. For anharmonic potentials the momentum-dependent quantum potential is cheaply estimated from the global Least-squares Fit to the trajectory momenta in the Taylor basis. Wave functions with nodes are described in the mixed coordinate space/trajectory representation at little additional computational cost. The nodeless wave function, represented by the trajectory ensemble, decays to the ground state. The mixed representation wave functions, with lower energy contributions projected out at each time step, decay to the excited energy states. The approach, illustrated by computing energy levels for anharmonic oscillators and energy level splitting for the double-well potential, can be used for the Boltzmann operator evolution.
The quantum trajectory dynamics in imaginary time with the momentum-dependent quantum potential [ J. Chem. Phys. 2010, 132, 014112], which can be used for calculations of low-lying energy eigenstates or for Boltzmann evolution, is extended to many dimensions. Evolution of a nodeless initial wave function represented in exponential form is described in terms of a trajectory ensemble. In the Lagrangian frame of reference the trajectories evolve according to the classical equations of motion under the influence of the inverted classical potential and the quantum potential defined by the gradient of the trajectory momenta. The initial wave function is sampled by the trajectories randomly, which makes the approach scalable to high dimensions. Another practical feature is approximate computation of the quantum potential and force from the global least squares fit of the trajectory momenta in a small basis. The importance sampling and the Eulerian frame-of-reference formulation, compensating for divergent dynamics in bound potentials, are introduced to reduce the number of trajectories in multidimensional applications. The zero-point energies are computed for model and chemical (H 2 O and SO 2 ) systems.
The Boltzmann evolution of a wavefunction can be recast as imaginary-time dynamics of the quantum trajectory ensemble. The quantum effects arise from the momentum-dependent quantum potential--computed approximately to be practical in high-dimensional systems--influencing the trajectories in addition to the external classical potential [S. Garashchuk, J. Chem. Phys. 132, 014112 (2010)]. For a nodeless wavefunction represented as ψ(x, t) = exp(-S(x, t)/ħ) with the trajectory momenta defined by ∇S(x, t), analysis of the Lagrangian and Eulerian evolution shows that for bound potentials the former is more accurate while the latter is more practical because the Lagrangian quantum trajectories diverge with time. Introduction of stationary and time-dependent components into the wavefunction representation generates new Lagrangian-type dynamics where the trajectory spreading is controlled improving efficiency of the trajectory description. As an illustration, different types of dynamics are used to compute zero-point energy of a strongly anharmonic well and low-lying eigenstates of a high-dimensional coupled harmonic system.
In this chapter we formulate quantum trajectory dynamics in imaginary time which seamlessly connects to Bohmian dynamics in real time. Approximate determination of the quantum potential and mixed representation of wavefunctions with nodes make this trajectory approach stable and scalable to high dimensions, while including the T&F Cat # K11224, Chapter 23, Page 362, 15-9-2010 leading quantum effects into wavefunction evolution. Reaction probabilities, energy eigenstates and reaction rate constants can be obtained by propagating an ensemble of quantum trajectories in real and/or imaginary time.
The quantum nature of nuclei often affects molecular structure and properties, which are associated with the motion of protons and other light nuclei at low temperature. However, incorporation of the nuclear quantum effects into theoretical studies of large molecular systems is an outstanding challenge in theoretical chemistry. In this paper, the de Broglie-Bohm formulation of the time-dependent Schrödinger equation is used as a framework for the development of approximate quantum corrections to the dynamics of a trajectory ensemble, representing a time-dependent wave function. Specifically, the trajectory-centered local Least Square Fit (L-LSF) to the trajectory attributes is introduced, formally analyzed and illustrated on model potentials, typically used to benchmark approximate and semiclassical dynamics methods. The L-LSF method improves the dynamics accuracy, compared to the globally defined approximations, developed previously, and is shown to converge to the exact (generally impractical for large systems) limit of quantum dynamics.
The computational cost of describing a general quantum system fully coupled by anharmonic interactions scales exponentially with the system size. Thus, an efficient basis representation of wave functions is essential, and when it comes to the large-amplitude motion of high-dimensional systems, the dynamic bases of Gaussian functions are often employed. The time dependence of such bases is determined from the variational principle or from classical dynamics; the former is challenging in implementation due to singular matrices, while the latter may not cover the configuration space relevant to quantum dynamics. Here we describe a method using Quantum Trajectory-guided Adaptable Gaussian (QTAG) bases "tuned"-including the basis position, phase, and width-to the wave function evolution, thanks to the continuity of the probability density in the course of the quantum trajectory dynamics. Thus, an efficient basis in configuration space is generated, bypassing the variational equations on the parameters of the Gaussians. We also propose a time propagator with basis transformation by projections which lends efficiency and stability to the QTAG dynamics, as demonstrated on standard tests and the ammonia inversion model.
We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.