Topological Foundations of Multi-Field Instabilities in Continua: Part 1: Foundations Part 2:Analytical Formulation for 1-D Spin Chains Part 3: Numerical Upscaling

Technical Summary: Topological Foundations of Multi-Field Instabilities in Continua: Part 2 — Analytical Formulation for 1-D Spin Chains

Problem Statement
Granular materials exhibit irreversible macroscopic deformation (dilatancy and failure) despite the time-reversibility of individual grain interactions. Conventional mechanics attributes this irreversibility primarily to frictional dissipation at grain contacts. This paper argues that such an attribution is incomplete, proposing that geometric coupling mechanisms at the grain-contact scale play a fundamental role independent of friction. The core problem addressed is the identification of a parameter-free, topological origin for configurational irreversibility in one-dimensional granular chains, specifically isolating the mechanism that drives dilatancy-driven matrix instability without relying on empirical energy landscapes or friction coefficients.

Methodology
The paper develops an analytical formulation based on the augmented Onsager framework ($A = D + L$), where $D$ is a symmetric dissipative block and $L$ is a skew-symmetric conservative block. The methodology proceeds through the following steps:

1. Topological Decomposition: The granular contact state is decomposed into a minimal Volumetric–Mechanical–Configurational (VMC) triad ($N=3$).
   • $V$: Volumetric channel (overlap/normal force).
   • $M$: Mechanical channel (tangential/shear transmission).
   • $C$: Configurational channel (fabric orientation/rotation of contact normals).
2. Graph-Theoretic Filtering: The 1-D granular chain is modeled as a path graph $P_N$. The paper utilizes the Hodge–Helmholtz decomposition to show that in 1-D, the first Betti number $\beta_1 = 0$. This topological constraint acts as a "cocycle filter," eliminating all solenoidal (curl-driven) cycle modes and restricting the system's evolution entirely to the cocycle (gradient) subspace.
3. Lie-Algebraic Mapping: The reversible coupling block $L$ for the $N=3$ triad is identified as an element of the Lie algebra $\mathfrak{so}(3)$. The configurational channel $C$ is reinterpreted as a classical spin degree of freedom. The operator $L$ acts as a rotation generator, where the state vector precesses about a coupling axis $\omega$.
4. Analytical Derivation: The paper derives closed-form analytical solutions for the linear evolution equation $\dot{q} = Lq + f_0$ (in the frictionless limit $D=0$) for both $N=3$ (odd) and $N=4$ (even) channel counts. These solutions are compared against direct numerical integration to validate the algebraic predictions.

Key Contributions

• The VMC Triad as a Minimal Gateway: The paper establishes that the three-channel VMC contact is the minimal topological unit capable of generating a "Gateway Layer." Unlike even-dimensional systems ($N=2, 4$) which form "Stable Layers" with closed symplectic orbits, the odd-dimensional ($N=3$) system possesses a structural zero eigenvalue mandated by the Parity Theorem ($\det(L)=0$).
• Spin-Reinterpretation of Irreversibility: By mapping the contact dynamics to $\mathfrak{so}(3)$, the paper demonstrates that the parity-mandated null mode corresponds exactly to the rotation axis of the coupling generator. Because a rotation operator exerts no torque along its own axis, any boundary forcing projected onto this axis cannot be dissipated or redirected by the conservative block. This results in a secular drift (unbounded growth) along the null direction, providing a geometric explanation for configurational irreversibility that persists even in the absence of friction.
• The Cocycle Filter: The paper rigorously defines the 1-D Cartesian restriction not as a simplification, but as an exact topological filter. It proves that the acyclic nature of the 1-D chain ($\beta_1=0$) mathematically forbids Schnakenberg cycle currents and curl-driven shear localisation. Consequently, the only available irreversible pathway in 1-D is the cocycle-borne null-mode drift, which manifests macroscopically as dilatancy.
• Frictionless Irreversibility: The analysis demonstrates that the onset of instability is a consequence of contact topology (the odd channel count $N=3$) rather than contact friction. The "arrow of time" in this system is driven by the non-retraceability of the forced path along the rotation axis, independent of the friction coefficient.

Results

• Parity-Driven Dynamics:
  • Odd $N$ (Gateway Layer): The system exhibits a secular drift along the null eigenvector $v_0 \propto (L_{MC}, 0, L_{VM})^T$. Under sustained axial boundary forcing, the null-mode amplitude grows linearly with time ($v_0^T q(t) \propto t$), representing an unbounded escape from the initial state.
  • Even $N$ (Stable Layer): The system exhibits harmonic containment. The dynamics are confined to a compact invariant torus, resulting in bounded, quasi-periodic oscillations about a shifted equilibrium.
• Closed-Form Solutions: The paper provides explicit analytical expressions for the state evolution $q(t)$ for $N=3$ and $N=4$. These solutions confirm that the odd-N system drifts secularly while the even-N system remains bounded, validating the theoretical contrast without numerical approximation.
• Activation Mechanism: The paper clarifies that in the frictionless limit ($D=0$), the drift is exact and unconditional. In the presence of dissipation ($D > 0$), the Gateway number $G = \|L\|_2^2 / (\sigma_{min}(D)\sigma_{max}(D))$ determines the magnitude of the transient overshoot before relaxation, but the existence of the drift mechanism is topological and independent of $G$.

Significance and Claims
The paper claims to provide a parameter-free, topological foundation for one-dimensional granular irreversibility. Its primary significance lies in decoupling the origin of dilatancy and configurational drift from frictional dissipation. By identifying the VMC triad as the irreducible generator of instability, the authors argue that granular irreversibility is an emergent property of contact-space connectivity, analogous to the preserved spin axis of a classical gyroscope.

The authors assert that this framework offers a rigorous alternative to empirical constitutive modeling, replacing curve-fitting with algebraic necessity derived from the parity of thermodynamic channel counts. The work serves as the analytical bridge between the general topological principles established in Part 1 and the numerical verification and DEM correspondence detailed in Part 3. The paper explicitly limits its scope to the linear operator evolution, noting that nonlinear saturation and mappings to nonlinear wave equations are reserved for future work.