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

Isothermal quantum dynamics: Investigations for the harmonic oscillator

Detlef Mentrup

Published May 26, 2003• 0 citations
OpenAlex

Abstract

Thermostated time evolutions are on a firm ground and widely used in classical molecular dynamics (MD) simulations. Hamilton´s equations of motion are supplemented by time-dependent pseudofriction terms that convert the microcanonical isoenergetic time evolution into a canonical isothermal time evolution, thus permitting the calculation of canonical ensemble averages by time averaging. However, similar methods for quantum MD schemes are still lacking. Given the rich dynamical behavior of ultracold trapped quantum gases depending on the value of the s-wave scattering length, it is timely to investigate how classical thermostating methods can be combined with powerful approximate quantum dynamics schemes to deal with interacting quantum systems at finite temperature. In this work, the popular method of Nose and Hoover to create canonically distributed positions and momenta in classical MD simulations is generalized to a genuine quantum system of infinite dimensionality. We show that for the quantum harmonic oscillator, the equations of motion in terms of coherent states may be modified in a Nose-Hoover manner to mimic the coupling of the system to a thermal bath and create a quantum canonical ensemble. The method is developed initially for a single particle and then generalized to the case of an arbitrary number of identical quantum particles, involving entangled distribution functions. The resulting isothermal equations of motion for bosons and fermions contain additional terms leading to Bose-attraction and Pauli-blocking, respectively. Questions of ergodicity are discussed for different coupling schemes. In the many-particle case, the superiority of the Nose-Hoover technique to a Langevin approach is demonstrated. In addition, the work contains an investigation of the Grilli-Tosatti thermostating method applied to the harmonic oscillator, and calculations for quantum wavefunctions moving with a time-invariant shape in a harmonic potential.

Related papers

ISOTHERMAL QUANTUM DYNAMICS: NOSÉ–HOOVER DYNAMICS FOR COHERENT STATES

Detlef Mentrup, Jürgen Schnack

The method of Nosé and Hoover1,2 to create canonically distributed positions and momenta in classical molecular dynamics simulations is frequently used. Hamilton's equations of motion are supplemented by time-dependent pseudofriction terms that convert the microcanonical isoenergetic time evolution into a canonical isothermal time evolution, thus permitting the calculation of canonical ensemble averages by time averaging. We show that for one quantum particle in an external harmonic oscillator, the equations of motion in terms of coherent states can easily be modified in an analogous manner to mimic the coupling of the system to a thermal bath and create a quantum canonical ensemble.3 The method is generalised to a system of two identical quantum particles. In the resulting equations of motion, one obtains an additional attractive term for bosons and a repulsive term for fermions in the dynamics of the pseudofriction coefficients, leading to a correctly sampled thermal weight.

International Journal of Modern Physics B•Nov 10, 2003•2 citations

Nosé-Hoover dynamics for coherent states

Mentrup, D., Schnack, J.

The popular method of Nosé and Hoover to create canonically distributed positions and momenta in classical molecular dynamics simulations is generalized to a genuine quantum system of infinite dimensionality. We show that for the quantum harmonic oscillator, the equations of motion in terms of coherent states can easily be modified in an analogous manner to mimic the coupling of the system to a thermal bath and create a quantum canonical ensemble. Possible applications to more complex systems, especially interacting Fermion systems, are proposed.

•Jan 1, 2001•9 citations

Nosé-Hoover sampling of quantum entangled distribution functions

Mentrup, D., Schnack, J.

While thermostated time evolutions stand on firm grounds and are widely used in classical molecular dynamics (MD) simulations (J. Phys. Chem. B 104 (2000) 159), similar methods for quantum MD schemes are still lacking. In the special case of a quantum particle in a harmonic potential, it has been shown that the framework of coherent states permits to set up equations of motion for an isothermal quantum dynamics (Physica A 297 (2001) 337). In the present article, these results are generalized to indistinguishable quantum particles. We investigate the consequences of the (anti-)symmetry of the many-particle wavefunction which leads to quantum entangled distribution functions. The resulting isothermal equations of motion for bosons and fermions contain new terms which cause Bose-attraction and Pauli-blocking. Questions of ergodicity are discussed for different coupling schemes.

•Jan 1, 2003•2 citations

Statistics and Nosé formalism for Ehrenfest dynamics

J L Alonso, A Castro, Jesús Clemente-Gallardo +3 more

