The pair correlation function of the Sine$_6$ process

Benedek Valk\'o; Lingfan Yuan; Shengqi Qiu; Spencer Venancio; Yahui Qu

REVIEW 2 major objections 5 minor 18 references

The pair correlation function of the Sine$_6$ process

T0 review · 2 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash

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

Pith's one-line read This paper derives an explicit single-variable formula for the pair correlation function of the Sine_6 process, in terms of Bessel functions and one integral.

desk verdict Explicit Bessel-function formula for the Sine_6 pair correlation; a real new result with a legitimate but addressable reliance on the authors' own ODE characterization. read the letter →

arxiv 2607.26223 v1 pith:KSCX5HUT submitted 2026-07-28 math.PR math-phmath.MP

Shengqi Qiu , Yahui Qu , Lingfan Yuan , Benedek Valkó , Spencer Venancio This is my paper · ORCID

classification math.PRmath-phmath.MP MSC 60B2060G5533C10

keywords Sine_6processGaussianbeta-ensemblepaircorrelationfunctionBesselfunctionspointrandommatrixtheoryODEreduction

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 aims to show that the pair correlation function of the Sine_6 point process—the bulk scaling limit of the Gaussian beta-ensemble at β=6—can be written as an explicit expression built from sines, powers, and Bessel functions of order ±1/6 and ±7/6, together with a single integral. This would be the first single-variable special-function representation of a bulk pair correlation for a beta-ensemble outside the classical values β=1,2,4, where determinantal and Pfaffian formulas have long been known. The derivation starts from an ODE characterization and reduces a third-order differential equation to a second-order one, whose solution is found by variation of parameters. If correct, the formula turns a six-dimensional integral representation into a one-dimensional expression, making asymptotics and numerical evaluation much more tractable.

What carries the argument

