These are preliminary reports that have not been peer-reviewed. They should not be regarded as conclusive, guide clinical practice/health-related behavior, or be reported in news media as established information. For more information, please see our FAQs.
Preprints are manuscripts made publicly available before they have been submitted for formal peer review and publication. They might contain new research findings or data. Preprints can be a draft or final version of an author's research but must not have been accepted for publication at the time of submission.
submitted on 07.09.2018 and posted on 10.09.2018by Marie-Madeleine Walz, Mohammad Ghahremanpour, Paul J. van Maaren, David van der Spoel
A longstanding goal of computational chemistry is to predict the state of materials in all phases with a single model. This is particularly relevant for materials that are difficult or dangerous to handle or compounds that have not yet been created. Progress towards this goal has been limited as most work has concentrated on just one phase, often determined by particular applications. In the framework of the development of the Alexandria force field we present here new polarizable force fields for alkali halides with Gaussian charge distributions for molecular dynamics simulations. We explore different descriptions of the Van der Waals interaction, like the commonly applied 12-6 Lennard-Jones (LJ), and compare it to \softer" ones, such as 8-6 LJ, Buckingham and a modified Buckingham potential. Our results for physico-chemical properties of the gas, liquid and solid phase of alkali halides, are compared to experimental data and calculations with reference polarizable and non-polarizable force fields. The new polarizable force field that employs a modified Buckingham potential predicts the tested properties for gas, liquid and solid phases with a very good accuracy. In contrast to reference force fields, this model reproduces the correct crystal structures for all alkali halides at low and high temperature. Seeing that experiments with molten salts may be tedious due to high temperatures and their corrosive nature, the models presented here can contribute significantly to our understanding of alkali halides in general and melts in particular.