pith. machine review for the scientific record. sign in

arxiv: 2512.17485 · v2 · submitted 2025-12-19 · 🧮 math.PR · stat.CO

Recognition: 2 theorem links

· Lean Theorem

Koenigs functions in the subcritical and critical Markov branching processes with Poisson probability reproduction of particles

Penka Mayster , Assen Tchorbadjieff
Authors on Pith no claims yet

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

classification 🧮 math.PR stat.CO
keywords Koenigs functionsMarkov branching processesPoisson distributionsubcritical branchingcritical branchingBell polynomialsEin functionKolmogorov equation
0
0 comments X

The pith

Koenigs functions yield explicit solutions for subcritical and critical Poisson branching processes via Bell polynomials and Ein(z).

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

The paper derives explicit forms for Koenigs functions in Markov branching processes where particle reproduction follows a Poisson law in continuous time. These forms are obtained by integrating the Kolmogorov equation and involve exponential Bell polynomials together with the modified exponential-integral function Ein(z). The solutions are presented for both subcritical and critical branching mechanisms. They lead to expressions for the limit conditional law in the subcritical regime and the invariant measure in the critical regime. Such closed-form results connect iteration theory of functions with probabilistic models of population dynamics.

Core claim

For subcritical and critical continuous-time Markov branching processes with Poisson offspring distribution, the associated Koenigs functions admit explicit representations in terms of exponential Bell polynomials and the function Ein(z). These representations follow from direct integration of the Kolmogorov backward equation, supplemented by graphical depictions of the iteration. Consequently, the limiting conditional distribution conditional on non-extinction can be written down for the subcritical case, while an invariant probability measure is obtained for the critical case.

What carries the argument

Koenigs functions, the iterative limits of normalized probability generating function iterates that encode the asymptotic behavior of the branching process.

Load-bearing premise

Particle reproduction occurs according to a Poisson probability distribution in continuous time under either subcritical or critical branching.

What would settle it

Direct substitution of the proposed explicit formula into the Kolmogorov equation for a chosen Poisson rate that fails to satisfy the functional equation would falsify the claimed solutions.

Figures

Figures reproduced from arXiv: 2512.17485 by Assen Tchorbadjieff, Penka Mayster.