The central mechanism is the reduction of the n=3 case of the ODE system (15) to a single third-order equation for q_3, and then to the second-order linear equation (28) via an explicit integrating-factor relation (λw' + (4−3iλ)w = 80). The reduced equation becomes a Bessel-type equation λr'' + (7/3)r' + (λ/4 + i/2)r = forcing after the change of variable r=λ² e^{-3iλ/2} q_3; its homogeneous solutions are the combinations r±(λ)=λ^{-1/6}(J_{±1/6}(λ/2) ± i J_{±7/6}(λ/2)), and variation of parameters gives the integral formula. The same idea is generalized in Lemma 4: a vector of Laguerre polynomials produces a solution z_n = e^{-inλ}λ^{n+1} y_n of the adjoint equation, which reduces the order-

What would settle it

Evaluate both sides of Theorem 1 at a specific nonzero λ, say λ=1, by high-precision numerical integration, and compare with a direct high-precision evaluation of the six-dimensional integral (10); any discrepancy beyond roundoff would disprove the formula. Alternatively, check that the q_3 constructed from the formula satisfies the third-order ODE (27) and the initial conditions q_3(0)=1, q_3'(0)=i, q_3''(0)=−9/8.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is Theorem 1: for λ>0, 4π² ρ_6^(2)(0,λ) = Re[1 − (4/27) λ^{-4} h(λ) + ℓ_+(λ)∫_0^λ b(s) r_-(s) ds − ℓ_-(λ)∫_0^λ b(s) r_+(s) ds], where h, b, r±, and ℓ± are explicit combinations of exponentials, powers, and Bessel functions of the first kind of orders ±1/6 and ±7/6. The authors prove this by taking a known ODE characterization of the correlation function for β=6, reducing the resulting third-order ODE for a coordinate q_3 to a second-order ODE, solving that equation in closed form using the two independent Bessel-type solutions r±, and then reconstructing the correlation function. They also show the same reduction idea works for every even β=2n,

Load-bearing premise

The load-bearing premise is that the ODE characterization quoted from the companion paper (Theorem 3, equations (15)–(17)) truly identifies ρ_6^(2); if that characterization is wrong or does not apply, the explicit formula does not describe the Sine_6 pair correlation.

Editorial extensions

If this is right

For β=6, the pair correlation is now computable by a single one-dimensional integral involving Bessel functions, rather than the six-dimensional integral representation in equation (10).
The formula yields a sharp large-λ expansion (Proposition 2): 4π²ρ_6^(2) = 1 − 2/(3λ²) + Γ(1/3)²[(2/9)cosλ/λ^{2/3} − (16/81)sinλ/λ^{5/3} − (64/729)cosλ/λ^{8/3} + (8/81)cos(2λ)/λ^{8/3}] + O(λ^{−11/3}).
For every even β=2n, the pair correlation can be obtained from a complex-valued ODE of order n−1, reducing the dimension of the previous description.
The small-λ behavior ρ_6^(2)(0,λ) = const·λ^6 + O(λ^8) follows from the integral representation and is consistent with the new formula.

Reading between the lines

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

If the same reduction can be pushed beyond n=3, explicit Bessel-type formulas may exist for β=10, 14, and other even integers; the general Lemma 4 gives a concrete route to test this by deriving the next-order ODE and searching for closed-form solutions.
The explicit formula makes it possible to check universality predictions for non-classical β numerically at finite λ, for example by comparing one-dimensional quadrature against the six-dimensional integral (10) at a few values of λ.
The Bessel-function structure of the β=6 bulk correlation is reminiscent of Painlevé-type representations seen for the soft-edge β=6 distribution, suggesting possible deeper connections between even-integer β correlations and hierarchies of special functions.
Because the formula is entire in λ, it could serve as a generating object for further identities involving Sine_6 correlation functions or related gap probabilities, though the paper does not pursue this.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

Main take: the paper proves a genuine new formula, Theorem 1, for the pair correlation of Sine_6 as a one-dimensional integral of Bessel functions. That is the first explicit single-variable special-function expression for a non-classical beta bulk correlation, and it goes beyond Forrester's multidimensional integral (9)-(10). The proof is explicit: from the 3x3 system (22) they reduce the resulting third-order ODE to a second-order ODE, identify homogeneous solutions as Bessel functions, and use variation of parameters. I checked the Wronskian sign and the variation-of-parameters step; the algebra is consistent. The paper also includes a general-n reduction to an (n-1)-order ODE, but only as a sketch, and an asymptotic expansion (Prop. 2) that is outlined rather than proved.

Where the soft spots are: the starting point, Theorem 3 of [14], is a preprint by two of the current co-authors. The paper neither proves nor numerically checks that theorem. So the derivation shows formula (8) follows from (22), but the identification of (22) with Sine_6 rests entirely on [14]. That is a genuine weakness, but it is not a circular argument—the target formula is not assumed. It is a reproducibility burden: a referee should verify Theorem 3 or ask the authors to include a check. The general-n reduction and the asymptotics section are also underdeveloped; they are auxiliary, not load-bearing for Theorem 1. Minor: the definition of fractional powers with principal branch is fine, but the integrals at lambda=0 need the limiting interpretation, which they state.

The citation pattern is appropriate: Forrester's integral, Rumanov's Painleve result, and the recent [3] and [5] are all credited. Self-citation of [14] is justified because the proof genuinely depends on it. No fitting, no free parameters, no invented objects.

Overall: a solid subfield contribution. The main theorem looks right, the proof is checkable, and the unverified input is pre-existing work by the same group rather than a hidden assumption. It deserves a serious referee; the referee's main job is to verify Theorem 3 of [14] and the reduction algebra. I would bring it to my reading group and cite it if I were working on Sine_beta correlations.

Referee Report

2 major / 5 minor

Summary. The paper derives an explicit formula for the pair correlation function ρ_6^{(2)}(0,λ) of the Sine_6 process in terms of Bessel functions of the first kind (Theorem 1, Eq. (8)). The proof specializes the ODE characterization of ρ_{2n}^{(2)} given in Theorem 3 (quoted from the preprint [14]) to n=3, reduces the resulting third-order ODE (27) to a second-order nonhomogeneous equation (28) via an explicit integrating-factor-like step, and solves it by variation of parameters using the Bessel-type homogeneous solutions (6). The paper also sketches a general reduction for β=2n (Section 4) and states large-λ asymptotics (Proposition 2), with an outline of how they follow from both the known multiple integral (10) and from Theorem 1.

Significance. If Theorem 3 of [14] is valid, formula (8) provides the first single-variable special-function representation of a Sine_β pair correlation function outside the classical cases β=1,2,4. The internal derivation is explicit and checkable, and the final formula is concrete enough for numerical evaluation and asymptotic analysis. The result is a worthwhile contribution to random matrix theory, and the connection to Bessel functions is elegant. However, the paper's sole connection to the Sine_6 process is an unproved same-author preprint theorem, which is a significant caveat on the otherwise interesting derivation.

major comments (2)

[Section 2, Theorem 3] The proof of Theorem 1 depends entirely on Theorem 3 of [14] (Eqs. (15)–(17)), which supplies both the ODE system (22) and the identification ρ_6^{(2)}=(4π²)^{-1}(1+2v^T Re q). This theorem is not proved or independently verified here, and [14] is a preprint by two of the co-authors. If there is a defect in the matrices (13), the vector (14), or the sign in (17), formula (8) may be a correct consequence of (22) but would not represent the Sine_6 pair correlation. I request that the authors make the n=3 input self-contained (e.g., by proving the needed part of Theorem 3 in an appendix) or provide a numerical comparison of (8) with the known integral representation (10) over a range of λ.
[Section 3, Eq. (28)] The reduction from (27) to (28) is stated as ‘direct differentiation’ without the computation, and the passage from (34)–(36) to (37) is similarly condensed. These algebraic manipulations are load-bearing: an error there would change the final formula. Since the proof is not otherwise computational, I ask that the intermediate expressions be provided, either in the text or in a supplementary appendix, so that the verification does not require the referee to redo the entire reduction.

minor comments (5)

[Eq. (32)] The leading asymptotic for r_-(λ) appears to have the wrong sign: from definition (6), J_{-7/6}(λ/2)~(λ/4)^{-7/6}/Γ(-1/6), so r_-(λ) ~ i 2^{7/3} Γ(-1/6)^{-1} λ^{-4/3}, with Γ(-1/6)<0. The displayed term in (32) has an additional minus sign. Please check and correct.
[Eq. (9)] The Selberg integral notation S_{2n}(-1+1/n, -1+1/n, 1/n) seems inconsistent with the standard convention S_n(α,β,γ)=∫∏ x_i^{α-1}(1-x_i)^{β-1}|Δ|^{2γ}. The integrand in (9) has exponents -1+1/n for the monomials, which would correspond to α=β=1/n. If this is intentional, please add a definition; otherwise correct the parameters.
[Section 4, after Eq. (39)] The general reduction for β=2n is only sketched: the reader is told that substituting the expressions for q_1,...,q_{n-1} into (39) leads to an ODE of order n-1 for q_n, with nonzero leading coefficient, but neither the ODE nor the nonvanishing is shown. Since this is a secondary claim, a more explicit statement (or a reference to a supplement) would be helpful.
[Section 3, around Eq. (34)] The initial conditions q_3(0)=1, q'_3(0)=i, q''_3(0)=-9/8 are quoted from Theorem 3's series recursion (16) without derivation. A one-line computation would make this easier to follow.
[Section 5, Proposition 2] The direct derivation of (11) from Theorem 1 is described only by an outline; if Proposition 2 is intended to be fully proved from the new formula, more details of the endpoint analysis are needed. The derivation from (10) via [8] is acceptable, but the claim that Theorem 1 also yields (11) should be either substantiated or labeled as a sketch.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the derivation is a direct algebraic reduction from an external ODE characterization, and the target formula is not assumed anywhere.

full rationale

The paper's main theorem (8) is derived by starting from Theorem 3 of [14], a self-cited preprint that characterizes the pair correlation function as the solution of an ODE system. The present paper then solves that ODE for n=3: it rewrites the 3x3 system (22) as a third-order ODE (27), reduces it to a second-order equation (28), solves the latter via Bessel functions and variation of parameters, and finally assembles the explicit expression (8). At no point is the target formula (8) assumed or fitted; it is produced by explicit computation. The self-citation is load-bearing in the sense that the connection to Sine_6 passes through [14], but that cited result is a distinct characterization (an ODE) rather than a disguised version of the conclusion. The paper also uses Forrester's independent integral representation (10) for asymptotics, providing an external checkpoint. No parameter is fitted and no equation is equivalent by construction to the final formula. The only potential concern is that [14] is a preprint by two of the current co-authors and its proof is not reproduced here, but that is a provenance/verifiability issue, not circularity: the cited theorem does not include (8), and the derivation from it is mathematically substantive.

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

There are no fitted free parameters: all constants in (5)-(8) are fixed by β=6, n=3. The axioms are standard Bessel-function theory, the definition of Sine_beta as a scaling limit, and the cited ODE characterization from [14]. No new entities are postulated; the auxiliary functions are explicit combinations of known functions.

assumptions (5)

domain assumption Sine_beta is the bulk scaling limit of the Gaussian beta-ensemble and equals the circular beta-ensemble limit.
Introduction, paragraphs 2-3; cites [16,17,12,13]. Needed to define the object whose pair correlation is computed.
standard math Forrester's integral representation (9) for rho_2n is valid.
Equation (9), citing [6,7,8]; used for asymptotics (Proposition 2) and as external context.
domain assumption Theorem 3 of [14]: rho_2n = (1/4π²)(1+2 v_n^T Re q_n), with q_n the entire solution of ODE (15).
Section 2, Theorem 3; this is the starting point of the proof of Theorem 1 and is not proved in the present paper. Self-citation to two current co-authors.
standard math Standard Bessel function identities (18)-(21), including the Wronskian (21).
Section 2, equations (18)-(21); used for variation of parameters and homogeneous solutions.
standard math Variation of parameters is valid on C\(−∞,0] and the integrals in (33) converge.
Section 3, after equation (32); relies on standard ODE theory and the asymptotics (32).

how reviews work

0 comments

Cite this review

Pith. "Pith review of The pair correlation function of the Sine$_6$ process." pith.science (2026). https://pith.science/paper/KSCX5HUT

@misc{pith2026260726223,
  author       = {Pith},
  title        = {Pith review of: The pair correlation function of the Sine$_6$ process},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KSCX5HUT}},
  note         = {Machine review of arXiv:2607.26223}
}

