Coagulation equations for non-spherical clusters

Technical Summary: Coagulation Equations for Non-Spherical Clusters

Problem Statement
This work investigates the long-time asymptotics of a coagulation model describing a system of particles characterized by two variables: volume ($v$) and surface area ($a$). Unlike standard one-dimensional coagulation models where particles are defined solely by volume, this model incorporates the geometric evolution of non-spherical particles. The aggregation mechanism is modeled in two distinct stages:

1. Collision: Two particles merge at a contact point, conserving both total volume and total surface area ($v_{new} = v_1 + v_2$, $a_{new} = a_1 + a_2$).
2. Fusion: Following collision, the particle geometry relaxes to reduce surface area while preserving volume, driving the particle toward a spherical shape. This process is governed by a "fusion term" in the evolution equation.

The particles are constrained by the isoperimetric inequality, $a \geq c_0 v^{2/3}$ (where $c_0 = (36\pi)^{1/3}$), ensuring physical realizability. The governing equation for the particle density $f(a, v, t)$ is:
$$\partial_t f + \partial_a [r(a, v)(c_0 v^{2/3} - a)f] = K[f]$$
where $K[f]$ is the classical Smoluchowski coagulation operator and $r(a, v)$ is the fusion rate. The study focuses on kernels with power-law scaling behaviors and analyzes the existence of self-similar solutions and the emergence of different asymptotic regimes based on the parameters of the fusion rate.

Methodology
The authors employ a rigorous mathematical analysis centered on the existence and properties of self-similar profiles. The methodology involves:

1. Scaling Analysis: The authors seek self-similar solutions of the form $f(a, v, t) = (1+t)^{-8/3\xi} g(a(1+t)^{-2/3\xi}, v(1+t)^{-\xi})$, where $\xi = 1/(1-\gamma)$ and $\gamma$ is the homogeneity degree of the coagulation kernel. This reduces the time-dependent partial differential equation to a stationary integro-differential equation for the profile $g$.
2. Truncation and Fixed-Point Arguments: To prove the existence of solutions, the authors introduce truncated versions of the coagulation kernel and the fusion term. They utilize a fixed-point theorem (specifically a variant of Tychonoff's theorem) on a compact set of measures equipped with the weak-$\ast$ topology.
3. Moment Estimates: A critical component of the proof involves deriving uniform moment estimates for powers of volume and area. The authors distinguish between cases where the fusion rate grows rapidly ($\mu > 0$) versus slowly or decays ($\mu \leq 0$).
   • For $\mu > 0$, the rapid growth of the fusion rate for large areas allows for the control of the total surface area, preventing particles from becoming arbitrarily non-spherical.
   • For $\mu \leq 0$, the fusion rate is insufficient to dominate coagulation for large particles, leading to different behaviors depending on initial data and the magnitude of the fusion constant.
4. Ramification Analysis: For the case $\mu < 0$ with specific initial conditions, the authors derive lower bounds on the total surface area to demonstrate "ramification," where the system evolves into a state where surface area grows significantly faster than volume, deviating from spherical shapes.

Key Contributions and Results

1. Existence of Self-Similar Profiles ($\mu > 0$):
   The paper proves the existence of volume-conserving self-similar profiles when the fusion rate parameter $\mu$ is positive. In this regime, the fusion term dominates coagulation for large particles, keeping the particle distribution close to the isoperimetric boundary ($a \approx c_0 v^{2/3}$). The authors establish that for kernels satisfying specific bounds (weak dependence on area, power-law dependence on volume), self-similar profiles exist with bounded moments for all powers of volume and area.

2. Existence of Self-Similar Profiles with Fast Fusion ($\mu \leq 0$):
   Contrary to the expectation that $\mu \leq 0$ precludes self-similarity, the authors show that if the prefactor of the fusion rate ($R_0$) is sufficiently large, self-similar profiles still exist. This result extends to the case $\mu = 0$ (previously studied in [14]) and demonstrates that a strong initial fusion rate can compensate for the lack of growth in the fusion kernel with respect to area.

3. Ramification Regime ($\mu < 0$):
   The paper identifies a distinct asymptotic regime termed "ramification." If $\mu < 0$ and the initial particle distribution has a sufficiently large surface area relative to its volume, the fusion term becomes negligible compared to coagulation. In this scenario, the total surface area grows exponentially in time ($\langle a \rangle(t) \sim e^{t/3}$), while the volume grows polynomially. This leads to a particle distribution where $a \gg v^{2/3}$, describing a system of highly ramified, fractal-like structures rather than spheres.

4. Moment Bounds and Regularity:
   The authors provide detailed estimates for moments $M_{n,k} = \int a^n v^k g(a,v) da dv$. They demonstrate that for $\alpha > 0$ (strong dependence of coagulation on small volumes), moments with negative powers of volume are bounded. For $\alpha = 0$, singularities may occur at small volumes, and bounds are restricted to moments with powers $k \geq \gamma$.

Significance and Claims
The paper claims to extend the mathematical theory of coagulation equations from one-dimensional (volume-only) models to two-dimensional models incorporating surface area. The primary significance lies in demonstrating that the interplay between collision and fusion rates fundamentally alters the long-time asymptotic behavior of the system.

• Geometric Evolution: The work highlights that the shape of particles (spherical vs. ramified) is not merely a static assumption but a dynamic outcome determined by the relative rates of fusion and coagulation.
• New Phenomena: The identification of the "ramification" regime represents a new mathematical phenomenon not present in standard one-dimensional coagulation models. It shows that under certain conditions, the system naturally evolves toward highly non-spherical, fractal-like structures.
• Mathematical Rigor: The paper provides rigorous proofs for the existence of self-similar profiles in a multi-variable setting, addressing technical challenges related to the non-compactness of the area variable and the specific scaling properties of the fusion term.

The authors explicitly state that their results are mathematical in nature, aiming to understand how much of the existing theory for one-dimensional coagulation can be carried over to the two-dimensional case and to observe the new phenomena arising from the inclusion of surface area dynamics. They do not propose new experimental setups but rather provide a theoretical framework for interpreting systems where particle geometry plays a critical role, such as in aerosol physics or the coalescence of viscous droplets.