[gmx-users] Minimisation

[Justin Lemkul]

On 5/22/17 7:25 AM, Pandya, Akash wrote:
> 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"
> #endif
>
> ; Include GLY Topology
> #include "GLY.itp"
> #ifdef POSRES
> #include "GLY_posre.itp"
> #endif
>
> [ moleculetype ]
> ; Name            nrexcl
> Protein_chain_H     3
>

A series of #include statements is not a useful diagnostic.  You should start by 
inspecting your coordinates, particularly whatever atom 6640 is, and its 
surroundings.  There may simply be a bad clash there that can be resolved 
easily.  Otherwise, follow 
http://www.gromacs.org/Documentation/Terminology/Blowing_Up#Diagnosing_an_Unstable_System 
  Glycine should come from the force field itself and likely is not 
topologically problematic, but if you constructed your own citrate parameters, 
you should check their correctness, as well.

-Justin

> Akash
>
> -----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
>
> --
> ==================================================
>
> 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.
>

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

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

==================================================