--- title: 'Directional Stability: Definition and Applications' url: https://www.emergentmind.com/topics/directional-stability type: topic --- # Directional Stability: Definition and Applications Directional stability denotes a class of stability concepts in which robustness, convergence, or sensitivity is evaluated relative to a prescribed direction, orientation, or marching axis rather than only in an undirected neighborhood. The term has no single field-independent definition. In variational analysis it refers to direction-specific metric regularity and coderivative criteria for multifunctions [1508.01930; 1712.04704]. In numerical PDEs it denotes stability inherited from representative ordinary differential equations posed along time or space [2506.06543]. In control and mechanics it describes stability of yaw, bearing formations, or basin escape with respect to particular physical directions [1503.07258; 1508.06961; 1504.04476]. In machine learning it refers to convergence of normalized parameter and gradient directions under gradient flow [2006.06657]. In directional statistics it concerns the robustness of density ridges on the sphere under perturbations of the underlying density and under iterative ridge-finding algorithms [2104.14977]. In interface physics and heterostructure energetics it is tied to crystallographic direction, planar-front stability, and orientation-dependent interface enthalpy and band lineup [1112.4091; 1403.5230]. ## 1. Direction-dependent regularity in variational analysis A central mathematical formulation of directional stability is directional Hölder or Lipschitz metric regularity for a multifunction \(F:X \Rightarrow Y\) with closed graph, reference point \((\bar x,\bar y)\), direction \((u,v)\), and exponent \(\alpha \in (0,1]\). The property requires constants \(\delta,\varepsilon,\eta>0\) such that \[ \operatorname{dist}(x,F^{-1}(y)) \le \kappa [\operatorname{dist}(y,F(x))]^\alpha \] for points \((x,y)\) lying simultaneously in a small neighborhood of \((\bar x,\bar y)\), in a conic directional neighborhood \((\bar x,\bar y)+\operatorname{cone} B((u,v),\varepsilon)\), and under the gauge condition \[ \operatorname{dist}(y,F(x)) \le \eta \|(x,y)-(\bar x,\bar y)\|^{1/\alpha}. \] When this is imposed only for \(y=\bar y\), the resulting property is directional Hölder metric subregularity [1508.01930]. The directional formulation is explicitly weaker than classical, undirected metric regularity when \((u,v)=(0,0)\), because the gauge condition is essential in the directional setting. It is also strictly stronger than directional Hölder or Lipschitz metric subregularity. In the Lipschitz case \(\alpha=1\), the classical relation to the Aubin property of \(F^{-1}\) and to Lyusternik–Graves-type error bounds remains the background model, while the directional theory replaces global criteria by direction-restricted slope and coderivative tests [1508.01930]. The slope-based characterization uses the lower semicontinuous envelope \[ \varphi(x,y):=\liminf_{u\to x}\operatorname{dist}(y,F(u)) \] together with local and nonlocal slopes. For all \(\alpha \in (0,1]\), directional Hölder metric regularity is equivalent to positivity of the directional nonlocal slope limit inferior \[ \liminf_{(x,y)\to_{(u,v)}(\bar x,\bar y),\, \varphi(x,y)>0,\, \varphi(x,y)^\alpha/\|(x,y)-(\bar x,\bar y)\|\to 0} |\Gamma (\varphi^\alpha)(\cdot,y)|(x)>0. \] For \(\alpha=1\), positivity of the corresponding local slope limit inferior is also equivalent; for \(\alpha\in(0,1)\), the local slope condition is sufficient but not necessary, as illustrated by \(F(x)=(x_1-x_2)^3\) [1508.01930]. A complementary finite-dimensional framework is provided by directional limiting normal cones, subdifferentials, and coderivatives. For a closed set \(\Omega \subset \mathbb R^n\), the directional limiting normal cone \[ N_\Omega(x;u):=\{\xi \mid \exists t_k \downarrow 0,\ u_k\to u,\ \xi_k\to \xi,\ \xi_k \in N_\Omega(x+t_k u_k)\} \] refines the nondirectional Mordukhovich normal cone by retaining the approach direction. Analogously, the directional coderivative of \(F\) at \((x,y)\) in direction \((u,v)\) is \[ D^*F((x,y);(u,v))(\eta):=\{\xi \mid (\xi,-\eta)\in N_{\operatorname{gph}F}((x,y);(u,v))\}. \] This directional calculus supports direction-specific tests for metric subregularity, metric regularity, subtransversality, calmness, and the Aubin property under very weak directional qualification conditions [1712.04704]. ## 2. Perturbation stability and sensitivity analysis Directional stability is especially prominent in perturbation theory. In the variational-analytic setting, exclusion of directional criticality yields persistence of regularity under smooth perturbations. In Asplund spaces, the directional critical limit set \(Cr^\alpha F((\bar x,\bar y),(u,v))\) encodes scaled coderivative behavior along \((u,v)\). If \((0,0)\notin Cr^\alpha F((\bar x,\bar y),(u,v))\), then \(F\) is directionally Hölder metrically regular of order \(\alpha\); moreover, if \(g:X\to Y\) is differentiable near \(\bar x\) and satisfies \[ \|Dg(x)\|\le c\, m_F(x,y), \qquad c\in(0,1), \] on a directional neighborhood, then \(F+g\) is directionally Hölder metrically regular at \((\bar x,\bar y+g(\bar x))\) in direction \((u,v+Dg(\bar x)u)\) [1508.01930]. For Lipschitz perturbations in the case \(\alpha=1\), the modulus transforms explicitly. If \(F\) is directionally metrically regular at \((\bar x,\bar y)\) in direction \((u,v)\) with modulus \(\tau\), and if \(g\) is locally Lipschitz near \(\bar x\) with Lipschitz constant \(\lambda\) and Hadamard differentiable at \(\bar x\) in direction \(u\), then \(F+g\) remains directionally metrically regular whenever \(\lambda \tau<1\), with new modulus \[ \kappa'=\frac{\tau}{1-\lambda \tau}. \] In both the Lipschitz and Hölder perturbation results, the output direction is transformed by the first-order rule \((u,v)\mapsto (u,v+Dg(\bar x)u)\) [1508.01930]. The same directional viewpoint enters parameterized optimization. For \[ (P_y)\quad \text{minimize } f(x,y)\ \text{ subject to }\ y\in g(x)-K, \] with feasible map \(\Phi(y)=\{x\in X\mid y\in g(x)-K\}\), solution map \(S(y)\), and value function \(v(y)\), directional metric regularity of \(G(x):=g(x)-K\) at \((\bar x,0)\) in direction \((h,d)\) yields constructive feasible-direction results and directional sensitivity bounds. If \(Dg(x_0)h-d\in T_K(g(x_0))\) and \(G\) is directionally metrically regular at \((x_0,0)\) in direction \((h,d)\), then \(h\) is a feasible direction. In finite-dimensional \(Y\), the same assumptions imply the directional Robinson condition \[ d\in Dg(x_0)X-T_K(g(x_0)). \] Under additional assumptions, the directional Hadamard derivatives of the value function satisfy \[ v'_-(0;d)\ge \inf_{x\in S(0)}\inf_{\lambda\in \Lambda(x)} D_yL(x,\lambda,0)d, \qquad v'_+(0;d)\le \inf_{x\in S(0)}\sup_{\lambda\in \Lambda(x)} D_yL(x,\lambda,0)d, \] and if each multiplier set \(\Lambda(x)\) is a singleton, then \[ v'(0;d)=\inf_{x\in S(0)} D_yL(x,\lambda(x),0)d. \] Directional regularity also yields feasibility estimates of the form \[ \operatorname{dist}(x,\Phi(y))\le \kappa [\operatorname{dist}(y,G(x))]^\alpha, \] which in turn imply directional Hölder or Lipschitz sensitivity of \(v(y)\) [1508.01930]. A related but distinct inverse-problem notion appears for the dissipative wave equation with boundary damping coefficient \(b\). For the measurement map \[ M(b)=\partial_\nu u_b|_{\Sigma_1}, \] the paper proves a directional Lipschitz stability estimate at the origin. For a fixed nontrivial \(h\in B_{1/2,1}(\Gamma_1)\), if the undamped solution satisfies \(h\,u_t^0\neq 0\) in \(L^2(\Sigma_1)\), then there exist \(K>0\) and \(p_0>0\) such that \[ K\|ph\|_{B_{1/2,1}(\Gamma_1)}\le \|M(ph)-M(0)\|_{L^2(\Sigma_1)} \] for \(00\) and \(\|u-v\|_{L^\infty}\le c-\rho\), then \[ \|m(u)-m(v)\|_{L^q(X)} \le \frac{1}{\rho}\|m(u)\|_{L^q(X)}\|u-v\|_{L^\infty(Y)}, \] \[ \|M(u)-M(v)\|_{L^q(X)} \le \frac{1}{\rho}\|M(u)\|_{L^q(X)}\|u-v\|_{L^\infty(Y)}. \] Under concavity hypotheses, the maximal solution map is Hadamard directionally differentiable; in the uniquely solvable case the derivative is the smallest solution of the linearized fixed-point equation \[ \delta=\Psi'((\mathbb S(u),u);(\delta,h)). \] The results are stated without sign restrictions on perturbation directions and apply simultaneously to elliptic and parabolic QVIs [2105.05895]. ## 3. Numerical and PDE notions of directional stability In discretization of advection–diffusion equations, directional stability is introduced through a directional-ODE framework that treats discrete PDE information as a representative ODE posed along a chosen marching direction, either time or space. For the general nonlinear ADE \[ \frac{\partial u}{\partial t}+v(x,t,u)\cdot \nabla u=\nabla\cdot(d(x,t,u)\nabla u)+f(x,t,u), \] the temporal-ODE scheme freezes space and marches in time, while the spatial-ODE scheme freezes time and marches in space [2506.06543]. For the 1D diffusion equation, the zeroth-order temporal representative ODE at node \(i\) is \[ \frac{du_i}{dT}=-A u_i(T)+B,\qquad A=\frac{2D}{(\Delta x)^2}, \] with exact update \[ u_i^{n+1}= \frac{B}{A}+\left(u_i^n-\frac{B}{A}\right)e^{-A\Delta t}. \] Higher-order temporal schemes approximate the neighbor sum by a polynomial and preserve the exact exponential damping \(e^{-2aT}\). The paper defines directional stability as stability inherited from the exact solution of the representative ODE along the chosen direction and assessed on the discrete interval. For diffusion-dominated ODEs, the temporal kernel \(e^{-2aT}\) and the spatial kernels \(e^{\pm \sqrt A x}\) enforce monotone damping or regularization in the marching direction [2506.06543]. The nonlinear diffusion representative ODE \[ \frac{du_i}{dT}=D(u_i(T))[A u_i(T)+B]+C,\qquad A=-\frac{2}{(\Delta x)^2}, \] admits a stationary-point analysis. Ignoring sources, the steady state satisfies \(u^*=-B/A\), and perturbation by \(\epsilon(t)\) gives \(d\epsilon/dt=F'(u^*)\epsilon\) with \(F'(u^*)=A D'(u^*)\). Since \(A<0\) and \(D(u^*)>0\), the stationary point is asymptotically stable. For the multi-stage temporal-ODE diffusion scheme, the paper proves unconditional stability in the asymptotic discrete sense \(\Delta t\to \infty\): the limit becomes the Explicit Gauss–Seidel update for \(P=0\) and the Fully-Implicit update for \(P\ge 1\), both known to be unconditionally convergent for diffusion-type problems. This is contrasted directly with the CFL restriction of explicit FTCS diffusion schemes [2506.06543]. Directional solidification supplies a different PDE meaning. In a capillary-wave description of binary-alloy directional solidification, planar-front stability is analyzed relative to the pulling direction and temperature gradient. The effective Hamiltonian is written in terms of the interface position \(Z(\mathbf x,t)\) and excess solute concentration \(C(\mathbf r,t)\), and the externally imposed temperature gradient \(S\) and pulling velocity \(V_P\) enter through the driving force. Linearization about a steadily moving planar interface leads to a universal dispersion relation for monochromatic perturbations. Instability corresponds to \(\operatorname{Re}(\omega)>0\) [1112.4091]. Within that dispersion relation, the gradient parameter \(m\) acts as a long-wavelength cutoff. The low-\(q\) branch recovers the Mullins–Sekerka instability, but for \(m>0\) the threshold shifts to a finite \(q_c\) determined by \(\omega_1(q_c)=\omega_1'(q_c)=0\). A second branch has a finite gap at \(q=0\), and the neutral line \(F_P'(v_P)=0\) marks the onset of a Cahn-type oscillatory instability. Beyond threshold, the planar-front displacement \(h(t)=Z(t)-Z_P(t)\) satisfies the nonlinear oscillator \[ M(\dot h(t))\,\ddot h(t)+R(\dot h(t))+m^2 h(t)=0. \] The model therefore connects directional stability of the planar front to both the Mullins–Sekerka branch and the oscillatory branch, and it links periodic interface motion to solute banding and, when low-velocity segments enter the Mullins–Sekerka unstable regime, to banded structures [1112.4091]. ## 4. Control, mechanics, and robustness to finite perturbations In aircraft dynamics, directional stability is the classical yaw or weathercock stability of a vehicle with respect to sideslip. It is quantified by the yawing moment derivative with respect to sideslip, \(C_{n_\beta}\), with static directional stability corresponding to \(C_{n_\beta}>0\). The complete loss of a vertical stabilizer drives \(C_{n_\beta}\to 0\) and \(C_{n_r}\to 0\), destabilizing the Dutch-roll mode and rendering the spiral mode near-neutral or unstable [1503.07258]. One control response is differential thrust. In a damaged Boeing 747-100 model, the lateral/directional state is \(x=[\phi,p,\beta,r]^T\), and the damaged aircraft uses inputs \(u=[\delta_a,\delta T]^T\). Differential thrust is mapped from a rudder-equivalent yawing moment by \[ qSb\, C_{n_{\delta_r}}\, \delta_r = \delta T\, y_e, \] which yields \[ \frac{\delta T}{\delta_r}=4.43\times 10^5 \frac{\text{lbf}}{\text{rad}} \] for the specified flight condition. A Lyapunov-based model reference adaptive controller restores a stable flight envelope despite engine time constant and time delay, and the reported simulations show tracking errors converging to zero in approximately \(15\) s under the stated scenario and under \(30\%\) additive uncertainty [1503.07258]. A related robust-control design uses \(H_\infty\) loop-shaping. The damaged aircraft plant is again expressed in the state \([\phi,p,\beta,r]\) with inputs \([\delta_a,\delta T]\), and the loop-shaping synthesis employs pre- and post-compensators \(W_1(s)\) and \(W_2(s)\). The reported maximum stability margin is \[ e_{\max}=0.2763, \] and the closed-loop responses settle in approximately \(15\) s with feasible actuator effort. In that usage, directional stability means practical recovery of weathercock stability and yaw damping by a “virtual rudder” implemented through differential thrust, despite the loss of the vertical tail [1508.02487]. In multi-agent systems, directional stability appears in bearing-based formation stabilization with directed interaction topologies. Each agent obeys single-integrator dynamics \(\dot p_i=u_i\), uses relative position measurements to its out-neighbors, and applies the distributed control law \[ u_i(t)=-\sum_{j\in \mathcal N_i} P_{g_{ij}^*}(p_i(t)-p_j(t)). \] The closed-loop stacked dynamics are \[ \dot p(t)=-L_B p(t), \] where \(L_B\) is the bearing Laplacian. For directed formations, stability and convergence depend on the equality \[ \operatorname{Null}(L_B)=\operatorname{Null}(R_B), \] which defines bearing persistence. If bearing persistence fails, undesired equilibria appear and global formation stability cannot be guaranteed. If the formation is both infinitesimally bearing rigid and bearing persistent, then \[ \operatorname{Null}(L_B)=\operatorname{span}\{\mathbf 1_n\otimes I_d,\ p^*\}, \] so the limit formation has the same shape as the target up to translation and scaling [1508.06961]. A dynamical-systems notion of directional stability is given by the stability threshold \(\mathrm{ST}\), the smallest perturbation magnitude capable of pushing a system out of an attractor’s basin. For an attractor \(\mathcal A\), \[ \mathrm{ST}=\min_{x\in \partial \mathcal B(\mathcal A)} d(x,\mathcal A). \] For a unit direction \(u\), the directional threshold is \[ \mathrm{ST}(u)=\inf\{a\ge 0 \mid x_0+a u \notin \mathcal B(\mathcal A)\}. \] The equality \[ \mathrm{ST}=\min_{\|u\|_M=1}\mathrm{ST}(u) \] identifies the most dangerous direction of perturbation. The paper contrasts this quantity with Basin Stability and emphasizes that random sampling is inefficient near threshold, deriving the asymptotic law \[ 1-S_B \sim (q-\sigma)^{\frac{N-1}{2}} \] for fixed perturbation amplitude \(q\) approaching the stability threshold \(\sigma\) from above [1504.04476]. ## 5. Directional convergence in learning and statistics In deep learning, directional stability refers to stabilization of geometric directions during optimization even when parameter norms diverge. For deep homogeneous networks trained by gradient flow on classification losses whose minimizers are “at infinity,” the parameter trajectory satisfies \[ \dot \theta(t)\in -\partial L(\theta(t)), \] and the paper proves that the normalized weights converge: \[ \hat \theta_\infty=\lim_{t\to \infty}\frac{\theta(t)}{\|\theta(t)\|}. \] The normalized trajectory has finite length, so weights diverge only in norm while their direction stabilizes [2006.06657]. Homogeneity then transfers directional convergence to predictions and margins. With \(L\)-positive homogeneity \(f(x;\alpha \theta)=\alpha^L f(x;\theta)\), the normalized margins \[ \gamma_i(t):=\frac{y_i f(x_i;\theta(t))}{\|\theta(t)\|^L} \] converge for all training examples, and the margin distribution converges. Under locally Lipschitz gradients, the gradient direction also stabilizes and asymptotically aligns with the flow: \[ \lim_{t\to\infty}\frac{\nabla L(\theta(t))}{\|\nabla L(\theta(t))\|}= -\lim_{t\to\infty}\frac{\dot \theta(t)}{\|\dot \theta(t)\|}. \] The analysis is based on unbounded nonsmooth Kurdyka–Łojasiewicz inequalities for definable functions, and it yields consequences for margin maximization and convergence of saliency maps [2006.06657]. For directional data on the sphere, stability concerns geometric features of an estimated density rather than optimization trajectories in parameter space. On the unit hypersphere \(\Omega_q\), the Riemannian gradient is \[ \operatorname{grad} f(x)=(I_{q+1}-x x^T)\nabla f(x), \] and the paper defines a \(d\)-dimensional directional density ridge by \[ \bar R_d=\{x\in \Omega_q:\ \bar V_d(x)^T \operatorname{grad} f(x)=0,\ \bar \lambda_{d+1}(x)<0\}. \] Here \(\bar V_d(x)\) contains the eigenvectors associated with the smallest \(q-d\) eigenvalues of the spherical Hessian [2104.14977]. The directional stability theorem states that if the regularity conditions \((\bar A1)-(\bar A3)\) hold for \(f\), and if \(\|f-g\|_{\infty,3}^*\) is sufficiently small, then the perturbed density \(g\) also satisfies the ridge regularity conditions and \[ \operatorname{Haus}(\bar R_d(f),\bar R_d(g)) = O(\|f-g\|_{\infty,2}^*). \] For the plug-in directional kernel density estimator \(f_h\), this becomes \[ \operatorname{Haus}(\bar R_d(f),\bar R_d(f_h)) = O(h^2)+O_P\!\left(\sqrt{\frac{|\log h|}{n h^{q+4}}}\right). \] The associated directional SCMS algorithm operates intrinsically on the sphere and converges linearly in geodesic distance under the stated assumptions, with contraction factor \[ \bar \Upsilon=\sqrt{1-\frac{\bar \beta_0 \bar \eta}{4}}. \] Thus, in this setting, directional stability combines perturbation robustness of ridge geometry with convergence stability of a ridge-finding iteration [2104.14977]. ## 6. Orientation, interfaces, and crystallographic direction In materials theory, directional stability is tied to crystallographic orientation and interface type. For Si/ZnS abrupt and polar-compensated interfaces, the two-interface supercell energetics are written as \[ E_{I1}+E_{I2}=E_\Omega-\sum_\alpha n_\alpha \mu_\alpha, \] and interface energies are compared across the \([111]\), \([100]\), and \([110]\) directions [1403.5230]. The work emphasizes that directions lacking generalized mirror symmetry generate artificial bulk fields in two-slab supercells and complicate both band-offset determination and orientation-dependent stability. The adopted strategy is to place each distinct interface in a variety of supercell environments, determine valence-band offsets by linear extrapolation of the local potential from slab centers to nominal interface planes, and extract interface energies from overdetermined linear systems across several supercells. For sufficiently large supercells, the spread of offsets for a given microscopic interface is typically \(0.1\)–\(0.3\) eV, much smaller than the total interfacial distribution [1403.5230]. The orientation dependence is strong. The nonpolar abrupt \([110]\) interface has absolute interface energy \(2.81\) eV/nm\(^2\). For \([100]\), polar-compensated interfaces are generally more stable than abrupt ones across Zn-rich and Zn-poor conditions. For \([111]\), reconstructions of the single-bond interfaces are stable over nearly the entire chemical-potential window \(3.7\) eV \(<\mu_x<5.4\) eV, while abrupt triple-bond interfaces are energetically disfavored. The paper states the inequalities \[ E_{[110]}