The structural properties of a two-Yukawa fluid: Simulation and analytical results

Standard Monte Carlo simulations are carried out to assess the accuracy of theoretical predictions for the structural properties of a model fluid interacting through a hard-core two-Yukawa potential composed of a short-range attractive well next to a hard repulsive core, followed by a smooth, long-range repulsive tail. Theoretical calculations are performed in the framework provided by the Ornstein-Zernike equation, solved either analytically with the mean spherical approximation (MSA) or iteratively with the hypernetted-chain (HNC) closure. Our analysis shows that both theories are generally accurate in a thermodynamic region corresponding to a dense vapor phase around the critical point. For a suitable choice of potential parameters, namely, when the attractive well is deep and/or large enough, the static structure factor displays a secondary low-Q peak. In this case HNC predictions closely follow the simulation results, whereas MSA results progressively worsen the more pronounced this low-Q peak is. We discuss the appearance of such a peak, also experimentally observed in colloidal suspensions and protein solutions, in terms of the formation of equilibrium clusters in the homogeneous fluid.     

Specific ion effects in solutions of globular proteins: comparison between analytical models and simulation

Monte Carlo simulations have been performed for ion distributions outside a single globular macroion and for a pair of macroions, in different salt solutions. The model that we use includes both electrostatic and van der Waals interactions between ions and between ions and macroions. Simulation results are compared with the predictions of the Ornstein-Zernike equation with the hypernetted chain closure approximation and the nonlinear Poisson-Boltzmann equation, both augmented by pertinent van der Waals terms. Ion distributions from analytical approximations are generally very close to the simulation results. This demonstrates that properties that are related to ion distributions in the double layer outside a single interface can to a good approximation be obtained from the Poisson-Boltzmann equation. We also present simulation and integral equation results for the mean force between two globular macroions (with properties corresponding to those of hen-egg-white lysozyme protein at pH 4.3) in different salt solutions. The mean force and potential of mean force between the macroions become more attractive upon increasing the polarizability of the counterions (anions), in qualitative agreement with experiments. We finally show that the deduced second virial coefficients agree quite well with experimental results.     

A State-dependent Division Scheme for the Potential of Simple Fluids

Abstract

A new division scheme for the pair potential into long-range and short-range parts is presented, which takes account of the dependence on the density and makes possible the extraction of the bridge function B(r) in the core region. The calculation of the correlation functions is carried out with the hybridized mean spherical approximation (HMSA). As attested by the comparison with the simulation results, HMSA used with this potential separation is suitable to produce accurate bridge function and pair correlation function for the Lennard-Jones fluid.     

Highly asymmetric electrolytes in the primitive model: Hypernetted chain solution in arbitrary spatial dimensions

The pair‐correlation functions for fluid ionic mixtures in arbitrary spatial dimensions are computed in hypernetted chain (HNC) approximation. In the primitive model (PM), all ions are approximated as nonoverlapping hyperspheres with Coulomb interactions. Our spectral HNC solver is based on a Fourier‐Bessel transform introduced by Talman (J. Comput. Phys. 1978, 29, 35), with logarithmically spaced computational grids. Numeric efficiency for arbitrary spatial dimensions is a commonly exploited virtue of this transform method. Here, we highlight another advantage of logarithmic grids, consisting in efficient sampling of pair‐correlation functions for highly asymmetric ionic mixtures. For three‐dimensional fluids, ion size and charge‐ratios larger than 1000 can be treated, corresponding to hitherto computationally not accessed micrometer‐sized colloidal spheres in 1‐1 electrolyte. Effective colloidal charge numbers are extracted from our PM results. For moderately large ion size and charge‐asymmetries, we present molecular dynamics simulation results that agree well with the approximate HNC pair correlations. © 2013 Wiley Periodicals, Inc.     

Electrokinetic properties of monovalent electrolytes confined in charged nanopores: Effect of geometry and ionic short-range correlations

The electrokinetic properties (such as capillary conductance, electroviscosity, and the streaming potential) are obtained for a restricted primitive model electrolyte confined in a slitlike nanopore made up of two infinite parallel plates and in a cylindrical cavity of infinite extension. The hypernetted chain/mean spherical approximation (HNC/MSA) is used to obtain the equilibrium ionic concentration profiles inside the pores, which in turn are used to calculate the electrokinetic properties via linear hydrodynamic equations. Our results are compared with those obtained via the classical Poisson–Boltzmann (PB) theory. Important quantitative and qualitative effects, attributed to geometry and to the proper consideration of short-range correlations by HNC/MSA, are discussed.     

