Effects of dielectric discontinuities on two charged plates

On Capacitance and Energy Storage of Supercapacitor with Dielectric Constant Discontinuity

The classical density functional theory (CDFT) is applied to investigate influences of electrode dielectric constant on specific differential capacitance Cd and specific energy storage E of a cylindrical electrode pore electrical double layer. Throughout all calculations the electrode dielectric constant varies from 5, corresponding to a dielectric electrode, to εwr= 10⁸ corresponding to a metal electrode. Main findings are summarized as below. (i): By using a far smaller value of the solution relative dielectric constant εr=10, which matches with the reality of extremely narrow tube, one discloses that a rather high saturation voltage is needed to attain the saturation energy storage in the ultra-small pore. (ii): Use of a realistic low εr=10 value brings two obvious effects. First, influence of bulk electrolyte concentration on the Cd is rather small except when the electrode potential is around the zero charge potential; influence on the E curve is almost unobservable. Second, there remain the Cd and E enhancing effects caused by counter-ion valency rise, but strength of the effects reduces greatly with dropping of the εr value; in contrast, the Cd and E reducing effects coming from the counter-ion size enhancing remain significant enough for the low εr value. (iii) A large value of electrode relative dielectric constant εrw always reduces both the capacitance and energy storage; moreover, the effect of the εrw value gets eventually unobservable for small enough pore when the εrw value is beyond the scope corresponding to dielectric electrode. It is analyzed that the above effects take their rise in the repulsion and attraction on the counter-ions and co-ions caused by the electrode bound charges and a strengthened inter-counter-ion electrostatic repulsion originated in the low εr value.

Charged nanorods at heterogeneously charged surfaces

We study the spatial and orientational distribution of charged nanorods (rodlike counterions) as well as the effective interaction mediated by them between two plane-parallel surfaces that carry fixed (quenched) heterogeneous charge distributions. The nanorods are assumed to have an internal charge distribution, specified by a multivalent monopolar moment and a finite quadrupolar moment, and the quenched surface charge is assumed to be randomly distributed with equal mean and variance on the two surfaces. While equally charged surfaces are known to repel within the traditional mean-field theories, the presence of multivalent counterions has been shown to cause attractive interactions between uniformly charged surfaces due to the prevalence of strong electrostatic couplings that grow rapidly with the counterion valency. We show that the combined effects due to electrostatic correlations (caused by the coupling between the mean surface field and the multivalent, monopolar, charge valency of counterions) as well as the disorder-induced interactions (caused by the coupling between the surface disorder field and the quadrupolar moment of counterions) lead to much stronger attractive interactions between two randomly charged surfaces. The interaction profile turns out to be a nonmonotonic function of the intersurface separation, displaying an attractive minimum at relatively small separations, where the ensuing attraction can exceed the maximum strong-coupling attraction (produced by multivalent monopolar counterions between uniformly charged surfaces) by more than an order of magnitude.

Electrolytes between dielectric charged surfaces: Simulations and theory

We present a simulation method to study electrolyte solutions in a dielectric slab geometry using a modified 3D Ewald summation. The method is fast and easy to implement, allowing us to rapidly resum an infinite series of image charges. In the weak coupling limit, we also develop a mean-field theory which allows us to predict the ionic distribution between the dielectric charged plates. The agreement between both approaches, theoretical and simulational, is very good, validating both methods. Examples of ionic density profiles in the strong electrostatic coupling limit are also presented. Finally, we explore the confinement of charge asymmetric electrolytes between neutral surfaces.

Distribution of counterions and interaction between two similarly charged dielectric slabs: roles of charge discreteness and dielectric inhomogeneity

