pith. sign in

arxiv: 2602.03228 · v4 · submitted 2026-02-03 · 🧬 q-bio.PE

Asymptotic Behavior of Integral Projection Models via Genealogical Quantities

Ryo Oizumi , Kensaku Kinjo , Yuki Chino This is my paper

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

classification 🧬 q-bio.PE
keywords integral projection modelsFredholm operatorsleading eigenvaluegenealogical quantitiesstable distributionreproductive valueBell polynomials
0
0 comments X

The pith

A reference-point operator turns continuous-state population kernels into explicit genealogical series for the leading eigenvalue and eigenfunctions.

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

The paper shows how to extract the dominant eigenvalue and the associated stable distribution and reproductive value directly from the kernel of an integral projection model by inserting a simple rank-one reference-point operator. This operator produces convergent series expansions whose terms are iterated kernels that admit a closed combinatorial expression via partial Bell polynomials. The resulting representations interpret population growth as an aggregation of ancestry-weighted contributions across generations, and they recover standard demographic quantities such as reproduction numbers and generation intervals without discretizing the state space.

Core claim

Under a dominant spectral separation condition, the reference-point construction yields series representations of the stable distribution and reproductive value in terms of iterated kernels that converge at the spectral radius; these series organize the leading eigensystem as a genealogical aggregation across generations and supply an Euler-Lotka-type equation through reference-point moments, all expressed combinatorially via ordinary partial Bell polynomials.

What carries the argument

The reference-point operator, a rank-one perturbation at the kernel level that makes point evaluation well-posed and induces a Markov-chain-style decomposition whose iterates are summed with partial Bell polynomials.

If this is right

  • Demographic indicators such as type reproduction numbers, generation intervals, and expected generation numbers are obtained directly from the continuous kernel without discretization.
  • The stable age or size distribution and reproductive value are expressed as infinite sums whose coefficients are explicit combinatorial functions of the iterated kernels.
  • The characteristic equation for the growth rate becomes an Euler-Lotka equation written in terms of reference-point moments.
  • The same construction applies to age-structured McKendrick-type equations and other positive-kernel Fredholm operators.

Where Pith is reading between the lines

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

  • The genealogical series may offer a new route to sensitivity analysis by differentiating term-by-term with respect to kernel parameters.
  • Because the construction avoids restrictive Hilbert-Schmidt assumptions, it could be applied to kernels with singularities that arise in certain spatial or size-structured models.
  • The combinatorial form via Bell polynomials suggests a possible link to branching-process interpretations that might simplify multi-type extensions.

Load-bearing premise

The dominant spectral separation condition holds for the kernels considered, including Hilbert-Schmidt, Doeblin-type, and rank-one perturbations.

What would settle it

Construct a specific continuous-state kernel satisfying the listed regularity conditions and compute numerically whether the partial-Bell-polynomial series for the leading eigenvalue converges to the value obtained by direct discretization or power iteration.

read the original abstract

We study the dominant eigenstructure of positive-kernel Fredholm operators arising in multi-state structured population models, including integral projection models and age-structured McKendrick-type equations. To obtain a determinant-free and interpretable characterization of the leading eigenvalue and eigenfunctions, we introduce a reference-point operator, a rank-one construction at the kernel level that renders point evaluation well posed and induces a Markov-chain-inspired decomposition in the continuous-state setting. This yields convergent series representations of the stable distribution and reproductive value in terms of iterated kernels, together with an Euler-Lotka-type characteristic equation expressed through reference-point moments. The iterates admit a closed combinatorial form via ordinary partial Bell polynomials, providing an explicit bridge from transition kernels to genealogical quantities. Under a dominant spectral separation condition, satisfied for a broad class of kernels including Hilbert-Schmidt, Doeblin-type, and rank-one perturbations, the expansion converges at the spectral radius and organizes the leading eigensystem as a genealogical aggregation across generations. As applications, we derive demographic indicators-type reproduction numbers, generation intervals, and expected generation numbers-directly from continuous-state kernels, without discretization and without restrictive Hilbert-Schmidt assumptions. The resulting framework clarifies how ancestry-weighted initial-state information accumulates across generations to determine population growth and composition.

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 paper introduces a reference-point operator (a rank-one perturbation of the transition kernel) for positive-kernel Fredholm operators arising in integral projection models. This construction yields series expansions for the leading eigenfunctions and eigenvalue via iterated kernels, expressed in closed form using ordinary partial Bell polynomials, together with an Euler-Lotka-type equation in terms of reference-point moments. Under a dominant spectral separation condition (simple leading eigenvalue r with |λ| < r for the rest of the spectrum), the expansions converge at the spectral radius and interpret the eigenstructure as a genealogical aggregation; the framework is claimed to apply to Hilbert-Schmidt, Doeblin-type, and rank-one perturbation kernels and to deliver demographic indicators without discretization.

