[gmx-developers] Position Restraints in Thermodynamic Integration

Carruthers, James j.carruthers2 at lancaster.ac.uk
Wed Sep 27 14:46:39 CEST 2017


In the online user guide for 2016.4 it states in mdp options that the restraint-lambdas directive only controls dihedral and pull code restraints. However, the 2016.4 manual states that position restraints can be interpolated for free energy calculations (pp. 137,138). I would like to be able to control position restraints as this will enable a convenient way to implement the Einstein crystal method of calculating solid state chemical potentials. This raises two questions, is there a technical reason why restraint-lambdas is restricted to controlling only those two parameters and is there a way to implement scaling position restraints outside of the *-lambdas framework that can be analysed by gmx bar?

Best regards,

More information about the gromacs.org_gmx-developers mailing list