Logo: to the web site of Uppsala University

Large scale quantum calculations for molar enthalpy of formation (Δ_fH⁰), standard entropy (S⁰), and heat capacity (C_V) are presented. A large data set may help to evaluate quantum thermochemistry tools in order to uncover possible hidden shortcomings and also to find experimental data that might need to be reinvestigated, indeed we list and annotate approximately 200 problematic thermochemistry measurements. Quantum methods systematically underestimate S⁰ for flexible molecules in the gas phase if only a single (minimum energy) conformation is taken into account. This problem can be tackled in principle by performing thermochemistry calculations for all stable conformations [Zheng et al., Phys. Chem. Chem. Phys. 13, 10885–10907 (2011)], but this is not practical for large molecules. We observe that the deviation of composite quantum thermochemistry recipes from experimental S⁰ corresponds roughly to the Boltzmann equation (S = RlnΩ), where R is the gas constant and Ω the number of possible conformations. This allows an empirical correction of the calculated entropy for molecules with multiple conformations. With the correction we find an RMSD from experiment of ≈13 J/mol K for 1273 compounds. This paper also provides predictions of Δ_fH⁰, S⁰, and C_V for well over 700 compounds for which no experimental data could be found in the literature. Finally, in order to facilitate the analysis of thermodynamics properties by others we have implemented a new tool obthermo in the OpenBabel program suite [O’Boyle et al., J. Cheminf. 3, 33 (2011)] including a table of reference atomization energy values for popular thermochemistry methods.

Data quality as well as library size are crucial issues for force field development. In order to predict molecular properties in a large chemical space, the foundation to build force fields on needs to encompass a large variety of chemical compounds. The tabulated molecular physicochemical properties also need to be accurate. Due to the limited transparency in data used for development of existing force fields it is hard to establish data quality and reusability is low. This paper presents the Alexandria library as an open and freely accessible database of optimized molecular geometries, frequencies, electrostatic moments up to the hexadecupole, electrostatic potential, polarizabilities, and thermochemistry, obtained from quantum chemistry calculations for 2704 compounds. Values are tabulated and where available compared to experimental data. This library can assist systematic development and training of empirical force fields for a broad range of molecules.

Gas-phase electric properties of molecules can be computed routinely using wave function methods or density functional theory (DFT). However, these methods remain computationally expensive for high-throughput screening of the vast chemical space of virtual compounds. Therefore, empirical force fields are a more practical choice in many cases, particularly since force field methods allow one to routinely predict the physicochemical properties in the condensed phases. This work presents Drude polarizable models, to increase the physical realism in empirical force fields, where the core particle is treated as a point charge and the Drude particle is treated either as a 1s-Gaussian or a ns-Slater (n = 1, 2, 3) charge density. Systematic parametrization to large high-quality quantum chemistry data obtained from the open access Alexandria Library (https://doi.org/10.5281/zenodo.1004711) ensures the transferability of these parameters. The dipole moments and isotropic polarizabilities of the isolated molecules predicted by the proposed Drude models are in agreement with experiment with accuracy similar to DFT calculations at the B3LYP/aug-cc-pVTZ level of theory. The results show that the inclusion of explicit polarization into the models reduces the root-mean-square deviation with respect to DFT calculations of the predicted dipole moments of 152 dimers and clusters by more than 50%. Finally, we show that the accuracy of the electrostatic interaction energy of the water dimers can be improved systematically by the introduction of polarizable smeared charges as a model for charge penetration.

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 toward 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 the 8–6 LJ, Buckingham, and a modified Buckingham potential. Our results for physicochemical properties of the gas, liquid, and solid phases of alkali halides are compared to experimental data and calculations with reference polarizable and nonpolarizable 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.

Spectroscopic analysis of compounds is typically combined with density functional theory, for instance, for assigning vibrational frequencies, limiting application to relatively small compounds. Accurate classical force fields could, in principle, complement these quantum-chemical tools. A relatively simple way to validate vibrational frequencies is by computing thermochemical properties. We present such a validation for over 1800 small molecules using the harmonic approximation. Two popular empirical force fields (GAFF and CGenFF) are compared to experimental data and results from Gaussian-4 quantum-chemical calculations. Frequency scaling factors of 1.035 (CGenFF) and 1.018 (GAFF) are derived from the zero-point energies. The force field calculations have larger deviation from experiment than the G4 method for standard entropy, but for heat capacity the results are comparable. For internal thermal energy and zero-point energy the deviations from G4 are relatively small. The work suggests that with some tuning force fields could indeed complement DFT in spectroscopical applications.