We have investigated SEI-formation reactions within a time-scale of few nanoseconds based on our rs@md simulation of the graphite/carbonate electrolyte model. Prior to the reactions the structural and dynamic properties of the model were consistently validated in terms of the baseline density and the diffusivity values. We have suggested the formation of negatively charged Li + /EC -structures as well as the flux of Li + s from the bulk electrolyte to the electrode as rate limiting steps for the EC reduction reaction. The effect of the Li + flux or "concentration polarization" was observed at high reaction rates. Since concentration polarization occurs in general at much longer times, this could be a consequence of the relatively large electron transfer rates used in our simulation. Nevertheless, our results give insights into the role of electro-neutrality and the structure of the EDL in determining the effective rates. Our results were based on the Li + injection scheme accounting for the charge transfer reactions at the interface to maintain electro-neutrality. Importantly, the surface charge at the graphite anode was assumed to be fixed in our simulations and a treatment accounting for charge fluctuations at the electrode surface due to the reactions could be considered in future work to account for the proper voltage dependence of the decomposition process. As the elementary electron-transfer rates were chosen based on the accessible time-range of our simulation, the derivation of these elementary rates from more rigorous frameworks (i.e Marcus theory) would be viewed as complementary to our work. Overall, our work shows how gas-phase energy barriers can be used to derive liquid-phase reaction rates based on the rs@md framework. The implications of the changed chemical environments on Li + diffusivity and coordination was also discussed. Next steps could be the consideration of more complex or more realistic electrolyte formulations, with the attempt of identifying the functionality of known additive compounds and their effect on the calculated growth rates. Using the derived structures to calculate observable mechanical properties (i.e bulk modulus) could also be of interest.