alphaFeedback

WorksMathematics

Slow Equilibration and Persistent Stationary Correlations in Spherical Spin Glasses

AI contributions
  • Anthropic Claude Opus 5.5
  • OpenAI ChatGPT 5.6
  • OpenAI ChatGPT 6
  • OpenAI ChatGPT 6.1
View details

Abstract

We study the Langevin dynamics of spherical mixed pp-spin glasses at fixed temperature, on time scales exponential in the dimension NN. Below the static transition temperature, the dynamics started from the uniform distribution on the sphere stays at total variation distance close to one from the Gibbs measure up to a time exponential in NN, with probability tending to one. We prove this for pure models of degree at least four at every such temperature, and for mixtures of minimum degree at least four that satisfy an explicit condition on a function of one variable. It extends to initial distributions that are independent of the disorder or depend weakly on it. Between the dynamical and static transition temperatures, we show for general finite mixtures that, with probability exponentially close to one, the generator has exponentially many exponentially small eigenvalues; a fixed part of the centered spin lies in the corresponding slow modes; the stationary spin correlation stays above every level below q∗2q_*^2 up to an exponential time, where q∗q_* is an explicit root of the plateau equation of the physics literature; the dynamics started from most points of the Gibbs measure stays far from equilibrium; and the two endpoints of the stationary dynamics share mutual information at least linear in NN. With probability tending to one, the correlation stays above every level below q∗q_*. Near the dynamical transition, exponentially many eigenvalues stay below e−cNe^{-cN}, where cc is proportional to an explicit power of the distance to the transition.

Authors’ Note

With the exception of this note, the entire document is machine-written. We have carefully supervised its preparation to ensure that the proofs are readable and that previous literature is appropriately cited. We are releasing this manuscript on Hexagon to disseminate the result as quickly as possible because it may be of interest to the mathematical community. We welcome suggestions concerning mathematical corrections, omitted citations, and attribution of ideas.

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

Introduction

Spherical spin glasses are mean-field models of disordered systems, in which the energy is a random Gaussian function on a high-dimensional sphere. Following the work of Kirkpatrick, Thirumalai, and Wolynes [100–102], pp-spin models have served as mean-field models of the glass transition, and spherical models are among those that can be solved: Crisanti and Sommers computed their statics [52], and Crisanti, Horner, and Sommers their dynamics [54]. For pure spherical pp-spin models with p≥3p \ge3, the standard picture involves a dynamical transition temperature TdT_d and a lower static transition temperature TsT_s. Between TsT_s and TdT_d, the Langevin dynamics is predicted to fall out of equilibrium, although the free energy keeps its annealed value: equilibrium dynamics is predicted to remain for very long times in one of exponentially many states that together carry the Gibbs measure [19, 54, 73]. Below TsT_s, the Gibbs measure itself condenses onto few states. Started instead from a random configuration, the dynamics below TdT_d is predicted to age and to stay away from the states that carry the Gibbs measure, on all time scales that survive the limit of infinite dimension [56]. Mixtures can behave differently: the two transitions may coincide (Example 18.10), and the long-time behavior can depend on the preparation [70].

In this article, we prove versions of these predictions in finite dimension NN, on time scales exponential in NN. We study the Langevin dynamics of spherical mixed pp-spin glasses at a fixed temperature and ask two questions. The first concerns temperatures below TsT_s: does the dynamics started from the uniform law on the sphere, which does not depend on the disorder, reach equilibrium, and on what time scale? The second concerns the interval Ts<T<TdT_s < T < T_d: what does the equilibrium dynamics remember of its starting point? For the first question, we show that the law of the dynamics started from the uniform law stays at total variation distance tending to one from the Gibbs measure up to a time exponential in NN. This holds at every temperature below TsT_s for pure models of degree at least four, and for the mixtures in a class defined by an explicit scalar condition (Theorem 1.2 and Corollary 1.3), and it extends to every initial law chosen independently of the disorder (Theorem 1.4). For the second question, we show for general finite mixtures that the generator has exponentially many exponentially small eigenvalues, that a fixed part of the spin lies in the corresponding slow modes, that the stationary spin correlation stays above a plateau up to an exponential time, that the dynamics started from most points of the Gibbs measure stays far from equilibrium, and that the two endpoints of the stationary dynamics share mutual information at least linear in NN (Theorems 1.5 to 1.7).

