Ising Supercritical Regime

• ISR is a regime characterized by smooth crossover ridges replacing the traditional phase boundary above the CEP, defined by peaks in specific heat and susceptibility.
• Universal scaling laws with critical exponents (β ≈ 0.326, γ ≈ 1.237) govern ISR, revealing an emergent Z2 symmetry across magnetic, quantum, and fluid systems.
• ISR leads to nontrivial phenomena such as divergent magnetic Grüneisen ratios and singular magnetocaloric effects, enabling advanced cooling protocols and deeper phase transition insights.

The Ising Supercritical Regime (ISR) generalizes the critical crossover structure above the critical endpoint (CEP) in systems exhibiting Ising universality, unifying supercriticality in magnetic, quantum, and fluid systems. In ISR, incipient phase boundaries above the CEP are replaced by symmetry-protected thermodynamic ridges—crossovers defined by maxima in susceptibilities and response functions—with universal scaling and emergent $Z_2$ (Ising) symmetry. These features yield highly nontrivial phenomena such as divergent magnetic Grüneisen ratios, singular magnetocaloric effects, and robust corrections induced by field-mixing asymmetry. The ISR has been confirmed in a wide range of contexts, from frustrated antiferromagnets (e.g., Nd$_3$BWO$_9$) to classical fluids, quantum spin chains, and dynamical nonequilibrium settings.

1. Emergence of ISR and Definition of Supercritical Crossovers

In systems of the Ising universality class, the traditional phase boundary between ordered and disordered phases terminates at a finite-temperature CEP (e.g., $(\mu_0 H_c, T_c)$ in Nd$_3$BWO$_9$ (Liu et al., 12 Jan 2026)). Beyond this point, the first-order transition evolves into a pair of smooth crossover ridges—the ISR—which in field–temperature (or pressure–temperature) phase space demarcate regimes of distinct fluctuational character. In the Ising context, these lines are loci of maxima of specific heat and susceptibility:

$$\text{ISR boundaries:}\quad C(T, H) = \max, \quad \chi(T, H) = \max,$$

which, in reduced variables $t\equiv(T-T_c)/T_c$, $h\equiv(H-H_c)/H_c$, obey the universal scaling

$h \propto t^{\beta+\gamma}$

with critical exponents $_3$0 of the 3D Ising class, yielding $_3$1 (Li et al., 2023, Liu et al., 12 Jan 2026, Cui et al., 21 Sep 2025). This structure is mirrored in classical fluids as the two Widom-type crossover lines $_3$2 and $_3$3 in the supercritical region, separating liquid-like, indistinguishable, and gas-like regimes with explicit $_3$4 symmetry about the critical isochore (Li et al., 2023, Cui et al., 21 Sep 2025, Li et al., 9 Dec 2025).

2. Universal Scaling Structure and Symmetry in ISR

The singular part of the free energy in the ISR region above the CEP takes the scaling form

$_3$5

where $_3$6 is the specific heat exponent. From this, all thermodynamic quantities inherit their scaling:

• Specific heat: $_3$7
• Susceptibility: $_3$8

Crossover ridges are located by maximizing $_3$9 or $_9$0 with respect to $_9$1 or $_9$2, yielding the ridge loci $_9$3 (Liu et al., 12 Jan 2026).

Importantly, the Ising $_9$4 symmetry enforces the presence of two symmetric crossovers—a consequence of the oddness of the equation of state and the self-reciprocal property between subcritical coexistence branches and supercritical ridges. This dual-ridge structure has been rigorously established in both mean-field (van der Waals) and Ising models, as well as verified empirically in supercritical fluids (Ar, CO$_9$5) (Li et al., 2023, Cui et al., 21 Sep 2025).

Supercritical Observable	Scaling at ISR	Physical Meaning
$_9$6, $_9$7 Peak	$_9$8	Crossover ridge location
Magnetocaloric Ratio $_9$9	$(\mu_0 H_c, T_c)$0	Divergence at CEP

3. Quantum and Dynamical Generalizations of ISR

In quantum Ising systems (e.g., TFIM), ISR emerges as a “quantum supercritical regime” (QSR), controlled by symmetry-breaking longitudinal field $(\mu_0 H_c, T_c)$1 in the $(\mu_0 H_c, T_c)$2–$(\mu_0 H_c, T_c)$3 plane. Crossover boundaries are given by

$(\mu_0 H_c, T_c)$4

where $(\mu_0 H_c, T_c)$5 is the dynamical exponent and $(\mu_0 H_c, T_c)$6 the correlation length exponent (Lv et al., 2024, Wang et al., 2024). The singular part of the free energy yields enhanced divergences in, e.g., the Grüneisen ratio:

$(\mu_0 H_c, T_c)$7

This supercritical region is accompanied by robust data collapse in scaling variables and distinctive dynamical phenomena including nonanalytic cusps in the Loschmidt rate function at quenches to the QSC boundary. The boundary is determined by the intersection of Lee-Yang-Fisher zeros with the real-time axis (Wang et al., 2024, Ouyang et al., 2023).

