Nagle–Kardar Model: Mixed-Range Ising System

• The Nagle–Kardar model is a one-dimensional Ising system with nearest-neighbor and mean-field couplings that produces continuous, first-order, and tricritical phase transitions.
• It employs transfer-matrix analysis, Landau expansion, and Fokker–Planck reduction to elucidate equilibrium properties, critical exponents, and Model A dynamics.
• The model reveals ensemble inequivalence, metastable switching via Arrhenius scaling, and insights into quantum extensions, serving as a testbed for nonequilibrium coarse-graining.

Searching arXiv for papers on the Nagle–Kardar model and related usage of the term. The Nagle–Kardar model denotes a one-dimensional Ising-like spin system with competing short-range and long-range interactions. In its standard classical form, it combines a nearest-neighbor coupling and a mean-field coupling on a ring of spin-$\tfrac12$ variables, and it is notable for exhibiting a line of continuous phase transitions, a first-order transition line, and a tricritical point despite its one-dimensional geometry (Kemmeter et al., 10 Nov 2025). Recent work has also used the same name for a quantum extension with a transverse field (Campa et al., 28 Apr 2026). In a separate strand of the literature, the label “Nagle–Kardar” has occasionally been used for the continuous Hwa–Kardar model of anisotropic self-organized criticality, but that usage refers to a different field-theoretic object and should be distinguished from the Ising-like model discussed here (Antonov et al., 2015).

1. Definition and Hamiltonian

The classical Nagle–Kardar model is defined on a one-dimensional ring of $N$ Ising spins,

$s_i=\pm 1,\qquad i=1,\dots,N,$

with periodic boundary conditions $s_{N+1}=s_1$ (Kemmeter et al., 10 Nov 2025). Its Hamiltonian is

$$H = -J \sum_{i=1}^N s_i s_{i+1} - \frac{K}{2N} \left(\sum_{i=1}^N s_i\right)^2 - h \sum_{i=1}^N s_i, \tag{1}$$

where $J$ is the nearest-neighbor coupling, $K$ is the mean-field coupling, and $h$ is an external field (Kemmeter et al., 10 Nov 2025).

This form makes explicit the coexistence of two interaction ranges. The nearest-neighbor term is sensitive to local domain structure: $J>0$ favors ferromagnetic alignment, while $J<0$ favors antiferromagnetic order. The mean-field term couples every spin to every other spin equally and favors global magnetization when $N$0 (Kemmeter et al., 10 Nov 2025). The model is therefore a paradigmatic mixed-range Ising system. A useful reformulation introduces the total magnetization

$N$1

and the number of defects $N$2, where defects are bonds with $N$3,

$N$4

In terms of $N$5 and $N$6, the Hamiltonian becomes

$N$7

This representation is central for both equilibrium analysis and dynamical coarse-graining (Kemmeter et al., 10 Nov 2025, Kemmeter et al., 24 Jun 2026).

The same mixed-range structure is used in finite-size thermodynamic studies, often with couplings denoted $N$8 and $N$9, or in reduced form $s_i=\pm 1,\qquad i=1,\dots,N,$0 and $s_i=\pm 1,\qquad i=1,\dots,N,$1 (Dantchev et al., 2024, Dantchev et al., 21 Aug 2025). In that notation, the model is a one-dimensional Ising chain with nearest-neighbor interaction and equivalent-neighbor ferromagnetic interaction of strength $s_i=\pm 1,\qquad i=1,\dots,N,$2 under periodic boundary conditions (Dantchev et al., 2024).

2. Equilibrium structure and tricriticality

The static properties are exactly accessible through transfer-matrix and mean-field methods (Kemmeter et al., 10 Nov 2025). Using a Hubbard–Stratonovich transformation, one obtains an effective free-energy density whose minimum determines the equilibrium magnetization $s_i=\pm 1,\qquad i=1,\dots,N,$3. The corresponding self-consistency equation is

$s_i=\pm 1,\qquad i=1,\dots,N,$4

The equilibrium defect density $s_i=\pm 1,\qquad i=1,\dots,N,$5 is then determined by $s_i=\pm 1,\qquad i=1,\dots,N,$6 (Kemmeter et al., 10 Nov 2025).

At zero field, $s_i=\pm 1,\qquad i=1,\dots,N,$7, the phase diagram in the $s_i=\pm 1,\qquad i=1,\dots,N,$8 plane contains a continuous transition line,

$s_i=\pm 1,\qquad i=1,\dots,N,$9

