---
title: 'Bessel Point Process: Hard-Edge Scaling in Random Matrices'
url: https://www.emergentmind.com/topics/bessel-point-process
type: topic
---

# Bessel Point Process: Hard-Edge Scaling in Random Matrices

The Bessel point process is the determinantal point process on $(0,\infty)$ governed by the Bessel kernel with parameter $\alpha>-1$; the same parameter is often denoted by $\nu$ or $a$ in the literature. It is the hard-edge scaling limit for eigenvalues near the origin in Laguerre or Wishart ensembles, and it is distinguished by projection-kernel determinantal correlations, rigidity, explicit conditional measures, and a substantial integrable-systems structure involving Fredholm determinants, Riemann–Hilbert problems, and Painlevé-type equations [1506.07581, 2109.09070, 2401.11213].

## 1. Definition and determinantal structure

Let $\mathbb R_+=(0,\infty)$ with Lebesgue measure. A standard form of the Bessel kernel is
\[
K^{\mathrm{Bessel}_\alpha}(x,y)
=
\frac{J_\alpha(\sqrt x)\,\sqrt y\,J_\alpha'(\sqrt y)-\sqrt x\,J_\alpha'(\sqrt x)\,J_\alpha(\sqrt y)}{2(x-y)},
\qquad \alpha>-1,
\]
where $J_\alpha$ is the Bessel function of the first kind of order $\alpha$. An equivalent Tracy–Widom form is
\[
K^{\mathrm{Bessel}_\alpha}(x,y)
=
\frac{J_\alpha(\sqrt x)\,\sqrt y\,J_{\alpha+1}(\sqrt y)-J_\alpha(\sqrt y)\,\sqrt x\,J_{\alpha+1}(\sqrt x)}{2(x-y)}.
\]
By the Macchi–Soshnikov theorem, the associated integral operator on $L^2(\mathbb R_+,dx)$ is an orthogonal projection of locally trace-class, and therefore defines a unique determinantal point process on $\operatorname{Conf}(\mathbb R_+)$ [1506.07581].

Its $n$-point correlation functions are
\[
\rho_n(x_1,\dots,x_n)=\det\!\bigl[K(x_i,x_j)\bigr]_{i,j=1}^n,
\]
in the sense that for test functions on $(\mathbb R_+)^n$ the joint intensities are given by determinants of the kernel. Equivalent integral representations also appear. One common hard-edge form is
\[
J_\nu(x,y)=\int_0^1 J_\nu(\sqrt x\,t)\,J_\nu(\sqrt y\,t)\,dt,
\]
while a scaled representation used in line-ensemble and field formulations is
\[
K_\alpha(x,y)=\int_0^{1/4} J_\alpha(\sqrt{zx})\,J_\alpha(\sqrt{zy})\,dz
\]
up to an innocuous gauge, or
\[
K^{\mathrm{Bes}_\alpha}(x,y)=\int_0^1 J_\alpha(2\sqrt{x\,u})\,J_\alpha(2\sqrt{y\,u})\,du
\]
in a different normalization [1904.04349, 2109.09070, 2109.09292].

The determinantal formulation immediately yields Fredholm-determinant expressions for gap events. For example, the probability of no points in a bounded set such as $[0,R]$ is expressed as the Fredholm determinant of the cut-down Bessel operator on $L^2([0,R])$ [2403.16219].

## 2. Hard-edge origin and dynamical realizations

The principal random-matrix origin is the Laguerre Unitary Ensemble. If $A$ is an $N\times(N+\alpha)$ complex Brownian matrix at time $1/2$, then $M=A A^\ast$ has the Wishart law, and its ordered eigenvalues $0\le \lambda_1<\cdots<\lambda_N$ have joint density
\[
C(N,\alpha)\,\prod_{j=1}^N \lambda_j^\alpha e^{-\lambda_j}\,\Delta(\lambda)^2\,1_{0\le \lambda_1<\cdots<\lambda_N}\,d\lambda.
\]
Under hard-edge scaling near the smallest eigenvalues, one lets $N\to\infty$ and sets $x=4N\lambda$, with time scaled as $t\mapsto 1+t/(4N)$. The resulting point process converges to the Bessel point process of index $\alpha$ [2109.09070].