4. ISR in Fluids: Crossover Loci and Emergent Symmetry

In classical fluids, ISR is manifested via two lines $(\mu_0 H_c, T_c)$8 and $(\mu_0 H_c, T_c)$9 in the ($_3$0) or ($_3$1) supercritical plane, defined by maxima of isothermal compressibility along isochoric paths parallel to $_3$2. These lines obey the same Ising scaling:

$_3$3

with $_3$4. The emergent $_3$5 symmetry in the scaling domain implies these branches are related by a simple inversion and their physical characterization (e.g., liquid-like/indistinguishable/gas-like) precisely mirrors the Ising magnetic case (Cui et al., 21 Sep 2025, Li et al., 9 Dec 2025).

Experimental validation comes from correlation between dynamic (inelastic X-ray/neutron scattering) and thermodynamic markers of the crossovers (compressibility, $_3$6 peaks), confirming the predicted exponents and symmetry structure (Li et al., 2023).

5. Asymmetric Corrections and Antisymmetric Scaling in ISR

Deviations from ideal Ising symmetry, especially in real fluids, are encapsulated by linear mixing of field variables (chemical potential, pressure, temperature) into the scaling fields. This produces universal antisymmetric corrections to the scaling of ISR boundaries and their “diameters” (mean lines between $_3$7 and $_3$8):

$_3$9

$_9$0

Higher-order cumulants (e.g., $_9$1, $_9$2 of the order parameter distribution) exhibit analogous scaling structure, with the “symmetry line” $_9$3 tracking the leading antisymmetric corrections and $_9$4 maxima coalescing near the ISR crossover lines (Li et al., 9 Dec 2025). Such corrections are tightly constrained by the field-mixing parameters and can be used to quantify the degree of asymmetry in real systems.

ISR Correction	Scaling Exponent	Parameter Dependence
Diameter in $_9$5	$_9$6	$_9$7, $_9$8, $_9$9
Diameter in $$\text{ISR boundaries:}\quad C(T, H) = \max, \quad \chi(T, H) = \max,$$0	$$\text{ISR boundaries:}\quad C(T, H) = \max, \quad \chi(T, H) = \max,$$1, $$\text{ISR boundaries:}\quad C(T, H) = \max, \quad \chi(T, H) = \max,$$2	$$\text{ISR boundaries:}\quad C(T, H) = \max, \quad \chi(T, H) = \max,$$3, $$\text{ISR boundaries:}\quad C(T, H) = \max, \quad \chi(T, H) = \max,$$4, $$\text{ISR boundaries:}\quad C(T, H) = \max, \quad \chi(T, H) = \max,$$5

6. ISR in Dynamical and Driven Ising Models

In nonequilibrium kinetic Ising models subjected to time-dependent fields, the supercritical regime is defined for $$\text{ISR boundaries:}\quad C(T, H) = \max, \quad \chi(T, H) = \max,$$6, period $$\text{ISR boundaries:}\quad C(T, H) = \max, \quad \chi(T, H) = \max,$$7 (slow driving). Under a nonzero bias field $$\text{ISR boundaries:}\quad C(T, H) = \max, \quad \chi(T, H) = \max,$$8, the system displays novel “sideband” peaks in the fluctuation spectrum and dynamic susceptibility, linked to resonant delay in magnetization switching governed by droplet nucleation dynamics. The scaling of the sideband positions and their dependence on system size, period, and effective interface tension have been quantitatively confirmed, both in the two-state Ising model and in the Blume–Capel generalization (where the presence of intermediate spin states reduces interface tension and modulates sideband structure) (Buendía et al., 11 Apr 2025).

7. Physical Consequences and Technological Implications

ISR has direct implications for magnetocaloric cooling and quantum thermodynamics. In Kagome antiferromagnets such as Nd$$\text{ISR boundaries:}\quad C(T, H) = \max, \quad \chi(T, H) = \max,$$9BWO$t\equiv(T-T_c)/T_c$0, the universality-protected divergence of the Grüneisen ratio within ISR yields record cooling rates, high entropy extraction, and sub-Kelvin refrigeration capabilities exceeding traditional paramagnetic salts and ferromagnets. The universality and robustness of ISR scaling point to broad utility in rare-earth–based Ising magnets and potentially quantum refrigerants where traditional $t\equiv(T-T_c)/T_c$1He-based systems are impractical (Liu et al., 12 Jan 2026, Lv et al., 2024).

ISR analysis enables unified understanding and predictive power across classical and quantum phase transitions, guiding experimental searches for new phases, critical cooling protocols, and advanced quantum simulation platforms. Emergent $t\equiv(T-T_c)/T_c$2 symmetry and scaling extend across diverse systems, resolving long-standing questions on asymmetry, singular corrections, and the boundaries of universality in critical phenomena.