Molecular dynamics (MD) simulations for cubic cells including 200 ionic pairs were performed with the AMOEBA-IL force field using the GPU version of TINKER 8.7, following a stepwise procedure for heating the system to the target temperature (the temperature with available experimental data for comparison). For each temperature, 1 ns, NPT simulations with an integration step of 1 fs, non-bonded cutoff of 10 Å and sampling each 10 ps were performed. For the heating process, the RESPA integrator, Bussi thermostat, and Berendsen barostat were employed. Long-range electrostatic effects were calculated via the smooth particle mesh Ewald (sPME) method, with a 10 Å real-space cutoff. The dimensions for the sPME reciprocal space grid were set as the closest integer to 1.2 times the axis of the cubic cell. Isotropic long-range corrections for the van der Waals term were used to account for van der Waals interactions beyond a cutoff distance of 10 Å. Test calculations with van der Waals cutoffs of 8 and 14 Å were performed to determine the influence of the cutoff.