<br><br><div class="gmail_quote">On Wed, Jan 19, 2011 at 4:40 PM, Xiaohu Li <span dir="ltr"><<a href="mailto:xiaohuli914@gmail.com">xiaohuli914@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br><div class="gmail_quote"><div class="im"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Message: 2<br>
Date: Wed, 19 Jan 2011 23:23:12 +0100<br>
From: Berk Hess <<a href="mailto:gmx3@hotmail.com" target="_blank">gmx3@hotmail.com</a>><br>
Subject: RE: [gmx-users] RE: Important: Bugs in NEMD calculation<br>
To: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
Message-ID: <COL113-W651FF8B51B6E10B91E0E4F8EF60@phx.gbl><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
<br>
<br>
<br>
It seems you are right.<br>
The density is calculated a few lines before 1/visc is stored, but it is not being used.<br>
Strange that this always went unnoticed, since it already seems to have been wrong in 4.0.<br>
The density of water is so close to 1000 that you wouldn't notice the difference<br>
<br>
But in a quick test I don't get the correct viscosity either way.<br>
I will have to find time for a thorough check.<br>
Do you get the right answer by dividing by the missing density term in src/mdlib/mdebin.c?<br>
<br>
Berk<br></blockquote></div><div>Well, my force field can not reproduce the experiment anyway. <br>However, the wrong way of calculating in NEMD(vol instead of dens) is getting visocisty much much lower than what I got from using g_energy -vis.<br>
If I'm understanding correctly, the one with visco.xvg(column 1 vs column 2) should be the one I'm supposed to look at, even with 12ns simulation, this does not converge, but roughly speaking the NEMD results with the vol replaced by rho(I wrote a outside script to calculate this) is on the same order of magnitude as the g_energy -vis using green-kudo(is that right?) expression.<br>
<font color="#888888">
<br>Xiaohu<br> </font></div><div><div></div><div class="h5"><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
Date: Wed, 19 Jan 2011 14:48:32 -0600<br>
From: <a href="mailto:xiaohuli914@gmail.com" target="_blank">xiaohuli914@gmail.com</a><br>
To: <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
Subject: [gmx-users] RE: Important: Bugs in NEMD calculation<br>
<br>
<br>
<br>
Date: Wed, 19 Jan 2011 20:43:20 +0100<br>
<br>
From: Berk Hess <<a href="mailto:gmx3@hotmail.com" target="_blank">gmx3@hotmail.com</a>><br>
<br>
Subject: RE: [gmx-users] Important: Bugs in NEMD calculation<br>
<br>
To: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
<br>
Message-ID: <COL113-W2019E7174AFFBAEF819A1A8EF60@phx.gbl><br>
<br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
> Date: Wed, 19 Jan 2011 19:13:12 +0100<br>
<br>
> From: <a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a><br>
<br>
> To: <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<br>
> Subject: Re: [gmx-users] Important: Bugs in NEMD calculation<br>
<br>
><br>
<br>
> On 2011-01-19 18.36, Xiaohu Li wrote:<br>
<br>
> > Hi, All,<br>
<br>
> > I've found a bug in the NEMD calculation for viscosity. This has<br>
<br>
> > been reviewed in /*Hess's paper at JCP 116 209 2002.*/<br>
<br>
> > The version of gromacs I'm using is the development version.<br>
<br>
> > Notice that this version correct a previous but of 4.5.3, where you uses<br>
<br>
> > NEMD, both the term V(eq. 21 on Hess's paper) and 1/eta(shear viscosity<br>
<br>
> > inverse) are supposed to be written to the *.edr file, however,<br>
<br>
> > the 4.5.3 versions does not have this. This version can be retrieved at<br>
<br>
> > <a href="http://repo.or.cz/w/gromacs.git/commit/c83de86d65ce7135be6cef15e9100d7516e6d9a7" target="_blank">http://repo.or.cz/w/gromacs.git/commit/c83de86d65ce7135be6cef15e9100d7516e6d9a7</a><br>
<br>
> > *However, even this version is buggy since eta=A*rho/(V*k**2)(eq. 20<br>
<br>
> > Hess's paper)*. *I have performed simulation and has found out that the<br>
<br>
> > V and eta which are written in *edr file does not match according to the<br>
<br>
> > formula, a little further check on the source code mdebin.c under the<br>
<br>
> > directory src/mdlib shows that it is actually calculating<br>
<br>
> > eta=A*Volume/(V*k**2) where density of rho should have been used. (This<br>
<br>
> > is at line 755 of mdebin.c ).<br>
<br>
> > I hope everyone who is using this can be aware of this, if you ever<br>
<br>
> > used this code to produce data, the V is correct from *edr, however, you<br>
<br>
> > need to manuelly get your eta using the above formula.<br>
<br>
> > For the GMX developers, I hope anyone of you can correct this bug.<br>
<br>
> ><br>
<br>
<br>
<br>
I fixed this bug recently is the git release-4-5-patches branch.<br>
<br>
You can get the fix from git and it will be in the 4.5.4 release<br>
<br>
(no date yet).<br>
Which one are you referring to? The one I got is the one you uploaded that fixed the zero viscosity and 2*cosZ*vel-x in edr file.<br>
This bug refers the wrong calculation of 1/eta.<br>
<br>
<br>
<br>
<br>
<br>
> Thanks for pointing that out. There also is a small issue in that the<br>
<br>
> volume is computed for a rectangular box<br>
<br>
> vol = box[XX][XX]*box[YY][YY]*box[ZZ][ZZ];<br>
<br>
> dens = (tmass*AMU)/(vol*NANO*NANO*NANO);<br>
<br>
><br>
<br>
> which would be incorrect for a non-rectangular box. You should however<br>
<br>
> use a rectangular box for this kind of calculations, although this is<br>
<br>
> not enforced by grompp.<br>
<br>
<br>
<br>
No.<br>
<br>
That formula is correct for any triclinic box!<br>
<br>
<br>
<br>
Berk<br>
<br>
<br>
<br>
><br>
<br>
> ><br>
<br>
> ><br>
<br>
> > Xiaohu<br>
<br>
> > *<br>
<br>
> ><br>
<br>
><br>
<br>
><br>
<br>
> --<br>
<br>
> David van der Spoel, Ph.D., Professor of Biology<br>
<br>
> Dept. of Cell & Molec. Biol., Uppsala University.<br>
<br>
> Box 596, 75124 Uppsala, Sweden. Phone: +46184714205.<br>
<br>
> <a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a> <a href="http://folding.bmc.uu.se" target="_blank">http://folding.bmc.uu.se</a><br>
<br>
> --<br>
<br>
> gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<br>
> <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
<br>
> Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
<br>
> Please don't post (un)subscribe requests to the list. Use the<br>
<br>
> www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
<br>
> Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
<br>
<br>
-------------- next part --------------<br>
<br>
An HTML attachment was scrubbed...<br>
<br>
URL: <a href="http://lists.gromacs.org/pipermail/gmx-users/attachments/20110119/1ffc9f59/attachment.html" target="_blank">http://lists.gromacs.org/pipermail/gmx-users/attachments/20110119/1ffc9f59/attachment.html</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
<br>
<br>
--<br>
<br>
gmx-users mailing list<br>
<br>
<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
<br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
<br>
<br>
<br>
End of gmx-users Digest, Vol 81, Issue 127<br>
<br>
******************************************<br>
<br>
<br>
<br>
<br>
<br>
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://lists.gromacs.org/pipermail/gmx-users/attachments/20110119/78703d32/attachment-0001.html" target="_blank">http://lists.gromacs.org/pipermail/gmx-users/attachments/20110119/78703d32/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Wed, 19 Jan 2011 17:29:02 -0500<br>
From: "Justin A. Lemkul" <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>><br>
Subject: Re: [gmx-users] CNT<br>
To: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
Message-ID: <<a href="mailto:4D3765AE.6080307@vt.edu" target="_blank">4D3765AE.6080307@vt.edu</a>><br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
<br>
<br>
trevor brown wrote:<br>
> Dear Justin,<br>
> Regarding our previous discussion under "CNT" topic, you want to know<br>
> exactly what I did to define C-C interaction. I have summarized it below.<br>
><br>
> 1. I copied oplsaa.ff folder in my working directory<br>
> 2. I added following lines to atomname2type.n2t<br>
> C opls_995 0 12.011 2 C 0.142 C 0.142<br>
> C opls_996 0 12.011 3 C 0.142 C 0.142 C 0.142<br>
> C opls_997 0 12.011 4 C 0.142 C 0.142 C 0.142 C 0.142<br>
> C opls_998 0 12.011 5 C 0.142 C 0.142 C 0.142 C 0.142<br>
> C 0.142<br>
> 3. I added those to atomtypes.atp<br>
> opls_995 12.01100 ;<br>
> opls_996 12.01100 ;<br>
> opls_997 12.01100 ;<br>
> opls_998 12.01100 ;<br>
><br>
> 4. I added those to ffbonded.itp<br>
> [ bondtypes ]<br>
> C C 1 0.14210 478900<br>
> [ angletypes ]<br>
> C C C 1 120.000 562.2<br>
> [ dihedraltypes ]<br>
> 5. I constructed .top file for CNT with the following, then I converted<br>
> it into .itp by vi editor<br>
> g_x2top -f cnt.gro -o cnt.top -nopairs -nexcl 5<br>
> 6. By pdb2gmx, I obtained .top for peptide then I converted it into .itp<br>
> 7. I combined CNT and peptide's pdbs in pymol and saved.<br>
> 8. I wrote a .top file given below,<br>
> ;<br>
> ; File 'topol.top' was generated<br>
> ; By user: onbekend (0)<br>
> ; On host: onbekend<br>
> ; At date: Mon Jan 10 02:51:19 2011<br>
> ;<br>
> ; This is a standalone topology file<br>
> ;<br>
> ; It was generated using program:<br>
> ; pdb2gmx - VERSION 4.5.3<br>
> ;<br>
> ; Command line was:<br>
> ; pdb2gmx -ignh -f vpgvg10.pdb<br>
> ;<br>
> ; Force field was read from the standard Gromacs share directory.<br>
> ;<br>
> ; Include forcefield parameters<br>
> #include "forcefield.itp"<br>
><br>
<br>
Be careful with this - if you move this topology between directories, it will<br>
call any forcefield.itp file that may be in the working directory. What did you<br>
remove the "oplsaa.ff" prefix?<br>
<br>
> ; Include topology for CNT<br>
> #include "cnt.itp"<br>
> ; Include topology for UW1<br>
> #include "uw1.itp"<br>
> ; Include water topology<br>
> #include "spc.itp"<br>
> [ system ]<br>
> ; Name<br>
> Protein and CNT in water<br>
> [ molecules ]<br>
> ; Compound #mols<br>
> Protein 1<br>
> CNT 1<br>
> SOL 4019<br>
> 9. When I make grompp for minimization it givea me the following error<br>
> WARNING 1 [file ffbonded.itp, line 2703]:<br>
> Overriding Bond parameters.<br>
> old: 0.151 292880 0.151 292880<br>
> new: C C 1 0.14210 478900<br>
><br>
> ERROR 1 [file ffbonded.itp, line 2707]:<br>
> Not enough parameters<br>
> Generated 334153 of the 334153 non-bonded parameter combinations<br>
> Generating 1-4 interactions: fudge = 0.5<br>
> Generated 334153 of the 334153 1-4 parameter combinations<br>
> -------------------------------------------------------<br>
> Program grompp, VERSION 4.5.3<br>
> Source code file: toppush.c, line: 1166<br>
> Fatal error:<br>
> Atomtype opls_996 not found<br>
><br>
> That is all what I exactly did. Is anything wrong or missing?<br>
><br>
<br>
You never defined nonbonded parameters for your new atom type. You've declared<br>
that it exists in the .atp file, but then you didn't add that atom type in<br>
ffnonbonded.itp.<br>
<br>
-Justin<br>
<br>
> best wishes<br>
> trevor<br>
><br>
><br>
><br>
><br>
<br>
--<br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br>
<br>
<br>
------------------------------<br>
<font color="#888888"><br>
--<br>
gmx-users mailing list<br>
<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
<br>
End of gmx-users Digest, Vol 81, Issue 130<br>
******************************************<br>
</font></blockquote></div></div></div><br>
</blockquote></div><br>