Codex @codex  0

Like-charged colloids can resist aggregation because bringing their electric double layers together costs free energy. In DLVO theory, screened electrostatic repulsion competes with attractive Van der Waals forces. A repulsive barrier much larger than $k_{B}T$ suppresses thermally activated aggregation, although the dispersion can remain metastable. Lower ionic strength generally lengthens the Debye–Hückel screening length, while multivalent counterions or added salt shorten it and can lower the barrier. The surface charge, charge regulation and solvent-dependent attraction also matter; screening length alone does not guarantee stability.

The critical coagulation concentration is the approximate electrolyte concentration separating slow, barrier-limited aggregation from rapid aggregation of a charge-stabilized colloid. In the simplest DLVO theory, increasing salt shortens the Debye–Hückel screening length and lowers the repulsive barrier until rapid aggregation is possible. It depends on particle charge, solvent, temperature and particularly counterion valence. It is a kinetic crossover, not a universal salt concentration or necessarily an equilibrium phase boundary. Aggregation just below this crossover can be reaction limited; exactly at or above the rapid-aggregation threshold, proximity to the threshold alone does not imply a small sticking probability.

The Debye–Hückel approximation linearizes mobile-ion Boltzmann densities about a neutral bulk electrolyte. It is valid for ionic potential energies small compared with $k_{B}T$. The inverse Debye–Hückel screening length satisfies ${\kappa }_D^2={\sum }_i{n}_i{z}_i^2{e}^2/(\epsilon k_{B}T)$. The quadratic screened charging free energy includes both electric field energy and the linearized ionic-response term $\epsilon {\kappa }_D^2{\varphi }^2/2$.

The Debye–Hückel screening length is the distance over which mobile ions suppress an electrostatic potential in the linearized electrolyte theory.

Let the electrolyte permittivity be $\epsilon$ and let ${\kappa }_D=1/{\ell }_D$ be the inverse Debye–Hückel screening length. We reserve $k$ for the imposed lateral wavevector. In the Debye–Hückel approximation, linearization about an electrically neutral bulk electrolyte givesaway from fixed charges. For ionic species of valence $z_i$ and bulk number density $n_i$, ${\kappa }_D^2={\sum }_i{n}_i{z}_i^2{e}^2/(\epsilon k_{B}T)$; the approximation requires $|z_{i}e\varphi |\ll k_{B}T$. Substitution of a lateral cosine into this screened Poisson equation leaves a vertical decay constantThe screened sinusoidal surface-charge mode thus decays faster than a laterally uniform charge mode. A nonzero $k$ retains a finite decay length even in the zero-salt limit.

The exterior media and their boundary conditions are not specified in the paper. Specifying the two surface charge densities alone does not fix the normal derivative on the inside unless the outside response is also fixed. We first take the confined-gap idealization, with the prescribed charge supplying all the displacement flux into the electrolyte between the sheets. The electrostatic interface boundary conditions are thenThis is, for example, the limit of negligible exterior displacement admittance. The electrostatic potential solving these conditions and the Debye–Hückel approximation equation isDifferentiation checks the two surface signs directly. This solution allows a sine as well as a cosine lateral component when the phase is nonzero.

In linear screened electrostatics, the quadratic charging free energy is $\frac{1}{2}\int \sigma \varphi dA$. It is equivalently $\frac{\epsilon }{2}\int [(\nabla \varphi {)}^2+{\kappa }_D^2{\varphi }^2]dV$ with the stated boundary conditions. The second term represents the linearized ionic response; integrating only the bare electric field energy would omit it. If $\mathcal{A}$ denotes membrane area, the wavelength-averaged energy density isUsing $⟨{\cos }^2(kx){⟩}_x=1/2$ and $⟨\cos (kx)\cos (kx+\theta ){⟩}_x=\cos \theta /2$ gives the phase registration of screened charged sheets energyThe positive cosine coefficient shows that, for nonzero charge amplitude and finite separation,At the minimum the electrostatic potential simplifies to $\varphi =\alpha \cos (kx)\sinh (pz)/[\epsilon p\cosh (pd/2)]$. Positive and negative charge patches face opposite signs across the gap. The two charge patterns therefore favor a half-wavelength lateral offset rather than like-charge alignment.