Statistical mechanical theories of the electric double layer

A comparison is made of present-day statistical mechanical theories of the diffuse part of the electric double layer in aqueous 1-1 electrolyte at a charged plane interface. These theories fall into three categories: (1) the modified Poisson-Boltzmann equation (MPB) based on the Kirkwood-Loeb charging of an ion; (2) the adaption to the electric double layer of the Bogoliubov-Born-Green-Yvon (BBGY) hierarchy of integral equations; (3) the use of the Ornstein-Zernike equation (OZ) for the direct correlation functions of the pair interfacial plane wall-molecular particle, as derived by Henderson, Abraham and Barker (HAB). The HAB-OZ equation is used in conjunction with the mean spherical approximation (MSA) or hypernetted chain approximation (HNC). All the theories make use of the primitive model of the bulk electrolyte, so that inhomogeneity of the dielectric permittivity next to the plane wall is neglected. Except perhaps for a variation of the BBGY theory, which uses a closure based on electrical neutrality, all the theories predict oscillatory behaviour in potential distribution as a function of distance at the higher electrolyte concentrations. The HAB-OZ model has the defect that electrostatic imaging is not consistent with the assumptions of central forces and pair-wise additivity of ionic interactions. It is found that the MPB provides the best overall agreement with Monte Carlo calculations.     

Electrolyte solutions that unmix. Hydrophobic ions in water

Simple models for the solvent-averaged ion-ion pair potentials for aqueous solutions of tetraalkylammonium halides were previously treated under the HNC (hypernetted chain) approximation to find the parameters needed to fit osmotic coefficient data of the corresponding real solutions. Here the model Pr ₄ NI is modified by changing the A_+– Gurney parameter to give the Pr ₄ NX model which exhibits the unmixing (separation into two solution phases of different concentrations) that has been reported for several real aqueous tetraalkylammonim salt solutions. The models used here are established at the McMillan-Mayer (MM) or solvent-averaged level, so careful attention is given to the choice of the thermodynamic potential from which we may derive the condition of material stability (stability with respect to separation into two phases of slightly different composition), and calculate the required thermodynamic coefficients from the MM pair correlation functions. The emphasis is on the study of the hydrophobic unmixing in terms of thermodynamic coefficients derived from the pair correlation functions calculated for the Pr ₄ NX model under HNC. In the temperature-concentration (T-c _S ) plane we can locate the solvability line (which separates states for which the HNC equation can be solved from the rest) and portions of the coexistence line. To locate the coexistence line in regimes in which the double tangent method is not effective we use a method based on the osmotic pressure and the solute chemical potential isotherms. Our results suggest that for the Pr ₄ NX model in the T-c _S plane the spinodal line lies within the solvability line which, in turn, lies within the coexistence line for the states where the latter could be determined. The role of the thermodynamic inconsistency implicit in the HNC correlation functions is given special attention, as is the remarkable role of the 1/r ⁴ cavity term in the model pair potential. Single-ion activity coefficients and DEL functions are calculated for some of the states studied.     

Analytical coarse-grained description for polymer melts

Starting from the Ornstein-Zernike equation the authors derive an analytical theory, at the level of pair correlation functions, which coarse grains polymer melts into liquids of interacting soft colloidal particles. Since it is analytical, the presented coarse-graining approach will be useful in developing multiscale modeling procedures to simulate complex fluids of macromolecules. The accuracy of the theory is tested by its capacity to reproduce the liquid structure, as given by the center-of-mass intermolecular total pair correlation function. The theory is found to agree well with the structure predicted by molecular dynamics simulations of the liquid described at the united atom level as well as by molecular dynamics simulations of the liquid of interacting colloidal particles. The authors perform simulations of the liquid of interacting colloidal particles having as input the potential obtained from their analytical total pair correlation function by enforcing the hypernetted-chain closure approximation. Tests systems are polyethylene melts of chains with increasing degrees of polymerization and polymer melts of chains with different chemical architectures. They also discuss the effect of adopting different conventional approximations for intra- and intermolecular monomer structure factors on the accuracy of the coarse-graining procedure, as well as the relevance of higher-order corrections to their expression.     

Exact Self-consistent Integral Equations for the Distribution Functions of Classical Fluids

Abstract

