alphaFeedback

WorksMathematics

An N−7/10N^{-7/10} Bound for the Finite-Size Bias of the Sherrington-Kirkpatrick Free Energy

AI contributions
  • OpenAI ChatGPT 5.6
  • OpenAI ChatGPT 6
  • OpenAI ChatGPT 6.1
  • Anthropic Claude Opus 5.5
  • Anthropic Claude Fable 5
  • Anthropic Claude Fable 5.1
View details
Ancillary data
Proof registry entry
View

Abstract

For every fixed inverse temperature β>1\beta>1, we prove that the expected free energy per spin of the zero-field Gaussian Sherrington-Kirkpatrick model differs from its Parisi limit by at most CβN−7/10C_\beta N^{-7/10}. This improves the exponent 2/32/3 of our earlier work (arXiv:2609.09103). We also compare the bound with the finite-size correction of order N−3/4N^{-3/4} predicted by Franz, Parisi, and Ricci-Tersenghi, and with the strict improvement over 2/32/3 that would follow from the ideas in OpenAI's recently released manuscripts.

Authors’ Note

With the exception of this note and small edits to the abstract, the entire document is machine written. We have taken care to supervise this writing to ensure that the proofs meet a minimum standard of readability and that previous literature is appropriately cited. However, the resulting document remains lacking in exposition and overall coherence. We are nonetheless releasing this manuscript on Hexagon to disseminate the result as quickly as possible, because we believe it may be of interest to the mathematics and physics communities. We welcome suggestions concerning mathematical corrections, omitted citations, and attribution of ideas.

P.L. was partially supported by NSF grant DMS-2450004.

Introduction

The Sherrington–Kirkpatrick model and the Parisi formula

Let (gij)1≤i<j≤N(g_{ij})_{1\le i<j\le N} be independent standard Gaussian variables. The Sherrington–Kirkpatrick (SK) model [76] at zero external field and inverse temperature β>0\beta> 0 has Hamiltonian

HN(σ)=1N∑1≤i<j≤Ngijσiσj,σ∈{−1,1}N.H_N(\sigma)=\frac{1}{\sqrt{N}}\sum_{1\le i<j\le N}g_{ij}\sigma_i\sigma_j,\qquad\sigma\in\{-1,1\}^N.

Its logarithmic partition function and expected pressure are

FN,β=log⁡∑σ∈{−1,1}NeβHN(σ),pN,β=EFN,βN.F_{N,\beta}=\log\sum_{\sigma\in\{-1,1\}^N}e^{\beta H_N(\sigma)},\qquad p_{N,\beta}=\frac{\mathbb{E}F_{N,\beta}}{N}.

For a probability measure μ\mu on [0,1][0,1] with distribution function α(q)=μ([0,q])\alpha(q)=\mu([0,q]), let uαu_\alpha be the solution of the Parisi equation on [0,1]×R[0,1]\times\mathbb{R}, and define the Parisi functional Pβ\mathcal{P}_\beta and the Parisi value PβP_\beta by

∂quα=−β22(∂xxuα+α(q)(∂xuα)2),uα(1,x)=log⁡(2cosh⁡x),Pβ(α)=uα(0,0)−β22∫01qα(q) dq,Pβ=inf⁡μ∈P([0,1])Pβ(α),(1)\begin{aligned} \partial_q u_\alpha&=-\frac{\beta^2}{2}\left(\partial_{xx}u_\alpha+\alpha(q)(\partial_xu_\alpha)^2\right),\qquad u_\alpha(1,x)=\log(2\cosh x),\\ \mathcal{P}_\beta(\alpha)&=u_\alpha(0,0)-\frac{\beta^2}{2}\int_0^1q\alpha(q)\,dq,\qquad P_\beta=\inf_{\mu\in\mathcal{P}([0,1])}\mathcal{P}_\beta(\alpha), \tag*{(1)} \end{aligned}

where P([0,1])\mathcal{P}([0,1]) is the set of probability measures on [0,1][0,1]. For a general μ\mu, the equation is understood in the weak sense of [43, Theorem 4].

Sherrington and Kirkpatrick introduced this model and studied it by a replica-symmetric calculation, which fails at low temperature [76]; a longer treatment, with numerical comparisons, is given in [45]. De Almeida and Thouless showed that the replica-symmetric saddle point becomes unstable [34]. Parisi proposed a hierarchical breaking of the replica symmetry in [62], and he developed a sequence of approximate solutions and the order parameter on [0,1][0,1] that leads to (1.1) in [63, 64]. He later interpreted this order parameter through the distribution of overlaps between states [65]. Mézard, Parisi, Sourlas, Toulouse, and Virasoro then described the organization of the states, with sample-dependent weights and ultrametric overlaps, in [50].

