Electrostatic interaction of a spherical particle in the vicinity of a circular orifice

The electrostatic interaction of a charged spherical particle in the vicinity of an orifice plane has been investigated in this paper.The particle can creep along the axis of the orifice and is immersed in a bulk electrolyte.By solving the Poisson-Boltzmann problem,we have obtained the effective electrostatic interaction for several values of reduced orifice radius h,including the cases of h > 1,h = 1 and h < 1.Two kinds of boundary conditions of the orifice plane are considered.One is the constant potential model corresponding to a conducting plane,the other is the constant charge model.In the constant potential model,there is an electrostatic attraction between the particle and the orifice plane when they get close to each other,while there is a pure electrostatic repulsion in the constant charge model.The interactions in both boundary models are sensitive to the parameters of the reduced orifice radius,the reduced particle-orifice distance,surface charge densities of the particle and orifice plane,and the reduced Debye screen constant corresponding to the salt-ion concentration and ion valence.     

Dipolar Poisson-Boltzmann equation: ions and dipoles close to charge interfaces

We present an extension to the Poisson-Boltzmann model where the dipolar features of solvent molecules are taken explicitly into account. The formulation is derived at mean-field level and can be extended to any order in a systematic expansion. It is applied to a two-plate system with oppositely charged surfaces. The ion distribution and profiles in the dipolar order parameter are calculated and can result in a large correction to the interplate pressure.     

A correlation-hole approach to the electric double layer with counter-ions only

ABSTRACT

We study a classical system of identically charged counter-ions near a planar wall carrying a uniform surface charge density. The equilibrium statistical mechanics of the system depends on a single dimensionless coupling parameter. A new self-consistent theory of the correlation-hole type is proposed which leads to a modified Poisson–Boltzmann integral equation for the density profile, convenient for analytical progress and straightforward to solve numerically. The exact density profiles are recovered in the limits of weak and strong couplings. In contrast to previous theoretical attempts of the test-charge family, the density profiles fulfil the contact-value theorem at all values of the coupling constant and exhibit the mean-field decay at asymptotically large distances from the wall, as expected. We furthermore show that the density corrections at large couplings exhibit the proper dependence on coupling parameter and distance to the charged wall. The numerical results for intermediate values of the coupling provide accurate density profiles which are in good agreement with those obtained by Monte Carlo simulations. The crossover to mean-field behaviour at large distance is studied in detail.     

Classical charged fluids at equilibrium near an interface: Exact analytical density profiles and surface tension

The structure of equilibrium density profiles in an electrolyte in the vicinity of an interface with an insulating or conductive medium is of crucial importance in chemical physics and colloidal science. The Coulomb interaction is responsible for screening effects, and in dilute solutions the latter effects give rise to universal leading corrections to nonideality, which distinguish electrolyte from nonelectrolyte solutions. An example is provided by the excess surface tension for an air-water interface, which is determined by the excess particle density, and which was first calculated by Onsager and Samaras. Because of the discrepancy between the dielectric constants on both sides of the interface, every charge in the electrolyte interacts with an electrostatic image, and the Boltzmann factor associated with the corresponding self-energy has an essential singularity over the length scalel from the wall. Besides Coulomb interactions, short-range repulsions must be taken into account in order to prevent the collapse between charges with opposite signs or between each charge and its image when the solvent dielectric constant is lower than that of the continuous medium on the other side of the interface. For a dilute and weaklycoupled electrolyte,l is negligible with respect to the bulk Debye screening length ξ_D. In the framework of the grand-canonical ensemble, systematic partial resummations in Mayer diagrammatics allow one to exhibit that, in this regime, the exact density profiles at leading order are the same as if they were calculated in a partially-linearized mean-field theory, where the screened pair interaction obeys an inhomogeneous Debye equation. In the latter equation the effective screening length depends on the distancex from the interface: it varies very fast over the lengthl and tends to its bulk value over a few ξDs. The equation can be solved iteratively at any distancex, and the exact density profiles are calculated analytically up to first order in the coupling parameter l/ξ_D. They show the interplay between three effects: (1) the geometric repulsion from the interface associated with the deformation of screening clouds, (2) the polarization effects described by the images on the other side of the interface, (3) the interaction between each charge and the potential drop created by the electric layer which appears as soon as the fluid has not a charge-symmetric composition. Moreover, the expressions allow us to go beyond Onsager-Samaras theory: the surface tension is calculated for charge-asymmetric electrolytes and for any value of the ratio between the dielectric constants on both sides of the interface. Similar diagrammatic techniques also allow one to investigate the charge renormalization in the dipolar effective pair interaction along the interface with an insulating medium.     

Towards a microscopic theory of wetting by ionic solutions. I. Surface properties of the semi-primitive model

