ABF question

From: Parisa Akhski (Parisa.Akhshi_at_chem.queensu.ca)
Date: Tue Jun 22 2010 - 11:43:06 CDT

Dear Jerome and other ABF method users,

I have a couple of questions using ABF method for an ion transport process:

(1) I have generated a PMF profile and I want to know whether the results are reliable or not? One criteria I think is the smoothness of the sampling. Also, according to the tutorial, history file is useful for assessing the convergence of the free energy profile. I am assuming that I need to plot the free energy versus time step for all individual bins to see whether they all converge or not. Is it exactly what it means? Is the force convergence another necessity for reliable results?

(2) I am a little confused about the sampling in constrained and unconstrained MD simulations. Since in constrained MD we freeze the reaction coordinate at discrete values of reaction coordinates (RC) and we have to specify the exact positions for RC and the number of bins and also the total number of "run" in configuration file (for example 5000000), we are determining the exact number of samples in each bin (equal number of samples in all bins) . And because we determine the time step as well, we have the control on the time that the particle spends on each bin. In unconstrained simulation, however, we have no control on the particle positions and the ion moves freely and the sampling is random. If these are true, can we conclude that in constrained MD the ion starts from the first bin and after spending the defined time, it moves to the next bin and...(the ion movement happens in just one direction starting form one end and ending at the other end) whereas, in unconstrained, it moves forth and back (without following just one direction)? Finally, can we say that the ns time scale in unconstrained MD (ABF) is just for the whole process not in each individual bin?

I apologize for these simple questions but I am really confused!



This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:54:16 CST