Open Access Article

This Open Access Article is licensed under a

Creative Commons Attribution 3.0 Unported Licence

Arghya
Majee
*^{ab},
Markus
Bier
^{ab} and
Rudolf
Podgornik
^{cd}
^{a}Max Planck Institute for Intelligent Systems, Heisenbergstr. 3, 70569 Stuttgart, Germany. E-mail: majee@is.mpg.de
^{b}Institute for Theoretical Physics IV, University of Stuttgart, Pfaffenwaldring 57, 70569 Stuttgart, Germany
^{c}Department of Theoretical Physics, J. Stefan Institute, Jamova c. 39, 1000 Ljubljana, Slovenia
^{d}Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, 1000 Ljubljana, Slovenia

Received
17th November 2017
, Accepted 2nd January 2018

First published on 2nd January 2018

The interaction between two chemically identical charge-regulated surfaces is studied using the classical density functional theory. In contrast to common expectations and assumptions, under certain realistic conditions we find a spontaneous emergence of disparate charge densities on the two surfaces. The surface charge densities can differ not only in their magnitude, but quite unexpectedly, even in their sign, implying that the electrostatic interaction between the two chemically identical surfaces can be attractive instead of repulsive. Moreover, an initial symmetry with equal charge densities on both surfaces can also be broken spontaneously upon decreasing the separation between the two surfaces. The origin of this phenomenon is a competition between the adsorption of ions from the solution to the surface and the interaction between the adsorbed ions already on the surface. These findings are fundamental for the understanding of the forces between colloidal objects and, in particular, they are bound to strongly influence the present picture of protein interaction.

Studies of the interaction between two charge-regulated surfaces have been performed for chemically identical surfaces with equal adsorption and desorption properties^{8–11} as well as for chemically non-identical surfaces.^{2,12–15} In all of the former cases, a certain basic symmetry of the problem was assumed a priori^{16–18} and the surface charge densities have been without exception constrained to be equal on both surfaces. However, the underlying physical reasoning for such an assumption is not generally applicable and is not based upon the detailed chemical nature of the surfaces bearing charge. The fact that the two interacting surfaces are chemically identical and, therefore, interact in the same way with the adjacent bathing solution, is not sufficient to infer the fundamental charge symmetry and to invoke equal surface charge densities in the application of the PB formalism. In fact, the charge distribution of the system should not be assumed a priori, but should follow from the minimization of the relevant total thermodynamic potential, yielding the equilibrium state in terms of the equilibrium electrostatic potential distribution between the surfaces as well as the equilibrium charge densities on the surfaces, without any additional assumptions. Whether this minimum implies an equal or unequal surface charge densities may and, as will be shown below, does depend on the parameters of the system under consideration.

Below we show that, depending on these parameters, a confined electrolyte in thermodynamic equilibrium with two chemically identical charge-regulated surfaces that can adsorb/desorb solution ions, can indeed adopt unequal surface charge densities, even at separations that are much larger than the Debye length. This happens due to an interplay of the adsorption of ions from the solution to the surface and the interaction between the already adsorbed ions at the surface. The model surface free energy is then related to the lattice fluid model and is composed of the surface entropy of mixing, the electrostatic energy of the adsorbed charges, the non-electrostatic energy penalty of adsorption and the change in the non-electrostatic interactions between the ions upon adsorption. The latter are assumed to be short-ranged, typically of the van der Waals type, hydrogen-bonding and/or of quantum-chemical origin, that allow for a nearest-neighbor-like description. Since in general the surface charge densities can differ in magnitude as well as in sign, an initial symmetry with equal charge densities on both surfaces can be spontaneously broken, and the surfaces can acquire different charge densities as they approach each other. At short separations, this implies surface charge densities differing even in their sign and consequently leading to an overall attractive interaction between the surfaces. An analytical treatment of the simpler system with only a single surface in contact with an electrolyte indicates that these findings are inherently related solely to the electrostatic interaction between the surfaces.

(1) |

(2) |

Fig. 2 Absolute value of the difference of the dimensionless charge densities at the two surfaces, |σ_{0}* − σ_{L}*|, as function of α and χ (see eqn (2)) for (a) κL = 0.1, (b) 1, (c) 1.5, and (d) 3. The white regions in each figure correspond to a symmetric situation where σ_{0}* = σ_{L}* whereas in the colored regions the equilibrium values of σ_{0}* and σ_{L}* differ from one another. |