As a first step towards a density functional theory (DFT) of wetting by ionic solutions we examine the density profiles of ions and solvent molecules confined near a charged wall, or between two walls, and the corresponding interfacial properties, including adsorption, surface tension, solvation force and electrostatic properties, within the semi-primitive model (SPM) of solutions made up of hard sphere solvent particles and charged hard spheres. Both monovalent and divalent cations with species-dependent diameters are considered. The density functional includes the best available Rosenfeld hard-sphere functional, as well as mean-field and Coulomb correlation contributions. The simpler mean-field functional is found to be adequate, at least for monovalent ions. The size differences lead to an interesting ‘fine structure’ of the density and charge density profiles. Cohesive interactions between all species are shown to lead to significant changes in the density profiles.     

The structure of the interface in the solvent-mediated interaction of dipolar surfaces

Summary Interaction of two dipolar surfaces separated by a polar medium is considered within the framework of nonlocal electrostatics. The dipolar-surface layers are modelled as regular lattices with fixed orientation of dipoles which are immersed into the solvent; solvent response is characterized by nonlocal dielectric function. The model is elaborated in order to reveal the role of the dipolar-layer discreteness in the electric field produced by one surface and the interaction between two surfaces (which gives rise to the so-called ?hydration? or ?structural? force acting between mineral surfaces and phospholipid bilayers). The discreteness effects are present only for commensurate lattices. Their special mutual arrangement then may lead to considerable reduction of structural forces,viz. the usual repulsion regime may change at short distances to attraction. Conditions are considered when repulsion is entirely replaced by attraction,i.e. the ?hydration barrier? disappears. In appended note we discuss the role of solvation of surface dipolar groups. We propose an explanation of why two modes of decay (one with oscillative fine structure) may be present in the dependence of the force upon distance, if the surface dipolar groups are immersed deep enough in the solvent, and how the long-range oscillative mode disappears when the surface is but weakly solvated. To speed up publication, the authors of this paper have agreed to not receive the proofs for correction.     

Electrostatic image effects for counterions between charged planar walls

We study the effect of dielectric inhomogeneities on the interaction between two planparallel charged surfaces with oppositely charged mobile charges in between. The dielectric constant between the surfaces is assumed to be different from the dielectric constant of the two semiinfinite regions bounded by the surfaces, giving rise to electrostatic image interactions. We show that on the weak-coupling level the image charge effects are generally small, making their mark only in the second-order fluctuation term. However, in the strong-coupling limit, the image effects are large and fundamental. They modify the interactions between the two surfaces in an essential way. Our calculations are particularly useful in the regime of parameters where computer simulations would be difficult and extremely time consuming due to the complicated nature of the long-range image potentials.     

Spin relaxation measurements of electrostatic bias in intermolecular exploration

We utilize the paramagnetic contribution to proton spin-lattice relaxation rate constants induced by freely diffusing charged paramagnetic centers to investigate the effect of charge on the intermolecular exploration of a protein by the small molecule. The proton NMR spectrum provided 255 resolved resonances that report how the explorer molecule local concentration varies with position on the surface. The measurements integrate over local dielectric constant variations, and, in principle, provide an experimental characterization of the surface free energy sampling biases introduced by the charge distribution on the protein. The experimental results for ribonuclease A obtained using positive, neutral, and negatively charged small nitroxide radicals are qualitatively similar to those expected from electrostatic calculations. However, while systematic electrostatic trends are apparent, the three different combinations of the data sets do not yield internally consistent values for the electrostatic contribution to the intermolecular free energy. We attribute this failure to the weakness of the electrostatic sampling bias for charged nitroxides in water and local variations in effective translational diffusion constant at the water-protein interface, which enters the nuclear spin relaxation equations for the nitroxide-proton dipolar coupling.     

Static polarizability effects on counterion distributions near charged dielectric surfaces: A coarse-grained Molecular Dynamics study employing the Drude model

Coarse-grained implicit solvent Molecular Dynamics (MD) simulations have been used to investigate the structure of the vicinal layer of polarizable counterions close to a charged interface. The classical Drude oscillator model was implemented to describe the static excess polarizability of the ions. The electrostatic layer correction with image charges (ELCIC) method was used to include the effects of the dielectric discontinuity between the aqueous solution and the bounding interfaces for the calculation of the electrostatic interactions. Cases with one or two charged bounding interfaces were investigated. The counterion density profile in the vicinity of the interfaces with different surface charge values was found to depend on the ionic polarizability. Ionic polarization effects are found to be relevant for ions with high excess polarizability near surfaces with high surface charge.     

Spectral broadening in suspensions of metallic spheres

