On Wed, Aug 7, 2013 at 12:49 PM, John Travers <jtravers70.yahoo.com> wrote:
> Thanks! I made the modifications to output nuclear cartesian forces in
> runmd.F90 and nmode.F90. For the force, the outputs from these two
> subroutines are slightly different. For example, I made an test on an
> isolated benzene molecule, the forces (for the 1st C) from runmd and nmode
> are " -0.1505921505018395E+01 -0.2234835466490895E+01
> 0.4347184130879028E+00" and " -0.1505920513395167E+01
> -0.2234834582828459E+01 0.4347184163951585E+00".
>
> Does anyone know the reason for this?
>
Computers ;). The force field is implemented differently for nmode and
sander, meaning that the order of operations will be different and errors
will accumulate -- especially in quantities like the forces where you're
adding contributions from a bunch of different force field terms.
Also, nmode and sander compute torsion terms fundamentally differently.
While both approaches are mathematically equivalent, the finite numerical
precision consequences are more ambiguous. These numbers appear to agree
quite well.
HTH,
Jason
--
Jason M. Swails
BioMaPS,
Rutgers University
Postdoctoral Researcher
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Aug 07 2013 - 12:00:04 PDT
