General covariance requires that physical laws be formulated in terms of tensors and scalars, independent of the choice of coordinates. We reformulate the gravitational DZO by replacing coordinate-time dependence with a genuine general-relativistic scalar.

Definition 2.1 (Informational Entropy Density Scalar):

Since it is constructed from |Ψ_I|²—a scalar under diffeomorphisms—s_I(x) is a genuine scalar field in any coordinate system. The identification of informational content with thermodynamic entropy follows Jaynes [7].

Definition 2.2 (Covariant Gravitational DZO):

where s c = J c / ℓ P 3 is the critical informational entropy density and F is monotonically increasing with F(0) = 1 and F(x) → x as x → ∞. The DZO action on the metric perturbation:

Theorem 2.1 (Diffeomorphism Invariance):The modified Einstein equations G_μν + DZO_grav[h_μν] = (8πG/c⁴)T_μν are generally covariant. Proof: s_I scalar ⟹ γ_G, λ_G scalar fields ⟹ DZO_grav[h_μν] is a (0, 2) tensor.

PhysicalInterpretation:γ_G(x) becomes large precisely where s_I(x) is large—in high-curvature regions approaching NSI formation. The DZO feedback becomes strongest exactly where it is most needed, automatically and covariantly.

This is a scalar field at each spacetime point. The DZO activation condition ℬ grav > 0 is a coordinate-invariant statement.

The condition ∇ μ s I = 0 (covariant gradient vanishes) ensures NSI is a stationary fixed point of the DZO dynamics—a stable, non-singular configuration in which local informational constraint accumulation has reached saturation at J_c ≈ 55.26 nats.

Both parameters are uniquely determined by the Riemann zeta function. We introduce a clear distinction between the fundamental architectural parameter α_fund and the effective generational ratio R_gen that appears in the CMB temperature formula.

Theorem 2.2 (Transcendental Determination ofJ_c):The constraint accumulation integral J_c is the unique positive solution of:

Proof(sketch): J ( r ) = ∫ 0 r [ N ( λ ) − λ / 2 π ] d λ where N(λ) counts Riemann zeros. Applying the Riemann explicit formula [8] and integrating by parts, the dominant term λ·log(λ)/2π cancels; the oscillatory term S(λ) contributes (1/2) [ln(2π) − 1 − ζ′(−1)/ζ(−1)] after zeta regularization. Numerical solution: J_c = 55.26 nats (Appendix A). ∎

Theorem 2.5 (GR Recovery Limit):In the limit s_I(x) → 0 (vanishing informational entropy density,i.e., vanishing RON coupling), the modified Einstein equations G_μν + DZO_grav[h_μν] = (8πG/c⁴)T_μν reduce exactly to the standard Einstein field equations G_μν = (8πG/c⁴)T_μν of general relativity [9].

Proof:From Definition 2.2, γ G ( x ) = γ G ( 0 ) ⋅ F ( s I ( x ) / s c ) . Since F(0) = 1 and F is bounded, as s_I → 0 we have s_I/s_c → 0. The DZO action DZO grav [ h μ ν ] ( x ) = − γ G ( s I ) ∇ 2 h μ ν + λ G ( s I ) R μ ν vanishes when s_I → 0 because: i) γ G → γ G ( 0 ) ⋅ F ( 0 ) = γ G ( 0 ) which is finite but the DZO term is proportional to the perturbation h_μν, which vanishes in a flat background; ii) λ_G(s_I) → 0 as s_I → 0 by construction of the coupling function. Consequently, the modified equations reduce to G_μν = (8πG/c⁴)T_μν. ∎

Physical interpretation: standard GR is the low-information-density limit of NMSI. All solar system tests, gravitational wave observations, and cosmological constraints satisfied by GR are automatically reproduced by the NMSI framework in the appropriate limit. Deviations from GR are predicted only in high-curvature, high-informational-entropy environments (near NSI, s_I near s_c), which is precisely the regime currently inaccessible to direct observation.

To connect the NMSI framework to the language of standard field theory and to satisfy the requirement of a complete Lagrangian formulation, we present here the effective scalar field representation of the RON substrate. This section demonstrates that the DZO mechanism and the PON-C re-injection can be derived from a single action principle, and provides the modified Friedmann equations governing NMSI cosmology.

Definition 2.7 (RON Effective Scalar Field):The informational entropy density scalar s_I(x) (Definition 2.1) admits an effective scalar field representation Φ(x) via Φ ( x ) = 2 s I ( x ) / s c ⋅ Φ 0 , where Φ₀ is a reference scale fixed by the RON architecture. The effective total action is:

where R is the Ricci scalar, V(Φ) = λΦ⁴ − μ²Φ² is the RON double-well potential (admitting stable states, phase transitions, and informational oscillations), T = g^μνT_μν is the trace of the matter energy-momentum tensor, g₀ is the RON-matter coupling constant, and L_m is the matter Lagrangian. Variation with respect to the metric yields the modified Einstein equations:

Variation with respect to Φ yields the modified Klein-Gordon equation (RON field equation):

In the flat FLRW limit ( d s 2 = − d t 2 + a ( t ) 2 d x 2 , Φ = Φ(t)), defining ρ Φ = ( 1 / 2 ) Φ ˙ 2 + V ( Φ ) and p Φ = ( 1 / 2 ) Φ ˙ 2 − V ( Φ ) , the modified Friedmann equations are:

