Pith. sign in

REVIEW 2 major objections 4 minor 23 references

Gelation and Positivity of Solutions to the Discrete Oort--Hulst--Safronov Coagulation Equation

T0 review · 2 major / 4 minor · reviewed 2026-08-01 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read For the discrete OHS coagulation equation, finite-initial-mass solutions lose mass in finite time under a kernel summability condition, and positive cluster sizes form a half-line from the smallest initial cluster.

desk verdict The gelation results are solid and citable; the positivity theorem has a real unproven step that likely sinks it as stated. read the letter →

arxiv 2607.26167 v1 pith:QNXCOKQZ submitted 2026-07-28 math.CA

Mashkoor Ali This is my paper
classification math.CA MSC 34A1234K30
keywords CoagulationequationGelationWeakformulationMassconservationPositivityofsolutionOort–Hulst–SafronovDiscreteKernelconditions
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper proves two structural results for the discrete Oort–Hulst–Safronov (OHS) coagulation equation, a model in which a cluster breaks into monomers when it collides with a larger cluster. First, if the coagulation kernel satisfies an explicit summability condition — essentially that suitably defined block intensities grow fast enough — then every solution with finite initial mass loses mass in finite time (gelation), and the gelation time is bounded by a constant involving the kernel, the initial mass, and the initial tail moment. The condition covers power-law, factorized, logarithmic, and mixed polynomial-logarithmic kernels, including cases not treated in earlier work on this discrete system. Second, if all collision rates are strictly positive, then for every positive time the set of cluster sizes with positive concentration is exactly the half-line {p*, p*+1, ...}, where p* is the smallest cluster present initially; this holds regardless of time, rates, or the arithmetic structure of the initial support. The method adapts a recent deterministic approach from the continuous Smoluchowski equation to the asymmetric flux structure of the discrete OHS system.

What carries the argument

The central analytical object is the weak formulation of Lemma 2.1: for any bounded test sequence (ψ_i), the evolution of Σ ψ_i ω_i(t) is governed by the discrete flux coefficient Φ_ψ(i,j) = (ψ_{i+1}−ψ_i)j − ψ_j for i≥j. Choosing ψ_i = i∧a makes Φ_ψ ≤ 0 and nonzero only for i≥a, yielding the key energy estimate ∫_0^t Δ(a)[Π_s(a)]² ds ≤ 2M_1^{in}, where Δ(a) is the block intensity and Π_t(a) the block occupation measure; summing over scales and applying Cauchy–Schwarz converts this into a uniform L² bound on the tail moment Υ_{2,s}, which forces mass loss in finite time. For positivity, the paper rewrites (1.1) as ˙ω_i = R_i − φ_i ω_i and uses the variation-of-constants identity ω_i(t)E_i(t)

What would settle it

Exhibit a solution satisfying Definition 2.1 for which ∫_0^t φ_i(s)ds = ∞ for some i and t>0 (φ_i non-integrable only where ω_i=0), making E_i(t) infinite and invalidating the division in (5.3) that proves Theorem 5.1.

Watch

Extended reading notes

Core claim

Under the assumptions of Theorem 3.1, whenever the kernel is symmetric, nonnegative, and satisfies ℓ = Σ_{a≥a0} [Δ(a)]^{−1/2} < ∞ for some a0 and r>1 (with Δ(a) = a·inf{Λ_{i,j}: a≤i,j≤⌊ra⌋}), and the initial data satisfy M_1^{in}<∞ and μ0 = Σ_{i>a0}(i−a0)ω_i^{in} > 0, every solution gels in finite time and T_gel ≤ 2ℓ² (r/(r−1))² M_1^{in}/μ0². Theorem 5.1 states that under the strict positivity assumption Λ_{i,j}>0 for all i,j, the positivity set J(t) = {i: ω_i(t)>0} equals {i≥p*} for all t>0, where p* = min{i: ω_i^{in}>0}, independently of t, of the kernel, and of the initial support beyond its minimum. The proof of the positivity theorem relies on a variation-of-constants representation and

Load-bearing premise

The proof of the positivity theorem assumes the integrated loss rate E_i(t)=exp(∫_0^t φ_i(s)ds) is finite for every finite time; the paper takes this to follow from the solution definition, but that definition only guarantees integrability of ω_i φ_i, not of φ_i itself, so a loss rate that blows up where ω_i vanishes would make the variation-of-constants step invalid.

