The DLVO theory (named after
Boris Derjaguin
Boris Vladimirovich Derjaguin (or Deryagin; russian: Бори́с Влади́мирович Деря́гин) (9 August 1902 in Moscow – 16 May 1994) was a Soviet and Russian chemist. As a member of the Russian Academy of Sciences, he laid the f ...
and
Lev Landau
Lev Davidovich Landau (russian: Лев Дави́дович Ланда́у; 22 January 1908 – 1 April 1968) was a Soviet- Azerbaijani physicist of Jewish descent who made fundamental contributions to many areas of theoretical physics.
His a ...
,
Evert Verwey and
Theodoor Overbeek
Jan Theodoor Gerard Overbeek (Groningen, January 30, 1911 – February 19, 2007) was a Dutch professor of physical chemistry at the Utrecht University.
Early life and education
Overbeek was born in Groningen. Overbeek's family moved to Rotterd ...
) explains the aggregation 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
counterion
160px, Polystyrene sulfonate, a cation-exchange resin, is typically supplied with as the counterion.">cation-exchange_resin.html" ;"title="Polystyrene sulfonate, a cation-exchange resin">Polystyrene sulfonate, a cation-exchange resin, is typical ...
s.
The electrostatic part of the DLVO interaction is computed in the
mean field approximation in the limit of low
surface potential
Surface charge is a two-dimensional surface with non-zero electric charge. These electric charges are constrained on this 2-D surface, and surface charge density, measured in coulombs per square meter (C•m−2), is used to describe the charge di ...
s - that is when the
potential energy
In physics, potential energy is the energy held by an object because of its position relative to other objects, stresses within itself, its electric charge, or other factors.
Common types of potential energy include the gravitational potentia ...
of an elementary charge on the surface is much smaller than the thermal energy scale,
.
For two spheres of radius
each having a charge
(expressed in units of the elementary charge)
separated by a center-to-center distance
in a fluid of
dielectric constant
The relative permittivity (in older texts, dielectric constant) is the permittivity of a material expressed as a ratio with the electric permittivity of a vacuum. A dielectric is an insulating material, and the dielectric constant of an insulat ...
containing a concentration
of monovalent ions,
the electrostatic potential takes the form of a screened-Coulomb or
Yukawa potential
In particle, atomic and condensed matter physics, a Yukawa potential (also called a screened Coulomb potential) is a potential named after the Japanese physicist Hideki Yukawa. The potential is of the form:
:V_\text(r)= -g^2\frac,
where is a m ...
,
:
where
is the
Bjerrum length
The Bjerrum length (after Danish chemist Niels Bjerrum 1879–1958 ) is the separation at which the electrostatic interaction between two elementary charges is comparable in magnitude to the thermal energy scale, k_\text T, where k_\text is the B ...
,
is the
Debye–Hückel screening length,
which is given by
, and
is the thermal energy scale at absolute temperature
.
Overview
DLVO theory is a theory of colloidal dispersion stability in which
zeta potential
Zeta potential is the electrical potential at the slipping plane. This plane is the interface which separates mobile fluid from fluid that remains attached to the surface.
Zeta potential is a scientific term for electrokinetic potential in coll ...
is used to explain that as two particles approach one another their ionic atmospheres begin to overlap and a repulsion force is developed. In this theory, two forces are considered to impact on colloidal stability:
Van der Waals force
In molecular physics, the van der Waals force is a distance-dependent interaction between atoms or molecules. Unlike ionic or covalent bonds, these attractions do not result from a chemical electronic bond; they are comparatively weak and th ...
s and
electrical double layer
A double layer (DL, also called an electrical double layer, EDL) is a structure that appears on the surface of an object when it is exposed to a fluid. The object might be a solid particle, a gas bubble, a liquid droplet, or a porous body. The D ...
forces.
The total
potential energy
In physics, potential energy is the energy held by an object because of its position relative to other objects, stresses within itself, its electric charge, or other factors.
Common types of potential energy include the gravitational potentia ...
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.
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
In chemistry and physics, activation energy is the minimum amount of energy that must be provided for compounds to result in a chemical reaction. The activation energy (''E''a) of a reaction is measured in joules per mole (J/mol), kilojoules p ...
, and subsequently passes through a shallow minimum, which is referred to as the secondary minimum.
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.
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
In physics, the kinetic energy of an object is the energy that it possesses due to its motion.
It is defined as the work needed to accelerate a body of a given mass from rest to its stated velocity. Having gained this energy during its accele ...
due to their velocity and mass.
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
colloid
A colloid is a mixture in which one substance consisting of microscopically dispersed insoluble particles is suspended throughout another substance. Some definitions specify that the particles must be dispersed in a liquid, while others extend ...
s can be considered to be trapped together by
Van der Waals force
In molecular physics, the van der Waals force is a distance-dependent interaction between atoms or molecules. Unlike ionic or covalent bonds, these attractions do not result from a chemical electronic bond; they are comparatively weak and th ...
s.
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. 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. Particles form weak attractions but are easily redispersed. Thus, the adhesion at secondary minimum can be reversible.
History
In 1923,
Debye and Hückel reported the first successful theory for the distribution of charges in ionic solutions.
The framework of linearized Debye–Hückel theory subsequently was applied to colloidal dispersions by Levine and Dube
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,
Derjaguin and
Landau
Landau ( pfl, Landach), officially Landau in der Pfalz, is an autonomous (''kreisfrei'') town surrounded by the Südliche Weinstraße ("Southern Wine Route") district of southern Rhineland-Palatinate, Germany. It is a university town (since 1990 ...
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.
Seven years later,
Verwey and
Overbeek independently arrived at the same result.
[
.
]
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.
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.
[M. Elimelech, J. Gregory, X. Jia, R. A. Williams, ''Particle Deposition and Aggregation Measurement: Modelling and Simulation'' (Boston: 1995).] 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
van der Waals force is actually the total name of dipole-dipole force, dipole-induced dipole force and dispersion forces,
[Jacob N. Israelacvili, ''Intermolecular and Surface Forces'' (London 2007).] 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
n, where C is a constant for interaction energy, decided by the molecule's property and n = 6 for van der Waals attraction.
[London, F. (1937), ''Trans Faraday Soc'', 33, 8–26.] 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.
So the net interaction energy for a molecule at a distance D away from the surface will therefore be
:
where
: ''w''(''r'') is the interaction energy between the molecule and the surface,
:
is the number density of the surface,
: ''z'' is the axis perpendicular to the surface and passesding across the molecule, with ''z'' = ''D'' at the point where the molecule is, and ''z'' = 0 at the surface,
: ''x'' is the axis perpendicular to the ''z'' axis, with ''x'' = 0 at the intersection.
Then the interaction energy of a large sphere of radius ''R'' and a flat surface can be calculated as
:
where
: ''W''(''D'') is the interaction energy between the sphere and the surface,
:
is the number density of the sphere.
For convenience,
Hamaker constant
The Hamaker constant ''A'' can be defined for a van der Waals (vdW) body–body interaction:
:A=\pi^2C\rho_1\rho_2,
where \rho_1 and \rho_2 are the number densities of the two interacting kinds of particles, and ''C'' is the London coefficient in ...
''A'' is given as
:
and the equation becomes
:
With a similar method and according to
Derjaguin approximation
The Derjaguin approximation (or sometimes also referred to as the proximity approximation), named after the Russian scientist Boris Derjaguin, expresses the force profile acting between finite size bodies in terms of the force profile between two ...
,
[Derjaguin B. V. (1934)''Kolloid Zeits'' 69, 155–164.] the van der Waals interaction energy between particles with different shapes can be calculated, such as energy between
: two spheres:
: sphere and surface:
: two surfaces:
per unit area.
Double layer force
A surface in a liquid may be charged by dissociation of surface groups (e.g. silanol groups for glass or silica surfaces) or by adsorption of charged molecules such as
polyelectrolyte
Polyelectrolytes are polymers whose repeating units bear an electrolyte group. Ion#Anions and cations, Polycations and polyanions are polyelectrolytes. These groups dissociation (chemistry), dissociate in aqueous solutions (water), making the pol ...
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
In thermodynamics, the Gibbs free energy (or Gibbs energy; symbol G) is a thermodynamic potential that can be used to calculate the maximum amount of work that may be performed by a thermodynamically closed system at constant temperature and pr ...
of EDL formation.
The thickness of the diffuse electric double layer is known as the
Debye screening length
A double layer (DL, also called an electrical double layer, EDL) is a structure that appears on the surface of an object when it is exposed to a fluid. The object might be a solid particle, a gas bubble, a liquid droplet, or a porous body. The D ...
. At a distance of two Debye screening lengths the electrical potential energy is reduced to 2 percent of the value at the surface wall.
:
with unit of ''m''
−1, where
:
is the number density of ion i in the bulk solution,
: ''z'' is the valency of the ion (for example, H
+ has a valency of +1, and Ca
2+ has a valency of +2),
:
is the
vacuum permittivity
Vacuum permittivity, commonly denoted (pronounced "epsilon nought" or "epsilon zero"), is the value of the absolute dielectric permittivity of classical vacuum. It may also be referred to as the permittivity of free space, the electric consta ...
,
is the
relative static permittivity
The relative permittivity (in older texts, dielectric constant) is the permittivity of a material expressed as a ratio with the vacuum permittivity, electric permittivity of a vacuum. A dielectric is an insulating material, and the dielectric co ...
,
: ''k''
B is the
Boltzmann constant
The Boltzmann constant ( or ) is the proportionality factor that relates the average relative kinetic energy of particles in a gas with the thermodynamic temperature of the gas. It occurs in the definitions of the kelvin and the gas constant, ...
.
The repulsive free energy per unit area between two planar surfaces is shown as
:
where
:
is the reduced surface potential,
,
:
is the potential on the surface.
The interaction free energy between two spheres of radius ''R'' is
:
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
:
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
The DLVO theory of colloidal stability has been extended to take into account the effect of
shear flow The term shear flow is used in solid mechanics as well as in fluid dynamics. The expression ''shear flow'' is used to indicate:
* a shear stress over a distance in a thin-walled structure (in solid mechanics);Higdon, Ohlsen, Stiles and Weese (1960) ...
in
fluid dynamic
In physics and engineering, fluid dynamics is a subdiscipline of fluid mechanics that describes the flow of fluids— liquids and gases. It has several subdisciplines, including ''aerodynamics'' (the study of air and other gases in motion) an ...
systems, which is relevant for many applications e.g.
microfluidics
Microfluidics refers to the behavior, precise control, and manipulation of fluids that are geometrically constrained to a small scale (typically sub-millimeter) at which surface forces dominate volumetric forces. It is a multidisciplinary field tha ...
,
chemical reactor
A chemical reactor is an enclosed volume in which a chemical reaction takes place. In chemical engineering, it is generally understood to be a process vessel used to carry out a chemical reaction, which is one of the classic unit operations in chem ...
s, atmospheric and
environmental flow
Environmental flows describe the quantity, timing, and quality of water flows required to sustain freshwater and estuarine ecosystems and the human livelihoods and well being that depend on these ecosystems. In the Indian context river flows requir ...
s, in the work of Alessio Zaccone and collaborators.
In this extended DLVO theory for sheared systems, the DLVO energy barrier for aggregation gets reduced by a negative contribution which is proportional to the
Péclet number
In continuum mechanics, the Péclet number (, after Jean Claude Eugène Péclet) is a class of dimensionless numbers relevant in the study of transport phenomena in a continuum. It is defined to be the ratio of the rate of advection of a physical ...
of the particles, i.e. proportional to the
shear rate
In physics, shear rate is the rate at which a progressive shearing deformation is applied to some material.
Simple shear
The shear rate for a fluid flowing between two parallel plates, one moving at a constant speed and the other one stationary ...
, to the
viscosity
The viscosity of a fluid is a measure of its resistance to deformation at a given rate. For liquids, it corresponds to the informal concept of "thickness": for example, syrup has a higher viscosity than water.
Viscosity quantifies the inte ...
of the medium, and to the cube of the colloidal particle size, while the proportionality coefficient depends on the flow geometry.
This result is obtained from an approximate solution to the governing
Smoluchowski convection-diffusion equation obtained via the
method of matched asymptotic expansions In mathematics, the method of matched asymptotic expansions is a common approach to finding an accurate approximation to the solution to an equation, or system of equations. It is particularly used when solving singularly perturbed differential eq ...
.
The theory explains a characteristic lag-time in the shear-induced aggregation of the particles, which decreases exponentially with the
shear rate
In physics, shear rate is the rate at which a progressive shearing deformation is applied to some material.
Simple shear
The shear rate for a fluid flowing between two parallel plates, one moving at a constant speed and the other one stationary ...
. It also explains the subsequent runaway (autocatalytic) regime of aggregation kinetics after the lag-time, as well as the characteristic bimodal cluster size distribution of aggregates found typically in shear-induced aggregating and self-assembling systems.
Furthermore, the theory has been verified across widely different flow conditions, in terms of different particles and microfluidic systems, and viscoelastic properties of the liquid phase.
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
A fullerene is an allotrope of carbon whose molecule consists of carbon atoms connected by single and double bonds so as to form a closed or partially closed mesh, with fused rings of five to seven atoms. The molecule may be a hollow sphere, ...
particles and
microorganisms
A microorganism, or microbe,, ''mikros'', "small") and ''organism'' from the el, ὀργανισμός, ''organismós'', "organism"). It is usually written as a single word but is sometimes hyphenated (''micro-organism''), especially in olde ...
.
Shortcomings
Additional forces beyond the DLVO construct have been reported to also play a major role in determining colloid stability.
DLVO theory is not effective in describing ordering processes such as the evolution of
colloidal crystal A colloidal crystal is an ordered array of colloid particles and fine grained materials analogous to a standard crystal whose repeating subunits are atoms or molecules. A natural example of this phenomenon can be found in the gem opal, where sphere ...
s in dilute dispersions with low salt concentrations. It also cannot explain the relation between the formation of colloidal crystals and salt concentrations.
[
N. Ise and I. S. Sogami,
''Structure Formation in Solution: Ionic Polymers and Colloidal Particles,''
(Springer, New York, 2005).
]
References
{{DEFAULTSORT:Dlvo Theory
Physical chemistry
Colloidal chemistry