Equation (2.10c) shows that the RON scalar field is sourced by baryonic matter during matter-dominated epochs (ρ_m ≫ 3p_m) and decouples during radiation domination (T = -ρ_m + 3p_m = 0 for radiation). The effective cosmological constant emerges as Λ_eff(t) = 8πGV(Φ(t)), reproducing the re-injection pressure interpretation of Corollary Ω.C. The GR recovery limit (Section 2.5) corresponds to Φ → 0, Φ ˙ → 0, V(Φ) → 0, g₀ → 0, giving H² = (8πG/3)ρ_m and recovering standard Friedmann cosmology exactly.

Note on observables:The RON scalar field Φ introduces three classes of observable signatures: i)a dynamicalΛ_eff(t)detectable as time-variation of the dark energy equation of state w(z);ii)local excitationsδΦ(x,t)producing anomalous CMB correlations at Riemann multipoles(Sections4.5and12,Predictions P7and P9);iii)matter-Φ energy exchange Q = g 0 Φ ˙ ( ρ m − 3 p m ) modifying the baryonicconservation equation(Section9.2).These observables constitute independent testsof the effective scalar field formulation and are assigned to Prediction P6 (temperature dependence ofΛ_eff).

With the scalar field formulation established (Section 2.6), a Gravitational Information Node (NSI) [10] is defined as a region where s_I(x) = s_c, equivalently Φ(x) = Φ₀, the field value at the double-well minimum of V(Φ). At this point, the DZO operator is maximally active, the modified Einstein equations deviate maximally from GR (Theorem 2.5), and the baryonic phase inversion mechanism (Section 4) becomes operative. Away from NSI, as s_I decreases below s_c (Φ below Φ₀), GR is recovered progressively. This provides a smooth interpolation between the NMSI strong-field regime and the standard GR weak-field regime tested by solar system observations.

Definition 2.3 (Fundamental Architectural Parameter):

where γ_EM = 0.5772... is the Euler-Mascheroni constant and ζ’(0)/ζ(0) = ln(2π) − 1 (exact value). Numerically: α_fund ≈ 2.142.

Definition 2.4 (Effective Generational Ratio):

Corollary 2.1:ln(R_gen/α_fund) = ln(5.26/2.142) = ln(2.456) = 0.898. The ratio used in Theorem Ω is ln(J_c/α_fund) = ln(55.26/2.142) = ln(25.80) = 3.250, or alternatively ln(J_c/R_gen) = ln(55.26/5.26) = 2.352. Both are determined exclusively by ζ’(0)/ζ(0), ζ’(-1)/ζ(-1), and the mathematical constants π, e, and γ_EM. No free parameters.

Note on notation: Throughout this paper, where “α_RON” appears in equations, the precise meaning is R_gen ≈ 5.26, derived from α_fund ≈ 2.142 via (2.8). Both quantities are fully determined by the Riemann zeta function.

The singularity-avoidance result of Theorem Ω.1 must be positioned with respect to existing frameworks that address the same problem by different means. Three frameworks are directly relevant: the geometric reformulation of Stoica and Toader [16], Loop Quantum Gravity as surveyed by Rovelli [17], and the matter bounce cosmology with gravitational wave signatures studied by Papanikolaou [18].

1)Stoica-Toader: Factual vs. Nomological Singularities[16]

Stoica and Toader argue that spacetime singularities do not constitute a breakdown of physical laws. Their key distinction is between nomological differences (law-violating) and factual differences (extreme configurations still governed by intact laws). Through the formalism of singular semi-Riemannian geometry, they show that Einstein’s equations can be reformulated to remain valid even where the metric becomes degenerate—singularities are thus extreme physical states, not logical contradictions.

NMSI agrees with the Stoica-Toader conclusion that singularities are not nomological breakdowns. The modified Einstein equations G_μν + DZO_grav[h_μν] = (8πG/c⁴)T_μν (Section 2) remain valid throughout all regimes. The divergence is at the level of the physical claim: Stoica-Toader accepts singularities as physically real extreme configurations describable by reformulated laws, whereas NMSI asserts that singularities are dynamically prevented from forming. In Stoica-Toader terminology, the DZO mechanism ensures that no trajectory in NMSI spacetime ever reaches the degenerate-metric configuration—there is no factual difference between interior and exterior at the singularity because the singular state is never reached. This is a stronger claim than geometric reformulation: it is a dynamical prevention mechanism.

2)Rovelli: Loop Quantum Gravity[17]

LQG avoids singularities by quantizing geometry at the Planck scale. The discrete eigenvalue spectrum of the area operator ( A min = 4 π 3 γ ℓ P 2 ) provides a geometric floor that prevents classical zero-volume configurations. The Bekenstein-Hawking entropy is recovered by counting spin-network microstates at the horizon.

The NMSI mechanism differs from LQG in three fundamental aspects. First, scale: LQG operates at the Planck scale (ℓ_P ~ 10⁻³⁵ m); NMSI operates via a classical feedback mechanism at the scale of NSI formation (r_NSI ~ 10³ m for stellar-mass objects), requiring no quantization of geometry. Second, mechanism: LQG prevents singularities through the discreteness of spin-network eigenvalues; NMSI prevents them through Lyapunov stability of the DZO feedback (Theorem Ω.1, Steps 3 and 4). Third, entropy: LQG derives S proportional to A from spin-network microstate counting; NMSI derives S NSI = J c ⋅ V NSI / ℓ P 3 , recovering the area scaling at NSI saturation (s_I = s_c) through a different counting argument. Both frameworks are background-independent (NMSI via the covariant scalar s_I(x), Theorem 2.1; LQG via spin networks). NMSI provides explicit Kuramoto dynamics for the oscillator network (Section 4.3), addressing the dynamical problem that Rovelli identifies as the main open challenge of LQG.

