Re: [AMBER] bond length energy

From: Mark Williamson <mjw.sdsc.edu>
Date: Wed, 17 Jun 2009 18:29:09 +0100

Thomas Lake wrote:
> Hi,
>
> Is there a function for calculating the bond length energy between a
> pair of atoms in a molecule? In my case I want to calculate the bond
> length energy of a disulfide (S-S) bond, as a function of time. I can
> calculate the bond distance using ptraj.
>
> Thanks
>
> Thomas Lake
> Imperial College London
>


Hi Thomas,

If you want to know the energy contribution from *only* the two body
bonding term, and know the length of this bond as a function of time and
its parameters, then you have all you need to calculate its potential
energy:


amber11/dat/leap/parm$ grep -e ^"S -S" parm99.dat
S -S 166.0 2.038 JCC,7,(1986),230; CYX (SCHERAGA)

$ cat /tmp/distance_output_from_ptraj.dat
2.2
2.1
2.3
2.4
etc....

$ awk < /tmp/distance_output_from_ptraj.dat '{print 166.0*($1-2.038)**2 }'
4.3565
0.638104
11.3949
21.7533
etc...


regards,

Mark


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Jul 06 2009 - 09:55:36 PDT
Custom Search