AMBER: Questions about pert.charge and atom type change.

From: Chunhu Tan <tanc.uci.edu>
Date: Wed, 30 Jun 2004 17:20:38 -0700

 Hi, all,

    I am trying to do TI with sander, but I am still confused with the
pert.charge
in xleap. My question is, when a H (with the charge 0.0791) atom becomes a
dummy, if
its pert.charge is 0-0.0791=-0.0791? If this is true, why the dummy H atoms
on
this page have no pert.charge?
http://web.njit.edu/all_topics/Prog_Lang_Docs/html/amber7/gibbs_box/HTML/xleap_edit_selected.html
Also on this page, the first CT's charge, in VAL, is 0.2985, in ALA,
is -0.1825, so its pert.charge
should be -0.1825-0.2985=0.4810, not the 0.105, am I right?

   Another question. When I perturbed CH3-OH to CH3-SH, if both the O
atom(type OH) and the H
atom attached to it (type HO) should be changed to type SH and HS
respectively? I tried to only
change the O, but when I saveamberparmpert, LEap told me :
     No bond parameter for: SH-HO
     Can't find angle parameter: CT-SH-HO
When I also changed the H, it is OK, but I am not sure it is right, and what
about the pert.charge of H?

Thanks in advance!!




-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Thu Jul 01 2004 - 01:53:00 PDT
Custom Search