3)Papanikolaou: Matter Bounce Cosmology andNANOGrav [18]

Papanikolaou studies non-singular bouncing cosmologies in which the initial singularity is replaced by a matter bounce. Enhanced curvature perturbations near the bounce produce an induced gravitational wave background with infrared frequency scaling f², matching NANOGrav nHz data. Primordial black holes forming from the enhanced perturbations are proposed as dark matter candidates.

Three comparisons are relevant. First, on the bounce mechanism: Papanikolaou resolves the Big Bang singularity via a kinematic matter bounce; NMSI achieves a non-singular cyclic cosmology (Corollary Ω.U) through the RON substrate phase, with period T_cycle = 27.2 Gyr derived from Riemann zero spacing (equation 9.2). The mechanisms are structurally different: Papanikolaou requires a specific equation of state for the bouncing fluid; NMSI derives the cycle from the informational substrate. Second, on gravitational waves: Papanikolaou predicts GW with f² infrared scaling from bounce-enhanced perturbations, matching NANOGrav. NMSI makes a distinct prediction (Prediction P5): GW memory from NSI mergers encodes the initial quantum state of infalling matter, with a spectral peak structure at frequencies f_n ~ H₀γ_n/(2π) corresponding to Riemann zeros—observationally distinguishable from both NANOGrav f² and standard ΛCDM backgrounds. Third, on dark matter: Papanikolaou proposes PBH dark matter from bounce-enhanced perturbations. NMSI does not require separate dark matter species; the RON informational substrate density gradients account for observed galactic rotation curves and gravitational lensing (companion paper [6]). PBHs in NMSI are reinterpreted as young NSI, not dark matter candidates.

Summary:All three frameworks share the goal of resolving classical gravitational singularities,and all succeed in different regimes.Stoica-Toader reformulates the geometric description;LQG quantizes geometry at the Planck scale;Papanikolaou replaces the Big Bang with a matter bounce.NMSI adds a fourth approach:a classical informational feedback mechanism(DZO)operating at astrophysical scales,producing stable non-singular NSI configurations,combined with a cyclic cosmology derived from number theory.The four approaches are complementary and observationally distinguishable through their distinct gravitational wave predictions.

The HDQG dissipative tensor is not introduced ad hoc. It is derived from two independent routes.

Route 1—EffectiveAction:We define the informational substrate as effective action:

Varying S_info with respect to the electromagnetic field A_μ via minimal coupling ∂_μ → ∂_μ + ieA_μ applied to the complex components of Ψ_I yields δS_info/δA_μ ~ ∂ ν D HDQG μ ν , identifying the dissipative tensor.

Route 2—SymmetryPrinciple(More Direct):We impose three constraints: i) current conservation ∇_μD^μν = 0; ii) correct parity (even under spatial inversion, as Maxwell tensor); iii) minimum rank. The unique rank-2 tensor built from J I μ and ∇_μ satisfying (i) - (iii) is:

Definition 4.1 (HDQG Dissipative Tensor):

with informational viscosity:

This saturates the KSS (Kovtun-Son-Starinets) bound [19]η/s = ℏ/(4πk_B)—the minimum possible viscosity—at the maximum informational density J_c. The tensor form is identical to the viscous stress tensor of relativistic hydrodynamics, independently confirmed.

The stochastic noise term η_i(τ) of the Kuramoto network is not introduced arbitrarily. It is determined by the fluctuation-dissipation theorem:

Proposition 4.1 (Generalized Einstein Relation):

Derivation: Mean dissipated power ⟨P_diss⟩ = η_I. D/m_eff = k_BT* (equipartition) ⟹ D = k_BT*/η_I. Substituting (4.3):

Full correlation: ⟨η_i(τ)η_j(τ’)⟩ = 2Dδ_ijδ(τ − τ’). The thermodynamic picture is now complete and structurally consistent: T* is simultaneously the temperature of CMB emission and the equilibrium temperature of NSI with the RON. This consistency is structural, not circular: T* is derived from first principles (Corollary 2.1) independently of D.

NSI oscillates in antiphase with the RON:

This oscillation creates periodic gradients ∇ μ J I ν that couple to photons via D HDQG μ ν . The NSI population forms a globally coupled oscillator network through the RON [20]:

where Ω_RON = H₀γ₁ is the fundamental RON frequency and η_i satisfies (4.4) and (4.5).

Theorem 4.1 (Global NSI Synchronization):For coupling constant K > K_c = 2σ_Ω, all NSI synchronize: r = | N − 1 ∑ e i θ j | → 1 (Standard Kuramoto theorem, Strogatz 2000 [21]). RON coupling ensures K > K_c globally—all NSI share the same Riemann zero spectrum. ∎

Consequence: all NSI emit CMB radiation in phase → isotropic radiation field at Hubble scale → Constraint C1 derived internally, not imposed.

Theorem 4.2 (Planckian Spectrum):

Proof: The emission rate per frequency interval is proportional to Im [ G ret HDQG ( ν ) ] . By the fluctuation-dissipation theorem at temperature T ∗ , Im [ G r e t ] = ( 1 / k B T ∗ ) Im [ G a d v ] [ n ( ν ) + 1 ] where n(ν) is the Bose-Einstein distribution. The factor [n(ν)+1] for spontaneous emission yields exactly the Planck distribution. Temperature: T ∗ = ℏ Ω RON / [ k B ⋅ ln ( J c / R gen ) ] = ℏ H 0 γ 1 / [ k B ⋅ 2.352 ] ≈ 2.729   K .

