All MD simulations were run using OpenMM 54 using its TIP4P/Ice force field implementation. The long-range electrostatics were treated using the Particle-Mesh Ewald (PME) algorithm, with an error tolerance parameter of 1e-5, a real space cutoff distance of 0.31 nm, and a switching distance of 0.279 nm. The cutoff (0.31nm) was chosen considering the shortest box length in supercells 211 and 221 of ice Ic, which was 0.6358 nm.