---
title: 'Smoluchowski Equation: Coagulation & Diffusion'
url: https://www.emergentmind.com/topics/smoluchowski-equation
type: topic
---

# Smoluchowski Equation: Coagulation & Diffusion

The Smoluchowski equation denotes two major, historically related classes of equations. In aggregation theory, the Smoluchowski coagulation equation is a mean-field population-balance model for the time evolution of cluster-size distributions under binary coalescence; in stochastic dynamics, the Smoluchowski equation is also the overdamped Fokker–Planck equation for Brownian motion in a potential. Across these usages, the central structures are gain–loss balance, diffusion or drift, weak formulations, and the competition between conservation laws and singular phenomena such as gelation, flux to infinity, trapping, or phase transition [2506.08017] [1507.03371].

## 1. Classical coagulation equation and weak formulations

In its continuous homogeneous form, the coagulation equation evolves a nonnegative size density \(f(t,m)\) or \(u(x,t)\) on \((0,\infty)\) by
\[
\partial_t f(t,m)
= \frac{1}{2}\int_0^m K(m',m-m')\,f(t,m')\,f(t,m-m')\,dm'
- f(t,m)\int_0^\infty K(m,m')\,f(t,m')\,dm',
\]
with symmetric nonnegative kernel \(K\) [1408.5228] [2506.08017]. In the discrete setting, with concentrations \(c_m(t)\) or \(f_m(x,t)\), the same gain–loss structure appears as
\[
\partial_t c_m(t)
= \frac{1}{2}\sum_{i+j=m} K(i,j)c_i(t)c_j(t)
- c_m(t)\sum_{j\ge 1}K(m,j)c_j(t).
\]
The first term counts formation of size \(m\) by mergers; the second removes size \(m\) by mergers with all partners [1401.0829] [1010.0598].

The measure-valued weak formulation is fundamental because the formal identity obtained by testing with \(\phi(x)=x\) does not automatically justify mass conservation. For a nonnegative Radon measure \(f_t\), the weak equation is
\[
\frac{d}{dt}\langle f_t,\phi\rangle
=
\frac12\iint_{(0,\infty)^2}
\big(\phi(x+y)-\phi(x)-\phi(y)\big)\,K(x,y)\,f_t(dx)\,f_t(dy),
\]
for suitable bounded measurable test functions \(\phi\) [1010.0598] [2501.03867]. In the generalized-moment framework of 2025, one instead writes truncated weighted moments
\[
m^b(r,t)=\int_0^r b(x)\,u(x,t)\,dx
\]
and requires
\[
m^b(r,t)=\int_0^r b(x)u_0(x)\,dx+\int_0^t\int_0^r b(x)Q[u](x,s)\,dx\,ds,
\]
for all continuous weights \(b\), all \(r>0\), and all \(t\in[0,T]\). This definition is proved equivalent to the Escobedo–Laurençot–Mischler–Perthame weak formulation [2506.08017].

Moment functionals organize much of the theory. For \(p\ge 0\),
\[
M_p(t)=\int_0^\infty x^p\,u(x,t)\,dx,
\]
so \(M_0\) is the total number of clusters and \(M_1\) is the total mass [2506.08017]. In broad classes of non-gelling kernels, \(M_1\) is conserved; in stronger-growth regimes, \(M_1\) may decrease because mass is transferred to an infinite cluster. The formal distinction between “mass-conserving” and “gelling” kernels is already visible at the level of homogeneous kernels \(K(ax,ay)=a^\gamma K(x,y)\), where the degree \(\gamma\) is a first classifier [2506.08017] [2501.03867].

## 2. Spatially inhomogeneous and measure-valued coagulation–diffusion equations

The spatially inhomogeneous Smoluchowski coagulation–diffusion equation augments the homogeneous coagulation operator with mass-dependent diffusion in physical space. In continuous mass variables,
\[
\partial_t f(t,x,m)=D(m)\,\Delta_x f(t,x,m)+Q^+(f)(t,x,m)-Q^-(f)(t,x,m),
\]
and in the discrete version,
\[
\partial_t f_m(x,t)=d(m)\Delta_x f_m(x,t)
+\frac12\sum_{i+j=m}K(i,j)f_i f_j
-f_m\sum_{j\ge1}K(m,j)f_j.
\]
At fixed spatial point \(x\), coagulation is exactly the homogeneous Smoluchowski mechanism; the novelty is the parabolic transport and the possibility of mass-dependent diffusivities [1401.0829] [1408.5228].