Figure 1
Figure 1. Figure 1: Graphics of h(s) for sub-, critical and supper- critical branching processes. The subcritical BP of [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Graphics of the function B(s) using eq. 8 for λ = 1/6 (black), λ = 5/6 (red) and λ = 17/18 (green). h (k) (0) = e −λλ k , k = 2, 3, ..., and the recurrence relation is specified as B (n+1)(0)e −λ = B (n) (0)e −λ {(λ + n − 1)e λ − nλ} −Xn k=2  n k  B (n+1−k) (0)e −λλ k . (8) We calculate the Taylor series expansion in the neighborhood of s = 0 with B(0) = 1 as B(s) = 1 +X∞ n=1 B (n) (0)s n n! , B′ (0) = (… view at source ↗
Figure 3
Figure 3. Figure 3: Histogram of fn = −B(n) (0) n! for different λ. For the subcritical Poisson reproduction law, 0 < λ < 1, the series expansion in the neighbourhood the point s = 1 h(s) = e λ(s−1) = X∞ k=0 λ k k! (s − 1)k , h(k) (1) = λ k . We start with B(1) = 0. Then, using the representation introduced in the following section B(s) = (1 − s)e G(s) , B′ (s) = −e G(s) + (1 − s)G ′ (s), we define B ′ (1) = −e G(1) < 0, G(1)… view at source ↗
Figure 4
Figure 4. Figure 4: Graphics of A(s) (left) and B(s) = e A(s) (right). −(λ − 1)e λX∞ k=0 (−1)kX k j=0  k j  (−1)k−j X k−j i=1  e λ λj i [k − j]i↓. The increasing and decreasing factorials, are defined by [x]k↓ = Γ(x + 1) Γ(x + 1 − k) , [x]k↑ = Γ(x + k) Γ(x) . Proof. In its explicit form, the function p(x) = e λx − (1 + xeλ ), (p(x))k = X k j=0 e λjx k j  (−1)k−j (1 + xeλ ) k−j . Then the integral log(B(s)) = (λ − 1)e λX… view at source ↗
Figure 5
Figure 5. Figure 5: Graphics of B(s) = (1 − s)e G(s) for λ = 1/3(blue), λ = 1/2(red) and λ = 2/3(green). = (1 − λ) log(1 − s) + (1 − λ) X∞ n=1 Z s x=0 (e λ(x−1) − 1)ndx (x − 1)n+1 . The indefinite integral for n = 1, 2, ... and Isaac Newton binomial expansion Z (e λ(x−1) − 1)ndx (x − 1)n+1 = (−1)n Z dx (x − 1)n+1 + Xn k=1  n k  (−1)n−k Z e λk(x−1)dx (x − 1)n+1 = (−1)n (−n)(x − 1)n + Xn k=1  n k  (−1)n−k    Xn j=1 −(λk)… view at source ↗
Figure 6
Figure 6. Figure 6: Graphics of G(s) (left) and B(s) = (1 − s)e G(s) (right). In particular, for n = 1 Z s x=0 g(x)dx (x − 1) = Z s x=0 e λ(x−1) − 1 − λ(x − 1) (1 − λ)(x − 1)2 ! dx = [PITH_FULL_IMAGE:figures/full_fig_p015_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Theorem 4.1. The explicit form of the function U(s) for λ = 1 is given by integration of the series expansion, U(s) = e X∞ k=0 (−1)kX k j=0  k j  e js(−1)k−j X k−j i=1  e j i (1 + se) k−j−i [k − j]i↓ −e X∞ k=0 (−1)kX k j=0  k j  (−1)k−j X k−j i=1  e j i [k − j]i↓. Compare with the geometric branching critical case, [14], and [13], m = 1 and 0 < s < 1. K f(x) = 1 (x − 1)2 + 1 1 − x , U(s) = − log(1 … view at source ↗
Figure 7
Figure 7. Figure 7: Graphics of Uλ(s) and U(s), U(0) = 0, U′ (0) = e (left) and respectively Cλ(s) = e Uλ(s) and C(s) = e U(s) , C(0) = 1 (right). The dashed line on right figure points out the minimal value at C(0) = 1. 4.2 In the neighborhood of the point s = 1 Consider the following decomposition in the neighborhood of s = 1, f(s) = K X∞ n=0 (s − 1)n n! − s ! = K [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
read the original abstract

Special functions have always played a central role in physics and in mathematics, arising as solutions of nonlinear differential equations, as well as in the theory of branching processes, which extensively uses probability generating functions. The theory of iteration of real functions leads to limit theorems for the discrete-time and real-time Markov branching processes. The Poisson reproduction of particles in real time is analysed through the integration of the Kolmogorov equation. These results are further extended by employing graphical representations of Koenigs functions under subcritical and critical branching mechanisms. The limit conditional law in the subcritical case and the invariant measure for the critical case are discussed, as well. The obtained explicit solutions contain the exponential Bell polynomials and the modified exponential-integral function $\rm{Ein} (z)$.

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

0 major / 2 minor

Summary. The manuscript derives explicit expressions for the Koenigs functions of subcritical and critical continuous-time Markov branching processes with Poisson offspring distribution. Integration of the Kolmogorov backward equation yields closed forms involving exponential Bell polynomials and the modified exponential-integral function Ein(z); substitution recovers the defining functional equation, from which the subcritical conditional limit law and critical invariant measure are obtained.

Significance. If the derivations hold, the explicit special-function representations connect standard branching-process iteration theory to Bell polynomials and Ein(z), offering concrete tools for further analytic work on limits and measures within the Poisson case. The approach remains inside the classical analytic framework for continuous-time branching processes.

minor comments (2)
  1. [Abstract] The abstract states that results are 'further extended by employing graphical representations of Koenigs functions' but provides no description of the figures, their construction, or the specific insights they convey; this should be clarified in the main text or figure captions.
  2. [Section on explicit solutions] The definition and normalization of the modified exponential-integral Ein(z) should be stated explicitly when first introduced, together with the precise relation to the exponential integral Ei or E_1 that is used in the derivation.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the positive assessment, including the recommendation for minor revision. The referee's summary correctly captures our derivation of explicit Koenigs functions via integration of the Kolmogorov backward equation, the use of exponential Bell polynomials and the modified exponential-integral function Ein(z), and the resulting limit conditional law and invariant measure.

Circularity Check

0 steps flagged

No significant circularity; derivations are self-contained within standard branching-process theory

full rationale

The manuscript integrates the Kolmogorov backward equation under a Poisson probability generating function, obtains explicit Koenigs functions expressed via exponential Bell polynomials and the modified exponential-integral Ein(z), and recovers the defining functional equation by direct substitution. These steps rely on classical iteration theory for real functions and the standard Kolmogorov forward/backward equations; no parameter is fitted to a data subset and then renamed as a prediction, no uniqueness theorem is imported from the authors' prior work, and no ansatz is smuggled via self-citation. The limiting conditional law and invariant measure follow directly from the same closed-form expressions without reduction to the inputs by construction. The derivation chain therefore remains independent of the target results.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper relies on standard axioms of continuous-time Markov branching processes and the specific Poisson reproduction assumption; no free parameters or invented entities are indicated in the abstract.

axioms (2)
  • standard math Kolmogorov forward equation governs the probability generating function of the branching process
    Standard differential equation for continuous-time Markov chains applied to branching processes.
  • domain assumption Reproduction of particles follows a Poisson probability distribution
    Central modeling choice for the offspring distribution in real time.

pith-pipeline@v0.9.0 · 5425 in / 1192 out tokens · 39513 ms · 2026-05-16T21:00:46.150491+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

24 extracted references · 24 canonical work pages

  1. [1]

    K. B. Athreya and P. E. Ney,Branching Processes, Springer, New York, 1972

    work page 1972

  2. [2]

    V . I. Dorogov and V . P. Chistyakov,Probabilistic Models of the Transformation of Particles, Nauka, Moscow, 1988 (in Russian)

    work page 1988

  3. [3]

    T. E. Harris,The Theory of Branching Processes, Springer, Berlin, 1963

    work page 1963

  4. [4]

    Pazsit and L

    I. Pazsit and L. P ´al,Neutron Fluctuations: A Treatise on the Physics of Branching Processes, Elsevier Science B.V ., Amsterdam, 2008

    work page 2008

  5. [5]

    B. A. Sevasyanov,Branching Processes, Nauka, Moscow, 1971 (in Russian; German transl.Verzwei- gungsprozesse, R. Oldenbourg Verlag, Munich, 1975)

    work page 1971

  6. [6]

    Schr ¨oder, ¨Uber iterierte Funktionen,Math

    E. Schr ¨oder, ¨Uber iterierte Funktionen,Math. Ann.3(1871), 296–322

    work page

  7. [7]

    M. Elin, V . Goryainov, S. Reich, and D. Shoikhet, Fractional iteration and functional equations for functions analytic in the unit disk,Comput. Methods Funct. Theory2(2002), no. 2, 353–366

    work page 2002

  8. [8]

    Karlin and J

    S. Karlin and J. McGregor, Embedding iterates of analytic functions with two fixed points into continuous semigroups,Trans. Amer. Math. Soc.132(1968), 137–145

    work page 1968

  9. [9]

    Koenigs, Recherches sur les int ´egrales de certaines ´equations fonctionnelles,Ann

    G. Koenigs, Recherches sur les int ´egrales de certaines ´equations fonctionnelles,Ann. Sci. ´Ecole Norm. Sup. (3) 1(1884), 2–41

    work page

  10. [10]

    V . V . Goryainov, Semigroups of analytic functions in analysis and applications,Russian Math. Surveys67 (2012), no. 6, 975–1021

    work page 2012

  11. [11]

    Kuczma, B

    M. Kuczma, B. Choczewski, and R. Ger,Iterative Functional Equations, Encyclopedia of Mathematics and its Applications, vol. 32, Cambridge Univ. Press, 1990

    work page 1990

  12. [12]

    Pakes, On the stationary measure for Markov branching processes,Mathematics13(2025), no

    A. Pakes, On the stationary measure for Markov branching processes,Mathematics13(2025), no. 11, 1802

    work page 2025

  13. [13]

    Tchorbadjieff and P

    A. Tchorbadjieff and P. Mayster, Geometric branching reproduction Markov processes,Modern Stochastics: Theory and Applications7(2020), no. 4, 357–378

    work page 2020

  14. [14]

    Tchorbadjieff and P

    A. Tchorbadjieff and P. Mayster, Factorial moments of the critical Markov branching process with geometric reproduction of particles,Modern Stochastics: Theory and Applications9(2022), no. 2, 229–244

    work page 2022

  15. [15]

    Tchorbadjieff and P

    A. Tchorbadjieff and P. Mayster, Wright function in the solution to the Kolmogorov equation of the Markov branching process with geometric reproduction of particles,Lithuanian Math. J.63(2023), 223–240

    work page 2023

  16. [16]

    Wang and T

    W. Wang and T. Wang, General identities on Bell polynomials,Comput. Math. Appl.58(2009), no. 1, 104–118

    work page 2009

  17. [17]

    Qi, D.-W

    F. Qi, D.-W. Niu, D. Lim, and Y .-H. Yao, Special values of the Bell polynomials of the second kind for some sequences and functions,J. Math. Anal. Appl.491(2020), 124382

    work page 2020

  18. [18]

    E. M. Wright, Iteration of the exponential function,Quart. J. Math.18(1947), no. 1, 228–235

    work page 1947

  19. [19]

    Jagers,Branching Processes with Biological Applications, John Wiley & Sons, 1975

    P. Jagers,Branching Processes with Biological Applications, John Wiley & Sons, 1975. 19

    work page 1975

  20. [20]

    Mode,Stochastic Processes in Demography and Their Computer Implementation, Springer, 1985

    C. Mode,Stochastic Processes in Demography and Their Computer Implementation, Springer, 1985

    work page 1985

  21. [21]

    Kimmel and D

    M. Kimmel and D. Axelrod,Branching Processes in Biology, Springer, 2002

    work page 2002

  22. [22]

    Haccou, P

    P. Haccou, P. Jagers, and V . Vatutin,Branching Processes: Variation, Growth and Extinction of Populations, Cambridge Univ. Press, 2007

    work page 2007

  23. [23]

    J. T. Conway, Indefinite integrals involving the exponential integral function,Integral Transforms Spec. Funct. 33(2022), no. 1, 1–15

    work page 2022

  24. [24]

    Comtet,Advanced Combinatorics: The Art of Finite and Infinite Expansions, D

    L. Comtet,Advanced Combinatorics: The Art of Finite and Infinite Expansions, D. Reidel Publishing Com- pany, Dordrecht, 1974. 20

    work page 1974