Finite Edwards-Anderson Order Parameter

• Finite Edwards-Anderson order parameter is defined through spatial and disorder averaging of squared local observables, capturing persistent frozen order in disordered phases.
• It distinguishes glassy states from homogeneous phases by measuring overlaps and local density fluctuations via experimental and simulation techniques.
• Finite-size corrections and scaling analyses are crucial, using methods like Monte Carlo sampling and field-theoretic approaches to assess the robustness of glassy order.

The finite Edwards-Anderson (EA) order parameter quantifies the extent of frozen, random order in disordered systems, particularly spin glasses and related models. It serves as a central diagnostic to distinguish true glassy phases from homogeneous or thermally fluctuating phases, both in classical and quantum contexts. Finite values of the EA order parameter reflect persistent local order that survives both disorder and ensemble averaging, marking regimes of replica symmetry breaking or glass rigidity. Techniques and definitions vary by physical system, but all constructions involve spatial averaging and disorder ensembles, directly linking q_{EA} to physical observables such as overlaps, local densities, or spectral linewidths.

1. Formal Definitions of the EA Order Parameter

The EA order parameter is generally defined in terms of the square of local observables, averaged over both thermal fluctuations and disorder. For an Ising spin glass of N spins with variables $S_i = \pm1$, the conventional finite-size EA parameter is

$$q_{\rm EA}(L,T) = \frac{1}{N}\sum_{i=1}^{N}\bigl[\langle S_i\rangle^2\bigr]_{\rm av}$$

but, operationally, it is often extracted as the position of the dominant peak of the overlap distribution $P(q)$, where the overlap between two replicas $\alpha$ and $\beta$ of the same disorder realization is

$q = \frac{1}{N}\sum_{i=1}^N S_i^\alpha S_i^\beta$

and

$q_{\rm EA}(L,T) = \arg\max_{q>0} P(q)$

In quantum or bosonic contexts, such as the disordered Bose-Hubbard model, the EA order parameter for the Bose glass is

$$q_{\rm EA} = \overline{\,\langle \hat n_i\rangle^2\,} - \overline{\langle \hat n_i\rangle}^2 \equiv \overline{\bigl(\langle \hat n_i\rangle - \overline{\langle \hat n_i\rangle}\bigr)^2}$$

where the overbar denotes a disorder average, and angular brackets denote a thermal or mean-field average (Thomson et al., 2016). In random-field Ising-nematic models relevant to FeSe, the EA parameter is

$$q_{\rm EA} = \frac{1}{N} \sum_{i} [ \langle \tau_i^z \rangle^2 ]_{\text{dis}} - \bar\phi^2$$

with $\tau_i^z$ a local Ising-nematic pseudospin, and $$q_{\rm EA}(L,T) = \frac{1}{N}\sum_{i=1}^{N}\bigl[\langle S_i\rangle^2\bigr]_{\rm av}$$0 the site-averaged nematic moment (Wiecki et al., 2021).

2. Physical Significance Across Systems

A finite $$q_{\rm EA}(L,T) = \frac{1}{N}\sum_{i=1}^{N}\bigl[\langle S_i\rangle^2\bigr]_{\rm av}$$1 signifies rigid local order that is locked to the randomness of the disorder rather than to a uniform global symmetry-broken state. In spin glasses, nonzero $$q_{\rm EA}(L,T) = \frac{1}{N}\sum_{i=1}^{N}\bigl[\langle S_i\rangle^2\bigr]_{\rm av}$$2 below the critical temperature $$q_{\rm EA}(L,T) = \frac{1}{N}\sum_{i=1}^{N}\bigl[\langle S_i\rangle^2\bigr]_{\rm av}$$3 demarcates the spin-glass phase and reflects frozen-in randomness in the spin configuration (Lewenstein et al., 2022). In the Bose glass, $$q_{\rm EA}(L,T) = \frac{1}{N}\sum_{i=1}^{N}\bigl[\langle S_i\rangle^2\bigr]_{\rm av}$$4 arises from localized superfluid puddles pinned by disorder, contrasting the uniform Mott insulator and clean superfluid, where $$q_{\rm EA}(L,T) = \frac{1}{N}\sum_{i=1}^{N}\bigl[\langle S_i\rangle^2\bigr]_{\rm av}$$5 (Thomson et al., 2016). In random-field Ising-nematic systems, a finite $$q_{\rm EA}(L,T) = \frac{1}{N}\sum_{i=1}^{N}\bigl[\langle S_i\rangle^2\bigr]_{\rm av}$$6 directly correlates with the NMR spectral width due to local nematic domains (Wiecki et al., 2021).

