<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;">Dear Justin and Erik,&nbsp;</p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><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>&gt; Thank you for the feedback. I have also tried the following:
&gt;
&gt; i) separating OH- and Na+ ions
&gt;
&gt; ii) placing the ions much closer to water droplet
&gt;
&gt; iii) increase the spacing between each ions, i.e. using a larger box
&gt;
&gt;
&gt; Perhaps I should try having the ions in a larger box,  then fill the box with
&gt; water molecules and then redo the energy minimisation and NVT?
&gt;

Well, like Erik said, there really isn't a smoking gun and what I'm observing 
might be either causative or symptomatic, there's just no way to know at this 
point.  I don't know what good any of those three options will do, because I 
don't understand what you're trying to accomplish.  What was the logic of 
putting a lattice of OH- and random Na+ in vacuo above a water droplet?  What is 
the goal of your simulations?</pre><pre><br></pre></div></div></blockquote><p style="font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;">Thank you for the valuable input, the goal is to simulate the ions in water droplet.</p><p style="font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;">Instead of having a box of water with the charged ions, I used the water/graphene system that has been NVT equilibrated for 20ns, and placed the ions on top of the droplet.</p><p style="font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><p style="font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;">If adding 20 OH- and 20 Na+ proves to be problematic, is there a way to add the ions incrementally without attracting a large force?</p><p style="font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><p style="font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;">The OH- ions were arranged in such a way "genconf -nbox 4 5 1" as I have tried random insertion, using gmx insert-molecules, and have found that the ions were too close together.</p><p style="font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><font face="monospace">From my other previous calculations, adding a cluster of H3O+ and Cl- ions (ranging from 40 to 200 ions) on top of the water droplet worked fine with NVE and NVT.<br>However, I have not been able to do so with OH- and Na+.<br><br><br><br>Regards,<br>Kester<br><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>
-Justin

