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 Calculation Of Equilibrium Cesium Ion-Pair Acidities In Tetrahydrofuran

DFT Calculation and Experimental Validation on the Interactions of Bis(trifluoromethylsulfonyl)imide and Hexafluorophosphate with Cesium

(Thio)urea Organocatalyst Equilibrium Acidities in Dmso

In silico calculation of acidity constants of carbonic acid conformers.

Thermodynamic Properties for the Aqueous Solutions of Cesium Borates at 298.15 K and 101 Kpa: Experimental and Correlation by Pitzer Ion-Interaction Model

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

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

Quantum Mechanical Prediction of Dissociation Constants for Thiazol-2-imine Derivatives

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

Ab Initio Computational Modeling on the Tautomerism of Monochalcogenocarboxylic Acids CH3C(O)XH (X = S, Se, and Te) in the Polar and Aprotic Solution

Local Acid Strength of Solutions and Its Quantitative Evaluation Using Excess Infrared Nitrile Probes

Unexpected Strong Acidity Enhancing the Effect in Protic Ionic Liquids Quantified by Equilibrium Acidity Studies: A Crucial Role of Cation Structures on Dictating the Solvation Properties.

A Theoretical Study of Electrical and Electrochemical Properties of Dicyanomethylene Derivatives of Squaric Acid

Thermodynamics of complexation of cesium ions with dibenzo-18-crown-6 in hydrophobic ionic liquids

Conformational distributions of tetrahydro-2-turoic acid in water at different pH values by their IR and vibrational circular dichroism spectra

Thermodynamics of cesium complexes formation with 18-crown-6 in ionic liquids

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

Theoretical Predictions of the Structure, Gas-Phase Acidity, and Aromaticity of Tetrathiosquaric Acid

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

Interaction of the cesium cation with calix[4]arene-bis(t-octylbenzo-18-crown-6): Extraction and DFT study