Progress toward chemical accuracy in the computer simulation of condensed phase reactions.

We describe a procedure for the generation of chemically accurate computer-simulation models to study chemical reactions in the condensed phase. The process involves (i) the use of a coupled semiempirical quantum and classical molecular mechanics method to represent solutes and solvent, respectively; (ii) the optimization of semiempirical quantum mechanics (QM) parameters to produce a computationally efficient and chemically accurate QM model; (iii) the calibration of a quantum/classical microsolvation model using ab initio quantum theory; and (iv) the use of statistical mechanical principles and methods to simulate, on massively parallel computers, the thermodynamic properties of chemical reactions in aqueous solution. The utility of this process is demonstrated by the calculation of the enthalpy of reaction in vacuum and free energy change in aqueous solution for a proton transfer involving methanol, methoxide, imidazole, and imidazolium, which are functional groups involved with proton transfers in many biochemical systems. An optimized semiempirical QM model is produced, which results in the calculation of heats of formation of the above chemical species to within 1.0 kcal/mol (1 kcal = 4.18 kJ) of experimental values. The use of the calibrated QM and microsolvation QM/MM (molecular mechanics) models for the simulation of a proton transfer in aqueous solution gives a calculated free energy that is within 1.0 kcal/mol (12.2 calculated vs. 12.8 experimental) of a value estimated from experimental pKa values of the reacting species.

Computation of Surface Electrical Potentials of Plant Cell Membranes : Correspondence to Published Zeta Potentials from Diverse Plant Sources

A Gouy-Chapman-Stern model has been developed for the computation of surface electrical potential (ψ0) of plant cell membranes in response to ionic solutes. The present model is a modification of an earlier version developed to compute the sorption of ions by wheat (Triticum aestivum L. cv Scout 66) root plasma membranes. A single set of model parameters generates values for ψ0 that correlate highly with published ζ potentials of protoplasts and plasma membrane vesicles from diverse plant sources. The model assumes ion binding to a negatively charged site (R− = 0.3074 μmol m−2) and to a neutral site (P0 = 2.4 μmol m−2) according to the reactions R− + IΖ ⇌ RIΖ−1 and P0 + IΖ ⇌ PIΖ, where IΖ represents an ion of charge Ζ. Binding constants for the negative site are 21,500 m−1 for H+, 20,000 m−1 for Al3+, 2,200 m−1 for La3+, 30 m−1 for Ca2+ and Mg2+, and 1 m−1 for Na+ and K+. Binding constants for the neutral site are 1/180 the value for binding to the negative site. Ion activities at the membrane surface, computed on the basis of ψ0, appear to determine many aspects of plant-mineral interactions, including mineral nutrition and the induction and alleviation of mineral toxicities, according to previous and ongoing studies. A computer program with instructions for the computation of ψ0, ion binding, ion concentrations, and ion activities at membrane surfaces may be requested from the authors.

Rock friction and its implications for earthquake prediction examined via models of Parkfield earthquakes.

The friction of rocks in the laboratory is a function of time, velocity of sliding, and displacement. Although the processes responsible for these dependencies are unknown, constitutive equations have been developed that do a reasonable job of describing the laboratory behavior. These constitutive laws have been used to create a model of earthquakes at Parkfield, CA, by using boundary conditions appropriate for the section of the fault that slips in magnitude 6 earthquakes every 20-30 years. The behavior of this model prior to the earthquakes is investigated to determine whether or not the model earthquakes could be predicted in the real world by using realistic instruments and instrument locations. Premonitory slip does occur in the model, but it is relatively restricted in time and space and detecting it from the surface may be difficult. The magnitude of the strain rate at the earth's surface due to this accelerating slip seems lower than the detectability limit of instruments in the presence of earth noise. Although not specifically modeled, microseismicity related to the accelerating creep and to creep events in the model should be detectable. In fact the logarithm of the moment rate on the hypocentral cell of the fault due to slip increases linearly with minus the logarithm of the time to the earthquake. This could conceivably be used to determine when the earthquake was going to occur. An unresolved question is whether this pattern of accelerating slip could be recognized from the microseismicity, given the discrete nature of seismic events. Nevertheless, the model results suggest that the most likely solution to earthquake prediction is to look for a pattern of acceleration in microseismicity and thereby identify the microearthquakes as foreshocks.