Nonreciprocal Blume-Capel Model

• The model introduces antisymmetric interactions and vacancy imbalance on a spin-1 lattice, yielding diverse phases including disorder, swap oscillations, and static order.
• Mean-field and Monte Carlo analyses reveal critical bifurcations, continuous transitions in 2D consistent with Ising exponents, and first-order endpoints at high coupling.
• The framework offers experimental knobs, such as tuning chemical potentials, to control dynamical phases in spin systems, relevant for colloidal and cold-atom setups.

The nonreciprocal Blume-Capel model (NR-BCM) extends the classical Blume-Capel framework by introducing nonreciprocal (antisymmetric) interactions and chemical-potential imbalance into a two-species, spin-1 lattice system. The model captures the interplay between nonequilibrium nonreciprocity and vacancy energetics, yielding a diverse phase structure not accessible in equilibrium spin systems. It provides a minimal setting for exploring the phenomenology of swapping limit cycles, restoration of static order via local chemistry, and criticality induced by defects in both two and three spatial dimensions (R et al., 21 Dec 2025).

1. Model Specification and Hamiltonian Structure

The NR-BCM is defined on a lattice where each site can be empty or occupied by one of two species, each carrying a spin $\sigma_i^\alpha \in \{-1, 0, +1\}$ with $\alpha \in \{A, B\}$. The system allows for vacancies ($\sigma = 0$) and is characterized by the following non-Hamiltonian, "selfish-energy" for spins of species $\alpha$: $$E_i^\alpha = -J\sum_{j\in\mathrm{nn}(i)}\sigma_i^\alpha\sigma_j^\alpha - K_{\alpha\beta}\sigma_i^\alpha\sigma_i^\beta + \Delta_\alpha(\sigma_i^\alpha)^2$$ where:

• $J>0$ is the intra-species ferromagnetic exchange.
• $K_{AB} = -K_{BA} \equiv K>0$ is an onsite, antisymmetric interspecies coupling introducing maximal nonreciprocity.
• $\Delta_\alpha$ is the single-ion anisotropy or chemical potential for vacancies of species $\alpha$.

For the symmetric case explored, $\Delta_A = -\Delta_B \equiv \Delta$, establishing a chemical-potential (vacancy) imbalance between species.

Glauber-type stochastic dynamics at temperature $T$ govern spin updates: $$\omega_i^\alpha(\sigma \to \sigma') = \frac{ \exp \left(-\beta [ E_i^\alpha(\sigma') - E_i^\alpha(\sigma) ] \right) }{ \sum_{p=-1}^{1} \exp\left(-\beta [ E_i^\alpha(p) - E_i^\alpha(\sigma) ] \right) } \,, \quad \beta = 1/(k_BT)$$

Key observables:

• Magnetizations $M_\alpha = \left<\sigma_i^\alpha\right>$
• Vacancy fractions $v_\alpha = \left<1_{\sigma_i^\alpha = 0}\right>$
• Combined amplitude $R = \sqrt{\frac{1}{2}(M_A^2 + M_B^2)}$
• "Angular momentum"-like oscillation strength $S = M_B\,\dot{M}_A - M_A\,\dot{M}_B$
• Control parameters: $\tilde{J} = 2dJ/k_BT$, $\tilde{K} = K/k_BT$, $\Delta$.

2. Mean-Field Dynamics and Phase Bifurcations

In the spatially uniform mean-field limit, evolution equations for magnetizations are: $$\tau \dot{M}_\alpha = -M_\alpha + \frac{2\sinh(\tilde{J} M_\alpha + \tilde{K}_{\alpha\beta} M_\beta)}{2\cosh(\tilde{J} M_\alpha + \tilde{K}_{\alpha\beta} M_\beta) + e^{\Delta_\alpha / k_BT}}$$ with antisymmetric $\tilde{K}_{AB} = -\tilde{K}_{BA} = \tilde{K}$ and $\Delta_A = -\Delta_B$.

The mean-field analysis identifies three principal dynamical regimes:

• Region I (Disorder): $(M_A, M_B) = (0,0)$
• Region II (Swap/Limit Cycle): persistent time-dependent oscillations of $M_A$, $M_B$ (only for $\tilde{K}\neq0$)
• Region III (Static Order): $(M_A, M_B) \neq (0,0)$, ferromagnetic-like order

Transitions among regions traverse several bifurcations (see (R et al., 21 Dec 2025)):

• Supercritical Hopf (I$\rightarrow$II)
• Saddle-node on invariant circle (SNIC, II$\rightarrow$III)
• Saddle-node of limit cycles (SNLC, oscillations terminate with large $\Delta$)
• Pitchfork (for $\tilde{K}=0$ equilibrium)
• Saddle-node (first-order) for equilibrium transitions
• A cusp-like point where saddle-node lines meet, analogous to a tricritical point

