AMBER Archive (2004)Subject: Re: AMBER: Amber Transition State Output
From: Johnson Agbo (agboj_at_chem.unr.edu)
Date: Wed Jun 09 2004 - 19:33:58 CDT
I am using nmode to find a transition state structure(coordinates) and the
energy barrier between my input structure or mininum structure and that of
the transition state. My input structure has 24.0 kcal/mol of energy. The
output file I am looking at has as its last part for the transition state
search as :
step = 56
F = -0.488232E+01 GRDMAX = 0.662498E-03 GNORM = 0.132406E-03
E-NONB E-ELE E-HBOND E-BOND
-0.20634E+01 -0.28338E+02 0.00000E+00
0.57860E+00
E-ANGLE E-DIHED E-NB14 E-EEL14
0.27568E+01 0.45998E+01 0.57571E+01
0.11826E+02
E-POL E-3BOD
0.00000E+00 0.00000E+00
Found transition state:
***************** ***************** *****************
step = 56
F = -0.488232E+01 GRDMAX = 0.662498E-03 GNORM = 0.132406E-03
E-NONB E-ELE E-HBOND E-BOND
-0.20634E+01 -0.28338E+02 0.00000E+00
0.57860E+00
E-ANGLE E-DIHED E-NB14 E-EEL14
0.27568E+01 0.45998E+01 0.57571E+01
0.11826E+02
E-POL E-3BOD
0.00000E+00 0.00000E+00
b1 = -0.07361 b2 = 0.06656 tlamba = -0.17361
rms of step length = 0.0577350278
ene @step k-1 = -4.882317 ene @k = 14.533882
ene @step k-1 = -4.882317 ene @k = 7.543089
ene @step k-1 = -4.882317 ene @k = 3.070845
ene @step k-1 = -4.882317 ene @k = 0.208641
ene @step k-1 = -4.882317 ene @k = -1.623461
ene @step k-1 = -4.882317 ene @k = -2.796260
ene @step k-1 = -4.882317 ene @k = -3.547019
ene @step k-1 = -4.882317 ene @k = -4.027604
ene @step k-1 = -4.882317 ene @k = -4.335235
ene @step k-1 = -4.882317 ene @k = -4.532149
ene @step k-1 = -4.882317 ene @k = -4.658191
Coordinate origin is now the molecular center of mass
b1 = 0.00151 b2 = 0.02023 tlamba = -0.09849
rms of step length = 0.4152882440
ene @step k-1 = -4.658191 ene @k = 58.765594
ene @step k-1 = -4.658191 ene @k = 34.889225
ene @step k-1 = -4.658191 ene @k = 19.807583
ene @step k-1 = -4.658191 ene @k = 10.320876
ene @step k-1 = -4.658191 ene @k = 4.384055
ene @step k-1 = -4.658191 ene @k = 0.693016
ene @step k-1 = -4.658191 ene @k = -1.582165
ene @step k-1 = -4.658191 ene @k = -2.968536
ene @step k-1 = -4.658191 ene @k = -3.799999
ene @step k-1 = -4.658191 ene @k = -4.287485
ene @step k-1 = -4.658191 ene @k = -4.563759
***************** ***************** *****************
I hope this helps.
Thanks.
On Wed, 9 Jun 2004, David A. Case wrote:
> On Wed, Jun 09, 2004, Johnson Agbo wrote:
>
>> How is the output file for the transition state interpreted. I cannot
>> locate the barrier height from the results of the output file.
>
> Can you be more specific about what sort of output you are looking at? Is
> this from nmode? Maybe you could post a bit of the output that you find
> hard to interpret, plus a little bit about what you are doing.
>
> ...dac
>
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber_at_scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
>
--
+++++++++++++++++++++++++++++++++++
+ Residential Address: +
+ 34 West 9th Steet, Apt. 8
Reno, NV. 89503 +
+ USA +
+ Tel. +1-775-322-5082 +
+ +
+ Office Address:
+ Dept. of Chemistry/216 +
+ Univ. of Nevada, Reno +
+ Reno, NV 89557-0020 +
+ USA +
+ Tel. +1-775-784-8002 +
+++++++++++++++++++++++++++++++++++
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|