AMBER Archive (2004)

Subject: Re: AMBER: Parameters for deuterium

From: David A. Case (
Date: Tue Jan 06 2004 - 10:32:14 CST

On Tue, Jan 06, 2004, Charlie Laughton wrote:
> I'd appreciate any thoughts/experiences on simulations incorporating
> deuterium (C-D or O-D bonds). Is it enough just to take corresponding H
> parameters and up the mass?

The force field is unchanged, and the masses do change. But you need to
remember that in classical mechanics, the distribution of configurations
is independent of the atomic masses, and so there are no "isotope effects" in
the usual sense. This means that the properties of MD simulations are not
affected by changing H to D. You would need to incorporate quantum effects
(especially zero-point motion) to see the real differences upon isotope

The one place where a simple classical model would work is for vibrational
frequencies: here if you change the atomic masses, you would see the
appropriate change in frequencies.

..hope this helps...dac


================================================================== David A. Case | e-mail: Dept. of Molecular Biology, TPC15 | fax: +1-858-784-8896 The Scripps Research Institute | phone: +1-858-784-9768 10550 N. Torrey Pines Rd. | home page: La Jolla CA 92037 USA | ================================================================== ----------------------------------------------------------------------- The AMBER Mail Reflector To post, send mail to To unsubscribe, send "unsubscribe amber" to