Sherrington-Kirkpatrick Spin Glass

• Sherrington-Kirkpatrick spin glass is a foundational mean-field model that captures disorder and frustration among Ising spins with infinite-range Gaussian interactions.
• The model utilizes replica symmetry breaking to describe a complex free energy landscape and phase transitions from paramagnetic to spin glass states.
• Quantum extensions incorporate transverse fields to explore phenomena like many-body localization and discrete time crystals in glassy systems.

The Sherrington-Kirkpatrick (SK) spin glass is the foundational mean-field model for spin glass behavior, encoding the interplay of disorder and frustration in a system of $N$ Ising spins with infinite-range, quenched random interactions. Introduced by Sherrington and Kirkpatrick in 1975 as a fully connected analogue of the Edwards-Anderson model, its rigorous mathematical structure and rich phenomenology have made it a canonical system in statistical mechanics, probability, condensed matter, and quantum information science (Sherrington et al., 30 May 2025). Central to the SK model are concepts such as replica symmetry breaking (RSB), complex free energy landscapes, glassy order, and, in quantum extensions, the effects of transverse fields, @@@@1@@@@, and non-trivial dynamical responses.

1. Formulation and Mean-Field Structure

The classical SK Hamiltonian for $N$ Ising spins $\sigma_i = \pm 1$ is

$$H[\sigma;J] = -\sum_{1\le i<j \le N} J_{ij} \sigma_i \sigma_j$$

where $J_{ij}$ are i.i.d. Gaussian variables with zero mean and variance $J^2/N$. This $1/N$ scaling ensures a non-trivial thermodynamic limit (Sherrington et al., 30 May 2025). The disorder-averaged thermodynamics are obtained from the quenched free energy: $$F = -\frac{1}{\beta} \; \mathbb{E}_J \ln Z[J], \qquad Z[J] = \sum_{\sigma} e^{-\beta H[\sigma;J]}$$

In the presence of a uniform transverse field $\Gamma$, the quantum SK model for $N$ spin-1/2's has Hamiltonian

$$\hat{\cal H} = -\sum_{1\le i<j\le N} J_{ij}\,\hat{\sigma}_i^z\,\hat{\sigma}_j^z - \Gamma\sum_{i=1}^N \hat{\sigma}_i^x$$

with $J_{ij}$ as above. Longitudinal fields or correlated random fields can also be included (Kiss et al., 2023, Hadjiagapiou, 2014).

2. Replica Method, Order Parameters, and RSB

The analytic solution uses the replica trick: $$\mathbb{E}\ln Z = \lim_{n\to 0} ( \mathbb{E}[Z^n] - 1 ) / n$$ (Sherrington et al., 30 May 2025). Disorder averaging introduces replica overlaps

$$q_{\alpha\beta} = \frac{1}{N} \sum_{i=1}^N \sigma_i^\alpha \sigma_i^\beta$$

and the free energy is given by a saddle point over these overlaps. Assuming replica symmetry (RS), all $q_{\alpha\beta} = q$, leading to the mean-field self-consistency: $$q = \int Dz \, \tanh^2(\beta J \sqrt{q}\, z), \qquad Dz = \frac{dz}{\sqrt{2\pi}} e^{-z^2/2}$$

However, below a critical temperature $T_c = J$ (Rodríguez-Camargo et al., 2021), the RS solution becomes unstable. Parisi's scheme introduces continuous replica symmetry breaking (full RSB): the overlap order parameter becomes a function $q(x)$, with $x \in [0,1]$ representing a hierarchy of ergodicity breaking. The Parisi functional is extremized with respect to $q(x)$, capturing the ultrametric structure of pure states and the nontrivial distribution of overlaps $P(q)$ (Sherrington et al., 30 May 2025, Talebi, 12 Jun 2025, Kiss et al., 2023).

3. Phase Diagram, Stability, and the de Almeida–Thouless Line

The SK model exhibits a high-temperature paramagnetic phase and a low-temperature spin glass phase. The transition is marked by the vanishing of the solution $q=0$, while the stability of the RS solution is controlled by the replicon eigenvalue: $\lambda_{\text{AT}} = 1 - \beta^2 J^2 \int Dz \, \sech^4(\beta J \sqrt{q} \, z)$ The line $\lambda_{\text{AT}}=0$ in $T$–$h$ space is the de Almeida–Thouless (AT) line, below which full RSB is required (Sherrington et al., 30 May 2025, Talebi, 12 Jun 2025, Hadjiagapiou, 2014). In quantum extensions, the phase boundary generalizes to a surface in $(T,\Gamma)$ (transverse field), with the quantum AT line (QuAT) marking the RSB onset at $T=0$ (Kiss et al., 2023, Young, 2017).

A table summarizing key phase boundaries in the classical/quantum SK model:

