Analytical expressions for the electrical potential near planar, cylindrical, and spherical surfaces for symmetric electrolytes

The conjecture of Tuinier (J. Colloid Interface Sci. 258 (2003) 45) for the electrical potentials near a cylindrical surface and near a spherical surface under the conditions of symmetric electrolyte and large scaled radius are derived by solving the corresponding Poisson-Boltzmann equation. The surface charge density-surface potential relations for these surfaces are also derived under the conditions of constant surface potential. We show that the level of surface charge density for planar, cylindrical, and spherical surfaces follows the order spherical surface > cylindrical surface > planar surface.     

Approximate analytical expressions for the electrical potential in a cavity containing salt-free medium

The electrical potential in a closed surface such as a cavity containing counterions only is derived for the cases of constant surface potential and constant surface charge density. The results obtained have applications in, for example, microemulsion-related systems in which ionic surfactants are introduced to maintain the stability of a dispersion and electroosmotic flow-related analysis. An analytical expression for the electrical potential is derived for a planar slit, and the methodology used is modified to derive approximate analytical expressions for spherical and cylindrical cavities. The higher the surface potential, the better the performance of these expressions. For the case where the surface potential is above ca. 50 mV, the performance of the approximate analytical expressions can further be improved by multiplying a correction function.     

Derivation of analytical expressions for the electrical potential distribution in lipid structures

The electrical potential inside a lipid structure, which is described by a modified Poisson-Boltzmann equation in the literature (Borukhov et al. Electrochim. Acta 2000, 46, 221), is solved, taking into account the effects of ionic sizes. Here, a micelle comprises an ionic surfactant layer and an aqueous core; the dissociation of the former yields a charged surface. The governing equation, which was solved numerically in a previous study for spherical geometry (Hsu et al. J. Phys. Chem. B 2003, 107, 14429), is solved analytically in this study for planar, cylindrical, and spherical geometries. The analytical results obtained are readily applicable for the evaluation of the spatial distributions of counterions inside a lipid structure. We show that if the linear size of a reverse micelle is fixed, the degree of dissociation of the surfactant layer follows the order planar > cylindrical > spherical.     

Contact angles of drops on curved superhydrophobic surfaces

Superhydrophobic surfaces have contact angles that exceed 150 degrees and are known to reduce surface fouling, protect surfaces, and improve liquid-liquid separations. Electrospun sub-micron fiber mats can perform as superhydrophobic surfaces. Superhydrophobic behavior is typically measured on planar surfaces, whereas applications may require curved surfaces. This paper discuses the measurement of water contact angles of fiber mats formed on cylindrical surfaces to create superhydrophobic behavior on curved surfaces. Equations are derived that relate the radius of curvature of spherical and cylindrical surfaces and drop size to the observed contact angle on the curved surfaces. Calculations from the equations agree well with experimental observations on spherical surfaces reported in literature and on cylindrical surfaces created in our lab.     

Solvation forces between silica bodies in supercritical carbon dioxide

We report Monte Carlo simulations of the solvation pressure between two planar surfaces, which represent the interface of spherical silica nanoparticles in supercritical carbon dioxide. Carbon dioxide (CO2) was modeled as an atomistic dumbbell or a spherical Lennard-Jones particle. The interaction between CO2 molecules and silica surfaces was characterized by the standard Steele potential with energetic heterogeneities representing the hydrogen bonds. The parameters for the solid-fluid interaction potentials were obtained by fitting our simulations to the experimental isotherms of CO2 sorption on mesoporous siliceous materials. We studied the dependence of the solvation force on the distance between planar silica surfaces at T = 318 K, at equilibrium bulk pressures p(bulk) ranging from 69 to 200 atm. At 69 atm, we observed a long-range attraction between the two surfaces, and it vanished when the pressure was increased to 102 and then 200 atm. The results obtained with different fluid models were consistent with each other. According to our observations, energetic heterogeneities of the surface have negligible influence on the solvation pressure. Using the Derjaguin approximation, we calculated the solvation forces between spherical silica nanoparticles in supercritical CO2 from the solvation pressures between the planar surfaces.     

Weak polyelectrolytes tethered to surfaces: Effect of geometry,acid–base equilibrium and electrical permittivity

