Numerical Solution of the Electrokinetic Equations for Multi-ionic Electrolytes Including Different Ionic Size Related Effects

Numerical study supplemented with simplified model on electrophoresis of a hydrophobic colloid incorporating finite ion size effects and ion-solvent interactions

We consider a modified electrokinetic model to study the electrophoresis of a hydrophobic particle by considering the finite sized ions. The mathematical model adopted in this study incorporates the ion steric repulsion, ion-solvent interactions as well as Maxwell stress on the electrolyte. The dielectric permittivity and viscosity of the electrolyte is considered to vary with the local ionic volume fraction. Based on this modified model for the electrokinetics we have analyzed the electrophoresis in a single as well as mixture of electrolytes of monovalent and non- z : z $z:z$ electrolytes. The dependence of viscosity on local ionic volume fraction modifies the hydrodynamic drag as well as diffusivity of ions, which are ignored in existing studies on electrophoresis. A simplified model for electrophoresis of a hydrophobic particle incorporating the ion steric repulsion and ion-solvent interactions is developed based on the first-order perturbation on applied electric field. This simplified model is established to be efficient for a Debye layer thinner than the particle size and a smaller range of slip length. This model can be implemented for any number of ionic species as well as non- z : z $z:z$ electrolytes. It is established that the ion steric interactions and dielectric decrement creates a counterion saturation in the Debye layer leading to an enhanced mobility compared to the standard model. However, experimental data for non-dilute cases often under predicts the theoretically determined mobility. The present modified model fills this lacuna and demonstrate that the consideration of finite ion size modifies the medium viscosity and hence, ionic mobility, which in combination lowers the mobility value.

Electrokinetic and dielectric response of a concentrated salt-free colloid: Different approaches to counterion finite-size effects

In the present work, a general model is developed for the electrokinetics and dielectric response of a concentrated salt-free colloid that takes into account the finite size of the counterions released by the particles to the solution. The effects associated with the counterion finite size have been addressed using a hard-sphere model approach elaborated by Carnahan and Starling [N. F. Carnahan and K. E. Starling, Equation of state for nonattracting rigid spheres, J. Chem. Phys. 51, 635 (1969)0021-960610.1063/1.1672048]. A more simple description of the finite size of the counterions based on that by Bikerman has also been considered for comparison. The studies carried out in this work include predictions on the effect of the finite counterion size on the equilibrium properties of the colloid and its electrokinetic and dielectric response when it is subjected to constant or alternating electric fields. The results show how important the counterion finite-size effects are for most of the electrokinetic and dielectric properties of highly charged and concentrated colloids, mainly for the static and dynamic electrophoretic mobilities. Furthermore, new insights are provided on the counterion condensation effect when counterions are allowed to have finite size. Focus is placed on the changes undergone by their concentration in the condensation layer for low-salt and highly charged colloids.

Applicability of the linearized Poisson-Boltzmann theory to contact angle problems and application to the carbon dioxide-brine-solid systems

In colloidal science and bioelectrostatics, the linear Poisson Boltzmann equation (LPBE) has been used extensively for the calculation of potential and surface charge density. Its fundamental assumption rests on the premises of low surface potential. In the geological sequestration of carbon dioxide in saline aquifers, very low pH conditions coupled with adsorption induced reduction of surface charge density result in low pH conditions that fit into the LPB theory. In this work, the Gouy–Chapman model of the electrical double layer has been employed in addition to the LPBE theory to develop a contact angle model that is a second-degree polynomial in pH. Our model contains the point of zero charge pH of solid surface. To render the model applicable to heterogeneous surfaces, we have further developed a model for the effective value of the point of zero charge pH. The point of zero charge pH model when integrated into our model enabled us to determine the point of zero charge pH of sandstone, quartz and mica using literature based experimental data. In this regard, a literature based thermodynamic model was used to calculate carbon dioxide solubility and pH of aqueous solution. Values of point of zero charge pH determined in this paper agree with reported ones. The novelty of our work stems from the fact that we have used the LPB theory in the context of interfacial science completely different from the classical approach, where the focus is on interparticle electrostatics involving colloidal stabilization.

Asymmetric double-layer charging in a cylindrical nanopore under closed confinement

This article presents a physical-mathematical treatment and numerical simulations of electric double layer charging in a closed, finite, and cylindrical nanopore of circular cross section, embedded in a polymeric host with charged walls and sealed at both ends by metal electrodes under an external voltage bias. Modified Poisson-Nernst-Planck equations were used to account for finite ion sizes, subject to an electroneutrality condition. The time evolution of the formation and relaxation of the double layers was explored. Moreover, equilibrium ion distributions and differential capacitance curves were investigated as functions of the pore surface charge density, electrolyte concentration, ion sizes, and pore size. Asymmetric properties of the differential capacitance curves reveal that the structure of the double layer near each electrode is controlled by the charge concentration along the pore surface and by charge asymmetry in the electrolyte. These results carry implications for accurately simulating cylindrical capacitors and electroactuators.