On Sun, 15 Jun 2008 12:17:01 +0300
Vitaly Chaban <[EMAIL PROTECTED]> wrote:
I've got a discrepancy calculating the diffusion coefficient via
Green-Kubo equation and via Einstein equation.
For methanol (200ps, dt=0.001):
frequency to write velocities = 5
frequency to write coordinates = 100
Any rationale behind this difference? That might well result in differences
of C(t) at short times.
$ g_msd
2.3 (and it is true)
How much of msd(t) did you use? The part of the curve you use to fit
to Einstein equation will affect enormously the result!
$ g_velacc -nonormalize -acflen 1001
$ g_analyze -f vac.xvg -integrate
1.3 (and it's almost two times lower...)
Here you use 1001 points to integrate! Which result to 1001*0.005fs=5 ps!
Again did you use the same length of the msd(t) to fit via Einstein equation?
Although the difference in time scales you probably use might explain
the difference, that difference is actually small (same units?)!
Do I miss any important keyword here? It seems the result of integration is
not correct in this case.
Thanks,
Vitaly V. Chaban
_______________________________________________
gmx-users mailing list [email protected]
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
-----------------------------------------------------
XAvier Periole - PhD
NMR & Molecular Dynamics Group
University of Groningen
The Netherlands
http://md.chem.rug.nl/~periole
-----------------------------------------------------
_______________________________________________
gmx-users mailing list [email protected]
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