A dynamical realization is given by the Dyson Bessel process: one starts with $N$ independent squared Bessel processes of index $\alpha$, begun at $0$ and conditioned never to collide. Under the same hard-edge scaling near time $1$, the top curves converge to a countable Bessel line ensemble whose one-time marginal is the Bessel point process. Its finite-dimensional distributions are governed by the extended Bessel kernel, and the limiting ensemble is stationary under horizontal translations because the extended kernel depends on time only through $s-t$ in the exponential factor [2109.09070].

A broader two-parameter object is the Bessel field $\mathcal B(\alpha,t)$. For fixed $(\alpha,t)$ it is the classical Bessel point process. Along any time-like path or space-like path in the $(\alpha,t)$ plane, the joint configuration is determinantal with an explicit extended kernel built from Bessel functions and exponential propagation factors. For fixed $t$, the field in the $\alpha$-direction is described as an exponential Gibbsian line ensemble; for fixed $\alpha$, the $t$-evolution is a squared Bessel Gibbsian line ensemble [2109.09292].

## 3. Rigidity and conditional measures

A point process on $\mathbb R_+$ is called rigid if for every bounded Borel set $B\subset \mathbb R_+$ the counting variable
\[
N_B(\xi)=|\xi\cap B|
\]
is almost surely determined by the restriction of the configuration to $\mathbb R_+\setminus B$. For every $\alpha>-1$, the Bessel point process is rigid. In particular, if $0<a<b<\infty$, then the number of points in $[a,b]$ is almost surely a measurable function of the configuration outside $[a,b]$ [1506.07581].

The proof uses additive statistics approximating the indicator of $[a,b]$. For an increasing sequence $T_n\to\infty$, define cutoff functions
\[
\psi_n^{(R)}(x)=
\begin{cases}
1, & 0<x\le R,\\[3pt]
1-\dfrac{\log(x/R)}{\log(T_n/R)}, & R<x\le T_n,\\[6pt]
0, & x>T_n,
\end{cases}
\]
and then set
\[
\phi_n(x)=\psi_n^{(b)}(x)-\psi_n^{(a)}(x),\qquad
A_n(\xi)=\sum_{x\in\xi}\phi_n(x).
\]
By construction, $\phi_n\equiv 1$ on $[a,b]$, $0\le \phi_n\le 1$, and $\operatorname{supp}\phi_n\subset[a,T_n]$. For bounded $\phi$ one has
\[
\operatorname{Var}\Bigl(\sum_{x\in\xi}\phi(x)\Bigr)
=
\iint_{\mathbb R_+\times\mathbb R_+} |\phi(x)-\phi(y)|^2\,|K^{\mathrm{Bessel}_\alpha}(x,y)|^2\,dx\,dy.
\]
Using the known small-argument and large-argument asymptotics of $J_\alpha$, together with the facts that for $x,y\gg1$ the Bessel kernel decays like $O((xy)^{-1/4}|x-y|^{-1})$ and that for $x$ or $y$ small the kernel is integrable, one obtains $\operatorname{Var}(A_n)\to0$, and the Ghosh–Peres criterion then yields rigidity [1506.07581].

Rigidity has a stronger conditional consequence. For a realization $X=(p_n)$ of the Bessel process and $R>0$, write $N(R)=\#\{p_n\le R\}$. The conditional distribution of the points inside $[0,R]$, given the configuration outside $[0,R]$, is almost surely an orthogonal polynomial ensemble on $[0,R]$ with exactly $N(R)$ points and weight
\[
\bar w_{X,\nu,R}(t)=t^\nu\prod_{p_n>R}\Bigl(1-\frac{t}{p_n}\Bigr)^2,\qquad t\in[0,R].
\]
Its joint density for ordered points $0\le t_1<\cdots<t_{N(R)}\le R$ is
\[
\frac1{Z_{X,\nu,R}}
\prod_{1\le i<j\le N(R)}(t_i-t_j)^2
\prod_{i=1}^{N(R)}\bar w_{X,\nu,R}(t_i),
\]
so the conditional law is again determinantal, now with a Christoffel–Darboux kernel [1904.04349].

