Dear all,
I went thru the valine to alanine Gibbs tutorial in two variations in the
charges for the resulting atoms; first I used the
0.105 and zeros according to the tutorial, second, I used charge for alanine
as in lib files.
Can you tell me why in the tutorial is this - to me arbitrary - charge
value? Is the standard procedure to change the
perturbed-atom (say C-to-H) charges to their actual values or zeroes?
I got two results both of which differ from those in the tutorial (I put
them below). Moreover, in the 2nd case (which I think
is the correct one) even the sign is reversed! Could you tell me what;s
wrong?
Thank you.
All the best,
Martin Lepsik
1st case:
Current Lambda = 0.000000
Last F.E. update: Lambda = 0.000000 Step = 11000 Method = F.E.P.
Accumulated "forward" quantities (Regular)
Lam+d_lam = 0.100000 F_energy = 5.59082
ELEC = 24.769 NONB = 1.324 14NB = 2.887
14EL = -25.201 BADH = -0.148
Accumulated "reverse" quantities (Regular)
Lam-d_lam = 0.000000 F_energy = -5.10742
ELEC = -25.078 NONB = -4.857 14NB = -3.274
14EL = 25.592 BADH = 0.120
2nd case:
Current Lambda = 0.000000
Last F.E. update: Lambda = 0.000000 Step = 11000 Method = F.E.P.
Accumulated "forward" quantities (Regular)
Lam+d_lam = 0.100000 F_energy = -3.58055
ELEC = 23.330 NONB = 0.275 14NB = 3.247
14EL = -32.252 BADH = -0.127
Accumulated "reverse" quantities (Regular)
Lam-d_lam = 0.000000 F_energy = 4.76917
ELEC = -23.653 NONB = -2.765 14NB = -3.258
14EL = 32.580 BADH = 0.115
----------------------------------------------------------------------------
----------------------------------------------------------------------------
----
Martin Lepsik, PhD student Phone: +420/2/20183-540, Fax:
+420/2/20183-292
Dept. of Theor Chem & Center for Complex molecular Systems and Biomolecules
Institute of Organic Chemistry and Biochemistry (IOCB)
Flemingovo nam 2,
Czech Academy of Sciences,
166 10, Prague 6,
Czech Rep.
URL: www.uochb.cas.cz/~teochem
Received on Mon Dec 16 2002 - 05:20:26 PST