The model and the two thresholds

Let SN={x∈RN:∣x∣2=N}S_N = \{x \in\mathbb{R}^N : \lvert x\rvert^2 = N\}, and let νN\nu_N be the uniform probability measure on SNS_N. For x,y∈RNx,y \in\mathbb{R}^N, we write R(x,y)=x⋅y/NR(x,y) = x \cdot y/N for the overlap and U(x)=x/NU(x) = x/\sqrt{N} for the normalized spin. A finite mixture is a nonzero polynomial

ξ(r)=∑p=2Pcprp,cp≥0.\xi(r) = \sum_{p=2}^{P} c_p r^p,\qquad c_p \ge0.

Its minimum degree dd is the smallest pp with cp>0c_p > 0, and ξ\xi is pure if exactly one coefficient is nonzero. The field HNH_N is the centered Gaussian process on SNS_N with covariance

EHN(x)HN(y)=Nξ(R(x,y)).\mathbb{E}H_N(x)H_N(y) = N\xi(R(x,y)).

For ξ(r)=rp\xi(r) = r^p, we may take HN(x)=N−(p−1)/2∑i1,…,ipgi1⋯ipxi1⋯xipH_N(x) = N^{-(p-1)/2} \sum_{i_1,\ldots,i_p} g_{i_1\cdots i_p}x_{i_1}\cdots x_{i_p} with independent standard Gaussian coefficients, and a mixture is a sum of independent fields of this form with weights cp\sqrt{c_p}. We allow a quadratic term and odd degrees. We refer to the randomness of HNH_N as the disorder. For an inverse temperature β≥0\beta\ge0, the Gibbs measure and the free energy are

π(dx)=eβHN(x)ZN,βνN(dx),ZN,β=∫SNeβHN dνN,F(β)=lim⁡N→∞1NElog⁡ZN,β.\pi(dx) = \frac{e^{\beta H_N(x)}}{Z_{N,\beta}}\nu_N(dx),\qquad Z_{N,\beta} = \int_{S_N}e^{\beta H_N}\,d\nu_N,\qquad F(\beta) = \lim_{N\to\infty}\frac{1}{N}\mathbb{E}\log Z_{N,\beta}.

The Gibbs measure favors large values of HNH_N; in the physics literature, the energy is usually −HN-H_N. The limit F(β)F(\beta) exists by the Parisi formula, which we recall below. The Langevin dynamics is the diffusion on SNS_N with generator

L=ΔSN+β∇SNHN⋅∇SN,L = \Delta_{S_N} + \beta\nabla_{S_N}H_N \cdot\nabla_{S_N},

that is, the solution of dXt=2PXt dWt+β∇SNHN(Xt) dt−(1−1/N)Xt dtdX_t = \sqrt{2}P_{X_t}\,dW_t + \beta\nabla_{S_N}H_N(X_t)\,dt - (1 - 1/N)X_t\,dt, where WW is a standard Brownian motion in RN\mathbb{R}^N and PxP_x is the orthogonal projection onto the tangent space of SNS_N at xx. It is reversible with respect to π\pi, and its Dirichlet form is E(f,f)=∫∣∇SNf∣2 dπ\mathcal{E}(f,f) = \int\lvert\nabla_{S_N}f\rvert^2\,d\pi. The operator −L-L has discrete spectrum 0=λ0<λ1≤λ2≤⋯0 = \lambda_0 < \lambda_1 \le\lambda_2 \le\cdots, where eigenvalues are counted with multiplicity. We write Pt=etLP_t = e^{tL}, Pt(x,⋅)P_t(x,\cdot) for its transition kernel, and μPt\mu P_t for the law at time tt of the diffusion with initial law μ\mu. All these objects depend on the realization of HNH_N, which we suppress from the notation. The total variation distance is ∥μ−ν∥TV=sup⁡A∣μ(A)−ν(A)∣\lVert\mu- \nu\rVert_{\mathrm{TV}} = \sup_A\lvert\mu(A) - \nu(A)\rvert. Physicists often write the generator as TΔ−∇E⋅∇T\Delta- \nabla E \cdot\nabla with E=−HNE = -H_N; this changes the time scale by a constant factor, which does not matter on the exponential time scales considered here.