Quantum dynamics (i.e. the Schrödinger equation) and classical dynamics (i.e. Hamilton equations) can both be formulated in equal geometric terms: a Poisson bracket defined on a manifold. In this paper, we first show that the hybrid quantum-classical dynamics prescribed by the Ehrenfest equations can also be formulated within this general framework, what has been used in the literature to construct propagation schemes for Ehrenfest dynamics. Then, the existence of a well-defined Poisson bracket allows us to arrive to a Liouville equation for a statistical ensemble of Ehrenfest systems. The study of a generic toy model shows that the evolution produced by Ehrenfest dynamics is ergodic and therefore the only constants of motion are functions of the Hamiltonian. The emergence of the canonical ensemble characterized by the Boltzmann distribution follows after an appropriate application of the principle of equal a priori probabilities to this case. Once we know the canonical distribution of an Ehrenfest system, it is straightforward to extend the formalism of Nosé (invented to do constant temperature molecular dynamics by a non-stochastic method) to our Ehrenfest formalism. This work also provides the basis for extending stochastic methods to Ehrenfest dynamics.

arXiv:1104.2154•Journal of Physics A Mathematical and Theoretical•Sep 7, 2011•29 citations

Canonical averaging in the second order quantized Hamilton dynamics by extension of the coherent state thermodynamics of the harmonic oscillator

Eric M. Heatwole, Oleg V. Prezhdo

A conceptually simple approximation to quantum mechanics, quantized Hamilton dynamics (QHD) includes zero-point energy, tunneling, dephasing, and other important quantum effects in a classical-like description. The hierarchy of coupled differential equations describing the time evolution of observables in QHD can be mapped in the second order onto a classical system with double the dimensionality of the original system. While QHD excels at dynamics with a single initial condition, the correct method for generating thermal initial conditions in QHD remains an open question. Using the coherent state representation of thermodynamics of the harmonic oscillator (HO) [Schnack, Europhys. Lett. 45, 647 (1999)], we develop canonical averaging for the second order QHD [Prezhdo, J. Chem. Phys. 117, 2995 (2002)]. The methodology is exact for the free particle and HO, and shows good agreement with quantum results for a variety of quartic potentials.

The Journal of Chemical Physics•May 28, 2007•21 citations

Quantum dynamics in classical constant-temperature baths.

Nkosinathi Dlamini

In this dissertation the formulation of various integration algorithms is studied, with a view to simulate quantum-classical systems in contact with a thermal bath. In particular focus is given to the constant temperature dynamics of the Nos e-Hoover, Nos e-Hoover Chain and Nos e-Hoover Power thermostat schemes. Through the use of the time symmetric Trotter factorisation of the Liouville operator, algorithms are derived that are both time-reversible and measure-preserving. The efficiency of these algorithms is tested via the constant temperature simulation of a low-dimensional harmonic system. In addition The Nos e-Hoover Power thermostat was then extended to the quantum-classical case. The damping of a tunnelling spin coupled to a thermalised harmonic mode was simulated and the results are presented.

