Theoretical Calculations of Acid Dissociation Constants: A Review Article

Acid dissociation constants, or pKa values, are essential for understanding many fundamental reactions in chemistry. These values reveal the deprotonation state of a molecule in a particular solvent. There is great interest in using theoretical methods to calculate the pKa values for many different types of molecules. These include molecules that have not been synthesized, those for which experimental pKa determinations are difficult, and for larger molecules where the local environment changes the usual pKa values, such as for certain amino acids that are part of a larger polypeptide chain. Chemical accuracy in pKa calculations is difficult to achieve, because an error of 1.36 kcal/mol in the change of free energy for deprotonation in solvent results in an error of 1 pKa unit. In this review the most valuable methods for determining accurate pKa values in aqueous solution are presented for educators interested in explaining or using these methods for their students.

Experimentation with different thermodynamic cycles used for pKa calculations on carboxylic acids using Complete Basis Set and Gaussian-n Models combined with CPCM Continuum Solvation Methods

Complete basis set and Gaussian-n methods were combined with Barone and Cossi's implementation of the polarizable conductor model (CPCM) continuum solvation methods to calculate pKa values for six carboxylic acids. Four different thermodynamic cycles were considered in this work. An experimental value of −264.61 kcal/mol for the free energy of solvation of H+, ΔGs(H+), was combined with a value for Ggas(H+) of −6.28 kcal/mol, to calculate pKa values with cycle 1. The complete basis set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. Thermodynamic cycles that include an explicit water in the cycle are not accurate when the free energy of solvation of a water molecule is used, but appear to become accurate when the experimental free energy of vaporization of water is used. This apparent improvement is an artifact of the standard state used in the calculation. Geometry relaxation in solution does not improve the results when using these later cycles. The use of cycle 1 and the complete basis set models combined with the CPCM solvation methods yielded pKa values accurate to less than half a pKa unit. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001

Accurate pKa Calculations for Carboxylic Acids Using Complete Basis Set and Gaussian-n Models Combined with CPCM Continuum Solvation Methods”

Complete Basis Set and Gaussian-n methods were combined with CPCM continuum solvation methods to calculate pKa values for six carboxylic acids. An experimental value of −264.61 kcal/mol for the free energy of solvation of H+, ΔGs(H+), was combined with a value for Ggas(H+) of −6.28 kcal/mol to calculate pKa values with Cycle 1. The Complete Basis Set gas-phase methods used to calculate gas-phase free energies are very accurate, with mean unsigned errors of 0.3 kcal/mol and standard deviations of 0.4 kcal/mol. The CPCM solvation calculations used to calculate condensed-phase free energies are slightly less accurate than the gas-phase models, and the best method has a mean unsigned error and standard deviation of 0.4 and 0.5 kcal/mol, respectively. The use of Cycle 1 and the Complete Basis Set models combined with the CPCM solvation methods yielded pKa values accurate to less than half a pKa unit.

Accurate relative pKa calculations for carboxylic acids using complete basis set and Gaussian-n models combined with continuum solvation methods

The complete basis set methods CBS-4, CBS-QB3, and CBS-APNO, and the Gaussian methods G2 and G3 were used to calculate the gas phase energy differences between six different carboxylic acids and their respective anions. Two different continuum methods, SM5.42R and CPCM, were used to calculate the free energy differences of solvation for the acids and their anions. Relative pKa values were calculated for each acid using one of the acids as a reference point. The CBS-QB3 and CBS-APNO gas phase calculations, combined with the CPCM/HF/6-31+G(d)//HF/6-31G(d) or CPCM/HF/6-31+G(d)//HF/6-31+G(d) continuum solvation calculations on the lowest energy gas phase conformer, and with the conformationally averaged values, give results accurate to ½ pKa unit. © 2001 American Institute of Physics.

AM1 and PM3 Calculations of the Potential Energy Surfaces for CH2OH Reactions with NO and NO2