The source particle method (SPM) due to Percus and Verlet of obtaining the single particle density, and the pair and triplet distribution functions of classical fluids (as well as in the variational theory of Bose liquids) is studied. Generalizations of hypernetted chain (HNC) equations are generated by holding fixed the coordinates of the definite group of source particles. Special attention is paid to the triplet distribution function and to the self-consistent calculation of the bridge function (elementary diagram) contribution in the pair distribution function. A comparison with the other exact integral equation theories including the BBGKY-equations is discussed.     

The electrical double layer for a fully asymmetric electrolyte around a spherical colloid: an integral equation study

The hypernetted chain/mean spherical approximation (HNC/MSA) integral equation for a totally asymmetric primitive model electrolyte around a spherical macroparticle is obtained and solved numerically in the case of size-asymmetric systems. The ensuing radial distribution functions show a very good agreement when compared to our Monte Carlo and molecular-dynamics simulations for spherical geometry and with respect to previous anisotropic reference HNC calculations in the planar limit. We report an analysis of the potential versus charge relationship, radial distribution functions, mean electrostatic potential, and cumulative reduced charge for representative examples of 1:1 and 2:2 salts with a size-asymmetry ratio of 2. Our results are collated with those of the modified Gouy-Chapman (MGC) and unequal radius modified Gouy-Chapman (URMGC) theories and with those of HNC/MSA in the restricted primitive model (RPM) to assess the importance of size-asymmetry effects. One of the most striking characteristics found is that, contrary to the general belief, away from the point of zero charge the properties of an asymmetric electrical double layer (EDL) are not those corresponding to a symmetric electrolyte with the size and charge of the counterion, i.e., counterions do not always dominate. This behavior suggests the existence of a new phenomenology in the EDL that genuinely belongs to a more realistic size-asymmetric model where steric correlations are taken into account consistently. Such novel features cannot be described by traditional mean-field theories such as MGC, URMGC, or even by enhanced formalisms, such as HNC/MSA, if they are based on the RPM.     

Maximum entropy formulation of the Kirkwood superposition approximation

Using a variational formulation, we derive the Kirkwood superposition approximation for systems at equilibrium in the thermodynamic limit. We define the entropy of the triplet correlation function and show that the Kirkwood closure brings the entropy to its maximal value. This approach leads to a different interpretation for the Kirkwood closure relation, usually explained by probabilistic considerations of dependence and independence of particles. The Kirkwood closure is generalized to finite volume systems at equilibrium by computing the pair correlation function in finite domains. Closure relations for high order correlation functions are also found using a variational approach. In particular, maximizing the entropy of quadruplets leads to the high order closure g(1234)=g(123)g(124)g(134)g(234)/[g(12)g(13)g(14)g(23)g(24)g(34)] used in the Born-Green-Yvon 2 equations which are a pair of integral equations for the triplet and pair correlation functions.     

Solvent phase behavior and the interaction of uniform and patterned solutes

Isotropic and anisotropic hypernetted-chain (HNC) integral equation theories are used to obtain the interaction of solutes both near and far from the solvent liquid-vapor coexistence. Spherically symmetrical and chemically patterned (patched) solutes are considered, and the influences of particle and patch sizes are investigated. Solvophilic and solvophobic solutes (or patches) are examined. Near coexistence, in the solvophobic case drying-like behavior occurs for solutes (patches) of sufficient size. This gives rise to relatively long ranged attractive forces that are strongly orientation dependent for the patched solute particles. We also report grand canonical Monte Carlo results for a pair of spherically symmetric solutes. This demonstrates that the anisotropic HNC theory gives qualitatively correct solvent structure in the vicinity of the solutes. Comparison with previous simulations also shows that the solute-solute potentials of mean force given by the anisotropic theory are more accurate (particularly at small separations) than those obtained using the isotropic method.     

Ab initio theoretical study of temperature and density dependence of molecular and thermodynamic properties of water in the entire fluid region: autoionization processes