Significance. If the central derivations are complete, the work supplies a determinant-free, combinatorially explicit bridge from continuous-state kernels to genealogical quantities (reproduction numbers, generation intervals, expected generation numbers) that organizes ancestry-weighted information across generations. The use of Bell polynomials for the iterates is a concrete technical strength that could make the expansions directly computable from kernel moments.

major comments (2)
  1. [Abstract and main convergence theorem (likely §3–4)] The manuscript asserts that the dominant spectral separation condition holds for the broad class of positive Hilbert-Schmidt kernels (abstract and the statement of the main convergence theorem). No explicit proof or counter-example verification is supplied. Block-diagonal kernels with two disjoint positive-measure supports are positive and Hilbert-Schmidt yet possess at least two eigenvalues of modulus equal to the spectral radius, violating separation; the Bell-polynomial series and reference-point Euler-Lotka equation lose their asymptotic justification in such cases.
  2. [Definition of the reference-point operator and subsequent series construction] The reference-point operator is defined as a rank-one construction that renders point evaluation well-posed, but the precise normalization and choice of reference point are not shown to be independent of the leading eigenfunction itself; this risks circularity in the genealogical aggregation when the operator is applied to general (non-strictly positive) kernels.
minor comments (2)
  1. [Notation and §2] The notation for the iterated kernels K_n and the partial Bell polynomials B_n should be introduced with an explicit recursive definition or generating-function identity in the main text rather than only in an appendix.
  2. [Numerical examples] Figure captions for any numerical illustrations of the series convergence should state the precise kernel, the value of r, and the truncation order used.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and constructive comments, which have helped us strengthen the clarity and precision of our claims. We address each major comment point by point below.

read point-by-point responses

  1. Referee: [Abstract and main convergence theorem (likely §3–4)] The manuscript asserts that the dominant spectral separation condition holds for the broad class of positive Hilbert-Schmidt kernels (abstract and the statement of the main convergence theorem). No explicit proof or counter-example verification is supplied. Block-diagonal kernels with two disjoint positive-measure supports are positive and Hilbert-Schmidt yet possess at least two eigenvalues of modulus equal to the spectral radius, violating separation; the Bell-polynomial series and reference-point Euler-Lotka equation lose their asymptotic justification in such cases.

    Authors: We appreciate the referee pointing out this counterexample, which correctly shows that not every positive Hilbert-Schmidt kernel satisfies dominant spectral separation. The manuscript presents spectral separation as an explicit hypothesis under which the series converge, while noting that it is satisfied for many kernels in the listed classes (including irreducible or strictly positive Hilbert-Schmidt kernels). In the revision we will (i) rephrase the abstract and Theorem 3.1 to state the separation condition more prominently as an assumption rather than an automatic property of all positive Hilbert-Schmidt kernels, (ii) add a brief paragraph in §2 citing standard results (Krein–Rutman and related positivity theorems) that guarantee a simple dominant eigenvalue when the kernel is irreducible or positive on a connected support, and (iii) explicitly exclude the block-diagonal case from the “broad class” claim. These changes will make the scope of the convergence result precise without altering the core derivations. revision: yes

  2. Referee: [Definition of the reference-point operator and subsequent series construction] The reference-point operator is defined as a rank-one construction that renders point evaluation well-posed, but the precise normalization and choice of reference point are not shown to be independent of the leading eigenfunction itself; this risks circularity in the genealogical aggregation when the operator is applied to general (non-strictly positive) kernels.

    Authors: The reference point x₀ is chosen a priori as any fixed point in the state space (independent of the unknown eigenfunctions) and the rank-one perturbation is defined directly from the original kernel K(·,x₀) and K(x₀,·). The normalization is fixed by requiring that the perturbed operator reproduces the original spectral radius while making point evaluation continuous; this construction does not presuppose knowledge of the leading eigenfunction. In the revision we will (i) rewrite the definition in §2.2 to state explicitly that x₀ is selected before any spectral computation, (ii) add a short lemma showing that the resulting Bell-polynomial series and Euler–Lotka equation are invariant under different admissible choices of x₀, and (iii) clarify that, for kernels that are not strictly positive, the genealogical interpretation still holds under the maintained spectral-separation hypothesis. These additions remove any appearance of circularity. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivations proceed from introduced reference-point operator