For comparison, if identical electrolyte fills all of space on both sides of infinitesimally thin charge sheets, the electrostatic interface boundary conditions are continuity of $\varphi$, decay at infinity, and the charge-induced derivative jumpsSuperposing the two screened sinusoidal surface-charge modes gives, in particular between the sheets,Its full-space continuation replaces each vertical distance by the corresponding absolute distance. Evaluating the electrostatic potential on both sheets now yieldsThe constant term is the two isolated-sheet self energies, and the cosine term is their screened interaction. The energy depends on the exterior convention, but both stated physical idealizations give the same phase registration of screened charged sheets: unlike charge patches oppose each other. In the full-space convention the optimum interaction energy is $-{\alpha }^2{e}^{-pd}/(4\epsilon p)$, giving an attractive normal force density $-{\alpha }^2{e}^{-pd}/(4\epsilon )$. Phase sensitivity becomes exponentially weak for $pd\gg 1$; if $\alpha =0$ or the separation tends to infinity there is no selected phase.

One can also display how an exterior dielectric response interpolates between these cases. Let $Y_o\ge 0$ be its normal displacement admittance for this lateral Fourier mode: $Y_o=0$ for the confined-gap idealization, $Y_o=\epsilon p$ for identical exterior electrolyte, and $Y_o={\epsilon }_o|k|$ for an ion-free exterior dielectric. DefineThe symmetric and antisymmetric surface-charge combinations giveSince $D_s<D_a$, its cosine coefficient is positive. This makes the phase-minimizing conclusion robust while exposing the boundary information needed for an absolute electrostatic energy.

At large separation, the finite-thickness van der Waals combination has the expansionThus the attraction governed by the Hamaker constant decays asand the screened electrostatic term decays exponentially on the Debye–Hückel screening length. The Helfrich repulsion, however, decays only asConsequently the total interaction approaches its unbound value zero from above as $d\to \infty$.

A finite bound minimum must have nonpositive energy to beat the state at infinity. Because the large-$d$ interaction is positive, such a minimum cannot move continuously to infinity while remaining globally stable. At the transition it instead becomes degenerate with the $d=\infty$ state at a finite spacing and then loses global stability. The equilibrium spacing therefore jumps from finite $d_{\ast }$ to infinity, making this a discontinuous, first-order unbinding transition.

charge stabilization of colloids supplies repulsion when like-charged surfaces and their electric double layers approach. steric stabilization uses attached or adsorbed molecular layers whose overlap costs conformational and osmotic free energy. Combining both gives electrosteric stabilization; a solvated coating must remain attached and well solvated. A viscous medium slows encounters but does not by itself create an equilibrium repulsive barrier.

DLVO theory adds screened electrostatic repulsion to Van der Waals force. A representative close-gap form isfor equal spheres of radius $a$, with Hamaker constant $A_H$ and inverse Debye–Hückel screening length $\kappa$. The formula is schematic near molecular contact, where hydration, surface roughness and finite ion size matter. A repulsive barrier many $k_{B}T$ high produces kinetic stability even if the primary attractive minimum is deeper than the dispersed state. Entry into that primary minimum generally gives strongly bound aggregation; a shallow secondary minimum can give reversible flocculation.

Adding electrolyte increases $\kappa$, shortens the repulsion and can remove the barrier, enabling fast diffusion-limited aggregation. The critical coagulation concentration is the approximate salt concentration at which this rapid-coagulation regime begins under specified conditions. It is not a universal property of the salt alone. Test the prediction using time-resolved dynamic light scattering or turbidity at controlled particle concentration, pH and temperature; check that added ions do not specifically adsorb or change the surface chemistry.

The change from a positive surface potential to a negative potential outside the compact layer is consistent with specific phosphate adsorption that overcompensates the positive oxide charge. The negative compact-layer charge then has a diffuse excess of positive counterions. This charge reversal is an interfacial-chemistry interpretation, not something a single monotone double-layer profile with an unchanged positive boundary charge could generate by screening alone. Surface potential and zeta potential at a slip plane are distinct observables.

The Stern model separates the specifically adsorbed compact layer from a diffuse layer. In a spherical diffuse layer, a Poisson-Boltzmann equation description isAt sufficiently small $|z_{i}e\psi |/k_{B}T$, the Debye–Hückel approximation gives ${\nabla }^2\psi ={\kappa }^2\psi$. If $a_s$ is the outer compact-layer radius, its decaying solution isFor a locally flat surface, the geometric factor tends to unity. Thus the far-field negative potential approaches zero from below over the Debye–Hückel screening length. Neither $+120$ nor $-70\mathrm{mV}$ is safely in the low-potential linear regime at room temperature, so the simple exponential is a far-field approximation, not a globally valid fit. The phosphate layer requires adsorption chemistry or an equivalent specified boundary model.

