&nbsp;Dear Gromacs users,<br> <br> &nbsp;I would like to compute force of interaction between 2 molecules as a function of relative position and, in particular, orientation.<br> <br> &nbsp;I was wondering if I could modify the pull code to include an angular harmonic<br> &nbsp;potential. Would that be possible, and if so, does anyone have a suggestion as<br> &nbsp;to where and how the modifications should be done?<br> <br> &nbsp;Secondly, but still connected to the above, does anyone have an idea to how<br> &nbsp;such an angle should be defined?<br> <br> &nbsp;All inputs are very welcome and highly appreciated.<br> <br> &nbsp;Best regards to all,<br> &nbsp;Soren Enemark<br>