Arxiv:1910.00465V2 [Cond-Mat.Soft] 29 Dec 2019 Xrsin 19.Rcnl,Hwvr Twssonthat Shown Was It However, Recently, Simple No [1–9]

arXiv:1910.00465v2 [cond-mat.soft] 29 Dec 2019 xrsin 19.Rcnl,hwvr twssonthat shown was it however, Recently, simple no [1–9]. are there expressions which for ob- potential to interaction PB-equation has an the tain one of solutions that numerical to was resort solutions to salt in double-layers lapping intact. remain de- will exponential repulsions condition double-layer this of under cay rea- that no of assume is to screening there within son and length, vanishes, other even Debye or each diminishes, the charges approach than smaller im- solutions distances surfaces salt double- charged in of when feature mersed However, generic a repulsions. be layer to exponential assumed Usually is screening [1–9]. potential (DLVO) Overbeek screened Derjaguin–Landau–Verwey– Together classical exponentially the form attractions, length. repulsions Waals screening der Debye van the with by decay set a with distance repulsions, double-layer expo- decaying to nentially leads eventually Poisson– that small, the equation, of (PB) is Boltzmann solution surfaces approximate poten- two an electrical the for allows between the which mid-plane that the is in implication tial The over- [1–9]. weakly lap only as- double-layers the electrical under evaluated that usually sumption is solution electrolyte in 2 ∗ † aoaoyo hsclCeity aut fCeia Engi Chemical of Faculty Chemistry, Physical of Laboratory E-mail: E-mail: o ogtm,tegnrloiinaotsrnl over- strongly about opinion general the time, long For immersed surfaces charged two between repulsion The [email protected] [email protected] ijn Stojimirović,Biljana itne nte leri eiei eoee.I hsre this In distance. recovered. separation di is square regime inverse algebraic an another assumes distances profile in pressure clear arise disjoining less the is interactions evidence algebraic the albeit that regimes, counterion-only show evidenc We experimental present interactions. we Here distance. separation codn ocnetoa idmeeti obelyrfo double-layer electric wisdom conventional to According idoe nvriyo ehooy OBx53 idoe 5 Eindhoven 513, Box PO Technology, of University Eindhoven .INTRODUCTION I. xeietlEiec o leri obeLyrForces Double-Layer Algebraic for Evidence Experimental 1 eateto nrai n nltclCeity Univers Chemistry, Analytical and Inorganic of Department cecsI,3 uiEns-nemt 25Gnv,Switzer Geneva, 1205 Ernest-Ansermet, Quai 30 II, Sciences ey nttt o aoaeil cec,UrctUniver Utrecht Science, Nanomaterials for Institute Debye 3 1 a tHffLbrtr o hscladClodChemistry, Colloid and Physical for Laboratory Hoff ’t Van akVis, Mark aula ,Urct38 H h Netherlands The CH, 3584 Utrecht 8, Padualaan 2, ∗ ec Tuinier, Remco Dtd aur ,2020) 1, January (Dated: ern n hmsr nttt o ope oeua Sys Molecular Complex for Institute & Chemistry and neering loosre hnivsiaiginincreainef- correlation ion-ion investigating when observed also h ltstepesr rfieas a nagbacde- algebraic on an charge has high 1 also of of profile limit pendence pressure the the in Langmuir plates by that the found calculated He first was [19]. bath counterion-only immersed a plates two in The between interaction exponentially. the for decays solution salt profile with counterion system decays the the plate where to charged contrast single in [16–18], the algebraically to concentra- adjacent counterion the profile the if systems tion such namely con- In only counterions, fluid cations. only adjacent tains the sur- of charged charged negatively bath is where a surface systems, in are are These faces systems. only 15]. the [14, modifies pressure additionally decaying situ- algebraically these effects In overspill charge finite surfaces. ations charged of non-uniformly surfaces and for size studied in- further an were overlapping as Strongly re- double-layers distance. algebraic decays separation an pressure of square is disjoining verse result calculated the The be where be- can pulsion 13]. repulsion solution 12, salt the [10, a limit analytically equilibrium in Donnan surfaces Donnan zero-field two the the tween as po- In known electrical also constant [11]. case spatially a but around high tential, roam possibly ions a poten- situation varying pendant in spatially DH the but in the weak whereas In a tial, [10]. in proximity elec- diffuse weak close ions a in case, of surfaces form for the in field limit tric pendant a has apart, far surfaces for Debye–Hückel limit weak-potential (DH) the ,3 2, leri eednisaeas on ncounterion- in found also are dependencies Algebraic o h omrrgm.I oho hs cases these of both In regime. former the for letP Philipse, P. Albert ietepesr easa h nes of inverse the as decays pressure the gime tnedpnec.A ml separation small At dependence. stance fagbaclydcyn double-layer decaying algebraically of e csnral ea xoetal with exponentially decay normally rces ohsrnl vrapn swl as well as overlapping strongly both /h 0 B h Netherlands The MB, 600 2 iia oe-a eedniswere dependencies power-law Similar . t fGeneva, of ity 3 sity, n rgrTrefalt Gregor and land 1, † tems, 2 fects [20, 21]. Cantilever (a) 100 pH 10 (b) While reports of experimentally measured exponential Silica Particle 500 mM KCl 50 double-layer forces are abundant [22–28], there have been 0 only a few articles showing algebraic double-layer forces -50 [29–31] in literature. In all the presented cases these -100 forces were measured in counterion-only regimes. On the vdW Fit -150 Silica Particle other hand algebraic interactions for the strongly over- Radius (µN/m) Force/Effective 0 10 20 30 lapping double-layers have only been studied theoreti- Separation Distance (nm) cally [10, 12, 13], and experimental instances of algebraic forces between surfaces in 1:1 salt solutions, as far as we FIG. 1. (a) Schematic representation of colloidal probe exper- know, have not been reported yet. The main aim of this iment. Van der Waals force measured in 500 mM KCl. The · −21 paper is to present direct and quantitative evidence for extracted Hamaker constant is H = 2.6 10 J [36]. the existence of algebraic interactions between surfaces immersed in simple salt solutions. We also explore different regimes where the algebraic forces are important 20 minutes. The substrate with particles was mounted experimentally. on the fluid cell. The geometry of the experiment is shown in Fig. 1a. The deflection of the cantilever was recorded when the particle on the cantilever was cen- II. MATERIALS AND METHODS tered above the selected one on the substrate with the precision of about 100 nm. For one pair of particles, A. Force Measurements the deflection is recorded in 150 approach-retract cycles with the cantilever velocity of 400 nm/s. The measure- Forces were measured between spherical silica parti- ment was done on 3–5 different pairs for each solution cles (Bangs Laboratories Inc, USA) with an average re- concentration. The spring constant of the cantilever was ported size of 5.2 µm. The colloidal probe technique determined by the Sader method [37], and the deflection based on atomic force microscopy was used [32–35]. A was converted to force using Hookes law. The approach single silica particle was glued on tip-less cantilever (Mi- part of the recorded curves is averaged and down-sampled croMasch, Tallin, Estonia) by first immersing the can- for increasing the force resolution and readability of the tilever in a small drop of glue (Araldite 2000+). The figures. These curves are then also used for theoretical substrate was made separately by spreading silica parti- analysis. cles on a quartz microscope slide (Plano GmbH, Wetzlar, Germany), which was previously cleaned in piranha so- Our analysis of the experimental forces only focuses lution (3:1 mixture of H2SO4 (98%) and H2O2 (30%)). on the electrostatic double-layer contribution. However Cantilevers with glued-on particles and substrate were in some cases the van der Waals forces become non- both heated at 1200◦C for 2 h to achieve firm attachment negligible. Therefore the van der Waals force for the and removal of the glue. During this sintering process, present silica–silica system was measured in 500 mM KCl particles shrink about 15%, so the average diameter is and the respective Hamaker constant was extracted, see 4.4 µm [36]. A root mean square (RMS) roughness of Fig. 1b. In order to be sure that we are measuring only 0.63 nm was measured by AFM imaging in liquid. Solu- the van der Waals force and the electrostatic interactions were made using KCl (Sigma Aldrich) and Milli-Q tions are completely screened at 500 mM, we have also water (Millipore). Solutions of pH 10 were made with preformed measurements in 1 M KCl. The 500 mM and addition of 1M KOH (Acros Organics), and solutions of 1 M curves overlap, and therefore we are confident that pH 3 with 1M HCl (Fisher Scientific). For experiments the electrostatic forces are completely screened. The van at pH 5.6 no adjustment was done. der Waals contribution was then subtracted from all the Force measurements were done at room temperature experimental force profiles. The measured forces F (h) 23 2◦C with a closed loop AFM (MFP-3D, Asylum between two spherical particles were then converted to ± Research) mounted on an inverted optical microscope the equivalent disjoining pressure between two plates, (Olympus IX70). Both cantilever and the substrate were Π(h), by calculating the derivative of the force versus cleaned in ethanol and water, and plasma treated for the separation distance, h, and by applying the Derjaguin 3 approximation [6] where h is the distance between the plates. Further details of the implementation of the full PB 1 dF (h) Π(h)= , (1) theory including constant charge regulation model are −πR · dh given in Ref. 36. where R is the particle radius.

2. Strongly Overlapping Double-Layers B. Theory For strongly overlapping double-layers analytical ap- 1. Poisson–Boltzmann Theory proximations of the PB equation can be found as was demonstrated recently [10, 12, 13]. These approximations The disjoining pressure between two charged plates is are applicable for symmetric z : z electrolytes and here calculated by solving the Poisson–Boltzmann (PB) equa- we will focus only on the 1:1 case. For 1:1 electrolytes, tion the disjoining pressure can be calculated by simplifying

2 2 Eq. (4) d Ψ(x) βe0 −ziΨ(x) i 2 = c e , (2) dx − εε0 i Π=2kBTcs[cosh(ΨM) 1]. (5) X − For strongly overlapping double-layers, the electrical po- where Ψ = βe0ψ is the rescaled dimensionless electrostatic potential, x is the coordinate normal to the plates, tential between the plates deviates only little from its i.e. ci is the bulk concentration of ion i, zi is the valence average value, the electric ﬁeld in the gap is approximately zero. The mid-plane potential in the equation of ion i, e0 is the elementary charge, ε0 is the dielectric permittivity of vacuum, ε is the relative dielectric permit- above can be then be replaced by a constant ΨM = Ψ, where Ψ turns out to be the dimensionless Donnan po- tivity of water and β = 1/(kBT ) is the inverse thermal tential. The value of the Donnan potential can be deter- energy, where kB is the Boltzmann constant and T is the temperature. Throughout, T = 298 K and ε = 80 are mined from the electro-neutrality condition for the two used as appropriate for dilute aqueous solutions. The plates and the intermediate salt solution PB equation is solved numerically with diﬀerent bound- c c− λ + − = , (6) ary conditions. One typically assumes a constant charge 2cs h

(CC) or constant potential (CP) on the plates at differ- where c+ and c− refer to the cation and anion number ent distances h, however the charge on the plates can be concentrations between the plates, respectively. We have also regulated upon approach [38–42]. Here we employ introduced characteristic length λ defined as constant regulation approximation, where we introduce σ 4 the regulation parameter p [39, 43]. The regulation pa- λ = = 2 , (7) cse0 κ ℓGC rameter is defined as where σ is the surface charge density of the plates and C p = dl , (3) κ = √8πℓBcs is the inverse Debye length, with ℓB = Cdl + Cin 2 βe0/(4πεε0) the Bjerrum length, and where C and C are double-layer and inner-layer capac- e0 dl in ℓGC = (8) itance, respectively. The regulation parameter interpo- 2πℓBσ lates between CP (p = 0) and CC (p = 1) boundary con- is the Gouy–Chapman length. At separation distances ditions and can be used for describing charge-regulation below λ the contribution of counterions to pressure is in general way. dominant, while at h > λ background salt contribution The solution of the PB equation yields the potential becomes more important. The Gouy–Chapman length profile ψ(x) between two charged plates positioned at is connected to the thickness of the double-layer in the x = h/2 and x = +h/2. The potential at the mid- counerion-only case and it is an analog to Debye length − plane, ψM = ψ(0), then permits to calculate the disjoin- in systems containing salt. ing pressure By combining Eq. (6) with the Boltzmann equilibrium for exchanging ions between bulk and the gap −ziΨM(h) Π(h)= k T ci e 1 , (4) B ∓Ψ i − c± = cse , (9) X 4

ΨM−ΨD we can calculate the value of the Donnan potential where γ = e ,ΨD being the diffuse-layer potential at the isolated plate surface. The disjoining pressure in λ Ψ= sinh−1 . (10) this case is equal to the pressure between the plates, since − h the bulk pressure vanishes without salt and is calculated The zero-field pressure can be finally calculated by in- by using Eq. (15) serting Eq. (10) into Eq. (5) 2εε Π= k Tc eΨM = 0 γ2. (16) B − β2e2ℓ2 λ 2 0 GC Π=2k Tcs 1+ 1 . (11) B  h −  2 s The equation above shows that γ is actually the rescaled   pressure. For the constant charge boundary conditions For λ/h 1, i.e. large separation and low surface charge, the relation between the separation distance and pressure ≪ the above equation reduces to becomes [30]

2 λ 2ℓGC −1 Π= k Tcs . (12) h = arctan γ . (17) B h γ The pressure for strongly overlapping double-layers At large separation distances the pressure tends to zero π therefore decays algebraically Π 1/h2. Another al- and one can approximate arctan(x) . With this we ∝ ≈ 2 gebraic dependence is recovered at small separations can rewrite the equation above as

λ 2σkBT 1 π 1 Π=2kBTcs = . (13) Π= 2 . (18) · h e0 · h 2βℓB · h This latter equation is similar to the ideal-gas equation of Here, the decay of the pressure is also algebraic. This state, P = NkBT/V , and we will refer to it as the ideal- equation was first derived by Langmuir [19] and we will gas equation [18]. Note that the ideal-gas analogy comes refer to it as the Langmuir equation. Note that the Lang- from the fact that the pressure between the plates in this muir equation features the same 1/h2 dependence as the regime can be calculated by the ideal-gas equation. zero-field result given in Eq. (12), albeit the pre-factor The disjoining pressure in the zero-field approximation and physical origin are different. is derived here for plates with constant charge densities. In the opposite limit, i.e. short separation distances, When one includes charge regulation effects, the resulting where arctan(x) x, Eq. (17) reduces to the ideal- ≈ pressure decay is not affected, but only the amplitude gas equation shown in Eq. (13). Therefore both the weakens [12, 13]. It has been also shown that the Donnan counterion-only and zero-field approximation lead to the potential in Eq. (10) results from the asymptotic solution same result for small separation distances. The ideal-gas of the PB-equation for vanishing electric field [10]. equation can be understood for the counterion-only case in the following way. When plates come sufficiently close together, all the coions are expelled and only counterions 3. Counterion-Only Double-Layers remain. Due to charge neutrality, the number of counterions in the gap must be equal to the number of charges In the counterion only case the disjoining pressure also on the plates. The distance between the plates and the decays algebraically [19, 30, 44]. Here the only ionic area of the plates therefore determine the concentration species are the counterions to the surfaces. If one as- of the ions in the gap as c = 2σ/(e0h); multiplying this sumes a positively charged surface and anions as counte- concentration by kBT leads to the disjoining pressure in rions, the PB equation reduces to Eq. (13). 2 2 For the case of asymmetric z : 1 salts, where z is the d Ψ(x) βe0c− Ψ(x) 2 = e . (14) dx εε0 valence of the coion, the multivalent coions get expelled from the gap between the charged plates upon approach This equation can be solved analytically for two charged [30, 36]. This situation results in a counterion-only sys- plates and the resulting potential profile is tem, however the salt is still present in the bulk. We can

2 xγ correct the Langmuir pressure for the bulk contribution. Ψ(x)=ΨM ln cos , (15) π 1 − ℓ If we add an additional term in the arctan(x) x GC ≈ 2 − 5

6 0 (a) (b) 10 6 (a) (b) 1.0 10 pH 10 -20 3 mM KCl 5 0.8 5 10 -40 10 0.6 exponential -60 4 4 10 10 -80 0.4 pH = 10 pH = 10

3 Pressure (Pa) Pressure (Pa) 10 3 -100 pH = 5.6 0.2 pH = 5.6 exponential 10 pH= 3.0 Regulation Parameter pH= 3.0

Diffuse-Layer Potential (mV) Diffuse-Layer -120 2 0.0 10 2 10 -4 -3 -2 -1 -4 -3 -2 -1 2 4 6 8 2 4 6 8 0 10 20 30 40 -9 -8 -7 10 10 10 10 10 10 10 10 10 10 10 Concentration (M) Concentration (M) Separation Distance (nm) Separation Distance (m)

FIG. 2. (a) Diffuse-layer potential and (b) regulation param- FIG. 3. Comparison of (a) log–lin and (b) log–log representa- eter as a function of KCl concentration for different pH conditions for disjoining pressure between two silica particles mea- tions. The diffuse-layer potentials and regulation parameters sured at pH 10 and 3 mM KCl. The algebraic 1/h and 1/h2, were extracted by fitting experimental force data with full PB as well as exponential curves are also shown. theory Eq. (2) using constant regulation approximation. Fur- ther details on this fitting procedure can be found in Ref. 36. charge at pH 10, where the fraction of charged silanol groups is the highest. On the other hand, regulation pa- expansion we arrive at [31] rameters are decreasing with increasing concentration at π 1 (z + 1)cs pH 10 and pH 5.6, while they are modestly increasing Π= 2 . (19) 2βℓB · (h +2ℓGC) − β at pH 3. These observations are consistent with earlier reports, where a similar decrease in the regulation param- We refer to this approximation as the corrected Langmuir eters was reported for pH 10, while at pH 4 the regulation equation. parameters were constant [24]. Note that regulation parameter of 1 represents constant charge conditions (CC), III. RESULTS AND DISCUSSION while at constant potential (CP) conditions the regulation parameter is 0 and the surface charge varies strongly We use the colloidal probe technique based on AFM when the two surfaces approach each other. Therefore at microscopy to experimentally study algebraic double- high pH and low concentration the silica particles be- layer forces between colloidal particles. Specifically, we have as CC surfaces, while they regulate more strongly measure the double-layer interactions between micron- at low pH and high salt concentrations. Surface proper- sized silica particles in the presence of KCl at different ties of silica reported in Fig. 2 are not the main inter- pH. We further analyze the forces in K4Fe(CN)6 to show est of the present paper, but they are used to determine experimental evidence of counterion-only induced alge- which asymptotic regimes are applicable at certain solu- braic double-layer interactions. tion conditions. The double-layer potentials are further First the van der Waals forces are subtracted from all converted by employing the Graham equation to surface the experimentally measured forces, which yields only the charge density and used to calculate characteristic length double-layer contribution to the force profile. The result- scales, such as Gouy–Chapman length. ing double-layer forces are then converted to disjoining The main focus of the present paper is the analysis pressures between plates by calculating the derivative of of the asymptotic laws in the experimental double-layer the force versus the separation distance. These experi- interactions. In Fig. 3, an example of disjoining pres- mental pressure profiles are then fitted with the full PB sure measured at pH 10 and 3 mM KCl is shown. The theory using the constant regulation approximation and exponential function e−κh, and two algebraic functions the silica surface properties are extracted from the fits. 1/h and 1/h2 are also shown in the figure. Note that Note that the concentrations are fixed to nominal values the amplitudes of these general functions are adjusted during fitting. The diffuse-layer potentials and regulation to fit the experimental data in order to show at which parameters for silica particles extracted from this fitting separation range such dependencies are applicable. The procedure are shown in Fig. 2. As expected for silica sur- data is shown in log–lin and log–log representations for faces in aqueous solutions, the double-layer potentials are clarity. This simple comparison reveals that the experi- increasing with increasing salt concentration due to elec- mental data has an exponential dependence only at sepa- trostatic screening. The particles acquire most negative ration distances above 15 nm, where the double-layers ∼ 6 from the two surfaces are only weakly overlapping. At When pH is increased to 10 (Fig. 4c), the 1/h2 ap- intermediate distances between 5 and 15 nm the double- proximation gets even less accurate since the magnitude layers overlap strongly and the 1/h2 dependence becomes of the charge again increases and consequently κλ further applicable, while at distances below 5 nm the pressure increases. On the other hand, the particles now follow the decreases as 1/h. CC behavior and the ideal-gas 1/h dependence is highly Let us now have a more detailed look at the appli- accurate at small separations. cability of the approximations introduced in the the- Looking at all the sub-figures in Fig. 4, one can de- ory section. First we will focus on strongly overlapping duce that the simple algebraic 1/h2 dependence from double-layers, where the zero-field approximation shown Eq. (12) describes silica–silica interactions at interme- in Eq. (11) is expected. In Fig. 4, the disjoining pres- diate separation distances at low pH and high salt con- sures for different pH conditions and different salt levels centrations, where κλ 1. Note that κλ 1 is not ∼ ≪ are shown. Note that pressures are normalized to bulk reached experimentally and therefore the agreement be- osmotic pressure while distance is normalized with De- tween Eq. (12) and experimental curves is not perfect. bye length. Together with experimental data we plot the The ideal-gas relation Eq. (13) is accurate at small sepa- algebraic zero-field approximations shown in Eqs. (11), ration distances, typically below κh 1, and at low con- ≤ (12), and (13). We further denote the value of κλ in centrations and high pH, where the regulation parameter the sub-figures, which represents the value of the char- is close to 1. Therefore experimental conditions where acteristic length in comparison to the Debye length and double-layers overlap strongly and algebraic interactions permits the evaluation of the accuracy of the zero-field can be reached in 1:1 electrolytes. approximations [10]. Let us now address the comparison of the zero- In the top row of Fig. 4, the data for pH 3 are shown. field approximation Eq. (12) with Langmuir expression, At pH 3 the silica surface exhibits the lowest surface Eq. (18). While the former is valid for strongly over- charge and the regulation parameter is around 1/2 as lapping double-layers in 1:1 electrolyte, the latter is de- evident from Fig. 2. By comparing the experimental rived for counterion-only double-layers. Both of these pressures to the calculated ones, one can deduce that approximations are valid at the intermediate separation the double-layers start to overlap strongly at distances distances and they both feature an inverse square dis- shorter than about 2-3 κh and the zero-field approxima- tance dependence 1/h2. The difference between them tion Eq. (11) becomes accurate. The simplified inverse is in the amplitude, i.e. strength of the interaction. In 2 square separation dependence 1/h shown in Eq. (12) Fig. 5 a comparison of the algebraic interactions is made represents a good approximation at intermediate separa- for weakly and highly charged systems, respectively. The tions at about κh 1 and gets more accurate at higher left sub-figure shows the disjoining pressure at pH 3 and ≈ 2 salt concentrations. The accuracy of the 1/h depen- 20 mM KCl, while the right subfigure shows the interac- dence, Eq. (12), is good at κλ 1, as discussed ear- ≤ tion at pH 10 and 1 mM KCl. In both cases the full PB lier [10], therefore one needs to go to high salt concentra- with constant regulation approximation is used to fit the tions and low pH for it to become applicable. Further- experimental data, and the PB theory perfectly describes more for pH 3 solutions the ideal-gas equation is not accu- the data in the whole distance range. At pH 10 the exper- rate since the surfaces strongly regulate and the ideal-gas imental pressure follows the 1/h ideal-gas dependence at equation is derived for constant charge surfaces. small separations, while at pH 3 this regime is not recov- At intermediate pH values of 5.6 shown in Fig. 4b, the ered due to stronger charge regulation of silica surfaces at zero-field approximation works the best at low concen- this condition. The zero-field 1/h2 dependence, Eq. (12), tration and small separation distances. The simple 1/h2 predicts the behavior rather well at pH 3, while it largely dependence, Eq. (12), gets again more accurate at higher overestimates the repulsion at pH 10. The opposite is concentrations, but it is less accurate compared to pH 3 true for the Langmuir equation, which works better at case, because the magnitude of the charge on silica is pH 10 and it overestimates the pressure at pH 3. These higher at higher pH. The ideal-gas equation works best results are consistent with the range of validity for both at low concentrations, since in this case the regulation pa- models. While the zero-field approximation is valid for rameter is close to 1 and the surface has approximately weakly charged surfaces, the Langmuir equation is appli- constant charge. cable for highly charged surfaces [10, 12, 18]. 7

100 pH 3.0, 1 mM KCl 100 pH 3.0, 10 mM KCl 100 pH 3.0, 20 mM KCl (a)

10 10 10

1 1 1

0.1 0.1 0.1

0.01 0.01 0.01

0.001 0.001 0.001

Normalized Pressure, 4 6 8 2 4 6 8 4 6 8 2 4 6 8 6 8 2 4 6 8 1 10 1 10 1 10 Normalized Separation Distance, Normalized Separation Distance, Normalized Separation Distance,

(b) 100 pH 5.6, 1 mM KCl 100 pH 5.6, 10 mM KCl 100 pH 5.6, 70 mM KCl

10 10 10

1 1 1 Eq. (10) 0.1 Eq. (11) 0.1 0.1 Eq. (12) 0.01 0.01 0.01

0.001 0.001 0.001

Normalized Pressure, 4 6 8 2 4 6 8 4 6 8 2 4 6 8 6 8 2 4 6 8 1 10 1 10 1 10 Normalized Separation Distance, Normalized Separation Distance, Normalized Separation Distance,

100 100 pH 10, 10 mM KCl 100 pH 10, 20 mM KCl (c)

10 10 10

1 1 1

0.1 0.1 0.1

0.01 0.01 0.01

0.001 pH 10, 1 mM KCl 0.001 0.001 4 6 8 2 4 6 8 4 6 8 2 4 6 8 Normalized Pressure, 0.1 1 10 1 10 1 10 Normalized Separation Distance, Normalized Separation Distance, Normalized Separation Distance,

FIG. 4. Zero-Field theory compared to experimental disjoining pressures measured at (a) pH 3.0, (b) pH 5.6, and (c) pH 10. Diﬀerent columns present diﬀerent salt levels. Experimental data is shown with symbols, while full, dashed, and dotted lines represent Eqs. (11), (12), and (13), respectively.

Low Charge 6 6 High Charge 10 10 In these systems, the multivalent ions have the same sign (a) pH 3.0 pH 10 (b) 20 mM KCl 1 mM KCl of charge as the surface. When two charged plates ap- 5 5 10 10 proach each other in these salts, at a certain distance,

4 4 10 10 Ideal Gas the multivalent coions are expelled from the gap between Langmuir

Pressure (Pa) the charged plates [30, 31, 36]. After the expulsion of Pressure (Pa) 3 3 Zero-Field 10 10 Full PB the coions, the systems behaves as counterion-only and Experimental 2 2 10 10 2 4 6 8 2 4 6 8 2 4 6 8 2 4 6 8 the pressure decays algebraically. An example of such -9 -8 -7 -9 -8 -7 10 10 10 10 10 10 Separation Distance (m) Separation Distance (m) interactions between two silica colloids in the presence of K4Fe(CN)6 salt at pH 10 is shown in Fig. 6. Again FIG. 5. Comparison of zero-ﬁeld and Langmuir approxima- the disjoining pressure proﬁle can be perfectly described tion for the (a) low charge and (b) high charge experimental with the full PB model with constant regulation approx- cases. imation. At separation distances above 80 nm, the ∼ pressure decays exponentially, and here both coions and counterions are still present in the gap. By lowering Another important experimental case where algebraic the separation distance, the experimental curve starts interactions are important is the multivalent coion case. 8

6 6 Exponential 10 (a) 10 (b) Algebraic 0.1 mM K4 Fe(CN)6 5 10 5 10 (a) 4 10 Langmuir Saturated Potential 4 10 3 10

Pressure (Pa) Pressure (Pa) 3 2 10 10 Ideal Gas Full PB Langmuir Corrected Langmuir 1 2 Ideal-Gas 10 10 High Charge 2 4 6 8 2 4 6 8 0 20 40 60 80 100 -9 -8 -7 10 10 10 Separation Distance (nm) Separation Distance (m)

FIG. 6. Experimentally measured disjoining pressure in the presence of 0.1 mM K4Fe(CN)6 at pH 10. (a) Log–lin and (b) lin–lin representations of experimental data. The full PB

model for the mixture of 1:z and 1:1 electrolyte is shown Low Charge together with ideal-gas (13), Langmuir (18), and corrected Zero-Field Debye-Hückel Langmuir (19) approximations.

Debye Length/Gouy-Chapmann Length, Normalized Separation Distance, to deviate from the exponential decay and this deviation Langmuir Saturated Potential (b) marks the transition into the counterion-only regime. At distances below 70 nm the four-valent [Fe(CN) ]4− pH = 10 ∼ 6 ions get expelled from the gap, and double-layer is com- prised of only K+ ions. Here the pressure decays algebraically. The Langmuir equation captures this 1/h2 dependence, but it overestimates the magnitude of the Ideal-Gas pH = 5.6 force. If one corrects the Langmuir equation by account- ing for the ions in the bulk and expanding the approximation with one more term, the resulting Eq. (19) describes the data very well in the separation range from pH = 3.0

5 to 80 nm. At distances below 10 nm, the ideal-gas Zero-Field Debye-Hückel equation offers a quantitative description of the data. In the multivalent coion systems where the pressure decays Debye Length/Gouy-Chapmann Length, algebraically, the counterion-only regime is typically im- Normalized Separation Distance, portant in a very wide separation range, and the exponential double-layer forces become evident only at very FIG. 7. Approximation validity map for interaction of charged plates across symmetric electrolytes. (a) Map show- large distances where they are weak. In the presence ing the relevant regimes with equations. (b) Experimentally of multivalent coions, the aggregation of colloids is also accessible ranges for silica in KCl solutions at pH 3.0, pH 5.6, inverse driven by the algebraic forces and leads to the and pH 10 are marked with different colors. Schultze-Hardy rule reported recently by some of us [45]. Similar algebraically decaying interactions were recently measured in the presence of like-charged polyelectrolytes add the ranges where experimental interactions between [31] and between stacked charged bilayers [46]. In the silica colloids were measured. In Fig. 7a, we present the case of polyelectrolytes, the transition between algebraic areas of validity for different approximations with the and exponential regimes is even more abrupt, due to the corresponding equations. The y-axis represents Debye larger charge of the polyelectrolytes as compared to mul- length divided by the Gouy–Chapman length and the tivalent coions. values are increasing with increasing surface charge den- Let us finally give an overview at what conditions the sity and decreasing salt concentration. On the x-axis, presented approximations are valid for symmetric elec- the separation distance relative to the Debye length is trolytes. Approximation validity maps are shown in given. The map can be divided into 5 regions. Note Fig. 7. Similar map was presented in [18] and here we that in reality the borders between the regions are not 9 sharp lines as drawn here, but are more fuzzy. The map IV. CONCLUSIONS can be divided into small separation (κh < 1) and large separation distances (κh > 1), where algebraic and ex- Algebraic double-layer interactions have been demon- ponential approximations are dominant, respectively. At strated experimentally for simple 1:1 salts. This type very small distances and intermediate charge, the ideal- of force is also present for systems containing highly- gas equation (13) is applicable. Note that the ideal-gas charged coions, where they come from the counterion- law is only valid for constant charge surfaces. At inter- only dominated regime. Both strongly overlapping mediate distances and for highly charged particles, the double-layer (zero-field) and counterion-only (Langmuir) Langmuir equation Eq. (18) is accurate. The range of va- regimes recover the inverse square distance dependence lidity can be extended to lower charges by using the cor- of disjoining pressures, albeit with different pre-factor. rected Langmuir approximation Eq. (19). At low charge The former is applicable for weakly charged surfaces and and intermediate distance, the zero-field approximation is harder to detect experimentally, while the latter is true Eq. (12) is applicable. At distances κh > 1, the pressure for highly charged surfaces. Independent of the surface profiles become exponential. Here we can approximate charge, at very small separation distances the 1/h ideal- the low charged case with the Debye–Hückel superpo- gas limit is recovered for non-regulating constant charge sition approximation, while for the highly charged sur- surfaces. In the other limit, at large separation distances, faces one has to replace the diffuse-layer potential in the the Debye–Hückel dependence is observed and pressures Debye–Hückel approximation with the effective saturated follow the conventional exponential profile. The assump- potential [9]. tion of exponential double-layer forces is therefore correct only at large separations, but at small separations the profiles are algebraic. The latter regime is also important in experimental situations, for example for colloidal aggregation in the presence of multivalent coions, or for systems containing like-charged polyelectrolytes. In Fig. 7b we highlight the regions where the measurements between silica colloids were done in the present paper. The position of the range, where the force measurements are located in the map, changes depending ASSOCIATED CONTENT on the salt concentration and pH of the solutions. The AFM measurements can be done reliably down to 1 nm ∼ Supporting Information separation distances, while the limit for large separations depends on the magnitude of the force, which in The Supporting Information is available free of charge turn depends on the charge of the surfaces. Therefore, on the ACS Publications website. the range where the measurements can be performed is moving down and to the right with increasing concen- Experimental force profiles for all pH conditions tration, by increasing κ. On the other hand, increasing and salt levels studied. Table with parameters extracted the pH shifts the ranges up due to increasing magnitude from fitting the experimental forces. of the surface charge, which decreases ℓGC. For all the pH conditions studied, the ideal-gas law is retrieved at low concentrations and sufficiently small separations for CC surfaces, which is consistent with the data presented AUTHOR INFORMATION in Fig. 4. For pH 10 and pH 3, the experimental data falls into the high and low charge regimes, respectively. Corresponding Authors Therefore, the Langmuir approximation is most accurate at pH 10 for intermediate distances, while the zero-field ∗E-mail: [email protected] approximation works best at pH 3. At pH 5.6 none of †E-mail: [email protected] the 1/h2 approximations is very accurate since the surface charge density is in this case in the intermediate Notes regime. The authors declare no competing financial interest. 10

ACKNOWLEDGMENTS Borkovec for providing access to the instruments in his laboratory and to Plinio Maroni for the help with AFM This research was supported by the Swiss National measurements. MV acknowledges the Netherlands Or- Science Foundation through grant 162420 and the Uni- ganisation for Scientiﬁc Research (NWO) for a Veni grant versity of Geneva. The authors are thankful to Michal (no. 722.017.005).

[1] Derjaguin, B. On the Repulsive Forces between Charged [14] Ghosal, S.; Sherwood, J. D. Repulsion Between Finite Colloid Particles and on the Theory of Slow Coagula- Charged Plates with Strongly Overlapped Electric Dou- tion and Stability of Lyophobe Sols. Trans. Faraday Soc. ble Layers. Langmuir 2016, 32, 9445–9450. 1940, 10, 333–346. [15] Ghosal, S.; Sherwood, J. D. Screened Coulomb Interac- [2] Derjaguin, B.; Landau, L. D. Theory of the Stability of tions with Non-Uniform Surface Charge. Proceedings of Strongly Charged Lyophobic Sols and of the Adhesion of the Royal Society A: Mathematical, Physical and Engi- Strongly Charged Particles in Solutions of Electrolytes. neering Sciences 2017, 473, 20160906. Acta Phys Chim 1941, 14, 633–662. [16] Gouy, M. Sur la constitution de la charge électrique àla [3] Verwey, E. J. W.; Overbeek, J. T. G. Theory of Stability surface d’un électrolyte. J. Phys. Theor. Appl. 1910, 9, of Lyophobic Colloids; Elsevier: Amsterdam, 1948. 457–468. [4] Kruyt, H. R.; Jonker, G. H.; Overbeek, J. T. G. Col- [17] Chapman, D. L. A Contribution to the Theory of Elec- loid Science, Edited by H.R. Kruyt. [Contributors: G.H. trocapillarity. Philos. Mag. 1913, 25, 475–481. Jonker and J. Th. G. Overbeek. English Translation by [18] Markovich, T.; Andelman, D.; Podgornik, R. Charged L.C. Jackson.] Vol. I. Irreversible Systems; Elsevier, Membranes: Poisson-Boltzmann Theory, DLVO 1952. Paradigm and Beyond. ArXiv160309451 Cond-Mat [5] Evans, D. F.; Wennerstrom, H. The Colloidal Domain; 2016, John Wiley: New York, 1999. [19] Langmuir, I. The Role of Attractive and Repulsive Forces [6] Russel, W. B.; Saville, D. A.; Schowalter, W. R. Col- in the Formation of Tactoids, Thixotropic, Gels Protein loidal Dispersions; Cambridge University Press: Cam- Crystals and Coacervates. J. Chem. Phys. 1938, 6, 873– bridge, 1989. 896. [7] Lyklema, J. Fundamentals of Interface and Colloid Sci- [20] Kjellander, R. Intricate Coupling between Ion-Ion and ence: Particulate Colloids; Elsevier, 2005. Ion-Surface Correlations in Double Layers as Illus- [8] Israelachvili, J. Intermolecular and Surface Forces, 3rd trated by Charge Inversion-Combined Effects of Strong ed.; Academic Press: London, 2011. Coulomb Correlations and Excluded Volume. J. Phys.- [9] Elimelech, M.; Gregory, J.; Jia, X. Particle Deposition Condens. Matter 2009, 21, 424101. and Aggregation: Measurement, Modelling and Simula- [21] Attard, P.; Mitchell, D. J.; Ninham, B. W. Beyond Pois- tion; Butterworth-Heinemann, 2013. son–Boltzmann: Images and Correlations in the Electric [10] Philipse, A. P.; Kuipers, B. W. M.; Vrij, A. Algebraic Double Layer. I. Counterions Only. J. Chem. Phys. 1988, Repulsions between Charged Planes with Strongly Over- 88, 4987–4996. lapping Electrical Double Layers. Langmuir 2013, 29, [22] Finlayson, S. D.; Bartlett, P. Non-Additivity of Pair In- 2859–2870. teractions in Charged Colloids. J. Chem. Phys. 2016, [11] Philipse, A.; Vrij, A. The Donnan Equilibrium: I. On the 145, 034905. Thermodynamic Foundation of the Donnan Equation of [23] Sainis, S. K.; Merrill, J. W.; Dufresne, E. R. Electro- State. J. Phys.: Condens. Matter 2011, 23, 194106. static Interactions of Colloidal Particles at Vanishing [12] Philipse, A. P.; Tuinier, R.; Kuipers, B. W. M.; Vrij, A.; Ionic Strength. Langmuir 2008, 24, 13334–13337. Vis, M. On the Repulsive Interaction Between Strongly [24] Smith, A. M.; Maroni, P.; Borkovec, M.; Trefalt, G. Mea- Overlapping Double Layers of Charge-Regulated Sur- suring Inner Layer Capacitance with the Colloidal Probe faces. Colloid and Interface Science Communications Technique. Colloids Interfaces 2018, 2, 65. 2017, 21, 10–14. [25] Israelachvili, J. N.; Adams, G. E. Measurement of Forces [13] Vis, M.; Tuinier, R.; M. Kuipers, B. W.; Vrij, A.; between Two Mica Surfaces in Aqueous Electrolyte Solu- P. Philipse, A. Interactions between Amphoteric Surfaces tions in Range 0-100 Nm. J. Chem. Soc.-Faraday Trans. with Strongly Overlapping Double Layers. Soft Matter I 1978, 74, 975–1001. 2018, 14, 4702–4710. [26] Prieve, D. C. Measurement of Colloidal Forces with TIRM. Adv. Colloid Interface Sci. 1999, 82, 93–125. 11

[27] Popa, I.; Sinha, P.; Finessi, M.; Maroni, P.; Papas- [36] Uzelac, B.; Valmacco, V.; Trefalt, G. Interactions be- tavrou, G.; Borkovec, M. Importance of Charge Regu- tween Silica Particles in the Presence of Multivalent lation in Attractive Double-Layer Forces between Dis- Coions. Soft Matter 2017, 13, 5741–5748. similar Surfaces. Phys Rev Lett 2010, 104, 228301. [37] Sader, J. E.; Chon, J. W. M.; Mulvaney, P. Calibration of [28] Smith, A. M.; Lee, A. A.; Perkin, S. The Electrostatic Rectangular Atomic Force Microscope Cantilevers. Rev. Screening Length in Concentrated Electrolytes Increases Sci. Instrum. 1999, 70, 3967–3969. with Concentration. J. Phys. Chem. Lett. 2016, 7, 2157– [38] Avni, Y.; Andelman, D.; Podgornik, R. Charge Regu- 2163. lation with Fixed and Mobile Charged Macromolecules. [29] Briscoe, W. H.; Horn, R. G. Direct Measurement of Sur- Current Opinion in Electrochemistry 2019, 13, 70–77. face Forces Due to Charging of Solids Immersed in a Non- [39] Behrens, S. H.; Borkovec, M. Electric Double Layer In- polar Liquid. Langmuir 2002, 18, 3945–3956. teraction of Ionizable Surfaces: Charge Regulation for [30] Montes Ruiz-Cabello, F. J. M.; Moazzami-Gudarzi, M.; Arbitrary Potentials. J. Chem. Phys. 1999, 111, 382– Elzbieciak-Wodka, M.; Maroni, P.; Labbez, C.; 385. Borkovec, M.; Trefalt, G. Long-Ranged and Soft Inter- [40] Kirkwood, J. G.; Shumaker, J. B. Forces between Pro- actions between Charged Colloidal Particles Induced by tein Molecules in Solution Arising from Fluctuations in Multivalent Coions. Soft Matter 2015, 11, 1562–1571. Proton Charge and Conﬁguration. Proc. Natl. Acad. Sci. [31] Moazzami-Gudarzi, M.; Kremer, T.; Valmacco, V.; U. S. A. 1952, 38, 863–871. Maroni, P.; Borkovec, M.; Trefalt, G. Interplay be- [41] Markovich, T.; Andelman, D.; Podgornik, R. Charge tween Depletion and Double-Layer Forces Acting be- Regulation: A Generalized Boundary Condition? EPL tween Charged Particles in Solutions of Like-Charged 2016, 113, 26004. Polyelectrolytes. Phys. Rev. Lett. 2016, 117, 088001. [42] Adˇzi´c, N.; Podgornik, R. Field-Theoretic Description of [32] Borkovec, M.; Szilagyi, I.; Popa, I.; Finessi, M.; Sinha, P.; Charge Regulation Interaction. Eur. Phys. J. E 2014, Maroni, P.; Papastavrou, G. Investigating Forces be- 37, 49. tween Charged Particles in the Presence of Oppositely [43] Trefalt, G.; Behrens, S. H.; Borkovec, M. Charge Regula- Charged Polyelectrolytes with the Multi-Particle Col- tion in the Electrical Double Layer: Ion Adsorption and loidal Probe Technique. Adv. Colloid Interface Sci. 2012, Surface Interactions. Langmuir 2016, 32, 380–400. 179-182, 85–98. [44] Briscoe, W. H.; Attard, P. Counterion-Only Electrical [33] Butt, H. J.; Cappella, B.; Kappl, M. Force Measurements Double Layer: A Constrained Entropy Approach. J. with the Atomic Force Microscope: Technique, Interpre- Chem. Phys. 2002, 117, 5452–5464. tation and Applications. Surf. Sci. Rep. 2005, 59, 1–152. [45] Cao, T.; Szilagyi, I.; Oncsik, T.; Borkovec, M.; Trefalt, G. [34] Butt, H. J. Measuring Electrostatic, van Der Waals, and Aggregation of Colloidal Particles in the Presence of Mul- Hydration Forces in Electrolyte Solutions with an Atomic tivalent Co-Ions: The Inverse Schulze–Hardy Rule. Lang- Force Microscope. Biophys. J. 1991, 60, 1438–1444. muir 2015, 31, 6610–6614. [35] Ducker, W. A.; Senden, T. J.; Pashley, R. M. Direct Mea- [46] Hishida, M.; Nomura, Y.; Akiyama, R.; Yamamura, Y.; surement of Colloidal Forces Using an Atomic Force Mi- Saito, K. Electrostatic Double-Layer Interaction between croscope. Nature 1991, 353, 239–241. Stacked Charged Bilayers. Phys. Rev. E 2017, 96, 040601.