---
title: Direction-Selective Criticality
url: https://www.emergentmind.com/topics/direction-selective-criticality
type: topic
---

# Direction-Selective Criticality

Direction-selective criticality denotes a class of critical phenomena in which instability, softening, or divergent response is not global across all relevant degrees of freedom, but is selected along particular directions, modes, symmetry sectors, or propagation channels. In recent arXiv literature, the term appears in several technically distinct settings: mean-field rotor Hamiltonians, where criticality is identified with the vanishing of curvature coefficients in a finite-dimensional collective sector of the microcanonical energy shell; spin-orbit-coupled magnets, where only one symmetry-inequivalent excitation branch becomes gapless; nemato-elastic systems, where Saint Venant compatibility suppresses incompatible nematic fluctuations and leaves only a compatible critical subspace; and directed networks and hypergraphs, where non-reciprocity, anchors, or feed-forward propagation reshape thresholds, scaling, and selected steady states [2603.29074, 2604.20173, 2507.23754, 2601.20726].

## 1. Core meaning and recurring definitions

Across these works, the “direction” relevant to criticality is not a single universal object. In the geometric rotor formulation, it is a unit vector in order-parameter space, probing perturbations $\delta m = u v$ of collective amplitudes $m = (m_1,\ldots,m_r)^{\mathsf T}$ [2603.29074]. In BNZS, it is a symmetry-distinct excitation sector associated with one of two Kramers-doublet-derived gaps, $\Delta_L$ and $\Delta_H$, under a field applied along a crystallographic direction [2604.20173]. In nemato-elasticity, it is a momentum-space direction $\hat q$ that determines whether an orbital nematic fluctuation projects onto a compatible or incompatible helical subspace [2507.23754]. In directed hypergraphs and non-reciprocal networks, it is embedded in asymmetric propagation itself: forward versus backward percolative channels, or the degree of reciprocity $\gamma$ in a non-Hermitian interaction matrix [2601.20726, 2312.12039].

A common contrast is with global criticality. BNZS explicitly distinguishes mode-selective or partial quantum criticality from “conventional, global quantum criticality,” where the entire low-energy spectrum softens collectively [2604.20173]. The rotor framework similarly replaces a thermodynamic-only characterization by a spectral geometric criterion that selects which collective mode loses quadratic rigidity first [2603.29074]. Nemato-elastic theory makes the same point in another language: only fluctuations lying in a compatible doublet become critical, while other symmetry components remain massive even at the transition [2507.23753]. These formulations differ in microscopic content, but they converge on a restricted-softening picture.

This suggests a broad conceptual distinction between two kinds of critical organization. In one, criticality is extensive across the low-energy sector. In the other, the system approaches instability through a constrained eigenspace, with the noncritical sector remaining stiff, gapped, ordered, or only weakly affected.

## 2. Geometric selection in mean-field rotor Hamiltonians

In finite-dimensional trigonometric mean-field rotor Hamiltonians, direction-selective criticality is formulated as a property of the extrinsic geometry of the microcanonical constant-energy shell
$$
\Sigma_E = \{x \in \Lambda \mid H(x)=E\},
$$
with $\Lambda \simeq \mathbb{R}^{2N}$. The key local observable is the trace of the Weingarten operator,
$$
W = \Pi [\nabla^2 H(x)] \Pi / \|\nabla H(x)\|,
$$
whose eigenvalues are the principal curvatures. In Euclidean gauge,
$$
\operatorname{Tr} W = \operatorname{div} n
= \frac{\Delta H}{\|\nabla H\|} - \frac{\nabla H^{\mathsf T} (\nabla^2 H)\nabla H}{\|\nabla H\|^3},
$$
and in the unit-normal gauge it directly controls microcanonical entropy derivatives through
$$
\partial_E S(E)=\langle \operatorname{Tr} W\rangle_E,
$$
with $S(E)=\ln \operatorname{area}(\Sigma_E)$ [2603.29074].

