Abstract:We introduce a novel approach to compute dissociation free energy and entropy values in simulations that employ a density functional theory description of the acidic moiety and of the solvent. The approach consists of utilizing an alchemical transformation of a weak acid A-COOH into the strong acid B-COOH, which makes it practical to employ alchemical free energy perturbation methods in the context of ab initio molecular dynamics simulations. The present alchemical transformation circumvents the need to tackle changes in the total number of electrons and atoms by replacing the chemical residue responsible for the change in acidity with an easily tunable external effective potential. Our investigation demonstrates that (1) a simple but effective class of external potentials that control acidity changes in the acetic/trifluoroacetic acid series can be achieved by replacing the methyl and trifluoromethyl substituents by screened dipoles. Using this dipole-field/quantum-mechanics (DF/QM) approach one can predict gas-phase geometries, proton dissociation energies, total dipole moments, and water binding energies in good agreement with full-QM values. (2) The resulting alchemical perturbation calculations are stable and well converged and allow one to compute absolute pK(a) values whose accuracy is limited primarily by the exchange-correlation functional employed: H-COOH=2.5+/-0.6 (full-QM calculation), 3.7 (exp); F(3)C-COOH=0.4+/-0.6 (DF/QM calculation), 0.5 (exp); H(3)C-COOH=3.1+/-0.7 (DF/QM calculation), 4.7 (exp); 3) Our DF/QM model predicts that the difference in acidity between H-COOH and H(3)C-COOH is dominated by solvent entropy effects, in excellent agreement with experimental observations. The calculated difference between the dissociation energies of these acids is DeltaDelta(d)U=0.0+/-0.26 kcal/mol while the experimental value is 0.0+/-0.1 kcal/mol.

First-Principle Predictions of Absolute Pka's of Organic Acids in Dimethyl Sulfoxide Solution

What Are the Pka Values of C–H Bonds in Aromatic Heterocyclic Compounds in DMSO?

Quantum-chemical Predictions of Redox Potentials of Organic Anions in Dimethyl Sulfoxide and Reevaluation of Bond Dissociation Enthalpies Measured by the Electrochemical Methods.

What are the pKa values of organophosphorus compounds?

Quantum-Chemical Predictions Of Pk(A)'S Of Thiols In Dmso

Theoretical Study on the Acidities of Chiral Phosphoric Acids in Dimethyl Sulfoxide: Hints for Organocatalysis

First-principle predictions of basicity of organic amines and phosphines in acetonitrile

Computation of Pk(A) Values of Substituted Aniline Radical Cations in Dimethylsulfoxide Solution

A Systematic Theoretical Study on the Acidities for Cations of Ionic Liquids in Dimethyl Sulfoxide.

Toward Prediction of the Chemistry in Ionic Liquids: an Accurate Computation of Absolute Pk(a) Values of Benzoic Acids and Benzenethiols.

pKa Prediction in Non-Aqueous Solvents

Computational Study on the Acidic Constants of Chiral Brønsted Acids in Dimethyl Sulfoxide.

Holistic Prediction of the Pka in Diverse Solvents Based on a Machine‐Learning Approach

First-Principle Calculation Of Equilibrium Cesium Ion-Pair Acidities In Tetrahydrofuran

Holistic Prediction of Pka in Diverse Solvents Based on Machine Learning Approach

Unexpected Solvation-Stabilisation of Ions in a Protic Ionic Liquid: Insights Disclosed by a Bond Energetic Study.

Quantum-chemical Predictions of Absolute Standard Redox Potentials of Diverse Organic Molecules and Free Radicals in Acetonitrile.

Simulating solvation and acidity in complex mixtures with first-principles accuracy: the case of CH$_3$SO$_3$H and H$_2$O$_2$ in phenol

Ab Initio Study of Different Acid Molecules Interacting with H2O

An Acidity Scale of 1,3-Dialkylimidazolium Salts in Dimethyl Sulfoxide Solution.

Combining ab initio quantum mechanics with a dipole-field model to describe acid dissociation reactions in water: first-principles free energy and entropy calculations