Molag and Stevens proved a universality statement for these conditional measures. If $X=(p_n)$ is any strictly increasing sequence in $(0,\infty)$ satisfying
\[
\lim_{n\to\infty}\frac{p_n}{n^2}=\pi^2,
\]
then for every compact $K\subset(0,\infty)^2$,
\[
\lim_{R\to\infty}
K_{N(R)}\!\bigl(x,y;\bar w_{X,\nu,R}\bigr)
=
J_\nu(x,y)
\]
uniformly on $K$. They also verify, via Soshnikov’s variance estimates, that a random Bessel configuration almost surely satisfies
\[
p_n=\pi^2 n^2+o(n^2),
\]
so almost every conditional orthogonal polynomial ensemble recovers the Bessel point process as $R\to\infty$ [1904.04349].

## 4. Fluctuations, counting statistics, and large gaps

For the counting function
\[
N(0,R)=\#\{\text{points in }(0,R)\},
\]
exponential-moment asymptotics are explicit. The moment-generating function
\[
M_\nu(\theta;[0,R])=\mathbb E[e^{\theta\,N(0,R)}]
\]
satisfies
\[
\log M_\nu(\theta;[0,R])
=
\theta\,\mu_\nu(R)
+\frac{\theta^2}{4}\,\log R
+\log\!\Bigl[G\!\bigl(1+\tfrac{\theta}{2\pi i}\bigr)\,
G\!\bigl(1-\tfrac{\theta}{2\pi i}\bigr)\Bigr]
+O(R^{-1/2}),
\]
where
\[
\mu_\nu(R)=\frac{\sqrt R}{\pi}-\nu
\]
and $G$ is Barnes’ $G$-function. Consequently,
\[
\mathbb E[N(0,R)]=\frac{\sqrt R}{\pi}-\nu+O(R^{-1/2}),
\]
and
\[
\operatorname{Var}(N(0,R))
=
\frac{\log R}{4\pi^2}
+\frac{1+\gamma_E+\log 4}{2\pi^2}
+O(R^{-1/2}).
\]
After normalization by the variance, $N(0,R)$ converges in distribution to $N(0,1)$ [1812.02188].

More general linear statistics also admit quantitative Gaussian asymptotics. For compactly supported $b$, define
\[
S_R=\sum_{x\in X} b(x/R).
\]
Under the centering condition $b(0)=0$, one has
\[
\operatorname{Var}(S_R)\sim 2\,c_3(b)\,\log R,
\qquad
\frac{S_R-\mathbb E[S_R]}{\sqrt{\operatorname{Var}(S_R)}}\Longrightarrow N(0,1),
\]
and, more precisely, if $F_{R,b}$ denotes the distribution function of
\[
\frac{S_R-\mathbb E[S_R]}{\sqrt{2\,c_3(b)\,\log R}},
\]
then
\[
\sup_{x\in\mathbb R}|F_{R,b}(x)-\Phi(x)|
\le
\frac{C}{\sqrt{\log R}}
\]
for all $R\ge1$, where $\Phi$ is the standard normal distribution function [2403.16219].

Large-gap asymptotics are likewise explicit. For disjoint consecutive intervals $I_j=(x_{j-1},x_j]$, piecewise thinning with retention probabilities $p_j$ yields a thinned process whose joint gap probability is again a Fredholm determinant. Writing $p_j=e^{u_j-u_{j+1}}$ with $u_{m+1}=0$, the exponential moment
\[
E_\nu(x,u)
=
\mathbb E\!\Bigl[\exp\Bigl(\sum_{j=1}^m u_j\,N(r x_{j-1},r x_j)\Bigr)\Bigr]
\]
has an asymptotic expansion involving
\[
\sigma(x)=\frac{\sqrt x}{2\pi^2},
\qquad
\Sigma(x,y)=\frac1{2\pi^2}\ln\!\frac{\sqrt x+\sqrt y}{|\sqrt x-\sqrt y|},
\]
together with Barnes $G$-factors and a quadratic form in the $u_j$ [1812.02188].

