The Smoluchowski approach to diffusion-controlled reactions is generalized to interacting substrate particles by including the osmotic pressure and hydrodynamic interactions of the nonideal particles in the Smoluchoswki equation within a local-density approximation. By solving the strictly linearized equation for the time-independent case with absorbing boundary conditions, we present an analytic expression for the diffusion-limited steady-state rate constant for small substrate concentrations in terms of an effective second virial coefficient B-2(*). Comparisons to Brownian dynamics simulations excluding hydrodynamic interactions show excellent agreement up to bulk number densities of B(2)(*)rho(0)less than or similar to 0.4 for hard sphere and repulsive Yukawa-like interactions between the substrates. Our study provides an alternative way to determine the second virial coefficient of interacting macromolecules experimentally by measuring their steady-state rate constant in diffusion-controlled reactions at low densities. (c) 2005 American Institute of Physics.

The permeation of hydrophobic, cylindrical nanopores by water molecules and ions is investigated under equilibrium and out-of-equilibrium conditions by extensive molecular-dynamics simulations. Neglecting the chemical structure of the confining pore surface, we focus on the effects of pore radius and electric field on permeation. The simulations confirm the intermittent filling of the pore by water, reported earlier under equilibrium conditions for pore radii larger than a critical radius R-c. Below this radius, water can still permeate the pore under the action of a strong electric field generated by an ion concentration imbalance at both ends of the pore embedded in a structureless membrane. The water driven into the channel undergoes considerable electrostriction characterized by a mean density up to twice the bulk density and by a dramatic drop in dielectric permittivity which can be traced back to a considerable distortion of the hydrogen-bond network inside the pore. The free-energy barrier to ion permeation is estimated by a variant of umbrella sampling for Na+, K+, Ca2+, and Cl- ions, and correlates well with known solvation free energies in bulk water. Starting from an initial imbalance in ion concentration, equilibrium is gradually restored by successive ion passages through the water-filled pore. At each passage the electric field across the pore drops, reducing the initial electrostriction, until the pore, of radius less than R-c, closes to water and hence to ion transport, thus providing a possible mechanism for voltage-dependent gating of hydrophobic pores. (C) 2005 American Institute of Physics.

An implicit solvent model is presented that couples hydrophobic, dispersion, and electrostatic solvation energies by minimizing the system Gibbs free energy with respect to the solvent volume exclusion function. The solvent accessible surface is the output of the theory. The method is illustrated with the solvation of simple solutes on different length scales and captures the sensitivity of hydration to the particular form of the solute-solvent interactions in agreement with recent computer simulations.

Recent studies on the solvation of atomistic and nanoscale solutes indicate that a strong coupling exists between the hydrophobic, dispersion, and electrostatic contributions to the solvation free energy, a facet not considered in current implicit solvent models. We suggest a theoretical formalism which accounts for coupling by minimizing the Gibbs free energy of the solvent with respect to a solvent volume exclusion function. The resulting differential equation is similar to the Laplace-Young equation for the geometrical description of capillary interfaces but is extended to microscopic scales by explicitly considering curvature corrections as well as dispersion and electrostatic contributions. Unlike existing implicit solvent approaches, the solvent accessible surface is an output of our model. The presented formalism is illustrated on spherically or cylindrically symmetrical systems of neutral or charged solutes on different length scales. The results are in agreement with computer simulations and, most importantly, demonstrate that our method captures the strong sensitivity of solvent expulsion and dewetting to the particular form of the solvent-solute interactions.