abstract

We derive an explicit formula for the pair correlation function of the Sine$_6$ process in terms of Bessel functions of the first kind. This provides the first 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$.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

18 extracted references

[14]

On the pair correlation function of the Sine β process

Yahui Qu and Benedek Valk ´o. On the pair correlation function of the Sine β process. https://arxiv.org/abs/2509.15446, 2025

arXiv 2025
[3]

Moments of CβE field partition function, Sineβ correlations and stochastic zeta.https://arxiv.org/abs/2602.08739, 2026

Theodoros Assiotis and Joseph Najnudel. Moments of CβE field partition function, Sineβ correlations and stochastic zeta.https://arxiv.org/abs/2602.08739, 2026

arXiv 2026
[5]

Long-Range Correlation of the Sineβ point Process

Laure Dumaz and Martin Malvy. Long-Range Correlation of the Sineβ point Process. https://arxiv.org/abs/2603.15289, 2026

arXiv 2026
[1]

Stegun.Handbook of mathematical functions with formulas, graphs, and mathematical tables.Wiley, 1984

Milton Abramowitz and Irene A. Stegun.Handbook of mathematical functions with formulas, graphs, and mathematical tables.Wiley, 1984

1984
[2]

Anderson, Alice Guionnet, and Ofer Zeitouni.Introduction to random matri- ces

