Hi everyone,<br><br>I'm trying to figure out the right commands for make_ndx so to generate the correct index.ndx file for g_bond analysis. From my PDB file:<br><br>ATOM 1 BN0 GGA X 1 107.255 17.133 66.368 1.00 0.00 P1 C<br>
ATOM 2 SC1 GGA X 1 110.304 14.942 66.116 1.00 0.00 P1 C<br>ATOM 3 SC2 GGA X 1 113.632 15.318 65.904 1.00 0.00 P1 C<br>ATOM 4 SC3 GGA X 1 114.063 12.624 65.310 1.00 0.00 P1 C<br>
ATOM 5 BN0 GGA X 2 108.171 20.156 68.028 1.00 0.00 P1 C<br>ATOM 6 SC1 GGA X 2 109.535 18.753 71.242 1.00 0.00 P1 C<br>ATOM 7 SC2 GGA X 2 112.445 17.547 72.396 1.00 0.00 P1 C<br>
ATOM 8 SC3 GGA X 2 111.301 16.592 74.758 1.00 0.00 P1 C<br>ATOM 9 BN0 GGA X 3 110.843 21.726 66.258 1.00 0.00 P1 C<br>ATOM 10 SC1 GGA X 3 111.984 24.201 68.852 1.00 0.00 P1 C<br>
ATOM 11 SC2 GGA X 3 114.049 24.264 71.496 1.00 0.00 P1 C<br>ATOM 12 SC3 GGA X 3 113.142 26.763 72.351 1.00 0.00 P1 C<br>...<rest of PDB file not included>...<br><br>Say I want to group together atoms SC1, SC2, and SC3 from all GGA residues in make_ndx:<br>
<br> 0 System : 19699 atoms<br> 1 GGA : 472 atoms<br> 2 DRG : 324 atoms<br> 3 W : 18669 atoms<br> 4 ION : 234 atoms<br><br>I typed in: <br>
<br>> r GGA & t SC1 | r GGA & t SC2 | r GGA & t SC3<br><br>Found 472 atoms with residue name GGA<br>Found 36 atoms with type SC1<br>Merged two groups with AND: 472 36 -> 0<br>Found 472 atoms with residue name GGA<br>
Merged two groups with OR: 0 472 -> 472<br>Found 24 atoms with type SC2<br>Merged two groups with AND: 472 24 -> 0<br>Found 472 atoms with residue name GGA<br>Merged two groups with OR: 0 472 -> 472<br>Found 0 atoms with type SC3<br>
<br> 5 GGA_&_SC1_GGA_&_SC2_GGA: 472 atoms<br><br>However, make_ndx "found 0 atoms with type SC3". I'm not sure if I'm doing this correctly. <br><br>If anyone has any insight, I would be most appreciative.<br>
<br>Thanks,<br>L<br><br>