a first-order transition line for sufficiently negative $s_{N+1}=s_1$0 and sufficiently large $s_{N+1}=s_1$1, and a tricritical point where the two meet,

$s_{N+1}=s_1$2

These results recur in later finite-size and dynamical studies (Kemmeter et al., 10 Nov 2025, Kemmeter et al., 24 Jun 2026).

A Landau expansion of the free energy around $s_{N+1}=s_1$3,

$s_{N+1}=s_1$4

shows that along the second-order line $s_{N+1}=s_1$5 but $s_{N+1}=s_1$6, whereas at the tricritical point both $s_{N+1}=s_1$7 and $s_{N+1}=s_1$8, leaving $s_{N+1}=s_1$9 as the leading stabilizing term (Kemmeter et al., 10 Nov 2025). This yields mean-field Ising static exponents on the critical line,

$$H = -J \sum_{i=1}^N s_i s_{i+1} - \frac{K}{2N} \left(\sum_{i=1}^N s_i\right)^2 - h \sum_{i=1}^N s_i, \tag{1}$$0

and tricritical exponents

$$H = -J \sum_{i=1}^N s_i s_{i+1} - \frac{K}{2N} \left(\sum_{i=1}^N s_i\right)^2 - h \sum_{i=1}^N s_i, \tag{1}$$1

The change in upper critical dimension from $$H = -J \sum_{i=1}^N s_i s_{i+1} - \frac{K}{2N} \left(\sum_{i=1}^N s_i\right)^2 - h \sum_{i=1}^N s_i, \tag{1}$$2 to $$H = -J \sum_{i=1}^N s_i s_{i+1} - \frac{K}{2N} \left(\sum_{i=1}^N s_i\right)^2 - h \sum_{i=1}^N s_i, \tag{1}$$3 reflects the shift from an effective $$H = -J \sum_{i=1}^N s_i s_{i+1} - \frac{K}{2N} \left(\sum_{i=1}^N s_i\right)^2 - h \sum_{i=1}^N s_i, \tag{1}$$4 to an effective $$H = -J \sum_{i=1}^N s_i s_{i+1} - \frac{K}{2N} \left(\sum_{i=1}^N s_i\right)^2 - h \sum_{i=1}^N s_i, \tag{1}$$5 description at tricriticality (Kemmeter et al., 10 Nov 2025).

A notable feature is that the one-dimensional geometry does not preclude bulk criticality because the infinite-range term changes the effective thermodynamics. This suggests that the model should not be interpreted as a conventional short-range 1D Ising chain perturbed weakly by finite-size corrections; rather, the mean-field term is structurally dominant in establishing long-range order (Kemmeter et al., 10 Nov 2025).

3. Fluctuations, large deviations, and finite-size behavior

The Nagle–Kardar model admits an explicit large-deviation description of magnetization fluctuations. The equilibrium magnetization distribution satisfies

$$H = -J \sum_{i=1}^N s_i s_{i+1} - \frac{K}{2N} \left(\sum_{i=1}^N s_i\right)^2 - h \sum_{i=1}^N s_i, \tag{1}$$6

with $$H = -J \sum_{i=1}^N s_i s_{i+1} - \frac{K}{2N} \left(\sum_{i=1}^N s_i\right)^2 - h \sum_{i=1}^N s_i, \tag{1}$$7 a quasi-potential whose minimum determines the equilibrium magnetization (Kemmeter et al., 10 Nov 2025). Near $$H = -J \sum_{i=1}^N s_i s_{i+1} - \frac{K}{2N} \left(\sum_{i=1}^N s_i\right)^2 - h \sum_{i=1}^N s_i, \tag{1}$$8, the form of $$H = -J \sum_{i=1}^N s_i s_{i+1} - \frac{K}{2N} \left(\sum_{i=1}^N s_i\right)^2 - h \sum_{i=1}^N s_i, \tag{1}$$9 changes across the phase diagram. Away from criticality, $J$0; on the critical line, $J$1; and at the tricritical point, $J$2 (Kemmeter et al., 10 Nov 2025).

These forms directly determine finite-size scaling of equilibrium fluctuations: $J$3 away from criticality,

$J$4

on the critical line, and

$J$5

at the tricritical point (Kemmeter et al., 10 Nov 2025). Equivalently, with $J$6,

$J$7

so that $J$8 on the critical line and $J$9 at tricriticality (Kemmeter et al., 10 Nov 2025).

The joint distribution of magnetization and defect density also satisfies a large-deviation principle,

$K$0

with explicit quasi-potential

