REVIEW 3 major objections 2 minor 24 references
The distribution function for the maximal height of $N$ non-intersecting Bessel paths
T0 review · 3 major / 2 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read For N non-intersecting Bessel paths, the probability that the outermost path stays below a level M is given exactly by a determinant series over the zeros of the Bessel function J_α.
desk verdict The formulas are likely new and may even be true, but the proof of the main theorem rests on a false determinant lemma, so the paper needs major revision before it is publishable. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The engine of the proof is the Karlin-McGregor ratio P = lim_{x→a, y→0} q_M(x,y)/q(x,y), where q and q_M are determinants of one-particle transition densities and q_M is built from the eigenfunction expansion of the Bessel generator on (0,M) with an absorbing wall at M; the eigenfunctions are Bessel functions evaluated at the zeros x_{n,α}. The asymptotic evaluation of the ratio when all starting points coalesce to a and all ending points to 0 is carried out with a determinant identity (Lemma 1) and a Schur-function expansion (Lemma 2) that extract the leading powers of Δ(x−a)Δ(y²) or Δ(x²)Δ(y²). The conversion to Hankel determinants is done by row operations using Bessel recurrence relations (Lemma 3) that replace the higher derivatives $J_α^{{(i−1)}}$(a x/M) with combinations of J_α and J_{α+1}, leaving the moment matrices of the discrete weights w_1, w_2, and w̃.
What would settle it
A decisive check is to simulate two independent Bessel processes with α = 1 started at the same small a, conditioned to end at 0 and not to collide, and compare the empirical frequency of max < M with formula (2.2); a mismatch beyond sampling error would refute the formula, as would any numerical disagreement between the α = 1/2 reduction of (2.4) and the classical Brownian-excursion formula (1.2).
Extended reading notes
Core claim
On its own terms, the paper establishes Theorem 1: for α > −1, N ≥ 1, M ≥ a, and a > 0, the probability that max_{0<t<1} b_N(t) < M equals c_N(α) $M^{{−N(3N+2α+1)/2}}$ times the determinant of the matrix whose (i,j) entry is ∑_{n=1}^∞ (−1)^{i−1} x_{n,α}^{i+2j+α−3} $J_α^{{(i−1)}}$(a x_{n,α}/M) $e^{{−x_{n,α}}$^2/($2M^{2}$)} / J_{α+1}^2(x_{n,α}), where x_{n,α} are the positive zeros of J_α and c_N(α) is an explicit constant. For a = 0, the formula is the analogous determinant (2.4) with powers x_{n,α}^{2i+2j+2α−4}. Theorem 2 rewrites both probabilities as Hankel determinants of moments of the discrete weights w_1, w_2, and w̃, so the distribution is expressed through multiple discrete orthogonal polynomials or discrete orthogonal polynomials. The paper also shows that the α = ±1/2 cases reproduce the known non-intersecting Brownian wall formulas, and that the single-path case N = 1, a = 0 reproduces Pitman and Yor's law of the Bessel bridge maximum.
Load-bearing premise
The argument assumes that the Karlin-McGregor determinant ratio (3.3) remains valid for Bessel processes with all N paths starting at the same point and ending at the same point, even though the Bessel process has a singular boundary at the origin and the infinite series over Bessel zeros is interchanged with the coalescing limits.
Editorial extensions
If this is right
- For every finite N and every α > −1, the maximum-height distribution is now an explicit series, so probabilities and their derivatives can be evaluated numerically without simulating the conditioned paths.
- In the special cases α = 1/2 and α = −1/2, the formulas reduce to the known maximum distributions for non-intersecting Brownian excursions and reflected Brownian motions, so the Bessel model is an exact interpolation between those two classical wall models.
- The Hankel-determinant writing brings the problem into the range of steepest-descent and orthogonal-polynomial asymptotics; the authors expect, under suitable scaling, convergence to the Tracy-Widom distribution of the Gaussian orthogonal ensemble as N → ∞.
- The N = 1, a = 0 reduction agrees with Pitman and Yor's known closed form, which anchors the general formula to a single-particle result.
- The formulas respect the natural large-M limit, tending to 1 as the ceiling M grows, which is consistent with the probability interpretation of the expression.
Reading between the lines
- Editorial inference: the same q_M/q determinant-ratio method should extend to joint survival probabilities for several levels, because the spectral expansion in Bessel zeros is the only model-dependent input and the determinant identities do not use the specific level M.
- Editorial inference: since squared Bessel processes are squared norms of multidimensional Brownian paths, the formula should translate into a statement about the radial components of N non-colliding d-dimensional Brownian particles confined to a ball of radius M, possibly connecting to Wishart-type eigenvalue statistics.
- Editorial inference: a natural testable extension is to let both N and α scale together; the exact finite-N series could be compared against large-deviation asymptotics before the conjectured GOE limit sets in.
- Editorial inference: the two starting-point cases a > 0 and a = 0 have different orthogonal-polynomial structures, so the transition as a → 0 may reveal a discrete-weight phenomenon not visible in the Brownian limits.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper considers an ensemble of N non-intersecting Bessel paths of order α, all starting at a≥0 and conditioned to end at 0 at time 1. The authors state two theorems: Theorem 1 gives the probability that the maximum of the outermost path stays below M as a determinant involving a sum over zeros of the Bessel function J_α (Eqs. (2.2) and (2.4)); Theorem 2 rewrites these probabilities as Hankel determinants associated with discrete (multiple) orthogonal polynomials (Eqs. (2.16) and (2.17)). The derivation proceeds via the Karlin-McGregor formula, a spectral expansion for the transition density with an absorbing wall at M, and asymptotic evaluations of the resulting determinants. Special cases α=-1/2 and α=1/2 are shown to reproduce the known non-intersecting Brownian motion formulas, and the N=1 case reproduces a result of Pitman-Yor.
Significance. If the stated formulas are correct, they provide the first explicit finite-N closed forms for the maximum-height distribution of non-intersecting Bessel paths, and the connection to discrete multiple orthogonal polynomials is a useful bridge to Riemann-Hilbert asymptotics. The paper is self-contained in the sense that it introduces no fitted parameters, and the agreement of the α=±1/2 special cases with the known Brownian formulas is a strong nontrivial check. The determinant representations in Theorem 2 are concrete and potentially computable. However, the validity of the proof is the critical issue.
major comments (3)
- [Section 3.2, Lemma 2 (Eq. (3.17))] Lemma 2 is false as stated. For N=2, a=1, y=(1,2), and f(z)=z^2, the left-hand side det_{i,j}[f(x_i y_j)] is identically zero for all x, while the right-hand side equals 8(x_2-x_1), which is nonzero. The flaw is in the proof: after applying Lemma 1, the step 'we order the index k=(k_1,...,k_N) such that {k_j} is decreasing' discards the summation over permutations that assign derivative orders to columns; for generic f the correct leading term is Δ(x-a) det_{i,j}[ y_j^{i-1} f^{(i-1)}(a y_j)/(i-1)! ] plus higher-order terms, not a product. Since this lemma is used directly in Eq. (4.13) of Proposition 1 to evaluate det[J_α(x_{n_j}/M x_i)], and Proposition 1 is the input to the proof of Theorem 1 in §4.2, the derivation of the central formulas (2.2) and (2.4) is not valid as written. The final formulas may be repairable, but the proof must be reworked.
- [Section 4.1, Proposition 1 (Eqs. (4.3)–(4.7))] The asymptotic expansions (4.3)-(4.7) are obtained by taking the limits x→a and y→0 inside the infinite sums over Bessel zeros and inside the Taylor-series representation of the modified Bessel function I_α, without a uniform-convergence or dominated-convergence justification. For example, in the derivation of (4.4), the series in (4.2) over n∈N^N is interchanged with the limit x→a and y→0; the terms contain oscillatory factors J_α(x_{n,α}/M ·) and the exponential e^{-x_{n,α}^2/(2M^2)}, and near y=0 the behavior J_α(x_{n,α} y/M) ∼ (x_{n,α} y/(2M))^α/Γ(α+1) makes the convergence non-uniform in y for α<0. This gap is load-bearing because Proposition 1 provides the main asymptotic estimates that feed Theorem 1.
- [Section 3.1, Eq. (3.3)] The Karlin-McGregor formula is invoked for the coalescing limit x→a, y→0, but the standard form of the formula requires distinct starting and ending points. The authors cite [13] for this limiting procedure, but no proof is given that the ratio of determinants converges to the non-intersecting probability in the presence of the singular boundary at 0. This is a gap in the derivation of the central identity (3.3); although the agreement of the final formulas with the known Brownian special cases is reassuring, it does not by itself supply the missing justification.
minor comments (2)
- [Section 3.2, proof of Lemma 2] The set N0 = {0,1,2,...} is used in the proof but is not defined before the lemma; please define it in the statement or in a preliminary remark.
- [Section 4.3, proof of Theorem 2] The row-reduction argument leading to (4.26) is described only in words; in particular, the parity properties of the polynomials P, Q, P~, Q~ in Lemma 3 are asserted but not explicitly used to justify the elimination of the J_α terms at each step. A more explicit induction would improve readability.
Circularity Check
No significant circularity: the maximum-height formula is derived self-containedly from the Bessel transition density, the Karlin-McGregor formula, and spectral theory; special cases are checked externally.
full rationale
The derivation starts with the Bessel transition density (1.9)-(1.10) and the absorbing-wall eigenfunction expansion (3.6)-(3.7), applies the Karlin-McGregor determinant ratio (3.3), and then evaluates the resulting determinants via Taylor expansion and Schur-function asymptotics. No parameter is fitted to the target probability, and no defining assumption is equivalent to the conclusion. The special cases a=0, N=1, and alpha=+/-1/2 are checked against the independent Brownian formulas (1.2)/(1.5) and Pitman-Yor [21], which provides external validation. The citation to Katori-Tanemura [13] for the Karlin-McGregor formula is an external theorem, not a self-citation, and is not used to forbid alternatives. The alleged failure of Lemma 2, if true, would be a correctness gap in the proof, not circularity; the analysis here therefore does not raise the circularity score on that basis.
Assumptions & free parameters
assumptions (3)
- domain assumption Karlin-McGregor formula for non-colliding Bessel processes with absorbing wall at M and coalescing endpoints
- domain assumption Spectral expansion (3.7) for the transition density with absorbing wall at M
- standard math Zeros x_{n,alpha} of J_alpha are real, simple, positive and form an orthonormal basis with weight 2/(M^2 J_{alpha+1}^2(x_{n,alpha}))
Cite this review
Pith. "Pith review of The distribution function for the maximal height of $N$ non-intersecting Bessel paths." pith.science (2026). https://pith.science/paper/OI66QYDA
@misc{pith2026190800736,
author = {Pith},
title = {Pith review of: The distribution function for the maximal height of $N$ non-intersecting Bessel paths},
year = {2026},
howpublished = {\url{https://pith.science/paper/OI66QYDA}},
note = {Machine review of arXiv:1908.00736}
}
abstract
In this paper, we consider $N$ non-intersecting Bessel paths starting at $x=a\geq 0$, and conditioned to end at the origin $x=0$. We derive the explicit formula of the distribution function for the maximum height. Depending on the starting point $a>0$ or $a=0$, the distribution functions are also given in terms of the Hankel determinants associated with the multiple discrete orthogonal polynomials or discrete orthogonal polynomials, respectively.
Figures
Reference graph
Works this paper leans on
-
[13]
M. Katori and H. Tanemura, Noncolliding Brownian motion and determinantal pro- cesses, J. Stat. Phys. , 129 (2007), no. 5-6, 1233-1277
work page 2007
-
[1]
Abraham, Reflecting Brownian snake and a Neumann-Dirichlet problem, Stochastic Process
R. Abraham, Reflecting Brownian snake and a Neumann-Dirichlet problem, Stochastic Process. Appl., 89 (2000), no. 2, 239-260
work page 2000
-
[2]
J. Arves´ u, J. Coussement and W. Van Assche, Some discrete multiple orthogonal poly- nomials, J. Comput. Appl. Math. , 153 (2003), 19-45
work page 2003
-
[3]
A. N. Borodin and P. Salminen, Handbook of Brownian Motion-Facts and Formulae , 2nd ed, Birkh¨ auser Verlag, Basel, 2002
work page 2002
-
[4]
S. Delvaux, A. B. J. Kuijlaars, P. Rom´ an and L. Zhang, Non-intersecting squared Bessel paths with one positive starting and ending point, J. Anal. Math. , 118 (2012), no. 1, 105-159
work page 2012
-
[5]
Feierl, The height of watermelons with wall, J
T. Feierl, The height of watermelons with wall, J. Phys. A , 45 (2012), no. 9, 095003, 26 pp
work page 2012
-
[6]
P. J. Forrester, S. N. Majumdar and G. Schehr, Non-intersecting Brownian walkers and Yang-Mills theory on the sphere, Nuclear Phys. B , 844 (2011), no. 3, 500-526
work page 2011
-
[7]
A. G¨ oing-Jaeschke and M. Yor, A survey and some generalizations of Bessel processes, Bernoulli, 9 (2003), no. 2, 313-349. 19
work page 2003
Show all 24 references
-
[8]
M. E. H. Ismail, Classical and Quantum Orthogonal Polynomials in One Variable , En- cyclopedia of Mathematics and its Applications, vol. 98, Cambridge University Press, 2005
2005
-
[9]
Itˆ o and H
K. Itˆ o and H. P. McKean,Diffusion Processes and their Sample Paths , Springer-Verlag, Berlin-New York, 1974
1974
-
[10]
Karatzas and S
I. Karatzas and S. E. Shreve, Brownian Motion and Stochastic Calculus, 2nd ed, Springer-Verlag, New York, 1991
1991
-
[11]
Katori, M
M. Katori, M. Izumi and N. Kobayashi, Two Bessel bridges conditioned never to collide, double Dirichlet series, and Jacobi theta function, J. Stat. Phys., 131 (2008), 1067-1083
2008
-
[12]
Katori and H
M. Katori and H. Tanemura, Symmetry of matrix-valued stochastic processes and non- colliding diffusion particle systems, J. Math. Phys. , 45 (2004), no. 8, 3058-3085
2004
-
[14]
Kobayashi, M
N. Kobayashi, M. Izumi and M. Katori, Maximum distributions of bridges of noncol- liding Brownian paths, Phys. Rev. E , 78 (2008), no. 5, 051102, 15 pp
2008
-
[15]
A. B. J. Kuijlaars, A. Mart´ ınez-Finkelshtein and F. Wielonsky, Non-intersecting squared Bessel paths and multiple orthogonal polynomials for modified Bessel weights, Comm. Math. Phys., 286 (2009), no. 1, 217-275
2009
-
[16]
A. B. J. Kuijlaars, A. Mart´ ınez-Finkelshtein and F. Wielonsky, Non-intersecting squared Bessel paths: critical time and double scaling limit, Comm. Math. Phys. , 308 (2011), no. 1, 227-279
2011
-
[17]
Liechty, Nonintersecting Brownian motions on the half-line and discrete Gaussian orthogonal polynomials, J
K. Liechty, Nonintersecting Brownian motions on the half-line and discrete Gaussian orthogonal polynomials, J. Stat. Phys. , 147 (2012), no. 3, 582-622
2012
-
[18]
I. G. Macdonald, Symmetric Functions and Hall Polynomials , Oxford Classic Texts in the Physical Sciences, The Clarendon Press, Oxford University Press, New York, 2015
2015
-
[19]
Motwani and P
R. Motwani and P. Raghavan, Randomized Algorithms, Cambridge University Press, Cambridge, 1995
1995
-
[20]
F. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller and B. V. Saunders, eds, NIST Digital Library of Mathematical Functions, http://dlmf.nist.gov/, Release 1.0.21 of 2018-12-15
2018
-
[21]
Pitman and M
J. Pitman and M. Yor, The law of the maximum of a Bessel bridge, Electron. J. Probab., 4 (1999), no. 15, 35 pp. 20
1999
-
[22]
Reif, Fundamentals of Statistical and Thermal Physics , McGraw-Hill, New York, 1965
F. Reif, Fundamentals of Statistical and Thermal Physics , McGraw-Hill, New York, 1965
1965
-
[23]
Revuz and M
D. Revuz and M. Yor, Continuous Martingales and Brownian Motion , 3rd ed, Springer- Verlag, Berlin, 1999
1999
-
[24]
Schehr, S
G. Schehr, S. N. Majumdar, A. Comtet and J. Randon-Furling, Exact distribution of the maximal height of p vicious walkers, Phys. Rev. Lett., 101 (2008), no. 15, 150601, 4 pp. 21
2008
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.