The Parisi formula states that pN,β→Pβp_{N,\beta}\to P_\beta as N→∞N\to\infty. Guerra proved the upper bound pN,β≤Pβp_{N,\beta}\le P_\beta for every NN [40], and Talagrand proved the matching lower bound in the limit [78]. The infimum in (1.1) is attained at a unique measure μβ\mu_\beta [10], the Parisi measure, whose distribution function we denote by αβ\alpha_\beta. For β>1\beta> 1, its support is an interval [0,qβ][0,q_\beta] with 0<qβ<10 < q_\beta< 1 [48, Theorem 1.1]; for β>1\beta> 1 close to 1, this was proved earlier by Zhou [86].

Finite-size corrections

In this paper, we study the rate of convergence in the Parisi formula at low temperature. By Guerra’s bound, the mean bias Pβ−pN,βP_\beta-p_{N,\beta} is nonnegative for every NN. Up to a factor 1/β1/\beta, it is the finite-size correction to the expected free energy per spin (§ 2.3).

At high temperature, the mean bias is small. For 0<β<10 < \beta< 1, Guerra and Toninelli showed by a quadratic coupling of two replicas that it is of order at most N−1N^{-1} [42, Section 3]. With a nonzero external field, Lopatto proved that the limiting pressure equals the replica-symmetric value up to and including the de Almeida–Thouless line [49, Theorem 1.1], and Kusuoka and Nakajima proved that the pressure is within O(N−1)O(N^{-1}) of this value on compact subsets of the strict de Almeida–Thouless region [47, Theorem 1.1].

Below the critical temperature, the mean bias is expected to decay more slowly. In [4 Theorem 1.1, (1.5)], we proved that 0≤Pβ−pN,β≤CβN−1/20 \le P_\beta-p_{N,\beta} \le C_\beta N^{-1/2} for every β>0\beta> 0 and N≥1N \ge1, with a numerical constant CC. For fixed β>1\beta> 1, we used the interval support of μβ\mu_\beta to prove the sharper upper bound CβN−2/3C_\beta N^{-2/3} and the lower bound cβN−11/12c_\beta N^{-11/12}, for all sufficiently large NN [4 Theorem 1.4]. In physics, the correction is expected to decay like N−ωN^{-\omega} with ω<1\omega< 1 at low temperature. Two values have been proposed for the full model, ω=2/3\omega= 2/3 and ω=3/4\omega= 3/4. A calculation retaining only longitudinal fluctuations gives ω=4/5\omega= 4/5 [72]. We review these predictions in § 2.4. We improve the exponent of the upper bound from 2/32/3 to 7/107/10.

Main result

Our main result is the following.

Theorem 1.1. For every β>1\beta> 1 there is a constant Cβ<∞C_\beta< \infty such that

Pβ−pN,β≤CβN−7/10(N≥2).P_\beta-p_{N,\beta} \le C_\beta N^{-7/10} \qquad(N \ge2).

Since Pβ−pN,β≥0P_\beta-p_{N,\beta} \ge0 for every NN [40], Theorem 1.1 is a two-sided estimate. In physical units, the expected free energy density is E[−FN,β/(βN)]\mathbb{E}[-F_{N,\beta}/(\beta N)], and its limit is −Pβ/β-P_\beta/\beta. By Theorem 1.1, we have

0≤E[−FN,ββN]−(−Pββ)=Pβ−pN,ββ≤CββN−7/10(N≥2).0 \le\mathbb{E}\left[-\frac{F_{N,\beta}}{\beta N}\right]-\left(-\frac{P_\beta}{\beta}\right)=\frac{P_\beta-p_{N,\beta}}{\beta}\le\frac{C_\beta}{\beta}N^{-7/10}\qquad(N \ge2).

For the expected logarithm of the partition function itself, the corresponding bound is

0≤NPβ−EFN,β≤CβN3/10.0 \le N P_{\beta} - \mathbb{E}F_{N,\beta} \le C_{\beta}N^{3/10}.

The constant CβC_{\beta} depends on β\beta, and we do not control it as β↓1\beta\downarrow1 or as β→∞\beta\to\infty. In particular, we do not obtain a corresponding bound for the ground-state energy.

