I like to apply harmonic constraints on all c-alpha atoms during MD-vacuo, with a condition apply to maintain inter c-alpha atom distance should be less than or equal to certain threshold (eg:9 Angs) during the integration.

Are there any method exist to call harmonic constraints for the above mentioned condition.

Appreciate your help