The two thresholds are defined through the functions

A(r)=(1−r)ξ′(r)r,Φβ(r)=β2ξ(r)+r+log⁡(1−r),0<r<1.A(r)=\frac{(1-r)\xi'(r)}{r},\qquad\Phi_\beta(r)=\beta^2\xi(r)+r+\log(1-r),\qquad0<r<1.

The dynamical threshold βd=βd(ξ)\beta_{\mathrm d}=\beta_{\mathrm d}(\xi) and the static threshold βs=βs(ξ)\beta_{\mathrm s}=\beta_{\mathrm s}(\xi) are given by

βd−2=sup⁡0<r<1A(r),βs2=inf⁡0<r<1−r−log⁡(1−r)ξ(r).(1)\beta_{\mathrm d}^{-2}=\sup_{0<r<1} A(r),\qquad\beta_{\mathrm s}^{2}=\inf_{0<r<1}\frac{-r-\log(1-r)}{\xi(r)}. \tag*{(1)}

For a pure model, we write βd(p)\beta_{\mathrm d}(p) and βs(p)\beta_{\mathrm s}(p) for the thresholds of ξ(r)=rp\xi(r)=r^p.

The static threshold is the point where the free energy leaves its annealed value. By Jensen’s inequality, we have F(β)≤lim⁡NN−1log⁡EZN,β=β2ξ(1)/2F(\beta)\leq\lim_{N}N^{-1}\log\mathbb{E}Z_{N,\beta}=\beta^2\xi(1)/2, and equality holds if and only if Φβ≤0\Phi_\beta\leq0 on (0,1)(0,1), that is, if and only if β≤βs\beta\leq\beta_{\mathrm s}. For even mixtures, Talagrand proved that the corresponding strict inequality is sufficient for this equality and that the non-strict inequality is necessary [146, Proposition 2.3]. The equivalence follows from the optimality criterion [146, Proposition 2.1] and the Parisi formula, which Chen proved for general mixtures [46] (see Lemma 3.18). Following Huang and Sellke [90, Definition 3], we say that (ξ,β)(\xi,\beta) is strictly replica symmetric if Φβ<0\Phi_\beta<0 on (0,1)(0,1) (Definition 3.16). For every β<βs\beta<\beta_{\mathrm s}, the pair is strictly replica symmetric and β2ξ′′(0)<1\beta^2\xi''(0)<1 (Lemma 3.18). We call the range β>βs\beta>\beta_{\mathrm s} the condensed phase.

The dynamical threshold enters through the derivative

Φβ′(r)=r1−r(β2A(r)−1).\Phi_\beta'(r)=\frac{r}{1-r}\left(\beta^2A(r)-1\right).

If β≤βd\beta\leq\beta_{\mathrm d}, then this identity implies that Φβ\Phi_\beta is nonincreasing, hence Φβ≤0\Phi_\beta\leq0. We conclude that βd≤βs\beta_{\mathrm d}\leq\beta_{\mathrm s}. For β>βd\beta>\beta_{\mathrm d}, the set {r∈(0,1):β2A(r)>1}\{r\in(0,1):\beta^2A(r)>1\} is nonempty and its closure does not contain 11, since A(1−)=0A(1^-)=0. We define its right endpoint

q∗=sup⁡{r∈(0,1):β2(1−r)ξ′(r)>r}.(2)q_*=\sup\{r\in(0,1):\beta^2(1-r)\xi'(r)>r\}. \tag*{(2)}

References

  1. [93]Huang, Montanari, Pham. [Bibliographic entry not printed on this page]

Paper details

vlnc:W_epvnaycc
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

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

Looking for related papers…
Contents