interactions between neutral dielectric bodies are traditionally viewed as being due predominantly to electromagnetic field fluctuations , or equivalently , dipole fluctuations that give rise to van der waals ( vdw ) fluctuation - induced interactions between them @xcite .
these interactions are attractive between identical bodies in vacuum or in a polarizable medium , such as water or an aqueous electrolyte ( or coulomb fluid ) .
they contribute one of the two key ingredients of the classical derjaguin - landau - verwey - overbeek ( dlvo ) theory of colloidal stability , the other one being the mean - field electrostatic interaction , which is repulsive for like - charged colloidal surfaces @xcite .
recent works have , however , highlighted the role of image - induced , ion - depletion effects in this scenario ( see , e.g. , refs . @xcite and references therein ) , leading to depletion of mobile solution ions from the vicinity of dielectric interfaces and , therefore , to an additional attractive force between apposed dielectric boundaries .
this is because most dielectric surfaces in the context of bio- and soft materials have a lower ( static ) dielectric constant than that of water and , therefore , solution ions experience repulsion from their same - sign image charges in the proximity of dielectric boundaries @xcite .
the recent advances in the study of image - induced , ion - depletion effects follow on the trail of the onsager - samaras framework formulated originally in the context of the surface tension problem of electrolytes @xcite
. such non - mean - field effects , which belong to the general class of depletion interactions @xcite , arise due to the discrete nature of mobile ions neglected in the collective mean - field description based on the standard poisson - boltzmann theory @xcite .
the studies of image - induced , ion - depletion effects have been focused exclusively on the case of strictly neutral ( charge - free ) dielectric surfaces . in this case , the ion - depletion interactions can be amplified in the presence of mobile multivalent ions in the solution @xcite due to stronger ion - image repulsions for these ions , even when the multivalent ions are present at small bulk concentrations around just a few mm @xcite .    in this paper
, we revisit the problem of interaction between neutral dielectric surfaces in a coulomb fluid by adding to it a novel feature : we relax the constraint of _ strict _ electroneutrality of surface boundaries , considered so far in the literature , by assuming that the surfaces are neutral only _ on the average _ , while microscopically they carry a _ quenched _ ( fixed ) random distribution of positive and negative charges .
we show that this seemingly simple generalization leads to significant qualitative changes in the distribution of ions and , consequently , also in the effective interactions between dielectric surfaces , especially when the intervening coulomb fluid contains mobile multivalent ions .
disordered charged systems are abundant in soft matter with examples ranging from polycrystalline surfaces with patchy surface potentials @xcite , dielectric contact surfaces @xcite , vapor - deposited amorphous films on solid substrates @xcite , surfactant - coated surfaces @xcite , dna microarrays @xcite , intrinsically disordered proteins @xcite , patchy colloids @xcite and random polyelectrolytes and polyampholytes @xcite .
surface charge disorder in these examples can exhibit highly random distributions as well as patchy and heterogeneous patterns , originating from different sources including specific electronic and/or structural properties of materials involved , surface grafting or adsorption of charged molecules and/or contaminants , synthetic and fabrication processes , etc
. the surface charge disorder can be highly sample specific and , at the same time , can depend strongly on the method of preparation .
it may be set and quenched for each sample ( which is the case of interest in this paper ) , annealed ( in which case the surface charges are mobile and in thermal equilibrium with the rest of the system ) , or partially quenched or partially annealed @xcite .
motivated by these examples , study of charge disorder and , in particular , effective interactions between random charge distributions has witnessed growing attention from theoreticians over the last several years @xcite , as well as from the simulation side where initial steps have been taken to include the effects of charge heterogeneity and disorder @xcite .
disorder effects have been studied extensively in the context of electromagnetic fluctuation - induced interactions between two apposed , randomly charged surfaces in vacuum @xcite , where they have been associated with anomalously long - ranged surface interactions observed in recent experiments @xcite .
they have also been studied in situations where a weakly coupled coulomb fluid , containing , e.g. , monovalent cation and anions , intervenes between the bounding surfaces @xcite .
this however leaves out the case of asymmetric coulomb fluids containing both monovalent ions and multivalent ( counter- ) ions .
these kinds of systems are in fact quite common in experiments within the biological context as in the case of viruses , dna condensates or other charged biopolymer aggregates @xcite and are expected to behave very differently since multivalent ions are known to couple strongly with fixed surface charges .
strong - coupling behavior of multivalent counterions at uniformly charged surfaces or surfaces with regular charge patterns has been studied extensively over the last decade and its connection to exotic phenomena such as like - charged attraction has been throughly discussed ( for recent reviews and a more exhaustive list of references , see refs .
@xcite ) . in disordered charged systems ,
such strong - coupling phenomena have been considered only in a few cases so far @xcite by assuming that bounding surfaces carry a _
finite _ mean surface charge density , to which multivalent counterions can couple strongly , in the same sense as considered in the context of uniformly charged surfaces as noted above .
yet , the presence of charge randomness on bounding surfaces was shown to give rise to novel phenomena such as strong surface attraction of multivalent counterions , characterized by a density profile that diverges at the surface , in clear violation of the contact - value theorem established for uniformly charged surfaces .
this kind of behavior originates from a singular , attractive , single - ion potential , which is created by the randomness in the distribution of surface charges and , as such , depends on the surface charge variance .
consequently , one can also show that the thermal entropy of counterions is diminished upon introducing a finite degree of charge randomness on the boundaries .
in other words , the system becomes more ` ordered ' as a direct outcome of the interplay between the thermal entropy of ions and the configurational entropy of charge randomness , entering through an ensemble average over various realizations of the disordered charge distribution .
this peculiar disorder - induced effect , which stands at odds with what one may expect intuitively , has been referred to as _ antifragility _
@xcite .
our analysis in this paper is focused on yet another facet of the disorder - induced effects by assuming , in a system of two plane - parallel dielectric slabs , that the randomly charged surfaces of the slabs bear _ no net charge_. this eliminates the direct strong coupling between multivalent ions and the mean surface charge and , thereby , also the ensuing strong - coupling interactions , considered in our previous papers @xcite , that would otherwise completely mask the vdw and image - induced , ion - depletion interactions between the slabs in an asymmetric coulomb fluid .
this allows us to address the question of how the presence of surface charge disorder affects the standard picture for the interaction of neutral bodies based on the vdw and ion - depletion effects .
the surface charge disorder has several different implications : first , it contributes a _ repulsive _ interaction between the slabs , which comes from self - interactions of disorder charges on the bounding surfaces with their image charges ; this contribution counteracts the _ attractive _ vdw interaction as discussed thoroughly elsewhere @xcite .
then , as noted above , the individual ions in the coulomb fluid experience an attractive disorder - induced potential , which strongly attracts them toward the bounding surfaces .
this effect counteracts the ion - image repulsions that tend to deplete ions from the slit region between the slabs and form the basis of the image - induced , ion - depletion mechanism for attraction between strictly neutral slabs .
thus , randomly charged bounding surfaces tend to accumulate more mono- and multivalent ions in the slit .
the system , however , responds differently to the increased number of ions : while the osmotic pressure due to monovalent ions increases and even becomes _ repulsive _ , consistent with the standard mean - field picture , the osmotic pressure due to multivalent ions , becomes ever more _ attractive _ !
this behavior is rooted in the combined effect of the surface charge disorder and the presence of mobile multivalent ions , with the latter creating strong inter - surface attractions upon further accumulation in the slit . as a result
, the effective total interaction pressure between randomly charged , net - neutral dielectric surfaces can differ qualitatively from what one expects based on the standard vdw @xcite and image - induced , ion - depletion interactions in the case of strictly neutral surfaces ( see , e.g. , refs . @xcite and references therein ) .
the net effect due to the interplay between disorder , mono- and multivalent - ion contributions yields a qualitatively different behavior for the effective surface - surface interactions , depending on the strength of disorder as quantified by the _ disorder coupling ( or strength ) parameter _ @xcite .
this feature bears some conceptual resemblance to the strong - weak coupling dichotomy that exists for net - charged surfaces , dependent on the _ electrostatic coupling parameter _ in that case @xcite . in weakly disordered systems , corresponding to a small disorder coupling parameter , one then discerns a distinct non - monotonic behavior for the interaction pressure as a function of the separation between the slabs , with a pronounced repulsive hump at intermediate separations ; conversely , in strongly disordered systems , corresponding to a large disorder coupling parameter , the effective interaction pressure becomes strongly attractive , with a range and magnitude larger than that of the vdw or the image - induced , ion - depletion interaction pressure as found between strictly neutral surfaces .
the organization of the paper is as follows : in section [ sec : model ] , we introduce our model and , in section [ sec : general ] , we briefly discuss the theoretical background and present the general results for the two - slab system .
the results of our analysis for the density profile of multivalent ions and the effective interactions between the slabs are presented in section [ sec : results ] .
we conclude our discussion in section [ sec : conclusion ] .     and randomly charged inner surfaces are immersed in a bathing ionic solution of dielectric constant @xmath0 .
the solution contains a mixture of monovalent and multivalent salts .
multivalent ions are confined in the slit region and are shown by large spheres ; monovalent salt anions and cations are shown by small red and blue spheres .
the slabs are assumed to be neutral on the average and their thickness is taken to be infinite.,width=321 ]
our model consists of two plane - parallel dielectric slabs of infinite surface area @xmath1 and dielectric constant @xmath2 with inner surfaces placed perpendicular to the @xmath3 axis at a separation distance of @xmath4 ( see fig .
[ f : schematic ] ) .
the slabs are immersed in an asymmetric coulomb fluid of dielectric constant @xmath0 at room temperature @xmath5 .
the coulomb fluid is a mixture of a monovalent @xmath6 salt of bulk ( reservoir ) concentration @xmath7 and a multivalent @xmath8 salt of bulk concentration @xmath9 with multivalent ions having charge valency of @xmath10 .
the multivalent ions are confined within the slit region @xmath11 , while the monovalent ions are dispersed throughout the space except in the region occupied by the dielectric slabs that are assumed to be impermeable to mobile ions @xcite .
we are interested only in the cases where the slab thickness is much larger than the debye ( or salt ) screening length and , thus , in the calculations to be presented later , we shall assume that the slab thickness is infinite .
the inner surfaces of the slabs are assumed to bear a _ quenched _ , random charge distribution @xmath12 , while they remain electrically _ net - neutral_. we assume that the random surface charges are distributed according to a gaussian probability distribution function , which is determined fully by its two moments , @xmath13 , and @xmath14 with @xmath15 , \label{eq : g}\ ] ] where @xmath16 is the surface disorder variance . by assumption , thus , the disordered charge distributions on the two slabs are statistically uncorrelated .
the effects of surface charge correlations ( or  patchiness " ) and the internal structure of multivalent ions , which are assumed to be point - like here , will be considered elsewhere @xcite .
in the present model , we have two types of mobile ions , namely , mono- and multivalent ions , that are expected to behave very differently in the presence of dielectric boundaries . on general grounds , and if the surfaces are assumed to be strictly neutral ( charge free ) , mobile ions interact only with their image charges .
this interaction scales with the second power of the ionic valency , being thus much larger in the case of multivalent ions with @xmath17 ( e.g. , trivalent and tetravalent ions ) . in this latter case ,
the ion - image interactions are dominated by the self - image interaction .
note that in most realistic examples , ions are dissolved in water , which is a highly polarizable medium with @xmath18 at room temperature ( @xmath19  k ) and , thus , we often have @xmath20 , giving image charges of the same sign and , hence , repulsive self - image interactions for each ion in the slit . such repulsive interactions lead to statistical correlations ( manifested as ion depletion ) between the mobile ions and the bounding surfaces that are expected to be weak for monovalent ions , but quite strong for multivalent ions .
this leads to a complicated problem in which different ionic species show distinctly different couplings to local electrostatic fields generated by the boundaries . while the monovalent ions can be therefore described appropriately by mean - field - type theories , such as the poisson - boltzmann or the debye - hckel theory @xcite , the multivalent ions require an altogether different description that should account for such large correlation effects .
the situation is , therefore , analogous to the one found in the case of asymmetric coulomb fluids confined between _ charged _ surfaces , where the dominant factor determining the behavior of multivalent counterions is again their electrostatic coupling to the local fields that are generated by the boundaries .
the difference is however that in the the case of a charged surface the coupling is determined by the so - called electrostatic coupling parameter @xcite , depending on the charge density of the surface boundaries , enabling the so - called _ strong - coupling approximation _ for multivalent counterions at charged surfaces ( see , e.g. , refs . @xcite and references therein ) , or more generally , the _ dressed multivalent - ion theory _
@xcite , which provides a very good approximation for the study of asymmetric coulomb fluids over a wide range of parameters as verified by explicit - ion and implicit - ion simulations .
the dressed multivalent - ion theory reproduces the mean - field debye - hckel theory and the standard strong - coupling theory for counterion - only systems @xcite as two limiting theories in the regime of large and small debye screening lengths and , therefore , bridges the gap between these two limits @xcite .
the key step in this latter approach is to integrate out the degrees of freedom associated with monovalent ions by means of a linearization approximation , justified only for highly asymmetric coulomb fluids @xmath17 , and yielding an effective debye - hckel ( dh ) interaction between the remaining multivalent ions and the surface charges ( if any ) @xcite .
then , one expands the partition function of the system in terms of the fugacity ( or bulk concentration , @xmath21 ) of multivalent ions or in terms of the inverse electrostatic coupling parameter for counterion - only systems @xcite .
the leading order of the virial expansion can be cast into a simple analytic theory because of its single - particle structure , which can successfully describe various features of strongly charged systems containing multivalent ions ( see , e.g. refs . @xcite and references therein ) .
the regime of applicability of this theory has been discussed extensively in recent literature @xcite , which we therefore do not reiterate here .    in the case of
strictly neutral ( charge - free ) surfaces , this strategy was shown to be effective as well @xcite , since , as is often the case in experimental systems containing asymmetric ionic mixtures @xcite , multivalent counterions of high valency ( for instance , cohex@xmath22 or polyamines such as spd@xmath22 , sp@xmath23 ) are present only in small bulk concentrations , of about a few mm , justifying fully the virial expansion scheme that underlies the dressed multivalent - ion theory .
this is the approach of choice that we adopt also in the present context , where the surface boundaries carry random charges and are neutral only on the average .
the only additional step here is that , due to the quenched disorder in the surface charge distribution , the free energy of the system follows by averaging over the whole ensemble of charge distributions , @xmath12 @xcite .      without delving further into details available in previous publications @xcite ,
we proceed by giving the general expression for the free energy of the two - slab system considered in this work , which has the form @xcite
@xmath24- \lambda_c k_{\mathrm{b}}t\int { \mathrm{d}}{{\bf r}}\ , \omega({{\bf r}})\ , \rme^{-\beta u({{\bf r}})}.\ ] ] here , @xmath25 and the green s function of the system is defined via @xmath26 where @xmath27 is the debye ( or salt ) screening parameter , which is non - zero only outside the region occupied by the dielectric slabs , i.e. , @xmath28 with @xmath29 being the total bulk concentration of monovalent ions .
the first term in eq . ( [ eq : exact_gpot ] ) represents the free energy of the system in the absence of multivalent ions .
it stems from _ self - interactions _ of disorder charges on the bounding surfaces of the slabs with their image charges ( note that the random charge distributions on the two slab are uncorrelated and , on the average , do not interact with each other ) .
this term depends only on the disorder variance , @xmath30 , and is non - vanishing only in inhomogeneous systems with a finite dielectric discontinuity at the bounding surfaces and/or a spatially inhomogeneous distribution of salt ions .
this contribution has been analyzed in the context of the fluctuation - induced forces between disordered surfaces in vacuum or in a weakly coupled coulomb fluid @xcite .    the second term in eq .
( [ eq : exact_gpot ] ) represents the contribution from multivalent ions on the leading ( virial ) order , in which @xmath31 is the _ effective _ single - particle interaction energy @xcite
@xmath32 ^ 2 .
\label{eq : u2}\ ] ] the first term in the above expression represents the self - image interaction of individual multivalent ions in the slit and the second term represents the contribution of the surface charge disorder to the effective single - particle interaction energy .
in the first term , @xmath33 gives the dielectric and/or salt polarization , or the image - charges effects ( corresponding to the generalized born energy ) , in which the free - space green s function , @xmath34 ( corresponding to the formation energy of individual ions in a homogeneous background ) , defined via @xmath35 , is subtracted from the total green s function .
the second term in eq .
( [ eq : u2 ] ) is found to be proportional to the disorder variance and shows an explicit temperature dependence and a quadratic dependence on the green s function and the multivalent - ion charge valency , @xmath36 ( these latter features can be understood by noting that the disorder term indeed stems from the sample - to - sample fluctuations , or variance , of the sample - dependent single - particle interaction energy as discussed in ref .
@xcite ) .    in the strong - coupling limit or within the multivalent dressed - ion theory @xcite , the number density of multivalent counterions can be obtained in terms of the effective single - particle interaction energy as @xcite @xmath37    in the specific example of two planar slabs , we can take advantage of the translational invariance of the green s function with respect to transverse ( in - plane ) direction coordinates @xmath38 and @xmath39 , since @xmath40 is only a function of @xmath41 , @xmath3 and @xmath42 , and write the free energy in terms of its fourier - bessel transform @xmath43 defined through @xmath44 for two semi - infinite slabs , one has @xmath45.\end{aligned}\ ] ] where @xmath46 , and @xmath47 in the absence of salt screening ( @xmath48 ) , @xmath49 reduces to the bare dielectric discontinuity parameter @xmath50 which gives a measure of the magnitude of  dielectric - image " charges ; while in a _ dielectrically homogeneous _ system , we have @xmath51 which gives a measure of the salt - induced image effects , or loosely speaking ,  salt - image " charges .
both image charge effects lead to depletion of ions from the vicinity of dielectric surfaces in a medium of relatively large polarizability , but they exhibit some fundamental differences , which have been discussed in the context of charged surfaces in our previous works @xcite .     ( a )
( b )    0.1 cm     ( c )    in the two - slab geometry , we can then re - express the ( number ) density profile of multivalent ions in the slit region , @xmath11 , as @xmath52 where @xmath53.\\\end{aligned}\ ] ] the _ interaction _ free energy of the system ( per @xmath54 and per unit area , @xmath1 ) can be written as @xmath55 where we have subtracted additive terms that are independent of surface separation @xmath4 and defined @xmath56
let us first focus on the case of a dielectrically homogeneous system with @xmath57 .
the density profiles of multivalent ions for this case are shown in fig .
[ f : neutral_dist]a , where we have rescaled the density profiles with their bulk value @xmath9 and the position in the slit @xmath58 with the bjerrum length @xmath59 . the inter - surface distance and the screening parameter
are rescaled in the same way . in the figure , we have fixed @xmath60 and @xmath61 , equivalent to mono- and multivalent salt concentrations of @xmath62  mm and @xmath63  mm when the bjerrum length is taken as @xmath64  nm ( appropriate for water at room temperature , i.e. , with @xmath19  k and @xmath18 ) and the multivalent ion valency is taken as @xmath65 ( note that throughout this paper we focus on the case of asymmetric coulomb fluids with tetravalent ions but the generalization of our results to other values of @xmath36 is straightforward ) .
the surface charge disorder variance is shown in the figure in terms of the dimensionless _ disorder coupling ( or strength ) parameter _
@xcite @xmath66 which is varied in the figure as @xmath67 and 2 ( corresponding to strictly neutral surfaces with @xmath68 , and randomly charged surfaces with @xmath69  nm@xmath70 and 0.04  nm@xmath70 , respectively ) .
as seen , for strictly neutral ( or disorder - free ) dielectrics , the density profile ( black solid curve ) shows a partial depletion of multivalent ions from the vicinity of the surface boundaries .
this is caused by the salt - image repulsion ( corresponding to the first term in eq .
( [ eq : u2 ] ) ) , which , as noted in the previous section , are caused by the discontinuity in the distribution of salt ions across the dielectric interfaces , giving a non - vanishing @xmath49 according to eq .
( [ eq : delta_s_homogen ] ) .
however , when the surfaces are randomly charged ( dashed curves ) , the situation turns out significantly different and the multivalent ions are attracted quite strongly towards the surfaces despite the  soft " salt - image repulsions ( moreover , one can note from the shown profiles that a larger amount of multivalent ions are pulled into the slit from the bulk solution ) .
in fact , the resulting disorder potential acting on individual multivalent ions ( corresponding to the second term in eq .
( [ eq : u2 ] ) ) exhibits a singular ( logarithmic ) behavior at the two surfaces and gives an algebraically diverging density as @xmath71 when @xmath72 .
this kind of phenomena has been discussed in detail in the context of disordered charged surfaces bearing a net charge density @xcite , where we show that the singular behavior of the density profile is closely connected with the _ anti - fragile _ behavior of multivalent counterions : adding a quenched disorder component to an otherwise uniform distribution of surface charges leads to a lowered entropy ( or thermal disorder ) for the counterions . in the present context , with net - neutral surfaces , the disorder - induced effects are qualitatively similar as they are regulated only by the variance of the charge distribution , @xmath30 , rather than the net charge density of the surface .
finite dielectric discontinuity at the bounding surfaces , @xmath73 , creates a much stronger depletion effect than in the case of dielectrically homogeneous system ( see fig .
[ f : neutral_dist]b ) . in fact
, the density profile of multivalent ions vanishes in the immediate vicinity of the surfaces because of the strong dielectric - image repulsion , provided again by the first term in eq .
( [ eq : u2 ] ) .
this term is singular itself and diverges on approach to the surface as @xmath74 when @xmath72 , thus dominating over the singular potential created by the charge randomness on the boundaries and , leading to the vanishing contact density of multivalent ions regardless of the disorder strength , @xmath75 , as seen in the figure ( even though a larger amount of multivalent ions are found in the slit at larger values of @xmath75 ) .
the competition between these two mechanisms of repulsion and attraction leads to a change in the shape of the density profile from unimodal to bimodal beyond a threshold value of the disorder strength , @xmath76 . the data in fig .
[ f : neutral_dist]c show @xmath76 as a function of the dielectric discontinuity parameter , @xmath77 , for fixed @xmath60 , @xmath63  mm and different values of @xmath78 and 60  mm ( corresponding to @xmath79 and 0.58 , respectively ) .
clearly , for larger salt screening and/or dielectric discontinuity parameter , a larger degree of charge disorder is required in order to counteract the image - charge repulsions .
these features of the density profiles are qualitatively similar to those found in the case of charged surfaces @xcite and , therefore , we shall not delve further into the details of the behavior of the density profile and proceed with the analysis of the effective interaction between net - neutral surfaces .
effective interactions between neutral dielectric slabs are standardly described in terms of the vdw interactions as , for instance , formulated within the lifshitz theory @xcite .
the vdw interaction pressure for two plane - parallel slabs can be expressed as @xcite @xmath80 where the first term comes from the thermal zero - frequency mode of the electromagnetic field - fluctuations and the second term comes from the higher - order matsubara frequencies , pertaining to quantum fluctuations .
@xmath81 is then the quantum part of the hamaker coefficient , which , based on an upper bound estimate in the case of hydrocarbon slabs interacting across an aqueous medium , can be taken as @xmath82 @xcite .
the vdw interaction is dominant mostly at the scale of a few nanometers for the inter - surface distance @xcite .    in the present model
, one needs to account for the electrostatic contribution to the inter - surface pressure , @xmath83 , as well .
this contribution can be written as the difference between the slit pressure of multivalent and monovalent ions and the bulk pressure @xmath84 , where @xmath29 is the total bulk concentration of monovalent ions , i.e. , @xmath85 .
the slit pressure can be calculated as @xmath86 , in which the first term is the contribution of multivalent ions that follows from the derivative of the free energy expression , eq .
( [ e : f_neutral ] ) , w.r.t .
the inter - surface distance when all other parameters are kept constant , while the second term is the contribution of monovalent ions that has been expressed in terms of the total mid - plane density of monovalent ions @xcite .
this latter quantity can be estimated here through the relation @xmath87\big|_{q=1}$ ] for @xmath88 , which has been shown by means of explicit monte - carlo simulations of a disorder - free system @xcite to give an accurate estimate of the distribution of monovalent ions in the slit .
the electrostatic pressure can thus be decomposed into its different components as @xmath89 where the components are obtained explicitly as @xmath90 -c_0\\ \label{eq : p_mon }   & & \beta p_{mon}=n(0)-n_b.\end{aligned}\ ] ] the contribution @xmath91 arises from self - interactions of random charges on the surfaces of the two slabs with their ( dielectric and/or salt ) image charges .
this contribution stems from the first term of the interaction free energy ( [ e : f_neutral ] ) and is present irrespective of multivalent ions , being typically comparable in range and strength with the vdw pressure .
it can , however , scale differently with the inter - surface distance and can be repulsive or attractive depending on the sign of dielectric discontinuity parameter , @xmath77 ( e.g. , it is repulsive for @xmath73 , applicable to the cases studied in this paper ) , resulting thus in a rather diverse behavior for the inter - surface interaction between two net - neutral slabs in the absence of multivalent ions @xcite .
the two terms @xmath92 and @xmath93 correspond to the disjoining or osmotic pressure contributions of multivalent and monovalent ions , respectively .
the former stems from the second term of the interaction free energy ( [ e : f_neutral ] ) , while , as noted above , the latter is included heuristically through eq .
( [ eq : p_mon ] ) .
, @xmath63  mm and @xmath94 and 35  mm as indicated on the graph ( colored dashed curves ) .
the corresponding vdw pressure for these parameter values closely overlap ( black dashed curves ) .
panel ( b ) shows the total pressure ( @xmath95 ) along with its electrostatic ( @xmath83 ) and vdw ( @xmath96 ) components for @xmath63  mm and @xmath62  mm .
the light - colored dashed curves show @xmath83 for @xmath97 and 35  mm ( from right to left ) . ,
title="fig:"](a )    0.3 cm    , @xmath63  mm and @xmath94 and 35  mm as indicated on the graph ( colored dashed curves ) .
the corresponding vdw pressure for these parameter values closely overlap ( black dashed curves ) .
panel ( b ) shows the total pressure ( @xmath95 ) along with its electrostatic ( @xmath83 ) and vdw ( @xmath96 ) components for @xmath63  mm and @xmath62  mm .
the light - colored dashed curves show @xmath83 for @xmath97 and 35  mm ( from right to left ) . ,
title="fig:"](b )    \(a )   + -0.4 cm   ( b )    \(c )   + -0.4 cm   ( d )    \(e )   + -0.4 cm   ( f )    the total pressure acting on the slabs is finally given by @xmath98 this quantity is plotted in fig .
[ f : neut_kappa]a as a function of the rescaled inter - surface distance for fixed @xmath99 , @xmath63  mm and four different values of @xmath94 and 35  mm , which correspond to @xmath100 and 0.45 , respectively ( colored dashed curves ) . for these same parameter values
, we also show the vdw pressure in fig .
[ f : neut_kappa]a ( black dashed curves ) , but these latter curves closely overlap and are not discernible at the implied resolution .
this clearly indicates that the dependence of the total pressure on the salt screening parameter enters mainly through the electrostatic contribution ( see below ) . in fact , as one can see from the figure , the total pressure , @xmath101 , is less attractive than the standard vdw pressure at large separations , but , as @xmath7 decreases , the total pressure becomes more strongly attractive as compared with the vdw contribution .
we show the electrostatic and vdw components of the total pressure separately in fig .
[ f : neut_kappa]b for @xmath63  mm and @xmath62  mm ( @xmath102 ) .
the electrostatic contribution , @xmath83 , shows a non - monotonic behavior : it is repulsive and decays to zero at large separations , while it becomes strongly attractive and diverges at small separations , showing thus a weakly repulsive maximum at intermediate inter - surface separations .
the light - colored dashed curves show the same quantity for larger values of @xmath97 and 35  mm ( from right to left ) . clearly , the extent of the repulsion and the non - monotonic behavior of the electrostatic pressure become more pronounced as @xmath7 is increased .
non - monotonic interaction pressures due to image - induced ion depletion have also been found in the absence of surface charge disorder but at relatively large bulk salt concentrations ( e.g. , above 250 mm ) @xcite .    before proceeding with a more detailed analysis of the electrostatic contribution
, we emphasize here that @xmath83 and @xmath96 show clearly different qualitative behaviors as one can see from fig .
[ f : neut_kappa]b .
the standard dlvo description of colloidal interactions @xcite in terms of the vdw interaction of neutral surfaces is therefore insufficient when the surfaces are neutral only _ on the average _ and otherwise carry random positive and negative charges , and/or when the system contains also an asymmetric coulomb fluid .
our results show that only a small degree of charge randomness with , e.g. , a surface charge variance of @xmath103  nm@xmath70 ( corresponding to @xmath99 ) , will be enough to generate a sizable deviation from the standard vdw prediction .      in order to gain further insight into the intriguing role of disorder - induced effects in the presence of multivalent ions
, we examine the behavior of the effective electrostatic pressure , @xmath83 , and its components in more detail .
first , we consider two strictly neutral ( disorder - free ) surfaces , in which case the disorder effects vanish and the electrostatic inter - surface pressure , @xmath83 , is given only by the osmotic pressures of multivalent and monovalent ions .
the results are shown by black solid curves in figs .
[ f : neutral_pressure](a , b ) and figs .
[ f : neutral_pressure](c , d ) for a dielectrically homogeneous ( @xmath57 ) and a dielectrically inhomogeneous system ( @xmath104 ) , respectively . for the latter case ,
the contributing components of the electrostatic pressure , i.e. , @xmath105 and @xmath106 , are plotted in figs .
[ f : neutral_pressure](e , f ) for comparison . the vdw pressure is shown by black dotted curves in panels a to d. in all cases , we have taken @xmath63  mm and @xmath62  mm .
it turns out that the electrostatic pressure in a disorder - free system ( @xmath107 ) , which is given by @xmath108 , is negative ( attractive ) for all inter - surface separations and tends to zero at large separations because of salt screening effects .
this means that the negative bulk pressure is stronger in magnitude than the slit pressure for both cases of @xmath57 and 0.95 .
this is because of the image - induced depletion of ions generated by salt and/or dielectric - image charges that are of the same sign when the surfaces are immersed in a medium of larger dielectric constant as assumed here . in the limit @xmath109 , @xmath83 in an inhomogeneous system
reduces to the bulk pressure , @xmath84 , due to complete depletion of ions from the slit under the strong repulsive force of dielectric - image charges , while in a dielectrically homogeneous system , complete depletion is not achieved ; hence , @xmath83 is less attractive in this latter case . on the other hand ,
for the given parameter values , the vdw pressure is more ( less ) attractive than the electrostatic pressure , @xmath83 , for @xmath104 ( @xmath57 ) as one can see from the figures .
next , we discuss how the introduction of charge disorder on the inner surfaces of the net - neutral dielectrics affects the effective electrostatic interaction between them . to this end , we vary the disorder strength parameter from @xmath107 up to @xmath110 ( corresponding to @xmath111  nm@xmath70 ) and divide the results into two categories of `` weakly '' and `` strongly '' disordered cases .    in the _ weak disorder regime _ ( typically @xmath112 ) , the electrostatic interaction pressure , @xmath83 , becomes gradullay more repulsive or positive ( especially at small separations ) as compared to that of strictly neutral surfaces when @xmath75 is increased up @xmath113 ( see colored dashed curves in figs .  [
f : neutral_pressure]a and c for @xmath57 and @xmath104 , respectively ) . the repulsive interaction pressure is stronger for larger values of the dielectric discontinuity parameter , @xmath77 ( compare panels a and c ) .
the most remarkable feature of our results is that @xmath83 in this regime adopts a non - monotonic behavior : for very small @xmath75 , it develops a point of zero electrostatic pressure followed by a shallow attractive minimum at intermediate separations ( see , e.g. , @xmath114 , red dashed curve , in panel a , or @xmath114 and 0.5 , red and blue dashed curves in panel c ) .
the depth of this minimum decreases and the interaction profile instead develops a repulsive ( positive ) hump at large values of @xmath75 ( e.g. , @xmath113 ) .
these behaviors are in clear contrast to what is observed in the case of strictly neutral ( disorder - free ) surfaces ( black solid curves ) .    in order to elucidate the origin of this difference
, we compare different components of @xmath83 in the case of an inhomogeneous system with @xmath104 and @xmath115 in fig .
[ f : neutral_pressure]e . for comparison , the light - blue solid curve and the black solid curve ( which is seen more clearly in the inset )
show the osmotic contribution of monovalent and multivalent ions , @xmath106 and @xmath105 , in the case of strictly neutral surfaces ( @xmath107 ) .
both contributions are attractive ( negative ) , while the former is the dominant one as expected since we have assumed that multivalent ions enter in small bulk concentrations ( here , @xmath63
mm ) . by introducing a weak charge randomness ( @xmath115 for the relevant curves in fig .
[ f : neutral_pressure]e ) , the negative pressure of both multivalent ( @xmath92 , red dashed curves , inset ) and monovalent ions ( @xmath93 , blue dashed curve , main set ) decrease in magnitude and , thus , become less attractive as compared with the disorder - free case .
it is important to note that the distribution of both monovalent and multivalent ions and , thus , their respective osmotic pressure , @xmath93 and @xmath92 , are affected by the surface charge disorder through the disorder - induced , single - ion potential , which is the second term in eqs .
( [ eq : u2 ] ) or ( [ eq : resu_neutral ] ) .
the decrease in the magnitude of attractive pressure components @xmath93 and @xmath92 in the weak disorder regime is because more mono- and multivalent ions are attracted into the slit from the bulk solution in the presence of surface charge disorder ( see fig .
[ f : neutral_dist ] and fig .
[ f : sum_multi ] below ) .
nevertheless , at very small separations , all ions are again totally depleted due to the repulsive forces of image charges and @xmath93 and @xmath92 reduce to the same bulk values as @xmath106 and @xmath105 do when @xmath116 , i.e. , @xmath117 and @xmath118 , respectively . the crucial role , however , is played by the repulsive disorder self - interaction component , @xmath91 ( green dashed curve ) , whose effect is amplified by increasing the charge disorder strength and/or the interfacial dielectric mismatch as follows from eq . .
it is the interplay between this repulsive disorder - induced pressure ( which dominates at small to intermediate separations ) and the attractive osmotic pressure of ions ( which dominate at larger separations ) that gives rise to the non - monotonic behavior mentioned above . for sufficiently small @xmath75 ,
the pressure of multivalent ions is typically too small ( as compared to the other two components ) to change this qualitative behavior .
this picture changes drastically in the _ strong disorder regime _
( typically @xmath119 ) , where the pressure due to multivalent ions becomes a key factor . in fig .
[ f : neutral_pressure]b and d , we plot @xmath83 for larger disorder strength parameters , @xmath120 up to 5 for @xmath57 and @xmath104 , respectively . in this regime
, we find a reverse trend caused by the surface charge disorder : the repulsive hump at intermediate separations now diminishes and eventually disappears when @xmath75 is increased to larger values ; one thus finds a highly attractive ( negative ) interaction pressure with a range that can be much larger than that of the vdw interaction pressure ( black dotted curves ) , or the mere image - induced , ion - depletion pressure ( black solid curves with @xmath107 ) , at sufficiently large disorder strengths .
the different components contributing to the electrostatic pressure for this case are shown in fig .
[ f : neutral_pressure]f for @xmath99 with other parameters being the same as figs .
[ f : neutral_pressure]c and d. as seen , while the disorder self - interaction contribution , @xmath91 , has become only slightly more repulsive ( green dashed curve ) , the other two components , @xmath92 and @xmath93 , show significant changes as compared to their counterparts in the weak disorder regime in fig .
[ f : neutral_pressure]e .
specifically , the attractive contribution due to the osmotic pressure of multivalent ions ( red dashed curve ) becomes much stronger than what we find in the weakly disordered or disorder - free cases , while the contribution from monovalent ions now becomes repulsive ( positive ) in contrast to its attractive ( negative ) behavior in these latter cases ( compare panels e and f ) .
the monovalent contribution becomes positive because of an increase in the number of monovalent ions in the slit that creates a larger osmotic pressure .
the average number of multivalent ions also increases in the slit as @xmath75 increases ; however , these ions correlate strongly with the disorder charges on the bounding surfaces of the slabs ( through the second term in eq .
( [ eq : resu_neutral ] ) which contributes directly to @xmath92 , eq . ( [ eq : p_c ] ) ) and , as a result , give rise to an effective attraction between the net - neutral slabs through the pressure component @xmath92 .
this effect is in clear contrast with the image - induced , ion - depletion mechanism , in which a _ decrease _ in the number of mono- and multivalent ions in the slit gives rise to attractive osmotic pressures on the slabs similarly from both types of ions .
this behavior originates from the ( singular ) attractive , single - ion potential , which is created by surface charge disorder .
it is thus one of the most fundamental aspects of the coupling between surface charge disorder and mobile multivalent ions that follows from our results .    , per rescaled surface area , @xmath121
, is plotted versus the rescaled inter - surface separation , @xmath122 , for a system similar to those considered in figs .
[ f : neutral_pressure]c and d. ]    , @xmath99 , @xmath123  mm and @xmath63 , 5 and 10  mm as indicated on the graph ( from top to bottom ) .
the corresponding electrostatic and vdw components of the total pressure , i.e. , @xmath83 and @xmath96 , are shown in rescaled units as well ( @xmath83 is shown by the light - colored dashed curves and @xmath96 by the black one , which closely overlap ) . panel ( b ) shows @xmath95 along with @xmath96 , @xmath83 and the three components contributing to the latter , i.e. , @xmath91 , @xmath92 and @xmath93 for fixed @xmath124  mm and other parameter values as in ( a ) . , title="fig:"](a )    0.3 cm    , @xmath99 , @xmath123  mm and @xmath63 , 5 and 10  mm as indicated on the graph ( from top to bottom ) .
the corresponding electrostatic and vdw components of the total pressure , i.e. , @xmath83 and @xmath96 , are shown in rescaled units as well ( @xmath83 is shown by the light - colored dashed curves and @xmath96 by the black one , which closely overlap ) .
panel ( b ) shows @xmath95 along with @xmath96 , @xmath83 and the three components contributing to the latter , i.e. , @xmath91 , @xmath92 and @xmath93 for fixed @xmath124  mm and other parameter values as in ( a ) . ,
title="fig:"](b )    the remarks in the preceding discussions on the attraction ( depletion ) of multivalent ions to ( from ) the slit region can be corroborated by considering the average total number of ions in the slit .
for multivalent ions , this quantity can be calculated from @xmath125 .
it shows a very different behavior as a function of the inter - surface separation in the weak disorder and strong disorder regimes as one can see in fig .
[ f : sum_multi ] ( we have used the same parameters as in figs .  [
f : neutral_pressure]c and d ) .
although , at fixed inter - surface distance , always more multivalent ions are attracted to the slit by increasing the disorder strength , the response of multivalent ions to the decrease in slab separation is quite different for different disorder regimes : for disorder - free and weakly disordered surfaces , multivalent counterions are quickly depleted from the slit by decreasing the slab separations , while , for sufficiently large disorder strengths , these ions ( as well as monovalent ions that are not shown here ) are more strongly attracted to the slit from the bulk solution due to stronger attraction experienced from the randomly charged , inner surfaces of the slabs . in the intermediate regime of disorder strengths ,
one find a non - monotonic behavior upon deceasing the inter - surface separation , with the average total number of multivalent ions first decreasing and then increasing at very small separations .
one should , however , note that in dielectrically inhomogeneous systems such as the one considered in fig .
[ f : sum_multi ] with @xmath104 , dielectric - image repulsions eventually win and , thus , even the curves for highly disordered systems eventually turn downward and we find @xmath126 for @xmath116 or , equivalently , @xmath127 for all cases in figs .  [ f : neutral_pressure]c and d ( these limiting behaviors occur at very small separations that are not physically meaningful , e.g. , below @xmath128   , and , for the sake of presentation , they have not been shown in the plots ) . the crossover from weak to strong disorder regimes
can be effected also by decreasing the salt screening parameter through a decrease in the bulk concentration of monovalent ions , @xmath7 ( and/or increasing the bulk concentration of multivalent ions , @xmath9 ) .
for instance , at sufficiently large @xmath7 such as @xmath123  mm , the electrostatic interaction pressure , @xmath83 , at the disorder strength value of @xmath99 and multivalent salt concentration of @xmath63  mm shows a weak - disorder behavior ( see fig .
[ f : kappa_c0]a ) , in contrast to what we found for @xmath99 at lower salt concentration of @xmath62  mm in fig .
[ f : neutral_pressure]d .
even though increasing @xmath9 up to @xmath124  mm enhances the attractive pressure as shown in the figure , the behavior of the effective interaction still remains in the weak - disorder regime as can be verified from the pressure components in fig .
[ f : kappa_c0]b ( e.g. , compare @xmath93 and @xmath92 with those in fig .
[ f : neutral_pressure]e ) .
this is clearly because of strong salt screening effects at large @xmath7 , in which case the strong - disorder behavior can be achieved only by taking even larger values of @xmath75 .
in this work , we have studied the effective interactions between _ net - neutral _ dielectric slabs that carry _ quenched _ ( fixed ) random charge distributions on their apposing surfaces while they are immersed in an asymmetric aqueous electrolyte ( coulomb fluid ) containing mobile monovalent and multivalent ions .
the effective interaction between quenched , random charge distributions have been considered in a series of recent works @xcite , and , while the role of debye screening due to a weakly coupled ( monovalent ) salt solution on these interactions have been analyzed , the role of multivalent ions that couple strongly with surface charges has received much less attention @xcite .    the interaction between ( electro- ) neutral dielectric slabs is traditionally described in terms of vdw or casimir - type forces @xcite that are always attractive between identical bodies in vacuum or in a polarizable medium such as water or an aqueous salt solution . in the dlvo context @xcite , the vdw attraction is counteracted by the mean - field electrostatic repulsion between like - charged surfaces , providing thus a mechanism for the stability of colloidal dispersions . in the case of strictly neutral surfaces , however , the electrostatic inter - surface repulsion is obviously absent but , in addition to the vdw interaction , there are image - induced , ion - depletion forces that are attractive as well and arise because the mobile solution ions are depleted from the vicinity of the dielectric interfaces as a consequence of the repulsion from their image charges ( these image charges are of the same sign when the surfaces are immersed in a medium of larger dielectric constant as compared with that of the slabs ) .
this effect is intrinsically non - mean - field and occurs because of the discrete nature of ions , neglected in the collective mean - field description based on the standard poisson - boltzmann theory @xcite .
image - induced , ion - depletion effects for strictly neutral dielectric surfaces have been studied extensively in recent years ( see , e.g. , refs . @xcite and references therein ) and
the role of multivalent ions , in particular , has been considered in refs .
@xcite .    in the present work ,
we add a new feature to the current understanding of the effective interaction between neutral surfaces and consider the situation , in which the surfaces are neutral only _ on the average _ but otherwise carry a quenched random distribution of positive and negative charges .
indeed , heterogeneously or randomly charged surfaces are commonplace in soft matter and have attracted a lot of attention in recent years ( see , e.g. , refs .
we have thus shown that the presence of _ quenched surface charge disorder _ in conjunction with _
mobile multivalent ions _ in the ionic solution leads to remarkable and significant qualitative changes in the standard picture commonly accepted for the interaction between neutral bodies .
the charge disorder has several manifestations in the present context : first , the self - interactions between the quenched random charges on the inner surfaces of the dielectric slabs and their image charges lead to a repulsive , short - ranged interaction pressure , @xmath91 , that tends to counteract the attractive vdw interaction pressure , @xmath96 as discussed thoroughly elsewhere @xcite . the surface charge disorder , on the other hand , modifies the single - ion interaction potentials in a way that creates a singular attractive potential acting between ions and the dielectric boundaries @xcite .
this effect works against the ion - image repulsions that tend to deplete ions from the slit region between the slabs .
thus , as a result of this attractive disorder - induced potential , more mono- and multivalent ions are accumulated in the slit when the confining dielectric boundaries are randomly charged .
however , we find qualitative differences in the ways this excess attraction affects the osmotic pressure of ions .
the osmotic pressure due to monovalent ions , @xmath93 , can become _ repulsive _ which , in fact , agrees with the standard mean - field picture since monovalent ions are only weakly coupled to the bounding surfaces and create a larger entropic pressure upon further accumulation in the slit . on the other hand ,
the osmotic pressure due to multivalent ions , @xmath92 , can become even more _ attractive _ !
this is in clear contrast with the image - induced , ion - depletion mechanism , in which a _ decrease _ in the number of mono- and multivalent ions in the slit gives rise to ( stronger ) attractive osmotic pressures on the slabs , similarly for both types of ions .
this rich behavior stems from combined effects of surface charge disorder and the presence of mobile multivalent ions , with the latter creating strong inter - surface attractions upon further accumulation in the slit . as a result , the effective total interaction pressure between randomly charged , net - neutral dielectric surfaces can differ qualitatively from what one expects based on the standard vdw and/or image - induced , ion - depletion interactions for strictly neutral surfaces , where they predict purely attractive interactions ( unless at relatively large bulk salt concentrations , e.g. , above 250 mm , as noted in ref .
@xcite ) . in particular , the net effect from the competing electrostatic components due to disorder , mono- and multivalent ions results in a distinct , disorder - induced _ non - monotonic _ behavior .    in the _ weak disorder regime _ , the net electrostatic pressure becomes gradually more repulsive ( positive ) at intermediate to large separations as the disorder coupling parameter ( or equivalently the disorder variance , @xmath30 ) on the bounding surfaces of the slabs is increased from zero ( strictly neutral or disorder - free system ) , exhibiting first a shallow minimum and then a pronounced repulsive hump at intermediate separations . as the disorder coupling parameter is further increased into the _ strong disorder regime _
( typically beyond @xmath129 ) , the behavior of the net electrostatic pressure is reversed ; the repulsive hump is diminished and one finds a strongly attractive ( negative ) interaction pressure with a range and magnitude larger than that of the vdw or the image - induced , ion - depletion interaction pressures between strictly neutral surfaces .
the crossover from weak to strong disorder regimes can be effected also by decreasing the salt screening parameter ( through decreasing the bulk concentration of monovalent ions ) and/or increasing the bulk concentration of multivalent ions .
the two different regimes of the disorder effects , the strong and the weak disorder regimes , bear also some resemblance to the strong - weak coupling dichotomy that exists for the net - charged surfaces , dependent on the electrostatic coupling parameter @xcite .
it is no small feat to be able to partition the behavior of this complicated system by two dimensionless coupling parameters that effectively govern its salient characteristics .
finally , we note the qualitative differences between the ion - mediated interactions we find between net - neutral surfaces in this work and those we reported between randomly charged surfaces carrying a finite mean surface charge density , @xmath130 , in refs .
one of the key differences is that the repulsion between mean surface charges on the inner surfaces of the slabs , @xmath131 , in the latter case contributes a dominant repulsive pressure to the net electrostatic pressure , which is stronger , both in magnitude and range , than @xmath91 and @xmath96 . on the other hand ,
when the surfaces carry a finite mean surface charge density , the osmotic pressure of monovalent ions , @xmath93 , turns out to be always repulsive and comparable in magnitude and range to @xmath131 ; these repulsive pressure components completely mask the short - range interaction pressures , @xmath91 and @xmath96 , whose roles in competition with the attractive pressure due to multivalent ions , @xmath92 , are brought up only in the case of net - neutral surfaces .
notably , we find that @xmath93 in the case of surfaces carrying a finite mean surface charge density hardly responds to changes in the disorder strength , while , in the case of net - neutral surfaces , it can change from a purely attractive component ( due to image - induced , ion - depletion effects ) for disorder - free and weakly disordered surfaces to a repulsive one ( due to disorder - driven ion accumulation in the slit ) for strongly disordered surfaces .
another difference between the two cases is that multivalent ( counter- ) ions in the case of surfaces carrying a finite mean surface charge density correlate strongly also with the mean charge densities of the confining boundaries of the slit , in line with the strong - coupling paradigm for electrostatics of charged surfaces known to generate strong like - charged surface attractions ( see , e.g. , refs .
such strong - coupling effects are clearly absent in the present case with net - neutral surfaces .    in summary ,
our analysis provides new insight into the intricate role of surface charge disorder in the context of net - neutral surfaces and its fingerprints on the effective ion - mediated interactions between them .
these interactions are predicted to occur with a range and magnitude comparable to or , for strongly disordered systems , much larger than the vdw and/or image - induced , ion - depletion effects considered previously only in the case of strictly neutral surfaces .
our predictions will thus be amenable to numerical and/or experimental verification .
a.n . acknowledges partial support from the royal society , the royal academy of engineering , and the british academy ( uk ) .
acknowledges support from the school of physics , institute for research in fundamental sciences ( ipm ) , where she stayed as a short - term visiting researcher during the completion of this work .
r.p . acknowledges support from the slovene agency of research and development ( arrs ) through grant nos .
p1 - 0055 and n10019 as well as the travel grant from the simons foundation and the aspen center for physics , supported by national science foundation grant phy-1066293 .
this constraint on multivalent counterions can be achieved by enclosing the outer regions ` behind ' the slabs in semi - permeable membranes , although , in the regime of parameters considered in this paper ( e.g. , with slab thicknesses much larger than the debye screening length ) , this constraint has no sizable impact on our results .