First, we consider the case κL = 1. As shown by Fig. 2(b), σ_{0}* and σ_{L}* are the same over a broad region (indicated by white) but not everywhere. In the dark blue region, the charge asymmetry is the highest and close to unity, implying that the two surfaces are oppositely charged. The line χ = −2α passes through the middle of this region. Along this line, the solution is σ_{0}* = σ_{L}* = 0 for α > α_{0} (≈−13.2 in this case) and for α ≲ α_{0} the dark blue region appears. For α ≳ α_{0} there are two more regions (one below the line χ = −2α and one above) where the charge asymmetry is present albeit with a lower contrast |σ_{0}* − σ_{L}*|. These two tails (light blue or greenish) are not inter-connected but with decreasing α they thicken, come close to each other, and finally merge with the dark blue region. The charge contrast in each of these two regions increases with decreasing α. Below the dark blue region and the lower tail, the equilibrium states are symmetric with σ_{0}* = σ_{L}* ≳ −½ whereas above the dark blue region and the upper tail, they are symmetric with σ_{0}* = σ_{L}* ≲ ½. In between the two tails the states are σ_{0}* = σ_{L}* ≈ 0. In other words, (σ_{0}*,σ_{L}*) changes from (≳−½,≳−½) to (≲½,≲½) across the dark blue region, from (≳−½,≳−½) to (≈0,≈0) across the lower tail, and from (≈0,≈0) to (≲½,≲½) across the upper tail.

With decreasing separation L between the two surfaces, the dark blue region in Fig. 2 broadens and starts to dominate over the tails, making them hardly visible (see Fig. 2(a)). With increasing separation all the regions shrink (see Fig. 2(c) and (d)) and the tails become increasingly difficult to be resolved numerically. In Fig. 2 both the interaction parameters α and χ are sampled with a tenth of the thermal energy k_{B}T = 1/β because ion adsorption is governed by a competition with the bulk solvation free energy, which can usually be measured within a similar accuracy.^{28,29} However, the dark blue region seems to be very stable and it remains present even at κL = 10. Upon increasing κL from 0.1 to 3 both α_{0} and the width of the dark blue region decrease relatively fast, whereas for κL between 3 and 10, they hardly change.

In order to explain these findings we consider a single surface in contact with an electrolyte. This problem is analytically solvable and the solution shows that on the line χ = −2α, there are two equally deep (local) minima of the grand potential Ω(σ*) at σ_{1}* and σ_{2}* = −σ_{1}*.^{26} For α ≳ −14, σ* = 0 corresponds to the single global minimum of the grand potential; see Fig. 3(a). However, for α ≲ −14, the global minimum shifts to states with σ_{1}* ≳ −½ and σ_{2}* = −σ_{1}* ≲ ½; see Fig. 3(b). For α ≲ −14, in the presence of a second surface, one surface acquires the charge density σ* = σ_{1}* and the other σ* = σ_{2}* = −σ_{1}* because the electrostatic attraction of two oppositely charged surfaces leads to a decrease of the grand potential of the system. Similarly, the one-surface problem shows equally deep global minima of Ω(σ*) at σ* ≳ −½ and σ* ≲ 0 for points in the upper part of the lower tail of Fig. 2(b) (e.g., see Fig. 3(c)). A second surface leads to charge densities of σ* ≳ −½ on one surface and of σ* ≳ 0 on the other such that the free energy cost in going upward the curve in Fig. 3(c) is balanced by a reduction of the free energy due to electrostatic attraction. As we go down the lower tail of Fig. 2(b), the one-surface problem shows two unequally deep minima at σ* ≳ −½ and σ* ≲ 0 for these points; see Fig. 3(d). Although for a single surface the minimum at σ* ≳ −½ is slightly deeper than the one at σ* ≲ 0, the combination σ_{0}* = σ_{L}* ≳ −½ for two surfaces would be too expensive due to strong electrostatic repulsion and σ_{0}* = σ_{L}* ≲ 0 is also a state of higher free energy. The balance for two surfaces is obtained for the combination (≳−½,≲0) by avoiding the repulsive interaction energy. As is shown in Fig. 3(e) and (f), a similar phenomenon occurs for the upper tail in Fig. 2(b) except for the fact that there the equilibrium states are at σ* ≳ 0 and ≲½. With increasing separation κL, the regions with charge asymmetry shrink because of a weaker electrostatic interaction due to screening.

Fig. 3 Variation of Δ(σ*) (see eqn (19) of the ESI, ref. 26), which is obtained after subtracting the bulk contribution from the grand potential functional per unit surface area for a system consisting of a single charge-regulated surface in contact with an electrolyte and expressed in the units of 1/(a^{2}β), as function of σ* for different combinations of the parameters α and χ. In all cases λ = βe^{2}/(4εκa^{2}) ≈ 21.3 is used. |

