Quantum Mechanics-Calibrated Classical Simulation of Earth-Abundant Divalent Metal Bis(trifluoromethanesulfonyl)imide Molar Conductivity in Organic Cosolvents with Onsager Transport Formalism

Molar conductivities are one of the fundamental transport properties of electrolytes that reflect the complex and dynamic interactions between ions and solvents comprehensively. The quantitative accuracy of experimental input-free simulations of molar conductivities is strongly influenced by the underlying interaction parameters employed in the model. When validated with experimental molar conductivities, the developed model could be used to reveal further atomistic level details about the solvation structures and correlated ion pair formation, providing in-depth knowledge about solution physical chemistry and shedding light on electrolyte-solvent system design rules. Divalent cations are more challenging to model than monovalent cations due to their higher charge densities and stronger interactions with the environment. Yet, they started attracting significant attention for next-generation energy storage purposes. In this work, we focus on two earth-abundant divalent cation electrolytes, Mg(TFSI)₂ and Ca(TFSI)₂ in a dimethylacetamide-tetrahydrofuran (DMA-THF) cosolvent system. We used quantum mechanical cluster models to optimize the force field parameters (including the pairwise nonbonded interaction parameters and atomic charges) to be applied in classical simulations. With the reliable force field model, we discussed the importance of including ion correlation explicitly in predicting the molar conductivities via the Onsager formalism and showed that the conventional Nernst–Einstein formula overestimates ionic mobilities due to its intrinsic independent and uncorrelated particle assumption. Further, we investigated the solvation structures and ion pair formations. We concluded that the suitability of the interaction potentials utilized in a classical model for particular systems needs to be assessed not solely by directly comparing the simulated molar conductivities with the measured ones but, more importantly, by using the correct formalism (Onsager) to deduce the simulated result from dynamics trajectories.