Isothermal quantum dynamics: Nosé-Hoover method for coherent states

D. Mentrup, J. Schnack

Detlef Mentrup, Jürgen Schnack

Detlef Mentrup, Jürgen Schnack

Detlef Mentrup

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.

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.

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.