Paper XL: Lattice Hopf Field Theory — Numerical Methods, Algorithm Design, and Initial Benchmarks

Abstract

We present the first lattice formulation of the Faddeev-Niemi (FN) Hopf field theory -- the sigma model $\mathbf{n}: \mathbb{R}^3 \to S^2$ with Dirichlet and quartic Skyrme terms, whose solitons carry integer-valued Hopf charge $H \in \pi_3(S^2) = \mathbb{Z}$. The continuum CP$^1$ parametrization $\mathbf{n} = z^\dagger \boldsymbol{\sigma} z$ is discretized on a cubic lattice with link variables $U_\mu(x) = z^\dagger(x) z(x{+}\hat{\mu}) / |z^\dagger(x) z(x{+}\hat{\mu})|$ that preserve the topological structure under discretization. We compute the lattice Hopf invariant via both the Whitehead integral and the topologically exact preimage linking number, validating integer-valuedness on the analytic Hopf map down to $N = 8$ (Appendix B.3). A hybrid Monte Carlo algorithm adapted to the CP$^1$ target space ($S^3 \subset \mathbb{C}^2$) is developed, incorporating geodesic leapfrog integration and microcanonical overrelaxation for reduced autocorrelation. We introduce a per-site topology constraint that preserves the Hopf charge $H_W = 1.07 \pm 0.002$ over 200 production sweeps with only $1.3\times$ overhead -- the first topology-preserving Monte Carlo with per-site local constraint in Hopf field theory. Initial production benchmarks at $N = 32$ across $\beta = 5, 10, 20$ show perfect Derrick virial balance $E_2 = E_4$ throughout. Crucially, vacuum subtraction ($H = 0$ reference runs) reveals that the physical soliton energy $E_{\text{sol}} = \langle E \rangle_{H=1} - \langle E \rangle_{H=0}$ vanishes with increasing $\beta$ ($E_{\text{sol}}/E_{\text{BS}} \to 0.03$ at $\beta = 20$): the 3D thermal MC preserves the approximate Hopf winding but does not stabilize a compact soliton. This provides direct computational evidence that the 4D lattice path integral is necessary for quantum soliton stabilization. The three-sector $F_2$ lattice extension is described at the conceptual level; the full lattice action and HMC force derivation are deferred to future work. Topological freezing at large $\beta$ is identified as a practical obstacle, and mitigation strategies are discussed. Strategies for resolving quantum stabilization of the $F_2$ saddle point via a 4D lattice path integral are outlined.