Is there a way within GROMACS to use both Buckingham and Lennard-Jones potentials in the same simulation? So for example water-water could be modelled with LJ and water-solute with Buckingham?