I am trying to add an additional harmonic restraint to a system to calculate the free energy of rotation. I am using the following approach: ( where kx and x0 are set to numerical values previously ) torsion = CustomTorsionForce('kx*(x - x0)^2') torsion.addPerTorsionParameter('x0') torsion.addPerTor...