•May 23, 2013•0 citations
Thermostated time evolutions are on a firm ground and widely used in classical molecular dynamics (MD) simulations. Hamilton´s equations of motion are supplemented by time-dependent pseudofriction terms that convert the microcanonical isoenergetic time evolution into a canonical isothermal time evolution, thus permitting the calculation of canonical ensemble averages by time averaging. However, similar methods for quantum MD schemes are still lacking. Given the rich dynamical behavior of ultracold trapped quantum gases depending on the value of the s-wave scattering length, it is timely to investigate how classical thermostating methods can be combined with powerful approximate quantum dynamics schemes to deal with interacting quantum systems at finite temperature. In this work, the popular method of Nose and Hoover to create canonically distributed positions and momenta in classical MD simulations is generalized to a genuine quantum system of infinite dimensionality. We show that for the quantum harmonic oscillator, the equations of motion in terms of coherent states may be modified in a Nose-Hoover manner to mimic the coupling of the system to a thermal bath and create a quantum canonical ensemble. The method is developed initially for a single particle and then generalized to the case of an arbitrary number of identical quantum particles, involving entangled distribution functions. The resulting isothermal equations of motion for bosons and fermions contain additional terms leading to Bose-attraction and Pauli-blocking, respectively. Questions of ergodicity are discussed for different coupling schemes. In the many-particle case, the superiority of the Nose-Hoover technique to a Langevin approach is demonstrated. In addition, the work contains an investigation of the Grilli-Tosatti thermostating method applied to the harmonic oscillator, and calculations for quantum wavefunctions moving with a time-invariant shape in a harmonic potential.
The method of Nosé and Hoover1,2 to create canonically distributed positions and momenta in classical molecular dynamics simulations is frequently used. Hamilton's equations of motion are supplemented by time-dependent pseudofriction terms that convert the microcanonical isoenergetic time evolution into a canonical isothermal time evolution, thus permitting the calculation of canonical ensemble averages by time averaging. We show that for one quantum particle in an external harmonic oscillator, the equations of motion in terms of coherent states can easily be modified in an analogous manner to mimic the coupling of the system to a thermal bath and create a quantum canonical ensemble.3 The method is generalised to a system of two identical quantum particles. In the resulting equations of motion, one obtains an additional attractive term for bosons and a repulsive term for fermions in the dynamics of the pseudofriction coefficients, leading to a correctly sampled thermal weight.
The popular method of Nosé and Hoover to create canonically distributed positions and momenta in classical molecular dynamics simulations is generalized to a genuine quantum system of infinite dimensionality. We show that for the quantum harmonic oscillator, the equations of motion in terms of coherent states can easily be modified in an analogous manner to mimic the coupling of the system to a thermal bath and create a quantum canonical ensemble. Possible applications to more complex systems, especially interacting Fermion systems, are proposed.
While thermostated time evolutions stand on firm grounds and are widely used in classical molecular dynamics (MD) simulations (J. Phys. Chem. B 104 (2000) 159), similar methods for quantum MD schemes are still lacking. In the special case of a quantum particle in a harmonic potential, it has been shown that the framework of coherent states permits to set up equations of motion for an isothermal quantum dynamics (Physica A 297 (2001) 337). In the present article, these results are generalized to indistinguishable quantum particles. We investigate the consequences of the (anti-)symmetry of the many-particle wavefunction which leads to quantum entangled distribution functions. The resulting isothermal equations of motion for bosons and fermions contain new terms which cause Bose-attraction and Pauli-blocking. Questions of ergodicity are discussed for different coupling schemes.
Quantum dynamics (i.e. the Schrödinger equation) and classical dynamics (i.e. Hamilton equations) can both be formulated in equal geometric terms: a Poisson bracket defined on a manifold. In this paper, we first show that the hybrid quantum-classical dynamics prescribed by the Ehrenfest equations can also be formulated within this general framework, what has been used in the literature to construct propagation schemes for Ehrenfest dynamics. Then, the existence of a well-defined Poisson bracket allows us to arrive to a Liouville equation for a statistical ensemble of Ehrenfest systems. The study of a generic toy model shows that the evolution produced by Ehrenfest dynamics is ergodic and therefore the only constants of motion are functions of the Hamiltonian. The emergence of the canonical ensemble characterized by the Boltzmann distribution follows after an appropriate application of the principle of equal a priori probabilities to this case. Once we know the canonical distribution of an Ehrenfest system, it is straightforward to extend the formalism of Nosé (invented to do constant temperature molecular dynamics by a non-stochastic method) to our Ehrenfest formalism. This work also provides the basis for extending stochastic methods to Ehrenfest dynamics.
A conceptually simple approximation to quantum mechanics, quantized Hamilton dynamics (QHD) includes zero-point energy, tunneling, dephasing, and other important quantum effects in a classical-like description. The hierarchy of coupled differential equations describing the time evolution of observables in QHD can be mapped in the second order onto a classical system with double the dimensionality of the original system. While QHD excels at dynamics with a single initial condition, the correct method for generating thermal initial conditions in QHD remains an open question. Using the coherent state representation of thermodynamics of the harmonic oscillator (HO) [Schnack, Europhys. Lett. 45, 647 (1999)], we develop canonical averaging for the second order QHD [Prezhdo, J. Chem. Phys. 117, 2995 (2002)]. The methodology is exact for the free particle and HO, and shows good agreement with quantum results for a variety of quartic potentials.
In this dissertation the formulation of various integration algorithms is studied, with a view to simulate quantum-classical systems in contact with a thermal bath. In particular focus is given to the constant temperature dynamics of the Nos e-Hoover, Nos e-Hoover Chain and Nos e-Hoover Power thermostat schemes. Through the use of the time symmetric Trotter factorisation of the Liouville operator, algorithms are derived that are both time-reversible and measure-preserving. The efficiency of these algorithms is tested via the constant temperature simulation of a low-dimensional harmonic system. In addition The Nos e-Hoover Power thermostat was then extended to the quantum-classical case. The damping of a tunnelling spin coupled to a thermalised harmonic mode was simulated and the results are presented.