<div>Hi, Mark:</div>
<div>&nbsp;</div>
<div>I made the index file and call two methanes as methane1 and methane2 , then I set energygrp_excl = methane1 methane2, but when I run grompp, it said methane1 is not an energy group. What do you understand this problem? Thanks. -Xianghong Qi<br>
<br></div>
<div class="gmail_quote">On Wed, Nov 19, 2008 at 2:17 AM, xianghong qi <span dir="ltr">&lt;<a href="mailto:xianghong001@gmail.com">xianghong001@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div>Yes, Mark. I now understand. I use united atom model. In my system, protein includes only methane pair.&nbsp; I still need those parameter to control protein/solvent, solvent/solvent interaction. Thanks so much for your great help.&nbsp; I will do it now. </div>

<div>&nbsp;</div>
<div>Best, -Xianghong Qi<br><br></div>
<div>
<div></div>
<div class="Wj3C7c">
<div class="gmail_quote">On Wed, Nov 19, 2008 at 1:48 AM, Mark Abraham <span dir="ltr">&lt;<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div>xianghong qi wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Thanks, Mark. &nbsp;That means I can set energygrp_excl = Protein Protein, then I can get rid of &nbsp;those parameter which are related to electrostatics, Vdw( like rcoulomb, vdw-type, rvdw, rvdw-switch, etc). &nbsp;Am I right?<br>
</blockquote><br></div>No. Section 7.3.19 indicates how this should be done. I think you want<br><br>energygrp_excl = Methane1 Methane2<br><br>so you will need to use make_ndx suitably to create these groups.<br><br>If &quot;Protein&quot; includes only your pair of methanes, then<br>
<br>energygrp_excl = Protein Protein<br><br>will have the same effect, only because there are no intramolecular nonbonded interactions in methane. It&#39;s thus not a generalizable solution... 
<div>
<div></div>
<div><br><br>Mark<br>_______________________________________________<br>gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br><a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.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 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/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></div></div>
<div>
<div></div>
<div class="Wj3C7c">-- <br>Some people make the world more special just by being in it.<br></div></div></blockquote></div><br><br clear="all"><br>-- <br>Some people make the world more special just by being in it.<br>