For Hamiltonians of the form
$$
H(\theta,p)=\sum_{i=1}^N \frac{p_i^2}{2}+Nv_0-\frac{N}{2}\sum_{a,b=1}^r J_{ab} m_a m_b,
$$
with
$$
m_a=\frac{1}{N}\sum_{i=1}^N \phi_a(\theta_i),
$$
the mean curvature per particle admits the universal collective expansion
$$
\operatorname{Tr}W(E,m)/N \simeq \frac{1}{c(E)} + m^{\mathsf T} C(E) m + O(\|m\|^3)+O(N^{-1}),
$$
where
$$
C(E)=\frac{1}{c(E)^2}\,[c(E)M-B],
$$
$$
M=\frac{1}{2}(JD+D^{\mathsf T}J), \qquad B=JQ_*J.
$$
The matrices $D$ and $Q_*$ encode closure of the trigonometric family and branch covariance on the reference branch, while $J$ contains the finite set of collective couplings [2603.29074].

The selection principle is spectral. If
$$
C(E) v_\alpha(E)=\lambda_\alpha(E) v_\alpha(E),
$$
then criticality occurs when the smallest curvature eigenvalue vanishes,
$$
\lambda_{\min}(E_c)=0,
$$
equivalently $\det C(E_c)=0$. The associated eigenvector $v_{\min}(E_c)$ is the critical channel. In this formulation, the energy shell loses quadratic geometric rigidity first along a distinguished collective direction, and phase transition onset is reinterpreted as a geometric instability intrinsic to $\Sigma_E$ [2603.29074].

The worked examples make the criterion explicit. In the Hamiltonian Mean-Field model,
$$
H=\sum_{i=1}^N \frac{p_i^2}{2}+\frac{N}{2}-\frac{N}{2}(m_x^2+m_y^2),
$$
with $v_0=1/2$, $J=I_2$, and $c(\epsilon)=2\epsilon-1$, the collective curvature form is
$$
C(\epsilon)=\frac{1}{c(\epsilon)^2}\left[c(\epsilon)I_2-\frac{1}{2}I_2\right],
$$
so the critical condition $c(\epsilon_c)=1/2$ yields $\epsilon_c=3/4$. In multimode diagonal sectors one has
$$
C_\alpha(E)=\frac{k_\alpha^2 J_\alpha}{c(E)^2}\left[c(E)-\frac{J_\alpha}{2}\right],
$$
hence $c(E_{c,\alpha})=J_\alpha/2$ and, on the disordered branch, $\epsilon_{c,\alpha}=v_0+J_\alpha/4$ [2603.29074].

A central limitation is explicit in the same framework: the geometric criterion identifies the mechanism and the critical energy, but not the order of the transition. First- versus second-order behavior still requires additional analysis of $S(E)$ and its derivatives, including microcanonical inflection-point analysis.

## 3. Sector-selective quantum and elastic criticality in condensed matter

In BaNd$_2$ZnS$_5$, direction-selective criticality takes the form of a field-induced, mode-selective quantum phase transition. Below the Néel temperature $T_N=2.9\,\mathrm{K}$, the Nd$^{3+}$ Kramers doublets produce two symmetry-inequivalent low-energy excitation sectors with gaps $\Delta_L$ and $\Delta_H$. For $H\parallel[110]$, the lower gap $\Delta_L$ softens continuously and collapses at $H_c \simeq 2\,\mathrm{T}$, while the higher gap $\Delta_H$ remains finite; neutron diffraction and thermodynamics indicate that the $q_2$ component collapses near $H_c$, whereas $q_1$ long-range order persists until much higher fields, near $12$–$12.7\,\mathrm{T}$. The intermediate partially critical phase appears only for $H\parallel[110]$; it is absent for $H\parallel[100]$ and $H\parallel[001]$, where $q_1$ and $q_2$ remain symmetry-equivalent under the field [2604.20173].

The thermodynamic signatures are correspondingly partial rather than global. At criticality, the ac susceptibility follows $\chi_{ac}\sim T^{-0.2}$ and collapses according to
$$
\chi_{ac}(T,H)=T^{-\gamma} f((H-H_c)/T^{1/(\nu z)}),
$$
with $\gamma=0.67$ and $\nu z=2.1$. The residual Sommerfeld coefficient $\gamma_0$ increases markedly and shows an apparent divergence near $H_c \approx 2\,\mathrm{T}$, indicating a dense set of gapless excitations confined to the critical symmetry sector. The transition is described as continuous within experimental resolution: Ehrenfest consistency holds for the heat-capacity anomalies, whereas Clausius-Clapeyron fails [2604.20173].