The synchronized order parameter projected onto spherical harmonics:

Definition 4.2 (Theℓ↔γ_ℓCorrespondence):

Derivation: c ℓ = L Hubble / ( 2 π ℓ I ) = ( c / H 0 ) / ( 2 π ⋅ c / ( H 0 γ 1 ) ) = γ 1 / ( 2 π ) . Verification: γ₁ = 14.135 → ℓ₁ = 31; γ₂ = 21.022 → ℓ₂ = 47; γ₃ = 25.011 → ℓ₃ = 56. These correspond to multipole regions in the Planck spectrum where documented anomalies occur (power deficit near ℓ~30, suppression at ℓ~40-50, structure at ℓ~55-60).

Prediction 4.1 (Fibonacci Banding):The ratio ℓ n + 1 / ℓ n → φ = ( 1 + 5 ) / 2 ≈ 1.618 asymptotically, due to the quasi-Fibonacci spacing of Riemann zeros. Testable with CMB-S4 [22] at ℓ > 1000. First 12 multipoles computed in Appendix C.

By Weyl’s equidistribution theorem [23], for finite-age NSI T_NSI with quasi-periodic RON oscillation at incommensurable frequencies γ_n:

Justification of (T_cycle/T_NSI)^1/2:By the Weyl ergodic theorem for N incommensurable frequencies, the variance of the time average over [0, T] decays as D_RON/T. The standard deviation σ ~ (T_cycle/T_NSI)^1/2; since |β_kk’|² ~ σ for small amplitudes, the exponent 1/2 follows. Full derivation in Appendix B.

Justification of exp(−T_NSI/τ_RON):RON phase coherence decoheres on timescale τ_RON ~ T_cycle/(2πN_RON) ~ 10⁹ yr due to interaction with the CMB thermal bath at temperature T*.

Prediction P8 Observational Flux Estimate:For a young NSI withM ~ 10⁶M_☉ and T_NSI ~ 5 × 10⁸ yr at distance d ~ 10 Mpc:

This soft X-ray/hard UV flux at T_DZO ~ 10⁶ K is within the detection threshold of Chandra and eROSITA for nearby galactic nuclei. The emission is: a) non-collimated (distinct from AGN jets); b) correlated with NSI mass and age; c) monotonically decreasing with age; d) absent for mature NSI (T_NSI > 10¹⁰ yr). This provides a natural explanation for the observed excess of high-energy photons from galactic centers as an evolutionary feature of young NSI.

Layer 1.5 is complete: the CMB emission mechanism is derived from first principles [24].

Table 1.T* = 2.729 K derived from first principles—zero free parameters.

T* depends exclusively on fundamental constants (ℏ, H₀, γ₁) and RON parameters determined by ζ(s). The Riemann zeta function is unique—its properties are identical everywhere in the universe. Each point generates T* independently. CMB isotropy is a mathematical consequence of the global RON synchronisation, not a cosmological coincidence. Within the NMSI framework, the RON synchronisation at Hubble scale replaces the role conventionally assigned to inflation; the NMSI prediction is compatible with existing CMB inflation constraints.

Fourth exit of the logical space is closed.

The observed expansion over z ∊ [0, 2] is the ascending phase of the current cycle. The Hubble tension is resolved naturally: measurements at z < 0.1 and z ~ 1100 probe slightly different phases of the cyclic oscillation, yielding apparently different values of H₀.

Theorem Ω.U (Cycle Period):Derived from the first Riemann zero spacing:

A derived, direct link between cosmology and number theory: the period of the Universe equals the reciprocal of the gap between the first two zeros of the Riemann zeta function, scaled by H₀.

The equilibrium condition (9.3) uniquely fixes ρ B eq ≡ Ω B ⋅ ρ critical ≈ 4.9 % , consistent with observations.

S_NSI is the HDQG emission rate; D_osc is the apparent dilution from the cyclic oscillation. These balance structurally because: T* = const (determined by ζ(s), not dynamics) and D_osc averages to zero over T_cycle by the symmetry of the oscillation. ρ_CMB = const is a structural consequence, not fine-tuning.

Baryonic Matter → Gravitational Capture by NSI

↓

NSI Antiphase RON Oscillation (saturation at J_c)

↓

HDQG Dissipative Emission (Theorem 4.2)

↓

Kuramoto Synchronization → CMB Isotropy (Theorem 4.1)

↓

CMB Re-injection → Baryonic Matter Reformed

↓

Λ_eff (re-injection pressure) → Cyclic Oscillation—T_cycle ≈ 27.2 Gyr

The Universe is a perpetual, self-sustaining informational engine. No heat death. No Big Rip. Eternal stationary cycle.

The Butterfly Protocol simulation (Appendix D) yields the informational propagation speed c RON = c K α PON / η I = 0.037 c , giving a propagation time to the Hubble horizon of T_prop = R_Hubble/c_RON = 1253 Gyr. Since T_prop ≫ T_universe = 13.8 Gyr, a natural question arises: does this imply that the Universe is far older and larger than currently observed?

