[gmx-users] On CHARMM angle bending potential

[Mr Bernard Ramos]

Hi David!
Thanks for the reply. I may need to check the details myself. Where in the Gromacs manual can I find them? Are there references I can check. I am only interested in knowing that the force field does not blow up. It would be a relief to know it does not diverege. Thanks 

Bernard 

    On Monday, July 4, 2016 9:28 PM, David van der Spoel <spoel at xray.bmc.uu.se> wrote:
 

 On 04/07/16 12:11, Mr Bernard Ramos wrote:
> Hi all!
>
> The potential for the angle-vibration potential in CHARMM is V = k(a - a0)^2. Its corresponding force field is proportional to 1/sin(a) which blows up when a approaches \pi or when the molecule is in linear arrangement. I am wondering how Gromacs calculates the angle-bending potential of a linear molecule such as carbon dioxide when it the force vector is singular at a equal to \pi?
>
> Thank you.
>
> Bernard
>
It depends on the force field. Gaff sets the reference angle to 178 or 
so. An alternative is to define the C as a virtual site, which is 
probably easiest to implement in most force fields, but only practical 
if the angle fluctuations are irrelevant for the problem at hand. My 
group are working on more systematic way of treating angle potentials, 
and one such linear angle is implemented in gromacs already, however it 
has not been parameterized rigorously and published.

-- 
David van der Spoel, Ph.D., Professor of Biology
Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:    +46184714205.
spoel at xray.bmc.uu.se    http://folding.bmc.uu.se
-- 
Gromacs Users mailing list

* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.