The classical forcefield used to model water molecules in the polymorphic forms of ice investigated was TIP4P/Ice: a four-body potential that, alongside the oxygen and two hydrogen atoms, employs an additional virtual atom (r M ). The dynamics of the virtual atom do not contribute any additional DOF to the system because its position can be reconstructed without uncertainty as: