Zeta potential of colloidal particle in solvent primitive model electrolyte solution: a density functional theory study

A systematic study of zeta potential for a spherical double layer (SDL) around a colloidal particle in electrolyte solutions, is performed using density functional theory and Monte Carlo simulation. The usual recipe under the solvent primitive model is employed to model the system, where macroion, counterions, and coions are represented by charged hard spheres of uniform charge density and the presence of solvent is taken into account by modelling it as neutral hard spheres. All the components of the system are embedded in a dielectric continuum in order to consider the electrostatic effect of the solvent. The density functional theory employs a suitable weighted density approximation to calculate the hard-sphere contribution, whereas the residual electrostatic interactions are calculated as a small perturbation around the uniform fluid. The zeta potential profiles of a SDL in the presence of a number of electrolytes have been calculated and are found to be considerably influenced in the presence of solvent with an increase in the concentration of the electrolyte. The theory successfully predicts the maxima and sign reversal of the zeta potential profiles at high macroion surface charge density and in the presence of multivalent counterions, as obtained from the Monte Carlo simulation.     

Effect of ionic size on the structure of spherical double layers: a Monte Carlo simulation and density functional theory study

The effect of ionic size on the diffuse layer characteristics of a spherical double layer is studied using Monte Carlo simulation and density functional theory within the restricted primitive model. The macroion is modelled as an impenetrable charged hard sphere carrying a uniform surface charge density, surrounded by the small ions represented as charged hard spheres and the solvent is taken as a dielectric continuum. The density functional theory uses a partially perturbative scheme, where the hard sphere contribution to the one particle correlation function is evaluated using weighted density approximation and the ionic interactions are calculated using a second-order functional Taylor expansion with respect to a bulk electrolyte. The Monte Carlo simulations have been performed in the canonical ensemble. The detailed comparison is made in terms of zeta potentials for a wide range of physical conditions including different ionic diameters. The zeta potentials show a maximum or a minimum with respect to the polyion surface charge density for a divalent counterion. The ionic distribution profiles show considerable variations with the concentration of the electrolyte, the valency of the ions constituting the electrolyte, and the ionic size. This model study shows clear manipulations of ionic size and charge correlations in dictating the overall structure of the diffuse layer.     

Depletion interactions between colloidal particles in polymer solutions: density functional approach

A density functional theory for colloid–polymer mixtures based on the weighted-density approximation has been developed to investigate the depletion effects acting between two colloids immersed in a bath of polymers and the depletion effects for a colloid near a planar hard wall. The theoretical results for the polymer-induced depletion interactions and the local polymer density distributions are in good agreement with the computer simulations. The calculation shows that the depletion interaction for a colloid near a planar hard wall is much stronger than that between two colloids in a polymer solution because of the strong confinement effect. The behaviour of the depletion interactions has been analysed as a function of the polymer density, the polymer chain length, and the colloid/polymer size ratio. Strong depletion effects appear in short-chain systems and with large colloid/polymer size ratios.     

Structural and thermodynamic properties of the restricted primitive model electrolyte in a mixture with uncharged hard spheres: a grand canonical Monte Carlo simulation and integral equation study

A restricted primitive model electrolyte in a mixture with uncharged hard spheres was studied at room temperature using grand canonical Monte Carlo computer simulation and Ornstein–Zernike integral equation theory in the hypernetted chain approximation (HNC). The mean spherical approximation results are also presented for a few cases. We obtained the pair distribution functions of species of the system, the dependencies of the total fluid density and the ionic fraction on the chemical potentials, the excess internal energy and the heat capacity at constant volume for a wide range of chemical potentials of the species from the simulations and HNC theory. In the majority of cases, good agreement between the theoretical predictions and simulation data is obtained. The composition of the mixture is determined by the chemical potentials of both species. The pair distribution functions have a Debye-like shape at low densities for various values of the ion fraction. By increasing the chemical potential of the uncharged component, weak trends for structuring of the solution are observed with the formation of ion-hard sphere-ion complexes. At high densities, a tendency for in-phase oscillations of ion–ion functions is observed similar to the pure electrolyte in the restricted primitive model. We analysed the chemical potential–density and the chemical potential–ion fraction projections of the equation of state in detail. Also, the heat capacity at constant volume has been calculated for the first time. The model and the results are useful for the development of the theory of inhomogeneous fluid mixtures.     

An association theory for the formation of ion oligomers and its application to 1-1 electrolytes

A theory for modelling electrolyte solutions which includes the formation of ion clusters of different size has been developed in the framework of the primitive model. In primitive models the solvent is described as a dielectric continuum and the solvent–solute interactions are neglected. For the dielectric constant the value of the pure solvent has been used. The ion cluster distribution is calculated from the mass action law. The association constants are related to integrals over the cluster distribution functions which are calculated with the Kirkwood superposition approximation from low density pair distribution functions. The ion clusters are defined by a certain distance which rules if two ions belong to the same cluster. This so-called Bjerrum distance is chosen according to fundamental investigations of the structure of ion cluster. All ion clusters are modelled as hard spheres. For the free ions and charged clusters the mean spherical approximation expression for the Coulomb interaction is added. The co-volumes of the ion clusters have been taken from the investigation of the ion cluster structures, and are consistent with the definition of an ion cluster chosen here.     

