[gmx-users] Velocity autocorrelation
Ignacio Fernández Galván
jellby at yahoo.com
Fri Mar 30 12:41:27 CEST 2012
--- On Thu, 29/3/12, Krzysztof Kuczera <kkuczera at ku.edu> wrote:
> You probably need to save velocities
> at least every 10 fs;
> look at your plot on a fine enough time scale - say
> 0-1 ps range
> probably running in NVE ensemble would help as well
Thanks for your suggestion. I was already saving velocities every 10 fs. I tried with a NVE ensemble, but the behaviour is the same.
I have uploaded a sample velocity ACF here:
That's the autocorrelation file for CN atoms in acetonitrile.
I have also uploaded the corresponding mdp file:
I still don't know how to properly analyze this, I would say that I want the amplitude to drop to almost zero, and that would mean a relaxation time of ~5-6 ps. But is this the correct/expected behaviour?
More information about the gromacs.org_gmx-users