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

Journal of Chemical Theory and Computation

Estimation of the Ground State Energy of an Atomic Solid by Employing Quantum Trajectory Dynamics with Friction

Bing Gu, Robert J. Hinde, Vitaly A. Rassolov, Sophya V. Garashchuk

Published May 20, 2015• 20 citations
Publisher

Abstract

Evolution with energy dissipation can be used to obtain the ground state of a quantum-mechanical system. This dissipation is introduced in the quantum trajectory framework by adding an empirical friction force to the equations of motion for the trajectories, which, as an ensemble, represent a wave function. The quantum effects in dynamics are incorporated via the quantum force derived from the properties of this ensemble. For scalability to large systems, the quantum force is computed approximately yet with sufficient accuracy to describe the strongly anharmonic ground state of solid (4)He represented by a simulation cell of 180 atoms.

DOI
10.1021/ct501176m

Related papers

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

The Schrödinger equation with friction from the quantum trajectory perspective

Sophya V. Garashchuk, Vaibhav Dixit, Bing Gu +1 more

Similarity of equations of motion for the classical and quantum trajectories is used to introduce a friction term dependent on the wavefunction phase into the time-dependent Schrödinger equation. The term describes irreversible energy loss by the quantum system. The force of friction is proportional to the velocity of a quantum trajectory. The resulting Schrödinger equation is nonlinear, conserves wavefunction normalization, and evolves an arbitrary wavefunction into the ground state of the system (of appropriate symmetry if applicable). Decrease in energy is proportional to the average kinetic energy of the quantum trajectory ensemble. Dynamics in the high friction regime is suitable for simple models of reactions proceeding with energy transfer from the system to the environment. Examples of dynamics are given for single and symmetric and asymmetric double well potentials.

The Journal of Chemical Physics•Feb 1, 2013•35 citations

Microscopic Quantum Friction

Pedro H. Pereira, F. Impens, C. Farina +2 more

We report on a microscopic theory of quantum friction. Our approach investigates the interplay between the dispersive response and the relative center-of-mass motion of two ground-state atoms. This coupling yields a quantum force, which can be expressed as a power series in the velocity. The significance of each contribution depends on its order parity: while even-order terms are reversible, odd-order terms are irreversible and only survive in the presence of an internal dissipation mechanism. In addition, we obtain general, model-independent properties for the work performed by these contributions for arbitrary scattering trajectories. These results enable an unambiguous identification of odd-parity terms with microscopic quantum friction. At room temperature, the dominant microscopic quantum friction is of first order in the velocity and presents a strong quantum character. Our microscopic theory reveals that several properties of quantum friction obtained in specific settings -- such as the cubic dependence on velocity at zero temperature -- are indeed universal features already present at the atomic scale.

quant-pharXiv:2601.13265•Jan 20, 2026

Quantum Langevin Dynamics

Mohammad Attrash, Roi Baer

Previous years researchers began to simulate open quantum system, taking into account the interaction between system and the environment. One approach to deal with this problem is to use the density matrix within the Liouville-von-Neumann formalism or the Markovian variant the Lindblad equations. Another way is to use a stochastic approach where a random force is added to the system. The benefit of the stochastic approach is to solve the dynamics of the system with less time and memory than the density matrix approaches. In this project we want to develop a stochastic approach that can deal with the stochastic wave functions approach. We did this on a 2-level system and found that it works well when comparing to a density matrix approach. Next, we tested a quantum particle connect to a bath of harmonic oscillators using the stochastic approach. We found that a friction term is necessary and applied it. Like in the classical Langevin equations the friction constant and the random force fluctuations are related by the fluctuation-dissipation constant. We showed that with friction the dynamics decays to an ensemble with energy of $E_{gs}+k_BT$. However, we also found here are problems. The system seems to absorb energy indefinitely if the temperature is higher than the zero point energy or if the system is a Morse oscillator. Thus more research is required to make this method work.

arXiv:2509.09743•Sep 11, 2025•0 citations

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

Time Evolution of Wigner Function in Dissipative Quantum Systems Using Entangled Trajectory Molecular Dynamics

Lifei Wang Lifei Wang, Yan Wang Yan Wang, Qin Zhang Qin Zhang +2 more

