HITS-MCM / gromacs-ramd

Random Acceleration Molecular Dynamics in GROMACS
GNU Lesser General Public License v2.1
33 stars 7 forks source link

Residence time #42

Open BerndDoser opened 4 months ago

BerndDoser commented 4 months ago

Additional atomic group pairs will be added to determine the residence time between the ligand and receptor. The atomic group pairs receptor_group_<n> and ligand_group_<n> must be defined in the index file (.ndx), e.g.

[ receptor_group_1 ]
   1    2
[ ligand_group_1 ]
   4    5

During the simulation, the distances between the center of mass (COM) of each pair will be calculated at every RAMD evaluation step. If the average of all distances exceeds the parameter 'ramd-residence-dist,' a message will be logged in the simulation file.

==== RAMD ==== Residence time obtained after <x> steps.

This message will only printed once.

Two new MDP options will be integrated:

The COM distances for the residence time will be saved in the last columns of the xvg file using the -ramd flag.