Radial distribution functions (RDFs) between the centers of mass of the ions, the solvent and between various atoms of the ions have been calculated in order to study the microstructure of the ILs. Thus, we focus here on the cation-cation, anion-anion, cation-anion, solvent-solvent distribution and the specific atoms comprising the two moieties in order to analyze and provide an insightful average picture of the IL/carbonate structure. For the [EMI] + cation ring centers, very weak but clear ordering can be observed. There seems to be a first shell of neighboring cations at around 250 pm for the pure IL (system A). This finding is different from other previously reported results about simulation studies of more polar and hydrophilic IL 1-ethyl-3-methylimidazolium acetate with water where no clear ordering was observed . Comparing the different systems A, B, C and D that have different levels of ionic liquid/carbonate concentration , the intensity (height) of the peak is increased from system A, B, C to D upon increasing the content of carbonate in the mixture indicating that the ordering between the cations is enhanced by the addition of carbonate in the mixture (Figure ). Our results in this study show that the first maxima of the ring -ring distance appear at smaller distances than for other aromatic systems, e.g. benzene, where the first maximum is found at to the atom H 1 is much more distinct and tighter than for the other two atoms, suggesting that H 1 is a stronger hydrogen bond donor than H 4 and H 5 for the pure IL system (system A). For panel (c), the first peak is much higher for coordination of the N 1 atom of [TFSI] -anion to the atom H 1 of the [EMI] + cation for the pure IL than for system D and the height of the first coordination shell is decreased when carbonate is added. In sharp contrast to panel (c), for panel (d), the first peak is much higher for coordination of the N 1 atom of [TFSI] -anion to the atom H4/5 of the