Re: [AMBER] Error in Minimization

From: Bill Ross <ross.cgl.ucsf.edu>
Date: Fri, 8 Mar 2019 02:21:21 -0800

Step 1: close/overlapping vdw for 2+ atoms, the most prominent being C,
1441: look at that atom and its neighbors, especially nonbonded
neighbors, especially any atoms you supplied params for.

Step 50: Now 1-4 (bonded) vdw's are impacted. C5 1446 == in same
neighborhood, might be the 'other atom', assuming only two are tangoing,
as is usually the case.


On 3/8/19 1:18 AM, Zhonghua Xia wrote:
> Dear all,
>
> I want to run conventional MD simulation on a protein complex with a small molecule. But it has been broken in the first step. The .in and .out file are shown below. I am wondering what is the problem and how to explain it. I appreciate any help from you:)
>
> Initial minimisation with Cartesian restrains for the solvent (hold the protein and Cl-)
> &cntrl
> imin=1, maxcyc=3000,
> ntr=1,
> restraint_wt=100,
> restraintmask=':1-90'
> &end
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 1 4.3472E+08 2.4011E+02 1.9793E+04 C 1441
>
> BOND = 9784.2491 ANGLE = 1786.4015 DIHED = 934.9579
> VDWAALS = ************* EEL = -62579.0562 HBOND = 0.0000
> 1-4 VDW = 1083.6816 1-4 EEL = 4202.4250 RESTRAINT = 0.0000
>
> .... RESTARTED DUE TO LINMIN FAILURE ...
>
> .... RESTARTED DUE TO LINMIN FAILURE ...
>
> .... RESTARTED DUE TO LINMIN FAILURE ...
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 50 -8.3706E+09 8.5391E+04 8.3864E+06 C5 1446
>
> BOND = 12553.8921 ANGLE = 2862.0259 DIHED = 952.0556
> VDWAALS = ************* EEL = -64213.1361 HBOND = 0.0000
> 1-4 VDW = ************* 1-4 EEL = 4191.5379 RESTRAINT = 3088.1670
> EAMBER = *************
>
> .... RESTARTED DUE TO LINMIN FAILURE ...
>
> .... RESTARTED DUE TO LINMIN FAILURE ...
>
>
> FINAL RESULTS
>
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 68 -8.3706E+09 8.5391E+04 8.3864E+06 C5 1446
>
> BOND = 12554.0865 ANGLE = 2862.0342 DIHED = 952.0557
> VDWAALS = ************* EEL = -64213.1344 HBOND = 0.0000
> 1-4 VDW = ************* 1-4 EEL = 4191.5380 RESTRAINT = 3088.1696
> EAMBER = *************
>
> ***** REPEATED LINMIN FAILURE *****
>
>
>
>
> Thank you very much in advance and kind regards,
> Zhonghua
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Mar 08 2019 - 02:30:03 PST
Custom Search