<div dir="ltr">Hi all,<br><br>I am trying to simulate a electroosmotic flow in a nanochannel. System consists of water + (Cl-) ions in a channel confined of silicon <111> walls. I start with taking water as only solvent, equilibrating the system for about 12-15 ps (gmxcheck runs clean), and then use genion to insert (Cl-) ions. Walls are charged appropriately to maintain electro-neutrality in the system. The problem is with Energy minimization after inserting ions using steep algorithm (haven't turned on the electric field yet). I am getting an error like this:<br>
<br>t = 0.038 ps: Water molecule starting at atom 2569 can not be settled.<br>Check for bad contacts and/or reduce the timestep.<br>Wrote pdb files with previous and current coordinates<br><br>and<br><br>Steepest Descents converged to machine precision in 88 steps,<br>
but did not reach the requested Fmax < 1000.<br>Potential Energy = 1.75456106043895e+06<br>Maximum force = 1.07111309902475e+11 on atom 4898<br>Norm of force = 1.48772609239112e+09<br><br>If it were only couple of such molecules, I could have removed them easily from the gro file. The problem is twofold, there are around 10 molecules which give this error message, plus even after removal of these molecules, energy minimization results in the same error on different sets of molecules. I have tried a number of ways to tackle this including tuning the time step. Is there any way to avoid this? Any help would be appreciated.<br>
<br>Thanks in advance.<br clear="all"><br>-- <br>Ravi Bhadauria<br>
</div>