Hi,&nbsp; Justin<br>I remove the del command. But I still got an error.<br><br>make_ndx -f a.pdb&nbsp; -o a.ndx &lt; make_ndx.input <br><br>make_ndx.input:<br><br>chain A and B&amp; !a H*<br>chain C&amp; !a H*<br>q<br><br><br><br>checking input for internal consistency...<br>
Opening library file /export/apps/share/gromacs/top//ffG53a6.itp<br>Opening library file /export/apps/share/gromacs/top//ffG53a6nb.itp<br>Opening library file /export/apps/share/gromacs/top//ffG53a6bon.itp<br>Opening library file /export/apps/share/gromacs/top//ff_dum.itp<br>
Generated 165 of the 1596 non-bonded parameter combinations<br>Opening library file /export/apps/share/gromacs/top//spc.itp<br>Opening library file /export/apps/share/gromacs/top//ions.itp<br>Excluding 3 bonded neighbours molecule type &#39;Protein_A&#39;<br>
Excluding 3 bonded neighbours molecule type &#39;Protein_B&#39;<br>Excluding 3 bonded neighbours molecule type &#39;Protein_C&#39;<br><br>NOTE 1 [file r-l_365130_G53a6.minim_traj_withH_0.6.top, line 39]:<br>&nbsp; System has non-zero total charge: -7.000002e+00<br>
<br><br><br>processing coordinates...<br>double-checking input for internal consistency...<br><br>WARNING 1 [file r-l_365130_G53a6.minim_traj_withH_0.6.top, line 39]:<br>&nbsp; With nstlist=0 atoms are only put into the box at step 0, therefore<br>
&nbsp; drifting atoms might cause the simulation to crash.<br>renumbering atomtypes...<br>converting bonded parameters...<br>initialising group options...<br>processing index file...<br>Making dummy/rest group for T-Coupling containing 4655 elements<br>
Making dummy/rest group for Acceleration containing 4655 elements<br>Making dummy/rest group for Freeze containing 4655 elements<br>Making dummy/rest group for Energy Mon. containing 1045 elements<br>Making dummy/rest group for VCM containing 4655 elements<br>
Number of degrees of freedom in T-Coupling group rest is 13962.00<br>Making dummy/rest group for User1 containing 4655 elements<br>Making dummy/rest group for User2 containing 4655 elements<br>Making dummy/rest group for XTC containing 4655 elements<br>
Making dummy/rest group for Or. Res. Fit containing 4655 elements<br>Making dummy/rest group for QMMM containing 4655 elements<br>T-Coupling&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; has 1 element(s): rest<br>Energy Mon.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; has 3 element(s): chAANDB_&amp;_!H* chC_&amp;_!H* rest<br>
Acceleration&nbsp;&nbsp;&nbsp;&nbsp; has 1 element(s): rest<br>Freeze&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; has 1 element(s): rest<br>User1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; has 1 element(s): rest<br>User2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; has 1 element(s): rest<br>VCM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; has 1 element(s): rest<br>XTC&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; has 1 element(s): rest<br>
Or. Res. Fit&nbsp;&nbsp;&nbsp;&nbsp; has 1 element(s): rest<br>QMMM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; has 1 element(s): rest<br>Checking consistency between energy and charge groups...<br><br>-------------------------------------------------------<br>Program grompp, VERSION 4.0.2<br>
Source code file: grompp.c, line: 150<br>Fatal error:<br>atoms 1 and 2 in charge group 1 of molecule type &#39;Protein_A&#39; are in different energy groups<br>-------------------------------------------------------<br><br>
<br><br><div class="gmail_quote">On Fri, Jan 30, 2009 at 5:11 PM, Justin A. Lemkul <span dir="ltr">&lt;<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</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="Ih2E3d"><br>
<br>
Liu Shiyong wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
-------------------------------------------------------<br>
Program grompp, VERSION 4.0.2<br>
Source code file: grompp.c, line: 150<br>
<br>
Fatal error:<br>
atoms 1 and 2 in charge group 1 of molecule type &#39;Protein_A&#39; are in different energy groups<br>
-------------------------------------------------------<br>
</blockquote>
<br></div>
Check your structure to make sure you are using the right coordinate file.<br>
<br>
Try to *not* delete every other group in the index file. &nbsp;Just because you don&#39;t need these groups in your energygrps does not mean they are not necessary for other functions.<br>
<br>
-Justin<br>
<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="Ih2E3d">
<br>
&quot;Can someone please tell Icarus that he&#39;s not the only one falling from the sky?&quot; (Urban Dance Squad)<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; :-) &nbsp;G &nbsp;R &nbsp;O &nbsp;M &nbsp;A &nbsp;C &nbsp;S &nbsp;(-:<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;GROningen Mixture of Alchemy and Childrens&#39; Stories<br>
<br>
<br>
I checked &nbsp;the a.ndx<br>
<br>
[ System ]<br>
 &nbsp; 1 &nbsp; &nbsp;2 &nbsp; &nbsp;3 &nbsp; &nbsp;4 &nbsp; &nbsp;5 &nbsp; &nbsp;6 &nbsp; &nbsp;7 &nbsp; &nbsp;8 &nbsp; &nbsp;9 &nbsp; 10 &nbsp; 11 &nbsp; 12 &nbsp; 13 &nbsp; 14 &nbsp; 15<br>
 &nbsp;16 &nbsp; 17 &nbsp; 18 &nbsp; 19 &nbsp; 20 &nbsp; 21 &nbsp; 22 &nbsp; 23 &nbsp; 24 &nbsp; 25 &nbsp; 26 &nbsp; 27 &nbsp; 28 &nbsp; 29 &nbsp; 30<br>
