alphaFeedback

WorksMathematics

Bloch's Law for Finite-Range Heisenberg Ferromagnets in Three Dimensions

AI contributions
  • unreleased internal OpenAI model
View details
Ancillary data
4 proof or workflow links
View

Abstract

We prove Bloch's T3/2 law for the spontaneous magnetization of three-dimensional quantum Heisenberg ferromagnets at every fixed positive quantum spin. The result holds for every nonnegative symmetric finite-range interaction whose support generates ℤ3, including spatially anisotropic couplings. The leading magnetization deficit has the exact coefficient determined by the determinant of the quadratic one-magnon dispersion. The thermodynamic limit is taken before the right field derivative at zero, and the low-temperature limit is taken last.

Introduction

Bloch’s law predicts that the spontaneous magnetization of a three-dimensional ferromagnet falls below its fully polarized value by an amount proportional to T3/2T^{3/2} at low temperature. We prove this law for the quantum Heisenberg ferromagnet at every fixed spin, allowing any nonnegative symmetric finite-range interaction whose support generates the lattice. The interaction may be spatially anisotropic and need not contain the ordinary nearest-neighbor bonds. Its quadratic dispersion determines the leading coefficient.

The model and the result

Fix S∈{12,1,32,…}S \in\left\{\frac{1}{2},1,\frac{3}{2},\ldots\right\}. The spin space at each site is HS=C2S+1\mathcal{H}_S=\mathbb{C}^{2S+1}, with orthonormal basis (em)m=−SS(e_m)_{m=-S}^{S} and operators

Szem=mem,S+em=S(S+1)−m(m+1) em+1(m<S),S+eS=0.S^z e_m = m e_m,\qquad S^+e_m=\sqrt{S(S+1)-m(m+1)}\,e_{m+1}\quad(m<S),\qquad S^+e_S=0.

Set S−=(S+)∗S^-=(S^+)^*, Sx=(S++S−)/2S^x=(S^++S^-)/2, and Sy=(S+−S−)/(2i)S^y=(S^+-S^-)/(2i). A subscript specifies the site where an operator acts. Let J:Z3→[0,∞)J:\mathbb{Z}^3\to[0,\infty) have finite support and satisfy

J(0)=0,J(z)=J(−z),⟨z:J(z)>0⟩Z=Z3.(1)J(0)=0,\qquad J(z)=J(-z),\qquad\langle z:J(z)>0\rangle_{\mathbb{Z}}=\mathbb{Z}^3. \tag*{(1)}

On ΛN=[−N,N]3∩Z3\Lambda_N=[-N,N]^3\cap\mathbb{Z}^3, with free boundary, define

HS,NJ=−12∑x,y∈ΛNJ(y−x)Sx⋅Sy,MS,N=∑x∈ΛNSxz.(2)H^J_{S,N}=-\frac{1}{2}\sum_{x,y\in\Lambda_N}J(y-x)S_x\cdot S_y,\qquad M_{S,N}=\sum_{x\in\Lambda_N}S_x^z. \tag*{(2)}

The sum counts each unordered bond twice; the factor 1/21/2 removes that duplication. We set the lattice spacing, Boltzmann constant, and ℏ\hbar equal to one, so T=β−1T=\beta^{-1}.

For β>0\beta>0 and h∈Rh\in\mathbb{R}, put

pβ(h)=lim⁡N→∞1β∣ΛN∣log⁡Tr⁡e−β(HS,NJ−hMS,N),mS,J(β)=∂h+pβ(0).(3)p_\beta(h)=\lim_{N\to\infty}\frac{1}{\beta|\Lambda_N|}\log\operatorname{Tr}e^{-\beta(H^J_{S,N}-hM_{S,N})},\qquad m_{S,J}(\beta)=\partial_h^+p_\beta(0). \tag*{(3)}

The finite-range pressure limit exists and is convex and SS-Lipschitz; Section 7 recalls the proof. Rotation by π\pi about the xx-axis makes the pressure even, so mS,J(β)∈[0,S]m_{S,J}(\beta)\in[0,S]. The derivative selects the magnetization in a field tending to zero through positive values after the thermodynamic limit.

Define the positive definite matrix

DS,J=S2∑z∈Z3J(z)zzT.(4)D_{S,J}=\frac{S}{2}\sum_{z\in\mathbb{Z}^3}J(z)zz^{\mathsf T}. \tag*{(4)}

It is the quadratic coefficient of the one-magnon dispersion:

εS,J(k)=S∑zJ(z)(1−cos⁡(k⋅z))=kTDS,Jk+O(∣k∣4).\varepsilon_{S,J}(k)=S\sum_z J(z)(1-\cos(k\cdot z))=k^{\mathsf T}D_{S,J}k+O(|k|^4).