The structural and thermodynamical properties of weak polyelectrolytes end-tethered to surfaces of arbitrary geometry are studied using a molecular theory. The theory is based on writing down the free energy functional of the system including all the basic interactions and the explicit acid–base equilibrium for the chargeable groups of the polymer. The theory explicitly includes the size, shape, conformations, and charge distribution of all the molecular species. The electrostatic interactions include a density-dependent dielectric function, modeled with the Maxwell–Garnett mixing formula, to account for the composition-dependent permittivity. The minimization of the free energy leads to the distribution of all molecular species and their dependence on bulk pH and salt concentration. We apply the theory to polymer chains end-tethered to planar, cylindrical, and spherical surfaces. The radius of the curved surfaces is small to enhance the curvature effect. We find that when the grafting surfaces are uncharged, the approximation of a constant dielectric function works very well for both structural and thermodynamic properties. The structure of weak polyelectrolytes tethered on cylindrical and spherical surfaces is different from that of polymers tethered on planar surfaces due to the available volume as a function of the distance from the surface. Specifically, the degree of dissociation increases with increasing curvature of the surface. This is a manifestation of the coupling between the local density of protons, counterions, and polymer segments. The results can be interpreted in terms of the local Le Chatelier principle for the acid–base equilibrium, with proper account of the three local contributions: counterions, protons, and chargeable groups. We find that one can achieve local changes of pH between one to two units within 1–2 nm. The thickness of the tethered layers as a function of bulk pH shows a large increase when the pH is equal to the bulk pK. However, the variation with salt concentration is different for the different geometries. The largest swelling is found for cylindrical surfaces. The predictions from scaling theories of a maximum in the thickness of the film as a function of salt concentration is found for planar films, but not for curved surfaces. Finally, the interactions between cylinders with tethered polyelectrolytes is very different from the equivalent planar surfaces. These results are important for the interpretation of force measurements with nanoscale AFM tips. The implications of the results for the rational design of responsive tethered polymer layers is discussed together with the limitations of the theoretical approach. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 2638–2662, 2006     

Electrical potential in a cylindrical double layer: a functional theory approach

Employing an iterative method in functional theory, the electrical potential distribution for the case of a cylindrical surface is solved. Although the analytical result derived is of an iterative nature, the second-order solution is found to be sufficiently accurate under conditions of practical significance. For the case of constant surface potential, the radius and the surface potential of a cylindrical surface can be estimated based on the extreme of the electrical potential distribution. The effects of the key parameters, including the number and the valence of the ions on a surface, the length of a particle, the relative permittivity of the liquid phase, the temperature, and the concentration of electrolyte on the surface potential, are examined. The general behavior of these effects is similar to that for a spherical surface, except that the surface potential of a cylindrical surface is independent of the electrolyte concentration. The present approach is also applicable to the case where a cylindrical surface remains at a constant charge density.     

Electrical potentials of two identical particles with fixed surface charge density in a salt-free medium

The electrical potentials of two identical planar, cylindrical, and spherical particles immersed in a salt-free dispersion are solved analytically by a perturbation approach for the case of constant surface charge density. The system under consideration simulates, for example, micelles, where the ionic species in the liquid phase come mainly from the dissociation of the functional groups on the droplet surface. We show that for planar particles, the present zero-order perturbation solution is exact, and for cylindrical and spherical particles, the first-order perturbation solution provides sufficiently accurate results, with an averaged percentage deviation on the order of 1% under typical conditions. In general, the higher the surface charge density, the higher the valence of counterions, the smaller the separation distance between two particles, and the smaller the curvature of particle surface, the better the performance of the perturbation solution.     

Extension of the Steele 10-4-3 potential for adsorption calculations in cylindrical, spherical, and other pore geometries

Simplified fluid-substrate interaction models derived from the Lennard-Jones potential are widely used in the simulation of gas physisorption phenomena. In this paper, we reinterpret the well known Steele 10-4-3 potential for a gas molecule interacting with a planar surface, and use the resultant scheme to derive new potentials for cylindrical and spherical pore geometries. These new potentials correctly recover the Steele result in the limit of infinite pore radius, a useful improvement over existing models. We demonstrate the new cylindrical Steele 10-4-3 potential in calculations of argon adsorption via fluid density functional theory. This potential yields markedly different adsorption behavior than existing cylindrical potentials, which follow from small but significant differences in both the strength and the shape of the fluid-surface interaction. These differences cannot be fully reconciled simply by reparameterizing (scaling) the existing models; the new potential is more realistic in design, and is especially to be preferred in studies where comparison with planar substrates is made. Finally, we discuss extensions of this approach to more complicated pore geometries, yielding a family of Steele-like potentials that all satisfy the correct planar limit.     

