Hi, I'm new to OpenMM. Could someone help on adding distance constrain between the mass centers of two molecules? e.g.

f= k(r-ro)^a, r>r0
f=0,r<=r0

r is the distance between the centers of two molecules.

How should define the external force in openmm?

Thank you very much!