This paper introduces a new computational method for molecular dynamics. The method combines the Backward-Euler scheme for the solution of stiff differential equations with a Langevin-equation approach to the establishment of thermal equilibrium. The method allows the user to choose a cutoof frequency _c. Vibrational modes with frequencies below _c will be fully excited (receive a mean energy of kT per mode), while modes with frequencies greater than _c will be effectively frozen by the method. By setting _c=kT/, one can obtain reasonable agreement with the quantum-mechanical energy distribution among the various modes, despite the classical character of the computation.

Click to go back to the publication list

Webpage design by Igor Zilberman