# Smoluchowski coagulation equation

The Smoluchowski coagulation equation is an integrodifferential equation that describes how a population of particles or clusters evolves when pairs merge irreversibly into larger ones, with the rate of each merger set by a coagulation kernel depending on the sizes of the two colliding particles. [Marian Smoluchowski](https://www.edgepedia.org/marian-smoluchowski) introduced it in 1916 to describe the clumping of colloidal particles through diffusion, and it has since been applied in physical chemistry, chemical engineering, atmospheric physics, planetary science, and economics.<sup>[1](https://journals.aps.org/pre/abstract/10.1103/PhysRevE.108.064110?ft=1)</sup> A review of coagulation–fragmentation models dates the derivation of the equation to 1917.<sup>[2](https://www.sciencedirect.com/science/article/pii/S0167278906003022)</sup>

| Key fact | Detail |
|---|---|
| Origin | Introduced by Marian Smoluchowski in 1916 for diffusion-driven clumping of colloidal particles<sup>[1](https://journals.aps.org/pre/abstract/10.1103/PhysRevE.108.064110?ft=1)</sup> |
| Governing input | The coagulation kernel K(x, y), the rate at which particles of sizes x and y merge<sup>[1](https://journals.aps.org/pre/abstract/10.1103/PhysRevE.108.064110?ft=1)</sup> |
| Discrete form | \( dc_g/dt = \tfrac{1}{2}\sum_{l=1}^{g-1} K(g-l,l)\,c_{g-l}c_l - c_g \sum_{l=1}^{\infty} K(g,l)\,c_l \)<sup>[3](https://www.sciencedirect.com/science/article/abs/pii/S016727890600306X)</sup> |
| Gelation | For the product kernel K = xy, mass conservation fails after a finite critical time, with the second moment diverging like \( 1/(t_c - t) \)<sup>[3](https://www.sciencedirect.com/science/article/abs/pii/S016727890600306X)</sup> |
| Gelation threshold | For kernels K = (xy)^γ, an infinitely large cluster forms in finite time when 1 ≥ γ > 1/2<sup>[1](https://journals.aps.org/pre/abstract/10.1103/PhysRevE.108.064110?ft=1)</sup> |
| Well-posedness | Norris (1998) gave existence and uniqueness conditions and a non-uniqueness example; Fournier and Laurençot (2006) proved uniqueness for homogeneous kernels of degree λ ∈ (−∞, 2] \ {0} given a finite moment of order λ<sup>[4](https://arxiv.org/abs/math/9801145)</sup><sup> • </sup><sup>[5](https://doi.org/10.1016/j.jfa.2005.07.013)</sup> |
| Main applications | Cloud droplet collision–coalescence, aerosol dynamics, polymerization, colloids<sup>[6](https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2020MS002301)</sup><sup> • </sup><sup>[7](https://arxiv.org/abs/2009.04436)</sup> |

## The equation

In the discrete form, particles come in integer sizes, and \( c_g(t) \) is the concentration of g-mers at time t. The equation reads<sup>[3](https://www.sciencedirect.com/science/article/abs/pii/S016727890600306X)</sup>

\[ \frac{dn_g}{dt} = \tfrac{1}{2}\sum_{j=1}^{g-1} K(g-j,\,j)\,n_{g-j}n_j \; - \; n_g\sum_{j=1}^{\infty} K(g,\,j)\,n_j. \]

The first term on the right is the gain in g-mers from coalescence of (g − l)-mers and l-mers; the second is the loss of g-mers sticking to particles of any other size. The kernel K(g, l) is the transition rate for the merger of a g-mer and an l-mer. Concentrations are often measured in units of the initial monomer concentration and time in units of \( 1/(c_1(0)K(1,1)) \).<sup>[3](https://www.sciencedirect.com/science/article/abs/pii/S016727890600306X)</sup>

The mass variable can also be continuous: f(x) then counts particles per unit volume in the mass range [x, x + dx], and the same gain-and-loss structure applies with x + y replacing the sum of indices. In cloud physics this continuous form is called the kinetic collection equation or the stochastic collection equation.<sup>[6](https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2020MS002301)</sup> A spatially extended variant, the Smoluchowski coagulation–diffusion PDE, models mass-bearing Brownian particles subject to short-range pairwise coagulation and has been derived as a kinetic limit from microscopic models of coagulating Brownian particles.<sup>[8](https://arxiv.org/abs/1401.0829)</sup>

## The coagulation kernel

The kernel K(x, y) encodes the physics of the collision process and is the main input separating one application from another.<sup>[6](https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2020MS002301)</sup> For cloud droplets it is built from a geometric swept volume, involving the drop radius and terminal fall velocity, multiplied by a collision efficiency tabulated from detailed fluid-dynamics trajectory calculations.<sup>[6](https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2020MS002301)</sup> For [diffusion-limited aggregation](https://www.edgepedia.org/diffusion-limited-aggregation) in colloidal systems, Smoluchowski himself derived the Brownian coagulation kernel.<sup>[9](https://pmc.ncbi.nlm.nih.gov/articles/PMC7597352/)</sup>

Mathematical work often restricts attention to homogeneous kernels, those satisfying K(ax, ay) = a^λ K(x, y), because the degree of homogeneity λ largely controls the behavior of solutions. Three kernels are exactly solvable: the constant kernel K = 2, the additive kernel x + y, and the multiplicative (product) kernel xy.<sup>[2](https://www.sciencedirect.com/science/article/pii/S0167278906003022)</sup><sup> • </sup><sup>[5](https://doi.org/10.1016/j.jfa.2005.07.013)</sup>

## Gelation

The equation is formally mass-conserving, but for sufficiently fast-growing kernels solutions lose mass in finite time. This failure, called gelation, corresponds physically to the formation of a giant structure, as in polymer gels, and is commonly described as a phase transition.<sup>[9](https://pmc.ncbi.nlm.nih.gov/articles/PMC7597352/)</sup> For the product kernel K(g, l) ∝ gl, an early exact solution showed that the total mass concentration ceases to be conserved after a finite critical time \( t_c \), while the second moment of the mass distribution diverges as \( \varphi_2(t) \propto 1/(t_c - t) \); at \( t = t_c \) itself, nothing singular happens to the mass spectrum or the particle number concentration.<sup>[3](https://www.sciencedirect.com/science/article/abs/pii/S016727890600306X)</sup>

Hendriks, Ernst, and Ziff, in a 1983 study of coagulation equations with gelation, derived criteria for the occurrence of gelation for classes of kernels, obtained critical exponents in the pre- and post-gelation stages, calculated bounds on the gelation time \( t_g \), and gave an exact post-gelation solution for \( K_{ij} = (ij)^\omega \) with ω > 1/2 and for \( K_{ij} = a^{i+j} \) with a > 1.<sup>[10](https://doi.org/10.1016/0167-2789(83)90245-2)</sup> They argued that \( K_{ij} \sim (ij)^\omega \) with ω ≈ 1 − 1/d, where d is dimensionality, effectively models the sol–gel transformation in polymerizing systems.<sup>[10](https://doi.org/10.1016/0167-2789(83)90245-2)</sup>

Modern results sharpen the threshold. Eggers and Fontelos showed that for kernels \( K = (xy)^\gamma \), the equation produces an infinitely large cluster in finite time when \( 1 \ge \gamma > 1/2 \), and computed the anomalous scaling dimensions of the dominant cluster growth; apart from the branch originating at the exactly solvable case γ = 1, they found a branch of solutions near γ = 1/2 that violates matching conditions for small cluster sizes that had been widely believed to hold universally.<sup>[1](https://journals.aps.org/pre/abstract/10.1103/PhysRevE.108.064110?ft=1)</sup> Fournier gave a short proof, revisiting work of Escobedo, Mischler, and Perthame (2002), that for homogeneous kernels of degree γ > 1 not vanishing on the diagonal, or of degree γ = 1 with a suitable logarithmic factor, any weak solution loses mass in finite time.<sup>[11](https://arxiv.org/abs/2501.03867)</sup> On the conservative side, a generalized moment approach by Kimura and Miyata showed that a wide class of kernels, including those with polynomial or partial exponential growth, still ensures mass conservation, and identified kernels growing like O(x(log x)²) that do induce gelation.<sup>[12](https://link.springer.com/article/10.1007/s13160-025-00744-z)</sup>

## Well-posedness

Because the equation is nonlinear and nonlocal, existence and uniqueness of solutions depend strongly on the kernel. James R. Norris, a probabilist then working on interacting particle systems, proved in 1998 sufficient conditions for existence and uniqueness for a wide class of kernels and initial mass distributions, constructed an example of non-uniqueness with two distinct conservative solutions, and showed that the stochastic coalescent, suitably normalized, converges weakly to the solution of the deterministic equation, giving a statistical derivation of Smoluchowski's equation.<sup>[4](https://arxiv.org/abs/math/9801145)</sup> Norris treated kernels bounded as K(x, y) ≤ φ(x)φ(y) for continuous sublinear φ, which also covers kernels that blow up as x → 0 or y → 0.<sup>[4](https://arxiv.org/abs/math/9801145)</sup>

Nicolas Fournier and Philippe Laurençot, mathematicians working on kinetic and coagulation equations, extended uniqueness in 2006 to measure-valued solutions for homogeneous kernels of degree λ ∈ (−∞, 2] \ {0}, assuming only that the moment of order λ of the initial datum is finite; earlier results of this type covered only the constant, additive, and product kernels.<sup>[5](https://doi.org/10.1016/j.jfa.2005.07.013)</sup> For the constant kernel they also obtained a comparison principle and a contraction property.<sup>[5](https://doi.org/10.1016/j.jfa.2005.07.013)</sup> A survey of coagulation equations for aerosol dynamics notes that well-posedness is now established for large classes of kernels, and studies in that setting also address stationary solutions when a source of small particles is present.<sup>[7](https://arxiv.org/abs/2009.04436)</sup>

## Probabilistic connections

The equation is the mean-field limit of stochastic particle systems. In the Marcus–Lushnikov process, a finite system of coagulating particles with a cut-off kernel, the possible limits as the number of particles grows solve either the Smoluchowski equation or the Flory equation, a variant that handles the post-gelation phase differently.<sup>[13](https://www.sciencedirect.com/science/article/pii/S0304414908000343)</sup> David J. Aldous's 1999 review of deterministic and stochastic coalescence models remains the standard probabilist's reference for this mean-field theory.<sup>[14](https://ftp.icm.edu.pl/packages/EMIS/journals/PS/viewarticlebbf4.html?id=263&layout=abstract)</sup> When spatial structure matters, the coagulation–diffusion PDE arises as the kinetic limit of coagulating Brownian particles.<sup>[8](https://arxiv.org/abs/1401.0829)</sup>

## Applications and limitations

**Cloud physics.** Bin microphysics schemes model collision–coalescence of cloud droplets by solving the Smoluchowski equation, known there as the stochastic collection equation.<sup>[15](https://journals.ametsoc.org/abstract/journals/bams/100/4/bams-d-18-0005.1.xml)</sup> Analytic solutions exist only for very simple collision kernels, so numerical solution is the rule, and the stiffness of the coupled equations forces short substeps in large-scale models whose dynamics time steps run several minutes to tens of minutes.<sup>[16](https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019MS001689)</sup>

The equation has a structural limitation in this setting: it is deterministic and mean-field, while droplet collision–coalescence is a stochastic process in which some droplets undergo unlikely collision sequences and grow faster than average, an effect proposed as important for rapid precipitation onset.<sup>[15](https://journals.ametsoc.org/abstract/journals/bams/100/4/bams-d-18-0005.1.xml)</sup> The kinetic collection equation is in fact only the mean-field approximation of a more general stochastic collection equation that accounts for fluctuations; Lagrangian particle-based methods, such as the superdroplet method, simulate the underlying stochastic process directly and can provide correct solutions even when the well-mixed volume is small.<sup>[6](https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2020MS002301)</sup><sup> • </sup><sup>[16](https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019MS001689)</sup>

**Aerosols and materials.** Binary coagulation, in which two particles merge into one, is a central process in aerosol dynamics, and the size distribution over time is described by the equation.<sup>[7](https://arxiv.org/abs/2009.04436)</sup> In polymerization, the kernel \( K_{ij} \sim (ij)^\omega \) with ω ≈ 1 − 1/d models the sol–gel transformation while approximately accounting for cross-linking and steric hindrance neglected in the classical Flory–Stockmayer theory, which corresponds to ω = 1.<sup>[10](https://doi.org/10.1016/0167-2789(83)90245-2)</sup> Extensions add fragmentation (as in the Becker–Döring equations of 1935), diffusion, sources of fresh monomers, and multicomponent particle composition.<sup>[2](https://www.sciencedirect.com/science/article/pii/S0167278906003022)</sup><sup> • </sup><sup>[7](https://arxiv.org/abs/2009.04436)</sup>

## References

1. [Jens Eggers and Marco Fontelos (2023). Anomalous dimensions of the Smoluchowski coagulation equation. Physical Review E 108, 064110.](https://journals.aps.org/pre/abstract/10.1103/PhysRevE.108.064110?ft=1)
2. [An introduction to mathematical models of coagulation–fragmentation processes, Chaos, Solitons & Fractals.](https://www.sciencedirect.com/science/article/pii/S0167278906003022)
3. [Gelation in coagulating systems, Chaos, Solitons & Fractals.](https://www.sciencedirect.com/science/article/abs/pii/S016727890600306X)
4. [James R. Norris (1998). Smoluchowski's coagulation equation: uniqueness, non-uniqueness and a hydrodynamic limit for the stochastic coalescent. arXiv:math/9801145.](https://arxiv.org/abs/math/9801145)
5. [Nicolas Fournier and Philippe Laurençot (2006). Well-posedness of Smoluchowski's coagulation equation for a class of homogeneous kernels. Journal of Functional Analysis 233(2), 351–379.](https://doi.org/10.1016/j.jfa.2005.07.013)
6. [Axel Seifert and colleagues (2020). Potential and Limitations of Machine Learning for Modeling Warm-Rain Cloud Microphysical Processes. Journal of Advances in Modeling Earth Systems.](https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2020MS002301)
7. [Marina A. Ferreira (2021). Coagulation equations for aerosol dynamics. arXiv:2009.04436.](https://arxiv.org/abs/2009.04436)
8. [Alan Hammond (2014). Coagulation and diffusion: a probabilistic perspective on the Smoluchowski PDE. arXiv:1401.0829.](https://arxiv.org/abs/1401.0829)
9. [The Smoluchowski Ensemble—Statistical Mechanics of Aggregation, PMC.](https://pmc.ncbi.nlm.nih.gov/articles/PMC7597352/)
10. [E. M. Hendriks, M. H. Ernst, and R. M. Ziff (1983). Coagulation Equations with Gelation. Journal of Statistical Physics 31(3).](https://doi.org/10.1016/0167-2789(83)90245-2)
11. [Nicolas Fournier (2025). On gelation for the Smoluchowski equation. arXiv:2501.03867.](https://arxiv.org/abs/2501.03867)
12. [Masato Kimura and H. Miyata (2025). Mass conservation and gelation for the Smoluchowski coagulation equation: a generalized moment approach. Japan Journal of Industrial and Applied Mathematics 42, 1609–1635.](https://link.springer.com/article/10.1007/s13160-025-00744-z)
13. [Marcus–Lushnikov processes, Smoluchowski's and Flory's models, Stochastic Processes and their Applications.](https://www.sciencedirect.com/science/article/pii/S0304414908000343)
14. [Probability Surveys Vol. 14 (2017), including Aldous's review of mean-field coalescence theory.](https://ftp.icm.edu.pl/packages/EMIS/journals/PS/viewarticlebbf4.html?id=263&layout=abstract)
15. [Hugh Morrison and colleagues (2020). Modeling of Cloud Microphysics: Can We Do Better? Bulletin of the American Meteorological Society.](https://journals.ametsoc.org/abstract/journals/bams/100/4/bams-d-18-0005.1.xml)
16. [Hugh Morrison and colleagues (2020). Confronting the Challenge of Modeling Cloud and Precipitation Microphysics. Journal of Advances in Modeling Earth Systems.](https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2019MS001689)

---
*Topic: Encyclopedia › Physical world and mathematics › Mathematics and statistics › Logic and discrete mathematics › General discrete mathematics and discrete structures › Combinatorics › Combinatorics in other fields › Combinatorics of gelation and aggregation*

*Initially written Sep 17, 2026 · Reviewed: Oct 15, 2026 · Edited: Oct 15, 2026 · Last review: Oct 15, 2026*

*Copyright 2026 EdgeChat AI, a subsidiary of Biostate AI.*

License: Edgepedia Community License 1.0, https://www.edgepedia.org/edgepedia/license