Editorial extensions

If this is right

  • For kernels of the form Λ_{i,j} ≥ C(i∧j)^γ with γ>1, every finite-mass solution gels in finite time, and the bound (3.13) gives T_gel ≤ 8C^{−1} a_0^{1−γ}(γ−1)^{−2}(r/(r−1))² M_1^{in}/μ0².
  • For the multiplicative kernel Λ_{i,j}=ij, choosing a0=1, r=2 gives T_gel ≤ 16 M_1^{in}/μ0² (equation (3.15)).
  • For symmetric polynomial kernels i^α j^β + i^β j^α, gelation occurs whenever α+β>1, a wider range than previously established; when α+β≥2 the result complements the known instantaneous-gelation conclusion with an explicit finite upper bound.
  • Under strict positivity of all rates, at any t>0 the positivity set is the half-line {i≥p*}: every cluster size at least the minimum initial size is present, and all smaller sizes remain absent for all times.
  • For slowly growing kernels Λ_{i,j} ≤ C(i+j) log(e+i∧j) with finite second moment, mass is conserved for all time, delineating the borderline of the gelation regime.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: The half-line positivity result suggests the discrete OHS model has a 'no gaps' property: initial holes in the cluster-size distribution are filled immediately, because species i is produced solely from i−1; this is qualitatively different from Smoluchowski coagulation, where positivity depends on the additive monoid generated by the initial support.
  • Editorial inference: The gelation bound depends on the kernel only through ℓ and on the data only through M_1^{in}/μ0², so it might be possible to estimate gelation times in applications from coarse block-averaged coagulation rates and the tail of the initial distribution, without knowing the full kernel.
  • Editorial inference: The proof of Theorem 5.1 assumes E_i(t)<∞, which the paper asserts follows from the solution definition; if that assertion can fail in borderline cases, the positivity conclusion may still hold but would require a separate argument — checking the validity of this step for the admitted solution class would settle the matter.
  • Editorial inference: Because the method requires Δ(a)>0 (positive rates within each block), kernels vanishing on the diagonal, such as Λ_{i,j}=|i−j|^γ, are excluded; testing numerically whether such kernels gel in finite time would map the boundary of the theory.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The paper studies the discrete Oort–Hulst–Safronov (OHS) coagulation equation. It introduces a weak formulation with an asymmetric flux coefficient and uses it to prove (Theorem 3.1) that every finite-initial-mass solution gels in finite time under the summability condition ℓ = Σ_a [Δ(a)]^{-1/2} < ∞, with an explicit gelation-time bound in terms of the initial tail moment. The theorem is applied to several kernel classes, including power-law, ratio-suppressed, and logarithmic kernels. Section 4 treats the critical logarithmic kernel, proving gelation for α>1 and a mass-conservation condition for slowly growing kernels. Section 5 proves a positivity result: under strict positivity of the kernel, the positivity set for t>0 is {i ≥ p*}, where p* is the minimal initially present cluster size. The paper also compares its gelation results with earlier work by [3] and [11].

Significance. If correct, Theorem 3.1 is a valuable discrete analogue of Fournier's continuous gelation theorem, giving explicit, simultaneous bounds on the gelation time for a broad class of kernels. The positivity theorem is surprisingly clean and sharply contrasts with the arithmetic monoid characterization for Smoluchowski coagulation. The weak formulation with the asymmetric flux is a useful tool. However, two load-bearing proof steps are currently unjustified: the symmetrization identity in Proposition 4.1 is false as printed, and the assertion in Section 5 that E_i(t) = exp(∫_0^t φ_i) < ∞ follows from Definition 2.1 is not correct. These gaps affect the critical-logarithmic gelation result and the entire positivity theorem, so the manuscript needs revision before the claims can be accepted as proved.