Positive definiteness follows because the interaction support spans R3\mathbb{R}^3.

Theorem 1.1 (Bloch’s law). For every fixed S∈{12,1,32,…}S\in\left\{\frac{1}{2},1,\frac{3}{2},\ldots\right\} and every interaction satisfying (1.1),

lim⁡β→∞β3/2(S−mS,J(β))=ζ(3/2)8π3/2det⁡DS,J.\lim_{\beta\to\infty}\beta^{3/2}\bigl(S-m_{S,J}(\beta)\bigr)=\frac{\zeta(3/2)}{8\pi^{3/2}\sqrt{\det D_{S,J}}}.

For J(±ej)=1J(\pm e_j)=1 and all other couplings zero, DS,J=SID_{S,J}=SI; the theorem gives the usual nearest-neighbor coefficient ζ(3/2)/(8π3/2S3/2)\zeta(3/2)/(8\pi^{3/2}S^{3/2}). The order of limits is part of the assertion: first volume tends to infinity at fixed β,h\beta,h, then h↓0h\downarrow0, and finally β→∞\beta\to\infty, with S,JS,J fixed.

Historical context and dependencies

Bloch derived the T3/2T^{3/2} reduction from long-wavelength spin waves [2]. The bosonic representation of Holstein and Primakoff [10] and Dyson’s analysis of spin-wave interactions [6, 7] developed the low-temperature expansion. In particular, Dyson identified the interaction contribution to the magnetization at order T4T^4; his estimates of the remaining higher-order terms were not a rigorous fixed-spin proof of the asymptotic law [7], Section 9. The leading coefficient in (1.6) is obtained from the low-temperature density of ideal magnons with quadratic energy kTDS,Jk\mathbf{k}^{\mathsf{T}}D_{S,J}\mathbf{k}. The issue is to justify this coefficient for the interacting spin system and the field derivative in (3).

Rigorous ordering theory distinguishes sharply between dimensions and between classical and quantum models. Mermin and Wagner excluded positive-temperature spontaneous magnetic order for the one- and two-dimensional isotropic Heisenberg models under their finite-range hypotheses [12]. Infrared bounds proved symmetry breaking for classical models [9], and Dyson, Lieb, and Simon established ordering for classes of quantum antiferromagnets [8]. These developments left the standard three-dimensional quantum ferromagnet as a separate problem, recorded by Lieb [11] and still discussed as open in recent accounts [16, 1]. The distinct companion manuscript [15], Theorem 1.1 proves low-temperature magnetized equilibrium states for the nearest-neighbor quantum ferromagnet. We use its finite-set algebra and deterministic lattice geometry, with their precise hypotheses stated at the points of use.

The free energy has a different rigorous history. Conlon and Solovej [4] and Tóth [17] obtained bounds for spin 1/21/2. Correggi, Giuliani, and Seiringer proved the leading low-temperature free-energy asymptotic for the nearest-neighbor model on Z3\mathbb{Z}^3 at every fixed spin [5], Theorem 2.1. A zero-field free-energy estimate does not by itself control the one-sided field derivative that defines spontaneous magnetization. Our proof therefore estimates the observable through a representation that remains uniform under arbitrary extra up-spin constraints, and only then passes to pressure increments.

The probabilistic representation belongs to the random-interchange approach of Tóth [17]; the symmetric-slot realization of higher spin is described by Nachtergaele [13]. The pin argument uses the stable-polynomial approach to negative dependence developed by Borcea, Brändén, and Liggett [3]. Diffusive smoothing uses the energy method of Nash [14]. The companion [15] provides the particular stable exchange symbols, presence inequality, local lattice inequality, and routed unit flows used below. We prove here the finite-range pin bootstrap, the quantitative long-return bound, and the comparison that recovers the full dispersion.

How returns determine the coefficient

Represent each spin by ℓ=2S\ell=2S spin-1/21/2 slots. Independent Poisson exchanges between slots, followed by a uniform permutation of the slots at each site, form directed cycles over a time period of length β\beta. A cycle has a single up or down color. Constraining all cycles through the boundary and through specified slot-time points to be up defines a pin law. If EE is its pin set and ii is an unpinned slot-time point, a change of weights gives

dE(i)=r(E,{i})1+r(E,{i}).d_E(i)=\frac{r(E,\{i\})}{1+r(E,\{i\})}.

where dE(i)d_E(i) is its down probability. Under the law with ii also pinned, r(E,{i})r(E,\{i\}) is the probability that the directed line from ii returns to ii at positive elapsed time before meeting EE. Section 2 makes this identity precise and represents rr by a walk on the slots left after exposing the cycles missing the augmented pin set. The exposed environment repeats over periods, while the walk uses fresh transition randomness on every traversal.