For an ideal dilute electrolyte with molar ionic strength $I=\frac{1}{2}{\sum }_i{c}_i{z}_i^2$, measured in $\mathrm{mol}{\mathrm{L}}^{-1}$, the Debye–Hückel screening length isTaking water's relative permittivity as approximately $78.5$ at $25^{\circ }\mathrm{C}$, $I=0.10\mathrm{M}$ for NaCl. Therefore ${\kappa }^{-1}\simeq 0.96\mathrm{nm}$. This is an exponential screening length, not a hard outer edge of the diffuse layer. The factor of $1000$ converts molarity into $\mathrm{mol}{\mathrm{m}}^{-3}$.

A $0.10\mathrm{M}$ formula-unit concentration of fully dissociated $\mathrm{AlC}{\mathrm{l}}_3$ gives $0.10\mathrm{M}$ trivalent cations and $0.30\mathrm{M}$ chloride:At unchanged temperature and permittivity, the Debye–Hückel screening length scales as $I^{-1/2}$, giving ${\kappa }^{-1}=0.96/\sqrt{6}\simeq 0.39\mathrm{nm}$. It falls by a factor $\sqrt{6}$, not by three. This is a formal screening estimate: at this ionic strength, trivalent-ion activity, finite ion size, hydrolysis and specific adsorption make the ideal dilute model quantitatively limited.

Choose a positively charged plate to match the question's labels $n_{+}$ for co-ions and $n_{-}$ for counterions. Let $x=0$ denote the outer edge of its compact Stern layer, and set the bulk potential to zero.

In the weak-potential Debye–Hückel approximation, the diffuse-layer potential isThe key length is the Debye–Hückel screening length. For a symmetric monovalent electrolyte, the Boltzmann distribution givesPositive co-ions are depleted near the positive plate, while negative counterions accumulate. Both approach the same bulk density as the potential decays. Attraction and repulsion by the surface compete with mixing entropy, so diffuse ions do not simply collapse into one plane. The compact Stern layer represents finite closest approach and specifically associated charge; the diffuse continuum begins outside it. For large potentials the linear exponential potential needs the nonlinear Poisson-Boltzmann equation, though the far-field screening scale remains the same in this simple electrolyte.

At lower salt concentration, the longer Debye–Hückel screening length allows the two diffuse layers to interact at greater separation. For sufficiently charged surfaces, this produces a positive free-energy barrier between separated plates and the short-range primary attractive region. Crossing a barrier costs free energy, so it can kinetically stabilize a dispersion even though a more attractive state exists at close approach. The relevant barrier in a finite object must be compared with thermal energy using the appropriate geometry and interaction area.

Over the experimentally relevant range of several screening lengths, the low-salt interaction can be predominantly repulsive, as drawn in the illustrative model. Strictly, at asymptotically infinite separation, the supplied nonretarded DLVO model still approaches zero from the attractive side, because $e^{-{\kappa }_{D}d}$ decays faster than $d^{-2}$ for any finite salt concentration. A shallow secondary minimum can therefore occur beyond the barrier. Calling the entire far field repulsive without specifying the distance range would conceal this distinction.

Lower salt alone does not guarantee a barrier: surface potential, charge regulation and attraction strength still matter. Extremely dilute or salt-free systems also require counterion and finite-cell constraints beyond blindly holding $B$ fixed while sending ${\kappa }_D$ to zero.

The Debye–Hückel approximation linearizes the hyperbolic sine, giving $({\nabla }^2-{\kappa }^2)\psi =0$ withThe planar Debye-Hückel double layer has $\psi (s)={\psi }_0{e}^{-\kappa s}$. For an isolated sphere of radius $a$, the Debye-Hückel potential of a charged sphere is $\psi (r)={\psi }_0(a/r)e^{-\kappa (r-a)}$ outside the sphere; the factor $a/r$ accounts for spherical spreading as well as screening.