Interaction thresholds for adsorption of quantum gases on surfaces and within pores of various shapes

Adsorption within pores and on surfaces occurs because of the attractive potential provided by the adsorbent. If the attraction is too weak, however, adsorption does not occur to any significant extent. This paper evaluates the criterion for such adsorption, at zero temperature, of the quantum gases 4He and H2. This criterion is expressed as a relationship between a threshold value of the well-depth (D) of the adsorption potential (on a semi-infinite planar surface) and the hard-core diameter (sigma) of the gas-surface pair potential. Six geometries are considered, of which two result in two-dimensional (2D) adsorbed phases, two result in one-dimensional (1D) phases, and two result in zero-dimensional phases. These are monolayer films on semi-infinite substrates or within a slit pore, linear or axial phases within cylindrical pores (within bulk solids) or cylindrical tubes, and single-particle adsorption within spherical pores or hollow spherical cavities, respectively. The criteria for film adsorption are consistent with analogous criteria for film wetting to occur, evaluated with a simple thermodynamic model.     

Critical coagulation concentration of a salt-free colloidal dispersion

Both exact and approximate analytical solutions of the Poisson-Boltzmann equation for two planar, parallel surfaces are derived for the case when a dispersion medium contains counterions only, and the results obtained are used to evaluate the critical coagulation concentration of a spherical dispersion. A correction factor, which is a function of the valence of counterions, the surface potential of a particle, and the potential on the midplane between two particles at the onset of coagulation, is derived to modify the classic Schulze-Hardy rule for the dependence of the critical coagulation concentration on the valence of counterions. The correction factor is found to increase with the increase in the valence of counterions and/or with the increase in the surface potential. However, it approaches a constant value of 0.8390 if the surface potential is sufficiently high.     

Counterion distributions and effective interactions of spherical polyelectrolyte brushes

We consider the electrosteric repulsion of colloidal particles whose surface carries a dense layer of long polyelectrolyte chains (spherical polyelectrolyte brushes). The theory of electrosteric repulsion of star polyelectrolytes developed recently is augmented to include particles with a finite core radius. It is shown that most of the counterions are confined within the brush layer. The strong osmotic pressure thus created within the brush layer dominates the repulsive interaction between two such particles. Because of this the pair interaction potential between spherical polyelectrolyte brushes can be given in terms of an analytic expression. The theoretical predictions are compared with available experimental data and semi-quantitative agreement between the two is found.     

Introducing interacting diffuse layers in TLM calculations: a reappraisal of the influence of the pore size on the swelling pressure and the osmotic efficiency of compacted bentonites

The truncation of the Gouy-Chapman diffuse part in compacted clay-rocks and bentonite is introduced into the electrical triple-layer model (TLM) recently developed by P. Leroy and A. Revil [J. Colloid Interface Sci. 270 (2004) 371]. The new model is used to explain the dependence of the osmotic efficiency and the swelling pressure as functions of the mean pore size of the medium, determined from the porosity and the specific surface. The truncation of the diffuse layer introduces a new variable in the system of equations to be solved, the electrical potential at the midplane between adjacent charged surfaces. This new variable is evaluated through a Taylor expansion of the electrical potential. The present model is able to capture the variation of the osmotic efficiency and the swelling pressure with the mean pore size. The partition of counterions between the Stern layer and the diffuse layer as a function of the pore size calculated by the TLM also shows a good consistency with the model. This implies that more than 90% of the counterions are located in the Stern layer.     

Study on osmotic pressure of non-ionic and ionic surfactant solutions in the micellar and microemulsion regions

