The selected solvents were [TFSI] ionic liquid, EC, DMC and 1M LiTFSI as lithium salt, with different mixing ratios as shown in Table . The bonded and non-bonded parameters for EC and DMC were obtained from the OPLS/AA force fields . The parameters of Lopes and Padua are based on the all atom optimized potential for liquid system force field (OPLS/AA) of Jorgensen et al. , which is one of the most general and widely used force fields used to describe ILs. Therfore, cations and anions of ILs ([EMI] + and [TFSI] -) were modeled using parameters of the systematic all-atom force field of Lopes and Padua . The Li ion parameters were adopted from ion-water potential calculations . Charges are downscaled to an absolute value of 0.8 to account for the charge transfer and polarizability within ILs. Charge scaling is an efficient way to approximately include polarizability and charge-transfer effects in classical simulations at no further computational cost .. Morrow and Maginn were the first to report a partial charge assignment for the [BMIM][PF 6 ] ion pair based on a noninteger molecular charge of ±0.904e. Subsequent work determined a molecular charge scaling of ±0.8e to be generally ideal for reproducing experimental ionic liquid-phase properties including enthalpy of vaporization, surface tension, and self-diffusion coefficients . In another work by Doherty et al, radial distribution functions (RDFs) were used to examine cation-anion interactions employing 0.8*OPLS-2009IL, and the scaled force field accurately reproduced RDFs from ab initio MD simulations .