Nonequilibrium molecular dynamics via Gauss's principle of least constraint

Gauss's principle of least constraint is used to develop nonequilibrium molecular-dynamics algo­ rithms for systems SUbject to constraints. The treatment not only includes "nonholonomic" constraints-those involving velocities-but it also provides a basis for simulating nonequilibrium steady states. We describe two applications of this new use of Gauss's principle. The first of these examples, the isothermal molecular dynamics of a three-particle chain, can be treated analytically. The second, the steady-state diffusiort~of a Lennard-Jones liquid, near its triple point, is studied nu­ merically. The measured diffusion coefficient agrees with inaependent estimates from eqUilibrium fluctuation theory and from Hamiltonian external fields. I. INTRODUCTION