A polarizable valence electron density based force field for high-energy interactions between atoms and molecules

High-accuracy molecular force field models suited for hot gases and plasmas are not as abundant as those geared toward ambient pressure and temperature conditions. Here, we present an improved version of our previous electron-density based force field model that can now account for polarization effects by adjusting the atomic valence electron contributions to match ab initio calculated Mulliken partial charges. Using a slightly modified version of the Hohenberg-Kohn theorem, we also include an improved theoretical formulation of our model when applied to systems with degenerate ground states. We present two variants of our polarizable model, fitted from ab initio reference data calculated at CCSD(T)/cc-pVTZ and CCSD(T)/CEP-31G levels of theory, that both accurately model water dimer interaction energies. Further improvements include the additional interaction components with fictitious non-spherically symmetric, yet atom-centered, electron densities and fitting the exchange and correlation coefficients against analytical expressions. The latter removes all unphysical oscillations that are observed in the previous non-polarizable variant of our force field.