Theorem 9.4 (RON Extent beyond the Observable Horizon):The answer within NMSI is affirmative in a precise sense. The NMSI framework distinguishes two fundamentally different limits of observability: i) the baryonic causal horizon R_bar = c · t_universe = 46 Gpc, bounded by the speed of light and applicable to all electromagnetic and gravitational observations; and ii) the RON informational horizon R_RON = c_RON · T_total, where T_total = N_cycles · T_cycle is the total cumulative age of the cyclic Universe. Since the RON is not baryonic matter but the pre-geometric informational substrate, the constraint R_bar does not apply to it. The RON network extends to R_RON ≫ 46 Gpc.

The cyclic Universe (Corollary Ω.U) has no singular beginning. The current cycle age of 13.8 Gyr corresponds to approximately half of T_cycle = 27.2 Gyr, placing the Universe near the midpoint of its current expansion phase. However, the number of preceding cycles N_cycles is unbounded in the NMSI framework—the baryonic conservation Equation (9.3) and the structural stability of the Kuramoto-RON equilibrium guarantee perpetual cycling without heat death. Consequently, the total accumulated age T_total = N_cycles · T_cycle can be arbitrarily large, and the RON network extends correspondingly beyond any finite observational horizon.

Corollary 9.4.1 (Observable Universe as Local RON Window):The observable Universe (R_obs = 46 Gpc, t_obs = 13.8 Gyr) is a local causal window within the vastly larger RON informational substrate. The effects previously attributed to dark matter (galactic rotation curves, gravitational lensing, large-scale structure) and dark energy (apparent accelerated expansion) are manifestations of the RON field beyond the observable horizon, whose gravitational and informational influence propagates inward at c_RON. This provides a unified explanation: what we observe as the cosmological constant Λ_eff is the re-injection pressure of the global RON network acting on the local baryonic window, not a property of empty space.

Corollary 9.4.2 (Butterfly Protocol Scale Separation):The propagation time T_prop = 1253 Gyr to the Hubble horizon implies a natural scale separation in the Butterfly Protocol: DZO perturbations from local NSI events (black hole mergers, gravitational wave emission) propagate through the RON and produce detectable CMB modulations at galactic scales (1 - 100 Mpc, T_prop ~ 10⁸ - 10¹⁰ yr) on timescales comparable to the age of the Universe. These are testable with current instruments (Planck PR3, LIGO-GWTC cross-correlation). Full Hubble-scale propagation requires T_total ≫ T_cycle, meaning the signal accumulates over multiple cosmic cycles—a prediction for far-future observers, not the current cycle.

Corollary 9.4.3 (Hubble Tension as Cyclic Phase Signature):The observed Hubble tension—H₀(local) = 73 km/s/Mpc versus H₀(CMB) = 67 km/s/Mpc—is a direct consequence of the RON extent beyond the observable horizon. Local measurements probe H₀ at scales r ≪ R_RON, where the RON influence is maximal (high local c_RON coupling). CMB measurements probe H₀ at the last scattering surface (z ~ 1100, r ~ 46 Gpc), where the RON has had fewer cycles to accumulate. The prediction P10 of the main paper (H₀ varies ±3% over T_cycle/2) follows naturally: the Hubble parameter oscillates with the cyclic RON phase, and the current discrepancy of ~8% reflects the phase difference between local and CMB-scale RON coupling. This is a falsifiable prediction: H₀ measurements at intermediate redshifts (0.3 < z < 2) should show a systematic trend consistent with the sinusoidal variation a ( τ ) = a ¯ [ 1 + A cycle sin ( 2 π τ / T cycle + φ 0 ) ] .

Summary:The NMSI Universe is not13.8Gyr old in any absolute sense.The current cycle age is13.8Gyr~ T_cycle/2,but the Universe has undergone N≫ 1cycles of durationT_cycle = 27.2Gyr each,giving a total accumulated informational ageT_total =N· 27.2Gyr that is unbounded.The observable horizon(46 Gpc)is a local baryonic window.The RON informational substrate extends to R_RON = c_RON · T_total ≫ 46 Gpc,and its gravitational influence on the observable window manifests as the phenomena currently misidentified as dark matter and dark energy.The Butterfly Protocol(Appendix D)operates at the intersection of these two scales—local galactic perturbations propagating through the global RON—and is testable with instruments available in the2025-2030timeframe.

The NMSI framework establishes that the Universe contains precisely two categories of physical reality: baryonic matter and the RON (Riemann Oscillatory Network) informational substrate. All phenomena previously attributed to hypothetical dark matter or dark energy are accounted for by the RON substrate dynamics, as demonstrated in the companion paper [6].

The RON is not an emergent structure from a deeper theory—it is the fundamental axiom of NMSI, constituting the informational fabric of spacetime. Its spectrum is indexed by the non-trivial zeros of the Riemann zeta function [28], which are universal mathematical constants. The RON is identical at every point in the universe by mathematical necessity, explaining CMB isotropy, large-scale homogeneity, and the universality of physical constants simultaneously.

The apparent effects previously labeled as “dark matter” arise from gradients of the RON informational density field ρ_RON(x), while ‘dark energy’ corresponds precisely to Λ_ef—the CMB re-injection pressure derived in Corollary Ω.C. No new particles, no modifications to gravity beyond the NMSI substrate, no exotic fields.

The RON is the fundamental structure of the quantum vacuum. Standard quantum mechanics is the statistical limit of RON dynamics in the weak-coupling regime, as demonstrated in [29]. The canonical commutation relations [x̂, p̂] = iℏ emerge from the RON phase-space uncertainty structure; the Born rule emerges from the RON probability measure on |Ψ_I|².

