Simulations were performed under NVT or NpT condititions at the temperatures, densities, or pressures indicated in Table . Various system sizes are used for orthorhombic box geometries with box length ratios of L z /L x = L z /L y ≈ 2.7933596497. All MD simulations were performed using GROMACS 5.0.6. The integration time step for all simulations was 2 fs, except for the rigid "coarse-grained" ILM2 model for [BMim][PF 6 ], where it was increased to 10 fs. The temperature of the simulated systems was controlled employing the Nosé-Hoover thermostat with a coupling time τ T = 1.0 ps. Constant pressure simulations were realized using a Rahman-Parrinello-barostat employing τ p = 2.0 ps. Note that the compressibility χ T was adjusted to each of the different systems. Both, the Lennard-Jones and electrostatic interactions were treated by smooth particle mesh Ewald summation. The Ewald convergence parameter was set to a relative accuracy of the Ewald sum of 10 -5 for the Coulomb interaction and 10 -3 for the LJ interaction. All bond lengths were kept fixed during the simulation runs and distance constraints were solved by means of the LINCS algorithm , except for TIP4P/2005 water, where the SETTLE procedure was used. The simulations were carried out in n X/W = 80 (except one simulation with n W = 320 for water) subsequent segments, indicated in Tables II III, and IV. All reported properties were then calculated for those segments separately in order to be able to estimate the uncertainty using standard statistical analysis procedures. Finally, we would like to shortly discuss how the "Or-thoBoXY" MD cells reported here were constructed: all simulation boxes listed in Tables III, and IV originated from cubic simulation cells containing either 256, 512, 1024, or 2048 molecules. These cells were created from scratch using Packmol. Then, three of those cells were stacked on top of each other, followed by a short simulation, employing Berendsen thermostats and barostats to robustly produce an initial system with a density close to the equilibrium value. The boxes were then morphed into their final form in a volume-preserving fashion by using short nonequilibrium MD simulation runs of 200 ps length by employing GROMACS' "deform" feature. The prepared systems were then equilibrated under proper NpT conditions for simulation times τ X = f X • ∆τ X . Here the ∆τ X are listed in Table , and f X = 2500. The simulation boxes used for all NVT simulations of water listed in Table were taken from Ref. V. RESULTS AND DISCUSSION