A self-consistent density-functional approach to the structure of electric double layer: charge-asymmetric electrolytes

A self-consistent density-functional approach has been employed to study the structure of an electric double layer formed from a charge-asymmetric (2:l) electrolyte within the restricted primitive model which corresponds to charged hard sphere ions and a continuum solvent. The particle correlation due to hard-core exclusions is evaluated by making use of the universality of the density functionals and the correlation function of the uniform hard sphere fluid obtained through the integral equation theory with an accurate closure relation whereas mean spherical approximation is employed for the electrical contribution. Numerical results on the diffuse layer potential drop, ionic density profile, and the mean electrostatic potential near the electrode surface at several surface charge densities are found to be in quantitative agreement with the available simulation data.     

Contact value relations and density functional theory for the electrical double layer

The recently developed contact-corrected density functional theory is applied to an electric double layer containing a valence and size symmetric electrolyte. The restricted primitive model is used to describe the electrolyte, and a non-polarisable, smooth, planar, uniformly charged hard surface is used to model the electrode. Comparison with Monte Carlo simulation results for this system shows that the theory satisfies the first contact value relation in line with the premise of the approach. Additionally, the theory is found to satisfy the semi-empirical second contact value relation quite closely. This is a welcome result, since no reference is made to this second relation nor is the relation used in the theoretical development. The second contact relation is closely related to the anomalous positive temperature derivative of the double layer capacitance at low temperature, low ion concentration, and low electrode charge, and indeed, the contact-corrected density functional theory predicts this effect.     

A new study of associating inhomogeneous fluids with classical density functional theory

ABSTRACT

A new density functional for the study of associating inhomogeneous fluids based on Wertheim's first-order thermodynamic perturbation theory is presented and compared to the most currently used associating density functionals. This functional is developed using the weighted density approximation in the range of association of hard spheres. We implement this functional within the framework of classical density functional theory together with modified fundamental measure theory to account for volume exclusion of hard spheres. This approach is tested against molecular simulations from literature of pure associating hard spheres and mixtures of non-associationg and associating hard spheres with different number of bonding sites close to a hard uniform wall. Furthermore, we compare and review our results with the performance of associating functionals from literature, one based on fundamental measure theory and the inhomogeneous version of Wertheim's perturbation theory. Results obtained with classical DFT and the three functionals show excellent agreement with molecular simulations in systems with one hard wall. For the cases of small pores where only one or two layers of fluid are allowed discrepancies between results with classical DFT and molecular simulations were found.     

Solvation force in hard ellipsoid fluids with HNW interaction

In present work, using density functional theory and extended restricted orientation model, the one particle density of hard Gaussian overlap fluid near the colloid walls is calculated. The hard needle–wall interaction between molecules and colloids are considered. Using non-linear equation, proposed by Grimson–Rickyazen, the solvation force of hard ellipsoidal molecular fluid with hard Gaussian overlap interaction is calculated. We could not find the exact or simulation results for comparison. The results in the case k = 2.0 are compared with the solvation force of one-dimensional hard rod fluids. The results are corresponded, qualitatively.     

Second-order Percus-Yevick theory for a confined hard-sphere fluid

A fluid of hard spheres confined between two hard walls and in equilibrium with a bulk hard-sphere fluid is studied using a second-order Percus-Yevick approximation. We refer to this approximation as second-order because the correlations that are calculated depend upon the position of two hard spheres in the confined fluid. However, because the correlation functions depend upon the positions of four particles (two hard spheres and two walls treated as giant hard spheres), this is the most demanding application of the second-order theory that has been attempted. When the two walls are far apart, this calculation reduces to our earlier second-order approximation calculations of the properties of hard spheres near a single hard wall. Our earlier calculations showed this approach to be accurate for the single-wall case. In this work we calculate the density profiles and the pressure of the hard-sphere fluid on the walls. We find, by comparison with grand canonical Monte Carlo results, that the second-order approximation is very accurate, even when the two walls have a small separation. We compare with a singlet approximation (in the sense that correlation functions that depend on the position of only one hard sphere are considered). The singlet approach is fairly satisfactory when the two walls are far apart but becomes unsatisfactory when the two walls have a small separation. We also examine a simple theory of the pressure of the confined hard spheres, based on the usual Percus-Yevick theory of hard-sphere mixtures. Given the simplicity of the latter approach the results of this simple (and explicit) theory are surprisingly good.     

On the relationship between the density functional formalism and the potential distribution theory for nonuniform fluids

It is shown that the variational principle for the grand potential of a nonuniform fluid as a functional of the singlet density yields the potential distribution theory for the equilibrium density. We derive the explicit form that the functional takes for a system of hard rods, and propose an approximate one for hard spheres. Attractive interactions are also considered in mean-field approximation. In all cases the pair direct correlation function of the nonuniform system is obtained and the density gradient expansion of the free energy is investigated.     