Parameter regime	Transition/instability	Criterion	Ref
$T=T_c=J$ (classical)	PM–SG (RS solution appears)	$q=0$ loses stability	(Sherrington et al., 30 May 2025)
$T<T_{AT}(h)$	SG: RS $\to$ full RSB	$\lambda_{\rm AT}=0$	(Talebi, 12 Jun 2025)
$\Gamma>\Gamma_c$	QPM–QSG (Quantum critical)	$\Gamma_c \approx 1.5J$	(Kiss et al., 2023)
$h>0$, $\Gamma>0$	No AT line (Quantum)	Only RS/ergodic glass	(Rajak et al., 2023)

4. Quantum SK Model: CTQMC, Parisi RSB, and Phase Diagram

For the quantum SK model in transverse field, the $N\to\infty$ solution is mapped to a self-consistent single-site effective theory, solvable via continuous-time quantum Monte Carlo (CTQMC) (Kiss et al., 2023). The quantum action incorporates both static and dynamic order parameters: $$S_{\rm eff}[y] = \int_0^\beta d\tau \big(y \sigma^z_\tau + \Gamma \sigma^x_\tau \big) - \frac{J^2}{2}\int_0^\beta d\tau d\tau' \, \tilde{\chi}(\tau-\tau') \, \sigma^z_\tau \sigma^z_{\tau'}$$ Replica symmetry breaking is implemented via continuous $Q(x)$ and corresponding Parisi flow equations for the field distribution $P(x,y)$ and scale-dependent free energy $\phi(x,y)$. The phase diagram features a quantum glass at small $\Gamma$, bounded by a critical line $\Gamma_c(T)$ with a continuous transition to a quantum paramagnet. In the quantum limit ($T=0$), $\Gamma_c \simeq 1.5J$ (Kiss et al., 2023, Young, 2017).

Inclusion of a longitudinal field $h > 0$ leads to restoration of replica symmetry and the disappearance of the AT boundary, with any nonzero $h$ rendering the system ergodic in the quantum glass phase (Rajak et al., 2023).

5. Free Energy Landscape and Dynamical Features

The structure of free energy minima and the landscape topology are central to glassy dynamics. The TAP (Thouless-Anderson-Palmer) equations define local magnetizations $m_i$, and every TAP minimum is paired with a nearby index-one saddle, with barriers $\Delta F$ distinguishing distinct classes of states:

• For $f > f_c$ (above the RSB threshold), barriers $\Delta F \sim N^{-2}$ vanish in large $N$, so only marginally stable states are relevant in large systems.
• For $f < f_c$, barriers grow as $\sim N^{1/3}$, sustaining the existence of pure states and glassy order in the thermodynamic limit (Aspelmeier et al., 2021).

This landscape underpins non-ergodic aging, slow relaxation, and the proliferation of metastable states.

6. Quantum Dynamics: Many-Body Localization and Discrete Time Crystals

In the quantum SK model, many-body localization (MBL) phenomena are realized despite infinite-range interactions. Numerical diagnostics (participation ratio, level statistics, Renyi entropy) reveal mobility edges in energy-density–$\Gamma$ space, separating non-ergodic, area-law-entangled glassy phases from ergodic, thermal regions (Mukherjee et al., 2017). Quantum spin glass order coincides with the MBL regime, with both destroyed above the quantum critical field $\Gamma_{CP} \simeq 1.5J$.

Periodically driven (Floquet) quantum SK models exhibit robust discrete time crystal (DTC) phases, even with long-range random interactions. The DTC order parameter follows the non-ergodic regime (as measured via the Shannon entropy), reinforcing the central organizing role of glassy non-ergodicity for exotic dynamical phases (Bothra et al., 27 Apr 2025).

7. Applications, Extensions, and Analytical Innovations

The SK model provides a framework for understanding quantum annealing performance, the effect of correlated disorder (e.g., joint Gaussian random fields with nonzero correlation), generalized models (diluted SK, multi-species, mixed $p$-spin), and spectral singularities near condensation (Das et al., 2024, Hadjiagapiou, 2014, Bates et al., 2018, Boettcher, 2019, Wang et al., 2024). New analytic approaches (e.g., distributional zeta-function method, rigorous path-integral-based variational characterizations) supplement the classical replica and Parisi frameworks, enabling exact results for classical and quantum phase diagrams, susceptibilities, finite-size corrections, and overlap distributions (Rodríguez-Camargo et al., 2021, Adhikari et al., 2019).

In the setting of quantum annealing, quantum tunneling introduces ergodicity: a small longitudinal field in the quantum SK model leads to instantaneous restoration of RS and elimination of the AT transition, yielding enhanced ground state preparation and fundamentally distinct quantum dynamics compared to the classical limit (Rajak et al., 2023).

---

The SK spin glass thus remains a central paradigm, exemplifying the interplay of disorder, frustration, and quantum dynamics in mean-field theory, and providing a blueprint for analytical, numerical, and experimental investigations into complex glassy systems (Sherrington et al., 30 May 2025, Kiss et al., 2023, Young, 2017).