pith. sign in

arxiv: 2601.04924 · v3 · submitted 2026-01-08 · ❄️ cond-mat.stat-mech · math.PR

Short-time statistics of extinction and blowup in reaction kinetics

Rotem Degany , Michael Assaf , Baruch Meerson This is my paper

Pith reviewed 2026-05-16 16:14 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech math.PR
keywords stochastic reactionsextinction timeblowup timeWKB approximationmaster equationshort-time tailessential singularityasymptotic analysis
0
0 comments X

The pith

A time-dependent WKB approximation to the master equation captures the essential singularity in short-time extinction and blowup statistics.

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

The authors show how to compute the short-time tail of the probability distribution for extinction or blowup times in stochastic reaction systems, where this tail has an essential singularity at zero. They apply a time-dependent WKB approximation directly to the master equation to capture the singularity, but this leaves a pre-exponential factor undetermined. By using WKB on the Laplace-transformed backward master equation and matching to an inner solution for moderate particle numbers, they determine the full asymptotic form. This is verified on three exactly solvable reaction examples. A sympathetic reader would care because these short-time statistics govern rare events that deterministic rate equations miss in small systems.

Core claim

The singularity at T=0 in P_m(T) is captured by a time-dependent WKB approximation applied directly to the master equation. Accurate results follow when the WKB solution from the Laplace-transformed backward master equation is matched to an inner solution valid for not too large m.

What carries the argument

Time-dependent WKB approximation applied to the master equation, matched to an inner solution via the Laplace-transformed backward equation.

If this is right

  • Accurate asymptotic expressions for the short-time tail of the extinction- or blowup-time distribution become available.
  • The pre-exponential factor in the asymptotic form can be calculated explicitly.
  • The method is verified to work on three reaction systems that are solvable exactly without approximations.
  • Both extinction and blowup events in well-mixed particle systems are covered by the same framework.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar WKB matching techniques could be applied to other stochastic processes with singular short-time statistics, such as in population genetics.
  • The approach might allow predictions for systems too complex for exact solutions by providing reliable short-time asymptotics.
  • If extended, it could inform models of epidemic outbreaks or chemical ignition where initial fluctuations matter.

Load-bearing premise

The WKB solution obtained from the Laplace-transformed backward master equation can be matched to an inner solution for not too large m such that the composite asymptotic accurately describes the short-time tail.

What would settle it

For one of the three solvable examples, compute the exact distribution numerically for small T and check if it agrees with the predicted WKB asymptotic form including the pre-factor.

Figures

Figures reproduced from arXiv: 2601.04924 by Baruch Meerson, Michael Assaf, Rotem Degany.

