I did not develop the force field of my molecule by myself. In stead I use the literature force field parameters. I guess not all the force field need the [pairs]. As long as I can get reasonable experimental observables, I can use such force field, right?<br>

<br><div class="gmail_quote">On Tue, Jan 26, 2010 at 6:50 PM, Mark Abraham <span dir="ltr">&lt;<a href="mailto:mark.abraham@anu.edu.au">mark.abraham@anu.edu.au</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 class="im"><br>
<br>
----- Original Message -----<br>
From: Yanmei Song &lt;<a href="mailto:ysong30@asu.edu">ysong30@asu.edu</a>&gt;<br>
Date: Wednesday, January 27, 2010 10:58<br>
Subject: Re: [gmx-users] Potential Energy<br>
To: <a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>, Discussion list for GROMACS users &lt;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>
<br>
&gt; Dear Justin:<br>
&gt;<br>
&gt; I finally figure out that the problem comes from the Coulomb 1-4 interactions defined by [pairs] part. When I remove this parts, the energy minimization ends with a big negative number. Thanks for your help.<br>
<br>
</div>That&#39;s not a solution, that&#39;s a random number generator. Whether or not you need such pairs is force-field dependent. Ad hoc changes are ridiculous. Messing with force fields is an expert topic, and parameterizing a new molecule is difficult and risky.<br>


<font color="#888888"><br>
Mark<br>
</font><div><div></div><div class="h5">--<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/search" target="_blank">http://www.gromacs.org/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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Yanmei Song<br>Ph.D. Candidate<br>Department of Chemical Engineering<br>Arizona State University<br>