Dear Amber community,
I have a question about the way that induced dipole is
implemented in the sander module of AMBER8.
When I looked in
$AMBERHOME/src/sander/prn_dipoles.f
I saw that
Indiptotal = sqrt((inddipx^2)+(inddipy^2)+(inddipz^2))
However, I have not seen any definition for
inddip(x,y,z) array in prn_dipoles.f .
Could someone please formulate inddip used in sander?
Isn't it proportional to field strength
(E*=E.COS(teta) where E is an electric field)and
polarizibility (alpha) [inddip = alpha.E*]?
If this is correct then what would be the value of E
for an atom group of a molecule immersed in POL3BOX.
One could also apply Debye's formula for a dilute
solution in a non-polar solvent to determine induced
dipole and polarizibility as a function of temperature
and other variables such as dielectric constant,
viscosity, molecular weight, etc., but I am
not exactly sure if inddip in sander works like that.
best regards,
Jenk.
__________________________________________________
Do You Yahoo!?
Tired of spam? Yahoo! Mail has the best spam protection around
http://www.rosswalker.co.uk/adsense_alternatives/
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Apr 22 2005 - 18:53:00 PDT