In nemato-elastic systems, the selection mechanism is instead enforced by compatibility. The strain tensor $\epsilon_{ij}=(\partial_i u_j+\partial_j u_i)/2$ must satisfy the Saint Venant relations,
$$
(\mathrm{Curl}\,\mathrm{Curl}\,\epsilon)_{ij}
\equiv \epsilon_{ikm}\epsilon_{jln}\partial_k\partial_l \epsilon_{mn}=0,
$$
or equivalently
$$
\epsilon_{ij,kl}+\epsilon_{kl,ij}-\epsilon_{ik,jl}-\epsilon_{jl,ik}=0.
$$
In a co-rotating helical basis, these constraints split the five-component traceless nematic fluctuation space into a compatible doublet and noncritical amplitudes. The effective Gaussian kernel can be written as
$$
S_{\mathrm{eff}}=\frac{1}{2V}\sum_q \Phi^\dagger(q)\left[(r+cq^2)I-\frac{g^2}{\mu}(I-\Gamma(\hat u))\right]\Phi(q),
$$
with masses
$$
r_{2,3}(q)=r+cq^2-\frac{g^2}{\mu}, \qquad
r_1(q)=r+cq^2-\frac{g^2}{\mu}\rho, \qquad
r_{4,5}(q)=r+cq^2.
$$
Thus only $\Phi_2$ and $\Phi_3$ become critical at $r_c=g^2/\mu$, while incompatible fluctuations remain gapped [2507.23754].

The directional selectivity follows from projection back to orbital channels. For an Ising nematic $Q_{2xy}$, the critical directions are $q\parallel[100]$ and $[010]$; for $Q_{x^2-y^2}$, they are $q\parallel[110]$ and $[1\bar{1}0]$. A companion formulation describes the same bifurcation as “compatible instability,” emphasizing that the critical modes are protected from pinning by defect strains, while defects generate long-ranged random longitudinal and transverse conjugate fields only in noncritical helical channels [2507.23753]. One consequence is that mean-field thermodynamics and widespread domain formation are not contradictory within this framework.

## 4. Directed propagation, non-reciprocity, and higher-order network criticality

A distinct use of direction-selective criticality appears in nonequilibrium directed systems. In a chain of adaptive excitable integrators, directionality is literal: coupling is strictly feed-forward, $i\to i+1$, the drive enters only at $i=1$, and the boundary at $i=N$ is open. The adaptive thresholds obey
$$
\Theta_i^{(t+1)}=\Theta_i^{(t)}(1-w)+|\tilde V_i^{(t)}|\,w,
$$
with $0<w\ll 1$, and the per-level gain
$$
g_i \equiv P(|\tilde V_i|>\Theta_i)
$$
self-organizes near marginal propagation, numerically $g_i\approx 0.43$–$0.47$ and $g_*\approx 0.4565$ in the fixed-threshold approximation. The model exhibits discrete scale invariance,
$$
P(S)\propto S^{-\tau}F\!\left(\frac{\ln S}{\ln \lambda}\right), \qquad \lambda\approx k_*\approx 1.48,
$$
with mixture and sum exponents near $-3$ and $-2$, respectively. Threshold and subthreshold spectra are Lorentzian with position-dependent corner frequencies, decreasing along the chain from $g_{0,1}\approx 0.318\,\mathrm{Hz}$ to $g_{0,7}\approx 0.003\,\mathrm{Hz}$ and from $f_{0,1}\approx 59.2\,\mathrm{Hz}$ to $f_{0,7}\approx 0.38\,\mathrm{Hz}$ [2201.07075].

In non-reciprocal neural networks, directionality is encoded by the reciprocity parameter $\gamma\in[-1,1]$ in the coupling statistics,
$$
\overline{J_{ij}J_{ji}-\frac{J_0^2}{N^2}}=\gamma\,\frac{J^2}{N}.
$$
Linear stability of the quiescent state is governed by
$$
\frac{1}{gJ}
=
\max\!\left(
1+\gamma,\,
\frac{J_0}{J}+\gamma\,\frac{J}{J_0}
\right).
$$
The system supports paramagnetic, ferromagnetic, and spin-glass-like regions. In the spin-glass region, reciprocal couplings produce marginal behavior, whereas decreasing reciprocity drives a smooth transition to chaos; the spin-glass region shrinks and disappears at $\gamma=-1$. Dynamic mean-field theory identifies a selected separatrix state through
$$
V(C_0^\star\,|\,C_0^\star,0)=0
$$
in the noiseless case and
$$
V(C_\sigma^\star\,|\,C_\sigma^\star,0)=-\frac{1}{2}\sigma^4
$$
with noise. In the ferromagnetic region, only fixed points are dynamically realizable; in the spin-glass region, the selected state is marginal in the ensemble description but single realizations generically display chaos, with $\Lambda>0$ for sufficiently large $N$ [2312.12039].

