DLVO theory

File:DLVO-Theorie.png

LVO-Theorie

In physical chemistry, the Derjaguin–Landau–Verwey–Overbeek (DLVO) theory explains the aggregation and kinetic stability of aqueous dispersions quantitatively and describes the force between charged surfaces interacting through a liquid medium. It combines the effects of the van der Waals attraction and the electrostatic repulsion due to the so-called double layer of counterions. The electrostatic part of the DLVO interaction is computed in the mean field approximation in the limit of low surface potentials - that is when the potential energy of an elementary charge on the surface is much smaller than the thermal energy scale, <math> k_\text{B} T</math>. For two spheres of radius <math>a</math> each having a charge <math>Z</math> (expressed in units of the elementary charge) separated by a center-to-center distance <math>r</math> in a fluid of dielectric constant <math>\epsilon_r</math> containing a concentration <math>n</math> of monovalent ions, the electrostatic potential takes the form of a screened-Coulomb or Yukawa potential,

<math display="block">\beta U(r) = Z^2 \lambda_\text{B} \, \left(\frac{e^{\kappa a}}{1 + \kappa a}\right)^2 \, \frac{e^{-\kappa r}}{r}, </math>

where

• <math>\lambda_\text{B}</math> is the Bjerrum length,
• <math>U</math> is the potential energy,
• <math>e</math> ≈ 2.71828 is Euler's number,
• <math>\kappa</math> is the inverse of the Debye–Hückel screening length (<math>\lambda_\text{D}</math>); <math>\kappa</math> is given by <math>\kappa^2 = 4 \pi \lambda_\text{B} n</math>, and
• <math>\beta^{-1} = k_\text{B} T</math> is the thermal energy scale at absolute temperature <math>T</math>

The DLVO theory is named after Boris Derjaguin and Lev Landau, Evert Verwey and Theodoor Overbeek who developed it between 1941 and 1948.

Overview

DLVO theory is a theory of colloidal dispersion stability in which zeta potential is used to explain that as two particles approach one another their ionic atmospheres begin to overlap and a repulsion force is developed.<ref>Template:Cite book</ref> In this theory, two forces are considered to impact on colloidal stability: Van der Waals forces and electrical double layer forces.