Adsorption isotherm of a Lennard-Jones nitrogen in a carbon slitlike pore

A density functional perturbation approximation based both on second-order perturbation theory and on the pore average density has been proposed to study the adsorption hysteresis of nitrogen in a carbon slit pore. The main advantage of the present approximation is that it is computationally much simpler than the original density functional approximation based on the second-order perturbation theory of liquids, and can be applied to several model fluids confined in a strong external field in order to study their structural and thermodynamic properties. The calculated adsorption hysteresis for the confined Lennard-Jones nitrogen is in very good agreement with computer simulation, even if its accuracy slightly deteriorates for the desorption branch. The calculated equilibrium particle density distributions also compare well with computer simulations, and are better than those of a density functional theory based on the so-called mean-field approximation.     

Specification of Density Functional Approximation by Radial Distribution Function of Bulk Fluid

A systematic methodology is proposed to deal with the weighted density approximation version of classical density functional theory by employing the knowledge of radial distribution function of bulk fluid.The present methodology results from the concept of universality of the free energy density functional combined with the test particle method.It is shown that the new method is very accurate for the predictions of density distribution of a hard sphere fluid at different confining geometries.The physical foundation of the present methodology is also applied to the quantum density functional theory.     

Forces between charged surfaces in a solvent primitive model electrolyte

Grand canonical Monte Carlo (GCMC) simulations have been performed to investigate the components of the force between parallel charged surfaces in an electrolyte. The solvent primitive model (SPM) was used to investigate the effect of neutral hard sphere solvent particles on the force between the surfaces. The effects of particle size, wall charge density, charge valency of the electrolyte, and the exclusion of neutral hard sphere are discussed. When solvent particles are considered, the total force between the charged surfaces is always repulsive, even for divalent counterions. This is different from the earlier conclusion reached with a restricted primitive model electrolyte. The repulsive force decreases in going from monovalent counterions to divalent counterions.     

Development of electron-proton density functionals for multicomponent density functional theory

We present a strategy for the development of electron-proton density functionals in multicomponent density functional theory, treating electrons and selected nuclei quantum mechanically without the Born-Oppenheimer approximation. An electron-proton functional is derived using an explicitly correlated electron-proton pair density. This functional provides accurate hydrogen nuclear densities, thereby enabling reliable calculations of molecular properties. This approach is potentially applicable to relatively large molecular systems with key hydrogen nuclei treated quantum mechanically.     

On the signs of fluorine spin coupling constants

From an exact expression for the free energy of a non-uniform fluid mixture a closure approximation for the inhomogeneous direct correlation functions is used to develop a theory of solvation forces in charged fluids based upon non-linear equations for the equilibrium ion number densities. In the limit of point ions, the expressions obtained reduce to those of the Poisson-Boltzmann theory of electrolytes. The numerical results obtained for a restricted primitive model electrolyte are compared with those of earlier work based on linear response theory and Poisson-Boltzmann theory with a simple Stern layer modification. At low electrolyte concentrations the agreement between all three theories is good. But at high electrolyte concentrations the Poisson-Boltzmann theory with a simple Stern layer correction fails to display the oscillations in the solvation force which characterize both the linear and non-linear theories.     

Associating fluids with four bonding sites against a hard wall: density functional theory

We present two new perturbation density functional theories to investigate non-uniform fluids of associating molecules. Each fluid molecule is modelled as a spherical hard core with four highly anisotropic square well sites placed in tetrahedral symmetry on the hard core surface. In one theory we apply the weighting from Tarazona's hard sphere density functional theory to Wertheim's bulk first-order perturbation theory. The other theory uses the inhomogeneous form of Wertheim's theory as a perturbation to Tarazona's hard-sphere density functional theory. Each theory approaches Tarazona's theory in the limit of zero association. We compare results from theory and simulation for density profiles, fraction of monomers, and adsorption of an associating fluid against a hard, smooth wall over a range of temperatures and densities. The non-uniform fluid theory which uses Tarazona's weighting of Wertheim's theory in the bulk is in good agreement with computer simulation results.     

Influence of surface affinity of planar walls on effective interaction between the wall and colloids

Using density functional theory, we investigate the effective interaction between a big colloid immersed in a sea of small colloids and a wall which has different affinity to the small colloids. Steele 10-4-3 potential is introduced to mimic both short-range repulsive and long-range attractive interactions between the wall and the small colloids. It is found that the surface affinity of the wall has a significant influence on the effective interaction. In the short-range repulsive case, the repulsion greatly enhances the big colloid-wall effective attraction, which sensitively depends on the concentration of small colloids, and is not sensitive to the repulsive strength. In the long-range attractive case, both the concentration of small colloids and the attractive strength have great effect on the effective interaction, and with an increase of the attractive strength, a strong repulsion may be induced when the big colloid is close to the wall. In low density limit of small colloids, the present results agree well with those of the Asakura and Oosawa(AO) approximation.