Figure 1
Figure 1. Figure 1: shows the universal extinction time probabil￾ity distribution P(T), obtained by a numerical inverse Laplace transform of R(s) from Eq. (9). Also shown are the large- and small-T asymptotics of P(T): Eqs. (10) and (12), respectively. Until now we have been discussing exact results for P(T) and asymptotics extracted from these exact results. Now we will show that the small-T tail of P(T) behav￾ior can be cap… view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: (a) Solutions of Eq. (6) for the annihilation: ex [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: The MTE [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: (a) Solutions of Eq. (33): Exact (solid line), [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: (a) Solutions of Eq. (45): exact (solid line), WKB [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
read the original abstract

We study the statistics of extinction and blowup times in well-mixed systems of stochastically reacting particles. We focus on the short-time tail, $T \to 0$, of the extinction- or blowup-time distribution $\mathcal{P}_m(T)$, where $m$ is the number of particles at $t=0$. This tail often exhibits an essential singularity at $T=0$, and we show that the singularity is captured by a time-dependent WKB (Wentzel-Kramers-Brillouin) approximation applied directly to the master equation. This approximation, however, leaves undetermined a large pre-exponential factor. We show how to calculate this factor by applying a leading- and a subleading-order WKB approximation to the Laplace-transformed backward master equation. Accurate asymptotic results can be obtained when this WKB solution can be matched to another approximate solution (the ``inner" solution), valid for not too large $m$. We demonstrate and verify this method on three examples of reactions which are also solvable without approximations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops an asymptotic method for the short-time tail (T→0) of the extinction- or blowup-time distribution P_m(T) in stochastic reaction networks. It applies a time-dependent WKB approximation directly to the master equation to capture the essential singularity at T=0; the undetermined prefactor is obtained from leading- and subleading-order WKB analysis of the Laplace-transformed backward master equation, which is then matched to an inner solution valid for not-too-large m. The procedure is demonstrated and verified on three exactly solvable reaction systems.

Significance. If the matching step can be placed on a firmer footing, the approach supplies a concrete route to short-time statistics for networks lacking closed-form solutions, extending standard WKB techniques to time-dependent master equations while providing explicit checks against exact benchmarks. The verification on three independent solvable models is a clear strength, as it confirms that the composite asymptotic reproduces both the exponential singularity and the prefactor in those cases.

major comments (2)
  1. [Description of the matching procedure (following the Laplace-transform WKB analysis)] The central claim that accurate asymptotics are obtained “when this WKB solution can be matched to another approximate solution (the inner solution)” rests on the existence of an overlap region whose validity is asserted only for the three solvable examples. No general criterion is given for when the WKB outer solution and the inner solution share a common range of m as T→0, nor is a remainder estimate supplied showing that the matched composite reproduces the exact short-time singularity to the claimed order for generic rates. This issue is load-bearing for the prefactor determination.
  2. [Verification on exactly solvable cases] The verification on the three solvable systems confirms reproduction of the known asymptotics, but the manuscript does not report a quantitative test of matching-point independence (e.g., variation of the composite prefactor when the matching m is varied within the putative overlap window). Such a test would directly address the skeptic’s concern about error control.
minor comments (2)
  1. The notation distinguishing the outer WKB solution, the inner solution, and the composite form should be introduced with explicit symbols and a brief statement of their respective domains of validity.
  2. A short table summarizing the three reaction networks, their exact short-time forms, and the WKB predictions would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We are grateful to the referee for the thorough review and constructive comments on our manuscript. The suggestions regarding the matching procedure and verification tests are helpful, and we will incorporate clarifications and additional analysis in a revised version. We address each major comment below.

read point-by-point responses

  1. Referee: [Description of the matching procedure (following the Laplace-transform WKB analysis)] The central claim that accurate asymptotics are obtained “when this WKB solution can be matched to another approximate solution (the inner solution)” rests on the existence of an overlap region whose validity is asserted only for the three solvable examples. No general criterion is given for when the WKB outer solution and the inner solution share a common range of m as T→0, nor is a remainder estimate supplied showing that the matched composite reproduces the exact short-time singularity to the claimed order for generic rates. This issue is load-bearing for the prefactor determination.

    Authors: We agree that a more explicit justification of the overlap region would strengthen the presentation. In the revised manuscript we will add a dedicated paragraph deriving the scaling of the overlap window as T→0 from the WKB eikonal and prefactor equations; this scaling is independent of the specific form of the rates provided the WKB ansatz remains valid. We will also state the conditions on the reaction rates under which the overlap is guaranteed to exist for sufficiently small T. A complete rigorous remainder estimate for completely arbitrary rates lies beyond the scope of the present work, but we will explicitly note this limitation and emphasize that the method is intended for networks where the inner solution can be constructed. revision: partial

  2. Referee: [Verification on exactly solvable cases] The verification on the three solvable systems confirms reproduction of the known asymptotics, but the manuscript does not report a quantitative test of matching-point independence (e.g., variation of the composite prefactor when the matching m is varied within the putative overlap window). Such a test would directly address the skeptic’s concern about error control.

    Authors: We thank the referee for this concrete suggestion. In the revised manuscript we will include, for each of the three exactly solvable models, a quantitative check in which the matching value of m is varied across the overlap window. We will report the resulting variation in the extracted prefactor (typically a few percent or less) and include a brief table or figure demonstrating this stability. This addition directly addresses error control and the robustness of the composite asymptotic. revision: yes

Circularity Check

0 steps flagged

No significant circularity; standard WKB derivation verified externally

full rationale

The paper applies time-dependent WKB directly to the master equation for the short-time tail of P_m(T), then uses leading/subleading WKB on the Laplace-transformed backward master equation to determine the pre-factor, followed by matching to an inner solution. This is demonstrated and verified on three exactly solvable reaction networks where exact results exist independently of the approximation. No step reduces the target singularity or pre-factor to a fitted parameter by construction, nor relies on self-citation chains for load-bearing justification. The verification on independent exact solutions provides external anchoring, keeping the derivation self-contained against benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard validity assumptions of the WKB approximation in the short-time limit and the existence of a matching region between outer and inner solutions; no free parameters are introduced and no new entities are postulated.

axioms (2)
  • domain assumption WKB approximation is valid for capturing the leading exponential behavior of the short-time tail of the extinction/blowup time distribution
    Invoked when applying the time-dependent WKB directly to the master equation
  • domain assumption Leading- and subleading-order WKB on the Laplace-transformed backward master equation yields the correct prefactor
    Required to determine the undetermined multiplicative factor left by the forward WKB

pith-pipeline@v0.9.0 · 5482 in / 1471 out tokens · 65441 ms · 2026-05-16T16:14:52.367836+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

17 extracted references · 17 canonical work pages

  1. [1]

    The solution is R(s, m) = 2−mm! 3−√1−8s 4 m/2 3+√1−8s 4 m/2 ,(7) where (a)k ≡Γ(a+k)/Γ(a) is the Pochhammer symbol, and we continue to assume thatmis an even number. In particular, we find R(s,2) = 1 s+ 1 , R(s,4) = 6 (s+ 1)(s+ 6) , R(s,6) = 90 (s+ 1)(s+ 6)(s+ 15) , R(s,8) = 2520 (s+ 1)(s+ 6)(s+ 15)(s+ 28) .(8) 3 In general, an increase ofmby 2 adds one po...

    work page

  2. [2]

    universal

    scales as lnz, therefore ¯Tm goes down as lnµ/µfor largeµ.) FIG. 3: The MTE ¯Tvsµfor the initial number of par- ticlesm= 1 (dashed red line),m= 10 (dashed blue line) andm= 100 (dashed green line). Also shown is the MTE in the universal limitm→ ∞(solid black line). Applying the Laplace transform to the backward mas- ter equation (29), we arrive at the equa...

    work page

  3. [3]

    Delbr¨ uck, J

    M. Delbr¨ uck, J. Chem. Phys.120, 120 (1940)

    work page 1940

  4. [4]

    M. S. Bartlett,Stochastic Population Models in Ecology and Epidemiology(Wiley, New York, 1960)

    work page 1960

  5. [5]

    Karlin and J

    S. Karlin and J. McGregor,Stochastic Models in Medicine and Biology(University of Wisconsin Press, Madison, WI, 1964)

    work page 1964

  6. [6]

    N. T. J. Bailey,The Mathematical Theory of Infectious Diseases and its Applications(Griffin, London, 1975)

    work page 1975

  7. [7]

    Ovaskainen and B

    O. Ovaskainen and B. Meerson, Stochastic models of pop- ulation extinction, Trends in Ecology & Evolution25, 643 (2010)

    work page 2010

  8. [8]

    Allen,Stochastic Population and Epidemic Models — Persistence and Extinction(Springer, Berlin, 2015)

    J. Allen,Stochastic Population and Epidemic Models — Persistence and Extinction(Springer, Berlin, 2015)

    work page 2015

  9. [9]

    Johanssen and D

    A. Johanssen and D. Sornette, Physica A294, 465 (2001)

    work page 2001

  10. [10]

    Assaf and B

    M. Assaf and B. Meerson, Extinction of metastable stochastic populations, Phys. Rev. E81, 021116 (2010)

    work page 2010

  11. [11]

    Assaf and B

    M. Assaf and B. Meerson, WKB theory of large devia- tions in stochastic populations, J. Phys. A: Math. Theor. 50, 263001 (2017)

    work page 2017

  12. [12]

    Elgart and A

    V. Elgart and A. Kamenev, Rare event statistics in reaction-diffusion systems, Phys. Rev. E70, 041106 (2004)

    work page 2004

  13. [13]

    Meerson, Short-time blowup statistics of a Brown- ian particle in repulsive potentialse, Phys

    B. Meerson, Short-time blowup statistics of a Brown- ian particle in repulsive potentialse, Phys. Rev. E112, 064110 (2025)

    work page 2025

  14. [14]

    D. A. McQuarrie, C. J. Jachimowski, and J. Russell, Ki- netics of Small Systems. II, J. Chem. Phys.40, 2914 (1964)

    work page 1964

  15. [15]

    Assaf and B

    M. Assaf and B. Meerson, Spectral formulation and WKB approximation for rare-event statistics in reaction systems, Phys. Rev. E74, 041115 (2006)

    work page 2006

  16. [16]

    C. W. Gardiner,Handbook of Stochastic Methods (Springer, Berlin, 2004) p. 259

    work page 2004

  17. [17]

    Meerson, Fluctuations of blowup time in a simple model of a super-malthusian catastrophe, J

    B. Meerson, Fluctuations of blowup time in a simple model of a super-malthusian catastrophe, J. Stat. Mech. (2025) 053201. Appendix A: Time-Dependent WKB Theory for Sec.III Here we employ the leading-order time-dependent WKB method, as applied directly to the master equa- tion, in the case of binary coalescence and linear decay, 9 2A→AandA→0, that we also...

    work page 2025