I am trying to implement a fraction-of-native-contacts (Qf) force to apply it in umbrella sampling calculations. First, I am using a CustomBondForce to add each native contact contribution as: energy_function = '(1/N)*(1/(1+exp(beta*(r-lambda*r0))))' Qf = CustomBondForce(energy_function) Where N is ...