⚗️ Colloidal Stability — DLVO Theory

DLVO (Derjaguin-Landau-Verwey-Overbeek) theory: total interaction energy V_T = V_EDL + V_vdW. Electrostatic repulsion vs van der Waals attraction determines colloidal stability.

ChemistryInteractive
DLVO interaction energy V_T vs separation h · Green = repulsion (EDL) · Red = attraction (vdW)

How it Works

DLVO theory calculates the total interaction energy between two spherical colloidal particles as a function of surface-to-surface separation h. The two contributions are: (1) electrostatic double-layer repulsion V_EDL which decays exponentially with the Debye length κ⁻¹, and (2) van der Waals attraction V_vdW which decays as 1/h at close range.

A colloid is stable when the energy barrier (maximum of V_T) is much greater than thermal energy kT ≈ 4.1 × 10⁻²¹ J at room temperature. Increasing salt concentration compresses the double layer (shorter Debye length), reducing the barrier until coagulation occurs.

V_EDL = 64πε₀εr·a·(kT/ze)²·tanh²(zeψ₀/4kT)·exp(−κh)
V_vdW = −A·a / (12h) [sphere-sphere, Derjaguin approx]
κ⁻¹ = √(ε₀εrkT / 2NAe²I) [Debye length]
V_T = V_EDL + V_vdW

Frequently Asked Questions

What is DLVO theory?

DLVO theory (Derjaguin-Landau-Verwey-Overbeek) describes the forces between charged particles in solution. The total interaction energy V_T = V_EDL + V_vdW is the sum of electrostatic double-layer repulsion and van der Waals attraction.

What is the electrical double layer?

The electrical double layer consists of a charged particle surface and a diffuse cloud of counterions in solution. It produces a repulsive force between particles. The characteristic decay length is the Debye length κ⁻¹, which decreases with increasing salt concentration.

What is the Debye length?

The Debye length κ⁻¹ = √(ε₀εrkT / 2e²NAI) is the characteristic screening length of the electrostatic interaction. At high salt concentration (high ionic strength I), the Debye length shrinks, reducing repulsion and promoting aggregation.

What is the Hamaker constant?

The Hamaker constant A characterizes the strength of van der Waals attraction between particles across a medium. For silica across water A ≈ 0.83×10⁻²⁰ J; for gold A ≈ 40×10⁻²⁰ J. Higher A means stronger attraction and less stable colloid.

What is zeta potential and why does it matter?

Zeta potential is the electric potential at the slipping plane of a particle moving in solution. Colloids with |ζ| > 30 mV are typically stable. As |ζ| decreases (e.g., by adding salt), aggregation becomes likely.

What is the critical coagulation concentration?

The critical coagulation concentration (CCC) is the salt concentration at which the energy barrier disappears and rapid coagulation occurs. The Schulze-Hardy rule states CCC ∝ 1/z⁶, where z is the valence of the counterion.

What is the primary minimum in DLVO energy?

The primary minimum is the deep attractive well at very close approach (less than 1 nm) where van der Waals attraction dominates. Particles that overcome the energy barrier and reach the primary minimum form irreversible aggregates.

What is the secondary minimum?

The secondary minimum is a shallow attractive well at larger separations (several nm) that can trap particles in reversible, loose aggregates called flocs. Unlike primary minimum aggregation, flocculation in the secondary minimum is reversible by gentle agitation.

Beyond DLVO: what else affects colloidal stability?

Beyond DLVO forces, steric stabilization (polymer brushes), depletion forces (non-adsorbing polymers), hydration forces (water structuring near hydrophilic surfaces), and hydrophobic attraction all affect colloidal stability.

How is colloidal stability measured experimentally?

Dynamic light scattering (DLS) tracks particle size over time to detect aggregation. Electrophoresis measures zeta potential. Turbidity measurements monitor sedimentation. Critical coagulation concentration is found by titrating salt and measuring aggregation onset.

About this simulation

This simulator numerically computes the full DLVO interaction-energy curve between two spherical colloidal particles across separations from 0.1 to 30 nm, summing electrostatic double-layer repulsion (V_EDL, decaying with the Debye length) and van der Waals attraction (V_vdW, decaying as 1/h). It scans the curve to find the energy barrier V_max in units of kT and classifies the colloid as stable, metastable, or actively coagulating.

🔬 What it shows

A live plot of V_T (gold), V_EDL (dashed green), and V_vdW (dashed red) against surface separation h. The barrier height V_max is marked with a dot and label, and the stats panel reports the Debye length, barrier height, primary-minimum depth, and an overall stability verdict.

🎮 How to use

Adjust particle radius a, surface potential ψ₀, ionic strength I, and the Hamaker constant A directly, or pick a material preset (silica, TiO₂, gold, polystyrene) to load textbook parameter combinations. Push ionic strength up to see the barrier collapse toward coagulation, or push |ψ₀| up to rebuild it.

💡 Did you know?

Gold's Hamaker constant (~40×10⁻²⁰ J) is roughly 50 times larger than silica's (~0.83×10⁻²⁰ J), which is exactly why gold nanoparticle suspensions are so much harder to keep stable without a strong stabilizing surface charge or coating.

Frequently asked questions

What does the V_max barrier value actually tell you about stability?

V_max is the height of the energy hump particles must climb over (in units of thermal energy kT) before they can fall into the deep primary-energy minimum and aggregate irreversibly. The simulator labels a colloid "Stable" above 25 kT, "Metastable" between 10-25 kT, and "Unstable / coagulating" below 10 kT, since particles constantly bump each other with roughly kT of thermal energy.

Why does raising the ionic strength slider destabilize the colloid?

Increasing ionic strength I shrinks the Debye length κ⁻¹ (it's computed directly from I in the code), which compresses the electrostatic double layer and makes V_EDL decay much faster with distance. That shortens the effective range of repulsion, shrinking or eliminating the energy barrier and letting van der Waals attraction win.

What's the difference between the "primary minimum" reading and the barrier?

The primary minimum is read directly from the first computed separation point (h = 0.1 nm), representing the deep attractive well at near-contact where van der Waals attraction dominates completely. The barrier V_max instead marks the highest point of V_T further out, which is what actually protects particles from reaching that deep well.

Why do different material presets need very different surface potentials and Hamaker constants?

Each preset (silica, TiO₂, gold, polystyrene) encodes real literature values because different materials genuinely differ in both their surface chemistry (which sets ψ₀) and their polarizability across water (which sets the Hamaker constant A). Gold's large A and modest |ψ₀| makes it intrinsically harder to stabilize than the more chemically stable, lower-A silica system.

Does the simulator account for anything beyond classic DLVO forces?

No — this simulation implements the two classic DLVO terms (V_EDL and V_vdW) exactly as derived from the Derjaguin approximation for spheres, with no steric, hydration, or depletion contributions. Real formulated colloids often add polymer-brush steric stabilization on top of these electrostatic and van der Waals forces to reach much higher practical stability.