To investigate the osmotic pressure of non-ionic and ionic surfactant solutions in the micellar and microemulsion regions, a potential of mean force including hard-core repulsion, van der Waals attraction and electric double layer repulsion is proposed to describe the interactions between micelles and between microemulsions. Both van der Waals attraction and electric double layer repulsion are represented using Yukawa tails. The explicit analytical expression of osmotic pressure derived from the first-order mean spherical approximation is implemented by accounting for the Donnan membrane effect. The proposed theory has been applied to micelle solutions of the non-ionic surfactant, n-dodecyl hexaoxyethylene monoether, the cationic surfactant, cetylpyridinium chloride, the anionic surfactant, sodium dodecyl sulfate, and spherical oil-in-water microemulsion system. Successful comparison is made between the proposed theory and the experimental osmotic pressure data for the studied surfactant solutions. Theoretical results show that the long-range electric double layer repulsion dramatically influences the osmotic pressure of both cationic and anionic surfactant solutions in the micellar region. The regressed model parameters such as effective micelle diameter, the mean aggregation number and effective micellar charge are in good agreement with those from static light scattering studies in the literature.     

Osmotic stress tensor in charge stabilized colloidal crystals

Osmotic stress tensor is introduced to describe the osmotic pressure in colloidal crystals within the framework of the theory of the Poisson-Boltzmann equation. The osmotic stress tensor is related to the fundamental stress tensor, which is associated with the Poisson-Boltzmann equation. It is shown that the osmotic stress tensor can be determined for colloidal crystals with arbitrary structures, as well as for media that are described by cell models. The general results are exemplified by spherical and cylindrical cell models.     

Stern potential near gold and chromium surfaces in contact with aqueous and alcoholic ammonium salt solutions

The repulsive force between two solid surfaces immersed in electrolyte solution is created by overlapping of the electrical double layers adjacent to the solid surfaces when the separation is large enough so that solvation forces can be neglected. When the repulsive force is known, one can calculate the Sternpotential at the solid surfaces. The repulsive interactions in electrolyte solutions were directly measured between two high-grade polished fused silica plates, one planar, the other spherically curved. The surfaces were covered with thin gold or chromium layers and separated by aqueous and alcoholic solutions of NH₄Cl, NH₄Br and NH₄SCN in concentrations 10^–5–10^–2 M. The absolute value of the negative Stern potential derived from the repulsive forces increased when regarding the anions in the order

Electrostatic double-layer interaction between spherical particles inside a rough capillary

Predictions of electrostatic double-layer interaction forces between two similarly charged spherical colloidal particles inside an infinitely long "rough" capillary are presented. A simple model of a rough cylindrical surface is proposed, which assumes the capillary wall to be a periodic function of axial position. The periodic roughness of the wall is characterized by the wavelength and amplitude of the undulations. The electrostatic double-layer interaction force between two spherical particles located axially inside this rough capillary is determined by solving the nonlinear Poisson-Boltzmann equation employing finite element analysis. The effect of surface roughness of the cylindrical enclosure on the interaction force between two particles is extensively studied on the basis of this model. The simulations are carried out for dimensionless amplitudes (amplitude/particle radii) ranging from 0.05 to 0.15 and scaled wavelengths (wavelength/particle radii) ranging from 0.4 to 4.0. The interaction force between the particles is significantly modified by the proximity of the rough capillary wall. Generally, the interaction force for rough capillaries oscillates around the corresponding interaction force in a smooth capillary depending on the magnitudes of the scaled amplitude and wavelength of the roughness. The influence of roughness on the electrostatic interactions becomes more pronounced when the surface potential of the cylinder wall is different from the sphere surface potentials. When the cylinder and the particle surfaces have large potential differences, the axial force experienced by a particle is dominated by the capillary roughness. There are dramatic oscillations of the force, which alternately becomes repulsive and attractive as the particle moves from the crest to the trough of the rough capillary wall. These results suggest that manipulation of colloidal particles in narrow microchannels may be subject to significant force variations owing to the roughness inherent in microfabricated channels etched on metal films.     

Water at hydrophobic substrates: curvature, pressure, and temperature effects

We studied the water density profile close to spherical and planar hydrophobic objects using molecular dynamics (MD) simulations. For normal pressure and room temperature, the depletion layer thickness of a planar substrate is approximately 2.5 Angstroms. Even for quite large spherical solutes with a radius of R = 18 Angstroms, the depletion layer thickness is reduced by 30%, which shows that substrate curvature and roughness is an experimentally important factor. Rising temperature leads to a substantial increase of the depletion layer thickness. The compressibility of the depletion layer is found to be surprisingly small and only approximately 5 times higher than that of bulk water. A high electrostatic surface potential of 0.5 V is found, which presumably plays an important role in the presence of charged solutes, since it can promote adsorption into the interfacial layer.