[gmx-users] Regarding index .ndx files

Kevin Boyd kevin.boyd at uconn.edu
Mon Sep 17 13:29:25 CEST 2018


Hi!

There are two key points here.

1) Index groups don't have to be mutually exclusive. So, you can have your
20+ indices, but also have an inclusive index that encompasses more or all
of your atoms, which you can then use for tc (or comm removal, etc.). The
only necessary thing is that the indices specified in tc-groups are
complete and exclusive.

2) Gromacs by default does just this - the "system" index Vytautas
mentioned is typically autogenerated and contains every atom in your
system. However, you should see if this is the correct option for you. A
lot of protocols e.g. couple protein separate from the solvent - your
initial setup seems to do something like this. If the only other particles
in your system are solvent, and you have an SOL index, you can just add
that column to your original setup. You might also have ions in your system
- these should be coupled in the same group as the solvent.

Kevin

On Mon, Sep 17, 2018 at 7:14 AM, Vytautas Rakeviius <vytautas1987 at yahoo.com>
wrote:

>  Without going into details and setting:
> tc-grps = System
> should allow to proceed without fatal errors.
>     On Monday, September 17, 2018, 1:58:50 PM GMT+3, Sebastian Muraru <
> sebastian.muraru at grabtop.upb.ro> wrote:
>
>  Hi there! I am a Masters student and this is my first email to the list.
> I would just like to ask whether someone can clarify something for me
> regarding the index .ndx files.
> I am using an nvt.mdp file in order to run an nvt equilibration step. The
> section of interest is the following (currently receiving an error):
>
> ###
> ; Temperature coupling is on
> tcoupl          = V-rescale                ; modified Berendsen thermostat
> tc-grps        = CX ssdna; two coupling groups - more accurate
> tau_t          = 0.1    0.1    ; time constant, in ps
> ref_t          = 300    300    ; reference temperature, one for each
> group, in K
> ###
>
> The current error says that a number of atoms are not part of any of the
> T-Coupling groups. Therefore I have the following questions:
> - Does this mean all atoms (including SOL) must be included in the groups
> specified using 'tc-grps'?
> - If my index.ndx file is very large (contains > 20 groups), do I have to
> manually write them down at 'tc-grps' in the mdp file or is there a way to
> automatically import the groups present in the index.ndx file?
>
> Best,
> Sebastian
> --
> Gromacs Users mailing list
>
> * Please search the archive at https://na01.safelinks.
> protection.outlook.com/?url=http%3A%2F%2Fwww.gromacs.org%
> 2FSupport%2FMailing_Lists%2FGMX-Users_List&data=02%7C01%
> 7Ckevin.boyd%40uconn.edu%7C8a6e0b22f16d450f328308d61c8ed732%
> 7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C636727797172451149&sdata=
> vBhFfQUHwCysT1RvQnt1o3k4UKQIIGhoRsOcFYNIwbw%3D&reserved=0 before posting!
>
> * Can't post? Read https://na01.safelinks.protection.outlook.com/?url=
> http%3A%2F%2Fwww.gromacs.org%2FSupport%2FMailing_Lists&
> data=02%7C01%7Ckevin.boyd%40uconn.edu%7C8a6e0b22f16d450f328308d61c8ed732%
> 7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C636727797172451149&sdata=
> w853EKi8mtyssqK934nZNPRtRzO6pS58FPeT9GKI6%2BQ%3D&reserved=0
>
> * For (un)subscribe requests visit
> https://na01.safelinks.protection.outlook.com/?url=
> https%3A%2F%2Fmaillist.sys.kth.se%2Fmailman%2Flistinfo%
> 2Fgromacs.org_gmx-users&data=02%7C01%7Ckevin.boyd%40uconn.edu%
> 7C8a6e0b22f16d450f328308d61c8ed732%7C17f1a87e2a254eaab9df9d439034
> b080%7C0%7C0%7C636727797172451149&sdata=%2BS%
> 2FOAXbJsBXoCKyLds6H1ABBC6ez1h4tTfGvjIibduI%3D&reserved=0 or send a mail
> to gmx-users-request at gromacs.org.
>
> --
> Gromacs Users mailing list
>
> * Please search the archive at https://na01.safelinks.
> protection.outlook.com/?url=http%3A%2F%2Fwww.gromacs.org%
> 2FSupport%2FMailing_Lists%2FGMX-Users_List&data=02%7C01%
> 7Ckevin.boyd%40uconn.edu%7C8a6e0b22f16d450f328308d61c8ed732%
> 7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C636727797172451149&sdata=
> vBhFfQUHwCysT1RvQnt1o3k4UKQIIGhoRsOcFYNIwbw%3D&reserved=0 before posting!
>
> * Can't post? Read https://na01.safelinks.protection.outlook.com/?url=
> http%3A%2F%2Fwww.gromacs.org%2FSupport%2FMailing_Lists&
> data=02%7C01%7Ckevin.boyd%40uconn.edu%7C8a6e0b22f16d450f328308d61c8ed732%
> 7C17f1a87e2a254eaab9df9d439034b080%7C0%7C0%7C636727797172451149&sdata=
> w853EKi8mtyssqK934nZNPRtRzO6pS58FPeT9GKI6%2BQ%3D&reserved=0
>
> * For (un)subscribe requests visit
> https://na01.safelinks.protection.outlook.com/?url=
> https%3A%2F%2Fmaillist.sys.kth.se%2Fmailman%2Flistinfo%
> 2Fgromacs.org_gmx-users&data=02%7C01%7Ckevin.boyd%40uconn.edu%
> 7C8a6e0b22f16d450f328308d61c8ed732%7C17f1a87e2a254eaab9df9d439034
> b080%7C0%7C0%7C636727797172451149&sdata=%2BS%
> 2FOAXbJsBXoCKyLds6H1ABBC6ez1h4tTfGvjIibduI%3D&reserved=0 or send a mail
> to gmx-users-request at gromacs.org.


More information about the gromacs.org_gmx-users mailing list