Dear Minnale,
As far as i know, at least 100ns simulation should be performed in order
to obtain the correct diffusion coefficient for lipids. I think, if you
performed 5 ns it is too short.
Hope this helps.
balazs
minnale wrote:
Hi all,
I have read the literature about Diffusion cofficient(Dynamics in
atomistic simulations of phospholipid membranes: Nuclear magnetic
resonance relaxation rates and lateral diffusion, Vol-125,
page-204703, Journal-THE JOURNAL OF CHEMICAL PHYSICS)
then I have given command like this
g_msd -f 5ns_popc.xtc -s .tpr -n p4.ndx -o msd_popc.xvg -lateral z
In the index file i have selceted p8 for the centre of mass and it
will be like this
[ P8 ]
8 60 112 164 ......................
I got the diffusion cofficient value like this
MSD gathered over 5000 ps with 501 restarts
# Diffusion constants fitted from time 500 to 4500 ps
# D[ P8] = 0.0131 (+/- 0.0190) (1e-5 cm^2/s)
The actual D=(1-10)*10-7 cm2/s from above mentioned article but iam
getting the value 1e-5 cm^2/s
I have tried in different ways with options -type z but I didnt get
the value near to actual value.
Could you please tell me where I am doing mistake
Thanks alot in advance.
Rediff Shopping
<http://adworks.rediff.com/cgi-bin/AdWorks/click.cgi/www.rediff.com/signature-default.htm/[EMAIL PROTECTED]/2206641_2199021/2201651/1?PARTNER=3&OAS_QUERY=null>
No virus found in this incoming message.
Checked by AVG - http://www.avg.com
Version: 8.0.138 / Virus Database: 270.4.11/1553 - Release Date: 15/07/2008 05:48
------------------------------------------------------------------------
_______________________________________________
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
_______________________________________________
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