The total potential energy is described as the sum of the attraction potential and the repulsion potential. When two particles approach each other, electrostatic repulsion increases and the interference between their electrical double layers increases. However, the Van der Waals attraction also increases as they get closer. At each distance, the net potential energy of the smaller value is subtracted from the larger value.<ref name=":0">{{#invoke:citation/CS1|citation |CitationClass=web }}</ref>

At very close distances, the combination of these forces results in a deep attractive well, which is referred to as the primary minimum. At larger distances, the energy profile goes through a maximum, or energy barrier, and subsequently passes through a shallow minimum, which is referred to as the secondary minimum.<ref name=":1">{{#invoke:citation/CS1|citation |CitationClass=web }}</ref>

At the maximum of the energy barrier, repulsion is greater than attraction. Particles rebound after interparticle contact, and remain dispersed throughout the medium. The maximum energy needs to be greater than the thermal energy. Otherwise, particles will aggregate due to the attraction potential.<ref name=":1" /> The height of the barrier indicates how stable the system is. Since particles have to overcome this barrier in order to aggregate, two particles on a collision course must have sufficient kinetic energy due to their velocity and mass.<ref name=":0" /> If the barrier is cleared, then the net interaction is all attractive, and as a result the particles aggregate. This inner region is often referred to as an energy trap since the colloids can be considered to be trapped together by Van der Waals forces.<ref name=":0" />

For a colloidal system, the thermodynamic equilibrium state may be reached when the particles are in deep primary minimum. At primary minimum, attractive forces overpower the repulsive forces at low molecular distances. Particles coagulate and this process is not reversible.<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref> However, when the maximum energy barrier is too high to overcome, the colloid particles may stay in the secondary minimum, where particles are held together but more weakly than in the primary minimum.<ref>Template:Cite journal</ref> Particles form weak attractions but are easily redispersed. Thus, the adhesion at secondary minimum can be reversible.<ref>{{#invoke:citation/CS1|citation |CitationClass=web }}</ref>

History

In 1923, Peter Debye and Erich Hückel reported the first successful theory for the distribution of charges in ionic solutions.<ref> Template:Citation </ref> The framework of linearized Debye–Hückel theory subsequently was applied to colloidal dispersions by S. Levine and G. P. Dube<ref> Template:Citation </ref><ref> Template:Citation </ref> who found that charged colloidal particles should experience a strong medium-range repulsion and a weaker long-range attraction. This theory did not explain the observed instability of colloidal dispersions against irreversible aggregation in solutions of high ionic strength. In 1941, Boris Derjaguin and Lev Landau introduced a theory for the stability of colloidal dispersions that invoked a fundamental instability driven by strong but short-ranged van der Waals attractions countered by the stabilizing influence of electrostatic repulsions.<ref> Template:Citation </ref> In 1948, Evert Verwey and Theodor Overbeek independently arrived at the same result.<ref name="Verwey"> Template:Citation </ref> This so-called DLVO theory resolved the failure of the Levine–Dube theory to account for the dependence of colloidal dispersions' stability on the ionic strength of the electrolyte.<ref> Template:Citation </ref>

Derivation

DLVO theory is the combined effect of van der Waals and double layer force. For the derivation, different conditions must be taken into account and different equations can be obtained.<ref name="Elimelech">M. Elimelech, J. Gregory, X. Jia, R. A. Williams, Particle Deposition and Aggregation Measurement: Modelling and Simulation (Boston: 1995).</ref> But some useful assumptions can effectively simplify the process, which are suitable for ordinary conditions. The simplified way to derive it is to add the two parts together.

van der Waals attraction

Template:Main article van der Waals force is actually the total name of dipole-dipole force, dipole-induced dipole force and dispersion forces,<ref name="Jacob">Jacob N. Israelacvili, Intermolecular and Surface Forces (London 2007).</ref> in which dispersion forces are the most important part because they are always present. Assume that the pair potential between two atoms or small molecules is purely attractive and of the form w = −C/rⁿ, where C is a constant for interaction energy, decided by the molecule's property and n = 6 for van der Waals attraction.<ref name="London">London, F. (1937), Trans Faraday Soc, 33, 8–26.</ref> With another assumption of additivity, the net interaction energy between a molecule and planar surface made up of like molecules will be the sum of the interaction energy between the molecule and every molecule in the surface body.<ref name="Jacob" /> So the net interaction energy for a molecule at a distance D away from the surface will therefore be <math display="block">w(D) = -2 \pi \, C \rho _1\, \int_{z=D}^{z= \infty \,} dz \int_{x=0}^{x=\infty \,}\frac{x \, dx}{(z^2+x^2)^3} = \frac{2 \pi C \rho _1}{4} \int_D^\infty \frac{dz}{z^4} = - \frac{ \pi C \rho _1 }{ 6 D^3 }</math>

where

• Template:Math is the interaction energy between the molecule and the surface,
• <math> \rho_1 </math> is the number density of the surface,
• Template:Math is the axis perpendicular to the surface and passesding across the molecule, with Template:Math at the point where the molecule is, and Template:Math at the surface,
• Template:Math is the axis perpendicular to the Template:Math axis, with Template:Math at the intersection.

Then the interaction energy of a large sphere of radius R and a flat surface can be calculated as <math display="block">W(D) = -\frac{2 \pi C \rho _1 \rho _2}{12} \int_{z=0}^{z=2R}\frac {(2R-z)zdz}{(D+z)^3} \approx -\frac{ \pi ^2 C \rho _1 \rho _2 R}{6D}</math>

where

• W(D) is the interaction energy between the sphere and the surface,
• <math>\rho_2</math> is the number density of the sphere.

For convenience, Hamaker constant A is given as <math display="block"> A = \pi^2C\rho_1\rho_2, </math> and the equation becomes <math display="block">W(D) = -\frac{AR}{6D}. </math>

With a similar method and according to Derjaguin approximation,<ref name="Derjaguin">Derjaguin B. V. (1934)Kolloid Zeits 69, 155–164.</ref> the van der Waals interaction energy between particles with different shapes can be calculated, such as energy between

• two spheres: <math>W(D) = -\frac{A}{6D} \frac{R_1 R_2}{(R_1 +R_2 )},</math>
• sphere and surface: <math>W(D) = -\frac{AR}{6D},</math>
• two surfaces: <math>W(D) = -\frac{A}{12 \pi D^2}</math> per unit area.

Double layer force

Template:Main article A surface in a liquid may be charged by dissociation of surface groups (e.g. silanol groups for glass or silica surfaces<ref>Template:Cite journal</ref>) or by adsorption of charged molecules such as polyelectrolyte from the surrounding solution. This results in the development of a wall surface potential which will attract counterions from the surrounding solution and repel co-ions. In equilibrium, the surface charge is balanced by oppositely charged counterions in solution. The region near the surface of enhanced counterion concentration is called the electrical double layer (EDL). The EDL can be approximated by a sub-division into two regions. Ions in the region closest to the charged wall surface are strongly bound to the surface. This immobile layer is called the Stern or Helmholtz layer. The region adjacent to the Stern layer is called the diffuse layer and contains loosely associated ions that are comparatively mobile. The total electrical double layer due to the formation of the counterion layers results in electrostatic screening of the wall charge and minimizes the Gibbs free energy of EDL formation.

The thickness of the diffuse electric double layer is known as the Debye screening length <math>1 / \kappa</math>. At a distance of two Debye screening lengths the electrical potential energy is reduced to 2 percent of the value at the surface wall.

<math display="block">\kappa = \sqrt{\sum_i \frac{\rho_{\infty i} e^2z^2_i}{\epsilon_r \epsilon_0 k_\text{B} T}}</math>

with unit of Template:Math, where

• <math>\rho_{\infty i}</math> is the number density of ion i in the bulk solution,
• Template:Math is the valency of the ion (for example, H⁺ has a valency of +1, and Ca²⁺ has a valency of +2),
• <math>\varepsilon_0</math> is the vacuum permittivity, <math>\epsilon_r</math> is the relative static permittivity,
• Template:Math is the Boltzmann constant.

The repulsive free energy per unit area between two planar surfaces is shown as <math display="block">W = \frac{64k_\text{B} T\rho_{\infty } \gamma ^2}{\kappa}e^{-\kappa D}</math> where

• <math>\gamma</math> is the reduced surface potential, <math>\gamma = \tanh\left(\frac{ze\psi_0}{4k_\text{B}T}\right)</math>,
• <math>\psi_0</math> is the potential on the surface.

The interaction free energy between two spheres of radius R is<ref> Template:Citation</ref> <math display="block">W = \frac{64\pi k_\text{B} TR\rho_{\infty} \gamma ^2}{\kappa ^2}e^{-\kappa D}.</math>

Combining the van der Waals interaction energy and the double layer interaction energy, the interaction between two particles or two surfaces in a liquid can be expressed as <math display="block">W(D) = W(D)_\text{A} + W(D)_\text{R},</math> where W(D)_R is the repulsive interaction energy due to electric repulsion, and W(D)_A is the attractive interaction energy due to van der Waals interaction.

Effect of shear flows

Alessio Zaccone and collaborators investigated the effects of shear-flow on particle aggregation<ref name=" PRE ">Template:Cite journal</ref> which can play an important role in applications e.g. microfluidics, chemical reactors, atmospheric and environmental flows. Their work showed a characteristic lag-time in the shear-induced aggregation of the particles, which decreases exponentially with the shear rate.<ref>Template:Cite journal</ref>

Application

Since the 1940s, the DLVO theory has been used to explain phenomena found in colloidal science, adsorption and many other fields. Due to the more recent popularity of nanoparticle research, DLVO theory has become even more popular because it can be used to explain behavior of both material nanoparticles such as fullerene particles and microorganisms. For example, DLVO theory has been widely applied to assess the degree of particle-particle interactions at controlled chemical conditions. For example, it has been used to investigate the colloidal stability of BaSO₄ (barium sulfate).<ref name="Gallardoetal2000">Template:Cite journal</ref><ref name="Seepmaetal2023a">Template:Cite journal</ref> and particle-particle interactions between magnesite, dolomite, quartz and serpentine.<ref name="Yaoetal2016">Template:Cite journal</ref>

Shortcomings

Additional forces beyond the DLVO construct have been reported to also play a major role in determining colloid stability.<ref name="McBride1997">Template:Cite journal</ref><ref name="Grasso2002">Template:Cite journal</ref> DLVO theory is not effective in describing ordering processes such as the evolution of colloidal crystals in dilute dispersions with low salt concentrations. It also cannot explain the relation between the formation of colloidal crystals and salt concentrations.<ref>N. Ise and I. S. Sogami, Structure Formation in Solution: Ionic Polymers and Colloidal Particles, (Springer, New York, 2005).</ref>