<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html;charset=GB2312" http-equiv="Content-Type">
  <title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
zjim wrote:
<blockquote cite="mid42B951B3.00000B.17584@m238.163.com" type="cite">dear
all:<br>
I wanted to change the opls torsion parameters into Ryckaert-Bellemans
parameters in order to apply them also in gromacs. According to the
manu, it can be done as follow:<br>
C0 = V0 + V2 + 1/2 (V1 + V3)<br>
C1 = 1/2(3V3 &#8722; V1)<br>
C2 = &#8722;V2<br>
C3 = &#8722;2V3<br>
C4 = 0<br>
C5 = 0<br>
yet, the reference i get doesn't have the value of coefficient V0, so
it was unable to achieve the change.So, my question how can i get the
value of V0 ? Is it a determinate value or not.does somebody know this,
please give some advice.<br>
<!-- footer --><br>
  <font style="font-size: 14.8px;">
<!--广告footer 开始--><!--广告footer 结束--> <br>
  </font></blockquote>
Hi Zjim,<br>
<br>
if you didn't find any value for V0, there is a high probability that
it should be put to 0. For instance, in the OPLS AA force field, V0 is
0.<br>
I would like to advise you to simply plot your torsional potential
function respect to the torsion angle as it is done un the manual and
check the shape of the curve. You should find minima and maxima at the
relevant conformations of your torsional angle system. Then you will be
sure that V0 is really 0.<br>
I hope this will help.<br>
<br>
Frederic, Tarragona Univ, Spain.<br>
<br>
</body>
</html>