&gt;
&gt;
&gt; Regards,
&gt;
&gt; Kester
&gt;
&gt;     On 11/10/14 5:04 AM, Erik Marklund wrote:
&gt;     &gt; Dear Kester,
&gt;     &gt;
&gt;     &gt; Hm. No obvious smoking gun as far as I can see. You could try a two-step minimisation where you first do one minimisation with soft-core potentials and then one without.
&gt;     &gt;
&gt;
&gt;     What's strange to me is the lattice of OH- at the top of the box and the sodium
&gt;     ions that are floating out in the middle of nowhere.  I don't understand the
&gt;     construction of this system, but having highly charged stuff in close contact
&gt;     and floating around in vacuum certainly looks like it could be trouble.
&gt;
&gt;     -Justin
&gt;
&gt;     &gt; Kind regards,
&gt;     &gt; Erik
&gt;     &gt;
&gt;     &gt; On 7 Nov 2014, at 05:36, Kester Wong&gt; wrote:
&gt;     &gt;
&gt;     &gt;
&gt;     &gt; Dear Erik,
&gt;     &gt;
&gt;     &gt; The starting structure consists of water on graphene that has already been energy minimised, and subsequently equilibrated (NVT) for 20ns.
&gt;     &gt; Using that (water droplet on graphene), I added the ions for energy minimisation.
&gt;     &gt; I followed a two-step minimisation:
&gt;     &gt; i) Minimisation with no constraints (-DFLEXIBLE)
&gt;     &gt; ii) Minimisation with constraints (SETTLES for TIPS3P water, and (-DCONSTRAINTS) for hydroxide)
&gt;     &gt;
&gt;     &gt; With the ions:
&gt;     &gt; I have tried placing the ions box lower toward water droplet, but the Max Force remains large upon energy minimisation.
&gt;     &gt;
&gt;     &gt; The input files can be accessed here:
&gt;     &gt; https://drive.google.com/folderview?id=0B7ym8d6G9-e2dG5pOFQ4SWVyZG8&amp;usp=sharing
&gt;     &gt;
&gt;     &gt; The hydroxide parameter is taken from Gerrit Groenhof, based on the supplementary material in this paper:
&gt;     &gt; http://dx.doi.org/10.1016/j.bpj.2014.04.062
&gt;     &gt;
&gt;     &gt; In short, this model was developed very similarly to the polarisable version of hydroxide in the SWM4-NDP force field.
&gt;     &gt; As for sodium and graphene, the topology parameters were taken from CHARMM27; TIPS3P was used for the water model.
&gt;     &gt;
&gt;     &gt; Visualisation after the energy minimisation showed a distorted virtual site of the hydroxide ions, where one of the four virtual sites no longer conform to the constraint. So, I am guessing whether the extremely large Max Force originates from the force constants in the hydroxide topology?
&gt;     &gt;
&gt;     &gt; Thank you for your time and assistance!
&gt;     &gt; -Kester
&gt;     &gt; --------- 원본 메일 ---------
&gt;     &gt; 보낸사람 : Erik Marklund&gt;
&gt;     &gt; 받는사람 : "&gt;"&gt;
&gt;     &gt; 받은날짜 : 2014년 11월 6일(목) 20:35:59
&gt;     &gt; 제목 : Re: [gmx-users] Very large Max Force [Energy minimisation] in water with ions
&gt;     &gt;
&gt;     &gt; Dear Kester,
&gt;     &gt;
&gt;     &gt; The potential energy is highly positive in the first case and the force is enormous in the second case, so no wonder that they blow up. How did you prepare these systems?
&gt;     &gt;
&gt;     &gt; Kind regards,
&gt;     &gt; Erik
&gt;     &gt;
&gt;     &gt; Erik Marklund, PhD
&gt;     &gt; Postdoctoral Research Fellow, Fulford JRF
&gt;     &gt;
&gt;     &gt; Department of Chemistry
&gt;     &gt; Physical &amp; Theoretical Chemistry Laboratory
&gt;     &gt; University of Oxford
&gt;     &gt; South Parks Road
&gt;     &gt; Oxford
&gt;     &gt; OX1 3QZ
&gt;     &gt;
&gt;     &gt; On 6 Nov 2014, at 04:31, Kester Wong &gt; wrote:
&gt;     &gt;
&gt;     &gt;
&gt;     &gt; Dear all,
&gt;     &gt;
&gt;     &gt; I have been trying to energy minimise ~1980 water (tips3p) molecules with 20 OH- anions (with virtual sites and dummy atoms) and 20 Na+ cations, however, I always ended up getting a very large Max Force.
&gt;     &gt;
&gt;     &gt; Using emtol = 1000 and emstep = 0.0001, I get:
&gt;     &gt;
&gt;     &gt; Potential Energy  =  3.79607808453231e+09
&gt;     &gt; Maximum force     =  9.86741557582003e+02 on atom 42394
&gt;     &gt; Norm of force     =  2.38846729392462e+01
&gt;     &gt;
&gt;     &gt; Using emtol = 700 and emstep = 0.0001, I get:
&gt;     &gt;
&gt;     &gt; Potential Energy  = -6.43566327273198e+13
&gt;     &gt; Maximum force     =  9.24568347101016e+25 on atom 42262
&gt;     &gt; Norm of force     =  5.74744822148669e+23
&gt;     &gt;
&gt;     &gt; As expected, both structures blew up during NVT equilibration.
&gt;     &gt;
&gt;     &gt; Could this be due to the arrangement of the ions? I tried separating the anions and cations in two boxes (placed on top of the water droplet), but that also yield a very large Max Force.
&gt;     &gt;
&gt;     &gt; FYI, the water droplet on graphene was equilibrated for 20ns.
&gt;     &gt; Any input is greatly appreciated. Thanks.
&gt;     &gt;
&gt;     &gt; - Kester
&gt;     &gt;
&gt;     &gt; --
&gt;     &gt; Gromacs Users mailing list
&gt;     &gt;
&gt;     &gt; * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
&gt;     &gt;
&gt;     &gt; * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
&gt;     &gt;
&gt;     &gt; * For (un)subscribe requests visit
&gt;     &gt; https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request@gromacs.org.
&gt;     &gt;
&gt;     &gt; --
&gt;     &gt; Gromacs Users mailing list
&gt;     &gt;
&gt;     &gt; * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
&gt;     &gt;
&gt;     &gt; * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
&gt;     &gt;
&gt;     &gt; * For (un)subscribe requests visit
&gt;     &gt; https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request@gromacs.org.
&gt;     &gt;
&gt;     &gt;
&gt;     &gt;
&gt;     &gt;
&gt;     &gt; --
&gt;     &gt; Gromacs Users mailing list
&gt;     &gt;
&gt;     &gt; * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
&gt;     &gt;
&gt;     &gt; * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
&gt;     &gt;
&gt;     &gt; * For (un)subscribe requests visit
&gt;     &gt; https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request@gromacs.org.
&gt;     &gt;
&gt;
&gt;     --
&gt;     ==================================================
&gt;
&gt;     Justin A. Lemkul, Ph.D.
&gt;     Ruth L. Kirschstein NRSA Postdoctoral Fellow
&gt;
&gt;     Department of Pharmaceutical Sciences
&gt;     School of Pharmacy
&gt;     Health Sciences Facility II, Room 629
&gt;     University of Maryland, Baltimore
&gt;     20 Penn St.
&gt;     Baltimore, MD 21201
&gt;
&gt;     jalemkul@outerbanks.umaryland.edu | (410) 706-7441
&gt;     http://mackerell.umaryland.edu/~jalemkul
&gt;
&gt;     ==================================================
&gt;     --
&gt;     Gromacs Users mailing list
&gt;
&gt;     * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
&gt;
&gt;     * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
&gt;
&gt;     * For (un)subscribe requests visit
&gt;     https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request@gromacs.org.
&gt;
&gt;
&gt;
&gt;

-- 
==================================================

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
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>
<img src='http://mail.ibs.re.kr/checkread/OTU1Mzcx/Z214LXVzZXJzQGdyb21hY3Mub3Jn/' width='1px' height='1px' />