[论文解读] Low energy excitations in a long prism geometry: computing the lower critical dimension of the Ising spin glass

Massimo Bernaschi, Luis A. Fernandez+4 more|arXiv (Cornell University)|Jan 12, 2026

Theoretical and Computational Physics被引用 0

一句话总结

论文引入矩形棱柱几何来研究伊辛自旋玻璃中的低能激发，计算在D=3中的下临界维数，并通过大规模模拟（带开放边界条件和Houdayer簇移动）分析相关函数的多重分形光谱。

ABSTRACT

We propose a general method for studying systems that display excitations with arbitrarily low energy in their low-temperature phase. We argue that in a rectangular right prism geometry, with longitudinal size much larger than the transverse size, correlations decay exponentially (at all temperatures) along the longitudinal dimension, but the scaling of the correlation length with the transverse size carries crucial information from which the lower critical dimension can be inferred. The method is applied in the particularly demanding context of Ising spin glasses at zero magnetic field. The lower critical dimension and the multifractal spectrum for the correlation function are computed from large-scale numerical simulations. Several technical novelties (such as the unexpectedly crucial performance of Houdayer's cluster method or the convenience of using open - rather than periodic - boundary conditions) allow us to study three-dimensional prisms with transverse dimensions up to $L=24$ and effectively infinite longitudinal dimensions down to low temperatures. The value that we find for the lower critical dimension turns out to be in agreement with expectations from both the Replica Symmetry Breaking theory and the Droplet model for spin glasses. We argue that our novel setting holds promise in clarifying which of the two competing theories more accurately describes three-dimensional spin glasses.

研究动机与目标

• Develop a general method to study systems with arbitrarily low-energy excitations in their low-temperature phase.
• Apply the method to the 3D Ising spin glass at zero field to extract the lower critical dimension.
• Compute the multifractal spectrum of the correlation function in the prism geometry.
• Demonstrate numerical strategies (open boundaries, Houdayer clusters) enabling large-scale equilibrium studies.
• Compare results with Replica Symmetry Breaking and Droplet model predictions.

提出的方法

• Use a rectangular right prism with longitudinal length M>>L to study low-energy excitations.
• Impose mutually incongruent boundary conditions across prism ends to extract excess free energy and correlation lengths.
• Relate ΔF scaling to L and M via ΔF ~ L^{D-1}/M^{b} with b conjectured to be D-independent (b=3/2 from MFT).
• Define and compute equilibrium correlation functions C^{(n)}(z1,z2) from replica overlaps Q(z).
• Employ a 1D toy model to interpret scaling: ξ_n ~ 1/E_k^{odd} and τ_n spectrum relations.
• Leverage Houdayer cluster moves and open boundary conditions to accelerate equilibration in long prisms.]
• research_questions: [

实验结果

研究问题

• RQ1 What is the lower critical dimension D_lc for the Ising spin glass at zero magnetic field in 3D?
• RQ2 How does the longitudinal correlation length ξ(L) scale with transverse size L in a long prism geometry?
• RQ3 What is the multifractal spectrum τ_n of the correlation functions in this geometry, and how does it compare to 1D toy models and DM/MFT predictions?
• RQ4 Do open boundary conditions provide a reliable and efficient route to equilibrium compared to periodic boundaries in this setting?
• RQ5 Can the results distinguish between Replica Symmetry Breaking and Droplet model descriptions of 3D spin glasses?

主要发现

L	M	BC-Z	ξ_{n=1}(T=0.7)	ξ_{n=1}(T≈T_c)
4	192	PBC	8.34(3)	4.03(5)
6	192	PBC	13.41(4)	5.87(3)
8	192	PBC	18.69(8)	7.83(3)
12	320	PBC	30.73(9)	11.68(4)
16	512	PBC	44.2(7)	15.8(3)
16	48	OBC	44.4(3)	15.74(11)
24	88	OBC	73.1(16)	24.1(5)

• ξ_{n=1} scales as a power law with L at low T and linearly with L at T_c, supporting a finite D_lc in the reported range.
• Best fit yields a_{3}=1.34(3)[3] and D_lc=2.49(3)[3], agreeing with mean-field predictions (a_3^{MFT}=4/3, D_lc^{MFT}=5/2).
• The multifractal spectrum τ_n is found to satisfy τ_n<n with τ_2≈1.43, τ_3≈1.71, τ_4≈1.90–1.92 across sizes, indicating strong multifractality.
• Open boundary conditions with L up to 24 and M large allow reliable extraction of ξ_{n=1} with manageable finite-size effects.
• Houdayer’s cluster moves yield a dramatic speed-up (~10^3) enabling equilibration of large prisms at T ≈ 0.63 T_c.
• Results favor mean-field-like scaling and are in good agreement with previous zero-temperature and finite-temperature estimates for D_lc.