[gmx-users] select an atom for "only" residue with make_ndx
jalemkul at vt.edu
Wed May 15 18:13:57 CEST 2019
On 5/15/19 11:03 AM, ABEL Stephane wrote:
> Hello all,
> I am currently doing simulations of a model of membranes with different types of lipids (LIPIDA, LIPIDB and LIPIDC) and I would like to select a group of atoms for only the same lipid with make_ndx. Below an example
> Three different lipids : LIPIDA, LIPIDB and LIPIDC contain few atoms with the same name (C31, C32 and C33). With make_ndx I would like to have the atoms C33 and C32 for "only" the LIPIDC
> So If I do
> rLIPIDC & aC32 | aC33 ---> I have the indices of the C32 and C33 of the 3 lipid types.
> So what is the exact command
Each different [moleculetype] already has its own default group, so
choosing atoms in LIPIDC should be as simple as:
N & a C31 | a C32 | a C33
where N is the group corresponding to LIPIDC.
Justin A. Lemkul, Ph.D.
Office: 301 Fralin Hall
Lab: 303 Engel Hall
Virginia Tech Department of Biochemistry
340 West Campus Dr.
Blacksburg, VA 24061
jalemkul at vt.edu | (540) 231-3129
More information about the gromacs.org_gmx-users