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

Modified quantum trajectory dynamics using a mixed wave function representation

Sophya V. Garashchuk, Vitaly A. Rassolov

Published Nov 2, 2004• 30 citations
Publisher

Abstract

Dynamics of quantum trajectories provides an efficient framework for description of various quantum effects in large systems, but it is unstable near the wave function density nodes where the quantum potential becomes singular. A mixed coordinate space/polar representation of the wave function is used to circumvent this problem. The resulting modified trajectory dynamics associated with the polar representation is nonsingular and smooth. The interference structure and the nodes of the wave function density are described, in principle, exactly in the coordinate representation. The approximate version of this approach is consistent with the semiclassical linearized quantum force method [S. Garashchuk and V. A. Rassolov, J. Chem. Phys. 120, 1181 (2004)]. This approach is exact for general wave functions with the density nodes in a locally quadratic potential.

DOI
10.1063/1.1804177

Related papers

Semiclassical nonadiabatic dynamics with quantum trajectories

Vitaly A. Rassolov, Sophya V. Garashchuk

Dynamics based on quantum trajectories with approximate quantum potential is generalized to nonadiabatic systems and its semiclassical properties are discussed. The formulation uses the mixed polar-coordinate space representation of a wave function. The polar part describes the overall time evolution of the wave-function components semiclassically using the single-surface approximate quantum potential. The coordinate part represents a complex``population'' amplitude, which in case of localized coupling can be solved for quantum mechanically in an efficient manner. In the high-energy regime this is accomplished by using a small basis determined by the coupling between surfaces. An illustration is given for a typical curve-crossing problem. The energy-resolved probabilities obtained from the time evolution of two wave packets for a wide range of energies are in excellent agreement with exact results for energies above the threshold of the diabatic reaction, including the case of total nonadiabatic transition.

Physical Review A•Mar 25, 2005•46 citations

Computation of correlation functions and wave function projections in the context of quantum trajectory dynamics

Sophya V. Garashchuk

The de Broglie-Bohm formulation of the Schrodinger equation implies conservation of the wave function probability density associated with each quantum trajectory in closed systems. This conservation property greatly simplifies numerical implementations of the quantum trajectory dynamics and increases its accuracy. The reconstruction of a wave function, however, becomes expensive or inaccurate as it requires fitting or interpolation procedures. In this paper we present a method of computing wave packet correlation functions and wave function projections, which typically contain all the desired information about dynamics, without the full knowledge of the wave function by making quadratic expansions of the wave function phase and amplitude near each trajectory similar to expansions used in semiclassical methods. Computation of the quantities of interest in this procedure is linear with respect to the number of trajectories. The introduced approximations are consistent with approximate quantum potential dynamics method. The projection technique is applied to model chemical systems and to the H+H(2) exchange reaction in three dimensions.

The Journal of Chemical Physics•Apr 17, 2007•5 citations

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

Sophya V. Garashchuk

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 Journal of Chemical Physics•Jan 7, 2010•30 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

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

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
Dynamics of quantum trajectories provides an efficient framework for description of various quantum effects in large systems, but it is unstable near the wave function density nodes where the quantum potential becomes singular. A mixed coordinate space/polar representation of the wave function is used to circumvent this problem. The resulting modified trajectory dynamics associated with the polar representation is nonsingular and smooth. The interference structure and the nodes of the wave function density are described, in principle, exactly in the coordinate representation. The approximate version of this approach is consistent with the semiclassical linearized quantum force method [S. Garashchuk and V. A. Rassolov, J. Chem. Phys. 120, 1181 (2004)]. This approach is exact for general wave functions with the density nodes in a locally quadratic potential.
Dynamics based on quantum trajectories with approximate quantum potential is generalized to nonadiabatic systems and its semiclassical properties are discussed. The formulation uses the mixed polar-coordinate space representation of a wave function. The polar part describes the overall time evolution of the wave-function components semiclassically using the single-surface approximate quantum potential. The coordinate part represents a complex``population'' amplitude, which in case of localized coupling can be solved for quantum mechanically in an efficient manner. In the high-energy regime this is accomplished by using a small basis determined by the coupling between surfaces. An illustration is given for a typical curve-crossing problem. The energy-resolved probabilities obtained from the time evolution of two wave packets for a wide range of energies are in excellent agreement with exact results for energies above the threshold of the diabatic reaction, including the case of total nonadiabatic transition.
The de Broglie-Bohm formulation of the Schrodinger equation implies conservation of the wave function probability density associated with each quantum trajectory in closed systems. This conservation property greatly simplifies numerical implementations of the quantum trajectory dynamics and increases its accuracy. The reconstruction of a wave function, however, becomes expensive or inaccurate as it requires fitting or interpolation procedures. In this paper we present a method of computing wave packet correlation functions and wave function projections, which typically contain all the desired information about dynamics, without the full knowledge of the wave function by making quadratic expansions of the wave function phase and amplitude near each trajectory similar to expansions used in semiclassical methods. Computation of the quantities of interest in this procedure is linear with respect to the number of trajectories. The introduced approximations are consistent with approximate quantum potential dynamics method. The projection technique is applied to model chemical systems and to the H+H(2) exchange reaction in three dimensions.
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 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.
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 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.