Pith. sign in

REVIEW 3 major objections 7 minor 23 references

Distributions of consecutive level spacings of circular unitary ensemble and their ratio: finite-size corrections and Riemann $\zeta$ zeros

T0 review · 3 major / 7 minor · reviewed 2026-08-06 · deepseek-v4-flash

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

Pith's one-line read For Haar-random unitary matrices, the gap-ratio distribution has no $O(N^{-2})$ finite-size correction; the first correction is $O(N^{-4})$, and this cancellation explains the $(\log(T/2\pi))^{-3}$ deviation of Riemann zeta zero gap ratios.

desk verdict The paper's claim that the O(N^{-2}) CUE gap-ratio correction vanishes is a strong, non-obvious result supported only numerically, and the zeta explanation is heuristic; still, it deserves a careful referee. read the letter →

arxiv 2507.10193 v1 pith:FP5OZEVT submitted 2025-07-14 math-ph cond-mat.dis-nnhep-thmath.MPmath.PR

Shinsuke M. Nishigaki This is my paper · ORCID
classification math-phcond-mat.dis-nnhep-thmath.MPmath.PR MSC 60B2015B5211M06
keywords circularunitaryensemblegap-ratiodistributionTracy-WidomequationsJanossydensityfinite-sizecorrectionsRiemannzetazerossinekernelconsecutivelevelspacings
topics The Riemann Hypothesis
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 establish a precise finite-size statement for the circular unitary ensemble: the ratio of consecutive level spacings, unlike the spacings themselves, has no correction of order $N^{-2}$ relative to the universal sine-kernel limit, so the first correction is $O(N^{-4})$. This matters because the gap ratio is the commonly used statistic for detecting quantum chaos, and a statistic that skips the first finite-size order can expose deeper corrections in a spectrum. The author further argues that this cancellation is exactly what makes Riemann zeta zero data show deviations from random-matrix predictions that scale as $(\log(T/2\pi))^{-3}$. A sympathetic reader would take the paper as giving the analytic CUE distributions and a new explanation for the zeta-zero finite-size scaling.

What carries the argument

The central object is the CUE Janossy density $\tilde J_1(0; [a_1,a_2])$, the probability that an interval contains no eigenphase except one fixed at $0$; it is a Fredholm determinant of the conditioned CUE kernel. The machinery is the Tracy-Widom system of nonlinear PDEs that gives the dependence of this determinant on the interval endpoints, supplemented by boundary conditions from the Neumann expansion. All spectral statistics studied here, including the gap-ratio distribution, are integrals of derivatives of this Janossy density, and the system's $N\mapsto -N$ symmetry after unfolding makes every derived distribution analytic in $1/N^2$, which is why the absence of the $N^{-2}$ term in $P_r$ is a meaningful cancellation.

What would settle it

A direct symbolic or very-high-precision numerical computation of the coefficient $P_r^{(2)}(r)$ from the large-$N$ expansion of the Tracy-Widom system: if $N^2(P_r(r)-P_r^{(0)}(r))$ has a nonzero scaling limit as $N\to\infty$, the claimed cancellation is false. For the zeta part, a gap-ratio histogram at larger $T$ whose deviation scales as $N_e^{-2}$ rather than $N_e^{-3}$ would falsify the kernel-based explanation.

Watch

Extended reading notes

Core claim

The paper derives, from the Tracy-Widom system for the CUE_J\'anossy density, the nearest-neighbor spacing distribution $P_{nn}(t)$, the joint distribution $P_c(a,b)$ of two consecutive spacings, and the gap-ratio distribution $P_r(r)$. Its central finding is that as $N\to\infty$, $P_c$ has the expected leading finite-size correction of order $N^{-2}$, but $P_r$ does not: the coefficient $P_r^{(2)}(r)$ vanishes, so $P_r(r)=P_r^{(0)}(r)+N^{-4}P_r^{(4)}(r)+O(N^{-6})$. The paper then uses the conjectured Riemann-zero kernel, which matches the CUE kernel to $O(N_e^{-2})$ and differs at $O(N_e^{-3})$, together with the effective size $N_e(T)\sim\log(T/2\pi)$, to show that the gap-ratio deviation of the zeta zeros should scale as $N_e^{-3}=(\log(T/2\pi))^{-3}$, consistent with histograms of hundreds of millions of zeros.

Load-bearing premise

