Re: [AMBER] drms comparison in sander during minimization?

From: David A Case <david.case.rutgers.edu>
Date: Thu, 22 Oct 2020 08:42:27 -0400

On Wed, Oct 21, 2020, Kenneth Huang wrote:
>
>A quick question- when we set drms=1E-4, what exactly is sander (or
>pmemd) comparing the value to in the output? Ie, given the final reported
>step from sander-
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 32599 -1.9473E+02 5.9628E-05 1.6881E-04 C2 20
>
> BOND = 0.5890 ANGLE = 4.4009 DIHED = 7.0049
> VDWAALS = -2.7607 EEL = 87.7876 EGB = -44.1792
> 1-4 VDW = 5.4083 1-4 EEL = -252.9773 RESTRAINT = 0.0000
>
>I'd been assuming it's the RMS value, but I realized the units don't
>quite make sense since drms is kcal/mol*A, whereas RMS should just be A
>unless it's not the strict Cartesian RMS that I've been thinking it was?

RMS means "root mean sqaure", and it doesn't per se specify what one is
taking the root mean square *of*.

The "RMS" value pointed to above is the root mean square of the elements
of the gradient, and indeed has units of kcal/mol*A. It should probably
have been labeled "DRMS" (derivative RMS) both for clarity and to match
the input variable drms.

....dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Oct 22 2020 - 06:00:03 PDT
Custom Search