[gmx-developers] Position Restraints in Thermodynamic Integration

[Carruthers, James]

Hello,

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,
James