full rationale

The paper introduces a new reference-point operator (rank-one construction) and derives series expansions for the stable distribution and reproductive value using iterated kernels and Bell polynomials. The dominant spectral separation condition is posited as an assumption for the class of kernels considered, but the leading eigenvalue and eigenfunctions are not defined in terms of themselves or fitted to the target quantities. No load-bearing self-citation reduces the central claims to prior unverified results by the same authors, and the genealogical aggregation follows from the operator decomposition rather than by construction. This is consistent with a self-contained derivation against the transition kernel.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central construction rests on introducing one new operator and invoking standard spectral properties of positive kernels; no fitted parameters or additional invented entities beyond the reference-point operator are stated.

axioms (2)
  • domain assumption Positive kernels admit a dominant eigenvalue separated from the rest of the spectrum under the stated conditions (Hilbert-Schmidt, Doeblin, rank-one perturbations).
    Invoked to guarantee convergence of the genealogical series at the spectral radius.
  • standard math Fredholm operators on suitable function spaces possess the spectral properties needed for the reference-point construction.
    Background from functional analysis used without proof in the abstract.
invented entities (1)
  • reference-point operator no independent evidence
    purpose: Rank-one construction at the kernel level that renders point evaluation well-posed and induces a Markov-chain decomposition.
    New object introduced to obtain the series representations and Euler-Lotka equation.

pith-pipeline@v0.9.0 · 5524 in / 1340 out tokens · 36548 ms · 2026-05-16T08:00:59.589431+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