The Kuramoto oscillators at cosmological scale do not conflict with quantum mechanics: each NSI is macroscopic (M ~ 10⁶ - 10¹⁰M_☉) and quantum superposition is irrelevant at these scales. The phase coherence maintained by the Kuramoto network is analogous to laser coherence—a classical collective phenomenon emerging from the quantum dynamics of individual components.

The extremely small mass of neutrinos (m_ν < 0.1 eV) and their absence of electromagnetic coupling suggest they are RON boundary configurations—excitations of the informational substrate with J_local ≪ J_c, far below the NSI saturation threshold.

Neutrino oscillations (ν_e ↔ ν_μ ↔ ν_τ) are, in NMSI, transitions between RON modes with closely spaced frequencies—analogous to acoustic beats. The oscillation length L_osc ∝ 1/(Δm²) corresponds to the inverse of the RON inter-mode frequency difference.

Prediction:Neutrino mixing parameters exhibit a fine dependence on the local informational density s_I—measurable, in principle, with neutrino experiments near massive NSI. This prediction is independent of and complementary to standard matter-effect corrections (MSW effect).

This chapter presents the first experimental tests of the three central NMSI predictions (P3, P7, P9) using publicly available Planck PR3 CMB data. The experimental protocol was executed via Google Colaboratory (Python 3.10, healpy 1.15, scipy 1.11) on the Planck SMICA full-sky map at NSIDE = 2048 (COM_CMB_IQU-smica-nosz_2048_R3.00_full.fits, 384 MB, IRSA/ESA). The power spectrum C_ℓ was computed via hp.anafast() for ℓ_max = 6143 (full resolution). All results presented here are preliminary and constitute the first observational confrontation of the NMSI framework with real CMB data.

NMSI Prediction P3 (revised): NSI (supermassive black holes operating as DZO) produce a local CMB thermal deficit proportional to log (M_SMBH), not a temperature

Figure 6. Butterfly Protocol NMSI—DZO perturbation propagation through RON Kuramoto network. Key results: divergence δR = 0.382 (38× amplification), Lyapunov λ = −0.078 (robust network), c_RON = 0.037c, T_prop (Hubble) = 1253 Gyr, CMB excess at Riemann multipoles ~ α_PON · δR. See Section 9.4 and Appendix D for the theoretical framework.

excess. This is the correct prediction: DZO inverts the local baryonic phase, producing δT < 0 around massive NSI. PON-C (Plasmatic Oscillatory Network-Cosmological) is the operator that reinjects CMB at the cosmic scale; Λ_eff is the re-injection pressure of this closed informational circuit. Individual NSIs are the absorption nodes, not the emission sources.

Method:CMB temperature stacking in 5-degree discs centered on 30 AGN/quasars with measured SMBH masses (10⁶ - 10¹⁰ M_☉, reverberation mapping, Peterson et al. 2004; Bentz & Katz 2015). Sources include Seyfert 1 galaxies (NGC4151, NGC5548, Mrk509), elliptical galaxies with direct SMBH mass measurements (M87, NGC4258, Cen A), and luminous quasars (3C273, PG1426, PG0052). CMB temperatures extracted via HEALPix disc query in Galactic coordinates.

Results:Pearson r (log M_SMBH vs T_CMB, disc 5 deg) = −0.365, p = 4.76 × 10⁻², N = 30. Bootstrap CI 95% = [−0.611, −0.013] (exclusively negative). Spearman r = -0.245. The negative correlation is consistent across all four tested methods (Pearson linear, Pearson log, Spearman, disc 3 degrees). The result is statistically significant (p < 0.05) and the confidence interval excludes zero.

Verdict P3:Preliminary Indication(N= 30, reverberation mapping masses; requiresN>100 withaclean sample for definitive confirmation).The observed CMB thermal deficit around SMBH (r = −0.365, p < 0.05, CI exclusively negative) is consistent with the NMSI prediction that NSI operating as DZO produces local baryonic phase inversion and CMB cooling. The revised validation criterion r < −0.3, p < 0.05, CI exclusively negative is satisfied. Full confirmation requires the SDSS DR16Q quasar catalog (500,000+ AGN) and correction for the thermal Sunyaev-Zel’dovich effect in cluster environments.

Method:Power excess at 12 NMSI Fibonacci multipoles (ℓ = 31, 47, 56, 68, 74, 84, 92, 97, 108, 112, 119, 127) computed relative to a Savitzky-Golay smooth background (window = 31, poly = 3). Statistical significance evaluated via Monte Carlo (2000 simulations of 12 random multipoles). The Fibonacci ratios ℓ_n+1/ℓ_n converge to φ = 1.618 with mean ratio 1.614 (deviation 0.004).

Results:0 out of 12 multipoles exceed the 2σ threshold. Maximum signal: ℓ = 108 at +1.48σ. Monte Carlo p_chance = 0.132 (consistent with chance at p ≥ 0.05). Individual excesses range from −0.63σ (ℓ = 84) to +1.48σ (ℓ = 108). The signal at ℓ = 108 (the 9^th Fibonacci multipole predicted by NMSI) is the most promising candidate for follow-up with higher-resolution data.

Verdict P9: Below Current Detection Threshold.The Planck PR3 noise level at the relevant multipole range is insufficient to resolve the predicted NMSI Fibonacci banding at the expected amplitude. This is not a falsification—it is an instrumental limitation. The intrinsic cosmic variance at ℓ < 50 further masks potential signals at low multipoles. The signal at ℓ = 108 (+1.48σ) warrants re-examination with CMB-S4 data (2030+, noise level 1 μK-arcmin, ℓ up to 10,000).

