Dear Rainer,<br><br>Thankyou very much. I read the abstract of this paper and it says torsional parameters fro alkanes have been updated. I do not have access to full contents of this paper, so could not get it.<br>Thanks again.<br>
<br>Sulatha<br><br><div class="gmail_quote">On Wed, Feb 9, 2011 at 2:27 PM, Rainer Boeckmann <span dir="ltr">&lt;<a href="mailto:rainer.boeckmann@biologie.uni-erlangen.de">rainer.boeckmann@biologie.uni-erlangen.de</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div style="word-wrap: break-word;"><div class="im">Dear Sulatha,<br><br>these are probably the more recent parameters from<br>
<br>Price, M. L. P., D. Ostrovsky, and W. L. Jorgensen. 2001. Gas-phase and liquid-state properties of<br>8esters, nitriles, and nitro compounds with the opls-aa force field. J. Comput. Chem. 22:1340–<br>1352.<br><br>Best<br>
Rainer<br></div><div><div class="im"><div><br></div><div>On Feb 9, 2011, at 6:48 AM, sulatha M. S wrote:</div><br></div><blockquote type="cite"><div><div></div><div class="h5">I looked at the the paper published in 1999, ( Jorgensen et al JACS, 121, 20, 4831, 1999), the aliphatic torsional parameters are the same as those from 1996. here are the values,<br>
 <br><br>                                   V1                    V2                    V3 (kcal/mol)<br>
CT-CT-CT-CT              1.740               -0.157                 0.279<br><br>HC-CT-CT-CT              0.000                0.000                0.366<br><br>HC-CT-CT-HC             0.000                 0.000               0.318<br>

<br><br>Values calculated from the equation given in p.62 of the manual in kJ/mol are <br><br>                                 c0                       c1                   c2                         c3<br>
CT-CT-CT-CT             3.56686            -1.88907            0.65688                -2.33467<br>
<br>
HC-CT-CT-CT              0.66526           1.99577            0.000                     -2.66102<br>
<br>
HC-CT-CT-HC             0.76567            -2.29702           0.000                     -3.06269<br><br>and the values given in ffoplsaabon.itp are:<br><br>CT-CT-CT-CT             2.9288            -1.4644                0.2092               -1.6736<br>



<br>

HC-CT-CT-CT             0.6276           1.88280               0.000                  -2.5104<br>

<br>

HC-CT-CT-HC             0.6276            1.8828               0.000                   -2.5104<br><br>So there is a difference. Which of these is correct ? Any help is highly appreciated.<br><br>Sulatha<br><br><br><br><div class="gmail_quote">

On Tue, Feb 8, 2011 at 7:48 PM, Andrew Paluch <span dir="ltr">&lt;<a href="mailto:apaluch@nd.edu" target="_blank">apaluch@nd.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

There is no error. The alkane dihedral parameters were updated in 1999, and differ from those originally published in 1996.<br><br>Andrew<br><br><div class="gmail_quote"><div><div></div><div>On Tue, Feb 8, 2011 at 5:20 AM, sulatha M. S <span dir="ltr">&lt;<a href="mailto:mssulatha@gmail.com" target="_blank">mssulatha@gmail.com</a>&gt;</span> wrote:<br>


</div></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div><div></div><div>Hi all,<br><pre>Hi</pre><pre> </pre><pre>I&#39;ve converted the OPLS-AA torsional potential parameters for the</pre>


<pre>alkane C-C-C-C (CT-CT-CT-CT in the gromacs parameter file notation),</pre><pre>C-C-C-H (CT-CT-CT-HC), and H-C-C-H (HC-CT-CT-HC) torsions from the</pre>
<pre>OPLS format given in Jorgensen et al, JACS 118, 11225 (1996) </pre><pre>to the Ryckaert-Bellemans format given in ffoplsaabon.itp and found</pre><pre>that the calculated values are different.</pre><pre> </pre><pre>A previous post to the gmx-users mailing list on March 27, 2008,</pre>



<pre>pointed out this issue for the H-C-C-H torsional potential but there was no response to that.</pre><pre> </pre><pre>Does anyone know if there is an error in the ffoplsaabon.itp file? Or</pre><pre>is there a newer set of OPLS-AA parameters?</pre>



<pre> </pre><pre>For the OPLS-AA parameters (in kcal/mol), I used:</pre><pre> </pre><pre>dihedral<span>              </span>V1<span>          </span>V2<span>             </span>V3</pre><pre>C-C-C-C<span>            </span>1.740<span>      </span>-0.157<span>         </span>0.279</pre>



<pre>C-C-C-H<span>             </span>0.0<span>          </span>0.0<span>            </span>0.366</pre><pre>H-C-C-H<span>             </span>0.0<span>          </span>0.0<span>            </span>0.318</pre>
<pre> </pre><pre>from which I calculated the Ryckaert-Bellemans parameters (in kJ/mol) as:</pre><pre> </pre><pre>dihedral<span>              </span>C0<span>    </span><span>      </span>C1<span>             </span>C2<span>           </span>C3</pre>



<pre>C-C-C-C<span>            </span>3.56686<span>   </span>-1.889076<span>   </span>0.65688     <span>  </span>-2.33467</pre><pre>C-C-C-H<span>             </span>0.66526<span>   </span>1.99577<span>     </span>0.0<span>           </span>-2.661024</pre>



<pre>H-C-C-H<span>             </span>0.76567<span>   </span>-2.29702<span>    </span>0.0<span>           </span>-3.06269</pre><pre> </pre><pre>the parameters in the ffoplsaabon.itp file are:</pre>
<pre> </pre><pre>dihedral<span>              </span>C0<span>          </span>C1<span>             </span>C2<span>           </span>C3</pre><pre>C-C-C-C<span>           </span>2.9288<span>      </span>-1.4644<span>      </span>0.2092<span>     </span>-1.6736</pre>



<pre>C-C-C-H<span>           </span>0.6276<span>       </span>1.88280<span>    </span>0.0<span>           </span>-2.5104</pre><pre>H-C-C-H<span>           </span>0.6276<span>       </span>1.8828<span>      </span>0.0<span>           </span>-2.5104</pre>



<br><br>Thankyou for any clarification.<br><font color="#888888"><br>Sulatha<br>
</font><br></div></div>--<br>
gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br></blockquote></div><br>
<br>--<br>
gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br></blockquote></div><br></div></div>
-- <br><div class="im">gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>Please don&#39;t post (un)subscribe requests to the list. Use the <br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a></div>
</blockquote></div><br><div class="im"><div>
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; font-size: medium;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<div>__________________________________________________________<br>Prof. Dr. Rainer Böckmann<br>Computational Biology<br>Universität Erlangen-Nürnberg<br>Staudtstr. 5<br>91058 Erlangen<br></div><div>Phone:        ++49 +9131 85-25409<span style="white-space: pre-wrap;">        </span>          FAX: ++49 +9131 85-25410<br>
E-Mail: <a href="mailto:rainer@bioinformatik.uni-saarland.de" target="_blank">rainer.boeckmann@biologie.uni-erlangen.de</a><br><a href="http://www.biotechnik.nat.uni-erlangen.de/research/boeckmann/" target="_blank">http://www.biotechnik.nat.uni-erlangen.de/research/boeckmann/</a></div>
<div>___________________________________________________________<br><br></div></div></span></div></span></div></span></div></span></div></span></div></span></div></span></span>
</div>
<br></div></div><br>--<br>
gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don&#39;t post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br></blockquote></div><br>