The distribution of counterions and the electrostatic interaction between two similarly charged dielectric slabs is studied in the strong coupling limit. Dielectric inhomogeneities and discreteness of charge on the slabs have been taken into account. It is found that the amount of dielectric constant difference between the slabs and the environment, and the discreteness of charge on the slabs have opposing effects on the equilibrium distribution of the counterions. At small interslab separations, increasing the amount of dielectric constant difference increases the tendency of the counterions toward the middle of the intersurface space between the slabs and the discreteness of charge pushes them to the surfaces of the slabs. In the limit of point charges, independent of the strength of dielectric inhomogeneity, counterions distribute near the surfaces of the slabs. The interaction between the slabs is attractive at low temperatures and its strength increases with the dielectric constant difference. At room temperature, the slabs may completely attract each other, reach to an equilibrium separation, or have two equilibrium separations with a barrier in between, depending on the system parameters.

Models of membrane electrostatics

Formulas are derived for the electrostatic potential of a charge in or near a membrane modeled as one or more dielectric slabs lying between two semi-infinite dielectrics. One can use these formulas in Monte Carlo codes to compute the distribution of ions near cell membranes more accurately than by using Poisson-Boltzmann theory or its linearized version. Here I use them to discuss the electric field of a uniformly charged membrane, the image charges of an ion, the distribution of salt ions near a charged membrane, the energy of a zwitterion near a lipid slab, and the effect of including the phosphate head groups as thin layers of high electric permittivity.

Simulation of an electrical double layer model with a low dielectric layer between the electrode and the electrolyte

We report Monte Carlo simulation results for double layers of 1:1 and 2:1 electrolytes near an electrode with an inner layer that has a dielectric constant, ε(2), smaller than that of the electrolyte, ε(3). The electrolyte is modeled in the implicit solvent framework (primitive model), while the electrode is a metal electrode in this study (ε(1) → ∞). The charged hard sphere ions are not allowed to enter into the inner layer. We show that the capacitance of the inner layer is C(δ) = ε(0)(ε(2) + ε(3))/2δ, where δ is the thickness of the inner layer. This result is different from that obtained from solutions of the Poisson-Boltzmann equation (ε(0)ε(2)/δ), indicating that interpretation of experimental data with a fitted ε(2) dielectric constant of the inner layer must be done using a different equation. We also show that the properties of the diffuse layer are not independent of the value of ε(2), which is a usual assumption of the Poisson-Boltzmann theory. This is mainly because the repulsive image charges repel both the counterions and the co-ions, while the electrode charge attracts the counterions and repels the co-ions.

Dressed counterions: polyvalent and monovalent ions at charged dielectric interfaces

We investigate the ion distribution and overcharging at charged interfaces with dielectric inhomogeneities in the presence of asymmetric electrolytes containing polyvalent and monovalent ions. We formulate an effective "dressed counterion" approach by integrating out the monovalent salt degrees of freedom and show that it agrees with results of explicit Monte Carlo simulations. We then apply the dressed counterion approach within the framework of the generalized strong-coupling theory, valid for polyvalent ions at low concentrations, which enables an analytical description for salt effects as well as dielectric inhomogeneities in the limit of strong Coulomb interactions. Limitations and applicability of this theory are examined by comparing the results with simulations.

Dressed counterions: strong electrostatic coupling in the presence of salt

We reformulate the theory of strong electrostatic coupling in order to describe an asymmetric electrolyte solution of monovalent salt ions and polyvalent counterions using field-theoretical techniques and Monte Carlo simulations. The theory is based on an asymmetric treatment of the different components of the electrolyte solution. The weak coupling Debye-Hückel approach is used in order to describe the monovalent salt ions while a strong coupling approach is used to tackle the polyvalent counterions. This combined weak-strong coupling approach effectively leads to dressed interactions between polyvalent counterions and thus directly affects the correlation attraction mediated by polyvalent counterions between like-charged objects. The general theory is specifically applied to a system composed of two uniformly charged plane-parallel surfaces in the presence of salt and polyvalent counterions. In the strong coupling limit for polyvalent counterions, the comparison with Monte Carlo simulations shows good agreement for large enough values of the electrostatic coupling parameter. We delineate two limiting laws that in fact encompass all the Monte Carlo data.

