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.

Comparison of Experimental and Theoretical Structures of a Transition State Analogue Used for the Induction of Anti-Cocaine Catalytic Antibodies

Comparison of the crystal structure of a transition state analogue that was used to raise catalytic antibodies for the benzoyl ester hydrolysis of cocaine with structures calculated by ab initio, semiempirical, and solvation semiempirical methods reveals that modeling of solvation is crucial for replicating the crystal structure geometry. Both SM3 and SM2 calculations, starting from the crystal structure TSA I, converged on structures similar to the crystal structure. The 3-21G(*)/HF, 6-31G*/HF, PM3, and AM1 calculations converged on structures similar to each other, but these gas-phase structures were significantly extended relative to the condensed phase structures. Two transition states for the hydrolysis of the benzoyl ester of cocaine were located with the SM3 method. The gas phase calculations failed to locate reasonable transition state structures for this reaction. These results imply that accurate modeling of the potential energy surfaces for the hydrolysis of cocaine requires solvation methods.

Prediction of Accurate Anharmonic Experimental Vibrational Frequencies for Water Clusters, (H2O)n, n = 2−5

Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2−5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.

Development of Measurement Methods for Testing of Hydrokinetic Devices to Evaluate the Environmental Effect on Local Substrate

Modifications and upgrades to the hydraulic flume facility in the Environmental Fluid Mechanics and Hydraulics Laboratory (EFM&H) at Bucknell University are described. These changes enable small-scale testing of model marine hydrokinetic(MHK) devices. The design of the experimental platform provides a controlled environment for testing of model MHK devices to determine their effect on localsubstrate. Specifically, the effects being studied are scour and erosion around a cylindrical support structure and deposition of sediment downstream from the device.

Microscopic and Spectroscopic Methods for Probing the Clay Water Interface

Clay minerals have a fundamental importance in many processes in soils and sediments such as the bioavailability of nutrients, water retention, the adsorption of common pollutants, and the formation of an impermeable barrier upon swelling. Many of the properties of clay minerals are due to the unique environment present at the clay mineral/water interface. Traditional techniques such as X-ray diffraction (XRD) and absorption isotherms have provided a wealth of information about this interface but have suffered from limitations. The methods and results presented herein are designed to yield new experimental information about the clay mineral/water interface.A new method of studying the swelling dynamics of clay minerals was developed using in situ atomic force microscopy (AFM). The preliminary results presented here demonstrate that this technique allows one to study individual clay mineral unit layers, explore the natural heterogeneities of samples, and monitor swelling dynamics of clay minerals in real time. Cation exchange experiments were conducted monitoring the swelling change of individual nontronite quasi-crystals as the chemical composition of the surrounding environment was manipulated several times. A proof of concept study has shown that the changes in swelling are from the exchange of interlayer cations and not from the mechanical force of replacing the solution in the fluid cell. A series of attenuated total internal reflection Fourier transform infrared spectroscopy (ATR-FTIR) experiments were performed to gain a better understanding of the organization of water within the interlayer region of two Fe-bearing clay minerals. These experiments made use of the Subtractive Kramers-Kronig (SKK) Transform and the calculation of difference spectra to obtain information about interfacial water hidden within the absorption bands of bulk water. The results indicate that the reduction of structural iron disrupts the organization of water around a strongly hydrated cation such as sodium as the cation transitions from an outer-sphere complex with the mineral surface to an inner-sphere complex. In the case of a less strongly hydrated cation such as potassium, reduction of structural iron actually increases the ordering of water molecules at the mineral surface. These effects were only noticed with the reduction of iron in the tetrahedral sheet close to the basal surface where the increased charge density is localized closer to the cations in the interlayer.