The AM1 and PM3 molecular orbital methods have been utilized to investigate the reactions of CH20H with NO and NO2 PM3 and AM1 calculated heats of formation differ from experimental values by 8.6 and 18.8 kcal mol-', respectively. The dominant reaction of CH20H with NO is predicted to produce the adduct HOCH2N0, supporting the hypothesis of Pagsberg, Munk, Anastasi, and Simpson. Calculated activation energies for the NO2 system predict the formation of the adducts HOCH2N02 and HOCH20N0. In addition, the PM3 calculations predict that the abstraction reaction producing CH20 and HN02 is more likely than one producing CH20 and HONO from reactions of CH20H with NO2.

A Metric to Evaluate and Synthesize Distributed Compliant Mechanisms

Compliant mechanisms with evenly distributed stresses have better load-bearing ability and larger range of motion than mechanisms with compliance and stresses lumped at flexural hinges. In this paper, we present a metric to quantify how uniformly the strain energy of deformation and thus the stresses are distributed throughout the mechanism topology. The resulting metric is used to optimize cross-sections of conceptual compliant topologies leading to designs with maximal stress distribution. This optimization framework is demonstrated for both single-port mechanisms and single-input single-output mechanisms. It is observed that the optimized designs have lower stresses than their nonoptimized counterparts, which implies an ability for single-port mechanisms to store larger strain energy, and single-input single-output mechanisms to perform larger output work before failure.

From Structure to Function: Characterization of Cu(I) Adducts in Leveler Additives by DFT Calculations

We present the first molecular model of the coordination complex formed by Cu(I) and imidazole-epichlorohydrin polymers. Our calculations show that the Cu(I) ion has linear coordination and the whole complex has neutral charge. Our model suggests salt couple pairing as the driving force for the formation of the surface-confined precipitation, which is crucial to obtain flat surfaces in industrial copper deposition processes, required for mass fabrication of state-of-the-art electronic and memory devices.

Vibronic Spectra of Jet-Cooled 2-Aminopurine center dot H2O Clusters Studied by UV Resonant Two-Photon Ionization Spectroscopy and Quantum Chemical Calculations

For understanding the major- and minor-groove hydration patterns of DNAs and RNAs, it is important to understand the local solvation of individual nucleobases at the molecular level. We have investigated the 2-aminopurine center dot H2O. monohydrate by two-color resonant two-photon ionization and UV/UV hole-burning spectroscopies, which reveal two isomers, denoted A and B. The electronic spectral shift delta nu of the S-1 <- S-0 transition relative to bare 9H-2-aminopurine (9H-2AP) is small for isomer A (-70 cm(-1)), while that of isomer B is much larger (delta nu = 889 cm(-1)). B3LYP geometry optimizations with the TZVP basis set predict four cluster isomers, of which three are doubly H-bonded, with H2O acting as an acceptor to a N-H or -NH2 group and as a donor to either of the pyrimidine N sites. The "sugar-edge" isomer A is calculated to be the most stable form with binding energy D-e = 56.4 kJ/mol. Isomers B and C are H-bonded between the -NH2 group and pyrimidine moieties and are 2.5 and 6.9 kJ/mol less stable, respectively. Time-dependent (TD) B3LYP/TZVP calculations predict the adiabatic energies of the lowest (1)pi pi* states of A and B in excellent agreement with the observed 0(0)(0) bands; also, the relative intensities of the A and B origin bands agree well with the calculated S-0 state relative energies. This allows unequivocal identification of the isomers. The R2PI spectra of 9H-2AP and of isomer A exhibit intense low-frequency out-of-plane overtone and combination bands, which is interpreted as a coupling of the optically excited (1)pi pi* state to the lower-lying (1)n pi* dark state. In contrast, these overtone and combination bands are much weaker for isomer B, implying that the (1)pi pi* state of B is planar and decoupled from the (1)n pi* state. These observations agree with the calculations, which predict the (1)n pi* above the (1)pi pi* state for isomer B but below the (1)pi pi* for both 9H-2AP and isomer A.