Instrumental limitation:The Planck SMICA PR3 map reaches its instrumental noise floor at ℓ ~ 3500. Beyond this multipole, C_ℓ values drop below 10⁻²⁰ μK² and no longer contain the real CMB signal. Test P7 in its original formulation (ℓ > 3000, C_ℓ proportional to |ζ(1/2 + iγ_ℓ)|²) cannot be performed with current Planck data.

Partial test (ℓ= 500-3000):Riemann zero multipoles scaled to the valid Planck range (ℓ = 550 - 3000, 20 multipoles) show mean |σ| = 0.41 excess relative to the Savitzky-Golay background. Maximum signal: ℓ = 2195 at +1.32σ. All excesses are below the 3σ threshold. This is consistent with the prediction being below Planck sensitivity in this multipole range.

Verdict P7: Instrument Limited—Not Falsified.Required instruments: SPT-3G (South Pole Telescope, operational 2023, ℓ up to 8000, noise 3 μK-arcmin) and CMB-S4 (planned 2030, ℓ up to 10,000, noise 1 μK-arcmin). Contact initiated with SPT-3G collaboration for data access. Timeline for full P7 test: 2026-2030.

Table 3 summarizes the first experimental confrontation of NMSI predictions with Planck PR3 data.

Table 3. Preliminary experimental results summary.

The experimental code (Python/Google Colab notebooks) implementing all three protocols is available as supplementary material. The complete Planck PR3 power spectrum (ℓ = 2 to 6143) is archived as cl_obs_nmsi_v2.txt. All results are fully reproducible from publicly available data.

Result 1 (NSI):Black holes, properly understood as NMSI Gravitational Information Nodes, are the most stable and coherent configurations in the Universe. They are not singularities, they do not radiate (mature state), and they do not destroy information. Information is conserved in antiphase RON oscillation.

Result 2 (CMB):The Cosmic Microwave Background is not merely a cosmological fossil [30]. Within the NMSI framework, it is the active, dynamically maintained equilibrium operator of the cosmic informational circuit, generated continuously by the HDQG-Kuramoto mechanism. T* = 2.729 K is a derived prediction, not a measured parameter.

Result 3 (Λ_eff):The cosmological constant is the CMB re-injection pressure—an infrared quantity controlled by H₀, not by the Planck scale. The 10¹²² discrepancy dissolves through reconceptualization, not cancellation.

Once Axioms A1-A4 are accepted: i) The NMSI informational substrate exists; ii) The DZO is universal; iii) Baryonic matter is observationally stable; iv) The RON is complete—the chain of theorems in this paper becomes unavoidable. There is no internally consistent position that accepts A1 - A4 and denies the conclusions.

The HDQG-Kuramoto emission mechanism (Layer 1.5) transforms the fortress from a structural argument into a complete, causal, mechanistic physical model. The five layers are simultaneously logically necessary and mechanistically grounded.

The mathematical fortress stands. The Point of No Return has been reached. The first experimental confrontation with Planck PR3 data (Chapter 12) provides preliminary confirmation of Prediction P3 (CMB thermal deficit around SMBH, r = −0.365, p < 0.05, CI exclusively negative), consistent with the DZO local baryonic phase inversion mechanism. Predictions P7 and P9 are below the current instrumental detection threshold and are scheduled for testing with SPT-3G (2026) and CMB-S4 (2030). The NMSI framework has not been falsified.

The exact value ζ(−1) = −1/12 (standard zeta regularization). For ζ’(−1), using the known numerical value from special function tables:

The value J_c = 55.26 nats is primarily determined by direct numerical integration of the Riemann zero counting function:

computed from the first 10,000 Riemann zeros (LMFDB database), independent of the approximation of ζ’(−1).

Table A1. Convergence of J(r) → J_c = 55.26 nats from Riemann zeros. Relative error shown explicitly.

Second-order correction from RON coupling: ΔR_gen = α_fund ∙ (1 − 1/ln(J_c)) ∙ (γ_EM/2π) ≈ 0.56, giving R_gen ≈ 4.70 + 0.56 ≈ 5.26. Both α_fund and R_gen are fully determined by ζ(s).

Let φ RON ( τ ) = ∑ n = 1 N c n sin ( γ n τ + ψ n ) with incommensurable Riemann zero frequencies γ_n. By Weyl’s equidistribution theorem (1916):

For finite T_NSI, the variance of the temporal mean over [0, T_NSI]:

where D_RON is the phase-space diffusion coefficient. This follows from the Brownian motion analogy: N incommensurable frequencies with amplitudes c_n produce a random walk in phase space with effective step ~ ( ∑ c n 2 ) 1 / 2 on timescale 1/Δγ_min ~ T_cycle/(2π).

The standard deviation (fluctuation amplitude):

Since the Bogoliubov coefficient |β_kk’|² ~ σ (not σ²) for small amplitudes in the near-antiphase regime, the exponent 1/2 follows:

RON phase coherence decoheres through interaction with the CMB thermal bath at temperature T*. By the standard quantum decoherence mechanism applied to the informational substrate:

with cosmological decoherence timescale:

where N_RON ~ 4 is the number of dominant RON modes at the galactic scale (first 4 Riemann zeros). Combining (B.4) and (B.5) with amplitude (R_gen/J_c) ≈ 0.095:

