Surface properties of fluids of charged platelike colloids

Surface properties of mixtures of charged platelike colloids and salt in contact with a charged planar wall are studied within density functional theory. The particles are modeled by hard cuboids with their edges constrained to be parallel to the Cartesian axes corresponding to the Zwanzig model [J. Chem. Phys. 39, 1714 (1963)] and the charges of the particles are concentrated at their centers. The density functional applied is an extension of a recently introduced functional for charged platelike colloids. It provides a qualitative approach because it does not determine the relation between the actual and the effective charges entering into the model. Technically motivated approximations, such as using the Zwanzig model, are expected not to influence the results qualitatively. Analytically and numerically calculated bulk and surface phase diagrams exhibit first-order wetting for sufficiently small macroion charges and isotropic bulk order as well as first-order drying for sufficiently large macroion charges and nematic bulk order. The asymptotic wetting and drying behaviors are investigated by means of effective interface potentials which turn out to be asymptotically the same as for a suitable neutral system governed by isotropic nonretarded dispersion forces. Wetting and drying points as well as predrying lines and the corresponding critical points have been located numerically. A crossover from monotonic to nonmonotonic electrostatic potential profiles upon varying the surface charge density has been observed. Nonmonotonic electrostatic potential profiles are equivalent to the occurrence of charge inversion. Due to the presence of both the Coulomb interactions and the hard-core repulsions, the surface potential and the surface charge do not vanish simultaneously, i.e., the point of zero charge and the isoelectric point of the surface do not coincide.     

Free isotropic-nematic interfaces in fluids of charged platelike colloids

Bulk properties and free interfaces of mixtures of charged platelike colloids and salt are studied within the density-functional theory. The particles are modeled by hard cuboids with their edges constrained to be parallel to the Cartesian axes corresponding to the Zwanzig model. The charges of the particles are concentrated in their center. The density functional is derived by functional integration of an extension of the Debye-Hückel pair distribution function with respect to the interaction potential. For sufficiently small macroion charges, the bulk phase diagrams exhibit one isotropic and one nematic phase separated by a first-order phase transition. With increasing platelet charge, the isotropic and nematic binodals are shifted to higher densities. The Donnan potential between the coexisting isotropic and nematic phases is inferred from bulk structure calculations. Nonmonotonic density and nematic order parameter profiles are found at a free interface interpolating between the coexisting isotropic and nematic bulk phases. Moreover, electrically charged layers form at the free interface leading to monotonically varying electrostatic potential profiles. Both the widths of the free interfaces and the bulk correlation lengths are approximately given by the Debye length. For fixed salt density, the interfacial tension decreases upon increasing the macroion charge.     

Effect of Discrete Macroion Charge Distributions on Electric Double Layer of a Spherical Macroion

The effect of replacing the conventional uniform macroion surface charge density with discrete macroion charge distributions on the structure of electric double layer (EDL) of a spherical macroion has been investigated by Monte Carlo (MC) simulations. Two discrete models have been investigated in addition to the central macroion charge: point charges localized on the macroion surface and finite-sized charges protruding into the solution. Both models have been studied with fixed and mobile macroion charges. The radial functions of local densities and electrostatic potential in EDL, are calculated and compared to the results obtained for the central macroion charge distribution. It is concluded that the model of charge distribution significantly affects the EDL structure close to the macroion, while the effect is much weaker at larger distances. With point charges localized on the macroion surface, counterions become stronger accumulated to the macroion, as a result the absolute values of surface potential ?₀ and zeta ξ potential are decreased. With protruding charges, the excluded volume effect dominates over the increased correlation ability; hence the counterions are less accumulated near the macroions and the absolute values of ?₀ and ξ potentials are increased.     

Effect of discrete macroion charge distributions in solutions of like-charged macroions

The effect of replacing the conventional uniform macroion surface charge density with discrete macroion charge distributions on structural properties of aqueous solutions of like-charged macroions has been investigated by Monte Carlo simulations. Two discrete charge distributions have been considered: point charges localized on the macroion surface and finite-sized charges protruding into the solution. Both discrete charge distributions have been examined with fixed and mobile macroion charges. Different boundary conditions have been applied to examine various properties. With point charges localized on the macroion surface, counterions become stronger accumulated to the macroion and the effect increases with counterion valence. As a consequence, with mono- and divalent counterions the potential of mean force between two macroions becomes less repulsive and with trivalent counterions more attractive. With protruding charges, the excluded volume effect dominates over the increased correlation ability; hence the counterions are less accumulated near the macroions and the potential of mean force between two macroions becomes more repulsive/less attractive.     

Effects of structural properties of the Stern layer on the electrophoretic migration of a highly charged spherical macroion

