AMBER Archive (2008)Subject: Re: AMBER: dihedral energy calculation inconsistent with gaussian's result
From: Jerome.GOLEBIOWSKI_at_unice.fr
Date: Tue Jun 24 2008  03:01:35 CDT
Hi,
the dihedral parameters are not the only term that is taken into account
for the total energy calculation. Indeed, the 14 terms are not zero for
such a system. The dihedral parameters are just a decomposition of the
remaining energy term once you have taken into account the coulomb and
the vdW interactions in your molecule.
If you want to check the energy barrier, you have to perform a
calculation with all the forcefield terms (especially 14 ones in your
case).
Hope this helps.
Jerome
On Tue, 20080624 at 10:43 +0800, Q733 wrote:
> Dear Amber users:
>
> I constructed a forcefield for an organic molecule containing 81 atoms
> using the Amber 2003 united atom forcefield, and checked the energy
> difference of XC2C2X
> dihedral between 0 and 180 degree .
>
> There are totally 3 sets of dihedral parameters:
> X C2C2X 1 1.60 0.0 3. Yang et
> al, 2005
> X C2C2X 1 0.60 180.0 2. Yang et
> al, 2005
> X C2C2X 1 1.00 180.0 1. Yang et
> al, 2005
>
> I added all of them according to the formula: Etors=(PK/IDIVF)*(1
> +cos(PN*phiphase)) and calculated the dihedral energy at 0 degree and
> 180 degree, which is 3.2 and 2 KJ respectively. The difference is only
> 1.2 kcal/mol, while in Amber8 manual on page 257 it is mentioned that
> the general energy barrier of H3CTCTH3 is about 3 kcal/mol, I can
> not understand why there is such a huge difference.
>
> I also calculated the energy difference of 2methylpetane with the
> dihedral of C2C2 at 0 degree and 180 degree using Gaussian03, and
> kept other parts of the molecule optimized and fixed, the result is
> about 14.76 kcal/mol, but the total energy difference of
> 2methylpentane calculated using 0 step of energy minimization using
> amber is about 9.76 kcal/mol. The energy does not match well.
>
> Is there anything wrong with my understanding or calculation ?
>
> Thanks very much in advance.
>
> Shanshan Qin
>
>

Jerome Golebiowski, PhD
LCMBA, team Chemometrics and Molecular Modeling
University of Nice, parc Valrose
06108 Nice Cedex2 France
tel : +33 (0)4 92 07 61 03
http://www.unice.fr/lcmba
http://www.unice.fr/lcmba/golebiowski

The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo_at_scripps.edu