Characteristic features:

• $$q_{\rm EA}(L,T) = \frac{1}{N}\sum_{i=1}^{N}\bigl[\langle S_i\rangle^2\bigr]_{\rm av}$$7 in clean, uniform phases or above $$q_{\rm EA}(L,T) = \frac{1}{N}\sum_{i=1}^{N}\bigl[\langle S_i\rangle^2\bigr]_{\rm av}$$8.
• $$q_{\rm EA}(L,T) = \frac{1}{N}\sum_{i=1}^{N}\bigl[\langle S_i\rangle^2\bigr]_{\rm av}$$9 only when disorder induces persistent inhomogeneities or glassy clusters.
• At criticality or near the transition, $P(q)$0 vanishes as a critical exponent.

3. Measurement Methodologies for Finite-Size q_{EA}

Spin Glasses

Monte Carlo simulation protocols for the 3D Edwards-Anderson model utilize parallel tempering with several thousand disorder samples. The overlap $P(q)$1 is computed from thermally independent replicas, the $P(q)$2 distribution is histogrammed, and $P(q)$3 is defined as the location of the main peak at $P(q)$4 (Yucesoy et al., 2012). The value of $P(q)$5 decreases with system size in the EA model, with reported values at $P(q)$6 for $P(q)$7 to $P(q)$8 ranging from $P(q)$90.65 to $\alpha$00.45.

Quantum/Bosonic Systems

In the Bose glass, quantum gas microscopes can detect the EA parameter through single-shot parity-resolved snapshots. By reconstructing local densities $\alpha$1 over many disorder patterns, one computes

$\alpha$2

Approximate $\alpha$3 estimates are possible using circular averaging for a single disorder realization, but with systematic biases (Thomson et al., 2016).

Field-Theoretic and Analytical Approaches

Saddle-point analysis in the Haake-Lewenstein-Wilkens framework expresses $\alpha$4 in terms of the squared local magnetization in an effective mean-field theory, yielding $\alpha$5 with $\alpha$6 from a bulk self-consistency equation. In three and four dimensions, $\alpha$7 for $\alpha$8, but in two dimensions the method overestimates $\alpha$9 due to strong boundary effects (Lewenstein et al., 2022).

Lee–Yang Zeros

The density and scaling of overlap-field Lee–Yang zeros provide an independent estimator: the initial slope of the integrated zero density $\beta$0 in the spin-glass phase yields $\beta$1, and this agrees within 5% with overlap-based estimates up to $\beta$2 (Baños et al., 2012).

Eigenstate Spin-Glass Order Parameter (ESG)

For many-body localized spin-glass order, reduced density matrices of two sites characterize $\beta$3 analogues in both static and dynamic regimes. The conventional $\beta$4 is compared with the ESG parameter derived from two-site reduced density matrices; both track each other in disorder-driven transitions and in finite-size scaling (Javanmard et al., 2018).

4. Finite-Size Corrections and Scaling

Finite-size $\beta$5 exhibits downward drift with increasing $\beta$6; no explicit scaling exponent for the correction is universally reported (Yucesoy et al., 2012). In glassy random energy models or 1/f noise models, the exact representation yields

$\beta$7

with $\beta$8 analytic in the ergodic phase and singular at the glass transition (Cao et al., 2016). In EA spin glasses, the broadening of the main $\beta$9 peak narrows and sharpens with $q = \frac{1}{N}\sum_{i=1}^N S_i^\alpha S_i^\beta$0, but the position $q = \frac{1}{N}\sum_{i=1}^N S_i^\alpha S_i^\beta$1 scales slowly; in the SK model, many additional peaks at lower $q = \frac{1}{N}\sum_{i=1}^N S_i^\alpha S_i^\beta$2 emerge characteristic of full RSB, but in the 3D EA model, $q = \frac{1}{N}\sum_{i=1}^N S_i^\alpha S_i^\beta$3 for peaks at small $q = \frac{1}{N}\sum_{i=1}^N S_i^\alpha S_i^\beta$4 remains small and flat, supporting a two-state scenario.

