The performance of lithium-ion batteries (LIB) using organic electrolytes strongly depends on the formation of a stable solid electrolyte interphase (SEI) film. Elucidating the dynamic evolution and spatial composition of the SEI can be very useful to study the stability of SEI components and help optimize the formation cycles of LIB. We propose a classical molecular dynamics simulation protocol for predicting the first stages of SEI formation using a reaction method involving the decomposition of EC and LiPF 6 molecules in the electrolyte. We accelerate the formation of SEI components near the anode surface by increasing the probability of reactions, implemented through a geometry matching scheme, followed by a force-field reconfiguration. We observe the formation of gases (C 2 H 4 ), inorganic (Li 2 CO 3 and LiF) and organic (LEDC) components. This protocol shows promise to be able to evaluate the effects of varying electrolyte compositions and additives on SEI layer structure and composition.

We quantify the effect of various sources of uncertainties in the prediction of thermo-physical properties of polymers using molecular dynamics simulations. We quantify how the choice of polymer builder, force field, molecular weight and data analysis affect predicted values of the glass transition temperature (Tg), room temperature density and coefficient of thermal expansion of poly(methyl-methacrylate) (PMMA) and polystyrene (PS). Interestingly, we find that the data analysis introduces significant uncertainties in Tg (approximately 5%) while the other properties are insensitive to it. The force field is the only variable that significantly affects the predictions of density. Polymer-consistent force field (PCFF) resulted in a higher density for PMMA than Dreiding and the opposite trend was observed in PS; in all cases the difference in density was less than 2%. Strongly correlated with density, we find that PCFF leads to a higher Tg than Dreiding for PMMA and both force fields predict similar Tg values for PS. The trends in Tg can be explained by differences in segmental mobility of the melts predicted by the two force fields. We find that the choice of amorphous polymer builder results in uncertainties in predictions comparable to those associated with the force field due to subtle, but persistent, differences in molecular structure. The results presented here provide insight into the physics behind molecular simulations of polymers and quantitative levels of uncertainties associated with individual sources that can help practitioners of molecular simulations interested in using their results in engineering applications.

No abstract