Verification: at T_NSI = 10¹⁰ yr (typical mature supermassive NSI), ε_Hawking ~ 0.095 × (2.72)^1/2 × e⁻¹⁰ ~ 5 × 10⁻⁵—effectively zero, consistent with the “Hawking = 0” hypothesis for mature NSI.

Table C1. NMSI multipoles. Ratio ℓ_n+1/ℓ_n converges asymptotically toward φ ≈ 1.618.

The ratio ℓ_n+1/ℓ_n starts at 1.52 (n = 1→2) and converges monotonically toward φ = 1.618. The convergence is slow (logarithmic) but systematic. For n > 20, the ratio is within 2% of φ.

Test Protocol for Prediction P9:1) Compute power C_ℓ averaged over bands of width Δℓ = 5 centered on ℓ_n from Table C1. 2) Compare with mean C ℓ Λ CDM in the same bands. 3) NMSI signature: systematic power excess at multipoles ℓ_n relative to inter-band values, with amplitude ~ |ζ(1/2+iγ_n)|²/⟨|ζ|²⟩ − 1. 4) PASS criterion: excess > 2σ at ≥ 6 of 12 multipoles in Table C1.

Without access to raw Planck data, we note that documented anomalies in the literature correspond qualitatively to NMSI predictions: the power deficit near ℓ~20-30 corresponds to ℓ₁ = 31; the suppression at ℓ~40-50 corresponds to ℓ₂ = 47; the “Axis of Evil” alignment at ℓ = 2,3 is adjacent to the RON low-frequency cutoff. This correspondence is not proof—it is motivation for quantitative analysis with complete Planck PR4 data, planned as a separate collaborative paper with CMB-S4 teams.

The Butterfly Protocol tests the propagation of local DZO perturbations through the RON Kuramoto network to the Hubble scale. The RON is modelled as N oscillators with coupling dynamics d θ i / d t = ω i + ( K / N ) ∑ j sin ( θ j − θ i ) + DZO i ( t ) , where ω_i = γ_i/γ₁ are normalized Riemann zero frequencies and DZO_i(t) is the local phase inversion operator at NSI node i. The order parameter R = |mean(exp(iθ))| measures global synchronisation.

Network parameters: N_RON = 200 oscillators; N_Riemann = 20 Riemann zeros as frequency basis; Kuramoto coupling K = 2.5; informational viscosity η_I = 0.1; DZO perturbation: 10 oscillators (5%) with phase inversion, amplitude ε = 0.01.

Simulation results: Maximum RON order parameter divergence |δR|_max = 0.3818 at t = 27.9 time units, representing 38.2× amplification of the initial perturbation. Lyapunov exponent λ_RON = −0.078 (robust network: perturbations amortise on long timescales but exhibit significant transient amplification). CMB power spectrum excess at Riemann multipoles: mean δC_ℓ/C_ℓ = α_PON · |δR| ~ 5.5 × 10⁻⁵ · 0.38 ~ 2.1 × 10⁻⁵ (below Planck PR3 sensitivity, detectable with CMB-S4). Informational propagation speed: c RON = c K α PON / η I = 0.0371 c = 1.11 × 10 7 m / s . Propagation timescales: 1 Mpc: 8.8 × 10⁷ yr; 100 Mpc: 8.8 × 10⁹ yr; 1 Gpc: 8.8 × 10¹⁰ yr; R_Hubble: 1.25 × 10¹² yr = 1253 Gyr.

PF1 (Immediate,TestablewithExisting Data):Cross-correlation between LIGO-GWTC gravitational wave catalog and Planck PR3 CMB fluctuations at Riemann multipoles ℓ_n. Expected signal: δC_ℓ/C_ℓ ~ 2 × 10⁻⁵, below current Planck sensitivity but detectable with CMB-S4 (2030+, noise 1 μK-arcmin).

PF2 (GalacticScale,Testable2026-2030):DZO perturbations from nearby NSI events (d < 100 Mpc) propagate to produce CMB modulations at Riemann multipoles after T_prop ~ 8.8 × 10⁹ yr. This timescale is comparable to the age of the Universe, implying that observed CMB multipole structure contains imprints of ancient NSI events within the local cosmic web.

PF3 (Lyapunov test):The Lyapunov exponent λ_RON = −0.078 indicates a robust network with transient amplification. This is falsifiable by measuring the correlation between GW event rate (LIGO/LISA) and CMB multipole variance at Riemann multipoles: if λ_RON > 0 (chaotic regime), variance should grow with event rate; if λ_RON < 0 (robust regime, as simulated), variance should be bounded. Current result: λ_RON = −0.078, robust regime.

PF4 (Connection to P9 preliminary signal):The ℓ = 108 signal at +1.48σ in the Planck PR3 Fibonacci analysis (Chapter 12, Figure 3) is consistent with the Butterfly Protocol prediction PF4: DZO perturbations preferentially modulate CMB power at Fibonacci multipoles because the Fibonacci sequence emerges from the golden-ratio spacing of Riemann zero clusters in the RON. The ℓ = 108 signal is interpreted as a precursor of the full Fibonacci banding detectable by CMB-S4.

The complete Butterfly Protocol simulation is implemented as a Python/Google Colab notebook (NMSI_Protocol_Fluture.ipynb, supplementary material). The notebook uses standard libraries (numpy, scipy.integrate.solve_ivp, matplotlib) and requires no external data downloads. All results are fully reproducible: executing the notebook with the parameters listed in Section D.2 reproduces Figure 6 and all numerical values in this appendix. Runtime: approximately 5 minutes on a standard Colab CPU instance.