Hi guys.<br><br>Well, I'm not "exactly" trying t reopen an old discussion. I'm trying to see a different aspect from it, as the subject indicates.<br><br>Since the message reproduced below is the last of the "fudges" questions, I'm worried about the fact of the fudgeLJ meaning nothing for the simulation.
<br><br>First, is it true, and there is no way around to get a working fudge LJ of 0.5 AND a fudgeQQ of 0.8333 in gromacs simulations?<br><br>Having this in mind, shouldn't it mean that the ffamber, the amber port to gromacs, have a "problem" in that point? Specially because it uses, up to the point I could see, the same "configurations" in the "default" section?
<br><br>I have deep concerns on this because I would like to move on to gromacs to make the next simulations I'm interested in. My interests in gromacs range from its speed to simplicity and range of analysis programs it already have (and, of course, the cost of the program). But this specific point is a big reason of concern to me, and I would like to see it clarified.
<br><br>Thanks a lot for any help in advance.<br><br>Sincerally yours,<br><br>Jones<br><br><div><span class="gmail_quote">On 12/7/06, <b class="gmail_sendername"><a href="mailto:chris.neale@utoronto.ca">chris.neale@utoronto.ca
</a></b> <<a href="mailto:chris.neale@utoronto.ca">chris.neale@utoronto.ca</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
> I am puzzled about this. In your opinion, when I use genpairs=no, fudgeLJ=<br>> 0.5 and fudgeQQ=0.8333, are these scaling factors, 0.5 and 0.8333, useful?<br>> On the other hand, when using genpairs=yes, how to set the scaling factor of
<br>> fudgeLJ and fudgeQQ?$<br><br>fudgeQQ is always used no matter what the value of genparis is.<br><br>(exerpt from<br><a href="http://www.gromacs.org/pipermail/gmx-users/2006-September/023743.html">http://www.gromacs.org/pipermail/gmx-users/2006-September/023743.html
</a>)<br>Each force-field has its own rules (e.g. gen-pairs and FudgeLJ/QQ), but these<br>apply to the information outlined above. For example, gen-pairs does NOT mean<br>"generate a [ pairs ] section for the molecule." Instead, it means "If LJ-14
<br>epsilon and sigma are not present in a [ pairs ] section entry, and<br>that type of<br>interaction is not explicitly formulated in [ pairtypes ], then it is<br>permissible to use the regular non-bonded parameters, and in that case scale
<br>them by FudgeLJ."<br><br>Therefore your settings indicate that coulombic 1-4 interactions will<br>be scaled by 0.8333 and the pairs must be taken directly from the [<br>pairs ] or [ pairtypes ] section (and they will NOT be scaled by
<br>fudgeLJ). In your case fudgeLJ does not matter (you could change the<br>value and it would not affect your simulation). However, I have always<br>hoped that it is set to 0.5 to indicate that this is what the<br>forcefield developers have done for you and included it in [ pairtypes
<br>] so that it's kind of a reference value for your piece of mind.<br><br>I was assuming that the previous messages were from an Amber<br>forcefield. Since Amber uses (to my knowledge) [ pairtypes ] instead<br>of genpairs=yes, I was confused by the combination of genpairs=yes and
<br>fudgeQQ=0.8333.<br><br><br>_______________________________________________<br>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a href="http://www.gromacs.org/mailman/listinfo/gmx-users">
http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please don'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't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php">http://www.gromacs.org/mailing_lists/users.php</a><br></blockquote></div><br>