Together with the lower bound in [4 Theorem 1.4], Theorem 1.1 implies that for every fixed β>1\beta> 1,

cβN−11/12≤Pβ−pN,β≤CβN−7/10.c_{\beta}N^{-11/12} \le P_{\beta} - p_{N,\beta} \le C_{\beta}N^{-7/10}.

for all sufficiently large NN. The exponent 7/107/10 lies between the exponent 2/32/3 of [4 Theorem 1.4] and the exponent 3/43/4 predicted by Franz, Parisi, and Ricci-Tersenghi [38], which we discuss in § 2.4. In § 2.6, we compare Theorem 1.1 with two conditional consequences of results stated in the OpenAI manuscripts [54, 56]. An uncentered comparison stated there would give an unspecified strict improvement over 2/32/3, and the variance asymptotic stated there would imply the lower bound Pβ−pN,β≥cβN−5/6P_{\beta} - p_{N,\beta} \ge c_{\beta}N^{-5/6}.

We now outline the proof; Section 2 gives a detailed overview. We follow the interpolation scheme of [4], which is a quantitative form of the Guerra interpolation [40] that Talagrand used for the positive fractional moments of the partition function [80]. For a tilt 0<λ≤λ00 < \lambda\le\lambda_{0}, let eN,β(λ)≥0e_{N,\beta}(\lambda) \ge0 be the error, in energy units, between the constrained Parisi value and the normalized logarithmic moment of order λ/β\lambda/\beta of the completed partition function. The constraint requires mass at least λ/β\lambda/\beta at zero (§ 4.1). Imposing this mass raises the Parisi value by at most Cβλ5C_{\beta}\lambda^{5} in energy units [4 Lemma 7.2]. Using this bound and the convexity of Gaussian logarithmic moments [26, Theorem 1], as in the proof of [4 Theorem 7.8, Step 1], we find

0≤Pβ−pN,ββ≤Cβλ5+2eN,β(λ).0 \le\frac{P_{\beta} - p_{N,\beta}}{\beta} \le C_{\beta}\lambda^{5} + 2e_{N,\beta}(\lambda).

For sufficiently large NN, we take λ=N−7/50\lambda= N^{-7/50}, so it suffices to show that eN,β(N−7/50)≤CβN−7/10e_{N,\beta}(N^{-7/50}) \le C_{\beta}N^{-7/10}.