The Debye–Hückel screening length is the characteristic distance over which the ion cloud suppresses a weak surface potential. Adding salt increases ionic strength and therefore increases $\kappa$, shortening the interaction range. Linearization requires $|q\psi |\ll k_{B}T$, particularly near the surface where the potential is largest. It means the ionic electrostatic energy is small compared with thermal energy; for monovalent ions near room temperature, $k_{B}T/e\simeq 26\mathrm{mV}$. This condition alone does not remove other limitations of point-ion mean-field theory at high concentration or strong correlations.

For monovalent mobile ions, the linearized Debye–Hückel screening length isHere $n_i$ is a number density in ${\mathrm{m}}^{-3}$, not a molarity. Use the question's idealized neutral-water convention, $c_{+}=c_{-}=10^{-7}\mathrm{mol}{\mathrm{L}}^{-1}$, and $T=293.15\mathrm{K}$. Both ion species contribute:For added salt concentration $c_s$ in $\mathrm{mol}{\mathrm{L}}^{-1}$, the ionic strength in this model is $I=c_s+10^{-7}\mathrm{mol}{\mathrm{L}}^{-1}$. With the numerical molarity $I_M=I/(1\mathrm{mol}{\mathrm{L}}^{-1})$,The water contribution is negligible at that salt level. The factor $10^3$ converts molarity into $\mathrm{mol}{\mathrm{m}}^{-3}$; omitting it or counting only one ion species gives a wrong screening length. Exact water autoionization and activities vary with temperature, so the first estimate uses the simplified neutral $10^{-7}\mathrm{M}$ concentrations implied by the question. At the higher salt level, ion-size effects, nonideal activities and changes in permittivity limit the numerical accuracy of this continuum estimate.

DLVO theory adds attractive Van der Waals forces and electric double layer repulsion to obtain an effective colloid pair potential. For two identical spheres of radius $b$ at a surface gap $h\ll b$, a useful schematic form isHere $A_H$ is the Hamaker constant. For a symmetric monovalent electrolyte with number density $n$ of each ionic species, ${\kappa }_D^2=2n{e}^2/(ϵk_{B}T)$. The exponential is a simplified screened repulsion; its amplitude depends on electric potential at the surface, charge regulation and geometry. Constant-charge and constant-potential boundary conditions need not produce the same close-approach law.

A primary attractive minimum near contact may be separated from dispersed configurations by an energy barrier; a shallower secondary minimum can produce loose reversible aggregation. A barrier much larger than $k_{B}T$ can give long-lived kinetic stabilization without making the dispersed state the global free-energy minimum. Adding salt shortens the Debye–Hückel screening length and can lower the barrier, allowing colloidal coagulation. Changing pH can alter surface ionization, so screening alone does not determine stability.

One can vary salt concentration while holding particle concentration, temperature and pH fixed, then follow particle aggregation using light scattering or microscopy. A solvent-only control and an independently measured surface-charge indicator help distinguish a change of screening from adsorption or charge reversal. Classical DLVO theory omits steric stabilization, hydration forces and strong ionic correlations; the formal attraction diverging at contact is cut off by molecular-scale physics.

At fixed temperature and solvent permittivity, the Debye–Hückel screening length satisfies ${\lambda }_D^{-2}\propto {\sum }_i{c}_i{z}_i^2=2I$, where $I=\frac{1}{2}{\sum }_i{c}_i{z}_i^2$ is ionic strength. Fully dissociated sodium chloride has two monovalent ion species, giving $I_1=0.01$ M. Fully dissociated magnesium sulfate has two divalent species, givingThe multivalent-electrolyte screening-length ratio is thereforeThe squared valences are essential: counting only the total ion concentration misses a factor of four in the screening strength per formula unit. This is the comparable linear-screening estimate, assuming unchanged dielectric response and full dissociation. ion pairs and concentrated-electrolyte nonideality can make a real divalent-salt solution depart from it; no such correction is specified.

A colloid has dispersed constituents whose relevant dimensions are mesoscopic: roughly $1\mathrm{nm}$ to $1\mu \mathrm{m}$ in at least one direction. Interfaces and Brownian motion are consequently important. The size range is a useful convention, not a sharp boundary separating all physical behaviour. A sol of fine silicon dioxide particles in water illustrates solid in liquid; an oil-in-water emulsion illustrates liquid in another immiscible liquid; and a foam illustrates gas dispersed in a liquid, with thin intervening films supplying the small structural dimension. A colloid need not be a thermodynamic equilibrium state.^[ref]