5. System-Dependent Behavior and Experimental Probes

System	Definition of $q = \frac{1}{N}\sum_{i=1}^N S_i^\alpha S_i^\beta$5	Experimental Proxy/Observable
3D Ising EA Spin Glass	$q = \frac{1}{N}\sum_{i=1}^N S_i^\alpha S_i^\beta$6 peak position	Overlap distribution from MC simulations
Disordered Bose-Hubbard	Density variance across sites	QGM snapshots, parity-resolved imaging
Random-Field Ising-Nematic	$q = \frac{1}{N}\sum_{i=1}^N S_i^\alpha S_i^\beta$7	NMR spectral width, local moments
Circular 1/f-Noise Model	Replica/Jack polynomial sum	Moments of Gibbs weights, finite-$q = \frac{1}{N}\sum_{i=1}^N S_i^\alpha S_i^\beta$8 scaling

NMR experiments measure $q = \frac{1}{N}\sum_{i=1}^N S_i^\alpha S_i^\beta$9, directly probing $q_{\rm EA}(L,T) = \arg\max_{q>0} P(q)$0 above the nematic transition (Wiecki et al., 2021). In ultracold atom experiments, the BG phase can be unambiguously detected by the rapid onset of $q_{\rm EA}(L,T) = \arg\max_{q>0} P(q)$1 with disorder, sharply discriminating it from MI and SF states (Thomson et al., 2016).

6. Controversies and Limitations

The Haake-Lewenstein-Wilkens approach, while deriving $q_{\rm EA}(L,T) = \arg\max_{q>0} P(q)$2 analytically for any dimension, overestimates the existence of a finite-$q_{\rm EA}(L,T) = \arg\max_{q>0} P(q)$3 spin-glass phase in $q_{\rm EA}(L,T) = \arg\max_{q>0} P(q)$4 due to neglect of dominant boundary or domain-wall effects, highlighting the limitations of mean-field approximations in low dimensions (Lewenstein et al., 2022). In 3D, both numerical and analytical work support finite $q_{\rm EA}(L,T) = \arg\max_{q>0} P(q)$5 below $q_{\rm EA}(L,T) = \arg\max_{q>0} P(q)$6 and a two-state (droplet/chaotic-pairs) scenario, in sharp contrast to the infinite-state structure of the SK model (Yucesoy et al., 2012). The finite-size behavior, slow drift with $q_{\rm EA}(L,T) = \arg\max_{q>0} P(q)$7, and system-specific corrections remain the subject of ongoing investigation. In quantum glass settings, the distinction between thermal and disorder-induced fluctuations in finite $q_{\rm EA}(L,T) = \arg\max_{q>0} P(q)$8 is subtle and demands high-precision, disorder-ensemble-resolved measurement protocols (Thomson et al., 2016, Javanmard et al., 2018).

7. Broader Implications and Extensions

A finite EA order parameter functions as a universal marker of quenched glassy order across a diverse range of systems, including spin glasses, bosonic glasses, random-field nematics, and log-correlated disordered landscapes (log-REM, 1/f-noise). In each context, $q_{\rm EA}(L,T) = \arg\max_{q>0} P(q)$9 encodes the persistent memory of frozen, disorder-pinned inhomogeneities and underpins the experimental identification of glass phases. In the quantum domain, $$q_{\rm EA} = \overline{\,\langle \hat n_i\rangle^2\,} - \overline{\langle \hat n_i\rangle}^2 \equiv \overline{\bigl(\langle \hat n_i\rangle - \overline{\langle \hat n_i\rangle}\bigr)^2}$$0 and its ESG generalizations provide scalable diagnostics that extend into non-equilibrium and many-body localized regimes, with potential application in quantum simulation platforms (Javanmard et al., 2018). The rigorous determination of $$q_{\rm EA} = \overline{\,\langle \hat n_i\rangle^2\,} - \overline{\langle \hat n_i\rangle}^2 \equiv \overline{\bigl(\langle \hat n_i\rangle - \overline{\langle \hat n_i\rangle}\bigr)^2}$$1 and its systematics remains central to the understanding of glassy matter, replica symmetry breaking, and the interplay between disorder, interactions, and topology in complex materials.