<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto">Hi Yu,<div><br></div><div>Yes. The reason for this is that tabulated potentials require you to provide the form of the repulsion and dispersion parts of the interaction separately. Otherwise it would never be possible to tabulated e.g. a switched or shifted interaction.</div><div><br></div><div>At first it sounds like a solution to rather tabulated</div><div><br></div><div>x^-12 - x^-6</div><div><br></div><div>And then use x=sigma/r, but then all cutoffs and shifts in the table would be applied to x rather than r, which would violate a bunch of physics.</div><div><br></div><div>So, in order to apply dispersion and repulsion separately we must have separate coefficients to apply for each interaction, and that way of specifying combination rules corresponds to normal c6/c12 with geometric combinations.</div><div><br></div><div>However, as Berk wrote there's a trick you can use: </div><div><br></div><div>If you read the topology section of the manual (in detail, sorry, it's somewhat condensed, but all the info is there), just as we have a section for pairrtypes (1,4 parameters) you can create a separate section for nonbonded parameters and explicitly specify c6 and c12 for each pair of atomtypes (have a look at the old Gromos force fields for an example).</div><div><br></div><div>You will have to translate from sigma/epsilon with a small script, and the section might be large if you have a ton of atomtypes, but it will work.</div><div><br></div><div>We might be able to automate this better when enabling the verlet tables, but that's still a bit away ;-)</div><div><br></div><div>Cheers,</div><div><br></div><div>Erik</div><div><br></div><div><br><div dir="ltr">--<div>Erik Lindahl <erik.lindahl@scilifelab.se></div><div>Professor of Biophysics</div><div>Science for Life Laboratory</div><div>Stockholm University & KTH</div><div>Office (SciLifeLab): +46 8 524 81567</div><div>Cell (Sweden): +46 73 4618050 </div><div>Cell (US): 1 267 307 8746</div><div><br></div></div><div dir="ltr"><br><blockquote type="cite">On Dec 1, 2019, at 02:42, Du, Yu <duyu@sioc.ac.cn> wrote:<br><br></blockquote></div><blockquote type="cite"><div dir="ltr">Hi Prof. Hess, <div><br></div><div>Thanks for your reply.</div><div><br></div><div>It is because user table requires comb-rule=1 that I have to change the comb-rule of AMBER and CHARMM </div><div>from 2 to 1. I realy have no special need for comb-rule!=1. So could you give some hints of "You can </div><div>list parameters for all combinations of atoms types you need/want explicitly."? I don't know the way </div><div>to explicitly list the combined values for all atom pairs.</div><div><br></div><div>Yes, I will take your caveat of larger errors.<br><br>Cheers, </div><div><br></div><div>Yu</div><div><br><blockquote name="replyContent" class="ReferenceQuote" style="padding-left:5px;margin-left:5px;border-left:#b6b6b6 2px solid;margin-right:0">-----Original Messages-----<br>
<b>From:</b><span id="rc_from">"Berk Hess" <hess@kth.se></span><br>
<b>Sent Time:</b><span id="rc_senttime">2019-12-01 02:26:54 (Sunday)</span><br>
<b>To:</b> gmx-developers@gromacs.org<br>
<b>Cc:</b> <br>
<b>Subject:</b> Re: [gmx-developers] comb-rule=2 support for forcetable<br><br><div dir="auto"><div>Hi,<div dir="auto"><br></div><div dir="auto">But you don't need to change the combination rule. You can list parameters for all combinations of atoms types you need/want explicitly.</div><div dir="auto"><br></div><div dir="auto">Note that changing the combination rule for a force field likely results in much larger errors than the differences between different force fields.</div><div dir="auto"><br></div><div dir="auto">Cheers,</div><div dir="auto"><br></div><div dir="auto">Berk</div><br><div class="gmail_extra"><br><div class="gmail_quote">On Nov 30, 2019 14:36, "Du, Yu" <duyu@sioc.ac.cn> wrote:<br type="attribution"><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Prof. Lindahl, <div><br></div><div>Thanks for your reply.</div><div><br></div><div>If user table is used e.g. coulombtype=PME-user and vdwtype=user, I think that the only difference </div><div>of geometric (comb-rule=1) and arithmetic (comb-rule=2) average of sigma that can't be eliminated </div><div>by editing the value in force field.</div><div><br></div><div>comb-rule=1 matches GROMOS force field which is deprecated for its physically incorrect scheme.</div><div>comb-rule=2 matches AMBER and CHARMM force fields which I want to use in the a compatible comb-rule</div><div>under circumstances of user table.</div><div><br></div><div>As you suggested, I indeed have copied and modified AMBER and CHARMM ff with a new name :), but in </div><div>comb-rule=1 for user specified nb potentials between groups. This is why I'm asking for a comb-rule=2</div><div>support.</div><div><br></div><div>By the way, I notice that forcetable in Verlet scheme has been under development. I'm excited about</div><div>it and hope it can move to GPU in the future.</div><div><br></div><div>Regards,</div><div>Yu<br><br><blockquote style="padding-left:5px;margin-left:5px;border-left:#b6b6b6 2px solid;margin-right:0">-----Original Messages-----<br>
<b>From:</b>"Erik Lindahl" <erik.lindahl@gmail.com><br>
<b>Sent Time:</b>2019-11-30 17:08:41 (Saturday)<br>
<b>To:</b> gmx-developers@gromacs.org<br>
<b>Cc:</b> <br>
<b>Subject:</b> Re: [gmx-developers] comb-rule=2 support for forcetable<br><br><div dir="ltr">Hi Yu,<div><br></div><div>That won't really require any changes in the code.</div><div><br></div><div>First, the reason we can't "support" it is that force field Lennard-Jones parameters have been determined for a specific set of combination rules and cutoffs, and when you change the combination rules you will change the result of the simulation. How good or bad that is depends on what you are doing, but it's something you'd have to determine and be ready to defend when publishing.</div><div><br></div><div>Having said that, if you absolutely want to change the setting for a particular force field, you can just edit the value in the force field file. I would strongly advise making a copy of the entire force field directory and changing its name to reflect this, so you don't forget and have an altered installed version of the force field ruin other simulations later.</div><div><br></div><div>Cheers,</div><div><br></div><div>Erik</div></div><br><div class="elided-text"><div dir="ltr">On Sat, Nov 30, 2019 at 3:11 AM Du, Yu <<a href="mailto:duyu@sioc.ac.cn" target="_blank">duyu@sioc.ac.cn</a>> wrote:<br></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb( 204 , 204 , 204 );padding-left:1ex">Dear GMX Developers, <div><br></div><div>The force table function (src/gromacs/tables/) in GROMACS is a magic, it helps to achieve a lot of work.</div><div>Could you please support it in comb-rule=2 for AMBER and CHARMM FF? </div><div><br></div><div>Sincerely, </div><div>Yu</div>-- <br>
Gromacs Developers mailing list<br>
<br>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List</a> before posting!<br>
<br>
* Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
* For (un)subscribe requests visit<br>
<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers</a> or send a mail to <a href="mailto:gmx-developers-request@gromacs.org" target="_blank">gmx-developers-request@gromacs.org</a>.</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Erik Lindahl <<a href="mailto:erik.lindahl@dbb.su.se" target="_blank">erik.lindahl@dbb.su.se</a>></div><div>Professor of Biophysics, Dept. Biochemistry & Biophysics, Stockholm University</div><div>Science for Life Laboratory, Box 1031, 17121 Solna, Sweden</div></div></div></div></div></div></div></div></div></blockquote></div></blockquote></div><br></div></div></div></blockquote><br><span>--<br>Du, Yu<br>PhD Student,<br>Shanghai Institute of Organic Chemistry<div>345 Ling Ling Rd., Shanghai, China. </div><div>Zip: 200032, <span style="font-size: 14px;">Tel: (86) 021 5492 5275</span></div></span></div><span>-- </span><br><span>Gromacs Developers mailing list</span><br><span></span><br><span>* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List before posting!</span><br><span></span><br><span>* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists</span><br><span></span><br><span>* For (un)subscribe requests visit</span><br><span>https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers or send a mail to gmx-developers-request@gromacs.org.</span></div></blockquote></div></body></html>