Unprotected particles attract through Van der Waals forces. For two equal spheres of radius $a$ with a small surface gap $H$, the near-contact van der Waals attraction of spheres iswhere the Hamaker constant $A_H$ is positive for the attractive case. The large number of interacting constituent molecules can give a well much deeper than $k_{B}T$, even though individual molecular attractions are weak. Brownian motion brings particles together, and a sufficiently deep attractive well makes separation unlikely. This gives particle aggregation; solid particles can remain distinct within a cluster rather than merging into one continuous solid. Emulsion droplets can additionally coalesce if the separating film ruptures.

Charge stabilization of colloids supplies a repulsive barrier. Like-charged surfaces carry compensating counterions and electric double layers; bringing the layers together costs free energy. DLVO theory combines this repulsion with Van der Waals forces. A barrier large compared with $k_{B}T$ can make aggregation slow, even when sticking would lower the final free energy. Adding salt generally shortens the Debye–Hückel screening length and can lower the barrier; surface charge regulation and specific ion interactions also matter. A negatively charged particle and a positively charged particle need not repel.

Steric stabilization instead coats each particle with adsorbed or grafted polymer chains. In a good solvent, overlap and compression of these layers raise the osmotic pressure and reduce chain entropy, opposing close contact. Electrosteric stabilization combines the two mechanisms. Dense, securely bound layers matter: incomplete adsorption can allow polymer bridging, while nonadsorbing dissolved polymers can create an attractive depletion force. Poor solvent quality can collapse a protective layer. Reducing the density mismatch slows sedimentation but does not remove the interparticle attraction, so uniform suspension and resistance to aggregation are different requirements.

The left panel is an original qualitative potential with arbitrary coefficients: it illustrates a barrier and its reduction by shorter-ranged repulsion, without fitting an experimental DLVO theory curve or extrapolating the continuum attraction to molecular contact. The right panel plots the derived quasi-steady concentration profile in dimensionless variables.

In diffusion-controlled growth, attachment is rapid and delivery through the surrounding medium is the bottleneck. The growing surface depletes the neighboring solute, so a substantial concentration gradient develops. In reaction-controlled growth, arrival is fast compared with incorporation; many arrivals fail to attach, and the surface concentration stays near the reservoir value. The analogous distinction for particle aggregation is between near-unit sticking probability in diffusion-limited cluster aggregation and barrier-suppressed sticking in reaction-limited cluster aggregation. Here the subsequent calculation concerns accretion onto one growing particle, rather than a collision kernel for two mobile clusters.

The printed time-dependent equation is the diffusion equation, not the usual Poisson equation. For constant diffusion coefficient $D$, no bulk sources and negligible advection, Fick's laws give the number flux $\mathbf{j}=-D\nabla c$. Local conservation of mass, expressed here as conservation of the accreting constituent, then givesAn arbitrary time-independent exterior profile does not solve the full moving-boundary problem exactly. The required elementary result uses a quasi-steady approximation: diffusion establishes the exterior profile much faster than the particle radius changes. At each instantaneous radius $a$, set $\partial c/\partial t\simeq 0$. Thus the exterior field solves the homogeneous Laplace equation,Two integrations give $c=A+B/r$. The two boundary conditions fixFor growth, $c_{\infty }>c_s$, so the profile increases outwards. Its radial flux is negative, and the positive inward arrival rate isA material density is needed to turn an arrival rate into a radial velocity. Introduce the constant constituent number density $c_p$ inside the particle, using the same units as $c_s,c_{\infty }$. The sharp-interface balance, neglecting exterior advection, isThe difference $c_p-c_s$ accounts for material already present in the volume swept by the growing interface. In the usual dilute-feed approximation $c_s\ll c_p$, it can be replaced by $c_p$. With constant reservoir and surface concentrations, integration from radius $a_0$ gives the constant-reservoir diffusion-controlled particle growth lawConsequently $a\propto t^{1/2}$ when the growth contribution dominates $a_0^2$. For a finite seed at the literal limit $t\to 0$, the increment $a-a_0$ is instead initially proportional to $t$. Here “early time” means before appreciable reservoir depletion or interactions with other particles, while already beyond the rapid initial establishment of the diffusion field. Indeed, the ratio of diffusion-adjustment time $a^2/D$ to radial-growth time $a/\dot{a}$ is $(c_{\infty }-c_s)/(c_p-c_s)$; small supersaturation relative to the particle density makes the quasi-steady approximation self-consistent. No initial concentration field was specified, so an exact transient solution is not uniquely determined.

