Proposed carbon dioxide sequestration scenarios in sedimentary reservoirs require investigation into the interactions between supercritical carbon dioxide, brines, and the mineral phases found in the basin and overlying caprock. Molecular simulations can help to understand the partitioning of metal cations between aqueous solutions and supercritical carbon dioxide where limited experimental dat...

We present a new methodology for computing solvation free energy, which is based upon the reference interaction site model (RISM)/hypernetted chain (HNC) solvation free energy expression, but which substitutes radial distribution functions taken from simulations for those calculated by simultaneous solution of the RISM and HNC equations. Consequently, solvation free energy can be obtained from ...

Three pulse stimulated photon echo peak shift ͑3PEPS͒ measurements were used to probe the solvation of a quadrupolar solute in three room temperature nondipolar solvents; benzene, CCl 4 , and CS 2 , and the results were compared with those for two polar solvents, methanol and acetonitrile, and one weakly polar solvent, toluene. Our data reveal three distinct solvent dynamical time scales; a sub-1...

A novel method for fast and accurate evaluation of the generalized Born radii in macromolecular solvation electrostatics calculations is proposed, based on the solvent accessibility of the first two solvation layers around an atom. The reverse generalized Born radii calculated by the method have correlation coefficient of 98.7% and RMSD of 0.031 A(-1) with the values obtained using a precise bu...

A local continuum solvation theory, exactly treating electrostatic matching conditions on the boundary of a cavity occupied by a solute particle, is extended to cover time-dependent solvation phenomena. The corresponding integral equation is solved with a complex-valued frequency-dependent dielectric function «~v!, resulting in a complex-valued v-dependent reaction field. The inverse Fourier tr...

A new solvation model that combines discrete and continuum descriptions of the solvent has been developed. The discrete solvent molecules are represented by effective fragment potentials (EFP), while the continuum is represented by the Onsager model. This (EFP+Onsager) model has been applied to the relative stabilities of the neutral and zwitterionic forms of glycine. Other supermolecule-contin...

This work extends our earlier two-domain formulation of a differential geometry based multiscale paradigm into a multidomain theory, which endows us the ability to simultaneously accommodate multiphysical descriptions of aqueous chemical, physical and biological systems, such as fuel cells, solar cells, nanofluidics, ion channels, viruses, RNA polymerases, molecular motors and large macromolecu...

We present a new set of geometry-based functional forms for parametrizing effective Coulomb radii and atomic surface tensions of organic solutes in water. In particular, the radii and surface tensions depend in some cases on distances to nearby atoms. Combining the surface tensions with electrostatic effects included in a Fock operator by the generalized Born model enables one to calculate free...

Thermodynamic measurements of the solvation of salts and electrolytes are relatively straightforward, but it is not possible to separate total solvation free energies into distinct cation and anion contributions without reference to an additional extrathermodynamic assumption. The present work attempts to resolve this difficulty using molecular dynamics simulations with the AMOEBA polarizable f...

Membrane channel proteins control the diffusion of ions across biological membranes. They are closely related to the processes of various organizational mechanisms, such as: cardiac impulse, muscle contraction and hormone secretion. Introducing a membrane region into implicit solvation models extends the ability of the Poisson-Boltzmann (PB) equation to handle membrane proteins. The use of late...