Simulations were performed in the canonical ensemble where the temperature was maintained by a Langevin thermostat 40 using a damping parameter of 100 fs. Since the origin of the correction term -C ij 4 /r 4 ij is electrostatics, hence, for our continuous, homogeneous, and isotropic dielectric medium, the correction term is scaled by the respective dielectric constant of the medium, ϵ.