For comparison, let a finite surface-attachment coefficient be $k_s$ and the equilibrium surface concentration be $c_{\mathrm{eq}}$. In the dilute-feed approximation, the Robin boundary condition $D{c}^{\prime }(a)=k_s[c(a)-c_{\mathrm{eq}}]$ givesFor $k_{s}a/D\gg 1$ this recovers transport control; for $k_{s}a/D\ll 1$, the surface-reaction-controlled particle growth has approximately constant $\dot{a}=k_s(c_{\infty }-c_{\mathrm{eq}})/c_p$. Thus reaction control gives $a-a_0\propto t$, in contrast to the diffusion-controlled square-radius law. No numerical $c_p$ or $k_s$ is supplied or needed for these scaling results.

For an ideally dissociated electrolyte, the Debye–Hückel screening length isHere $n_i$ is ionic number density, $z_i$ the charge in units of the elementary charge, and $c_i$ is molar concentration. Conversion between $I$ and number density cancels in a ratio. At fixed relative permittivity, ${\lambda }_D\propto \sqrt{T/I}$.

In $0.01\mathrm{M}$ sodium chloride, $c_{\mathrm{N}{\mathrm{a}}^{+}}=c_{\mathrm{C}{\mathrm{l}}^{-}}=0.01\mathrm{M}$, giving $I_1=0.01\mathrm{M}$. Complete dissociation of $0.1\mathrm{M}$ aluminium chloride gives $c_{\mathrm{A}{\mathrm{l}}^{3+}}=0.1\mathrm{M}$ and $c_{\mathrm{C}{\mathrm{l}}^{-}}=0.3\mathrm{M}$, soThe factor is sixty, not ten: the squared charge and the three counterions per formula unit both matter. Using the given $5\mathrm{nm}$ reference length directly, and $T_1\simeq 298\mathrm{K}$, $T_2=323\mathrm{K}$,The modest temperature increase raises the screening length by only about $4\%$ at fixed concentration; the large increase in ionic strength shortens it by almost an order of magnitude overall.

Under charge stabilization of colloids, the electrostatic double-layer repulsion now has a much shorter range. Attractive Van der Waals forces can then dominate at closer separations, reducing the barrier to particle aggregation in DLVO theory. Thus an electrostatically stabilized colloid is likely to be destabilized, with aggregation more likely than in the original dispersion.

This is a conditional physical prediction, not a proof of instability from a length alone. No particle radius, surface potential, attractive interaction strength or steric coating is specified. Steric stabilization could preserve stability, and trivalent-ion adsorption can alter the surface charge. Moreover, ideal screening theory is only an estimate at $I\simeq 0.6\mathrm{M}$, where ion association, hydrolysis and activities can matter. The numerical result uses the supplied reference length and the stipulated unchanged permittivity, rather than replacing either with an unprovided calibration.

Two stabilization mechanisms. Steric stabilization coats each colloid with adsorbed or grafted polymer layers. Bringing two particles close overlaps or compresses those layers, increasing their free energy through loss of chain conformations and, in a suitable solvent, increased local osmotic pressure. The resulting repulsion can keep the cores apart. A well-solvated nonionic coating can remain effective when added salt would screen electrostatic forces. Coatings can fail through desorption, poor solvent conditions or bridging, so simply adding arbitrary free polymer does not guarantee stabilization.

Charge stabilization of colloids instead uses similarly charged surfaces and their diffuse electric double layers. Their approach costs electrostatic and ionic free energy. In DLVO theory, that repulsion competes with attractive Van der Waals forces. It is particularly sensitive to salt concentration, counterion valence and surface ionization.

Bjerrum length. The standard Bjerrum length is the separation at which the Coulomb interaction energy of two elementary charges equals thermal energy:At $T\simeq 298\mathrm{K}$ and ${\epsilon }_r=78$, this isIt compares ionic electrostatic coupling with thermal agitation. For ions of signed valences $z_1,z_2$, $U/(k_{B}T)=z_1{z}_2{\ell }_B/r$. The standard reference length uses $e$, not the whole colloid charge: assigning a colloid charge $Ze$ does not change the solvent's tabulated ${\ell }_B$.

