Structural and thermodynamic properties of ionic fluids are related to those of a simpler "mimic" system with short ranged intermolecular interactions in a spatially varying effective field by use of Local Molecular Field (LMF) theory, already successfully applied to nonuniform simple fluids. By consistently using the LMF approximation to describe only the slowly varying part of the Coulomb interaction, which we view as arising from a rigid Gaussian charge distribution with an appropriately chosen width σ, exceptionally accurate results can be found. In this paper we study a uniform system of charged hard spheres in a uniform neutralizing background, where these ideas can be presented in their simplest form. At low densities the LMF theory reduces to a generalized version of the Poisson-Boltzmann approximation, but the predicted structure factor satisfies the exact StillingerLovett moment conditions, and with optimal choice of σ the lowest order approximation remains accurate for much stronger couplings. At high density and strong couplings the pair correlation function in the uniform mimic system with short ranged interactions is very similar to that of the full ionic system. A simple analytic formula can then describe the difference in internal energy between the ionic system and the associated mimic system.

Strong, short-ranged positional correlations involving counterions can induce a net attractive force between negatively charged strands of DNA and lead to the formation of ion pairs in dilute ionic solutions. However, the long range of the Coulomb interactions impedes the development of a simple local picture. We address this general problem by mapping the properties of a nonuniform system with Coulomb interactions onto those of a simpler system with short-ranged intermolecular interactions in an effective external field that accounts for the averaged effects of appropriately chosen long-ranged and slowly varying components of the Coulomb interactions. The remaining short-ranged components combine with the other molecular core interactions and strongly affect pair correlations in dense or strongly coupled systems. We show that pair correlation functions in the effective short-ranged system closely resemble those in the uniform primitive model of ionic solutions and illustrate the formation of ion pairs and clusters at low densities. The theory accurately describes detailed features of the effective attraction between two equally charged walls at strong coupling and intermediate separations of the walls. Analytical results for the minimal coupling strength needed to get any attraction and for the separation at which the attractive force is a maximum are presented.effective short-ranged model ͉ ion pairing ͉ mean field theory ͉ Poission-Boltzmann S trong Coulomb interactions in crowded, nonuniform environments have important experimental consequences in a wide variety of biophysical applications ranging from DNA packaging in viruses to transport in ion channels (1-4). These interactions present major challenges to theory and computer simulations not only because of their characteristic long range but also because they can be very strong at short distances. Here, we present a local molecular field (LMF) theory (5) that averages over particular long-ranged and slowly varying components of the Coulomb interactions (6) while still maintaining an accurate description of the short-ranged components. Our model provides a general and physically suggestive theory for strongly coupled Coulomb systems and reduces exactly to the classical Poisson-Boltzmann (PB) approximation for dilute, weakly coupled systems.We consider a general starting point where a molecule of species i, described by a rigid body frame with center at r i , interacts with an external field, fi (r i ), that comes from fixed charged solutes, or walls, or particular fixed molecules of a mobile species, as illustrated below. The subscript f indicates the source of the field, which we treat as a special fixed species f. The interaction between a pair of molecules of species i and j is assumed to have the general form w ij (r ij ) ϭ w s,ij (r ij ) ϩ w q,ij (r ij ), where r ij ϵ r j Ϫ r i . The w s,ij (r ij ) denote general (repulsive core and other), short-ranged intermolecular interactions. There are angular coordinates expressing orientations of the body fram...

Effective attraction between like-charged walls mediated by counterions is studied using local molecular field (LMF) theory. Monte Carlo simulations of the "mimic system" given by LMF theory, with short-ranged "Coulomb core" interactions in an effective single particle potential incorporating a mean-field average of the long-ranged Coulomb interactions, provide a direct test of the theory, and are in excellent agreement with more complex simulations of the full Coulomb system by Moreira and Netz [Eur. Phys. J. E 8, 33 (2002)]. A simple, generally applicable criterion to determine the consistency parameter sigma(min) needed for accurate use of the LMF theory is presented.

Hydration free energy of hydrophobic solutes studied by a reference interaction site model with a repulsive bridge correction and a thermodynamic perturbation method A simple weighted-density-functional method: Test and its application to hard sphere fluid in spherical cavity This paper determines the excess free energy associated with the formation of a spherical cavity in a hard sphere fluid. The solvation free energy can be calculated by integration of the structural changes induced by inserting the cavity using a number of different exact thermodynamic pathways. We consider three such pathways, including a new density route derived here. Structural information about the nonuniform hard sphere fluid in the presence of a general external field is given by the recently developed hydrostatic linear response ͑HLR͒ integral equation. Use of the HLR results in the different pathways gives a generally accurate determination of the solvation free energy for cavities over a wide range of sizes, from zero to infinity. Results for a related method, the Gaussian field model, are also discussed.

Percus showed that approximate theories for the structure of nonuniform hard sphere fluids can be generated by linear truncations of functional expansions of the nonuniform density rho(r) about that of an appropriately chosen uniform system. We consider the most general such truncation, which we refer to as the shifted linear response (SLR) equation, where the density response rho(r) to an external field phi(r) is expanded to linear order at each r about a different uniform system with a locally shifted chemical potential. Special cases include the Percus-Yevick (PY) approximation for nonuniform fluids, with no shift of the chemical potential, and the hydrostatic linear response (HLR) equation, where the chemical potential is shifted by the local value of phi(r). The HLR equation gives exact results for very slowly varying phi(r) and reduces to the PY approximation for hard core phi(r), where generally accurate results are found. We show that a truncated expansion about the bulk density (the PY approximation) also gives exact results for localized fields that are nonzero only in a "tiny" region whose volume V(phi) can accommodate at most one particle. The SLR equation can also exactly describe a limit where the fluid is confined by hard walls to a very narrow slit. This limit can be related to the localized field limit by a simple shift of the chemical potential, leading to an expansion about the ideal gas. We then try to develop a systematic way of choosing an optimal local shift in the SLR equation for general phi(r) by requiring that the predicted rho(r) is insensitive to small variations about the appropriate local shift, a property that the exact expansion to all orders would obey. The resulting insensitivity criterion (IC) gives a theory that reduces to the HLR equation for slowly varying phi(r) and is much more accurate than HLR both for very narrow slits, where the IC agrees with exact results, and for fields confined to "tiny" regions, where the IC gives very accurate (but not exact) results. However, the IC is significantly less accurate than the PY and HLR equations for single hard core fields. Only a small change in the predicted reference density is needed to correct this remaining limit.