Hi, I am building my own routine to use the Gromacs API to do MS-EVB. I need to treat O-H bonds of 1 hydronium ion with a Morse Potential. How can I specify in the *.mdp file, to just describe bonding in hydronium using this potential? Thanks for you help, Jackson