On the zeta side, the load-bearing premise is that Riemann zero spacings are described by the conjectured finite-size kernel (41) with effective size $N_e(T)=(\log(T/2\pi))/\sqrt{12\Lambda}$; if this identification is wrong the $(\log T)^{-3}$ explanation fails, while the CUE cancellation claim itself stands independently.

Editorial extensions

If this is right

  • For CUE_N, $P_r^{(2)}(r)=0$ means the gap-ratio statistic is insensitive to the leading finite-size part of the kernel, so comparisons of spectra with theory can be pushed to $O(N^{-4})$.
  • Riemann zeta zero gap-ratio deviations from the sine kernel should be dominated by $N_e^{-3}\propto(\log(T/2\pi))^{-3}$ rather than by $N_e^{-2}$, which the paper's histograms support and larger data sets can test more sharply.
  • The joint distribution $P_c(a,b)$ retains the $O(N^{-2})$ correction, so comparing $P_c$ and $P_r$ between a given spectrum and CUE_N distinguishes kernels that match to first order from those matching to second order.
  • The Tracy-Widom/J\'anossy method yields explicit finite-$N$ expressions for $P_{nn}(t)$ and $P_c(a,b)$, providing benchmarks for numerical spectra of quantum-chaotic systems.
  • Because the would-be $N^{-2}$ term vanishes, any observed $N^{-2}$ deviation in a gap-ratio measurement would signal a kernel genuinely different from CUE at that order, rather than mere finite-size noise.

Reading between the lines

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

  • The mechanism of the cancellation is left open by the paper; a likely target for a proof is an algebraic identity in the Tracy-Widom system showing that $P_r^{(2)}(r)$ integrates to zero for all $r$, which would turn a numerical observation into a theorem.
  • The same ratio-based cancellation may occur for other circular ensembles with integrable kernels sharing the $N\mapsto -N$ symmetry, giving a testable prediction for circular orthogonal and symplectic ensembles and for finite-$N$ spectra of many-body systems.
  • Because the zeta kernel differs from CUE at $O(N_e^{-3})$, the gap-ratio statistic is a natural place to search for arithmetic effects in zeta-zero data at higher heights, provided samples are large enough to resolve the $N_e^{-3}$ shape.
  • The method could be extended to ratios of next-nearest spacings or higher-order gap ratios, where analogous cancellations may appear at higher orders and would show how generic the $N^{-2}$ suppression is.
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

3 major / 7 minor

Summary. The paper develops a Tracy-Widom (TW) PDE framework for the Jánossy density of the circular unitary ensemble CUE_N and uses it to compute the joint distribution of two consecutive spacings P_c(a,b), the nearest-neighbor spacing distribution P_nn(t), and the gap-ratio distribution P_r(r). The central claim is that the leading finite-N correction to P_r relative to the sine-kernel limit is O(N^{-4}), because the would-be O(N^{-2}) correction vanishes exactly; this is supported numerically by an N^4 scaling collapse for N=8..16. The paper then applies the same logic to the Riemann zeta zeros, arguing that the observed (log(T/2π))^{-3} scaling of the gap-ratio deviation follows from the conjectured finite-size kernel (41)–(42) for zeta zeros combined with the CUE_N cancellation. The manuscript includes a detailed derivation of the TW system, high-precision cross-checks against Nyström quadrature, and a Mathematica notebook.

Significance. If the claimed cancellation is correct, the result is of genuine interest: the gap ratio becomes a sharper probe of finite-size effects than the spacing distributions, since its leading correction is pushed to O(N^{-4}). The paper also provides a plausible mechanism linking this cancellation to the observed zeta-zero gap-ratio scaling. Strengths include the systematic TW derivation, the internal cross-check against Nyström quadrature at the 10^{-36} level for N=10, the availability of a reproducible notebook, and the clean scaling collapse in Fig. 3. The main weakness is that the load-bearing cancellation P_r^{(2)}(r)=0 is asserted from numerical evidence only, without an analytic derivation or error bound; the zeta conclusion is therefore conditional on an unproved structural fact, and the fitted exponent in Fig. 6 is not statistically tied to the claimed value -3.