$K$1

Minima of $K$2 define locally stable macrostates, while saddle points define activation barriers between them (Kemmeter et al., 10 Nov 2025, Kemmeter et al., 24 Jun 2026).

This two-variable structure is one of the model’s distinctive features. Magnetization $K$3 is the thermodynamic order parameter, but defect density $K$4 carries the domain-wall information required to track local rearrangements and activated transitions (Kemmeter et al., 24 Jun 2026).

4. Glauber dynamics, Fokker–Planck reduction, and Model A criticality

The canonical dynamics most commonly studied for the Nagle–Kardar model is single-spin-flip Glauber dynamics with non-conserved magnetization (Kemmeter et al., 10 Nov 2025). At each update, one site is selected, a spin flip $K$5 is proposed, the associated energy change $K$6 is computed from the Hamiltonian, and the move is accepted with probability

$K$7

The microscopic dynamics obeys detailed balance with respect to the canonical Gibbs distribution (Kemmeter et al., 10 Nov 2025).

Because the Hamiltonian depends only on magnetization $K$8 and defects $K$9, the full Markov process can be coarse-grained to a master equation for the macrostate probability $h$0. In the large-$h$1 limit, one introduces

$h$2

and obtains a two-dimensional Fokker–Planck equation,

$h$3

The drift terms $h$4 govern deterministic relaxation of magnetization and defect density, while $h$5 describe finite-size fluctuations (Kemmeter et al., 10 Nov 2025, Kemmeter et al., 24 Jun 2026).

The associated Langevin equation is

$h$6

with $h$7 and $h$8 equal to the diffusion matrix derived from the Fokker–Planck coefficients (Kemmeter et al., 10 Nov 2025). This mesoscopic description is sufficient to recover critical slowing down, metastable switching, and first-passage asymptotics.

Near criticality, center manifold theory reduces the deterministic dynamics to an effective one-dimensional evolution for the slow magnetization mode $h$9. Along the second-order line,

$J>0$0

which implies

$J>0$1

whereas at the tricritical point the cubic coefficient vanishes and the leading term is quintic,

$J>0$2

implying

$J>0$3

These decay laws yield the dynamic decay exponents

$J>0$4

on the critical line and

$J>0$5

at the tricritical point (Kemmeter et al., 10 Nov 2025).

Using the dynamic scaling relation

$J>0$6

together with the static exponents, one finds in both cases

$J>0$7

This establishes that the critical dynamics of the Nagle–Kardar model belongs to the Model A universality class in the Hohenberg–Halperin classification: purely relaxational, Markovian dynamics with a non-conserved scalar order parameter (Kemmeter et al., 10 Nov 2025, Kemmeter et al., 24 Jun 2026).

A common misconception is that the mean-field interaction should force mean-field dynamic behavior beyond Model A. The dynamical analyses instead show that, although the static critical structure is modified by the long-range term and tricriticality appears, the dynamic critical exponent remains $J>0$8 and the macroscopic relaxation mechanism is still that of non-conserved order-parameter dynamics (Kemmeter et al., 10 Nov 2025).

5. Metastability, first-passage times, and Arrhenius behavior

Because the quasi-potential $J>0$9 may have multiple minima, the Nagle–Kardar model exhibits bistability and metastability in the ordered region (Kemmeter et al., 10 Nov 2025). For $J<0$0, there are typically two symmetric minima at $J<0$1; for $J<0$2, symmetry is broken and one minimum becomes metastable (Kemmeter et al., 10 Nov 2025, Kemmeter et al., 24 Jun 2026).

The average first-passage time $J<0$3 for transitions between locally stable macrostates obeys an Arrhenius-type large-deviation law,

$J<0$4

where $J<0$5 is the barrier height between the initial minimum and the relevant saddle point in the quasi-potential landscape (Kemmeter et al., 10 Nov 2025). This is the natural finite-system analogue of activated barrier crossing, but with system size $J<0$6 multiplying the barrier because $J<0$7 is a large-deviation rate function.

Microscopic Glauber simulations and mesoscopic Langevin simulations both reproduce this scaling (Kemmeter et al., 24 Jun 2026). For symmetric bistable cases, left-to-right and right-to-left passage times coincide; with an external field they differ because the relevant barriers differ (Kemmeter et al., 10 Nov 2025). This direct correspondence between microscopic dynamics, macroscopic quasi-potential, and activated switching is one of the reasons the model is considered a useful testbed for nonequilibrium coarse-graining (Kemmeter et al., 24 Jun 2026).

