What are the atoms that are causing the problems actually doing, where are they located etc?
Catch ya, Dr. Dallas Warren Drug Delivery, Disposition and Dynamics Monash Institute of Pharmaceutical Sciences, Monash University 381 Royal Parade, Parkville VIC 3010 dallas.war...@pharm.monash.edu.au +61 3 9903 9167 --------------------------------- When the only tool you own is a hammer, every problem begins to resemble a nail. From: gmx-users-boun...@gromacs.org [mailto:gmx-users-boun...@gromacs.org] On Behalf Of Ravi Bhadauria Sent: Thursday, 8 April 2010 3:11 PM To: gmx-users@gromacs.org Subject: [gmx-users] energy minimization in charged systems Hi all, 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: t = 0.038 ps: Water molecule starting at atom 2569 can not be settled. Check for bad contacts and/or reduce the timestep. Wrote pdb files with previous and current coordinates and Steepest Descents converged to machine precision in 88 steps, but did not reach the requested Fmax < 1000. Potential Energy = 1.75456106043895e+06 Maximum force = 1.07111309902475e+11 on atom 4898 Norm of force = 1.48772609239112e+09 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. Thanks in advance. -- Ravi Bhadauria -- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/mailing_lists/users.php