A schematic mean-field phase diagram, with $\Delta$ on the horizontal and $\tilde{J}$ on the vertical axis, demarcates these regimes and bifurcation loci.

3. Monte Carlo Analysis in Two Dimensions

Monte Carlo simulations were conducted on square lattices with periodic boundaries ($L=20$–$120$). At $\Delta=0$, finite systems display droplet-mediated oscillations, but as $L\to\infty$, spiral topological defects proliferate, leading to the destruction of global swap order: $R, S \to 0$ and the system remains disordered.

Introducing a chemical-potential imbalance, $\Delta \gtrsim 1$, promotes robust static ferromagnetic order ($R>0$, $S=0$ for infinite $L$). The disorder–order transition is continuous, with susceptibility and Binder cumulant analysis indicating critical scaling in the 2D Ising universality class ($\nu \approx 1.01$, $\gamma \approx 1.78$, $\beta \approx 0.12$). The susceptibility peak scales as $\chi_{\max} \sim L^{\gamma/\nu}$ with $\gamma/\nu \approx 1.756\pm0.007 < 2$, confirming second-order character.

Within the static ordered phase, a smooth crossover (for moderate $\tilde{J}$) evolves into a first-order transition for $\tilde{J} \gtrsim 7.0$ ($\chi_{\max}\sim L^2$), terminating at a critical point reminiscent of the liquid–gas endpoint.

4. Monte Carlo Analysis in Three Dimensions

Simulations on cubic lattices ($L^3$) with identical update protocols confirm the existence of a stable swap (limit-cycle) phase at small $\Delta$ and intermediate $\tilde{J}$, with sustained oscillatory order ($S>0$) and no spiral defects. Upon increasing $\Delta$, the pathway for ordering follows swap $\rightarrow$ disorder $\rightarrow$ static order, in contrast to the mean-field scenario where a direct swap-to-static transition is allowed via SNIC or SNLC bifurcations. The expected critical exponents in 3D are anticipated to cross over toward mean-field XY values for Hopf-driven (I$\rightarrow$II) transitions, though they were not directly extracted.

5. Role of Single-Ion Anisotropy and Vacancies

The antisymmetric anisotropy $\Delta_A = -\Delta_B$ functions as a local chemical-potential imbalance, strongly biasing vacancy formation into one species (typically $A$). This local vacancy asymmetry deprives the "predator" species of on-site coupling targets, weakening nonreciprocal dynamical effects and suppressing swap oscillations. A sufficiently large $\Delta$ restores static, equilibrium-like order, even in the presence of strong nonreciprocal interactions and in low dimensions where global swapping would otherwise be destabilized by defects.

Control of vacancy energetics, such as via light-tunable adsorption in colloids or chemical potentials in cold-atom spinor systems, offers a practical experimental knob to stabilize or suppress dynamical, nonreciprocal phases (R et al., 21 Dec 2025).

6. Related Models: Directed Small-World Substrates

Directed small-world rewiring introduces nonreciprocal couplings by randomizing the directionality of links in the lattice (rewiring probability $q$). The Hamiltonian for the nonreciprocal spin-1 Blume-Capel model on such substrates reads: $$H_{\mathrm{directed}} = -\sum_{i,j} J_{ij} S_i S_j + \Delta \sum_i S_i^2, \quad J_{ij} = J A_{ij}$$ Transitions are classified by $q$:

• $q < q_c \approx 0.35$: Second-order, but with exponents that depart from regular 2D Blume-Capel values, indicating a change in universality class.
• $q > q_c$: First-order, with a discontinuity in the magnetization and Binder cumulant signature, indicative of latent-heat behavior.

The presence of a tricritical point at $q=q_c$ and continuously varying exponents for $q<q_c$ highlights the rich phase structure introduced by nonreciprocal disorder (Fernandes et al., 2010).

7. Significance and Outlook

The NR-BCM illustrates how competing local and nonlocal interactions—specifically, vacancy energetics and nonreciprocal couplings—shape collective dynamics, criticality, and order in spin systems. The model provides a solvable framework for exploring nonequilibrium phase transitions, defect-driven instability of nonreciprocal phases, and the restoration of equilibrium universality under suitable vacancy biasing. These findings establish vacancy energetics as a central mechanism for controlling nonreciprocal dynamical phases, extending the reach of equilibrium critical phenomena into fundamentally nonequilibrium regimes (R et al., 21 Dec 2025, Fernandes et al., 2010).