Pith. sign in

REVIEW 4 major objections 3 minor 14 references

Integral representation of probabilities in Kingman coalescent

T0 review · 4 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read This paper claims that the block-counting process of the Kingman coalescent satisfies a small-time local central limit theorem: for $n=\lfloor(2+\sqrt{t}\,v)/t\rfloor$, $d_\theta^n(t)\sim \frac{3}{2\sqrt{2\pi}}\sqrt{t}\,e^{-3v^3/4}$…

desk verdict The integral representation idea is worth a look, but the main local CLT is not proved: Lemma 3.3's quadratic coefficient is wrong and the exponents are inconsistent. read the letter →

arxiv 1908.01474 v2 pith:3AV6SS3Z submitted 2019-08-05 math.PR

Youzhou Zhou This is my paper · ORCID
classification math.PR MSC 60F1060C05
keywords Kingmancoalescentblock-countingprocessfinite-timedistributionintegralrepresentationlocalcentrallimittheoremsteepestdescentsmall-timeasymptoticsJacobithetafunction
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's aim is to convert the unwieldy explicit distribution of the Kingman coalescent into a form that supports asymptotic analysis. The author starts from the known infinite-series formula for $d_\theta^n(t)=P(D_t=n)$, rewrites the combinatorial coefficients and exponential factors as contour integrals, and sums the resulting geometric series to get a clean integral representation. A second, deformed representation then feeds a steepest-descent argument whose payoff is a local central limit theorem in the small-time regime: the block count $D_t$ fluctuates around $2/t$ on the $\sqrt{t}$ scale with a Gaussian-type law. If the argument is correct, it is the first local CLT for this process at small times, and the same integral representation is proposed as a route to large and moderate deviations.

What carries the argument

The argument hangs on the phase function $\psi_t(z)=-2z^2-(2n+\theta)t\log\cos z$, which appears as $\exp(\psi_t((z+\pi)/2)/t)\sin((z+\pi)/2)\varphi_t(z)$ in the integrand. The real part of $\psi_t$ decreases monotonically along the deformed contour $L_1$ away from the critical point $z=-\pi$, so the integral localizes near $-\pi$; the Jacobi triple product identity converts the $\theta$ series into the explicit infinite product $\varphi_t(z)$, and the oddness of the full integrand allows the two contour pieces to be combined. On the short segment $L_{1,t}$ the integral is reduced to a Fourier-type integral in a variable $y$, which produces the claimed $\sqrt t$ scale and exponential factor.

What would settle it

Expand $\frac1t\psi_t((z(y)+\pi)/2)$ at $z(y)=-\pi+t^{1/4}e^{i\pi/4}\sqrt y$ using $\psi_t(z)=(-2+(2n+\theta)t/2)z^2+\cdots$; compare the coefficient of $y^2$ with the paper's $-(2n+\theta)/(12\cdot24)$. The direct expansion gives $-(2n+\theta)/192$, which on the scale $(2n+\theta)t\to4$ is $-1/48$ rather than $-1/72$. Evaluating the resulting Gaussian integral and comparing it with the original series for small $t$ and fixed $v$ would show whether the stated constant $3/(2\sqrt{2\pi})$ is correct.

Watch

Extended reading notes

Core claim

The central claim is that a contour-integral representation makes the finite-time distribution of the Kingman coalescent analytically accessible. Theorem 2.1 writes $d_\theta^n(t)$ as a one-dimensional complex integral whose integrand is $e^{-w^2/2t}\,i\sin(w/2)/(\cos(w/2))^{2n+\theta}$ times an explicit prefactor. Theorem 3.1 deforms the integration contour to a unit circle and uses the Jacobi triple product to absorb the $\theta$ sum into an infinite product $\varphi_t(z)$. The main theorem then claims that for $n=\lfloor(2+\sqrt t v)/t\rfloor$, $d_\theta^n(t)\sim \frac{3}{2\sqrt{2\pi}}\sqrt t\,e^{-3v^3/4}$ as $t\to0$, describing the fluctuation of the block count around $2/t$ on the $\sqrt t$ scale. The theorem as printed displays the exponent $-3v^3/4$, while the closing line of the proof writes $-3v^2/4$; this exponent discrepancy is one of the load-bearing points checked by the calculation in the falsifier below.

Load-bearing premise

The argument rests on the Taylor expansion of the exponent function along the deformed contour having exactly the quadratic coefficient stated in Lemma 3.3; a direct expansion of that same function appears to give a different coefficient, and if the direct expansion is right the final constant does not follow.

Editorial extensions

