Solvation free energies of anions: from curated reference data to predictive models
Authors:
- Thomas NevolianisTNThomas Nevolianis
- Jonathan ZhengJZJonathan Zheng
- Simon MülLerSMSimon MülLer
- Matthias BaumannMBMatthias Baumann
- Sofja TshepelevitshSTSofja Tshepelevitsh
- Ivari KaljurandIKIvari Kaljurand
- Ivo LeitoILIvo Leito
- Irina SmirnovaISIrina Smirnova
- William GreenWGWilliam Green
- Kai LeonhardKLKai Leonhard
Preprint
- 27 June 2025
- preprint
- Published by American Chemical Society (ACS) in ChemRxiv
Abstract
Predicting the physicochemical properties of ionizable solutes, including solubility and lipophilicity, is of broad significance. Such predictions rely on the accurate determination of solvation free energies for ions. However, the limited availability of high-quality reference data poses a challenge in developing accurate, inexpensive computational prediction methods. In this study, we address both issues of data quality and availability. We present three databases and models related to ionic phenomena: 1 (1) 8,241 pK a datapoints across 8 solvents, (2) 5,536 gas-phase acidities from DLPNO- CCSD(T) QM calculations, and (3) 6,090 solvation free energies of anions across 8 solvents obtained from a thermodynamic cycle. We also report 6,088 solvation free energies of neutral conjugate solutes computed using the COSMO-RS method. The pKa data were obtained from the iBonD database, cleaned, and combined with a separate compilation of trustworthy reference pKa data. Gas-phase acidities were computed for most of the acids present in the pKa corpus. Leveraging these data, we compiled values for solvation free energies of anions. We then trained several graph neural network models, which can be used as an alternative to QM approaches to quickly estimate these properties. The pKa and gas-phase acidity models accept reaction SMILES strings of the acid dissociation as inputs, whereas the solvation energy model accepts the SMILES string of the anion. Our microscopic pKa model achieves good accuracy, with an overall test mean average error of 0.58 units on unseen solutes and 0.59 on the SAMPL7 challenge (the lowest error so far among multi-solvent models). Our gas-phase acidity model had mean absolute errors slightly above 2 kcal/mol when evaluated against experimental data. The anionic solvation free energy model had mean absolute errors of less than 3 kcal/mol in several test evaluations, comparable to (though less reliable than) several widely-used QM-based solvation models. The models and data are free and publicly available.Keywords
All Related Versions
Funding Information
- Deutsche Forschungsgemeinschaft (191948804)
This publication has 85 references indexed in Scilit:
- Cathode Design for Aqueous Rechargeable Multivalent Ion Batteries: Challenges and OpportunitiesAdvanced Functional Materials, 2021
- Systematic microsolvation approach with a cluster‐continuum scheme and conformational samplingJournal of Computational Chemistry, 2020
- Advances and challenges in modeling solvated reaction mechanisms for renewable fuels and chemicalsWIREs Computational Molecular Science, 2019
- On the Basicity of Organic Bases in Different MediaEuropean Journal of Organic Chemistry, 2019
- Software update: the ORCA program system, version 4.0WIREs Computational Molecular Science, 2017
- Materials Science and Materials Chemistry for Large Scale Electrochemical Energy Storage: From Transportation to Electrical GridAdvanced Functional Materials, 2012
- The ORCA program systemWIREs Computational Molecular Science, 2011
- Anchor Points for the Unified Brønsted Acidity Scale: The rCCC Model for the Calculation of Standard Gibbs Energies of Proton Solvation in Eleven Representative Liquid MediaChemistry – A European Journal, 2011
- Measurement of real free energies of transfer of individual ions from water to other solvents with the jet cellJournal of Solution Chemistry, 1991
- Gas‐phase equilibrium affinity scales and chemical ionization mass spectrometryMass Spectrometry Reviews, 1989