Nonreciprocal Disorder Prevents Zero-Temperature Freezing in a Ferromagnet

Technical Summary: Nonreciprocal Disorder Prevents Zero-Temperature Freezing in a Ferromagnet

Problem Statement
While nonreciprocal (NR) interactions are known to drive systems out of equilibrium and generate cyclic currents, the specific consequences of quenched (spatially fixed, random) nonreciprocity in extended, finite-dimensional systems remain largely unexplored. In equilibrium statistical mechanics, the competition between quenched disorder and ferromagnetism (e.g., in Random-Field or Random-Bond Ising models) leads to phenomena such as rare-region effects, pinning, and zero-temperature freezing, where the system becomes trapped in a metastable state. The central question addressed is how these phenomena are modified when the disorder is nonreciprocal. Specifically, does quenched nonreciprocity destroy ferromagnetic order, and does it prevent the system from freezing at zero temperature ($T=0$)?

Methodology
The authors introduce a minimal nonreciprocal analogue of the $\pm J$ Random-Bond Ising Model (RBIM) on a 2D square lattice. The model consists of Ising spins $\sigma_i \in \{-1, +1\}$ with nearest-neighbor couplings $(J_{ij}, J_{ji})$ drawn independently:

• With probability $1-p$: Reciprocal ferromagnetic bonds $(+1, +1)$.
• With probability $p/2$: Nonreciprocal bonds $(+1, -1)$.
• With probability $p/2$: Nonreciprocal bonds $(-1, +1)$.

The system evolves via single-spin-flip Glauber dynamics at temperature $T$. At $T=0$, spins flip deterministically if the energy cost $\Delta E_i < 0$, do not flip if $\Delta E_i > 0$, and flip with probability $1/2$ if $\Delta E_i = 0$.

The study employs three complementary approaches:

1. Numerical Simulation: Large-scale Monte Carlo simulations ($N \times N$ lattices up to $512^2$) to characterize magnetization, spatial/temporal correlation functions, and domain growth.
2. Mean-Field Theory: A single-site mean-field treatment to derive the phase diagram and critical density $p_c$.
3. Analytical Arguments: A gauge-invariance argument to establish an upper bound on the critical density and a field-theoretic coarse-graining to analyze the relevance of nonreciprocity near the critical point.

Key Contributions and Results

1. Destruction of Ferromagnetic Order:
   The system exhibits a continuous nonequilibrium phase transition from an ordered ferromagnetic phase to a disordered paramagnetic phase. The critical temperature $T_c(p)$ decreases monotonically with the density of nonreciprocal bonds $p$, reaching zero at a finite critical density $p_c \approx 0.20$.

   • Upper Bound: Using a gauge-invariance argument, the authors prove that long-range order requires the percolation of reciprocal bonds. Since reciprocal bonds must form a connected cluster to sustain order, and the bond percolation threshold on a 2D square lattice is $1/2$, they establish the rigorous bound $p_c(T) \le 1/2$.
   • Phase Diagram: The phase diagram resembles that of the $\pm J$ RBIM at finite temperatures, with ordered and disordered phases separated by a second-order transition. Mean-field theory predicts $p_c^{MF} = 1/3$, qualitatively consistent with simulations.

2. Prevention of Zero-Temperature Freezing:
   Unlike equilibrium disordered models (such as the RBIM), which freeze into a spin-glass or metastable state at $T=0$ due to "pinning" of domain walls, the disordered phase in this nonreciprocal model remains dynamically active at $T=0$.

   • Mechanism: The authors identify an asymmetry between local and global energy. A spin flip can decrease the local "selfish" energy ($E_i$) deterministically while increasing or leaving unchanged the global energy ($E$). This allows for "Escher-like" infinite descent in local energy without violating global constraints.
   • Persistent Dynamics: Specific configurations of nonreciprocal bonds (e.g., $2 \times 2$ plaquettes with zero boundary field) support deterministic cycles. These act as sources of localized, athermal driving, preventing any spin from becoming permanently frozen. Consequently, the system does not reach an absorbing state.

3. Rare-Region Effects and Coarsening:

   • Rare-Region Reversals: In the disordered phase at $T=0$, the autocorrelation function $C_t(t)$ exhibits stretched-exponential tails ($C_t(t) \sim \exp[-(t/\tau)^{\alpha}]$ with $\alpha < 1$). These tails arise from rare regions with locally low nonreciprocal density that sustain transient local order. Crucially, unlike in equilibrium where such regions freeze at $T=0$, here they continue to reverse due to nonreciprocal driving.
   • Activated Coarsening: In the ordered phase, domain growth is initially curvature-driven ($L(t) \sim t^{1/2}$). At later times, nonreciprocal bonds act as pinning sites, slowing growth to an activated form ($L(t) \sim \ln t$). However, unlike the RBIM where pinning halts growth entirely at $T=0$, the nonreciprocal driving ensures that fully absorbing configurations are exponentially unlikely, allowing coarsening to proceed indefinitely.

Significance
The paper demonstrates that quenched nonreciprocity fundamentally alters the low-temperature behavior of ferromagnetic systems. While finite-temperature behavior remains similar to equilibrium disordered models (due to the irrelevance of nonreciprocity near the clean-Ising fixed point in the renormalization group sense), the zero-temperature limit is qualitatively distinct. The authors show that nonreciprocal interactions supply a "run-and-chase" cyclic source of deterministic athermal driving that prevents freezing. This suggests that the destruction of fixed points and the emergence of persistent dynamics, previously observed in fully-connected (mean-field) nonreciprocal systems, also persist in finite-dimensional spatially extended systems. The work establishes that spatially quenched nonreciprocity can frustrate large-scale order while simultaneously preventing the system from settling into a static, frozen state.