In this work, MD simulations were performed using LAMMPS (Large Scale Atomic/Molecular Massively Parallel Simulator) . All potential models used to describe intermolecular interactions were based on the OPLS/AA force field (Optimized Potentials for Liquid Simulations/All Atom) . In the non-polarizable OPLS/AA force field formalism, the total energy for each ionic liquid system are evaluated as a sum of individual energies for the harmonic bond stretching and angle bending terms, a cosine series for torsional energetics, and Coulomb and 12-6 Lennard-Jones terms for the nonbonded interactions. Geometric combining rules for the Lennard-Jones coefficients were utilized. Nonbonded interactions were evaluated intermolecularly and for intramolecular atom pairs separated by three or more bonds. To apply the same parameters for both intra-and intermolecular interactions, the 1,4-intramolecular interactions were reduced by a factor of 2. A 12-6 Lennard-Jones (LJ) pair interaction model with 15.0 Å cut-off was used to model the van der Waal's interactions.