🔢 mathematics

The pair correlation function of the Sine $_6$ process

This paper derives the first explicit single-variable special function representation of the pair correlation function for the bulk limit of a beta-ensemble beyond the classical values of $\beta=1,2,$ and $4$, specifically expressing it in terms of Bessel functions for the Sine $_6$ process.

Original authors: Shengqi Qiu, Yahui Qu, Lingfan Yuan, Benedek Valkó, Spencer Venancio

Published 2026-07-30

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Shengqi Qiu, Yahui Qu, Lingfan Yuan, Benedek Valkó, Spencer Venancio

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a crowded dance floor where the dancers are invisible, but their movements leave a distinct pattern on the floor. In the world of mathematics and physics, scientists study "ensembles," which are like massive collections of these invisible dancers (often representing energy levels in atoms or eigenvalues in complex matrices). For decades, researchers have been obsessed with understanding how these dancers avoid bumping into each other. They've found that for three specific types of dance floors—called $\beta = 1, 2,$ and $4$—the rules are crystal clear. The dancers follow predictable patterns that can be written down using familiar mathematical tools, like sine waves and simple integrals. It's like knowing the exact steps to a waltz or a tango.

But what happens when the dance floor changes its rules to something in between, or something entirely new? Specifically, what if the "repulsion" between dancers is set to a value called $\beta = 6$ ? This is the "Sine6 process." For a long time, this was a mystery. While we knew the dancers existed and moved in a certain way, no one could write down a single, neat formula to describe exactly how they spaced themselves out. It was like watching a chaotic mosh pit and trying to predict the exact distance between every pair of people without a clear rulebook. Understanding this isn't just about abstract math; it helps physicists and mathematicians understand the fundamental limits of randomness and order in the universe, from the behavior of electrons to the distribution of prime numbers.

This paper, written by Shengqi Qiu, Yahui Qu, Lingfan Yuan, Benedek Valkó, and Spencer Venancio, finally cracks the code for the $\beta = 6$ dance floor. The authors have derived the first explicit, single-variable formula that describes the "pair correlation function" for this specific case. In plain English, they found a mathematical recipe that tells us exactly how likely it is to find two dancers at a specific distance from each other in this chaotic crowd.

Before this discovery, the best we had for $\beta = 6$ was a massive, complicated multi-dimensional integral (a giant, multi-layered math equation involving six different variables at once). It was like trying to describe a symphony by listing every single note played by every instrument simultaneously—it was technically correct but impossible to use or understand intuitively. The authors' breakthrough is that they managed to shrink this giant, six-dimensional monster down into a much more manageable form. They expressed the answer using "Bessel functions," which are special mathematical waves that often appear when dealing with circular or oscillating systems, combined with a single integration (a single sweep of a mathematical curve).

The team didn't just guess this formula; they proved it. They started with a known system of complex differential equations (mathematical rules describing how things change) that governs the $\beta = 6$ process. This system was originally a third-order equation, which is like a car with three gears that are all fighting each other. The authors' clever trick was to show that this complex system could be "reduced" to a simpler, second-order equation. By doing this, they were able to solve it explicitly. Their result is a precise formula involving Bessel functions of the first kind (specifically with indices like $1/6$ and $7/6$ ), sine, cosine, and power functions.

The paper also checks its own work by looking at what happens when the distance between dancers ( $\lambda$ ) gets very large. They show that their new formula matches the expected behavior: the dancers eventually space themselves out evenly, with small ripples of randomness fading away. They even provide a detailed expansion of how these ripples behave, showing terms that oscillate like $\frac{\cos(\lambda)}{\lambda^{2/3}}$ . This confirms that their new, simpler formula isn't just a lucky guess; it captures the true, deep behavior of the system.

In short, this paper takes a problem that was previously stuck in a fog of complex, multi-dimensional math and clears the air with a sharp, single-variable lens. It's the first time anyone has been able to write down the "dance steps" for the Sine6 process in a way that relies on standard special functions, opening the door for easier calculations and deeper understanding of this specific type of mathematical randomness.

