Is Gromacs able to calculate the work and forces of a restrained bond (a NMR restrain) moved with a constant rate between two distances to create Jarzinsky profiles just as Amber does? If not, which subroutine you recommend I should modify to do this?