The electrophoretic migration of a highly charged spherical macroion suspended in an aqueous solution of NaCl is studied using the molecular dynamic method. The objective is to examine the effects of the colloidal surface charge density on the electrophoretic mobility (μ) of the spherical macroion. The bare charge and the size of the macroion are varied separately to induce changes in the colloidal surface charge density. Our results indicate that μ depends on colloidal surface charge density in a nonmonotonic manner, but that this relationship is independent of the way the surface charge density is varied. It is found that an increase in colloidal surface charge density may lead to the formation of new sublayers in the Stern layer. The μ profile is also found to have a local maximum for a bare charge at which a new sublayer is formed in the Stern layer, and a local minimum for a bare charge at which the outer sublayer becomes relatively dense. Finally, the electrophoretic flow caused by the migration of the spherical macroion is studied to find that one decisive factor causing the electrophoretic flow is the ability of the macroion to carry anions in the electrolyte solution.     

Long-range many-body polyelectrolyte bridging interactions

We investigate polyelectrolyte bridging interactions mediated by charged, flexible, polyelectrolyte chains between fixed cylindrical macroions of opposite charge in a two-dimensional hexagonal crystalline array. We show that in the asymptotic regime of small macroion density, the polyelectrolyte-mediated attraction is long range, falling off approximately linearly with the macroion array density. We investigate the polyelectrolyte free energy as a function of the macroion density and derive several analytic limiting laws valid in different regimes of the parameter space.     

Structure of inhomogeneous polymer solutions: a density functional approach

The structure of polymer solutions confined between surfaces is studied using a density functional theory where the polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres. The present theory uses the concept of universality of the free energy density functional to obtain the first-order direct correlation function of the nonuniform system from that of the corresponding uniform system, calculated through the Verlet-modified type bridge function. The uniform bulk fluid direct correlation function required as input has been calculated from the reference interaction site model integral equation theory using the Percus-Yevick closure relation. The calculated results on the density profiles of the polymer as well as the solvent are shown to compare well with computer simulation results.     

Influence of atomistic physics on electro-osmotic flow: an analysis based on density functional theory

Molecular density profiles and charge distributions determined by density functional theory (DFT) are used in conjunction with the continuum Navier-Stokes equations to compute electro-osmotic flows in nanoscale channels. The ion species of the electrolyte are represented as centrally charged hard spheres, and the solvent is treated as a dense fluid of neutral hard spheres having a uniform dielectric constant. The model explicitly accounts for Lennard-Jones interactions among fluid and wall molecules, hard sphere repulsions, and short range electrical interactions, as well as long range Coulombic interactions. Only the last of these interactions is included in classical Poisson-Boltzmann (PB) modeling of the electric field. Although the proposed DFT approach is quite general, the sample calculations presented here are limited to symmetric monovalent electrolytes. For a prescribed surface charge, this DFT model predicts larger counterion concentrations near charged channel walls, relative to classical PB modeling, and hence smaller concentrations in the channel center. This shifting of counterions toward the walls reduces the effective thickness of the Debye layer and reduces electro-osmotic velocities as compared to classical PB modeling. Zeta potentials and fluid speeds computed by the DFT model are as much as two or three times smaller than corresponding PB results. This disparity generally increases with increasing electrolyte concentration, increasing surface charge density and decreasing channel width. The DFT results are found to be comparable to those obtained by molecular dynamics simulation, but require considerably less computing time.     

Molecular Dynamics Simulation of the Electrical Double Layer of Spherical Macroion

The structure of the electrical double layer (EDL) of a spherical macroion with a total charge of 60 elementary charges is studied by molecular dynamics methods. In calculations we used two models: continuous and discrete. In the continuous model, the total charge was concentrated in the center of the macroion; in the discrete model, elementary charges were randomly distributed over the surface of the macroion. The radial profiles of local densities and electric potential in EDL, as well as the degree of counterion binding by the macroion, are calculated with allowance for the Lennard-Jones and electrostatic interactions. It is established that the character of charge distribution significantly affects the EDL structure near the macroion, whereas its effect is much weaker at larger distances. The results obtained are compared with the experimental data on the surface potential and the diffuse part of EDL of sodium dodecyl sulfate micelles in aqueous solution, as well as on the micelle-bound charge. It is shown that even weak specific interaction between counterions and a macroion can substantially influence the structure of its EDL.     

Structures and correlation functions of multicomponent and polydisperse hard-sphere mixtures from a density functional theory