major comments (2)
  1. [Section 5, after Eq. (5.2)] The statement 'E_i(t) := exp(∫_0^t φ_i(s) ds), which satisfies E_i(t) ≥ 1 and E_i(t) < ∞ for all finite t by Definition 2.1' is false. Definition 2.1 only requires ∫_0^T ω_i(s) φ_i(s) ds < ∞, not ∫_0^T φ_i(s) ds < ∞. Since φ_i ≥ 0, weighted integrability does not control the unweighted integral when ω_i is small or zero. Consequently the variation-of-constants formula (5.3) is not justified, and Lemma 5.1 and Lemma 5.3, which divide by E_i(t) and use E_i(t) < ∞ to conclude positivity, are not valid as written. For example, a scalar trajectory y(t) = (T-t)^a with φ(t) = a/(T-t), a>1, satisfies ∫φ = ∞, ∫φ y < ∞, is absolutely continuous, and hits zero at T. The proof needs either an additional argument establishing local integrability of φ_i under the stated assumptions, or extra hypotheses on the kernel/initial data that ensure it. This is load-bearing because Theorem 5.1 depends entirely
  2. [Section 4, Proposition 4.1, before Eq. (4.4)] The 'symmetrization identity' printed as 2Σ_{i≥j} ijω_iω_j = Σ_{i,j≤a} ijω_iω_j is false. For example, with a=2 the left side equals 2ω_1^2 + 4ω_1ω_2 + 8ω_2^2 while the right side equals (ω_1 + 2ω_2)^2. The step from the lower bound −Φψ_a Λ ≥ C C_α ij to the estimate (4.4) requires a correct treatment of the summation region: one must restrict the lower bound to pairs i,j ≤ a and use the inequality 2Σ_{i≥j, i,j≤a} ijω_iω_j ≥ (Σ_{i≤a} iω_i)^2. As printed, the proof of Proposition 4.1 is invalid. This is load-bearing for the claimed gelation in the critical logarithmic case α>1.
minor comments (4)
  1. [Theorem 3.1 statement and proof] The theorem and its proof refer to H(a) without defining it; the correct object is Δ(a). The phrase 'H(a)>0' should read 'Δ(a)>0', and the proof step 'Multiplying (3.7) by [H(a)]^{-1/2}' should use [Δ(a)]^{-1/2}. The same notation error appears in Section 4 and in Example 3.1(5).
  2. [Example 3.1(3)] The constant in the bound for ℓ appears to be off by a factor of 2 relative to the standard integral comparison ∫_a0^∞ dx/[x log^{1+ε/2}(e+x)] = 2/(ε log^{ε/2}(e+a0)). Please check the displayed constants.
  3. [Section 5, Lemma 5.3] The base case of the induction uses 'Since E_{p*}(t) < ∞' to divide; this relies on the same unjustified integrability of φ_{p*}. Once the major issue is resolved, the wording should be updated to make the finiteness of E_i explicit or avoid the integrating factor.
  4. [Section 4, Proposition 4.2] The proof is said to adapt Steps 1–3 of Proposition 7 in [16], but the adaptation is sketched rather than fully detailed. In particular, the derivation of the Gronwall inequality and the dominated-convergence argument would benefit from more detail to be fully checkable.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: gelation and positivity proofs are self-contained from the weak formulation; the Section 5 E_i-finiteness assertion is a non-circular technical gap.

full rationale

The central derivation chain is self-contained. Theorem 3.1 derives the gelation-time bound directly from Lemma 2.1 (the weak formulation) with test sequence ψ_i = i ∧ a; the block intensity Δ(a) and the initial tail moment μ0 are hypotheses, not outputs, and the explicit bound (3.3) is proved rather than assumed. Proposition 4.1 similarly constructs a test sequence and derives an a priori L²-in-time bound on M_1 from the kernel lower bound, with no fitted parameter. The positivity theorem is a structural consequence of the variation-of-constants identity (5.3) and equivalence (5.4) under assumption (5.1); it does not assume the conclusion. Citations to the author's earlier works [2], [3] are contextual background on well-posedness and prior gelation criteria; the new theorems are conditional on any solution and are proven from equations, so these citations are not load-bearing. The only flagged issue is non-circular: after Eq. (5.2), the paper asserts that E_i(t) = exp(∫_0^t φ_i(s)ds) is finite 'by Definition 2.1', but Definition 2.1 only guarantees ∫_0^T ω_i φ_i < ∞, not ∫_0^T φ_i < ∞. If φ_i were non-integrable on the zero set of ω_i, E_i could be infinite, breaking Lemma 5.3's division step. This is an omitted justification in the positivity proof, not a reduction of the theorem to its inputs, and it does not affect the gelation results.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

