AMOEBA Polarizable Force Field for Molecular Dynamics Simulations of Glyme Solvents

J Chem Theory Comput. 2023 Feb 14;19(3):1023-1034. doi: 10.1021/acs.jctc.2c00926. Epub 2023 Jan 24.

Abstract

Classical molecular dynamics (MD) simulations of electrolyte systems are important to gain insight into the atom-scale properties that determine the battery-relevant performance. The recent Tinker-HP software release enables efficient and accurate MD simulations with the AMOEBA polarizable force field. In this work, we developed a procedure to construct a universal AMOEBA model for the solvent family of glymes (glycol methyl ethers), which involves a refinement scheme for valence parameters by fitting the AMOEBA-derived atomic forces to those computed at the DFT level. The refined AMOEBA model provides a good description of both local and nonlocal properties in terms of the spectroscopic response of glyme molecules, as well as the liquid glyme density and dielectric constant. In addition, the complexation energies of alkali and alkaline-earth metal cations with tetraglyme molecules obtained from AMOEBA calculations are in good agreement with DFT results, demonstrating the suitability of the developed AMOEBA model for an accurate simulation of glyme-based battery electrolytes. We also expect the procedure to be transferable to the development of AMOEBA models for other battery electrolyte systems.