**From:** Boyang Wang (*pkuwangboyang_at_yahoo.com.cn*)

**Date:** Wed Mar 01 2006 - 12:18:30 CST

**Next message:**Philip Blood: "Re: Re: NAMD with Myrinet MX on MacOS X"**Previous message:**jz7_at_duke.edu: "question about FEP calculation"**Next in thread:**Boyang Wang: "Re: About Langevin Dynamics."**Maybe reply:**Boyang Wang: "Re: About Langevin Dynamics."**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

Hi Jerome and all,

I found your previous email. That is really helpful. But my question is:

The Gamma in Langevin equation should have a unit of (s)^-1, as "mass*gamma*v" has the unit of force.

But in NAMD configuration files, the Langevin damping coefficient has a unit of (s) or (ps), which puzzles me. I thought this Langevin damping coefficient was the reciprocal of the "gamma" in the Langevin equation.

Thanks for all your time.

Boyang.

Tue Mar 02 2004 - 15:39:24 CST

J¨¦rôme wrote:

"

When doing Langevin dynamics, NAMD adds two forces to those deriving from the

force field :

* a frictional force, equal to - mass * gamma * velocity (gamma is the damping

coefficient)

* a stochastic force, that is normally distributed with zero average and a rms

value sigma.

As I said previously, the frictional term constantly drains the system's

kinetic energy, while the stochastic force gives back a certain amount of

energy.

To use this as a thermostat, the force rms, sigma, has to be adjusted using

the fluctuation-dissipation theorem. In this case, the FDT yields :

sigma = sqrt( 2 * kT * gamma * mass / delta_t), where delta_t is the timestep.

So NAMD uses the value of gamma provided by the user and computes accordingly

the value of sigma that will lead to the desired temperature. You can see

that if gamma is very small, sigma will be small too, so that the overall

Langevin force will be negligible compared to "physical" ones from the force

field.

All this is likely to be explained with much detail in nonequilibrium

statistical mechanics books.

Cheers,

Jerome

"

---------------------------------

ÑÅ»¢1GÃâ·ÑÓÊÏä°Ù·Ö°Ù·ÀÀ¬»øÐÅ

ÑÅ»¢ÖúÊÖ-ËÑË÷¡¢É±¶¾¡¢·ÀÉ§ÈÅ

**Next message:**Philip Blood: "Re: Re: NAMD with Myrinet MX on MacOS X"**Previous message:**jz7_at_duke.edu: "question about FEP calculation"**Next in thread:**Boyang Wang: "Re: About Langevin Dynamics."**Maybe reply:**Boyang Wang: "Re: About Langevin Dynamics."**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

*
This archive was generated by hypermail 2.1.6
: Wed Feb 29 2012 - 15:41:40 CST
*