Experimental techniques such as infrared spectroscopy (IR) and raman spectroscopy have often been the methods of choice to study ion-ion and ion-solvent coordination structures. However, these spectroscopic techniques have several limitations for multicomponent solvent systems. These experimental techniques, for example, have difficulties in quantitatively deconvoluting the overlapping and the overtone peaks of the individual component systems . On the other hand, theoretical methods such as ab initio molecular dynamics (AIMD) and classical molecular dynamics (MD) simulations have proven useful in modeling the solvation and transport properties of multicomponent solvent systems. MD simulations have proven useful in characterizing microscopic structural features of Li + -anion coordination. MD simulations can help connect microscopic structural features of Li +anion coordination to the underlying dynamical properties such as Li mobility. To this end, ab initio molecular dynamics (AIMD) simulations may add valuable insight into solvation mechanisms and energetics of ionic liquid systems, but this comes at the cost of relatively higher computational demands . The high computational cost of AIMD often limits simulations to short timescales of shorter simulation lengths which consequently prohibits description of the medium-long structural order present in ionic liquids . Classical molecular dynamics (MD), therefore, has largely been the method of choice for investigations of IL structure . On the basis of the above considerations, therefore, the purpose of the present study is to provide some general insight into the solvation structures of both [EMI] + and Li + ions in pure and mixed IL/carbonate solvent systems, by employing classical molecular dynamics simulations. This study can help to better understand the solvation phenomena in these electrolytes, significantly improving the rational design of electrolytes for battery applications.