<br>
[ chAANDB_&amp;_!H* ]<br>
 &nbsp; 1 &nbsp; &nbsp;5 &nbsp; &nbsp;6 &nbsp; &nbsp;7 &nbsp; &nbsp;8 &nbsp; &nbsp;9 &nbsp; 10 &nbsp; 11 &nbsp; 12 &nbsp; 14 &nbsp; 15 &nbsp; 16 &nbsp; 17 &nbsp; 18 &nbsp; 19<br>
 &nbsp;23 &nbsp; 24 &nbsp; 25 &nbsp; 26 &nbsp; 27 &nbsp; 28 &nbsp; 29 &nbsp; 30 &nbsp; 31 &nbsp; 32 &nbsp; 34 &nbsp; 35 &nbsp; 36 &nbsp; 37 &nbsp; 38<br>
<br>
[ chC_&amp;_!H* ]<br>
2571 2575 2576 2577 2578 2580 2581 2582 2583 2584 2585 2586 2587 2588 2589<br>
2590 2591 2592 2593 2594 2595 2597 2598 2599 2600 2601 2602 2605 2606 2607<br>
<br>
<br>
<br>
-- <br>
Shiyong Liu<br>
Postdoc<br>
center for bioinformatics in the university of kansas<br>
Lab: (785)864-1962<br></div>
Email: <a href="mailto:syliu@ku.edu" target="_blank">syliu@ku.edu</a> &lt;mailto:<a href="mailto:syliu@ku.edu" target="_blank">syliu@ku.edu</a>&gt; (<a href="mailto:shiyongliu@ku.edu" target="_blank">shiyongliu@ku.edu</a> &lt;mailto:<a href="mailto:shiyongliu@ku.edu" target="_blank">shiyongliu@ku.edu</a>&gt; or <a href="mailto:liushiyong@ku.edu" target="_blank">liushiyong@ku.edu</a> &lt;mailto:<a href="mailto:liushiyong@ku.edu" target="_blank">liushiyong@ku.edu</a>&gt;)<div class="Ih2E3d">
<br>
Homepage: <a href="http://www.people.ku.edu/%7Esyliu" target="_blank">http://www.people.ku.edu/~syliu</a><br>
Lab: &nbsp; &nbsp;<a href="http://vakser.bioinformatics.ku.edu/people" target="_blank">http://vakser.bioinformatics.ku.edu/people</a><br>
Phone: &nbsp; &nbsp; &nbsp;(785) 864-1962<br>
<br>
<br></div>
------------------------------------------------------------------------<br>
<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>
</blockquote>
<br>
-- <br>
========================================<br>
<br>
Justin A. Lemkul<br>
Graduate Research Assistant<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<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>
</blockquote></div><br><br clear="all"><br>-- <br>Shiyong Liu<br>Postdoc<br>center for bioinformatics in the university of kansas<br>Lab: (785)864-1962<br>Email: <a href="mailto:syliu@ku.edu">syliu@ku.edu</a> (<a href="mailto:shiyongliu@ku.edu">shiyongliu@ku.edu</a> or <a href="mailto:liushiyong@ku.edu">liushiyong@ku.edu</a>)<br>
Homepage: <a href="http://www.people.ku.edu/~syliu">http://www.people.ku.edu/~syliu</a><br>Lab: &nbsp; &nbsp;<a href="http://vakser.bioinformatics.ku.edu/people">http://vakser.bioinformatics.ku.edu/people</a><br>Phone:         &nbsp; &nbsp; &nbsp;(785) 864-1962<br>