<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 all,&nbsp;</p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;">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.</p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;">Using emtol = 1000 and emstep = 0.0001, I get:</p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><p style="border: 0px; padding: 0px; margin: 0px; cursor: text;">Potential Energy &nbsp;= &nbsp;3.79607808453231e+09</p><p style="border: 0px; padding: 0px; margin: 0px; cursor: text;">Maximum force &nbsp; &nbsp; = &nbsp;9.86741557582003e+02 on atom 42394</p><p style="border: 0px; padding: 0px; margin: 0px; cursor: text;">Norm of force &nbsp; &nbsp; = &nbsp;2.38846729392462e+01</p><div><br></div><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><span style="font-size: 10pt;">Using emtol = 700 and emstep = 0.0001, I get:</span></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><p style="border: 0px; padding: 0px; margin: 0px; cursor: text;">Potential Energy &nbsp;= -6.43566327273198e+13</p><p style="border: 0px; padding: 0px; margin: 0px; cursor: text;">Maximum force &nbsp; &nbsp; = &nbsp;9.24568347101016e+25 on atom 42262</p><p style="border: 0px; padding: 0px; margin: 0px; cursor: text;">Norm of force &nbsp; &nbsp; = &nbsp;5.74744822148669e+23</p><div style="font-family: Gulim; font-size: 10pt;"><br></div><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;">As expected, both structures blew up during NVT equilibration.</p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;">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.&nbsp;</p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;">FYI, the water droplet on graphene was equilibrated for 20ns.&nbsp;</p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;">Any input is greatly appreciated. Thanks.</p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;"><br></p><p style="font-family: Gulim; font-size: 10pt; border: 0px; padding: 0px; margin: 0px; cursor: text;">- Kester</p></p>
<img src='http://mail.ibs.re.kr/checkread/OTM1NzIz/Z214LXVzZXJzQGdyb21hY3Mub3Jn/' width='1px' height='1px' />