REVIEW 3 major objections 4 minor 38 references
CLT for Circular beta-Ensembles at High Temperature
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read At high temperature, the circular beta-ensemble's linear statistics satisfy a central limit theorem with an explicit covariance that interpolates between white noise and random-matrix noise.
desk verdict First explicit CLT at the beta/N crossover for circular beta ensembles, with an honest proof chain; but the advertised W2 rate is not obtained by the argument as written. 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 load-bearing object is the operator $L=A+2\pi\beta W$ acting on the Hilbert space $H=\{\psi\in L^2(\mathbb{T}):\psi'\in L^2(\mathbb{T}),\int\psi\,d\mu_\beta^V=0\}$ with inner product $\langle\phi,\psi\rangle_H=\int \phi'\psi'\,d\mu_\beta^V$. Here $A\phi=-\phi''-(\log\mu_\beta^V)'\phi'$ is a Sturm-Liouville operator and $W\phi=-H(\phi'\mu_\beta^V)$ involves the Hilbert transform $H$. Proposition 4.3 supplies a complete orthonormal eigenbasis $(\varphi_j)$ of $H$ with eigenvalues $\kappa_j\sim\alpha j^2$ and Lipschitz bounds $\|\varphi_j^{(k)}\|_{\mathrm{Lip}}\le C_k\kappa_j^{(k+1)/2}$. This spectral structure is what lets the proof decompose an arbitrary test function into eigenmodes, truncate the expansion with controllable error, and apply a normal-approximation theorem to the finite-dimensional projections.
What would settle it
For $V=0$, simulate the particle system at fixed $\beta$ and measure the variance of $N^{-1/2}\sum_i\cos(x_i)$; Theorem 1.2 predicts it approaches $1/(2(1+\beta))$, so a clear mismatch at large $N$ would refute the CLT. More directly, one can numerically diagonalize a discretization of $L$ for a non-constant smooth $V$ and check whether the eigenvalue ratios $\kappa_j/j^2$ converge to a positive constant and whether the stated eigenfunction Lipschitz bounds hold.
Extended reading notes
Core claim
The central discovery is that, in the high-temperature scaling where the inverse temperature is $2\beta/N$, the random measure $\nu_N=\sqrt{N}(\mu_N-\mu_\beta^V)$ has Gaussian macroscopic fluctuations with an explicit covariance operator. For $\beta>0$, $V\in C^{3,1}(\mathbb{T})$ and $\psi\in C^{2\gamma+1}(\mathbb{T})$, $\gamma\ge 2$, with $\int\psi\,d\mu_\beta^V=0$, Theorem 1.2 states that $\nu_N(\psi)=N^{-1/2}\sum_i\psi(x_i)$ converges in law to $\mathcal{N}(0,\sigma_\beta^V(\psi)^2)$ with $\sigma_\beta^V(\psi)^2=\langle\psi,L^{-1}\psi\rangle_H=\int \psi'(L^{-1}\psi)'\,d\mu_\beta^V$, and $W_2(\nu_N(\psi),\mathcal{N}(0,\sigma^2))\le C\sqrt{\log N}\,N^{-(\gamma-1)/(\gamma+1)}$. Here $\mu_\beta^V$ is the unique minimizer of $\beta E(\mu)+K(\mu|\mu_0^V)$, the operator $L$ is defined by $-L\phi=\phi''+2\pi\beta H(\mu_\beta^V\phi')+(\log\mu_\beta^V)'\phi'$, and $H$ is the Hilbert transform. In the unweighted case $V=0$, the variance is $\sigma_\beta^0(\psi)^2=2\sum_{k\ge1}\frac{1}{1+\beta/k}|\widehat{\psi}_k|^2$, which visibly interpolates between the $L^2$ norm and the $H^{1/2}$ semi-norm.
Load-bearing premise
The load-bearing assumption is that the operator $L$, including the Hilbert-transform part $W$, is self-adjoint on the Sobolev space $H$ and has a complete eigenbasis with eigenvalues growing like the square of their index and with derivatives of eigenfunctions controlled by powers of those eigenvalues. The rate proof depends on this spectral picture, so if any of it fails, the stated CLT and convergence rate do not follow.
Editorial extensions
If this is right
- At high temperature, the macroscopic fluctuations of the circular beta-ensemble are Gaussian with a covariance operator that interpolates continuously between the $L^2$ white-noise covariance at $\beta=0$ and the $H^{1/2}$ random-matrix covariance at $\beta=\infty$.
- The rate of convergence in the Wasserstein $W_2$ metric is $O(\sqrt{\log N}\,N^{-(\gamma-1)/(\gamma+1)})$ for $\psi\in C^{2\gamma+1}$, so smoother test functions yield faster convergence; for $C^\infty$ test functions the rate is $O(\sqrt{\log N}/N)$.
- For $V=0$, the limiting variance has the explicit Fourier formula $2\sum_{k\ge1}\frac{1}{1+\beta/k}|\widehat{\psi}_k|^2$, reducing to the classical formulas in the limits $\beta\to0$ and $\beta\to\infty$.
- The equilibrium measure $\mu_\beta^V$ is the unique minimizer of $\beta E(\mu)+K(\mu|\mu_0^V)$, has a density bounded above and below by positive constants, and is as smooth as the potential $V$.
- The concentration estimate of Theorem 1.4 gives exponential tail bounds for $W_1(\mu_N,\mu_\beta^V)$, which imply almost sure convergence of the empirical measure whenever $\beta\gg N^{-1}$.
Reading between the lines
- The same operator-level strategy would likely extend to Gaussian beta-ensembles on the real line and to general potentials, giving an explicit limiting variance in settings where only implicit central limit theorems were previously available.
- The variance formula suggests a two-parameter family of Gaussian fields parametrized by $\beta$; at finite $\beta$ the fluctuations are neither pure white noise nor pure random-matrix noise, which could be tested numerically as a systematic interpolation between the two known regimes.
- The paper leaves open whether the condition $\beta\inf\mu_\beta^V\to\infty$ can fail for some smooth non-constant potential; if it does fail, the $H^{1/2}$ universality at $\beta=\infty$ would break and the limiting variance would retain a dependence on $V$.
- A direct numerical check of the unweighted formula $\sigma_\beta^0(\psi)^2=2\sum_{k\ge1}\frac{1}{1+\beta/k}|\widehat{\psi}_k|^2$, mode by mode, would isolate the operator's spectrum from other sources of error and sharpen confidence in the whole proof.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the Circular beta-Ensemble on the torus in the high-temperature regime where the inverse temperature scales as beta/N, with an external potential V. The main object is the signed fluctuation field nu_N = sqrt(N)(mu_N - mu^V_beta) around the equilibrium measure mu^V_beta, which minimizes the functional beta E + relative entropy. Theorem 1.2 states a CLT for linear statistics nu_N(psi) with explicit variance sigma^V_beta(psi)^2 = <psi, L^{-1} psi>_H, where L is the Sturm-Liouville-type operator in (15), together with a quantitative Wasserstein-2 rate sqrt(log N) N^{-(gamma-1)/(gamma+1)}. The proof combines a concentration estimate in W1 (Theorem 1.4), obtained from a Coulomb transport inequality and an energy regularization, with the Stein-type normal approximation of Lambert-Ledoux-Webb and a spectral decomposition of L. Sections 2-8 develop the equilibrium measure properties, the concentration proof, the eigenbasis regularity, and the beta -> 0/infty behavior of the variance.
Significance. The CLT at high temperature is a natural and interesting result: it gives an explicit crossover covariance structure interpolating between the L2 noise of independent particles at beta = 0 and the H^{1/2} noise of the usual beta-ensemble at beta = infinity, and it provides quantitative convergence in W2. The overall strategy is coherent and builds on solid external tools, and the spectral analysis of the operator L is a genuine contribution. The concentration inequality of Theorem 1.4 is also of independent interest. However, the advertised W2 rate in Theorem 1.2 is not supported by the estimates as written, and there are gaps in the spectral self-adjointness argument; these affect a stated quantitative result, though the CLT itself appears salvageable.
major comments (3)
- [Section 4, proof of Theorem 1.2, Eqs. (58)-(60)] The stated W2 rate in Theorem 1.2 is not what the proof delivers. Taking d to be the integer part of N^{1/(4(gamma+1))}, estimate (58) gives W2(nu_N(psi), nu_N(psi[d])) <= C sqrt(log N) d^{-2(gamma-1)} = C sqrt(log N) N^{-(gamma-1)/(2(gamma+1))}. The dominant term in (60) is of the same order, or slower if the displayed Proposition 4.4 bound is used literally. Hence the proof establishes only W2 <= C sqrt(log N) N^{-(gamma-1)/(2(gamma+1))}, not the claimed N^{-(gamma-1)/(gamma+1)}. The weak convergence is unaffected, but the quantitative rate in the theorem and the related claim in Remark 1.1 must be corrected.
- [Section 6, Proposition 6.3] Self-adjointness of W on H is asserted without proof. The text states 'Since W is non-negative and self-adjoint on H', but only positivity of the quadratic form is shown in Lemma 6.1. Positivity of a quadratic form does not imply symmetry or self-adjointness for an unbounded operator, and this property is load-bearing for the spectral theorem and the min-max comparison that yield the eigenvalue asymptotics. Please provide a domain and a proof of self-adjointness, for instance by an integration by parts using H* = -H and the equilibrium equation, or cite a precise reference.
- [Section 4, Proposition 4.4 and Eq. (60)] There is an inconsistency in the d-dependence of the second error term. Proposition 4.4 states the second contribution as sqrt(log N) sqrt(sum kappa_j^2) sum kappa_j, which for kappa_j ~ alpha j^2 is of order sqrt(log N) d^{11/2}. Equation (60), however, writes sqrt(log N) d^4. The square-root of the bound in (55) naturally gives sqrt(log N) sqrt(sum kappa_j^2) sqrt(sum kappa_j), which is d^4; the displayed Proposition 4.4 is weaker. The two forms are not interchangeable as written, and this discrepancy affects the resulting rate in Theorem 1.2. Please reconcile the statement of Proposition 4.4 with its proof and with equation (60).
minor comments (4)
- [Proposition 4.3(c)] There is a stray period in the statement: 'kappa_j ~ alpha j^2 as j -> infinity. for every j >= 1' should be 'kappa_j ~ alpha j^2 as j -> infinity, for every j >= 1'.
- [Abstract and Introduction] There is a typo in the abstract: 'equilibrium mesure' should be 'equilibrium measure'. Similar small typos appear elsewhere and should be cleaned up.
- [Remark 1.1] The remark that the rate for C^infty test functions is sqrt(log N)/N is based on the currently stated exponent; if the correct rate is the weaker one identified above, this remark must be updated accordingly.
- [Section 3, proof of Theorem 1.4] The citation to [Chafai, Hardy, Maida, Theorem 1.1] and 'the discussion below' could be made more precise by specifying the exact form of the Coulomb transport inequality used in (41).
Circularity Check
No significant circularity: the CLT is derived from independent external tools (Sturm-Liouville theory, Lambert-Ledoux-Webb normal approximation, Chafai-Hardy-Maida transport inequality) and from internally proved spectral estimates, with no fitted parameter or target assumption smuggled into the inputs.
full rationale
The derivation chain is self-contained and does not reduce to its inputs. The equilibrium measure mu_V_beta is defined as the minimizer of F_beta, which is an independent variational object justified by classical large-deviation and convexity results; the CLT variance is then expressed in terms of the operator L built from this equilibrium measure, but L is not defined using the target fluctuation law. The spectral Proposition 4.3 is proven, not assumed: (a)-(c) follow from Sturm-Liouville theory for A plus positivity of W (Lemma 6.1, Lemma 6.2, Proposition 6.3), and (d) follows from an internal bootstrap on the eigenfunction equation; no eigen-gap or eigenvalue-growth condition equivalent to the CLT is imposed. The concentration estimate (Theorem 1.4) is derived from the Coulomb transport inequality of Chafai-Hardy-Maida and a regularization lemma, not from the CLT being proven. The normal approximation step (Theorem 4.5) is cited to Lambert-Ledoux-Webb, an external published result whose stated assumptions do not include the high-temperature CLT; the citation is load-bearing but it is legitimate independent support, not a self-referential uniqueness or ansatz citation. The identity eta_infty = sigma_V_beta(psi)^2 is verified by the spectral computation in Proposition 6.3, so the variance formula is not assumed by construction. There are no fitted parameters, no quantity called a prediction that was used as an input, and no renaming of a known empirical pattern. The advertised W2 exponent in Theorem 1.2 may or may not follow from the displayed estimates (a possible correctness concern), but that is unrelated to circularity: nothing in the proof assumes the target convergence or variance.
Assumptions & free parameters
assumptions (7)
- standard math The logarithmic energy E and relative entropy K have compact level sets and are strictly convex, so F^V_beta has a unique minimizer.
- domain assumption Large deviation principle for the empirical measure at speed beta N with rate F^V_beta - F^V_beta(mu^V_beta), as in Berman (2018) and Garcia-Zelada (2018).
- domain assumption Coulomb transport inequality of Chafai-Hardy-Maida (2018): W1(mu,nu)^2 <= 4 pi E(mu-nu) for probability measures with finite logarithmic energy.
- domain assumption Stein normal approximation theorem of Lambert-Ledoux-Webb (2017), quoted as Theorem 4.5: W2(F,N(0,Id)) <= sqrt(E||F+K^{-1}LF||^2) + sqrt(E||I-K^{-1}Gamma||^2).
- standard math Sturm-Liouville spectral theory and Weyl's law for A = -(phi' mu)'/mu on T: orthonormal eigenbasis of H with eigenvalues lambda_j ~ alpha j^2.
- standard math Hilbert transform properties: H(e^{ikx}) = i sgn(k) e^{ikx}, H* = -H, H is an isometry on L2_0(T), and (Hf)' = H(f').
- standard math Sobolev embedding H^{m+1}(T) subset C^{m,1/2}(T).
Cite this review
Pith. "Pith review of CLT for Circular beta-Ensembles at High Temperature." pith.science (2026). https://pith.science/paper/OLBURSUI
@misc{pith2026190901142,
author = {Pith},
title = {Pith review of: CLT for Circular beta-Ensembles at High Temperature},
year = {2026},
howpublished = {\url{https://pith.science/paper/OLBURSUI}},
note = {Machine review of arXiv:1909.01142}
}
read the original abstract
We consider the macroscopic large N limit of the Circular beta-Ensemble at high temperature, and its weighted version as well, in the regime where the inverse temperature scales as beta/N for some parameter beta>0. More precisely, in the large N limit, the equilibrium measure of this particle system is described as the unique minimizer of a functional which interpolates between the relative entropy (beta=0) and the weighted logarithmic energy (beta=\infty). More precisely, we provide subGaussian concentration estimates in the W1 metric for the deviations of the empirical measure to this equilibrium mesure. The purpose of this work is to show that the fluctuation of the empirical measure around the equilibrium measure converges towards a Gaussian field whose covariance structure interpolates between the Lebesgue L^2 (beta=0) and the Sobolev H^{1/2} (beta=\infty) norms. We furthermore obtain a rate of convergence for the fluctuations in the W_2 metric. Our proof uses the normal approximation result of Lambert, Ledoux and Webb [2017] the Coulomb transport inequality of Chafai, Hardy, Maida [2018] and a spectral analysis for the operator associated with the limiting covariance structure.
Reference graph
Works this paper leans on
-
[1]
G. Akemann and S.-S. Byun. The high temperature crossover for general 2D C oulomb gases. J. Stat. Phys., 175 0 (6): 0 1043--1065, 2019
work page 2019
-
[2]
R. Allez , J.-P. Bouchaud , and A. Guionnet . Invariant Beta Ensembles and the Gauss-Wigner Crossover . Physical Review Letters, 109 0 (9): 0 094102, Aug. 2012. URL http://dx.doi.org/10.1103/PhysRevLett.109.094102
-
[3]
G. W. Anderson, A. Guionnet, and O. Zeitouni. An Introduction to Random Matrices, volume 118. Cambridge Studies in Advanced Mathematics, 2010
work page 2010
-
[4]
F. Benaych-Georges and S. P \'e ch \'e . Poisson statistics for matrix ensembles at large temperature. J. Stat. Phys., 161 0 (3): 0 633--656, 2015
work page 2015
-
[5]
R. J. Berman. On large deviations for G ibbs measures, mean energy and gamma-convergence. Constructive Approximation, 48 0 (1): 0 3--30, 2018
work page 2018
-
[6]
A. Borodin and S. Serfaty. Renormalized energy concentration in random matrices. Comm. Math. Phys., 320 0 (1): 0 199--244, 2013
work page 2013
-
[7]
B. M. Brown, M. S. P. Eastham, and K. M. Schmidt. Periodic differential operators, volume 230 of Operator Theory: Advances and Applications. Birkh\" a user/Springer Basel AG, Basel, 2013. ISBN 978-3-0348-0527-8; 978-3-0348-0528-5. doi:10.1007/978-3-0348-0528-5. URL https://mathscinet.ams.org/mathscinet-getitem?mr=2978285
-
[8]
D. Chafa\" i , A. Hardy, and M. Ma\" i da. Concentration for C oulomb gases and C oulomb transport inequalities. Journal of Functional Analysis, 275 0 (6): 0 1447 -- 1483, 2018. ISSN 0022-1236. doi:https://doi.org/10.1016/j.jfa.2018.06.004. URL http://www.sciencedirect.com/science/article/pii/S0022123618302209
Show all 38 references
-
[9]
Chatterjee
S. Chatterjee. Fluctuations of eigenvalues and second order P oincar \'e inequalities. Probab. Theory Related Fields, 143: 0 1--40, 2009
2009
-
[10]
Deift, A
P. Deift, A. Its, and I. Krasovsky. Toeplitz matrices and T oeplitz determinants under the impetus of the I sing model: some history and some recent results. Comm. Pure Appl. Math., 66 0 (9): 0 1360--1438, 2013. ISSN 0010-3640. doi:10.1002/cpa.21467. URL https://mathscinet.ams...
2013 doi
-
[11]
Dembo and O
A. Dembo and O. Zeitouni. Large deviations techniques and applications, volume 38 of Stochastic Modelling and Applied Probability. Springer-Verlag, Berlin, 2010. URL http://dx.doi.org/10.1007/978-3-642-03311-7. Corrected reprint of the second (1998) edition
1998 doi
-
[12]
D \"o bler and M
C. D \"o bler and M. Stolz. Stein's method and the multivariate CLT for traces of powers on the classical compact groups. Electron. J. Probab., 16: 0 2375--2405, 2011
2011
-
[13]
D \"o bler and M
C. D \"o bler and M. Stolz. A quantitative central limit theorem for linear statistics of random matrix eigenvalues. J. Theor. Probab., 27: 0 945--953, 2014
2014
-
[14]
Dumitriu and A
I. Dumitriu and A. Edelman. Matrix models for beta ensembles. J. Math. Phy., 43 0 (11): 0 5830--5847, 2002
2002
-
[15]
J. Fulman. Stein's method, heat kernel, and traces of powers of elements of compact lie groups. Electron. J. Probab., 17 0 (66): 0 16 pp., 2012
2012
-
[16]
Garc \' a-Zelada
D. Garc \' a-Zelada. A large deviation principle for empirical measures on polish spaces: Application to singular gibbs measures on manifolds. Preprint arXiv:1703.02680. To appear in Annales de l'Institut Henri Poincar \'e . , 2018
2018 arXiv
-
[17]
Guionnet and T
A. Guionnet and T. Bodineau. About the stationary states of vortex systems. Ann. Inst. H. Poincar\'e Probab. Statist., 53 0 (2): 0 205--237, 1999
1999
-
[18]
Hiai and D
F. Hiai and D. Petz. The semicircle law, free random variables and entropy, volume 77 of Mathematical Surveys and Monographs. 2000
2000
-
[19]
Johansson
K. Johansson. On S zeg \"o 's asymptotic formula for T oeplitz determinants and generalizations. Bull. Sci. Math., 112 0 (3): 0 257--304, 1988
1988
-
[20]
T. Johnson. Exchangeable pairs, switchings, and random regular graphs. Electron. J. Combin., 22 0 (1): 0 1--33, 2015
2015
-
[21]
T. Kato. Perturbation theory for linear operators. Reprint of the 1980 edition, Springer-Verlag, Berlin, 1995
1980
-
[22]
G. Lambert. Mesoscopic central limit theorem for the circular beta-ensembles and applications. Preprint arXiv:1902.06611, 2019
1902 arXiv
-
[23]
Lambert, M
G. Lambert, M. Ledoux, and C. Webb. Quantitative normal approximation of linear statistics of -ensembles. Ann. of Probab. Vol. 47, No. 5, 2619--2685, 2019
2019
-
[24]
Ma \" i da and \'E
M. Ma \" i da and \'E . Maurel-Segala. Free transport-entropy inequalities for non-convex potentials and application to concentration for random matrices. Probab. Theory Related Fields, 159 0 (1-2): 0 329--356, 2014. URL http://dx.doi.org/10.1007/s00440-013-0508-x
2014 doi
-
[25]
V. A. Marchenko. Sturm- L iouville operators and applications . AMS Chelsea Publishing, Providence, RI, revised edition, 2011. ISBN 978-0-8218-5316-0. doi:10.1090/chel/373. URL https://mathscinet.ams.org/mathscinet-getitem?mr=2798059
2011 doi
-
[26]
Nakano and K
F. Nakano and K. D. Trinh. Gaussian beta ensembles at high temperature: eigenvalue fluctuations and bulk statistics. J. Stat. Phys., 173 0 (2): 0 296--321, 2018
2018
-
[27]
Nakano and K
F. Nakano and K. D. Trinh. Poisson statistics for beta ensembles on the real line at high temperature. Preprint arXiv:1910.00766, 2019
1910 arXiv
-
[28]
C. Pakzad. Poisson statistics at the edge of gaussian beta-ensembles at high temperature. Preprint arXiv:1804.08214, 2018
2018 arXiv
-
[29]
C. Pakzad. Large deviations principle for the largest eigenvalue of the gaussian --ensemble at high temperature. J. Theor. Probab., pages 1--19, 2019 a
2019
-
[30]
C. Pakzad. Extremes of chi triangular array from the gaussian βββ-ensemble at high temperature. Preprint arXiv:1903.02103, 2019 b
1903 arXiv
-
[31]
Reed and B
M. Reed and B. Simon. Methods of modern mathematical physics, volume IV. Analysis of operators. Academic Press [Harcourt Brace Jovanovich, Publishers], New York-London, 1978
1978
-
[32]
N. Ross. Fundamentals of stein's method. Probab. Surveys, 8: 0 210--293, 2011
2011
-
[33]
Rougerie and S
N. Rougerie and S. Serfaty. Higher-dimensional C oulomb gases and renormalized energy functionals. Comm. Pure Appl. Math., 69 0 (3): 0 519--605, 2016. URL http://dx.doi.org/10.1002/cpa.21570
2016 doi
-
[34]
E. B. Saff and V. Totik. Logarithmic potentials with external fields, volume 316 of Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Springer-Verlag, Berlin, 1997. URL http://dx.doi.org/10.1007/978-3-662-03329-6. Appendix B by Th...
1997 doi
-
[35]
B. Simon. Orthogonal polynomials on the unit circle. P art 1 and Part2 , volume 54 of American Mathematical Society Colloquium Publications. American Mathematical Society, Providence, RI, 2005. ISBN 0-8218-3675-7. doi:10.1090/coll/054.2/01. URL https://mathscinet.ams.org/maths...
2005 doi
-
[36]
H. Spohn. Generalized gibbs ensembles of the classical toda chain. Journal of Statistical Physics, May 2019. ISSN 1572-9613. doi:10.1007/s10955-019-02320-5. URL http://dx.doi.org/10.1007/s10955-019-02320-5
2019 doi
-
[37]
K. D. Trinh. Global spectrum fluctuations for gaussian beta ensembles: a martingale approach. J. Theor. Probab., pages 1--18, 2017
2017
-
[38]
C. Webb. Linear statistics of the circular -ensemble, stein's method, and circular dyson brownian motion. Electron. J. Probab., 21 0 (25): 0 16 pp., 2016
2016
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.