If this is right

  • If the theorem holds, the block count has the small-time fluctuation description $D_t\approx 2/t+\sqrt t\,Z$ with $Z$ distributed according to the limiting density given by the formula, refining the known deterministic limit $tD_t\to2$.
  • Because $d_\theta^n(t)$ are the convex coefficients in the transition functions of Fleming-Viot and infinitely-many-neutral-alleles diffusion models, the local CLT provides a concrete small-time approximation to those transition probabilities.
  • The same integral representation is proposed as a basis for large-deviation and moderate-deviation estimates for $D_t$ at small times, going beyond the central-limit scale.
  • The explicit prefactor in the theorem gives a quantitative benchmark that numerical evaluations or other asymptotic methods should reproduce in the same regime.

Reading between the lines

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

  • If the quadratic-coefficient discrepancy is real, the local-CLT shape may survive but the variance and prefactor would change; checking the $y^2$ coefficient directly is the cheapest way to test the stated constant.
  • A corrected Gaussian factor would read as $\exp(-c v^2)$ for some constant $c$, not the cubic exponent printed in the theorem; a careful re-derivation of Lemma 3.3 would determine $c$.
  • The same contour machinery could be tried on $\Lambda$-coalescents with a Kingman component, where small-time fluctuations have so far been studied with probabilistic rather than analytic methods; the analytic route would be a new comparison.
  • Controlling the remainder terms in the uniform convergence step (3) more carefully should yield higher-order asymptotic expansions in powers of $\sqrt t$, extending the local CLT to a full small-time expansion.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 3 minor

Summary. The paper starts from Tavaré's infinite-series expression for the probabilities d_θ^n(t) = P(D_t = n) of the block-counting process of the Kingman coalescent and rewrites it as a contour integral (Theorem 2.1). It then deforms the contour and uses the Jacobi triple product to obtain a second representation (Theorem 3.1). Finally, a steepest-descent analysis of this representation is used to claim a small-time local central limit theorem (Theorem 3.2): for n = floor((2 + √t v)/t), d_θ^n(t) ∼ (3/(2√(2π))) √t e^{-3v^3/4}. The main novelty claimed is the integral representation suited to asymptotic analysis and the new local CLT.

Significance. The approach is appealing: an explicit integral representation could become a useful tool for small-time asymptotics, and the derivation is non-circular, starting from Tavaré's independent series rather than assuming the conclusion. The paper also contains concrete, falsifiable asymptotic predictions. However, the core asymptotic calculation contains arithmetic inconsistencies that currently invalidate the proof; the result as stated is not established. If the constants and exponents are corrected and the contour deformation is justified, the method may still yield a publishable local CLT.

major comments (4)
  1. [§3.2, Lemma 3.3, Eq. (2)] The quadratic coefficient in the steepest-descent expansion is miscalculated. Substituting w = (z(y)+π)/2 = 2^{-1} t^{1/4} e^{iπ/4} √y into ψ_t(w) = -2w^2 - (2n+θ)t log cos w and using log cos w = -w^2/2 - w^4/12 + O(w^6) gives (1/t)ψ_t(z+π/2) = ((2n+θ)t - 4)/(8√t) i y - (2n+θ)t/192 y^2 + O(√t log^6(1/t)). The printed coefficient -(2n+θ)t/(12·24) = -(2n+θ)t/288 differs from -(2n+θ)t/192, so the limiting Gaussian variance is -1/72 instead of -1/48. This changes the Gaussian integral and the final constant in Lemma 3.3. The displayed Bernoulli expansion in Eq. (2) also has the wrong sign for the z^4 term under the standard convention B_4 = -1/30: ψ_t should have coefficient +(2n+θ)t/12 for z^4, and the y^2 term in Lemma 3.3 is negative only because w^4 = -t y^2/16.
  2. [§3.2, Lemma 3.3 and Theorem 3.2] The theorem statement and its proof disagree on the exponential argument: Theorem 3.2 displays e^{-3v^3/4}, while the last line of the proof uses e^{-3v^2/4}. A probability approximation of the form e^{-3v^3/4} is invalid because it is unbounded as v → -∞; the correct Gaussian form must be e^{-3v^2/4}. This is not a purely cosmetic typo because the asymptotic statement as printed is not a probability density in v.
  3. [Proof of Theorem 3.2] Even after correcting the quadratic coefficient, the constant 3/(2√(2π)) does not follow from the stated Lemma 3.3 and Eq. (1). With the corrected expansion, Lemma 3.3 gives ∫_{C+}(K_t(z)-K_t(-z))φ_t(z)dz ∼ -i√(3π) √t e^{-3v^2/4}. Using Eq. (1), which has the factor i/√(2πt), and the Stirling estimate (2n-1+θ choose n)/2^{2n-1+θ} ∼ √t/√(2π), one obtains d_θ^n(t) ∼ (√3/(2√π)) √t e^{-3v^2/4}, not 3/(2√(2π))√t e^{-3v^2/4}. The intermediate line in the proof of Theorem 3.2 and the displayed 'Therefore' are arithmetically incompatible; the proof must be recomputed.
  4. [§3.2, contour deformation after Eq. (1)] The displayed identity ∫_{C+}(K_t(z)-K_t(-z))φ_t(z)dz = -∫_{L1}K_t(z)φ_t(z)dz + ∫_{L1}K_t(z)φ_t(z)dz has two identical integrals on the right-hand side, so as printed it asserts that the left-hand side is zero. The intended relation among L1, L2, and the map z ↦ -z is not stated precisely, and the assertion that contributions along the branch cuts cancel is not proved. A complete contour-deformation argument is needed before Lemma 3.3 can be applied.