We study the effective dielectric constant of dilute suspensions of metallic spheres in a transparent background, e.g. glass or water. In the low density limit the effective dielectric constant as a function of frequency exhibits a strong absorption peak corresponding to the dipolar surface plasmon of isolated spheres. We investigate how the absorption peak is modified at higher density due to electrostatic interactions between spheres. We find appreciable deviations from the Clausius-Mossotti formula which may be seen most clearly in a Cole-Cole plot of the effective dielectric constant.     

Electrophoretic mobility without charge driven by polarisation of the nanoparticle–water interface

Polarisation of the interface, spontaneously occurring when water is in contact with hydrophobic solutes or air, couples with the uniform external field to produce a non-zero force acting on a suspended particle. This force exists even in the absence of a net particle charge, and its direction is affected by the first-order, dipolar and the second-order, qudrupolar orientational order parameters of the interfacial water. The quadrupolar polarisation gives rise to an effectively negative charge. The corresponding surface charge density is inversely proportional to the area of the shear surface. As a result, the overall contribution from the quadrupolar polarisation to the particle mobility becomes negligible compared to experimentally reported values for particles exceeding a few nanometres in size. In contrast, the contribution of the dipolar order of the interface to the effective surface charge scales inversely with the particle size and dominates the zero-charge mobility of submicron particles. The corresponding electrokinetic charge is determined by the preferential orientation of interfacial dipoles relative to the surface normal.     

Polyelectrolyte adsorption

The problem of charged polymer chains (polyelectrolytes) as they adsorb on a planar surface is addressed theoretically. We review the basic mechanisms and theory underlying polyelectrolyte adsorption on a single surface in two situations: adsorption of a single charged chain, and adsorption from a bulk solution in θ solvent conditions. The behavior of flexible and semi-rigid chains is discussed separately and is expressed as function of the polymer and surface charges, ionic strength of the solution and polymer bulk concentration. We mainly review mean-field results and briefly comment about fluctuation effects. The phenomenon of polyelectrolyte adsorption on a planar surface as presented here is of relevance to the stabilization of colloidal suspensions. In this respect we also mention calculations of the inter-plate force between two planar surfaces in presence of polyelectrolyte. Finally, we comment on the problem of charge overcompensation and its implication to multi-layers formation of alternating positive and negative polyelectrolytes on planar surfaces and colloidal particles.     

Dielectric alignment of a non-polar molecule substituting a polar molecule in a lattice of molecular dipoles

The dielectric alignment of a dilute non-polar component in a polar solvent may be determined by NMR spectroscopy. In this paper a series expansion of the alignment up to second order terms in the dipolar interaction is presented for a non-polar molecule replacing a single polar molecule in a rigid lattice of molecular dipoles. Assuming isotropic polarizabilities for all molecules Van Vleck earlier applied the lattice model in a theory of the dielectric constant of a dilute solution.In the present calculation the polarizabilities of the non-polar molecule and the dipoles are assumed to be anisotropic and isotropic, respectively. The anisotropy of the non-polar molecule is relevant, since the value of the alignment is zero for the isotropic case. Different results are obtained for cubic lattices and for a lattice points are uniformly distributed.In absence of a rigid dipole moment of the solvent molecules the latter formula may be compared with that implicity deduced by Buckingham for the Kerr effect. In the limit of strong dipoles the orientation imposed by the rigid dipole moments appears to be the major contribution to the alignment of the non-polar molecule.     

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.     

Spin Exchange Between Charged Paramagnetic Particles in Dilute Solutions

Kinetic equations for the spin density matrix which take into account binary collisions and a method of calculating the spin exchange effective radius have been generalized to the case of dilute solutions of charged paramagnetic particles. The effective radius of the spin exchange and rate constant of the bimolecular spin exchange between charged paramagnetic particles in solutions have been calculated numerically. Calculations have been performed under the assumption that the exchange interaction is isotropic and decays exponentially with the increase in the distance between radicals, and the solution has a given dielectric permittivity and Debye screening radius. Dependences of the spin exchange rate constant on the mutual diffusion coefficient, exchange and electrostatic interactions parameters have been found numerically. The theory has been applied to experimental results taken from the literature. The rate constant of the spin exchange between radicals of like charge found from the experiment and calculated within the developed theory are in good qualitative agreement .     

An integral equation theory for the dense dipolar hard-sphere fluid

An integral equation theory based upon a simple linearization of the hypernetted-chain closure approximation is solved numerically for dipolar hard spheres. The pair-correlation function and thermodynamic properties are found to be in excellent agreement with Monte Carlo results. The dielectric constant in this approximation is larger than that predicted by previous theories.     

Poisson-Nernst-Planck Equations for Simulating Biomolecular Diffusion-Reaction Processes I: Finite Element Solutions

