Quantum Stochastic Molecular Dynamics Simulations of the Viscosity of Superfluid Helium

Decoherent quantum equations of motion are derived that yield the trajectory of an open quantum system. The viscosity of superfluid Lennard-Jones helium-4 is obtained with a quantum stochastic molecular dynamics algorithm. The momentum state occupancy entropy is counted with a continuous representation of boson number and averages are obtained with umbrella sampling. Instantaneous snapshots of the Bose-Einstein condensed system show multiple highly occupied momentum states. The viscosity is obtained from the Onsager-Green-Kubo relation with the time correlation function modified in the quantum case. On the saturation curve, at higher temperatures the viscosities of the classical and quantum liquids are equal. With decreasing temperature the viscosity of the classical liquid increases whereas that of the quantum liquid decreases. Below the $\lambda$-transition the viscosity lies significantly below the classical value, being small but positive due to the mixture of condensed and uncondensed bosons. The computed trajectories give a physical explanation of the molecular mechanism for superfluidity.