Once σ*(r) is known, the grand potential per unit surface area of the system can be obtained by evaluating [σ*]; see eqn (2). The dependence of as function of the separation κL describes an effective interaction potential between the surfaces and is shown in Fig. 4 for different combinations of the parameters α and χ. In each case, increases initially with increasing κL and shows a maximum at some finite separation, typically well above the molecular length scale. Upon increasing κL further, the electrostatic interaction vanishes exponentially ∼exp(−κL) (see Fig. 5) and the osmotic (or entropic) contribution (=−2IL/β ∼ L) of the ions to dominates. The effective force per unit surface area −∂/∂L is negative up to the distance L of the maximum in Fig. 4 and therefore, the interaction is attractive. This implies that the electrostatic attraction is sufficiently strong to overcome the repulsive osmotic pressure. At larger separations, however, the electrostatic interaction weakens and in the limit L → ∞, the effective force per unit surface area equals the constant osmotic pressure 2I/β. The occurrence of the maximum in Fig. 4(a) at a larger separation compared to the cases in Fig. 4(b) and (c) is related to enhanced electrostatic interactions due to the stronger asymmetry in the surface charge densities. Note that an effective interaction potential is a mesoscopic concept, which incorporates the energy and entropy balance of all microscopic degrees of freedom, e.g., the surface charge densities and the ion number density profiles, by minimization of the microscopic grand potential under the constraint of fixed mesoscopic degrees of freedom, e.g., the wall separation. In an even more microscopic (atomistic) approach, one could attempt to replace the parameters α and χ in favor of free energy contributions of the corresponding processes. In that sense an effective interaction energy always contains both energetic and entropic contributions.

Fig. 4 Variation of the effective interaction potential per unit surface area between two charge-regulated surfaces in the units of k_{B}T per nm^{2} as function of the scaled separation κL for three different state points in Fig. 2(b) from (a) the dark blue region, (b) the lower tail and (c) the upper tail. As shown by the plots, the interaction energy increases initially, shows a maximum, and ultimately decreases monotonically. The initial increase in the interaction energy corresponds to a negative effective force at short distances L, implying that the interaction is attractive there. The value of κL at the maximum of is indicated by arrows. As expected, the attraction in the dark blue region of Fig. 2 is stronger than in the other colored regions due to a larger contrast in the surface charge densities. |

Fig. 5 Electrostatic part _{el} of the effective interaction potential per unit surface area between two charge-regulated surfaces in units of k_{B}T per nm^{2} as function of the scaled separation κL for three different state points as in Fig. 4. For α = −17.0 and χ = 34.0, the interaction is attractive everywhere implied by the opposite signs of the charge densities at the two surfaces. For the other two cases, the electrostatic interaction energies increase initially, show a maximum, and then decay to zero. Both the curves show kinks within the repulsive part of the interaction which are related to the discontinuities of the surface charge densities as functions of the wall separation κL (see Fig. 2). As expected, the interaction decays exponentially to zero which is confirmed by the semi-logarithmic plot in the inset. |

The electrostatic part _{el} of the total interaction energy after subtracting the ideal osmotic (or entropic) contribution (= −2IL/β ∼ L) of the ions and the surface tensions acting at the two solid–liquid interfaces, is shown as a function of the separation κL in Fig. 5 for the same values of the α and the χ parameters as in Fig. 4. For α = −17.0 and χ = 34.0, which is a point on the line χ = −2α in Fig. 2, the interaction is attractive due to the opposite charge densities at the two surfaces everywhere within the range of κL shown here. For the other two cases, the interactions are attractive at short distances κL and they become repulsive with increasing separation. These two curves show kinks corresponding to the discontinuities of the surface charge densities as functions of the wall separation κL (see Fig. 2). For example, if one considers the curve corresponding to α = −13.0 and χ = 24.1, the two surfaces are oppositely charged up to κL ≈ 1.5, then, from κL ≈ 1.5 to κL ≈ 1.9, the two surfaces carry equal charge densities (σ_{0}* = σ_{L}* ≲ ½), from κL ≈ 1.9 to κL ≈ 3.4 the surfaces adopt charge densities which are different in magnitude but have equal signs (σ_{0}* ≈ 0, σ_{L}* ≲ ½), and finally, beyond κL ≈ 3.4, both surfaces become equally charged (σ_{0}* = σ_{L}* ≈ 0). A similar phenomenon occurs for α = −12.9 and χ = 26.5 where the interaction changes from attractive to repulsive at κL ≈ 2.9, and at κL ≈ 3.2 the charge densities at the two surfaces become equal. As expected, at large separations, the electrostatic interaction in all cases decay exponentially ∼exp(−κL); see the inset of Fig. 5. Please note that the equilibrium states are characterized by a minimum of the total interaction potential as function of the wall separation L. As (L) corresponds to the minimum of the grand potential functional Ω[σ*,n_{±}] under the constraint of a fixed wall separation L, it is necessarily continuous with respect to L, but its derivatives with respect to L may be discontinuous at first-order phase transitions. In the present work no first-order bulk phase transitions are considered but the observed spontaneous symmetry breaking of the surface charge densities corresponds to first-order surface phase transitions. Hence, kinks can occur only in the surface contribution _{el}(L), and they are hardly visible in the total interaction (L) (note the widely different scales in Fig. 4 and 5).

