[gmx-users] Minimisation

Pandya, Akash akash.pandya.15 at ucl.ac.uk
Mon May 22 13:25:17 CEST 2017

What could be wrong with the topology could you please elaborate? I added glycine and citrate molecules  to the simulation box randomly. Could it be the fact I added topology for glycine and citrate in the wrong way?

This is how my topol.top file looks:

; Include Citrate Topology
#include "Citrate.itp"
#ifdef POSRES
#include "Citrate_posre.itp"

; Include GLY Topology
#include "GLY.itp"
#ifdef POSRES
#include "GLY_posre.itp"

[ moleculetype ]
; Name            nrexcl
Protein_chain_H     3


-----Original Message-----
From: gromacs.org_gmx-users-bounces at maillist.sys.kth.se [mailto:gromacs.org_gmx-users-bounces at maillist.sys.kth.se] On Behalf Of Justin Lemkul
Sent: 22 May 2017 12:13
To: gmx-users at gromacs.org
Subject: Re: [gmx-users] Minimisation

On 5/22/17 6:47 AM, Pandya, Akash wrote:
> Hi all,
> During Minimisation I get the following output that the simulation ended prematurely.
> Steepest Descents:
>    Tolerance (Fmax)   =  1.00000e+03
>    Number of steps    =       500000
> Step=    0, Dmax= 1.0e-02 nm, Epot=  7.39413e+26 Fmax=         inf, atom= 6640
> Step=   14, Dmax= 1.2e-06 nm, Epot=  7.39413e+26 Fmax=         inf, atom= 6640
> Energy minimization has stopped, but the forces have not converged to 
> the requested precision Fmax < 1000 (which may not be possible for your system).
> It stopped because the algorithm tried to make a new step whose size 
> was too small, or there was no change in the energy since last step. 
> Either way, we regard the minimization as converged to within the 
> available machine precision, given your starting configuration and EM parameters.
> Double precision normally gives you higher accuracy, but this is often 
> not needed for preparing to run molecular dynamics.
> You might need to increase your constraint accuracy, or turn off 
> constraints altogether (set constraints = none in mdp file)
> How do I increase my constraint accuracy? Which file do I have to change?

That's not your problem.  You have infinite force, which means you have some catastrophic problem with either your coordinates or topology.



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 at 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 at gromacs.org.

More information about the gromacs.org_gmx-users mailing list