Resolution of the Cosmic-Ray Muon Puzzle: Non-Linear Substrate Shear Bifurcation, Neutral Pion Channel Suppression, and Excess Muon Production in Ultra-High-Energy Air Showers

Combined observations from major extensive air shower (EAS) arrays—including the Pierre Auger Observatory, IceCube, Telescope Array, and KASCADE-Grande—have established an acute anomaly in ultra-high-energy astroparticle physics known as the Cosmic-Ray Muon Puzzle. When primary cosmic rays exceeding E_0 ~ 10^17 eV (√s ≥ 14 TeV, beyond the Large Hadron Collider regime) collide with atmospheric nuclei, the measured ground-level muon count (N_μ) exceeds the predictions of state-of-the-art QCD hadronic interaction models (EPOS-LHC, QGSJET-II-04, SIBYLL 2.3d) by 30% to 80% (>8σ significance). This excess grows logarithmically with primary energy and cannot be resolved by shifting primary composition toward heavy iron nuclei without violating independent fluorescence measurements of the shower maximum depth (X_max). Building upon our foundational series on the physical vacuum as an elastically pre-stressed continuum, we show that the muon excess is an intrinsic elastodynamic consequence of non-linear substrate shear bifurcation at trans-LHC collision strain rates. In our ontology, neutral pions (π^0 -> 2γ) are non-topological breathing dilatation modes, whereas charged pions (π^±) and muons (μ^±) are radially excited torsional vortex solitons stabilized by the holding pre-stress P_0 ≈ 7.48 × 10^24 Pa. Above the critical onset threshold E_th ≈ 1.0 × 10^17 eV, non-linear bulk hardening resists volumetric breathing modes, bifurcating collision strain energy preferentially into long-lived torsional vortex solitons. We derive the analytical muon enhancement factor R_μ(E_0) = 1 + 0.097 ln(E_0 / E_th), accurately reproducing the empirical 22% excess at 10^18 eV and 45%–60% excess at 10^19–5×10^19 eV while preserving X_max concordance under Occam's razor without exotic quark fireballs or Lorentz invariance violation.