Directed hypergraph percolation generalizes the same logic to higher-order interactions with asymmetric functional dependencies. A directed hyperedge maps an input set $N^{(-)}(\alpha)$ to an output set $N^{(+)}(\alpha)$, while anchor nodes encode indispensable participants. The renormalized availability is
$$
\pi_N = 1-\theta+\theta p_N,
$$
and the percolation threshold is controlled by
$$
\hat\Lambda
=
p_N
\frac{\langle q^{\rm in}q^{\rm out}\rangle}{\langle q^{\rm out}\rangle}
\frac{\langle p_H^{[\mathbf m]}\,\pi_N^{m-2}\,m^{\rm out}m^{\rm in}\rangle}{\langle m^{\rm out}\rangle}
=1.
$$
The Hypergraph Giant In Component, Hypergraph Giant Out Component, and Hypergraph Giant Strongly Connected Component emerge simultaneously, but post-critical scaling differs. In finite-moment regimes,
$$
\beta^{(+)}=\beta^{(-)}=1, \qquad \beta=2.
$$
In maximally correlated heavy-tailed regimes, anomalous exponents depend on whether node or hyperedge percolation is considered, and anchor-free systems can exhibit modified composition rules such as $\beta_R=\beta_R^{(+)}+\beta_R^{(-)}-1$ in the reported $\gamma<3$ regimes [2601.20726].

Taken together, these systems show that “direction-selective” need not refer only to spatial anisotropy. It can refer to asymmetric causal structure, forward/backward reachability, or non-reciprocal spectral selection of attractors and exponents.

## 5. Strain direction and optimization trade-offs as selectors of instability channels

In hole-doped manganite La$_{0.75}$Ca$_{0.25}$MnO$_3$, direction-selective criticality is implemented by uniaxial strain as a crystallographically resolved tuning field. The structural response is decomposed into breathing and Jahn-Teller modes,
$$
Q_{1,i}=\frac{X_i+Y_i+Z_i}{\sqrt 3}, \qquad
Q_{2,i}=\frac{X_i-Y_i}{\sqrt 2}, \qquad
Q_{3,i}=\frac{2Z_i-X_i-Y_i}{\sqrt 6},
$$
with $\rho_i\equiv \sqrt{Q_{2,i}^2+Q_{3,i}^2}$, alongside site-average and site-selective combinations such as $Q_{1,\mathrm{avg}}$, $Q_{1,\mathrm{diff}}$, $Q_{\rho,\mathrm{avg}}$, and $Q_{\rho,\mathrm{diff}}$. Extreme uniaxial strain up to nearly $8\%$ along $(100)$, $(010)$, $(110)$, and $(\bar1\bar10)$ stabilizes qualitatively distinct responses rather than different strengths of one phase [2605.27979].

The selected channel depends on direction. Along $(100)$, the response is predominantly cooperative Jahn-Teller: $Q_{\rho,\mathrm{avg}}$ grows, $Q_{\rho,\mathrm{diff}}$ remains vanishingly small, and the orbital pattern tends toward staggered $|3y^2-r^2\rangle/|3x^2-r^2\rangle$ above about $4\%$. Along $(010)$, both Jahn-Teller and breathing amplitudes increase, but $Q_{\rho,\mathrm{diff}}$ and $Q_{1,\mathrm{diff}}$ become much larger above about $4\%$, and both Mn sites evolve toward $|x^2-y^2\rangle$. Along the diagonals, strain is applied along Mn–O bonds, suppresses the in-plane rotation $\phi$, and produces strong site selectivity in both $\rho_i$ and $Q_1$, stabilizing the CF phase with ferro-orbital order within one Mn sublattice and C-type charge disproportionation between sublattices. Under biaxial strain the FM$\to$A-type AFM boundary appears near $4\%$, whereas under diagonal uniaxial strain the FM metallic ground state remains stable to at least $8\%$ [2605.27979].