Greg W. Anderson, Alice Guionnet, and Ofer Zeitouni.Introduction to random matri- ces. Cambridge University Press, 2009. 10

2009
[4]

Universality of generalβ- ensembles.Duke Math

Paul Bourgade, L ´aszl´o Erd ˝os, and Horng-Tzer Yau. Universality of generalβ- ensembles.Duke Math. J., 163(6):1127–1190, 2014

2014
[6]

Forrester

Peter J. Forrester. Selberg correlation integrals and the1/r 2 quantum many-body system.Nuclear Physics B, 388(3):671–699, 1992

1992
[7]

Forrester

Peter J. Forrester. Addendum to ‘Selberg correlation integrals and the1/r 2 quantum many body system’.Nuclear Physics B, 416(1):377–385, 1994

1994

Show all 18 references

[8]

Forrester.Log-gases and random matrices, volume 34 ofLondon Mathematical Society Monographs Series

Peter J. Forrester.Log-gases and random matrices, volume 34 ofLondon Mathematical Society Monographs Series. Princeton University Press, Princeton, NJ, 2010

2010
[9]

Forrester

Peter J. Forrester. Differential identities for the structure function of some random matrix ensembles.Journal of Statistical Physics, 183(2), May 2021

2021
[10]

On the Tracy-Widomβ distribution forβ= 6.Symmetry, Integrability and Geometry: Meth- ods and Applications, November 2016

Tamara Grava, Alexander Its, Andrei Kapaev, and Francesco Mezzadri. On the Tracy-Widomβ distribution forβ= 6.Symmetry, Integrability and Geometry: Meth- ods and Applications, November 2016