The temperature and density dependence of the molecular and thermodynamic properties of water is investigated theoretically by means of the ab initio electronic structure theory combined with the reference interaction site model method, so-called RISM-SCF. We consider the autoionization process (H2O + H2O right harpoon over left harpoon H3O+ + OH-) by regarding H2O, H3O+, and OH- as "solute" molecules in an aqueous solution and evaluate molecular geometry, electronic structure, solvation structure, and the ionic product of water (pKw) of these species as functions of thermodynamic conditions. In our previous paper, we calculated these properties by using essentially the same method in a wide range of density values (0.6-1.4 g/cm3). However, the calculation was limited at rather higher density (>0.6 g/cm3) due to the difficulty of convergence, which is inherent to the hypernetted-chain (HNC) closure. The problem is overcome in this study by employing the Kovalenko-Hirata (KH) closure which hybridizes the HNC and the mean-spherical approximation (MSA). Here, we present the results for the thermodynamic range of densities from 0.025 to 1.0 g/cm3 and for temperatures from 300 to 800 K including the supercritical point.     

Autocorrelation functions of a fluid of molecules interacting through steeply repulsive and attractive forces

A new method is presented for an extension of Enskog's approximation for the evaluation of the autocorrelation functions of a fluid, and this approach is used to evaluate these functions when the interaction between the molecules includes both steeply repulsive and steeply attractive forces. Consequently the correlation functions depend upon the temperature in a nontrivial way. As an example, the method is applied to calculate the velocity and force autocorrelation functions of a fluid when the molecules interact through the specific potential, V(r)=4epsilon[(sigma/r)2n-(sigma/r)n] when the parameter n is large. There is a relationship between this model and the "sticky sphere" one which is exploited in the theoretical computations. The results obtained from the theory are compared with molecular dynamics simulation for n=72 and 144 and for a range of temperatures from T=epsilon/kB down to epsilon/3kB. The two approaches agree very well for a range of state points, especially at short times. At later times the theory predicts a more oscillatory behavior than the simulation especially at very low reduced temperatures.     

A local density functional approximation for the ion distribution near a charged electrode in the restricted primitive model electrolyte

A local HNC/MSA approximation is developed and applied to the 1:1 restricted primitive model electrolyte. Improvement of ion density profiles in front of a charged electrode is achieved by employing ion—ion direct correlation functions from homogeneous systems of non-neutral composition as found locally in the inhomogeneous double layer. This approximation is related to the density functional approach for inhomogeneous fluids. The local HNC/MSA method predicts, at higher surface charges, layering of counterions and charge inversion as seen in Monte Carlo data, a strong increase of the surface potential with charging and a maximum in the double layer capacitance.     

Pair diffusion, hydrodynamic interactions, and available volume in dense fluids

We calculate the pair diffusion coefficient D(r) as a function of the distance r between two hard sphere particles in a dense monodisperse fluid. The distance-dependent pair diffusion coefficient describes the hydrodynamic interactions between particles in a fluid that are central to theories of polymer and colloid dynamics. We determine D(r) from the propagators (Green's functions) of particle pairs obtained from molecular dynamics simulations. At distances exceeding ~3 molecular diameters, the calculated pair diffusion coefficients are in excellent agreement with predictions from exact macroscopic hydrodynamic theory for large Brownian particles suspended in a solvent bath, as well as the Oseen approximation. However, the asymptotic 1/r distance dependence of D(r) associated with hydrodynamic effects emerges only after the pair distance dynamics has been followed for relatively long times, indicating non-negligible memory effects in the pair diffusion at short times. Deviations of the calculated D(r) from the hydrodynamic models at short distances r reflect the underlying many-body fluid structure, and are found to be correlated to differences in the local available volume. The procedure used here to determine the pair diffusion coefficients can also be used for single-particle diffusion in confinement with spherical symmetry.     

Small angle neutron scattering study of sodium dodecyl sulfate micellar growth driven by addition of a hydrotropic salt

The structures of aggregates formed in aqueous solutions of an anionic surfactant, sodium dodecyl sulfate (SDS), with the addition of a cationic hydrotropic salt, p-toluidine hydrochloride (PTHC), have been investigated by small angle neutron scattering (SANS). The SANS spectra exhibit a pronounced peak at low salt concentration, indicating the presence of repulsive intermicellar interactions. Model-independent real space information about the structure is obtained from a generalized indirect Fourier transformation (GIFT) technique in combination with a suitable model for the interparticle structure factor. The interparticle interaction is captured using the rescaled mean spherical approximation (RMSA) closure relation and a Yukawa form of the interaction potential. Further quantification of the geometrical parameters of the micelles was achieved by a complete fit of the SANS data using a prolate ellipsoidal form factor and the RMSA structure factor. The present study shows that PTHC induces a decrease in the fractional charge of the micelles due to adsorption at the micellar surface and consequent growth of the SDS micelles from nearly globular to rodlike as the concentration of PTHC increases.