In this paper we developed accurate finite element methods for solving 3-D Poisson-Nernst-Planck (PNP) equations with singular permanent charges for electrodiffusion in solvated biomolecular systems. The electrostatic Poisson equation was defined in the biomolecules and in the solvent, while the Nernst-Planck equation was defined only in the solvent. We applied a stable regularization scheme to remove the singular component of the electrostatic potential induced by the permanent charges inside biomolecules, and formulated regular, well-posed PNP equations. An inexact-Newton method was used to solve the coupled nonlinear elliptic equations for the steady problems; while an Adams-Bashforth-Crank-Nicolson method was devised for time integration for the unsteady electrodiffusion. We numerically investigated the conditioning of the stiffness matrices for the finite element approximations of the two formulations of the Nernst-Planck equation, and theoretically proved that the transformed formulation is always associated with an ill-conditioned stiffness matrix. We also studied the electroneutrality of the solution and its relation with the boundary conditions on the molecular surface, and concluded that a large net charge concentration is always present near the molecular surface due to the presence of multiple species of charged particles in the solution. The numerical methods are shown to be accurate and stable by various test problems, and are applicable to real large-scale biophysical electrodiffusion problems.     

Simulations of counterions at charged plates

Using Monte Carlo simulations, we study the counterion distribution close to planar charged walls in two geometries: i) when only one charged wall is present and the counterions are confined to one half-space, and ii) when the counterions are confined between two equally charged walls. In both cases the surface charge is smeared out and the dielectric constant is the same everywhere. We obtain the counterion density profile and compare it with both the Poisson-Boltzmann theory (asymptotically exact in the limit of weak coupling, i.e. low surface charge, high temperature and low counterion valence) and the strong-coupling theory (valid in the opposite limit of high surface charge, low temperature and high counterion valence) and with previously calculated correction terms to both theories for different values of the coupling parameter, thereby establishing the domain of validity of the asymptotic limits. Gaussian corrections to the leading Poisson-Boltzmann behavior (obtained via a systematic loop expansion) in general perform quite poorly: At coupling strengths low enough so that the Gaussian (or one-loop) correction does describe the numerical deviations from the Poisson-Boltzmann result correctly, the leading Poisson-Boltzmann term by itself matches the data within high accuracy. This reflects the slow convergence of the loop expansion. For a single charged plane, the counterion pair correlation function indicates a behavioral change from a three-dimensional, weakly correlated counterion distribution (at low coupling) to a two-dimensional, strongly correlated counterion distribution (at high coupling), which is paralleled by the specific-heat capacity which displays a rounded hump at intermediate coupling strengths. For the case of counterions confined between two equally charged walls, we analyze the inter-wall pressure and establish the complete phase diagram, featuring attraction between the walls for large enough coupling strength and at intermediate wall separation. Depending on the thermodynamic ensemble, the phase diagram exhibits a discontinuous transition where the inter-wall distance jumps to infinity (in the absence of a chemical potential coupling to the inter-wall distance, as for charged lamellae in excess solvent) or a critical point where two coexisting states with different inter-wall distance become indistinguishable (in the presence of a chemical potential, as for charged lamellae with a finite fixed solvent fraction). The attractive pressure decays with the inter-wall distance as an inverse cube, similar to analytic predictions, although the amplitude differs by an order of magnitude from previous theoretical results. Finally, we discuss in detail our simulation methods and compare the finite-size scaling behavior of different boundary conditions (periodic, minimal image and open). Received 6 November 2001     

Ground State and Charge Renormalization in a Nonlinear Model of Relativistic Atoms

We study the reduced Bogoliubov-Dirac-Fock (BDF) energy which allows to describe relativistic electrons interacting with the Dirac sea, in an external electrostatic potential. The model can be seen as a mean-field approximation of Quantum Electrodynamics (QED) where photons and the so-called exchange term are neglected. A state of the system is described by its one-body density matrix, an infinite rank self-adjoint operator which is a compact perturbation of the negative spectral projector of the free Dirac operator (the Dirac sea). We study the minimization of the reduced BDF energy under a charge constraint. We prove the existence of minimizers for a large range of values of the charge, and any positive value of the coupling constant α. Our result covers neutral and positively charged molecules, provided that the positive charge is not large enough to create electron-positron pairs. We also prove that the density of any minimizer is an L ¹ function and compute the effective charge of the system, recovering the usual renormalization of charge: the physical coupling constant is related to α by the formula α_phys ≃ α(1 + 2α/(3π) log Λ)⁻¹, where Λ is the ultraviolet cut-off. We eventually prove an estimate on the highest number of electrons which can be bound by a nucleus of charge Z. In the nonrelativistic limit, we obtain that this number is  ≤  2Z, recovering a result of Lieb. This work is based on a series of papers by Hainzl, Lewin, Séré and Solovej on the mean-field approximation of no-photon QED.