On Tue, Aug 01, 2006, Ilyas Yildirim wrote:
>
> Is there any way to set the bondii radii to zero or any value for a
> specific atom type using xleap?
No.
> Another question is the following: When the bondii radii for a specific
> type is set to zero, and igb=2 is used, I am getting NaN results in the
> output file (when bondii radii is not zero - even when it is very close
> to zero - the result is ok). I was wondering why this is creating a
> singularity in the system when implicit solvent is used.
You can't use a zero radius for any atom with a non-zero partial charge: that
would create real problems with the electrostatic energies. Even for atoms
with zero partial charge, using a zero dielectric radius is probably dangerous
(as you have found), since the code was certainly never carefully checked for
this possibility.
...dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Aug 02 2006 - 06:07:13 PDT
