Hi<div><br></div><div>I think this is what you need to do:</div><div><br></div><div>(1) Ensure that the topol.top and index.ndx files are set up in the manner described on the wikki.</div><div><br></div><div>(2) Put in the mdp file the following:</div>
<div><br></div><div><span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; ">energygrps = A B C<br>energygrp_table = A C B B B C</span></div><div><br></div>
<div>(3) Now create a set of files as follows</div><div><br></div><div>table.xvg - containing the A B interactions</div><div>table_A_C.xvg - containing the A C interactions</div><div>table_B_B.xvg - containing the B B interactions</div>
<div>table_B_C.xvg - containing the B C interactions</div><div><br></div><div>You now run with</div><div><br></div><div>mdrun -table table.xvg</div><div><br></div><div>gromacs will read in all these files as energygrp_table tells it that it must read table_A_C.xvg, table_B_B.xvg and table_B_C.xvg as well as table.xvg.</div>
<div><br></div><div>To put this all another way, when you use the command</div><div><br></div><div>vdw-type = user</div><div><br></div><div>(incidentally I think you said that there are no electrostatics so you can use <span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; ">coulombtype = cut-off as it will make no difference)</span></div>
<div><br></div><div>Then you have to provide a table file with definitions of f(r), g(r) and h(r). Then gromacs it will use this table for all the non-bonded interactions (so in your case AA, AB, AC, BB, BC and CC interactions). However, when you add:</div>
<div><br></div><div><span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; ">energygrp_table = A C</span></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Then you are saying, "oh but please do something different for A C interactions, do what is in the table_A_C.xvg file instead." </span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">As a side note, correct me if I am wrong, but there are no AA or CC interactions in your molecule? With the set up I describe above they would have a functional form equivalent to the A B interaction. This is obviously not important if their parameters are equal to zero.</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Finally, if your still not sure that stuff has been read in you can look in the log file, which should say what tables have been read in. It will also give you warnings (although these maybe in standard error I can't remember) if your tables are not set up properly.</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">ciao</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Gareth</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></font></div><div><br></div><div><br><div class="gmail_quote">On Tue, Apr 20, 2010 at 7:38 PM, Martin Vartorelli <span dir="ltr"><<a href="mailto:mrvartorelli@gmail.com">mrvartorelli@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi Gromacs users,<br>
<br>
previously I've posted my situation with a simulation and I'm still confused:<br>
<br>
I have 3 groups of different atoms: A, B and C, and tabulated bonding<br>
and non-bonding potentials.<br>
<br>
The tabulated bonded potentials are specified in the files:<br>
<br>
table_b1.xvg<br>
table_a1.xvg<br>
table_a2.xvg<br>
table_a3.xvg<br>
<br>
Because the systems is simply this molecule: A--(B)n--C, the *.mdp<br>
file contains the following:<br>
<br>
coulombtype = user<br>
vdw-type = user<br>
energygrps = A B C<br>
energygrp_table = A B A C B B B C<br>
<br>
>From the manual and tutorials I understand that, for the non bonding<br>
part of U, I can suply a file called "table.xvg" that contains the 7<br>
non-zero columns for r, f, f', g, g', h and h'. The file "table.xvg"<br>
must contain the functions for AB interactions and the remaining<br>
interactions must be in the files "table_A_C.xvg", "table_B_B.xvg" and<br>
"table_B_C.xvg".<br>
<br>
<br>
So, in the command line I can write: mdrun .... -table table.xvg -tableb table<br>
<br>
I'm right?<br>
<br>
As a test I have also included the file "table_A_B.xvg" (the same<br>
potential functions that "table.xvg") in the same directory, and when<br>
I look at the log file it shows that this file is readed too. I'm<br>
still unclear... how must I specify the files in an univocal way for<br>
the non bonding interactions that I want?<br>
<br>
Gromacs version: 4.0.7<br>
<br>
Any help will be appreciated.<br>
<br>
Martin.<br>
<br>
PD: I've received some private e-mails explaining to me some issues<br>
about spam, mailing lists, private mails, etc. Thanks to all and I<br>
apologize if I've bother someone, it wasn't my intention. And for<br>
those who spontaneously shows me their support, I don't want it and I<br>
don't need it!<br>
<font color="#888888">--<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'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" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</font></blockquote></div><br></div>