No free parameters are fitted to data; the paper is proof-based. The main extra-logical inputs are the solution concept, the well-posedness background from earlier papers, and the E_i-finiteness assertion in the positivity section.

assumptions (3)
  • domain assumption The weak formulation Lemma 2.1 is valid: term-by-term summation and integration interchange is justified by Definition 2.1's integrability condition.
    The proof asserts the interchange is justified by the integrability condition of Definition 2.1, but does not give a full justification for the infinite double-sum.
  • domain assumption Global solutions in the sense of Definition 2.1 exist for the kernel classes considered in Examples 3.1 and Propositions 3.2–3.3.
    The theorems are conditional on any solution; existence is not proved here and is cited to prior works [2,3,5,11,12], which may not cover all new kernel classes (e.g., ratio-suppressed kernels).
  • ad hoc to paper E_i(t) = exp(∫_0^t φ_i) < ∞ for all finite t under Definition 2.1 (stated after (5.2)).
    Load-bearing for Theorem 5.1's division step in (5.3); not implied by the stated condition ∫ω_iφ_i<∞, which permits non-integrable φ_i on sets where ω_i=0.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Gelation and Positivity of Solutions to the Discrete Oort--Hulst--Safronov Coagulation Equation." pith.science (2026). https://pith.science/paper/QNXCOKQZ

@misc{pith2026260726167,
  author       = {Pith},
  title        = {Pith review of: Gelation and Positivity of Solutions to the Discrete Oort--Hulst--Safronov Coagulation Equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QNXCOKQZ}},
  note         = {Machine review of arXiv:2607.26167}
}
read the original abstract

Motivated by the recent deterministic approach of Fournier~\cite{F2025} to gelation for the continuous Smoluchowski coagulation equation, we adapt his method to the discrete Oort--Hulst--Safronov (OHS) coagulation system. We show that under a suitable condition on the coagulation kernel, every solution with finite initial mass loses mass in finite time, and we give an explicit bound on the gelation time. We also prove gelation in the critical logarithmic case and provide a sufficient condition for mass conservation. Finally, we study the positivity of solutions and show that, for any positive time, a cluster size has positive concentration if and only if it is at least as large as the smallest cluster present initially.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

23 extracted references · 1 linked inside Pith

  1. [3]

    K., Global existence of solutions to the discrete Safronov-Dubovski ˇi coagulation equations and failure of mass-conservation,J

    Ali, M., Giri, A. K., Global existence of solutions to the discrete Safronov-Dubovski ˇi coagulation equations and failure of mass-conservation,J. Math. Anal. Appl.,519(1), 126755, 2023

    2023

  2. [11]

    Das, A., Saha, J., The discrete Safronov-Dubovski ˇi aggregation equation: instantaneous gelation and nonexistence theorem,J. Math. Anal. Appl.,514(1), 126310, 2022

    2022

  3. [1]

    Aldous, D., Deterministic and stochastic models for coalescence (aggregation and coagulation): a review of the mean-field theory for probabilists,Bernoulli,5, 3–48, 1999

    1999

  4. [2]

    K., On the discrete Safronov-Dubovski ˇi coagulation equations: well- posedness, mass conservation, and asymptotic behavior,Math

    Ali, M., Rai, P., Giri, A. K., On the discrete Safronov-Dubovski ˇi coagulation equations: well- posedness, mass conservation, and asymptotic behavior,Math. Methods Appl. Sci.,47(6), 5191–5206, 2024. 18

    2024

  5. [4]

    Andreis, L., Iyer, T., Magnanini, E., Gelation in cluster coagulation processes,arXiv:2308.10232, 2023

    arXiv 2023

  6. [5]

    Methods Appl

    Bagland, V., Convergence of a discrete Oort-Hulst-Safronov equation,Math. Methods Appl. Sci., 28(13), 1613–1632, 2005

    2005

  7. [6]

    Bagland, V., Lauren¸ cot, Ph., Self-similar solutions to the Oort-Hulst-Safronov coagulation equation, SIAM J. Math. Anal.,39, 345–378, 2007

    2007

  8. [7]

    Banasiak, J., Lamb, W., Lauren¸ cot, Ph.,Analytic Methods for Coagulation-Fragmentation Models, Volumes 1 & 2, CRC Press, Boca Raton, 2019

    2019

