Error in NAMD Controller Class Reference- Adaptive Tempering

From: Sachin Natesh (snatesh_at_uchicago.edu)
Date: Thu Sep 10 2015 - 19:13:41 CDT

Hello all,

In the NAMD controller class reference, which can be found here
<http://www.ks.uiuc.edu/Research/namd/doxygen/classController.html>, I
notice a possible error in the computation of the adaptive tempering
temperature update (lines 1875-1878; see below)

01875 //dT is new temperature
01876 BigReal
<http://www.ks.uiuc.edu/Research/namd/doxygen/common_8h.html#a40> dT =
((potentialEnergy-potEnergyAverage)/BOLTZMANN
<http://www.ks.uiuc.edu/Research/namd/doxygen/common_8h.html#a1>+adaptTempT)*adaptTempDt;
01877 dT += random
<http://www.ks.uiuc.edu/Research/namd/doxygen/classController.html#p83>->gaussian
<http://www.ks.uiuc.edu/Research/namd/doxygen/classRandom.html#a6>()**sqrt(2.*adaptTempDt)**adaptTempT;
01878 dT += adaptTempT;

On line 1877, the stochastic portion of the temperature update is added to
what has been computed for the new temperature thus far, though the
integration time-step, adaptTempDt, is placed within the sqrt function.
This seems particularly erroneous to me as *equation 1* of the paper upon
which the NAMD implementation is based (here
<http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2905458/>) calls for the
integration time-step to fall outside of the radical, which can be realized
through simple cross multiplication of the dt differential.

Does the NAMD implementation of the adaptive tempering algorithm actually
contain this error, is it the online documentation that is flawed, or am I
missing something?

Any information in this regard would be much appreciated.

Best,

Sachin Natesh

This archive was generated by hypermail 2.1.6 : Thu Dec 31 2015 - 23:22:03 CST