The structures of nonuniform binary hard-sphere mixtures and the correlation functions of uniform ternary hard-sphere mixtures were studied using a modified fundamental-measure theory based on the weight functions of Rosenfeld [Rosenfeld, Phys. Rev. Lett. 63, 980 (1989)] and Boublik-Mansoori-Carnahan-Starling-Leland equation of state [Boublik, J. Chem. Phys. 53, 471 (1970); Mansoori et al., J. Chem. Phys. 54, 1523 (1971)]. The theoretical predictions agreed very well with the molecular simulations for the overall density profiles, the local compositions, and the radial distribution functions of uniform as well as inhomogeneous hard-sphere mixtures. The density functional theory was further extended to represent the structure of a polydisperse hard-sphere fluid near a hard wall. Excellent agreement was also achieved between theory and Monte Carlo simulations. The density functional theory predicted oscillatory size segregations near a hard wall for a polydisperse hard-sphere fluid of a uniform size distribution.     

Density-functional theory of spherical electric double layers and zeta potentials of colloidal particles in restricted-primitive-model electrolyte solutions

A density-functional theory is proposed to describe the density profiles of small ions around an isolated colloidal particle in the framework of the restricted primitive model where the small ions have uniform size and the solvent is represented by a dielectric continuum. The excess Helmholtz energy functional is derived from a modified fundamental measure theory for the hard-sphere repulsion and a quadratic functional Taylor expansion for the electrostatic interactions. The theoretical predictions are in good agreement with the results from Monte Carlo simulations and from previous investigations using integral-equation theory for the ionic density profiles and the zeta potentials of spherical particles at a variety of solution conditions. Like the integral-equation approaches, the density-functional theory is able to capture the oscillatory density profiles of small ions and the charge inversion (overcharging) phenomena for particles with elevated charge density. In particular, our density-functional theory predicts the formation of a second counterion layer near the surface of highly charged spherical particle. Conversely, the nonlinear Poisson-Boltzmann theory and its variations are unable to represent the oscillatory behavior of small ion distributions and charge inversion. Finally, our density-functional theory predicts charge inversion even in a 1:1 electrolyte solution as long as the salt concentration is sufficiently high.     

Viscosity and electrophoretic mobility of cesium fullerenehexamalonate in aqueous solutions--comparing experiments and theories on nanometer-sized spherical polyelectrolyte

The viscosity of aqueous solutions of cesium fullerenehexamalonate T h -C 66(COOCs) 12, a rigid spherical nanometer-sized polyvalent salt, was measured by the Ubbelohde-type viscometer. The measurements were performed without added salt at 25 degrees C in the concentration range between 7 and 320 g/dm (3). THe concentration dependence of the obtained reduced viscosity was compared with the theoretical prediction, taking into account contributions stemming from the intrinsic viscosity, hydrodynamic perturbations of the hypothetically bare fullerenehexamalonate macroion, the primary electroviscous effect, and the secondary electroviscous effect. Using the geometric radius of the bare macroion from the previous measurements of the estimated effective charge of the macroion and from the small-angle X-ray scattering data of the estimated thickness of the compact shell of counterions electrostatically bound to the macroion, a good agreement between theory and the experiment was obtained in the range of the lowest and of the highest concentrations. Electrostatic interactions are identified as the main cause of the increased reduced viscosity at the lowest measured concentrations. At the highest concentrations, electrostatic interactions are effectively screened, and the influence of binary hydrodynamic interactions and perturbations of the hypothetical bare macroion prevails over electrostatic contributions to the increased viscosity. The electrophoretic mobility of the fullerenehexamalonate ion in aqueous salt-free medium was computed with the same value for the radius of the fullerenehexamalonate macroion as that used in the calculation of viscosity. The numerical solution of Ohshima's equation agreed well with the experimental values.     

A macroion electrokinetics algorithm

A numerical algorithm is presented for the standard model of macroion electrokinetics and certain generalizations of it. The macroion consists of a cylindrical section with identical, hemispheroidal endcaps, each piece having arbitrary length. The system of one macroion and adjoining salt solution is subjected to an arbitrary sequence of pulsed electrical fields and pulsed translational and rotational velocities. Numerical solutions are obtained for the time dependent electrostatic and mobile ion concentration fields and the solvent velocity. From these fields the dielectric response, force, and torque are calculated. Generalizations of the standard model include the diffusive motion of macroion surface charges, partial slip of solvent motion at the macroion surface, and a simple model for the reactive exchange of surface charge with solution ions. The primary illustrative application is to recent measurements of electric birefringence versus applied field frequency for poly-(tetrafluorothylene) colloidal particles, but a few results are presented for the dielectric response of DNA fragments and of spherical colloidal particles. The source code and additional details are provided as supplementary documentation.     

On the parallel-perpendicular transition for a nematic phase at a wall

