Re: [AMBER] NaN VDWAALS Error During Minimization

From: David A Case <david.case.rutgers.edu>
Date: Tue, 5 Jan 2021 22:44:39 -0500

On Tue, Jan 05, 2021, Jatin Kashyap wrote:
>
>I am getting an error[1] upon performing minimization of drug-protein
>complex by using the AMBER tutorial[2].

You have two atoms on top of one another in your starting structure:
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 1 NaN NaN 9.2217E+09 O1 18
>
> BOND = 13750.1241 ANGLE = 10976.5308 DIHED = 2934.9414
> VDWAALS = NaN EEL = Infinity EGB = 60228.2741
> 1-4 VDW = 381037.7111 1-4 EEL = 14860.5376 RESTRAINT = 0.0000

Dave Cerutti pointed you to mdgx, but (for me) it's much simpler just
to use the cpptraj "check" command: it will print out short non-bonded
contacts and bad bonds (the two most common problems with initial
structures.)

cpptraj <prmtop file>
trajin <inpcrd file>
check
go

....dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Jan 05 2021 - 20:00:02 PST
Custom Search