Theory of Cation Solvation in the Helmholtz Layer of Li-ion Battery Electrolytes

The solvation environments of Li$^+$ in conventional non-aqueous battery electrolytes, such as LiPF$_6$ in mixtures of ethylene carbaronate (EC) and ethyl methyl carbonate (EMC), are often used to rationalize the transport properties of electrolytes and solid electrolyte interphase (SEI) formation. In the SEI, the solvation environments in the compact electrical double layer (EDL) next to the electrode, also known as the Helmholtz layer, determine (partially) what species can react to form the SEI, with bulk solvation environments often being used as a proxy. Here we develop and test a theory of cation solvation in the Helmholtz layer of non-aqueous Li-ion battery electrolytes. First, we validate the theory against bulk and diffuse EDL atomistic molecular dynamics (MD) simulations of LiPF$_6$ EC/EMC mixtures as a function of surface charge, where we find the theory can capture the solvation environments well. Next we turn to the Helmholtz layer, where we find that the main effect of the solvation structures next to the electrode is an apparent reduction in the number of binding sites between Li$^+$ and the solvents, again where we find good agreement with our developed theory. Finally, by solving a simplified version of the theory, we find that the probability of Li$^+$ binding to each solvent remains equal to the bulk probability, suggesting that the bulk solvation environments are a reasonable place to start when understanding new battery electrolytes. Our developed formalism can be parameterized from bulk MD simulations and used to predict the solvation environments in the Helmholtz layer, which can be used to determine what could react and form the SEI.