We use an Onsager-level density functional theory to investigate the behaviour of the nematic phase in contact with a solid wall. The nematic consists of hard rigid rods having perfect uniform alignment and uniform spatial density. In the absence of any particle-wall interactions besides excluded-volume forces, we predict a director orientation parallel to the wall. We show that this preference for parallel alignment is due to the entropy associated with the larger volume available to the particles in their parallel orientation. An adsorption energy favouring normal alignment gives rise to a transition from a high temperature parallel orientation to a low temperature normal orientation. We derive expressions for the temperature of this transition, relating it explicitly to the wall adsorption energy, particle axial ratio, and nematic density. Effects such as layering near the wall and imperfect nematic order are argued not to be necessary for the existence of this transition.     

Effect of attractions on the structure of polymer solutions confined between surfaces: a density functional approach

A density functional theory is presented to study the effect of attractions on the structure of polymer solutions confined between surfaces. The polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres, both having Yukawa-type attractions and the mixture being confined between attractive Yukawa-type surfaces. The present theory treats the ideal gas free energy functional exactly and uses weighted density approximation for the hard chain and hard sphere contributions to the excess free energy functional. The attractive interactions are calculated using the direct correlation function obtained from the polymer reference interaction site model theory along with the mean spherical approximation closure. The theoretical predictions on the density profiles of the polymer and the solvent molecules are found to agree quite well with the Monte Carlo simulation results for varying densities, chain lengths, wall separations, and different sets of interaction potentials.     

Renormalized charge in a two-dimensional model of colloidal suspension from hypernetted chain approach

The renormalized charge of a simple two-dimensional model of colloidal suspension was determined by solving the hypernetted chain approximation and Ornstein-Zernike equations. At the infinite dilution limit, the asymptotic behavior of the correlation functions is used to define the effective interactions between the components of the system and these effective interactions were compared to those derived from the Poisson-Boltzmann theory. The results we obtained show that, in contrast to the mean-field theory, the renormalized charge does not saturate, but exhibits a maximum value and then decays monotonically as the bare charge increases. The results also suggest that beyond the counterion layer near to the macroion surface, the ionic cloud is not a diffuse layer which can be handled by means of the linearized theory, as the two-state model claims, but a more complex structure is settled by the correlations between microions.     

Computer simulation of macroion layering in a wedge film

The layering of macroions confined to a wedge slit formed by two uncharged hard walls is studied using a canonical Monte Carlo method combined with a simulation cell that contains both wedge-shaped slit and bath regions. The macroion solution is modeled within a one-component fluid approach that in an effective way incorporates the double layer repulsion due to simple electrolyte ions as well as the discrete nature of an aqueous solvent. The layer formation under a wedge confinement is analyzed by carrying out separate simulation runs for a set of consecutive wedge segments designed to represent a single wedge slit. As the wedge thickness progressively increases, the sequence of regions along the wedge film with distinct features of macroion layering has been established. This sequence comprises (i) a wedge region of the thickness smaller than the macroion diameter that is free of macroions; (ii) a region with a one-dimensional macroion chain along the wedge corner at a wedge thickness of a one macroion diameter; (iii) a region comprising a low-ordered macroion monolayer that extends until the wedge thickness slightly above two macroion diameters; (iv) a region comprising a pair of well-defined two-dimensional configurations of macroions segregated on each of the wedge walls; and (v) a free-of-macroions wedge region between two surface monolayers that now originates from an electrostatic repulsion imposed by the surface macroions, which is followed by (vi) a well-defined macroion monolayer film between two surface monolayers, a less defined bilayer film, a three-layer film, and so on up to the bulk solution. Once formed, the macroion surface monolayers persist for all remaining wedge thicknesses up to the bulk, forming in such a way effective charged wedge boundaries. Such a formation of the macroion surface monolayers on the uncharged confining walls is related to the haloing mechanism for regulating the stability in colloidal suspensions [Tohver et al. Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 8951] and is discussed as well. Finally, the estimated boundary of the free of macroion region between surface monolayers correlates well with the location of the boundary of the so-called "vacuum" phase that has been observed experimentally in an aqueous suspension of charged polystyrene spheres bounded by electrostatically repulsive glass walls [Pieransky et al. Phys. Rev. Lett. 1983, 50, 900].     

Application of the Mermin entropy principle to the “apparatus” density functional theory

An “apparatus” operator approach is presented for the extension of the density functional theory of Hohenberg, Kohn, and Mermin. Using the Mermin entropy principle, a one-to-one correspondence is established between the density matrix for the system and the electron charge density for a finite-temperature system in the presence of an apparatus. In the zero-temperature limit in the absence of an apparatus, the Hohenberg–Kohn theory is recovered. The central aspect of this new density functional theory is that the principle of maximum entropy is applied to the system plus its surroundings under the additional constraint that the electron charge density is given. The system is treated as a subsystem of a composite system and is not necessarily in the equilibrium state as in the Mermin theory. As an example, it is shown how, in principle, excited states are encompassed by the theory.