Calculation of Diffusion coefficient

From: Souvik Sinha (
Date: Mon Jan 29 2018 - 08:01:13 CST


Currently, I am estimating a peptide insertion (into membrane) free energy
profile using ABF. I want to calculate the local diffusion coefficient
along the reaction coordinate for which I need the estimation of 'time
autocorrelation function of the fluctuations of the instantaneous forces' *[
J. Phys. Chem. B, 2011, 115 (49), pp 14704–14712 ]*. Unfortunately, I
forgot to use the "*outputSystemForce*" option in the *'colvar' *block
during the ABF simulations.

>From a previous thread
I found that taking the assumption of similar force distribution and
autocorrelation time across different window will certainly narrow the
rerun option to a single window along a specific RC (for the purpose of
error analysis). But I think, the autocorrelation function of instantaneous
force will differ along RC because of the different environment (water to
lipid bilayer) and that is what, I guess, will be reflected in the values
of diffusion coefficient along RC. In that case, do I need to rerun all the
windows with "*outputSystemForce*" option enabled (which is pretty
troublesome right now)? Is there any post-processing option from which I
can extract the system force trajectory for all the windows?


Souvik Sinha
Research Fellow
Bioinformatics Centre (SGD LAB)
Bose Institute
Contact: 033 25693275

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2018 - 23:20:49 CST