For unions of widely separated intervals, large-$r$ gap probabilities acquire a theta-function structure. If
\[
I=(0,x_1)\cup(x_2,x_3)\cup\cdots\cup(x_{2g},x_{2g+1}),
\qquad 0<x_1<\cdots<x_{2g+1},
\]
then
\[
F(x_1,\dots,x_{2g+1})
=
\mathbb P\{\text{no points in }I\}
=
\det(I-K_\alpha|_{L^2(I)}).
\]
As $r\to\infty$, $\log F(r x_1,\dots,r x_{2g+1})$ admits an expansion of the form
\[
c\,r-d_1\sqrt r+\log\theta(\nu(r))
-\frac1{32}\sum_{j=1}^{2g+1}\int_M^r B(-x_j,\nu(t))\,dt
+C+O(r^{-1/2}),
\]
where the coefficients are described through a genus-$g$ Riemann surface, normalized holomorphic one-forms, and a linear flow on a $g$-dimensional torus. In the ergodic case, Birkhoff’s ergodic theorem identifies the leading logarithmic term explicitly [2101.09216].

## 5. Fredholm determinants, Riemann–Hilbert problems, and Painlevé structures

A central integrable observable is the joint generating function of occupancy numbers on disjoint intervals. For
\[
0=x_0<x_1<\cdots<x_k<\infty,
\qquad
N_j=\#\{\text{particles in }(x_{j-1},x_j)\},
\]
set
\[
F(\mathbf x,\mathbf s)=\mathbb E\Bigl[\prod_{j=1}^k s_j^{N_j}\Bigr].
\]
This function is a Fredholm determinant,
\[
F(\mathbf x,\mathbf s)=\det[I-\mathbf K_{\vec s}]_{L^2(0,x_k)},
\]
with deformed kernel
\[
\mathbf K_{\vec s}(u,v)
=
\sum_{j=1}^k (1-s_j)\,\chi_{(x_{j-1},x_j)}(u)\,K_{\rm Be}(u,v).
\]
Charlier and Doeraene derive a $2\times2$ Riemann–Hilbert problem with simple poles at $0,-x_1,\dots,-x_k$; from its Lax pair they obtain a system of $k$ coupled Painlevé V equations for functions $q_j$. Their main representation writes
\[
\ln F(x\vec r,\vec s)
=
\sum_{j=1}^k (s_{j+1}-s_j)\int_0^x \ln(\xi)\,q_j(\xi)^2\,d\xi,
\]
with hard-edge boundary conditions
\[
q_j(x)=\sqrt{s_{j+1}-s_j}\,J_a(\sqrt{x\,r_j})\,(1+O(x)),
\qquad x\to0.
\]
For $k=1$, this reduces to the Tracy–Widom description of the one-interval Bessel gap probability in terms of a single Painlevé V equation [1709.07365].

A more general framework treats multiplicative statistics depending on two parameters. For a bounded measurable weight $\sigma$, define
\[
Q[\sigma](x,t)=\det[I-M_{\sigma,x,t}K_\nu M_{\sigma,x,t}]_{L^2(0,\infty)},
\qquad
(M_{\sigma,x,t}\psi)(\lambda)=\sqrt{\sigma(x^{-2}\lambda+t)}\,\psi(\lambda).
\]
Introduce
\[
v(x,t)=-\partial_x\log Q[\sigma](x,t)-\frac{4\nu^2-1}{8x}.
\]
Then $v$ satisfies the nonlinear integrable PDE
\[
(2\,v_{xt}\,v_x)+v_t\,v_{xxt}-v_x\,v_{xtt}=\frac{\nu^2}{4}.
\]
This equation is equivalent to the compatibility of the linear system
\[
L f=\lambda f,
\qquad
(\lambda-t)f=Bf,
\]
with
\[
L=-\partial_x^2+2\,v_{xt}(x,t),
\qquad
B=-v_x\partial_x+\frac{v_{xt}}2,
\]
or, in Lax form,
\[
\partial_t L=[B,L].
\]
Thus the $t$-deformation is isospectral for a Sturm–Liouville problem. When $\sigma$ is specialized to a step function, the PDE reduces to the coupled Painlevé V system above; in the single-interval case it reduces to the Tracy–Widom Painlevé V equation [2401.11213].