There are two geometric tasks. For rough estimates, choose three linearly independent interaction vectors. They generate a sublattice of finite index; on each of its cosets, the chosen bonds form a copy of the ordinary cubic lattice. Section 3 transfers local Nash and flow estimates to these copies and compares their energies with the full interaction energy. Section 4 then proves a sparse-pin estimate by downward induction on the pin set. This induction supplies enough regular time for diffusion without assuming a finite-range ordering theorem.

Sharp asymptotics require all interaction bonds. The ideal walk has site jump rates SJ(z)S J(z), and its heat kernel ptS,Jp_t^{S,J} satisfies

ptS,J(x,x)∼κt−3/2,κ=(4π)−3/2(det⁡DS,J)−1/2.p_t^{S,J}(x,x) \sim\kappa t^{-3/2}, \qquad\kappa= (4\pi)^{-3/2}(\det D_{S,J})^{-1/2}.

After a full traversal the slot index is uniform, so the ideal nn-period return probability to a specified slot is ℓ−1pnβS,J(x,x)\ell^{-1}p_{n\beta}^{S,J}(x,x). Summing over the ℓ\ell slots at a site cancels this factor, leaving the finite sums κ∑n≤Kn−3/2\kappa\sum_{n\le K} n^{-3/2} that approach the Bloch coefficient. The selected sublattice is used only for estimates; the full jump kernel determines this coefficient.

The main analytic difficulty is controlling arbitrarily late returns in the same repeated environment. Section 5 truncates a quadratic hitting-probability identity after KK periods. A mass already spread for time comparable to KβK\beta has a flow to infinity with energy at most C(Kβ)−1/2C(K\beta)^{-1/2}. Averaging over one further period gives a remainder bounded by CK−1/2β−3/2+oK(β−3/2)C K^{-1/2}\beta^{-3/2}+o_K(\beta^{-3/2}), with CC independent of KK. Its K=1K=1 case first improves the exposed-slot density to O(β−3/2)O(\beta^{-3/2}).

This improved density permits the comparison in Section 6. The ideal path rarely encounters an exposed strand over a fixed number of periods. Splitting a return at its midpoint combines that small encounter probability with diffusive smoothing, giving an error oK(β−3/2)o_K(\beta^{-3/2}) in each return entry. The comparison retains the same exposed environment in both halves; independence of their environments is unnecessary. The finite-return comparison and the long-return remainder give matching upper and lower down-probability bounds, with the lower bound requiring a pin-free neighborhood of the starting point.

Finally, Section 7 expresses a positive field as a mixture of extra pins. The probability of selecting any given slot is at most 1−e−βh1-e^{-\beta h}, so the required neighborhood is free of extra pins with probability tending to one as h↓0h\downarrow0 at fixed β,K\beta,K. Integrating the finite-volume derivative bounds before taking the thermodynamic limit proves (1.6) in the stated order. The uniform pin estimates, the finite-index geometric reduction, and the energy bound for a spread source are the reusable parts of the argument.

Exchange pictures and pinned returns

We express the magnetization through colored exchange cycles, then turn pinning probabilities into one-particle return probabilities. Two changes of law are involved. Pinning removes the color factor from cycles that meet the pins. Exposing the cycles that avoid the pins leaves a conditional exchange process on the remaining slots. The last part of this section identifies its first return moment with an unconditioned Markov walk in the exposed environment. These constructions are finite-dimensional; the spatial estimates begin in Section 3.

Spin slots and the exchange law

Fix β>0\beta>0 and a cube Λ=ΛN\Lambda=\Lambda_N. Set ℓ=2S\ell=2S and

b0=max⁡{∣z∣∞:J(z)>0},∂Λ={x∈Λ:dist⁡∞(x,Z3∖Λ)≤b0},Λ∘=Λ∖∂Λ.b_0=\max\{|z|_\infty:J(z)>0\},\qquad\partial\Lambda=\{x\in\Lambda:\operatorname{dist}_\infty(x,\mathbb{Z}^3\setminus\Lambda)\le b_0\},\qquad\Lambda^\circ=\Lambda\setminus\partial\Lambda.

No interaction bond joins Λ∘\Lambda^{\circ} to the exterior of Λ\Lambda. The slots are the elements of VΛ=Λ×{1,…,ℓ}V_{\Lambda}=\Lambda\times\{1,\ldots,\ell\}. Above each unordered bond {x,y}⊂Λ\{x,y\}\subset\Lambda with J(y−x)>0J(y-x)>0, include all ℓ2\ell^2 unordered slot pairs. Write EΛ\mathcal{E}_{\Lambda} for the resulting slot edges and set