A different but related selector appears in the spatiotemporal TDANN model of primate MT. There the control parameter is the balance between a contrastive objective and a spatial regularizer,
$$
\mathcal{L}_{\mathrm{total}}
=
\mathcal{L}_{\mathrm{contrast}}
+
\alpha \sum_k \mathcal{L}_{\mathrm{spatial},k},
$$
with
$$
\mathcal{L}_{\mathrm{spatial},k}
=
\frac{1-\mathrm{corr}(\mathbf r_k,\mathbf D_k)}{2}.
$$
Varying $\alpha$ induces distinct map regimes: fragmented maps with many defects at $\alpha=0$; an intermediate regime near $\alpha\approx 0.5$ with smooth direction slabs, pinwheels, fraction of direction-selective units $(\mathrm{DSI}>0.5)$ of about $72\%$, median DSI about $0.68$, median circular variance about $0.73$, primary FWHM of $68.2^\circ$, and pinwheel density near $6\,\mathrm{mm}^{-2}$; then defect proliferation near $\alpha\approx 1.25$; and oversmoothing near $\alpha\approx 2.5$ [2605.11718].

The paper explicitly stops short of claiming established criticality. It reports bifurcation-like behavior, re-entrant defect density, and “qualitative hallmarks” suggestive of competing-energy landscapes, while noting that no finite-size scaling, critical exponents, or precise $\alpha_c$ are provided. This is therefore a near-critical or phase-like use of direction-selective organization rather than a demonstrated universality class [2605.11718].

## 6. Universality, distinctions, and open problems

The literature assigns different meanings to universality. In the rotor case, universality refers to the collective geometric expansion of $\operatorname{Tr}W$ within a broad class of finite-dimensional trigonometric mean-field interactions, with model dependence reduced to a finite set of collective couplings and closure data [2603.29074]. In nemato-elasticity, universality refers to compatibility itself: the gauge constraints of elasticity suppress incompatible fluctuations independently of crystalline anisotropy in the ideal medium [2507.23754]. In directed hypergraph percolation, by contrast, universality can break down: anomalous exponents depend on heavy tails, maximal correlations, anchors, and whether node or hyperedge percolation is performed [2601.20726]. In BNZS, the extracted exponents are reported to be smaller than in conventional BEC- or Ising-like transitions, consistent with criticality confined to a restricted sector rather than global softening [2604.20173].

Several misconceptions are explicitly corrected by these works. Direction-selective criticality is not synonymous with ordinary anisotropy in real space. It may be selection in order-parameter space, as in the eigenvectors of $C(E)$ in rotor Hamiltonians; symmetry-sector selectivity, as in $\Delta_L\to 0$ while $\Delta_H$ stays finite in BNZS; compatibility-selected momentum directions in nematicity; forward/backward branching channels in directed hypergraphs; or crystallographic strain selection of Jahn-Teller versus breathing instabilities in manganites [2603.29074, 2604.20173, 2605.27979, 2601.20726]. It also does not, by itself, fix transition order. The rotor framework requires additional thermodynamic analysis to distinguish first- from second-order transitions, and the MT study does not establish a critical point despite reporting phase-like behavior [2603.29074, 2605.11718].

The outstanding problems are correspondingly heterogeneous. Rotor models require branch-dependent analysis when multiple competing phases are present, including expansions around symmetry-broken branches [2603.29074]. Nemato-elastic theory leaves open dynamic criticality, explicit lattice-anisotropic implementations, and nonlinear fluctuation effects beyond the Gaussian sector [2507.23754]. BNZS motivates angle-dependent field studies, neutron spectroscopy of branch-selective softening, and further tests of partially critical phases in spin-orbit-coupled rare-earth magnets [2604.20173]. The directed-chain model leaves open continuous-time extensions, within-level spatial structure, and the behavior under strongly autocorrelated inputs [2201.07075]. The MT model identifies dense $\alpha$ sweeps, finite-area scaling, and defect statistics as the next step if a true critical point is to be claimed [2605.11718].

A plausible synthesis is that direction-selective criticality names a constrained route to instability: a system may approach criticality not by uniformly softening all relevant fluctuations, but by reorganizing its accessible fluctuation space so that only a selected subset becomes soft. The selector can be geometric curvature, compatibility, anisotropic exchange, non-reciprocal propagation, higher-order functional asymmetry, uniaxial strain, or an optimization trade-off. What unifies these otherwise disparate cases is the replacement of global softening by restricted critical channels.

Source: https://www.emergentmind.com/topics/direction-selective-criticality