2016
[11]

On theβ= 6Tracy-Widom distribution and the second Calogero-Painlev ´e system.https://arxiv.org/abs/2010.06733, 2020

Alexander Its and Andrei Prokhorov. On theβ= 6Tracy-Widom distribution and the second Calogero-Painlev ´e system.https://arxiv.org/abs/2010.06733, 2020

2010
[12]

Eigenvalue statistics for CMV matrices: from Pois- son to clock via random matrix ensembles.Duke Math

Rowan Killip and Mihai Stoiciu. Eigenvalue statistics for CMV matrices: from Pois- son to clock via random matrix ensembles.Duke Math. J., 146(3):361–399, 2009

2009
[13]

Level statistics for one-dimensional Schr ¨odinger operators and Gaussian beta ensemble.Journal of Statistical Physics, 156(1):66–93, April 2014

Fumihiko Nakano. Level statistics for one-dimensional Schr ¨odinger operators and Gaussian beta ensemble.Journal of Statistical Physics, 156(1):66–93, April 2014

2014
[15]

Painlev ´e representation of Tracy–Widom β distribution forβ= 6

Igor Rumanov. Painlev ´e representation of Tracy–Widom β distribution forβ= 6. Communications in Mathematical Physics, 342(3):843–868, November 2015. 11

2015
[16]

Continuum limits of random matrices and the Brownian carousel.Inventiones Math., 177:463–508, 2009

Benedek Valk ´o and B ´alint Vir ´ag. Continuum limits of random matrices and the Brownian carousel.Inventiones Math., 177:463–508, 2009

2009
[17]

The Sine β operator.Inventiones mathematicae, 209(1):275–327, 2017

Benedek Valk ´o and B ´alint Vir ´ag. The Sine β operator.Inventiones mathematicae, 209(1):275–327, 2017

2017
[18]

Wong.Asymptotic approximations of integrals, volume 34 ofClassics in Applied Math- ematics

R. Wong.Asymptotic approximations of integrals, volume 34 ofClassics in Applied Math- ematics. Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2001. Shengqi Qiu, DEPARTMENT OFMATHEMATICS, UNIVERSITY OFWISCONSIN– MADISON, 480 LINCOLN DR, MADISON, WI 537...

2001

Pith tools

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

[1] [14]

On the pair correlation function of the Sine β process

Yahui Qu and Benedek Valk ´o. On the pair correlation function of the Sine β process. https://arxiv.org/abs/2509.15446, 2025

arXiv 2025

[2] [3]

Moments of CβE field partition function, Sineβ correlations and stochastic zeta.https://arxiv.org/abs/2602.08739, 2026

Theodoros Assiotis and Joseph Najnudel. Moments of CβE field partition function, Sineβ correlations and stochastic zeta.https://arxiv.org/abs/2602.08739, 2026

arXiv 2026

[3] [5]

Long-Range Correlation of the Sineβ point Process

Laure Dumaz and Martin Malvy. Long-Range Correlation of the Sineβ point Process. https://arxiv.org/abs/2603.15289, 2026

arXiv 2026

[4] [1]

Stegun.Handbook of mathematical functions with formulas, graphs, and mathematical tables.Wiley, 1984

Milton Abramowitz and Irene A. Stegun.Handbook of mathematical functions with formulas, graphs, and mathematical tables.Wiley, 1984

1984

[5] [2]

Anderson, Alice Guionnet, and Ofer Zeitouni.Introduction to random matri- ces

Greg W. Anderson, Alice Guionnet, and Ofer Zeitouni.Introduction to random matri- ces. Cambridge University Press, 2009. 10

2009

[6] [4]

Universality of generalβ- ensembles.Duke Math

Paul Bourgade, L ´aszl´o Erd ˝os, and Horng-Tzer Yau. Universality of generalβ- ensembles.Duke Math. J., 163(6):1127–1190, 2014

2014

[7] [6]

Forrester

Peter J. Forrester. Selberg correlation integrals and the1/r 2 quantum many-body system.Nuclear Physics B, 388(3):671–699, 1992

1992

[8] [7]

Forrester

Peter J. Forrester. Addendum to ‘Selberg correlation integrals and the1/r 2 quantum many body system’.Nuclear Physics B, 416(1):377–385, 1994

1994