After equilibration, the systems were subjected to 10 ns of MD in the NPT ensemble. In all cases, the RESPA integrator, Bussi thermostat, and Montecarlo barostat were used for all IL systems studied. The integration step used was 1 fs and sampling was done each 1 ps. A non-bonded interaction cutoff distance of 10 Å was used, coupled with PME. PME summation was set as the closest integer to 1.2 times the axis of the cubic cell.