Norris formulated a general measure-solution theory on a type space \(E\), with measurable mass \(m(y)\), diffusivity \(a(y)\), and symmetric finite coagulation kernel \(K(y,y',dz)\) satisfying \(m(z)=m(y)+m(y')\) \(K(y,y',\cdot)\)-a.e. The mild equation is
\[
\mu_t+\int_0^t P_{t-s}K^-(\mu_s)\,ds
=
P_t\mu_0+\int_0^t P_{t-s}K^+(\mu_s)\,ds,
\]
where \(P_t\) is the heat semigroup with type-dependent diffusivity [1408.5228]. Under the structural bound \(K(y,y',E)\le w(y)w(y')\), local boundedness assumptions on \(a\), and weighted integrability of the initial data, there exists a maximal strong solution, uniqueness holds within the strong-solution class, and mass is conserved when the initial mass–weight moment is finite [1408.5228].

A particularly important example is the Einstein–Smoluchowski model for Brownian coagulation of spherical colloids in \(\mathbb{R}^3\). There,
\[
D(m)=\frac{k_B T}{6\pi\eta\,r(m)},\qquad r(m)\propto m^{1/3},
\]
and the Brownian kernel is
\[
K(m,m')=4\pi\big(D(m)+D(m')\big)\big(r(m)+r(m')\big).
\]
With \(a(y)=2D(m(y))\) and a suitable choice of weight \(\phi(m)=m^{1/6}+m^{5/6}\), the Norris framework accommodates the singular behavior as \(m\downarrow0\), yields global existence, and yields mass conservation when
\[
\int_{\mathbb{R}^3\times(0,\infty)} \big[m^{2/3}+m^{4/3}\big]\,\mu_0(dx,dm)<\infty
\]
[1408.5228].

The coagulation–diffusion PDE also arises as a kinetic limit of many Brownian particles undergoing short-range coagulation. In \(d\ge3\), with \(N\varepsilon^{d-2}=Z\), Hammond and Rezakhanlou show that the empirical measures converge to weak solutions with macroscopic kernel
\[
\beta(n,m)=\alpha(n,m)\int_{\mathbb{R}^d}V(x)\,[1-u_{n,m}(x)]\,dx,
\]
where \(u_{n,m}\) solves
\[
-\Delta u_{n,m}(x)=\frac{\alpha(n,m)}{d(n)+d(m)}V(x)\,[1-u_{n,m}(x)],
\qquad u_{n,m}(x)\to0\ \text{as}\ |x|\to\infty.
\]
The factor \(1-u_{n,m}\) is the survival correction accounting for short-range depletion of near-collision pairs [1401.0829].

## 3. Mass conservation, fluxes, and gelation

A central issue in coagulation theory is that the equation is formally mass-conserving but may lose finite-size mass through gelation. The 2025 generalized-moment analysis introduces a mass flux
\[
J[u](x,t)=\int_0^x\int_{x-y}^\infty y\,K(y,z)\,u(y,t)\,u(z,t)\,dz\,dy,
\qquad J[u](0,t)=0,
\]
and proves the weak conservation law
\[
x\,\partial_t u(x,t)=-\partial_x J[u](x,t)
\]
in the sense of distributions [2506.08017]. Moreover,
\[
\int_0^t J[u](x,s)\,ds
=
\int_0^x y\big(u_0(y)-u(y,t)\big)\,dy,
\]
so that formally
\[
M_1(t)-M_1(0)=-\int_0^t J[u](\infty,s)\,ds.
\]
This identifies gelation as flux of mass to infinity [2506.08017].

The generalized-moment method yields sharp sufficient conditions for conservation and gelation. If a continuous weight \(b\) satisfies \(M^b(0)<\infty\) and the kernel obeys inequalities of the form
\[
(x+y)K(x,y)\le C_1(b(x)+x+1)(b(y)+y+1),
\]
\[
\big(b(x+y)-b(x)-b(y)\big)K(x,y)\le C_2 K_2^b(x,y),
\]
then any weak solution conserves mass:
\[
M_1(t)=M_1(0),\qquad M^b(t)<\infty.
\]
If instead \(b\) is nondecreasing, the kernel is nondegenerate at large sizes, and
\[
\big(b(x+y)-b(x)-b(y)\big)K(x,y)\le -\lambda xy+\mu(x+y+1),
\]
then gelation occurs; if \(\mu=0\), finite-time gelation is guaranteed whenever
\[
T>\frac{2M^b(0)}{\lambda M_1(0)^2},
\]
and one may write
\[
M_1(t)=M_1(0)-g(t),\qquad g(t)=\int_0^t J[u](\infty,s)\,ds\ge0.
\]
These criteria apply to possibly inhomogeneous kernels, not only to homogeneous ones [2506.08017].

For homogeneous kernels
\[
K(x,y)=x^\alpha y^\beta+x^\beta y^\alpha,\qquad 0\le\alpha\le\beta\le1,
\]
the generalized-moment framework recovers the classical threshold: if \(\alpha+\beta\le1\), then mass is conserved; if \(\alpha+\beta>1\), the regime is consistent with gelation [2506.08017]. The additive kernel \(K(x,y)=x+y\) sits at the critical line \(\alpha+\beta=1\) and conserves mass; the multiplicative kernel \(K(x,y)=xy\) lies beyond it and is a canonical gelling example [2506.08017].

A separate 2025 analysis revisits the Escobedo–Mischler–Perthame proof and gives a short deterministic criterion for finite-time gelation. If
\[
H(a):=a\inf\{K(x,y):x,y\in[a,ra]\},\qquad
k:=\int_{x_0}^\infty [H(a)]^{-1/2}\,da<\infty,
\]
then any weak solution with \(M_1(f_0)<\infty\) and \(f_0((x_0,\infty))>0\) loses mass in finite time [2501.03867]. For continuous homogeneous kernels of degree \(\gamma>1\) with \(K(1,1)>0\), this condition holds automatically, so any weak solution gels [2501.03867]. At the critical homogeneity \(\gamma=1\), the power of the logarithmic correction required for gelation depends on the kernel shape: for \(K_1(x,y)=x+y\), the threshold is \(\alpha_c=1\), whereas for the near-diagonal band kernel
\[
K_0(x,y)=(x\wedge y)\mathbf{1}_{x/y\in[1/2,2]},
\]
the threshold is \(\alpha_c=2\) for monodisperse data [2501.03867].

## 4. Self-similar, stationary, and forced regimes

Self-similarity organizes long-time behavior in non-gelling regimes and singular behavior near gelation. For the multiplicative family
\[
K(x,y)=(xy)^s,\qquad s<\frac12,
\]
the kernel is homogeneous of degree \(\lambda=2s<1\), so gelation is absent and mass is conserved [2212.12581]. Global-in-time self-similar solutions are sought as
\[
n(x,t)=t^\alpha f(t^\beta x),
\]
with
\[
\beta=-\frac{1}{1-2s},\qquad \alpha=2\beta.
\]
Equivalently,
\[
n(x,t)=S(t)^{-2}f\!\left(\frac{x}{S(t)}\right),\qquad S(t)\sim t^{1/(1-2s)}.
\]
For \(s<0\), the profile has a three-region structure consisting of very-fast decay near the origin, a lognormal intermediate region, and a gamma-distribution tail; for \(s\in(0,1/2)\), it has a singular origin, a Pareto-like intermediate region, and the same gamma tail [2212.12581].

Near gelation, self-similarity becomes more delicate. For kernels
\[
K(x,y)=(xy)^\gamma,\qquad 1\ge\gamma>\frac12,
\]
the equation gels in finite time, and the similarity ansatz
\[
n(x,t)\sim t'^{\,\beta(2\gamma+1)-1}\,\psi(x\,t'^\beta),\qquad t'=t_g-t,
\]
leads to an anomalous drift exponent
\[
\nu=\frac{1-(2\gamma-1)\beta}{\beta}.
\]
The Laplace-transform formulation yields a nonlocal similarity equation for \(\Phi\) and its fractional integral \(F\), and the analysis identifies a stable lower branch continued from the exactly solvable \(\gamma=1\) product kernel together with an upper branch near \(\gamma\approx1/2\) that violates scaling relations often assumed universal [2212.12714]. In particular, the linear laws \(\sigma=(2\gamma-1)/2\) and \(\nu=\gamma-1/2\) are not generally valid [2212.12714].

Open and forced versions of the coagulation equation replace exact mass conservation by flux balance. With a constant flux of dust particles entering at zero size, one obtains weak flux solutions satisfying
\[
\int_{(0,z]}x\,f_t(dx)-\int_{(0,z]}x\,f_0(dx)
=
-\int_0^t J_{f_s}(z)\,ds+Jt.
\]
For a broad non-gelling class of kernels with homogeneity \(\gamma<1\) and \(|\gamma+2\lambda|<1\), there are time-dependent flux solutions with linearly increasing mass
\[
M(t)=M(0)+Jt,
\]
and size-averaged time-integrated densities are bounded above by the explicit constant-flux profile
\[
f(x)\propto x^{-(\gamma+3)/2}.
\]
For the constant kernel, the Bernstein transform solves a Riccati ODE and converges to the stationary flux measure \(f_\infty(x)\propto x^{-3/2}\) [2412.07745].

In discrete forced coagulation with source \(S_k\) and removal \(r_k c_k\),
\[
\frac{d}{dt}c_k
=
\frac12\sum_{\ell=1}^{k-1}a_{k-\ell,\ell}c_{k-\ell}c_\ell
-c_k\sum_{\ell\ge1}a_{k,\ell}c_\ell
+S_k-r_kc_k,
\]
well-posedness holds for a large class of kernels \(a_{k,\ell}\le A_*(k^\alpha \ell^\beta+k^\beta \ell^\alpha)\) and removal laws \(r_k\ge R_*k^\gamma\), and under a suitable smallness condition solutions converge exponentially to a unique equilibrium [1801.03083]. For aggregation–fragmentation steady states with steeply decaying tails, a fast two-stage numerical algorithm constructs a small “seed system” and extends it through a quadratic recurrence for \(n_{k+1}\), yielding orders-of-magnitude speedups over standard schemes [1503.02421].

## 5. Particle systems, stochastic approximations, and numerical methods

A principal finite-particle approximation to coagulation is the Marcus–Lushnikov process. If
\[
\mu_t^n=\frac1n\sum_{k=1}^{N(t)}\delta_{X_k(t)},
\]
then coalescence of masses \(x\) and \(y\) replaces \(\delta_x+\delta_y\) by \(\delta_{x+y}\), and the empirical measure satisfies a martingale formulation converging to the weak Smoluchowski equation [1010.0598]. For homogeneous-like kernels with degree up to \(1\), Cepeda and Fournier introduce a Wasserstein-type distance \(d_\lambda\) and prove an explicit convergence rate
\[
O(n^{-1/2})
\]
uniform on compact time intervals, together with constructive discrete approximations of the initial data achieving the same rate [1010.0598].

For large discrete systems, recent computational work exploits low-rank structure of the kernel. Low-rank Monte Carlo factorizes or upper-bounds the collision kernel as a sum of rank-one terms and samples collisions with segment trees. The per-collision cost is
\[
O(R\log M),
\]
where \(R\) is the factorization rank and \(M\) the number of size classes [2312.02788]. The Brownian benchmark kernel
\[
C_{ij}=\left(i^{1/3}+j^{1/3}\right)\left(i^{-1/3}+j^{-1/3}\right)
\]
is exactly rank \(3\), the linear kernel \(C_{ij}=i+j\) is rank \(2\), and the ballistic kernel admits a low-rank upper bound with acceptance probability at least \(1/\sqrt2\) [2312.02788]. The method remains exact after acceptance–rejection correction and preserves total mass at the particle level [2312.02788].

A related low-rank strategy applies to coagulation combined with Ostwald ripening in volume space:
\[
\frac{\partial n(v,t)}{\partial t}
=
\frac12\int_0^v K(u,v-u)n(u,t)n(v-u,t)\,du
-n(v,t)\int_0^\infty K(v,u)n(u,t)\,du
-\partial_v J(v,t),
\]
with
\[
J(v,t)=G(v,t)n(v,t)-D_v(v,t)\partial_v n(v,t),\qquad D_v=d_0 G.
\]
In the dimensionless variables of the 2026 study, the coupled model admits an exact parametric solution of the form
\[
\Phi(\xi,b)=\frac{n^2(b)}{V(b)}\exp\!\left(-\frac{n(b)}{V(b)}\xi\right),
\]
so the scaled profile is universally exponential,
\[
f(z)=e^{-z},
\]
independent of the initial volume distribution [2601.14011].

The same coagulation formalism has been carried into astrophysical applications. For primordial black-hole clusters, the merger kernel is obtained by velocity-averaging the gravitational-wave capture cross section, yielding
\[
\mathcal{K}_{ij}\propto
\left(\frac{m_i m_j}{m_0^2}\right)^{15/14}
\left(\frac{m_i+m_j}{m_0}\right)^{9/14},
\]
hence an effective homogeneous degree
\[
\lambda=\frac{39}{14}>1.
\]
This implies finite-time runaway in the Smoluchowski sense. The simulations use an event-driven full-conditioning Monte Carlo scheme with waiting time
\[
\tau=\frac{V}{\mathcal{K}_0}\ln\!\left(\frac1{r_1}\right),
\]
and show that mass segregation further shortens the runaway time [2604.01684].

## 6. Overdamped Fokker–Planck, geometric, and active-matter variants

In another established usage, the Smoluchowski equation is the overdamped Fokker–Planck equation for the probability density \(P(x,t)\) of a Brownian particle in a potential \(U(x)\):
\[
\partial_t P(x,t)=D\,\partial_x\!\left[e^{-\beta U(x)}\partial_x\big(e^{\beta U(x)}P(x,t)\big)\right].
\]
For constant \(D\), this is equivalently
\[
\partial_t P=D\,\partial_x^2 P + D\beta\,\partial_x\!\big(U'(x)P\big)
\]
[1507.03371]. With localized time-dependent sinks, exact Laplace-transform methods reduce the PDE to an integral equation for the single function \(P(0,t)\); once \(P(0,t)\) is known, the full density follows by convolution with the Green function. This program has been carried out for a flat potential, a piecewise linear potential, and a parabolic potential, with exact solutions for constant, linear, inverse, and exponential sink strengths in the cases treated [1408.3288] [1408.3290] [1507.03371].

On curved manifolds, the intrinsic Smoluchowski equation reads
\[
\partial_t p = D\,\Delta_M p + \nabla_M\!\cdot(\mu p \nabla_M U).
\]
On the sphere \(S^2\),
\[
\partial_t p = D\,\Delta_{S^2}p + \nabla_{S^2}\!\cdot(\mu p \nabla_{S^2}U),
\]
with conserved probability and flux
\[
J=-D\nabla_{S^2}p-\mu p\nabla_{S^2}U.
\]
A geometry-based algorithm evolves particles in the tangent plane and projects them back to the sphere; it converges in the weak sense, captures curvature-induced effects in both transient and stationary regimes, and reproduces the stationary Boltzmann density \(p_*\propto e^{-U/k_BT}\) [2104.02843].

The same overdamped framework also supports nonlinear orientational dynamics. For the modified Smoluchowski equation on the unit sphere with dipolar potential,
\[
\partial_t f
=
-\nu(|J_f|)\,\nabla_\omega\!\cdot\!\big(f\,\nabla_\omega(\omega\cdot\Omega_f)\big)
+\tau(|J_f|)\,\Delta_\omega f,
\]
where
\[
J_f=\int_{S^{d-1}}\omega f(\omega,t)\,d\omega,\qquad \Omega_f=\frac{J_f}{|J_f|},
\]
equilibria are von Mises–Fisher states \(f_\infty=\rho M_{\kappa\Omega}\) solving the self-consistency condition
\[
\sigma(\kappa)=\rho\,c_d(\kappa).
\]
The theory identifies stable and unstable equilibria, exponential convergence rates, second-order transitions when \(h(|J|)/|J|\) is non-increasing, and first-order transitions with hysteresis when the ratio grows in part of the state space [1212.3920].

For active Brownian swimmers in \(\mathbb{R}^2\), the phase-space Smoluchowski equation becomes
\[
\partial_t P + U_0\hat{\mathbf u}(\varphi)\cdot\nabla P
=
D_B\nabla^2P + D_\Omega \partial_\varphi^2 P.
\]
After marginalization over the angle \(\varphi\), the long-time position density is Gaussian with effective diffusivity
\[
D_{\rm eff}=D_B+\frac{U_0^2}{2D_\Omega},
\]
and the exact mean-square displacement is
\[
\langle r^2(t)\rangle
=
4D_B t + 2U_0^2\left(\frac{t}{D_\Omega}-\frac{1-e^{-D_\Omega t}}{D_\Omega^2}\right).
\]
The kurtosis quantifies non-Gaussianity generated by orientational persistence; for \(\mathrm{Pe}\lesssim0.1\), the deviation from Gaussian scales as \(t^{-2}\) at short times and \(t^{-1}\) asymptotically [1501.07237].

Taken together, these strands show that “Smoluchowski equation” is not a single model but a family of structurally related equations. In coagulation theory it governs mass exchange across size classes, weak solutions, flux to infinity, and gelation; in stochastic transport it governs overdamped diffusion, reactive sinks, curvature, and active matter. The modern literature extends both lines through generalized moments, measure solutions, kinetic limits, low-rank computation, and explicit singular or stationary asymptotics [2506.08017] [1401.0829]

Source: https://www.emergentmind.com/topics/smoluchowski-equation