References

  1. [1]J. E. Björnberg and K. Ryan. Exponential decay in 𝔸(n)-invariant quantum spin systems. Annales Henri Poincaré (2026), published online March 25, 2026. doi:10.1007/s00023-026-01690-3.DOI
  2. [2]F. Bloch. Zur Theorie des Ferromagnetismus. Zeitschrift für Physik 61 (1930), 206–219. doi:10.1007/BF01339661.DOI
  3. [3]J. Borcea, P. Brändén, and T. M. Liggett. Negative dependence and the geometry of polynomials. Journal of the American Mathematical Society 22 (2009), 521–567. doi:10.1090/S0894-0347-08-00618-8.DOI
  4. [4]J. G. Conlon and J. P. Solovej. Upper bound on the free energy of the spin 1/2 Heisenberg ferromagnet. Letters in Mathematical Physics 23 (1991), 223–231. doi:10.1007/BF01885500.DOI
  5. [5]M. Correggi, A. Giuliani, and R. Seiringer. Validity of the spin-wave approximation for the free energy of the Heisenberg ferromagnet. Communications in Mathematical Physics 339 (2015), 279–307. doi:10.1007/s00220-015-2402-0.DOI
  6. [6]F. J. Dyson. General theory of spin-wave interactions. Physical Review 102 (1956), 1217–1230. doi:10.1103/PhysRev.102.1217.DOI
  7. [7]F. J. Dyson. Thermodynamic behavior of an ideal ferromagnet. Physical Review 102 (1956), 1230–1244. doi:10.1103/PhysRev.102.1230.DOI
  8. [8]F. J. Dyson, E. H. Lieb, and B. Simon. Phase transitions in quantum spin systems with isotropic and nonisotropic interactions. Journal of Statistical Physics 18 (1978), 335–383. doi:10.1007/BF01106729.DOI
  9. [9]J. Fröhlich, B. Simon, and T. Spencer. Infrared bounds, phase transitions and continuous symmetry breaking. Communications in Mathematical Physics 50 (1976), 79–95. doi:10.1007/BF01608557.DOI
  10. [10]T. Holstein and H. Primakoff. Field dependence of the intrinsic domain magnetization of a ferromagnet. Physical Review 58 (1940), 1098–1113. doi:10.1103/PhysRev.58.1098.DOI
  11. [11]E. H. Lieb. Long range order for the quantum Heisenberg model. Open Problems in Mathematical Physics, contribution dated January 27, 1999, Problem A. Author’s problem statement.
  12. [12]N. D. Mermin and H. Wagner. Absence of ferromagnetism or antiferromagnetism in one- or two-dimensional isotropic Heisenberg models. Physical Review Letters 17 (1966), 1133–1136; erratum, 17 (1966), 1307. doi:10.1103/PhysRevLett.17.1133.DOI
  13. [13]B. Nachtergaele. Quasi-state decompositions for quantum spin systems. In B. Grigelionis et al. (eds.), Probability Theory and Mathematical Statistics, Proceedings of the Sixth Vilnius Conference, VSP/TEV, Utrecht–Tokyo–Vilnius, 1994, pp. 565–590. arXiv:cond-mat/9312012.arxiv.org/abs/cond-mat/9312012
  14. [14]J. Nash. Continuity of solutions of parabolic and elliptic equations. American Journal of Mathematics 80 (1958), 931–954. doi:10.2307/2372841.DOI
  15. [15]OpenAI. Spontaneous magnetization in the quantum Heisenberg ferromagnet. OpenAI Math Release preprint OAI:Spontaneous-magnetization-in-the-quantum-Heisenberg-ferromagnet-September-24-2026, 2026.
  16. [16]R. Seiringer. The Heisenberg ferromagnet: a dilute Bose gas in disguise. In Mini-Workshop: New Directions in Correlated Quantum Systems, Oberwolfach Reports 22 (2025), 357–358. doi:10.4171/OWR/2025/7.
  17. [17]B. Tóth. Improved lower bound on the thermodynamic pressure of the spin 1/2 Heisenberg ferromagnet. Letters in Mathematical Physics 28 (1993), 75–84. doi:10.1007/BF00739568.DOI

Paper details

vlnc:W_d7r7u4nt
v1 · latestHistory (1)

published · Fresh — published within the last 30 days

Usage

views · 1

downloads · 0

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

LicenseRef-OpenAI-public

DOI · 10.83295/628k-3s2y

Original source

  • Mathematics

unreleased internal OpenAI model · View statement →

Looking for related papers…
Contents