Monte Carlo determination of mixed electrolytes next to a planar dielectric interface with different surface charge distributions

Employing canonical ensemble Monte Carlo simulations, we report a calculation of the distribution of small ions next to a planar negatively charged surface in the presence of mixed electrolytes of monovalent and trivalent salt ions within the framework of the primitive model under more realistic hydrated ion size conditions. The effects of surface charge discreteness and dielectric breakdown on charge inversion are discussed based on increasing concentration of both monovalent and trivalent salt. Moreover, a comparison of the simulation results for different discretization models is made along with the case of uniformly distributed charge in terms of the ionic density profiles as well as the integrated charge distribution function. For finite size charged groups located inside the lower dielectric region, a complete equivalence with the case of uniform distribution is observed if the quantities of interest are exclusively analyzed as a function of the distance to the charged interface. With protruding head groups into the aqueous solution, the excluded volume dominates over the correlation effect, therefore the ions are less accumulated in the vicinity of the charged surface, inducing that the onset position of charge inversion experiences an evident shift toward the aqueous environment. Overall, the effect of repulsive image forces on the diffuse double layer structure can be significant at low surface charge density irrespectively of surface charge distributions.

The role of multipoles in counterion-mediated interactions between charged surfaces: strong and weak coupling

We present general arguments for the importance, or lack thereof, of structure in the charge distribution of counterions for counterion-mediated interactions between bounding symmetrically charged surfaces. We show that on the mean field or weak coupling level, the charge quadrupole contributes the lowest order modification to the contact value theorem and thus to the intersurface electrostatic interactions. The image effects are non-existent on the mean field level even with multipoles. On the strong coupling level the quadrupoles and higher order multipoles contribute additional terms to the interaction free energy only in the presence of dielectric inhomogeneities. Without them, the monopole is the only multipole that contributes to the strong coupling electrostatics. We explore the consequences of these statements in all their generality.

One-dimensional counterion gas between charged surfaces: exact results compared with weak- and strong-coupling analyses

We evaluate exactly the statistical integral for an inhomogeneous one-dimensional (1D) counterion-only Coulomb gas between two charged boundaries and from this compute the effective interaction, or disjoining pressure, between the bounding surfaces. Our exact results are compared to the limiting cases of weak and strong couplings which are the same for 1D and three-dimensional (3D) systems. For systems with a large number of counterions it is found that the weak-coupling (mean-field) approximation for the disjoining pressure works perfectly and that fluctuations around the mean-field in 1D are much smaller than in 3D. In the case of few counterions it works less well and strong-coupling approximation performs much better as it takes into account properly the discreteness of the counterion charges.

Weak- and strong-coupling electrostatic interactions between asymmetrically charged planar surfaces

We compare weak- and strong-coupling theory of counterion-mediated electrostatic interactions between two asymmetrically charged plates with extensive Monte Carlo simulations. Analytical results in both weak- and strong-coupling limits compare excellently with simulations in their respective regimes of validity. The system shows a surprisingly rich structure in terms of interactions between the surfaces as well as fundamental qualitative differences in behavior in the weak- and the strong-coupling limits.

Strong-coupling electrostatics in the presence of dielectric inhomogeneities

We study the strong-coupling (SC) interaction between two like-charged membranes of finite thickness embedded in a medium of higher dielectric constant. A generalized SC theory is applied along with extensive Monte Carlo simulations to study the image charge effects induced by multiple dielectric discontinuities in this system. These effects lead to strong counterion crowding in the central region of the intersurface space upon increasing the solvent-membrane dielectric mismatch and change the membrane interactions from attractive to repulsive at small separations. These features agree quantitatively with the SC theory at elevated couplings or dielectric mismatch where the correlation hole around counterions is larger than the thickness of the central counterion layer.