This integrable description is also compatible with the Its–Izergin–Korepin–Slavnov theory of integrable operators and a Bessel-model Riemann–Hilbert problem. The resulting identities connect Fredholm determinants, boundary-value problems for the associated Sturm–Liouville equation, and nonlocal or coupled Painlevé-type reductions [2401.11213].

## 6. Related processes, limits, and variants

The infinite Bessel process arises in the study of infinite Pickrell measures. For $s>-1$, one defines a modified kernel
\[
J^{(s)}(x,y)=\int_0^1 J_s(\sqrt{t\,x})\,J_s(\sqrt{t\,y})\,dt,
\]
equivalently written in a Tracy–Widom-type Christoffel–Darboux form. The corresponding determinantal process is obtained as the hard-edge scaling limit of radial parts of finite-dimensional Pickrell measures after the rescaling $\lambda\mapsto n^2\lambda$. Different parameters give mutually singular laws, and the ergodic decomposition analysis shows that the additional “Gaussian parameter” vanishes almost surely [1610.04646].

A distinct but related object is the finite-temperature deformation of the discrete Bessel point process on
\[
\mathbb Z'=\{\dots,-3/2,-1/2,1/2,3/2,\dots\}.
\]
For a decay weight $\sigma:\mathbb Z'\to[0,1]$ and $L>0$, the kernel
\[
K_\sigma(a,b)=\sum_{l\in\mathbb Z'} \sigma(l)\,J_{a+l}(2L)\,J_{b+l}(2L)
\]
defines a determinantal point process on $\mathbb Z'$. Its largest-particle distribution
\[
Q_\sigma(L,s)=\Pr[a_{\max}\le s]=\det(1-P_sK_\sigma P_s)
\]
satisfies a reduction of the $2$D Toda hierarchy,
\[
\frac{\partial^2}{\partial L^2}\ln q_s(L)+L\frac{\partial}{\partial L}\ln q_s(L)+4
=
4\,\frac{q_{s+1}(L)\,q_{s-1}(L)}{q_s(L)^2},
\]
together with a discrete integro-differential Painlevé II system. Under a formal continuum scaling, this deformation converges to the finite-temperature Airy framework and its KdV/integro-Painlevé II structure [2207.01421].

Another direction replaces the determinantal kernel by a random differential operator. For the stochastic Bessel operator
\[
\mathfrak G^{\beta,a}
=
-\frac1{m(x)}\,d\Bigl[\frac1{s(x)}\,d\,\cdot\,\Bigr],
\]
with random weights
\[
m(x)=\exp\!\bigl(-(a+1)x-(2/\sqrt\beta)B(x)\bigr),
\qquad
s(x)=\exp\!\bigl(ax+(2/\sqrt\beta)B(x)\bigr),
\]
the low-lying eigenvalues define a point process at the hard edge. In the high-temperature limit $\beta\to0$, the rescaled eigenvalues
\[
\mu^\beta(k)=\beta\ln(1/\Lambda^{\beta,a}(k))
\]
converge to a limiting point process $M_0$ characterized through coupled stochastic differential equations and an alternating reflected Brownian-motion construction. No closed-form joint density, Fredholm determinant, or Pfaffian kernel is known for $M_0$, and the limiting process is not Poisson [2411.13972].

These variants clarify a structural boundary. The classical $\beta=2$ hard-edge Bessel point process is determinantal and admits exact Fredholm and Painlevé descriptions; discrete, finite-temperature, infinite-dimensional, and high-temperature analogues preserve some of its hard-edge or Bessel-function features, but they generally modify the kernel structure, the integrable hierarchy, or the determinantal character itself [2207.01421, 2411.13972].

Source: https://www.emergentmind.com/topics/bessel-point-process