minor comments (3)
  1. [Abstract and Introduction] Citation numbers are inconsistent between the abstract and the body: Kingman is cited as [7] in the abstract but [9] in the text, and Tavaré is cited as [12] in the abstract but [13] in the text. There are also typographical errors such as 'bench-mark' and 'ininitely'.
  2. [Lemma 3.2] The statement contains a redundant and confusing double equality: |∫_{L^c_{1,t}} K_t φ dz| = |∫_{L^c_{1,t}} K_t φ dz| ≤ ... . Only one integral should appear on the left-hand side.
  3. [Throughout] The phrase 'steep descent' should be 'steepest descent' throughout Sections 3.1 and 3.2, and the figures should be referred to at the points where the contours are first described.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the derivation is self-contained from Tavaré's explicit series; only a minor non-load-bearing self-citation appears.

full rationale

The paper's derivation chain is self-contained: it starts from Tavaré's explicit finite-time series, rewrites the combinatorial coefficients and exponential factors using Cauchy's integral formula and Fourier transforms, sums a geometric series, and then applies the Jacobi triple product to obtain the alternative integral representation in Theorem 3.1. The subsequent steepest-descent estimates in Lemmas 3.2 and 3.3 are carried out on that representation, and Theorem 3.2 follows by combining these estimates with Stirling's formula. No parameter is fitted to a subset of the probabilities being predicted, and no target asymptotic is assumed as an input. The only self-citation is reference [14], mentioned in the introduction as an application context for the probabilities; it carries no load in the proof. The skeptical reviewer's noted discrepancy in Lemma 3.3 — the quadratic coefficient appears as -(2n+θ)t/72 rather than the direct expansion value -(2n+θ)t/48, and the final exponent is printed as v^3 in the theorem but v^2 in the proof — is a potential mathematical error and a correctness risk, not a circularity: the proof's final constant is not obtained by assuming the theorem. No circular step can be exhibited, so the paper should not receive a high circularity score.

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

The paper introduces no new physical or mathematical entities. It relies on known results (Tavaré's formula, Jacobi triple product) and on an asserted but under-proved contour deformation. The main unproved assumption is the correctness of the steepest descent expansion in Lemma 3.3.

assumptions (4)
  • domain assumption Tavaré's explicit series for d_θ^n(t) is correct and is the starting point.
    The paper relies on the explicit finite-time distribution from Tavaré [13] without re-deriving it. This is an external known result.
  • standard math The Jacobi triple product identity (Lemma 3.1) holds for the relevant q and x.
    Used in Theorem 3.1 to rewrite the theta-type function h_t(z|θ) as an infinite product.
  • ad hoc to paper The contour deformation from the horizontal lines to circles around odd multiples of π is valid, with contributions on branch cuts canceling by symmetry.
    The paper asserts this without a detailed proof. It is load-bearing for Theorem 3.1 and is not a standard textbook result.
  • domain assumption The principal branch of (cos w/2)^{2n+θ} is used consistently and the integrand decays sufficiently at infinity on the chosen contours.
    Required for Theorem 2.1 and for the contour deformation; the paper only sketches the decay estimate.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Integral representation of probabilities in Kingman coalescent." pith.science (2026). https://pith.science/paper/3AV6SS3Z

@misc{pith2026190801474,
  author       = {Pith},
  title        = {Pith review of: Integral representation of probabilities in Kingman coalescent},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3AV6SS3Z}},
  note         = {Machine review of arXiv:1908.01474}
}
read the original abstract