major comments (3)
  1. [§4.4, Fig. 3 and Eqs. (39)–(40)] The central claim P_r^{(2)}(r)=0 is not derived. The text states that the cancellation mechanism is 'obscure from the structure of the PDEs', and the only support is the N^4 scaling collapse of five curves with N=8,10,12,14,16. Since P_r is defined in Eq. (40) as an integral of P_c, and since Fig. 2 already shows a well-defined N^{-2} correction to P_c, it should be possible to compute P_r^{(2)}(r) = ∫_0^{2π/(1+r)} b P_c^{(2)}(-rb,b) db from the N^{-2} expansion of the TW system (32)–(33). Please provide that computation, or a rigorous argument showing the integral vanishes, and quantify the contamination from N^{-4} and N^{-6} terms at the plotted values of N.
  2. [§5.2, Fig. 6] The empirical support for the claimed exponent -3 is weaker than stated. The optimal fit is reported as 0.1896 N_e^{-3.081}, i.e. a fitted exponent of -3.081, with no uncertainties and only six points. This does not distinguish -3.081 from -3, and it is not a test of the predicted exact exponent. Please fit with the slope fixed at -3 and report the residuals or a goodness-of-fit statistic, or provide error bars; otherwise the conclusion that the data confirm (log T)^{-3} scaling is not quantitatively supported.
  3. [§5.2, Fig. 5 and the footnote on p. 14] The zeta argument is presented as a consequence of the CUE cancellation plus the conjectured kernel (41)–(42), but no predicted coefficient for the O(N_e^{-3}) gap-ratio correction is computed from that kernel. Fig. 5 shows only noisy histograms, and the footnote simultaneously claims a match to CUE_N up to O(N_e^{-4}) while incorporating an O(N_e^{-3}) term, which is confusing. Please derive the explicit N_e^{-3} contribution to P_{RZ,r} from (41)–(42) and compare it to the histograms with a quantitative test, and clarify the order to which the adjusted kernel is matched.
minor comments (7)
  1. [§2, Eq. (2)] The unfolded variables are denoted by the same symbols x,y as the original eigenphases; please state explicitly that after the map x ↦ (2/N)x the variables are rescaled, to avoid ambiguity.
  2. [§4.1] The text 'WorkingPrecision → 5 MachinePrecision' is unclear; specify the actual precision setting used in the NDSolve integrations.
  3. [§4.4] In the sentence introducing the 1/N^2 expansion of P_r, the variable t is used for P_r^{(2)}(t) and P_r^{(4)}(t); the argument should be r throughout.
  4. [§5.2 and Fig. 6] The fit in Fig. 6 is described as an 'optimal linear fit'; on a log-log plot this is a power-law fit, not a linear fit, and the wording should be corrected.
  5. [§5.2, footnote] The footnote appears internally inconsistent: it claims the kernel can be adjusted to match the CUE_N expansion up to O(N_e^{-4}) while also incorporating an O(N_e^{-3}) term with ar α = 1 + Q/(√3 Λ^{3/2} N_e); please clarify whether the matching order is O(N_e^{-3}) or O(N_e^{-4}).
  6. [References [8] and [9]] Ref. [8] is given as Phys. Rev. E 110, 084101 (2013), but the cited article appears to be Phys. Rev. Lett. 110, 084101 (2013); Ref. [9] lists 'P. Vivo, and E. Vivo' as authors, which looks like an error. Please verify the author lists and journal names.
  7. [Throughout] There are several minor typos, including 'taylored' in §4.1 and 'deviations of from' in §5.2; a careful proofreading pass is needed.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular derivation: the CUE gap-ratio calculation is self-contained and cross-checked; the Riemann-zeta scaling analysis depends on an external conjectured kernel, not on fitted outputs of this paper.

full rationale

The central CUE result is obtained by solving the Tracy-Widom system (32)-(33), whose derivation is reproduced in Sections 3 and 4.1 from the Tracy-Widom theorem [14], with boundary conditions (33) from Neumann expansion. The numerical Jánossy densities are independently verified against Nyström quadrature (34), so the calculation does not reduce to a fitted input. The claimed O(N^{-2}) cancellation in the gap-ratio distribution is presented as an observed numerical scaling (Fig. 3) with the paper explicitly stating that the mechanism remains obscure; this is an admitted limitation of analytic support, not a disguised use of the target result as an assumption. Self-citations [12,13] supply the TW-for-Jánossy framework and the sine-kernel limit, but the framework equations are rederived here and the sine-kernel limit is a standard external result; therefore the self-citations are not load-bearing. The Riemann-zeta section imports the conjectured kernel (41) and effective size (42) from Bogomolny et al. [3] and Forrester-Mays [4], and compares with independent Odlyzko/LMFDB zero data; the N_e^{-3} scaling is confirmed by a free power-law fit (0.1896 N_e^{-3.081}), not forced by any parameter fitted in this paper. No equation is defined in terms of the quantity it purports to predict, and no fitted parameter is renamed as a prediction. The only concerns are the absence of an analytic derivation of P_r^{(2)}=0 and reliance on the conjectured zeta kernel, which are correctness/rigor risks, not circularity.

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

