On Wed, Apr 21, 2004, david.evans.ulsop.ac.uk wrote:
>
> I am trying to run an MD simulation where the non-bonded interactions of a
> small group of atoms are set to zero...
>
> The system however gives NaN halfway through an MD simulation when two atoms
> come close to each other.
I don't think you can avoid this without changing the code. Amber still
computes r**-12 for every pair of atoms, then multiplies by the appropriate
LJ parameters. What is probably required is to test _first_ for zero
parameters, and skip the calculation if appropriate. Another approach would
be to hack the excluded atom list, so that the dummy atoms were excluded from
all other atoms. Yet another idea would be to actually remove the atoms, and
add an artificial bond or NMR-like restraint to "tether" the parts of the
molecule together.
....good luck...dac
--
==================================================================
David A. Case | e-mail: case.scripps.edu
Dept. of Molecular Biology, TPC15 | fax: +1-858-784-8896
The Scripps Research Institute | phone: +1-858-784-9768
10550 N. Torrey Pines Rd. | home page:
La Jolla CA 92037 USA | http://www.scripps.edu/case
==================================================================
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Apr 21 2004 - 17:53:00 PDT