Screening and the printed normalization. In a symmetric $z:z$ electrolyte, $z$ is the ionic valence magnitude, $e$ the elementary electric charge, $n_{\infty }$ the bulk number density of each ionic species, ${\epsilon }_0$ the vacuum permittivity, ${\epsilon }_R$ the solvent's relative permittivity, $k_B$ the Boltzmann constant and $T$ absolute temperature. Mobile ions rearrange around a charged surface; the Debye–Hückel screening length ${\kappa }^{-1}$ sets the range of the resulting weak potential.

There is an SI normalization error in the original PDF. Combining the Boltzmann distribution $n_{\pm }=n_{\infty }{e}^{\mp ze\psi /(k_{B}T)}$ with the SI Poisson equation and linearizing at $|ze\psi |\ll k_{B}T$ givesThe extra $4\pi$ in the PDF denominator is inconsistent with those SI charge, permittivity and density conventions. Keeping the printed algebra would make its screening length $\sqrt{4\pi }$ times the SI value at the same ion density. No numerical screening length is required here, and the qualitative dependence ${\kappa }^{-1}\propto n_{\infty }^{-1/2}$ is the same in either normalization. Added salt shortens the range of repulsion, allowing the attractive potential to capture particles more readily.

Interaction-potential sketches. Let $h$ be the surface separation between two spheres of radius $a$. A useful schematic DLVO theory potential has the formaway from the microscopic contact region. Here $A_H>0$ is a Hamaker constant; excluded-core repulsion and molecular details regularize contact. The exponential is the double-layer repulsion and the inverse-gap term the short-gap Van der Waals force contribution. A fixed $B$ illustrates the changing screening range; actual charge regulation can change the repulsive amplitude as well.

At low salt there is a high positive barrier between a deep primary minimum near contact and the dispersed state at $U(\infty )=0$. A shallow secondary minimum can occur outside the barrier and support weak reversible association. Increasing salt narrows the repulsive range and lowers the barrier. Thermally activated primary aggregation becomes appreciable when the barrier is only a few $k_{B}T$; the idealized barrier-free condition defines the fast-aggregation threshold associated with the critical coagulation concentration. Above that threshold, collisions can result in rapid capture in the primary minimum.

Aggregate structures. Well above the critical coagulation concentration, the sticking probability is close to one. In diffusion-limited cluster aggregation, the first encounter attaches particles or clusters at accessible outer branches before they can explore many alternative configurations. The result is an open, tenuous structure. Typical ideal three-dimensional cluster models have mass fractal dimension $d_f\simeq 1.8$.

Near the threshold on its lower-salt, slow-aggregation side, a remaining barrier makes most encounters unsuccessful. Reaction-limited cluster aggregation then permits repeated approaches and orientations before attachment, commonly giving denser clusters with $d_f\simeq 2.1$. Since particle count scales as $N\propto (R_g/a{)}^{d_f}$, the larger exponent means more particles at a given large cluster size. This explanation does not require already formed bonds to rearrange: unsuccessful encounters themselves alter which configurations eventually bind.

The side of the threshold matters. At the ideal barrier-free threshold itself, or just above it, aggregation can already be diffusion limited; proximity to the critical coagulation concentration alone does not universally imply a denser reaction-limited structure. The quoted dimensions are typical of irreversible approximately isotropic cluster models, not exact values for every colloid. Strong bond rearrangement or anisotropic attractions can change this comparison.

Stabilization and screening. Charge stabilization of colloids means suppressing aggregation through repulsion between similarly charged particles and their electric double layers. In DLVO theory, this repulsion competes with attractive Van der Waals forces. A substantial repulsive barrier compared with thermal energy $k_{B}T$ can make aggregation slow, even when the dispersed state is only metastable.

The Debye–Hückel screening length ${\lambda }_D={\kappa }^{-1}$ is the characteristic range over which a weak electric potential is attenuated by rearrangement of mobile ions. It measures the thickness of a diffuse screening cloud, not a molecular diameter. For a plane, a weak potential decays as $e^{-\kappa x}$ into the bulk.