No invented entities. The central CUE claim rests on the standard TW/PDE machinery plus the author's prior framework, cross-checked internally. The zeta explanation rests on an external conjectured kernel and effective N_e, which are the main domain assumptions. The only fitted quantity is the diagnostic line in Fig. 6, not used to derive the claimed scaling.

free parameters (1)
  • Fit exponent and coefficient for mean gap-ratio deviation = 0.1896 N_e^{-3.081}
    Fig. 6 optimal linear fit to six zeta-zero data points; the claimed scaling exponent is -3, so this is a diagnostic fit, not an input parameter.
assumptions (5)
  • standard math Tracy-Widom theorem for integrable kernels (Eqs. 5-6) as stated in [14], with the paper's correction to the exponential variant.
    Central machinery; quoted from literature and modified; correctness supported by cross-check with Nyström quadrature (Sec. 4.1).
  • standard math The CUE_N kernel (2) and its unfolding expansion in powers of 1/N^2.
    Known DPP kernel for Haar-distributed U(N); the expansion is used to analyze the scaling of corrections.
  • standard math Gaudin-Mehta theorem connecting conditional gap probabilities to Fredholm determinants (Eq. 9).
    Basis for expressing the Janossy density as a Fredholm determinant.
  • domain assumption Conjectured Riemann zeta kernel KRZ (Eq. 41) and effective matrix size N_e(T) (Eq. 42), taken from [3,5].
    The explanation of the (log T)^{-3} scaling relies on this conjectured kernel; if the kernel identification is wrong, the explanation fails. Stated in Sec. 5.1-5.2.
  • standard math The distribution identities P_c(a1,a2) = -∂^2 J1/∂a1∂a2 (Eq. 39) and P_r(r) integral (Eq. 40).
    Standard representations used to derive the gap-ratio distribution from the Janossy density.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Distributions of consecutive level spacings of circular unitary ensemble and their ratio: finite-size corrections and Riemann $\zeta$ zeros." pith.science (2026). https://pith.science/paper/FP5OZEVT

@misc{pith2026250710193,
  author       = {Pith},
  title        = {Pith review of: Distributions of consecutive level spacings of circular unitary ensemble and their ratio: finite-size corrections and Riemann $\zeta$ zeros},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FP5OZEVT}},
  note         = {Machine review of arXiv:2507.10193}
}
abstract

We compute the joint distribution of two consecutive eigenphase spacings and their ratio for Haar-distributed $\mathrm{U}(N)$ matrices (the circular unitary ensemble) using our framework for J\'{a}nossy densities in random matrix theory, formulated via the Tracy-Widom system of nonlinear PDEs. Our result shows that the leading finite-$N$ correction in the gap-ratio distribution relative to the universal sine-kernel limit is of $\mathcal{O}(N^{-4})$, reflecting a nontrivial cancellation of the $\mathcal{O}(N^{-2})$ part present in the joint distributions of consecutive spacings. This finding suggests the potential to extract subtle finite-size corrections from the energy spectra of quantum-chaotic systems and explains why the deviation of the gap-ratio distribution of the Riemann zeta zeros $\{1/2+i\gamma_n\}, \gamma_n\approx T\gg1$ from the sine-kernel prediction scales as $\left(\log(T/2\pi)\right)^{-3}$.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

