HI,
I have been learning Tutorial on Dr. Case's thermodynamic
integration.
When I calculate charging free energy in Vacuum, in my md print outs
EELEC and 1-4 EEL terms are zero, which to me makes sense since I have
set the charges to zero.
But when I simulate the system in solvent, the EELEC and 1-4 EEL terms
are no longer zero (even though I have set the charges of the atoms on
toluene to zero). While I can understand partly why EELEC is non-zero
(this is due to the charge to charge interactions of water molecules), I
do not know why 1-4 EEL is NOT zero. 1-4EEL is the interaction
separated by 3 bonds and all charges on toluene are zero. I am
wondering where else this 1-4 EEL could have come from. Or could I be
making mistakes in preparing the parm file?
Thanks,
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Aug 06 2006 - 06:07:11 PDT