The inter-surface force is usually measured by using a surface forces apparatus (SFA), atomic force microscopy (AFM), or optical tweezers.^{27,30} In order to observe the anomalous attraction discussed in the preceding paragraph for a surface with an appropriate charge regulation behavior, one can either fix the distance between the surfaces and change the ionic strength of the solution or fix the ionic strength and vary the separation between the surfaces. For relatively small separations compared to the Debye length, L ≲ 1/κ, the system is expected to exhibit a broad parameter range of surface charge asymmetry (see, e.g., the colored region in Fig. 2(a)). Charge asymmetry is also present for larger separations L, but the corresponding parameter range is smaller and more difficult to find (see Fig. 2). In order to avoid possible additional effects occurring at short separations L in a real experimental setup, it is advisable to use low ionic strengths, i.e., large Debye lengths 1/κ. For example, I = 0.1 mM in water leads to 1/κ ≈ 30 nm, so that κL ≈ 0.1 for a separation length L ≈ 3 nm, which is much larger than molecular dimensions and, therefore, a mean-field-like theory as the one presented here is expected to work well.^{30} Moreover, it is not necessary to go to such small values of κL: between κL = 0.1 and 1 (e.g., κL = 0.5), one can expect to have a sufficiently broad parameter range of surface charge asymmetry (see the colored regions in Fig. 2). Possible candidates for the type of surfaces described here are biomolecules like lipids or proteins^{27} or solid colloidal particles (e.g., made of silica) grafted with particular surface groups (e.g., −NH_{2} or −COOH) (see ref. 31 and 32). The parameter χ can be adjusted by means of an appropriate arrangement and density of surface groups. On the other hand, the parameter α can be tuned by changing counterion concentration, e.g., the pH, in the solvent. As mentioned earlier, the parameters α and χ can be obtained, e.g., by fitting experimentally measured profiles of the effective force. For example, for the synthetic cationic double-chain surfactant didodecyldimethylammonium (DDA^{+}) the values α = −7.4, χ = 14.75 for bromide (Br^{−}) and α = −3.4, χ = 14.75 for chloride (Cl^{−}) counterions have been obtained in ref. 22. This demonstrates that the parameter ranges for α and χ, for which surface charge asymmetry is predicted here, are experimentally accessible, in particular in the case of low ionic strengths.

We finish our discussion by briefly commenting on the importance of the possible electrostatic attraction due to surface charge asymmetry in comparison with the van der Waals (vdW) attraction present in the system. The vdW interaction is usually estimated in terms of the Hamaker coefficient.^{33} For a pair of parallel planar silica surfaces interacting across water, the Hamaker coefficient is A ≈ 4.8 × 10^{−21} J (see ref. 34), so that the vdW attraction energy per unit surface area −A/(12πL^{2}) ≈ 0.03k_{B}T per nm^{2} for L = 1 nm, which corresponds roughly to the thickness of the lines in Fig. 4. Hence the vdW interaction is qualitatively and quantitatively irrelevant for the effective interaction potentials considered in Fig. 4. The same can be expected for biological molecules, where the Hamaker coefficients are typically similar or smaller than 1k_{B}T (see ref. 35).

It is important to note that our findings are not restricted to the case Θ = 2. For an asymmetric charge interval corresponding to Θ ≠ 2, the colored regions of Fig. 2 shift in the α−χ plane but the qualitative features remain the same. In fact, we obtain asymmetric equilibrium states σ_{0}* ≠ σ_{L}* even for Θ = 1, where both surfaces can acquire only negative charges or remain uncharged; there the origin of asymmetry is similar to the greenish regions in Fig. 2.