23 extracted references · 22 canonical work pages

  1. [1]

    F. J. Dyson, Statistical theory of the energy levels of complex systems. I , J. Math. Phys., 3, 140 (1962)

    work page 1962

  2. [2]

    J. P. Keating and N. C. Snaith, Random Matrix Theory and ζ(1/2 + it), Commun. Math. Phys., 214, 57 (2000)

    work page 2000

  3. [3]

    Bogomolny, O

    E. Bogomolny, O. Bohigas, P. Leboeuf, and A. G. Monastra, On the spacing distribution of the Riemann zeros: corrections to the asymptotic result , J. Phys. A, 39, 10743 (2006)

    work page 2006

  4. [4]

    P. J. Forrester and A. Mays, Finite size corrections in random matrix theory and Odlyzko’s data set for the Riemann zeros , Proc. R. Soc. A 471, 20150436 (2015)

    work page 2015

  5. [5]

    Bornemann, P

    F. Bornemann, P. J. Forrester, and A. Mays, Finite size effects for spacing distributions in random matrix theory: circular ensembles and Riemann zeros , Stud. Appl. Math. 138, 401 (2017)

    work page 2017

  6. [6]

    P. J. Forrester and B.-J. Shen, Finite size corrections in the bulk for circular β ensembles, arXiv: 2505.09865 [math-ph]

    arXiv

  7. [7]

    Oganesyan and D

    V. Oganesyan and D. A. Huse, Localization of interacting fermions at high temperature , Phys. Rev. B, 75, 155111 (2007)

    work page 2007

  8. [8]

    Y. Y. Atas, E. Bogomolny, O. Giraud, and G. Roux, Distribution of the ratio of consecutive level spacings in random matrix ensembles , Phys. Rev. E, 110, 084101 (2013)

    work page 2013

Show all 23 references
  1. [9]

    Y. Y. Atas, E. Bogomolny, O. Giraud, P. Vivo, and E. Vivo, Joint probability densities of level spacing ratios in random matrices , J. Phys. A, 46, 355204 (2013)

    2013

  2. [10]

    S. M. Nishigaki, M. Giordano, T. G. Kov´ acs, and F. Pittler, Critical statistics at the mobility edge of QCD Dirac spectra, PoS LATTICE 2013, 018 (2013)

    2013

  3. [11]

    P. J. Forrester and A. M. Odlyzko, Gaussian unitary ensemble eigenvalues and Riemann ζ function zeros: A nonlinear equation for a new statistic , Phys. Rev. E, 54, R4493 (1996)

    1996

  4. [12]

    S. M. Nishigaki, Tracy-Widom method for J´ anossy density and joint distribution of extremal eigenvalues of random matrices , Prog. Theor. Exp. Phys. 2021, 113A01 (2021)

    2021

  5. [13]

    S. M. Nishigaki, Distributions of consecutive level spacings of Gaussian unitary ensemble and their ratio: ab initio derivation , Prog. Theor. Exp. Phys. 2024, 081A01 (2024)

    2024

  6. [14]

    C. A. Tracy and H. Widom, Fredholm determinants, differential equations and matrix models, Commun. Math. Phys., 163, 33 (1994)

    1994

  7. [15]

    P. J. Forrester, Log-gases and random matrices , Princeton University Press (2010), chap. 9.1

    2010

  8. [16]

    S. M. Nishigaki, Level spacing distribution of critical random matrix ensembles, Phys. Rev. E, 58, R6915 (1998)

    1998

  9. [17]

    S. M. Nishigaki, Level spacings at the metal-insulator transition in the Anderson Hamiltonians and multifractal random matrix ensembles , Phys. Rev. E, 59, 2853 (1999)

    1999

  10. [18]

    Bornemann, On the numerical evaluation of Fredholm determinants , Math

    F. Bornemann, On the numerical evaluation of Fredholm determinants , Math. Comp., 79, 871 (2010)

    2010

  11. [19]

    Rudnick and P

    Z. Rudnick and P. Sarnak, Zeros of principal L-functions and random matrix theory , Duke Math. J. 81, 269 (1996)

    1996

  12. [20]

    H. L. Montgomery, The pair correlation of zeros of the zeta function , in Analytic Number Theory. Proceedings of the 24th Symposium on Pure Mathematics, ed. H. G. Diamond (Amer. Math. Soc., Providence, RI, 1973), p. 181

    1973

  13. [21]

    A. M. Odlyzko, On the distribution of spacings between zeros of the zeta function , Math. Comp. 48, 273 (1987)

    1987

  14. [22]

    A. M. Odlyzko, The 1022-nd zero of the Riemann zeta function , in Dynamical, spectral, and arithmetic zeta functions , ed. M. van Frankenhuisen and M.I. Lapidus, Contemporary Math. 290 (Amer. Math. Soc., Providence, RI, 2001), p. 139

    2001

  15. [23]

    (Available at: https://www.lmfdb.org/zeros/zeta/, date last accessed: July, 2024)

    The LMFDB Collaboration, The L-functions and modular forms database (2024). (Available at: https://www.lmfdb.org/zeros/zeta/, date last accessed: July, 2024). 15/15

    2024

Pith tools

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