References

  1. [1]Michael Aizenman, Joel Louis Lebowitz, and David Ruelle. Some rigorous results on the Sherrington–Kirkpatrick spin glass model. Communications in Mathematical Physics, 112(1):3–20, 1987. doi: 10.1007/BF01217677.DOI
  2. [2]Michael Aizenman, Robert Sims, and Shannon Lee Starr. Extended variational principle for the Sherrington–Kirkpatrick spin-glass model. Physical Review B, 68(21):214403, 2003. doi: 10.1103/PhysRevB.68.214403.DOI
  3. [3]Theodore Wilbur Anderson. The integral of a symmetric unimodal function over a symmetric convex set and some probability inequalities. Proceedings of the American Mathematical Society, 6(2):170–176, 1955. doi: 10.1090/S0002-9939-1955-0069229-1.DOI
  4. [4]P. M. Aronow and Patrick Lopatto. Quantitative Parisi formulas and fluctuations in the Sherrington–Kirkpatrick model, 2026. URL https://arxiv.org/abs/2609.09103v1. Statement and equation locators refer to arXiv:2609.09103v1.
  5. [5]P. M. Aronow and Patrick Lopatto. The Parisi formula for the Edwards–Anderson model in the high-dimensional limit. Hexagon, 2026. hexagon:2610.00155v1.
  6. [6]Timo Aspelmeier. Bond chaos in the Sherrington–Kirkpatrick model. Journal of Physics A: Mathematical and Theoretical, 41(20):205005, 2008. doi: 10.1088/1751-8113/41/20/205005. Section locators refer to arXiv:0802.3204v1.DOI
  7. [7]Timo Aspelmeier. An exact relation between free energy fluctuations and bond chaos in the Sherrington–Kirkpatrick model. Journal of Statistical Mechanics: Theory and Experiment, 2008(4):P04018, 2008. doi: 10.1088/1742-5468/2008/04/P04018. Equation locators refer to arXiv:0802.3007v1.DOI
  8. [8]Timo Aspelmeier, Alain Billoire, Enzo Marinari, and Michael A. Moore. Finite-size corrections in the Sherrington–Kirkpatrick model. Journal of Physics A: Mathematical and Theoretical, 41(32):324008, 2008. doi: 10.1088/1751-8113/41/32/324008.DOI
  9. [9]Antonio Auffinger and Wei-Kuo Chen. On properties of Parisi measures. Probability Theory and Related Fields, 161(3–4):817–850, 2015. doi: 10.1007/s00440-014-0563-y.DOI
  10. [10]Antonio Auffinger and Wei-Kuo Chen. The Parisi formula has a unique minimizer. Communications in Mathematical Physics, 335(3):1429–1444, 2015. doi: 10.1007/s00220-014-2254-z. Statement locators refer to arXiv:1402.5132v3.DOI
  11. [11]Antonio Auffinger, Wei-Kuo Chen, and Qiang Zheng. The SK model is infinite step replica symmetry breaking at zero temperature. Communications on Pure and Applied Mathematics, 73(5):921–943, 2020. doi: 10.1002/cpa.21886.
  12. [12]JinHo Baik and Ji Oon Lee. Fluctuations of the free energy of the spherical Sherrington–Kirkpatrick model. Journal of Statistical Physics, 165(2):185–224, 2016. doi: 10.1007/s10955-016-1610-0. Statement locators refer to arXiv:1505.07349v3.DOI
  13. [13]Sergei Natanovich Bernstein. Démonstration du théorème de Weierstrass fondée sur le calcul des probabilités. Communications de la Société mathématique de Kharkov. Deuxième série, 13(1):1–2, 1912. URL https://www.mathnet.ru/eng/khmo107.
  14. [14]Alain Bill oire. Numerical estimate of the finite-size corrections to the free energy of the Sherrington–Kirkpatrick model using Guerra–Toninelli interpolation. Physical Review B, 73(13):132201, 2006. doi: 10.1103/PhysRevB.73.132201.
  15. [15]Stefan Boettcher. Extremal optimization for Sherrington–Kirkpatrick spin glasses. The European Physical Journal B, 46(4):501–505, 2005. doi: 10.1140/epjb/e2005-00280-6. Section locators refer to arXiv:cond-mat/0407130v2.DOI
  16. [16]Stefan Boettcher. Simulations of ground state fluctuations in mean-field Ising spin glasses. Journal of Statistical Mechanics: Theory and Experiment, 2010(7):P07002, 2010. doi: 10.1088/1742-5468/2010/07/P07002. Section locators refer to arXiv:0906.1292v2.DOI
  17. [17]Stefan Boettcher. Ground state properties of the diluted Sherrington–Kirkpatrick spin glass. Physical Review Letters, 124(17):177202, 2020. doi: 10.1103/PhysRevLett.124.177202. Table locators refer to arXiv:1912.04974v3.DOI
  18. [19]Jean-Philippe Bouchaud, Florent Krzakala, and Olivier C. Martin. Energy exponents and corrections to scaling in Ising spin glasses. Physical Review B, 68(22):224404, 2003. doi: 10.1103/PhysRevB.68.224404. Section locators refer to arXiv:cond-mat/0212070v2.DOI
  19. [20]Anton Bovier and Anton Klimovsky. The Aizenman–Sims–Starr and Guerra’s schemes for the SK model with multidimensional spins. Electronic Journal of Probability, 14(8):161–241, 2009. doi: 10.1214/EJP.v14-611. Section locators refer to arXiv:0802.3467v2.DOI
  20. [21]Herm Jan Brascamp and Elliott Hershel Lieb. On extensions of the Brunn–Minkowski and Prékopa–Leindler theorems, including inequalities for log concave functions, and with an application to the diffusion equation. Journal of Functional Analysis, 22 (4):366–389, 1976. doi: 10.1016/0022-1236(76)90004-5.DOI
  21. [22]Sourav Chatterjee. Disorder chaos and multiple valleys in spin glasses, 2009. Statement locators refer to arXiv:0907.3381v4.arxiv.org/abs/0907.3381
  22. [23]Hong-Bin Chen. FRSB in the SK spin glass: convergence to full-interval support at zero temperature, 2026. Statement locators refer to arXiv:2607.18032v1.arxiv.org/abs/2607.18032

Paper details

vlnc:W_g8cxf4vw
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

OpenAI ChatGPT 5.6, OpenAI ChatGPT 6, OpenAI ChatGPT 6.1, Anthropic Claude Opus 5.5, Anthropic Claude Fable 5, Anthropic Claude Fable 5.1 · View statement →

Looking for related papers…
Contents