<html>
<head>
</head>
<body>
Hi Dallas,<br>
<br>
A lot of MD force fields have this problems with aliphatic dihedrals. I would
guess that OPLS is slightly better than Gromos96, which in turn should be
better than Gromacs. (The "gromacs" force field is really Gromos87 with some
modifications mentioned in the manual - we don't try to develop our own force
field, it is too<br>
much religion in it :-)<br>
<br>
If you want to stick to united atom I would recommend either Ryckaert-Bellemans
potentials or The Kuwajima alternative. The R-B potential is already defined<br>
in Gromacs; just change the dihedral type to "3" and make sure you REMOVE
the 1,4 interactions since these are included in this functional form.<br>
<br>
Cheers,<br>
<br>
Erik<br>
<br>
<br>
Dallas Warren wrote:<br>
<blockquote type="cite" cite="mid:5.1.0.14.2.20021115131802.00b2a728@pop.ozhosting.com">
 I was just wondering if someone here knows the reason for the value used 
for the CH2-CH2 dihedral in the GROMACS force field?<br>
  <br>
 I have been doing some simulations containing bile salts, which have a short
hydrocarbon chain attached to essentially cholesterol.&nbsp; The chain was too
straight, so I started to look into why this would be the case.&nbsp; A simulation
of 100 decane molecules showed chains that are very straight with only one,
may be two gauche confirmations per chain.&nbsp; Something isn't looking right
here ;-)<br>
  <br>
 So, then looked for a paper on the energies of a CH2-CH2 dihedral, pulled 
Allinger et al. J. Chem. Phys., 106(12) (1997) which is "The torsional conformations
of butane: definitive energetics from ab initio methods".&nbsp; The difference
in energy between the anti and gauche conformations is 0.62 kcal/mol.&nbsp; I
performed an energy minimisation of butane in the anti and gauche confirmations
and found that the energy difference using the GROMACS forcefield is 1.58
kcal/mol!&nbsp; The energy for the dihedral is over twice of what it really should
be, no wonder the chains are straight!<br>
  <br>
 How was the value decided on?&nbsp; It just seems to me that is a huge difference,
even if the force field isn't meant for hydrocarbons, proteins/amino acids
still contain chains of CH2-CH2 and those sections are going to be too rigid.<br>
  <br>
 Catch ya,<br>
  <x-sigsep></x-sigsep>
  <p><b>Dr. Dallas Warren<br>
  </b><i>Research Fellow<br>
  </i>Department of Pharmaceutical Biology and Pharmacology<br>
 Victorian College of Pharmacy<br>
 Monash University<br>
 381 Royal Parade<br>
 Parkville VIC 3010<br>
  <font color="#0000ff"><u><a class="moz-txt-link-abbreviated" href="mailto:dallas.warren@vcp.monash.edu.au">dallas.warren@vcp.monash.edu.au</a><br>
  </u></font>+61 3 9903 9076</p>
  </blockquote>
  <br>
  </body>
  </html>