Drude for GROMACS

But in this case I would also need parameters for PD1A and OD2C2B cause they’re not present in current release and neither in the charmm version

PD1A has been part of the force field for a long time. PD1AN was simply removed and supplanted by PD1A. For OD2C2B, ε = -0.19 kcal/mol, Rmin/2 = 1.97 Å. If you need further parameters or information about the implementation, again, please contact Alex’s group. They’re responsible for that publication and the information it contains.