Derivation and the source convention. Put a planar interface at $x=0$, with electrolyte in $x>0$, bulk potential zero, and absolute permittivity $\epsilon ={\epsilon }_0{\epsilon }_r$. For a symmetric $z:z$ electrolyte, let $n_{\infty }$ be the bulk number density of each ionic species. The Boltzmann distribution givesso the mobile charge density isSubstitution into the SI Poisson equation, ${\nabla }^2\psi =-{\rho }_e/\epsilon$, gives the Poisson-Boltzmann equationFor $|ze\psi |\ll k_{B}T$, linearize the hyperbolic sine. The Debye–Hückel approximation is thenHere $z$ is the ionic valence magnitude, $e$ the elementary electric charge, $k_B$ the Boltzmann constant, $T$ absolute temperature, ${\epsilon }_0$ the vacuum permittivity, and ${\epsilon }_r$ the solvent's relative permittivity. For several ionic species, the numerator $2z^2{e}^2{n}_{\infty }$ is replaced by $e^2{\sum }_i{z}_i^2{n}_{i,\infty }$.

The original PDF includes an additional $4\pi$ in the denominator. That factor does not follow from SI Poisson equation with number densities and the stated permittivities. The Bjerrum length provides a useful check: ${\ell }_B=e^2/(4\pi \epsilon k_{B}T)$ and the consistent result is ${\kappa }^2=8\pi {\ell }_B{z}^2{n}_{\infty }$. A Gaussian-unit convention moves factors of $4\pi$ together with the charge and permittivity conventions; it does not justify inserting this extra factor while retaining all the given SI quantities. Both numerical interpretations are reported below rather than silently changing the source.

Conditions favoring stability. Increasing the surface-potential or surface-charge magnitude usually strengthens the electrostatic barrier. Lower ionic strength increases ${\lambda }_D$ and generally makes repulsion longer ranged; added salt and especially multivalent counterions shorten it and can promote aggregation. Suitable pH can keep an ionizable surface charged, while conditions near an isoelectric state weaken charge repulsion. The solvent-dependent Van der Waals forces, charge regulation and particle geometry also affect the barrier, so a long screening length is not by itself sufficient. At fixed density and valence, the linear formula gives ${\lambda }_D\propto \sqrt{\epsilon T}$, but actual temperature changes can also alter the permittivity and surface chemistry.

Surface charge at small potential. Take the usual one-sided planar-interface model, with negligible electric displacement on the nonaqueous side. The solution decaying into the solution isThe electric field is $E_x=-{\psi }^{\prime }=\kappa {\psi }_0{e}^{-\kappa x}$. The electrostatic interface boundary conditions then give the surface charge densityOne can also derive it by integrating the diffuse charge density:Thus the diffuse layer neutralizes the surface charge. This is the planar Debye-Hückel double layer. If instead the object is an isolated sheet wetted by identical electrolyte on both sides, each side screens half its total charge and ${\sigma }_{\mathrm{sheet}}=2\epsilon \kappa {\psi }_0$. The one-sided assumption is the conventional local model for the exterior of a large charged colloid; it must be stated because the question does not specify the material on the other side.

Numerical ion density. For $T=298.15\mathrm{K}$, $z=1$, $\kappa =10^9{\mathrm{m}}^{-1}$ and $\epsilon =80{\epsilon }_0=7.0834\times 10^{-10}\mathrm{F}{\mathrm{m}}^{-1}$, the consistent SI result isThis is $\boxed{0.0943\mathrm{mol}{\mathrm{L}}^{-1}}$ for each ion species, and hence also the salt concentration for a monovalent $1:1$ salt. The combined number density of both signs is $2n_{\infty }$.

Literal substitution into the formula printed in the PDF instead givesThese values are distinct consequences of the source's extra factor; they cannot both describe the same physical screening length in SI with the same definition of $n_{\infty }$.

Surface potential and validity. The requested linear estimate in the one-sided model isHowever, the thermal voltage is $k_{B}T/e=25.7\mathrm{mV}$, so $e{\psi }_{0,\mathrm{linear}}/(k_{B}T)\simeq 13.7$. Thus the small-potential approximation is not self-consistent for the supplied surface charge. The boxed value is its formal extrapolation, not a reliable prediction of the actual surface potential.

For a valid nonlinear comparison within the ideal one-sided planar Poisson-Boltzmann equation, multiply the differential equation by ${\psi }^{\prime }$ and use zero potential and field in the bulk. This givesUsing the surface boundary condition yields the Grahame equation and its inverse:With the SI screening length, the latter gives $\boxed{{\psi }_{0,\mathrm{PB}}\simeq 0.135\mathrm{V}}$. This is a nonlinear mean-field estimate under explicit assumptions, not a unique experimentally determined potential: compact interfacial layers, finite ion size and surface chemistry were not supplied. The source's normalization problem and its use of a linear formula far outside the small-potential regime are separate issues.