Problem Statement
The paper addresses the challenge of deriving an explicit, single-variable special function representation for the pair correlation function of the Sine $_\beta$ process in the bulk scaling limit of Gaussian $\beta$ -ensembles. While explicit formulas exist for the classical values $\beta = 1, 2, 4$ (expressed via determinants, Pfaffians, and sinc functions), no such closed-form expression in terms of standard special functions was known for non-classical values of $\beta$ . Specifically, the authors focus on the case $\beta = 6$ (the Sine $_6$ process), aiming to move beyond integral representations involving multidimensional integrals or generalized hypergeometric functions.

Methodology
The authors build upon the characterization of the pair correlation function $\rho^{(2)}_{2n}(0, \lambda)$ provided in prior work [14], which describes the function as the solution to a complex-valued system of $n$ first-order differential equations. For the case $n=3$ (corresponding to $\beta=6$ ), this system yields a third-order ordinary differential equation (ODE) for the third component of the solution vector, $q_3(\lambda)$ .

The core of the methodology involves a reduction of this third-order ODE to a second-order ODE. The authors demonstrate that the third-order equation can be transformed into a second-order equation with a specific non-homogeneous term involving a polynomial function $h(\lambda)$ . By introducing a change of variables $r(\lambda) = \lambda^2 e^{-3i\lambda/2} q_3(\lambda)$ , the problem is mapped to a second-order linear ODE with variable coefficients. The homogeneous part of this equation is solved using Bessel functions of the first kind, specifically $J_{\pm 1/6}$ and $J_{\pm 7/6}$ . The particular solution is constructed using the method of variation of parameters, leading to an integral representation involving these Bessel functions and the non-homogeneous term.

Key Contributions and Results

Explicit Formula for Sine $_6$ : The primary result (Theorem 1) is an explicit formula for the pair correlation function $\rho^{(2)}_6(0, \lambda)$ . The formula expresses the function in terms of Bessel functions of the first kind ( $J_{\pm 1/6}, J_{\pm 7/6}$ ), elementary functions (exponentials, powers), and a single definite integral. This is the first such representation for a bulk limit of a $\beta$ -ensemble where $\beta$ is not 1, 2, or 4.
ODE Reduction: The paper establishes that the $n$ -dimensional ODE system characterizing $\rho^{(2)}_{2n}$ can be reduced to a single complex-valued ODE of order $n-1$ . This is demonstrated generally in Section 4 using properties of generalized Laguerre polynomials and a specific vector-valued polynomial solution $y_n(\lambda)$ .
Asymptotic Expansion: The authors derive the large $\lambda$ asymptotic behavior of $\rho^{(2)}_6(0, \lambda)$ (Proposition 2). They provide an expansion up to order $O(\lambda^{-11/3})$ , refining previous decay estimates. This expansion includes oscillatory terms involving $\cos(\lambda)$ and $\sin(\lambda)$ with fractional powers of $\lambda$ .
Analytic Properties: The paper notes that while the formula is derived for real $\lambda$ , the function $\rho^{(2)}_6(0, \lambda)$ can be extended to an entire function on the complex plane.

Significance
The paper claims significance primarily for providing the first explicit single-variable special function representation of the bulk pair correlation function for a non-classical $\beta$ value. This breaks the reliance on multidimensional integral representations (such as the Selberg-type integrals found in Forrester's work) or abstract characterizations involving random entire functions. The result bridges the gap between the classical cases ( $\beta=1,2,4$ ) and general $\beta$ , offering a concrete analytical tool for studying the Sine $_6$ process. The authors also note that while similar explicit results exist for the soft edge scaling limit (e.g., Tracy-Widom distributions for $\beta=6$ via Painlevé equations), this work addresses the bulk limit, which had previously lacked such a representation.

More like this