The selected solvents were 1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl) imide [TFSI] ionic liquid, ethylene carbonate (EC), diethyl carbonate (DEC) and 1M LiFSI as lithium salt. The bonded and non-bonded parameters for EC and DMC were obtained from the OPLS-AA force fields . The cations and anions of ILs ([EMI] + , [TFSI] -) and [FSI] - were modeled using parameters of the systematic all-atom force field of opes and adua he i ion parameters were adopted from ion-water potential calculations Compositions of the simulated ionic liquids are shown in (1.52 g/cm 3 ). Periodic boundary conditions were applied on all three sides of the cubic simulation box to represent the bulk solution. The system was simulated at temperature, T = 298 K and pressure, P = 1 atm. We used the velocity erlet integrator with a timestep of 1 fs to integrate the e uations of motion he ose-oover thermostat and barostat were used to keep the temperature and pressure constant. We used PACKMOL to create initial configurations. The Ewald method was used for treating electrostatic interactions.