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.

An Approach to the Estimation of Solute Hydrogen-Bond Donor Acidity

Proposed Hydrogen-Bonding Index of Donor or Acceptor Reflecting Its Intrinsic Contribution to Hydrogen-Bonding Strength

A 19F Solvatomagnetic Comparison Method for the Determination of the α1 Scale of Solvent Hydrogen‐Bond Donation and Its Application to Deep Eutectic Solvents

Determination of the Hydrogen-Bond Basicity Descriptor by Reversed-Phase Liquid Chromatography

Density, Acidity and Theoretical Calculation of Hydrophobic Deep Eutectic Solvents Formed by Amide and Lactic Acid

An Approach to Computing Solvent Reorganization Energy.

pKa Prediction in Non-Aqueous Solvents

Experimental quantification of electrostatics in X-H···π hydrogen bonds

Aerosol pH and Ion Activities of HSO 4 - and SO 4 2- in Supersaturated Single Droplets

H-Bond Acceptor Parameters for Anions

Computation of standard equilibrium acidity of C–H acids in ionic media: shedding light on predicting changes of chemical behavior by switching solvent system from molecular to ionic

Surface tensions and surface potentials of acid solutions

Quantification and Application of Lewis Acidity and Basicity of Benzene and Its Derivates

Estimation on the Individual Hydrogen-Bond Strength in Molecules with Multiple Hydrogen Bonds.

Effect of Solvent Dielectric Constant and Acidity on the OH Vibration Frequency in Hydrogen-Bonded Complexes of Fluorinated Ethanols

Aerosol pH and Ion Activities of HSO4- and SO42- in Supersaturated Single Droplets

Determination of Pka Values of C-H Bonds in Polar Fluorinated Arenes Referred to a New CF3SO2-Substituted Anchor Compound.

Accurate Description of Solvent-Exposed Salt Bridges with a Non-polarizable Force Field Incorporating Solvent Effects

Heterolytic (R-H)+· Bond Dissociation Free Energies in Solution Derived from Electrochemically-monitored Acid-Base Titration (EABT)

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

Real single ion solvation free energies with quantum mechanical simulation