Re: AMBER: Help in radii

From: FyD <fyd.q4md-forcefieldtools.org>
Date: Wed, 16 Jan 2008 21:31:07 +0100

Quoting snoze pa <snoze.pa.gmail.com>:

> I am getting following error message after running a small Gaussian
> calculation.
> *
> No NMR shielding tensors so no spin-rotation constants.
> Leave Link 601
> (Enter g03/l602.exe)
> Merz-Kollman atomic radii used.
> Atom Element Radius
>
> [...]
> 14 26 0.00
> GetVDW: no radius for atom 14 atomic number 26.
> Error termination via Lnk1e in g03/l602.exe
> *
> How can I get the GetVDW for atomic no 14 in Gaussian. Do I need to
> mention it in the beginning of the input Gaussian file
> like Pop=(ReadRadii)! Any help. Thanks in advance

In R.E.D., this problem is automatically handled & the input is
automatically generated... You can also modify the source code, and
change the Fe radii if you do not like the "1.8" value set by default.

> What is the actual radii of Fe shall i use? in amber archive i saw somebody
> was suggesting 0.49
> http://archive.ambermd.org/200408/0096.html
> while FyD notebook it is recommended to use larger one like 2.0 etc. The
> program like vmd suggest 2.01 while in gammes it is 1.8 Which one is the
> right one. thanks in advance

Why not doing some testing: Are you sure this radii strongly affects
the Fe charge value ? Did you test that ?

regards, Francois


-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Jan 20 2008 - 06:07:09 PST
Custom Search