Show all 23 references
  1. [8]

    K., Rai, P., Giri, A

    Barik, P. K., Rai, P., Giri, A. K., Mass-conserving weak solutions to the Oort-Hulst-Safronov coagu- lation equation with singular rates,J. Differ. Equ.,11(5), 1125–1138, 2022

    2022

  2. [9]

    P., On the positivity of solutions to the Smoluchowski equations,Mathematika,42, 406–412, 1995

    Da Costa, F. P., On the positivity of solutions to the Smoluchowski equations,Mathematika,42, 406–412, 1995

    1995

  3. [10]

    P., Mathematical aspects of coagulation-fragmentation equations, inMathematics of Energy and Climate Change, Springer International Publishing, Cham, 83–162, 2015

    Da Costa, F. P., Mathematical aspects of coagulation-fragmentation equations, inMathematics of Energy and Climate Change, Springer International Publishing, Cham, 83–162, 2015

    2015

  4. [12]

    Davidson, J., Existence and uniqueness theorem for the Safronov-Dubovski ˇi coagulation equation, Z. Angew. Math. Phys.,65(4), 757–766, 2014

    2014

  5. [13]

    Drake, R. L., A general mathematical survey of the coagulation equation, inTopics in Current Aerosol Research, Part 2, International Reviews in Aerosol Physics and Chemistry, Pergamon Press, Oxford, 203–376, 1972

    1972

  6. [14]

    B., Structural stability of disperse systems and finite nature of a coagulation front,J

    Dubovski ˇi, P. B., Structural stability of disperse systems and finite nature of a coagulation front,J. Experim. Theor. Phys.,89(2), 384–390, 1999

    1999

  7. [15]

    B., A ‘triangle’ of interconnected coagulation models,J

    Dubovski ˇi, P. B., A ‘triangle’ of interconnected coagulation models,J. Phys. A: Math. Gen.,32(5), 781–793, 1999

    1999

  8. [16]

    Math´ ematique,363(G6), 583–591, 2025

    Fournier, N., On gelation for the Smoluchowski coagulation equation,Comptes Rendus. Math´ ematique,363(G6), 583–591, 2025

    2025

  9. [17]

    P., Theoretical analysis of a discrete population balance model with sum kernel,Port

    Kaushik, S., Kumar, R., da Costa, F. P., Theoretical analysis of a discrete population balance model with sum kernel,Port. Math.,80(3), 343–367, 2023

    2023

  10. [18]

    Lachowicz, M., Lauren¸ cot, Ph., Wrzosek, D., On the Oort-Hulst-Safronov coagulation equation and its relation to the Smoluchowski equation,SIAM J. Math. Anal.,34, 1399–1421, 2003

    2003

  11. [19]

    Lauren¸ cot, Ph., Convergence to self-similar solutions for a coagulation equation,Z. Angew. Math. Phys.,56, 398–411, 2005

    2005

  12. [20]

    Lauren¸ cot, Ph., Self-similar solutions to a coagulation equation with multiplicative kernel,Physica D,222, 80–87, 2006

    2006

  13. [21]

    H., Van de Hulst, H

    Oort, J. H., Van de Hulst, H. C., Gas and smoke in interstellar space,Bull. Astronom. Inst. Nether- lands,10, 187–210, 1946

    1946

  14. [22]

    Probab.,41(3B), 1806–1830, 2013

    Rezakhanlou, F., Gelation for Marcus-Lushnikov process,Ann. Probab.,41(3B), 1806–1830, 2013

    2013

  15. [23]

    S.,Evolution of the Protoplanetary Cloud and Formation of the Earth and the Planets, Israel Program for Scientific Translations, Jerusalem, 1972

    Safronov, V. S.,Evolution of the Protoplanetary Cloud and Formation of the Earth and the Planets, Israel Program for Scientific Translations, Jerusalem, 1972. Jindal Global Business School, O.P. Jindal Global University, Sonipat–131001, Haryana, India Email address:mashkoor.al...

    1972

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.