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.

The Equivalent Dipole Potential of Water for the Electronic Structure of Threonine

Ab Initio Investigation Of Structure, Stability, Thermal Behavior, Bonding, And Infrared Spectra Of Ionized Water Cluster (H2o)(6)(+) (Vol 145, 154307, 2016)

Ab Initio Investigation of Possible Candidate Structures and Properties of Water Cluster (H2O)7+ Via Particle Swarm Optimization Method

Ab Initio Investigation of the Lower Energy Candidate Structures for (h2o)(5)+ Water Cluster

Application of the flexible and polarisable hydration ion model to study Th(IV) aqueous solutions

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

Delicate balance of hydrogen bonding forces in D-threoninol

Dynamics and structural changes of small water clusters on ionization

Structural and Dynamic Properties of Solvated Hydroxide and Hydronium Ions in Water from Ab Initio Modeling.

Reactivity of Metaphosphate and Thiometaphosphate in Water: A DFT Study

The Structure of Hydrated Electron. Part 1. Magnetic Resonance of Internally Trapping Water Anions: A Density Functional Theory Study

Quantifying the Structure of Water and Hydrated Monovalent Ions by Density Functional Theory-Based Molecular Dynamics.

Theoretical Study on Tautomerism and Proton Transfer of 2,4-Dithiothymine

The electronic structure of liquid water within density functional theory

Infrared Spectra of Small Anionic Water Clusters from Density Functional Theory and Wavefunction Theory Calculations

Density Functional Study on the Structural and Thermodynamic Properties of Aqueous DNA-electrolyte Solution in the Framework of Cell Model

A Systematical Comparison of DFT Methods in Reproducing the Interaction Energies of Halide Series with Protein Moieties.

Ab Initio Study of Different Acid Molecules Interacting with H2O

First-Principles Study Of The Infrared Spectrum In Liquid Water From A Systematically Improved Description Of H-Bond Network

Density Functional Theory Studies on the Real and Apparent Water-Exchange Reaction Kinetics of Al3+ in Aqueous Solution

First-principles study on structural and chemical asymmetry of a biomimetic water-splitting dimanganese complex