<p style="border:0; padding:0; margin:0; font-family:'Gulim'; font-size:10pt; cursor: text;"><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><span style="font-family: monospace; font-size: large; line-height: 24px;">Dear Justin and all,</span></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><div style="font-family: 돋움, arial; font-size: 12px; margin-top: 30px; margin-left: 0.8em; color: rgb(0, 102, 204); font-weight: bold;">-------- 원본 메일 ---------</div><p style="font-family: monospace; font-size: 12pt; line-height: 1.5; margin-bottom: 0px; margin-top: 0px;"><blockquote style="font-family: Gulim; font-size: 12px; border-left-style: solid; border-left-width: 2px; margin: 0pt 0pt 0pt 0.8em; padding-left: 1em;"><div style="font-family:arial,돋움;line-height:1.5"><b>보낸사람</b> : Justin Lemkul <jalemkul@vt.edu><br><b>받는사람</b> : <gmx-users@gromacs.org><br><b>받은날짜</b> : 2014년 9월 29일(월) 21:32:13<br><b>제목</b> : Re: [gmx-users] Problem with constraints in NVT calculations.<br><!-- original content --><div style="margin-top:5px;"><pre>On 9/29/14 8:01 AM, Kester Wong wrote:
> Dear all,
>
>
> I have been using the SWM4-NDP water model that works well with CHARMM27 force
> field.
>
Time out. SWM4-NDP is a polarizable water model, and CHARMM27 is an additive
force field. The two are not compatible. If you're using them together, you're
basically just gambling that they work. You shouldn't do this. </pre></div></div></blockquote></p><p style="font-family: monospace; font-size: 12pt; line-height: 1.5; margin-bottom: 0px; margin-top: 0px;"><font face="monospace">Thank you for the information. </font><span style="font-size: 12pt; line-height: 1.5;">Actually, I took the data from the ff-drude-water-ions.itp file and included them in the CHARMM27 itp files. Now I realised that it does not work like that.</span></p><p style="margin-bottom: 0px; margin-top: 0px;"><br></p><p style="font-family: monospace; font-size: 12pt; line-height: 1.5; margin-bottom: 0px; margin-top: 0px;">If the ffnonbonded.itp file contains the appropriate nonbonded parameters for the polarisable water model (i.e. the atomtypes and dummies for the Drude water and ions model), will it work? </p><p style="font-family: monospace; font-size: 12pt; line-height: 1.5; margin-bottom: 0px; margin-top: 0px;">Below is my input topologies:</p><p style="font-family: monospace; font-size: 12pt; line-height: 1.5; margin-bottom: 0px; margin-top: 0px;"><span style="font-size: medium;">https://drive.google.com/file/d/0B7ym8d6G9-e2dmdQOHNaU3c3Qkk/view?usp=sharing</span><font face="monospace"><br></font><font face="monospace"><span style="font-size: 12px; line-height: 18px; white-space: pre;"><br></span></font><blockquote style="font-family: Gulim; font-size: 12px; border-left-style: solid; border-left-width: 2px; margin: 0pt 0pt 0pt 0.8em; padding-left: 1em;"><div style="font-family:arial,돋움;line-height:1.5"><div style="margin-top:5px;"><pre>> The force field (most recent version) was obtained from virtualchemistry, and
> consists of polarisable ions and the SWM4-NDP water model from Lamoureux et al.
>
> The water on graphene (energy minimisation and NVT) calculations are stable, and
> have produced a few NVT runs of up to 10ns.
>
> I also compared the NVT runs with other previously calculated structures,
> namely, the visualised NVT runs are very much comparable to that obtained by
> CHARMM27 force field with TIPS3P water model.
>
>
> However, as the polarisable H3O+ and Cl- ions are inserted into the simulation
> box, the NVT run tends to crash as a result of too many LINCS warnings.
>
> I suspect that the constraints in the H3O is the cause of the LINCS warnings,
> and after removing the [ constraints ] in the topology file, the calculation
> seems to be running.
>
>
> ;[ constraints ]
>
> ;; i funct doh dhh
>
> ;1 3 1 0.102
>
> ;1 4 1 0.102
>
> ;1 5 1 0.102
>
> ;3 4 1 0.169124
>
> ;4 5 1 0.169124
>
> ;3 5 1 0.169124
>
>
>
> Can anyone please tell me if what I am doing is correct?
>
> As I am not sure if removing the H3O constraints in the NVT run is appropriate.
>
>
> With regard to the massless Drude and standard oxygen mass (thank you Justin for
> providing the information), I would like to know why are they not appropriate
> (as mentioned in the previous email) for the SWM4-NDP model?
>
I wasn't making a declarative statement that it won't work. I was asking
whether anyone had verified that using a massless Drude in that model actually
reproduced all of the physical properties of the SWM4-NDP model. We do all of
our development on the Drude FF with Drudes of mass = 0.4 amu and the extended
Lagrangian. I don't know how a massless Drude behaves in this model. It may
work fine, but you should very rigorously verify this yourself before doing
anything by doing simulations of pure water.
-Justin
</pre></div></div></blockquote><font face="monospace">I have Read the paper by Jochen Hub et al [Chem. Sci., 5, 1745 2013], that their calculated data (water with ions) are in line with experimental results.<br>Upon further inspection, I found more papers that reported the use of 0.4 amu mass and an extended Lagrangian for the Drude particle. e.g. by Yu et al. [J. CHem. Theory Comput. 6, 774-786 2010]. <br>Could you please advise if there is anything I could do, to make the additive and polarisable FFs usable after modification?</font></p><p style="margin-bottom: 0px; margin-top: 0px;"><font face="monospace"><font size="4"><span style="line-height: 1.5;">I am more than willing to go through the </span><span style="line-height: 24px;">necessary</span><span style="line-height: 1.5;"> changes/implementation if it can be done as this is the only avenue that I have left. </span></font></font></p><p style="margin-bottom: 0px; margin-top: 0px;"><br></p><p style="margin-bottom: 0px; margin-top: 0px;"><br></p><p style="margin-bottom: 0px; margin-top: 0px;"><span style="font-family: monospace; font-size: large; line-height: 24px;">Thank you in advance.</span></p><p style="margin-bottom: 0px; margin-top: 0px;"><br></p><p style="margin-bottom: 0px; margin-top: 0px;"><span style="font-family: monospace; font-size: large; line-height: 24px;">Regards,</span></p><p style="margin-bottom: 0px; margin-top: 0px;"><span style="font-family: monospace; font-size: large; line-height: 24px;">Kester</span><font face="monospace"><br><br><span style="font-family: monospace; font-size: 12px; line-height: 18px; white-space: pre;"><br></span></font><blockquote style="font-family: Gulim; font-size: 12px; line-height: 1.5; border-left-style: solid; border-left-width: 2px; margin: 0pt 0pt 0pt 0.8em; padding-left: 1em;"><div style="font-family:arial,돋움;line-height:1.5"><div style="margin-top:5px;"><pre><br></pre><pre><br></pre><pre><br></pre><pre>--
==================================================
Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow
Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201
jalemkul@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul
==================================================
--
Gromacs Users mailing list
* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request@gromacs.org.
</pre></div><!-- original content --><br></div></blockquote></p><p style="font-family: Gulim; font-size: 10pt;"></p><p style="font-family: Gulim; font-size: 10pt;"></p><p style="font-family: Gulim; font-size: 10pt;"></p><p style="font-family: Gulim; font-size: 10pt;"></p><p style="font-family: Gulim; font-size: 10pt;"></p></p>
<img src='http://mail.ibs.re.kr/checkread/ODIxMDMz/Z214LXVzZXJzQGdyb21hY3Mub3Jn/' width='1px' height='1px' />