27 extracted references · 27 canonical work pages

  1. [1]

    Chung, K. L. (1960).Markov Chains with Stationary Transition Probabilities. Springer, New York

    work page 1960

  2. [2]

    (1974).Advanced Combinatorics: The Art of Finite and Infinite Expansions

    Comtet, L. (1974).Advanced Combinatorics: The Art of Finite and Infinite Expansions. D. Reidel Publishing Company, Dordrecht

    work page 1974

  3. [3]

    and Tuljapurkar, S

    Coulson, T. and Tuljapurkar, S. (2008). Dynamics of a variable-sized population in a stochastic environment.The American Naturalist, 172, 599–612

    work page 2008

  4. [4]

    Coulson, T., Tuljapurkar, S., and Childs, D. Z. (2010). Using evolutionary demography to link life history theory, quantitative genetics and population ecology.Journal of Animal Ecology, 79, 1226–1240

    work page 2010

  5. [5]

    Doak, D. F. et al. (2021). A critical comparison of integral projection and matrix projection models.Ecological Monographs, 91, e01447

    work page 2021

  6. [6]

    Elderd, B. D. and Miller, T. E. X. (2016). Quantifying variation in population growth. Ecological Monographs, 86, 125–144

    work page 2016

  7. [7]

    Ellner, S. P. and Rees, M. (2006). Integral projection models for populations in temporally varying environments.The American Naturalist, 167, 410–428

    work page 2006

  8. [8]

    Feller, W. (1941). The integral equation of renewal theory.Annals of Mathematical Statistics, 12, 243–267

    work page 1941

  9. [9]

    Fredholm, I. (1903). Sur une classe d’´ equations fonctionnelles.Acta Mathematica, 27, 365–390

    work page 1903

  10. [10]

    Krein, M. G. and Rutman, M. A. (1948). Linear operators leaving invariant a cone in a Banach space.American Mathematical Society Translations

    work page 1948

  11. [11]

    Heesterbeek, J. A. P. and Roberts, M. G. (2007). The type-reproduction number in structured epidemic models.Mathematical Biosciences, 206, 3–10

    work page 2007

  12. [12]

    (1981).Geometric Theory of Semilinear Parabolic Equations

    Henry, D. (1981).Geometric Theory of Semilinear Parabolic Equations. Springer

    work page 1981

  13. [13]

    and Herv´ e, L

    Hennion, H. and Herv´ e, L. (2001).Limit Theorems for Markov Chains and Stochastic Properties of Dynamical Systems. Springer

    work page 2001

  14. [14]

    Inaba, H. (2013). On the definition and computation of the type-reproduction number. Journal of Mathematical Biology, 66, 1065–1097

    work page 2013

  15. [15]

    (2017).Age-Structured Population Dynamics in Demography and Epidemiology

    Inaba, H. (2017).Age-Structured Population Dynamics in Demography and Epidemiology. Springer

    work page 2017

  16. [16]

    McKendrick, A. G. (1925). Applications of mathematics to medical problems.Proceedings of the Edinburgh Mathematical Society, 44, 98–130

    work page 1925

  17. [17]

    Merow, C. et al. (2014). Advancing population ecology with integral projection models. Methods in Ecology and Evolution, 5, 99–110

    work page 2014

  18. [18]

    Mode, C. J. (1971).Multitype Branching Processes. American Elsevier

    work page 1971

  19. [19]

    (1877).Untersuchungen ¨ uber das logarithmische und Newton’sche Potential

    Neumann, C. (1877).Untersuchungen ¨ uber das logarithmische und Newton’sche Potential. Teubner. Nicol` e, F. et al. (2011). Interdependent effects of habitat quality and climate.Journal of Ecology, 99, 1211–1218

    work page 2011

  20. [20]

    (2022).Population Dynamics Based on Individual Stochasticity

    Oizumi, R. (2022).Population Dynamics Based on Individual Stochasticity. Springer. ASYMPTOTIC BEHAVIOR OF IPMS 35

    work page 2022

  21. [21]

    Oizumi, R., Inaba, H., Takada, T., Enatsu, Y., and Kinjo, K. (2022). Sensitivity analysis on declining population in Japan.PLOS ONE, 17, e0273817

    work page 2022

  22. [22]

    Oizumi, R., Kuniya, T., and Enatsu, Y. (2016). Reconsideration of eigenvalue problems. PLOS ONE, 11, e0157715

    work page 2016

  23. [23]

    and Takada, T

    Oizumi, R. and Takada, T. (2013). Optimal harvesting strategies.Journal of Theoretical Biology, 323, 76–89

    work page 2013

  24. [24]

    (1958).An Introduction to Combinatorial Analysis

    Riordan, J. (1958).An Introduction to Combinatorial Analysis. John Wiley & Sons

    work page 1958

  25. [25]

    and Wenska, T

    Schumitzky, A. and Wenska, T. (1975). An operator approach to population dynamics. SIAM Journal on Mathematical Analysis, 6, 229–235

    work page 1975

  26. [26]

    Sharpe, F. R. and Lotka, A. J. (1911). A problem in age-distribution.Philosophical Magazine, 21, 435–438. Verduyn Lunel, S. M. (1989).Exponential Type Calculus for Linear Delay Equations. Verduyn Lunel, S. M. (1990). Series expansions for functional differential equations.Journal of Differential Equations, 85, 17–53. von Foerster, V. (1959). Some remarks ...

    work page 1911

  27. [27]

    White, J. W. et al. (2016). Fitting state-space integral projection models.Ecological Applications, 26, 2677–2694. Email address:ooizumi-ryou@ipss.go.jp Email address:ft41140@ns.kogakuin.ac.jp Email address:chino@nycu.edu.tw

    work page 2016