The dissipative quantum systems are treated using Klein-Kramers equation, combined with the Gaussian kernel trajectory ensemble, for time evolution of Wigner function $ρ_ω(q, p, t)$ in phase space. The entangled trajectory molecular dynamics approach is used to obtain trajectory solutions for the Klein-Kramers equation with three models: free particle, damped harmonic oscillator and metastable potential. It is found that the performance of semiclassical Wigner propagation is effectively for the relaxation of damped harmonic oscillator and dissipative decay of a metastable state. In addition, the energy of trajectory ensemble decays faster with smaller friction value and changes slightly with variable temperature parameters.

The Journal of Atomic and Molecular Sciences•Jul 1, 2016•0 citations
Evolution with energy dissipation can be used to obtain the ground state of a quantum-mechanical system. This dissipation is introduced in the quantum trajectory framework by adding an empirical friction force to the equations of motion for the trajectories, which, as an ensemble, represent a wave function. The quantum effects in dynamics are incorporated via the quantum force derived from the properties of this ensemble. For scalability to large systems, the quantum force is computed approximately yet with sufficient accuracy to describe the strongly anharmonic ground state of solid (4)He represented by a simulation cell of 180 atoms.
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.
Similarity of equations of motion for the classical and quantum trajectories is used to introduce a friction term dependent on the wavefunction phase into the time-dependent Schrödinger equation. The term describes irreversible energy loss by the quantum system. The force of friction is proportional to the velocity of a quantum trajectory. The resulting Schrödinger equation is nonlinear, conserves wavefunction normalization, and evolves an arbitrary wavefunction into the ground state of the system (of appropriate symmetry if applicable). Decrease in energy is proportional to the average kinetic energy of the quantum trajectory ensemble. Dynamics in the high friction regime is suitable for simple models of reactions proceeding with energy transfer from the system to the environment. Examples of dynamics are given for single and symmetric and asymmetric double well potentials.
We report on a microscopic theory of quantum friction. Our approach investigates the interplay between the dispersive response and the relative center-of-mass motion of two ground-state atoms. This coupling yields a quantum force, which can be expressed as a power series in the velocity. The significance of each contribution depends on its order parity: while even-order terms are reversible, odd-order terms are irreversible and only survive in the presence of an internal dissipation mechanism. In addition, we obtain general, model-independent properties for the work performed by these contributions for arbitrary scattering trajectories. These results enable an unambiguous identification of odd-parity terms with microscopic quantum friction. At room temperature, the dominant microscopic quantum friction is of first order in the velocity and presents a strong quantum character. Our microscopic theory reveals that several properties of quantum friction obtained in specific settings -- such as the cubic dependence on velocity at zero temperature -- are indeed universal features already present at the atomic scale.
Previous years researchers began to simulate open quantum system, taking into account the interaction between system and the environment. One approach to deal with this problem is to use the density matrix within the Liouville-von-Neumann formalism or the Markovian variant the Lindblad equations. Another way is to use a stochastic approach where a random force is added to the system. The benefit of the stochastic approach is to solve the dynamics of the system with less time and memory than the density matrix approaches. In this project we want to develop a stochastic approach that can deal with the stochastic wave functions approach. We did this on a 2-level system and found that it works well when comparing to a density matrix approach. Next, we tested a quantum particle connect to a bath of harmonic oscillators using the stochastic approach. We found that a friction term is necessary and applied it. Like in the classical Langevin equations the friction constant and the random force fluctuations are related by the fluctuation-dissipation constant. We showed that with friction the dynamics decays to an ensemble with energy of Egs+kBTE_{gs}+k_BTEgs​+kB​T. However, we also found here are problems. The system seems to absorb energy indefinitely if the temperature is higher than the zero point energy or if the system is a Morse oscillator. Thus more research is required to make this method work.
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 dissipative quantum systems are treated using Klein-Kramers equation, combined with the Gaussian kernel trajectory ensemble, for time evolution of Wigner function ρω(q,p,t)ρ_ω(q, p, t)ρω​(q,p,t) in phase space. The entangled trajectory molecular dynamics approach is used to obtain trajectory solutions for the Klein-Kramers equation with three models: free particle, damped harmonic oscillator and metastable potential. It is found that the performance of semiclassical Wigner propagation is effectively for the relaxation of damped harmonic oscillator and dissipative decay of a metastable state. In addition, the energy of trajectory ensemble decays faster with smaller friction value and changes slightly with variable temperature parameters.