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) . Although MD simulations using nonpolarizable force fields provide important insight into molecular level correlations and structure in ionic systems . However, nonpolarizable force fields do have poor dynamics , which has led to critical evaluations. As an alternative to the nonpolarizable force fields, methods such as AIMD , and polarizable force fields have varied advantages over their nonpolarizable counterparts, but this comes at a cost of higher computational demands.