- T. Markovich, D. Andelman and R. Podgornik, in Handbook of Lipid Membranes, ed. C. Safynia and J. Raedler, Taylor & Francis, 2018 Search PubMed.
- I. Popa, P. Sinha, M. Finessi, P. Maroni, G. Papastavrou and M. Borkovec, Phys. Rev. Lett., 2010, 104, 228301 CrossRef PubMed.
- M. Lund and B. Jönsson, Q. Rev. Biophys., 2013, 46, 265 CrossRef CAS PubMed.
- B. W. Ninham and V. A. Parsegian, J. Theor. Biol., 1971, 31, 405 CrossRef CAS PubMed.
- H. Diamant and D. Andelman, J. Phys. Chem., 1996, 100, 13732 CrossRef CAS.
- G. Longo, M. O. de la Cruz and I. Szleifer, ACS Nano, 2013, 7, 2693 CrossRef CAS PubMed.
- A. C. Maggs and R. Podgornik, EPL, 2014, 108, 68003 CrossRef.
- S. H. Behrens and M. Borkovec, J. Phys. Chem. B, 1999, 103, 2918 CrossRef CAS.
- P. M. Biesheuvel, J. Colloid Interface Sci., 2004, 275, 514 CrossRef CAS PubMed.
- M. Borkovec and S. H. Behrens, J. Phys. Chem. B, 2008, 112, 10796 CrossRef PubMed.
- G. Trefalt, S. H. Behrens and M. Borkovec, Langmuir, 2016, 32, 380 CrossRef CAS PubMed.
- D. Chan, T. W. Healy and L. R. White, J. Chem. Soc., Faraday Trans., 1976, 72, 2844 RSC.
- D. McCormack, S. L. Carnie and D. Y. C. Chan, J. Colloid Interface Sci., 1995, 169, 177 CrossRef CAS.
- S. H. Behrens and M. Borkovec, Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top., 1999, 60, 7040 CrossRef CAS.
- D. Y. C. Chan, T. W. Healy, T. Supasiti and S. Usui, J. Colloid Interface Sci., 2006, 296, 150 CrossRef CAS PubMed.
- J. E. Sader and D. Y. C. Chan, J. Colloid Interface Sci., 1999, 213, 268 CrossRef CAS PubMed.
- J. C. Neu, Phys. Rev. Lett., 1999, 82, 1072 CrossRef CAS.
- E. Trizac and J.-L. Raimbault, Phys. Rev. E: Stat. Phys., Plasmas, Fluids, Relat. Interdiscip. Top., 1999, 60, 6530 CrossRef CAS.
- T. Hiemstra, W. H. V. Riemsdijk and G. H. Bolt, J. Colloid Interface Sci., 1989, 133, 91 CrossRef CAS.
- T. Hiemstra, J. C. M. D. Wit and W. H. V. Riemsdijk, J. Colloid Interface Sci., 1989, 133, 105 CrossRef CAS.
- M. Kosmulski, Langmuir, 2005, 21, 7421 CrossRef CAS PubMed.
- D. Harries, R. Podgornik, V. A. Parsegian, E. Mar-Or and D. Andelman, J. Chem. Phys., 2006, 124, 224702 CrossRef PubMed.
- T. Markovich, D. Andelman and R. Podgornik, Europhys. Lett., 2016, 113, 26004 CrossRef.
- N. Adžić and R. Podgornik, Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2015, 91, 022715 CrossRef PubMed.
- M. Dubois, T. Zemb, N. Fuller, R. P. Rand and V. A. Parsegian, J. Chem. Phys., 1998, 108, 7855 CrossRef CAS.
- See the ESI for details of the calculations.
- J. N. Israelachvili, Intermolecular and Surface Forces, Academic Press, Amsterdam, 2011 Search PubMed.
- Y. Marcus, Pure Appl. Chem., 1983, 55, 977 CrossRef CAS.
- H. D. Inerowicz, W. Li and I. Persson, J. Chem. Soc., Faraday Trans., 1994, 90, 2223 RSC.
- G. Trefalt, T. Palberg and M. Borkovec, Curr. Opin. Colloid Interface Sci., 2017, 27, 9 CrossRef CAS.
- N. Solin, L. Han, S. Che and O. Terasaki, Catal. Commun., 2009, 10, 1386 CrossRef CAS.
- L. Yuan, L. Chen, X. Chen, R. Liu and G. Ge, Langmuir, 2017, 33, 8724 CrossRef CAS PubMed.
- H. C. Hamaker, Physica, 1937, 4, 1058 CrossRef CAS.
- L. Bergström, Adv. Colloid Interface Sci., 1997, 70, 125 CrossRef.
- D. Leckband and J. Israelachvili, Q. Rev. Biophys., 2001, 34, 105 CrossRef CAS PubMed.

## Footnote |

† Electronic supplementary information (ESI) available. See DOI: 10.1039/c7sm02270k |

This journal is © The Royal Society of Chemistry 2018 |