Kingman Coalescent was first proposed by Kingman [7] in population genetics to describe population's genealogical structure. Now it becomes a bench-mark model for coalescent process. Extensive studies have been conducted on Kingman coalescent. In particular, its explicit finite time distribution was obtained by Tavar\'e [12]. However, very few people use this explicit distribution to do analysis for it is an intractable infinite series. In this article, we are going to establish a complex integral representation for the finite time distribution, then we use steepest descent method to analyze this integral representation to obtain local central limit theorem at small time regime.

Figures

Figures reproduced from arXiv: 1908.01474 by the authors.

Figure 1
Figure 1. Deformed Rectangular Contour Then 2 Z ∞−iA −∞−iA f(w)dw = Z ∞−iA −∞−iA f(w)dw + Z −∞+iA ∞+iA f(w)dw = Z L f(w)dw where contour L are two parallel horizontal lines. One can easily deform the contour L to contours in [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Steep Descent Contour = Z C+ gt(z)φt(z)dz + Z C+ gt(z)φt(z)dz =2 Z C+ gt(z)φt(z)dz where C+ is the upper half circle and it has counter clockwise orientation. Hence (1) c θ n (t) = Y∞ k=1 [1 − e −4kπ2 t ]e (1−θ) 2t 8 i √ 2πt Z C+ (Kt(z) − Kt(−z))φt(z)dz. Now we consider two paths ( [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

14 extracted references · 11 canonical work pages

  1. [1]

    Milton Abramowitz and Irene A. Stegun. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. Dover, New York, ninth dover printing, tenth gpo printing edition, 1964

    1964

  2. [2]

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

    work page 1999

  3. [3]

    Sociedade Brasileira de Matem´ atica, Rio de Janeiro, 2009

    Nathana¨ el Berestycki.Recent progress in coalescent theory, volume 16 of Ensaios Matem´ aticos [Mathemat- ical Surveys]. Sociedade Brasileira de Matem´ atica, Rio de Janeiro, 2009

    work page 2009

  4. [4]

    Some large deviations in Kingman’s coalescent

    Andrej Depperschmidt, Peter Pfaffelhuber, and Annika Scheuringer. Some large deviations in Kingman’s coalescent. Electron. Commun. Probab., 20, no. 7:1–14, 2015

    work page 2015

  5. [5]

    S. N. Ethier and R. C. Griffiths. The transition function of a Fleming-Viot process. Ann. Probab. , 21(3):1571–1590, 1993

    work page 1993

  6. [6]

    Steven N. Evans. Kingman’s coalescent as a random metric space. In Stochastic models (Ottawa, ON, 1998), volume 26 of CMS Conf. Proc., pages 105–114. Amer. Math. Soc., Providence, RI, 2000

    work page 1998

  7. [7]

    R. C. Griffiths. Asymptotic line-of-descent distributions. J. Math. Biol. , 21(1):67–75, 1984

    work page 1984

  8. [8]

    R. C. Griffiths. Coalescent lineage distributions. Adv. Appl. Prob., 38:405–429, 2006

    work page 2006

Show all 14 references
  1. [9]

    J. F. C. Kingman. The coalescent. Stochastic Process. Appl., 13(3):235–248, 1982

    1982

  2. [10]

    Diffusion limits at small times for Λ-coalescents with a Kingman com- ponent

    Vlada Limic and Anna Talarczyk. Diffusion limits at small times for Λ-coalescents with a Kingman com- ponent. Electron. J. Probab., 20,no. 45:1–20, 2015

    2015

  3. [11]

    Second-order asymptotics for the block counting process in a class of regularly varying Λ-coalescents

    Vlada Limic and Anna Talarczyk. Second-order asymptotics for the block counting process in a class of regularly varying Λ-coalescents. Ann. Probab., 43(3):1419–1455, 2015

    2015

  4. [12]

    Coalescents with multiple collisions

    Jim Pitman. Coalescents with multiple collisions. Ann. Probab., 27(4):1870–1902, 1999

    1902

  5. [13]

    Line-of-descent and genealogical processes, and their applications in population genetics models

    Simon Tavar´ e. Line-of-descent and genealogical processes, and their applications in population genetics models. Theoret. Population Biol., 26(2):119–164, 1984

    1984

  6. [14]

    Ergodic inequality of a two-parameter infinitely-many-alleles diffusion model

    Youzhou Zhou. Ergodic inequality of a two-parameter infinitely-many-alleles diffusion model. J. Appl. Probab., 52(1):238–246, 2015. Department of Mathematical Science, Xi’an Jiaotong-Liverpool University, 111 Renai Road, Suzhou, China 215 123 E-mail address : youzhou.zhou@xjtlu.edu.cn

    2015

Pith tools

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