A plausible implication is that the Nagle–Kardar model occupies an intermediate position between exactly solvable equilibrium mixed-range spin systems and stochastic metastability models: it is simple enough to allow explicit quasi-potentials and coarse-grained drifts, but rich enough to display tricritical scaling and activated phase switching (Kemmeter et al., 24 Jun 2026).

6. Finite-size fluctuation-induced forces and ensemble dependence

Finite systems of the Nagle–Kardar model display strong ensemble dependence because the long-range interaction makes the model non-additive (Dantchev et al., 21 Aug 2025). In the grand-canonical ensemble with periodic boundary conditions, the finite-size excess free energy defines a Casimir force,

$J<0$8

while in the canonical ensemble with fixed magnetization, the corresponding excess Helmholtz free energy defines a Helmholtz force,

$J<0$9

Here $N$00 for the one-dimensional ring (Dantchev et al., 21 Aug 2025).

Exact asymptotics show that on the critical line the critical Casimir amplitude is

$N$01

while at the tricritical point it is

$N$02

In both cases the corresponding critical Casimir force is repulsive near the singular regime (Dantchev et al., 2024, Dantchev et al., 21 Aug 2025). This is unusual because periodic boundary conditions are normally associated with attractive critical Casimir forces. The Nagle–Kardar model therefore violates the commonly cited “boundary condition rule” according to which equivalent boundary conditions imply attraction and conflicting ones imply repulsion (Dantchev et al., 2024, Dantchev et al., 21 Aug 2025).

The Helmholtz force behaves differently. In the canonical ensemble with fixed magnetization, the contribution of the long-range term cancels exactly from the excess free energy, so the Helmholtz force is independent of $N$03 and coincides with that of the standard one-dimensional Ising model at fixed magnetization (Dantchev et al., 21 Aug 2025). It nevertheless changes sign as a function of temperature and magnetization (Dantchev et al., 21 Aug 2025).

This ensemble dependence is not a minor finite-size correction. It directly reflects the non-additivity of the long-range interaction and the inequivalence of thermodynamic potentials in finite systems (Dantchev et al., 21 Aug 2025). The Nagle–Kardar model has therefore become a controlled setting for comparing fluctuation-induced forces across conjugate ensembles.

7. Quantum extension and related nomenclature

A quantum extension of the Nagle–Kardar model adds a transverse field and is defined by

$N$04

with $N$05 Pauli matrices (Campa et al., 28 Apr 2026). This quantum Nagle–Kardar model retains the competition between nearest-neighbor and mean-field interactions while introducing quantum fluctuations. The thermodynamic analysis shows that ensemble inequivalence present in the classical model is removed above a threshold transverse field $N$06: for $N$07, only second-order transition lines remain and ensemble equivalence is restored (Campa et al., 28 Apr 2026). This suggests that quantum fluctuations can smooth out the first-order structure responsible for inequivalence in the classical model.

A separate terminological issue concerns the phrase “Nagle–Kardar” in the context of anisotropic stochastic growth and self-organized criticality. Some papers describe the continuous Hwa–Kardar model as “Hwa–Kardar (Nagle–Kardar)” (Antonov et al., 2015). That model is a stochastic height equation with anisotropic diffusion, nonlinear transport along a preferred direction, and white noise, used for self-organized critical systems such as driven sandpiles (Antonov et al., 2015). It is not the same object as the Ising-like mixed-range spin model governed by Hamiltonian (1). The overlap in naming is therefore historical or contextual rather than structural.

This distinction matters because the two model families belong to different research programs. The Ising-like Nagle–Kardar model concerns competing short- and long-range interactions, tricriticality, metastability, ensemble inequivalence, and Glauber dynamics (Kemmeter et al., 10 Nov 2025, Kemmeter et al., 24 Jun 2026, Dantchev et al., 21 Aug 2025). The Hwa–Kardar model concerns anisotropic stochastic growth, self-organized criticality, and renormalization-group analysis under turbulent advection (Antonov et al., 2015). Conflating the two obscures both.

The Nagle–Kardar model, in its classical and quantum forms, is therefore best understood as an exactly tractable mixed-range Ising system whose combination of local domain-wall structure and global ferromagnetic coupling produces tricritical equilibrium behavior, Model A critical dynamics, activated metastability, and pronounced ensemble dependence (Kemmeter et al., 10 Nov 2025, Kemmeter et al., 24 Jun 2026, Dantchev et al., 21 Aug 2025, Campa et al., 28 Apr 2026).