Preprint
Article

This version is not peer-reviewed.

Quantum Informational Relativity: Foundational Structure, Microphysical Dynamics, Cosmological Predictions, and Numerical Validation (Complete Derivations and Unified Theoretical Framework)

Submitted:

21 November 2025

Posted:

25 November 2025

You are already at the latest version

Abstract
Quantum Informational Relativity (QIR) proposes a unified formalism in which microphysical dynamics, spacetime structure, and cosmological evolution are consistently derived from an underlying informational field-theoretic framework. This monograph consolidates the complete formulation of QIR from the operator-level action and classical limit to the Hamiltonian structure, consistency relations, and effective constants and organizes it into a coherent, self-contained theory. At the microscopic level, QIR predicts informational solitons, canonical fluctuation modes, confinement, mass generation, and small but measurable shifts in Standard Model observables. A Wilsonian mapping links these scales to cosmology, producing a smooth continuation across 30+ orders of magnitude. Cosmologically, QIR generates percent-level deviations in the growth function and lensing potential, easing the S8 tension while leaving background distances unchanged. Fully nonlinear N-body simulations performed with a modified SWIFT engine confirm these signatures and reproduce realistic cosmic structures. Altogether, this work provides the first complete exposition of QIR as a mathematically consistent, observationally testable, and scale-continuous alternative to standard microphysical and gravitational dynamics.
Keywords: 
;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  

Notations and Conventions

This monograph employs both classical and operator-valued structures. To ensure full clarity, we fix from the outset a strict and unambiguous distinction between classical fields, quantum operators, canonical variables, effective quantities, and geometric objects.

Classical and Quantum Quantities

We use the following global convention:
  • Classical objects carry no hat:
    p ( x ) , Z ( p ) , I μ ν [ p ] , T μ ν , g μ ν .
  • Operator-valued objects carry a hat:
    p ^ ( x ) , Z ^ ( p ^ ) , I ^ μ ν [ p ^ ] , L ^ , π ^ ( x ) .
  • The spacetime metric g μ ν and curvature tensors R μ ν , R , G μ ν are treated as classical throughout the theory and therefore never carry a hat.
  • Expectation values in the QIR Hilbert space H QIR are written
    O ^ = Ψ | O ^ | Ψ .
Normal ordering is assumed for all operator expressions involving p ^ unless stated otherwise.

Canonical Quantization

The conjugate momentum operator is defined by
π ^ ( x ) δ L ^ δ ( 0 p ^ ( x ) ) .
Equal-time canonical commutation relations are
[ p ^ ( t , x ) , p ^ ( t , y ) ] = 0 , [ π ^ ( t , x ) , π ^ ( t , y ) ] = 0 ,
[ p ^ ( t , x ) , π ^ ( t , y ) ] = i δ ( 3 ) ( x y ) .

Canonical Normalization of Fluctuations

Fluctuations around a background configuration p 0 are written
p ( x ) = p 0 ( x ) + δ p ( x ) , p ^ ( x ) = p 0 ( x ) + δ p ^ ( x ) .
The canonically normalized fields are
χ ( x ) = Z ( p 0 ) δ p ( x ) , χ ^ ( x ) = Z ( p 0 ) δ p ^ ( x ) ,
with canonical equal-time relations
[ χ ^ ( t , x ) , 0 χ ^ ( t , y ) ] = i δ ( 3 ) ( x y ) .

Mode Expansion

The normalized quantum fluctuation admits the standard decomposition
χ ^ ( x ) = d 3 k ( 2 π ) 3 1 2 ω k a ^ k e i k · x + a ^ k e i k · x ,
with
[ a ^ k , a ^ k ] = ( 2 π ) 3 δ ( 3 ) ( k k ) .

Internal Multiplets and Topological Sectors

When internal structure is required, we extend
p A ( x ) , A = 1 , , N ,
with operator counterpart p ^ A ( x ) . Phase maps and winding structures used for particle classification are always defined using the normalized multiplet.

Spacetime and Geometry

We work on a Lorentzian manifold ( M , g μ ν ) with signature ( , + , + , + ) . Indices are raised and lowered with g μ ν . Curvature tensors follow
G μ ν = R μ ν 1 2 g μ ν R .

Informational Modulation

Classically,
Z ( p ) = 1 + β p 2 .
The operator version is
Z ^ ( p ^ ) = 1 + β p ^ 2 .
The effective local light speed is
c ( x ) = c 0 Z ( p ( x ) ) .

Physical and Effective Metrics

The physical metric is g μ ν . An auxiliary effective metric is occasionally used:
g μ ν eff = Z ( p ) g μ ν ,
but all gravitational dynamics are defined with g μ ν .

Energy–Momentum Tensors

Matter: T μ ν (classical). Informational:
I μ ν [ p ] and I ^ μ ν [ p ^ ] .
Total conservation:
μ ( T μ ν + I μ ν ) = 0 .

Einstein Equation (QIR Form)

G μ ν = 8 π G c ( x ) 4 Z ( p ) T μ ν + I μ ν .

Units

Natural units = 1 are assumed. The factor c ( x ) is kept explicit to track informational modulation.

Fourier Transform

f ( x ) = d 3 k ( 2 π ) 3 f ˜ ( k ) e i k · x , f ˜ ( k ) = d 3 x f ( x ) e i k · x .

Functional Derivatives

δ S δ g μ ν , δ S δ p , Z ( p ) = d Z d p .
These conventions remain fixed throughout the monograph unless explicitly overridden for specific calculations.

1. Introduction

The search for a unified description of microphysical dynamics, spacetime geometry, and cosmological structure formation remains one of the central challenges in contemporary theoretical physics. Quantum field theory (QFT) successfully describes local excitations and particle interactions, while General Relativity (GR) governs the large-scale curvature of spacetime. Yet the conceptual foundations of these frameworks differ profoundly: QFT is rooted in operator-valued fields defined on a fixed background, whereas GR treats geometry itself as dynamical. Reconciling these perspectives has proven difficult, not only at the level of mathematical consistency, but also in terms of identifying a common set of underlying physical principles.
A recurring idea in several approaches to unification is that spacetime, fields, or interactions may arise from deeper informational, statistical, or emergent structures. These perspectives appear in contexts ranging from entanglement-based derivations of geometry to thermodynamic formulations of gravitational dynamics. While diverse in implementation, they share the notion that traditional geometric and quantum concepts might both originate from an informational substrate.
Quantum Informational Relativity (QIR) develops this line of thought in a specific and operationally well-defined manner. The central dynamical entity is an informational field p ( x ) whose modulation factor Z ( p ) affects inertial, interaction, and gravitational responses while preserving relativistic covariance. In this view, particles, interactions, and geometric curvature are not independent ingredients but manifestations of a single informational degree of freedom. The resulting structure retains the full tensorial and geometric content of GR, yet introduces a controlled and covariant deformation through the informational sector.
A distinctive feature of QIR is that informational fluctuations, their canonical normalization, and their topological properties naturally give rise to microphysical structures such as stable localized excitations, phase sectors, and emergent mass scales. At the opposite end of the spectrum, the same informational coupling modifies the growth of cosmic structures and the lensing potential at late times, providing a possible explanation for the mild yet persistent tension between weak-lensing surveys and Λ CDM predictions. Because the underlying modulation is the same at all scales, QIR establishes a continuous bridge between microphysical and cosmological regimes without introducing additional free parameters.
Beyond analytic calculations, the informational modification can be implemented directly in N-body simulations, allowing nonlinear structure formation to be tested in fully dynamical settings. This provides an essential consistency check: any viable extension of GR must recover realistic halo distributions and large-scale structures while remaining compatible with observed growth suppression.
The purpose of this monograph is twofold. First, it presents a complete and self-contained construction of QIR, from the operator-level action and canonical quantization to geometric consistency and scale transitions. This includes full derivations and explicit intermediate steps, ensuring mathematical transparency. Second, it synthesizes the observational and phenomenological consequences across microphysical, cosmological, and nonlinear regimes, demonstrating how a single informational field can produce coherent physics over many orders of magnitude.
The document is organized as follows. Section 2 develops the foundational operator structure of QIR, including the fundamental action, the informational stress tensor, the unified Einstein equation, the Hamiltonian formulation, canonical normalization, and the derivation of effective propagation and coupling scales. Section 3 presents the quantization of the informational field, including operator canonical quantization, mode expansion on curved backgrounds, the propagator, spectral decomposition, renormalized quantum observables, and the vacuum structure. Section 4 confronts the microphysical predictions of QIR with high-energy data, examining informational corrections, decay processes, scattering cross sections, spectral signatures, informational mixing, and combined collider constraints. Section 5 develops the cosmological sector, including background evolution, linear growth, weak lensing, CMB constraints, nonlinear structure formation, and unified cosmological tests. Section 6 presents the numerical validation using SWIFT, including the simulation pipeline, growth extraction, matter power-spectrum evolution, cosmic-web classification, velocity-field diagnostics, and global nonlinear validation. Section 7 investigates astrophysical and propagation signatures across gravitational potentials, halo dynamics, radiation transport, time-delay effects, and high-energy propagation. Section 8 provides a general discussion and comparative assessment of QIR, including structural synthesis, physical interpretation, internal consistency, comparisons with GR, Λ CDM, and modified gravity, limitations, and falsifiable predictions. Section 9 concludes the monograph.

2. Foundational Action and Operator Structure

The Quantum Informational Relativity (QIR) framework is built upon a single informational operator field p ^ ( x ) whose local amplitude modulates inertial, geometric, and coupling strengths through a dimensionless informational modulation function Z ( p ^ ) > 0 . The goal of this section is to establish rigorously the full operator-level dynamical structure: the fundamental action, its variations with respect to the metric and the informational field, the unified Einstein equation, and the covariant canonical quantization of p ^ ( x ) . All quantities are written in the physical g-frame unless explicitly stated.

2.1. Fundamental Operator Action

The Quantum Informational Relativity (QIR) framework is grounded in the postulate that the fundamental dynamical entity is an operator-valued informational field p ^ ( x ) defined on a Lorentzian manifold ( M , g μ ν ) . Its amplitude encodes local informational density, and its variations renormalize inertial and geometric responses through a positive definite, dimensionless modulation function Z ( p ^ ) .
In this subsection we construct the fundamental action, justify its operator structure, and examine its mathematical properties. We proceed step by step to make explicit the assumptions, domains of definition, and the covariant identities needed for the derivations that will follow in subsequent subsections.

2.1.1. Principles Guiding the Action

The informational field p ^ is assumed to satisfy the following general principles:
  • Covariance: all dynamical quantities must transform as well-defined geometric objects under diffeomorphisms.
  • Positivity: the informational modulation Z ( p ^ ) must be everywhere positive to ensure a well-posed kinetic structure and a hyperbolic operator equation of motion.
  • Hermiticity: the action must be self-adjoint in the Hilbert space H QIR ; this ensures that energies, momenta, and stress tensors are real-valued expectation values.
  • Minimal coupling: the covariant derivative D μ includes the Levi-Civita connection of g μ ν , and possibly internal gauge connections acting on multiplet extensions of p ^ .
  • Locality: the action depends on p ^ and its first derivatives only, consistent with standard field-theoretic structure.
These principles uniquely restrict the admissible form of the action.

2.1.2. Construction of the Operator Action

Let p ^ ( x ) be an operator-valued scalar field acting on the QIR Hilbert space H QIR . The most general diffeomorphism-invariant, Hermitian, local action built from p ^ , D μ p ^ , and g μ ν up to mass dimension four is
S [ p ^ , g μ ν ] = M d 4 x g 1 2 A ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) U ( p ^ ) ,
where A ( p ^ ) and U ( p ^ ) are operator-valued functions defined as power series in p ^ . Hermiticity requires A ( p ^ ) to be a real, positive, even function.
The essential feature of QIR is the identification
A ( p ^ ) = Z ( p ^ ) ,
where Z ( p ^ ) is the informational modulation governing inertial and coupling strength.

2.1.3. Informational Modulation and Its Properties

The modulation function Z ( p ^ ) encodes how the local informational density renormalises all dynamical responses. For consistency, it must satisfy:
  • Z ( p ^ ) > 0 for all admissible states, ensuring positivity of the kinetic operator.
  • Z ( p ^ ) is a self-adjoint operator function, i.e. Z ( p ^ ) = Z ( p ^ ) .
  • Z ( p ^ ) must commute with scalars but not necessarily with p ^ itself unless p ^ is taken to be Hermitian.
A widely used form is
Z ( p ^ ) = 1 1 + β p ^ 2 ,
which is positive and even. However, our derivations remain valid for any differentiable operator function Z ( p ^ ) .

2.1.4. Final Form of the Operator Action

With the above identifications, the fundamental QIR action is
S [ p ^ , g ] = d 4 x g 1 2 Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) V ( p ^ ) .
The potential V ( p ^ ) is chosen to be a real, bounded-from-below, self-adjoint operator function. A common example is
V ( p ^ ) = 1 2 μ 2 p ^ 2 + 1 4 λ p ^ 4 ,
though nothing in the subsequent derivations depends on this specific form.

2.1.5. Hermiticity and Operator Ordering

Because Z ( p ^ ) and p ^ may not commute, care must be exercised in the definition of ( D μ p ^ ) Z ( p ^ ) D ν p ^ . Two choices are common:
  • Symmetric ordering:
    1 2 Z ( p ^ ) g μ ν ( D μ p ^ ) D ν p ^ + D ν p ^ g μ ν ( D μ p ^ ) Z ( p ^ ) .
  • Left ordering:
    Z ( p ^ ) g μ ν ( D μ p ^ ) D ν p ^ .
QIR chooses the symmetric ordering to guarantee strict Hermiticity. Thus we interpret the kinetic term in (2) as
1 2 Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) 1 2 Z ( p ^ ) g μ ν ( D μ p ^ ) D ν p ^ + D ν p ^ g μ ν ( D μ p ^ ) Z ( p ^ ) .
In the semiclassical limit p ^ p , the two orderings coincide, ensuring consistency with classical dynamics.

2.1.6. Domain of the Covariant Derivative

The operator derivative D μ acts as
D μ p ^ = μ p ^ + Γ μ p ^ ,
where Γ μ is either:
  • the Levi-Civita connection acting on tensors (trivial for scalars),
  • an internal gauge connection acting on multiplets p ^ A , A = 1 , , N .
Thus for a scalar operator,
D μ p ^ = μ p ^ ,
but when internal indices are present,
D μ p ^ A = μ p ^ A + i ( A μ ) B A p ^ B .
All subsequent variations of D μ p ^ assume these transformation rules.

2.1.7. Classical Limit and Emergent Interpretation

Taking expectation values in a suitable state | Ψ gives
p ( x ) = Ψ | p ^ ( x ) | Ψ , Z ( p ) = Ψ | Z ( p ^ ) | Ψ ,
and the action reduces to its classical version:
S cl [ p , g ] = d 4 x g 1 2 Z ( p ) g μ ν μ p ν p V ( p ) .
This classical limit recovers all the expressions used in microphysical and cosmological analyses.
This concludes the full construction, motivation, and operator-level justification of the fundamental QIR action. The next Section 2.2 performs the complete metric variation, including the explicit treatment of
δ g μ ν , δ g μ ν , δ g , δ Γ μ ν α , δ ( D μ p ^ ) ,
and leads to the fully derived informational stress tensor I μ ν and the geometric sector of the unified field equations.

2.2. Metric Variation and the Informational Stress Tensor

In this section we perform the complete variation of the operator action (2) with respect to the metric g μ ν . This derivation is essential because it produces the informational energy–momentum tensor I ^ μ ν and establishes its role in the unified Einstein equation.
The variation involves multiple tensorial identities, variations of the determinant, inverse metric, and Christoffel symbols, as well as the metric dependence of the covariant derivative acting on p ^ . For full transparency, every required identity is derived explicitly.

2.2.1. Variation Identities for the Metric

We begin with standard results from differential geometry.

Inverse Metric

Since g μ α g α ν = δ μ ν , variation gives
δ g μ ν = g μ α g ν β δ g α β .

Determinant

Starting from δ g = g g μ ν δ g μ ν , one obtains
δ g = 1 2 g g μ ν δ g μ ν = 1 2 g g μ ν δ g μ ν .

Christoffel Symbols

The Levi-Civita connection Γ μ ν α = 1 2 g α λ ( μ g λ ν + ν g λ μ λ g μ ν ) varies as
δ Γ μ ν α = 1 2 g α λ μ δ g λ ν + ν δ g λ μ λ δ g μ ν .

Variation of the Covariant Derivative Acting on p ^

For a scalar operator-valued field p ^ ,
D μ p ^ = μ p ^ , δ ( D μ p ^ ) = 0 .
For a multiplet p ^ A with gauge connection A μ ,
D μ p ^ A = μ p ^ A + i ( A μ ) A B p ^ B , δ ( D μ p ^ A ) = 0 ,
because A μ does not depend on the metric.

2.2.2. Varying the Kinetic Term

The kinetic term of the action is
K = 1 2 g Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) .
Its variation is
δ K = 1 2 ( δ g ) Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) + 1 2 g Z ( p ^ ) δ g μ ν ( D μ p ^ ) ( D ν p ^ ) .
Using (6) and (5), this becomes
δ K = 1 4 g g α β δ g α β Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) + 1 2 g Z ( p ^ ) g μ α g ν β δ g α β ( D μ p ^ ) ( D ν p ^ ) .
Lowering indices on ( D μ p ^ ) yields
( D μ p ^ ) ( D ν p ^ ) g μ α g ν β = ( D α p ^ ) ( D β p ^ ) .
Thus
δ K = 1 2 g 1 2 g α β Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) Z ( p ^ ) ( D α p ^ ) ( D β p ^ ) δ g α β .

2.2.3. Variation of the Potential Term

The potential term is
U = g V ( p ^ ) .
Since V ( p ^ ) does not depend on the metric,
δ U = ( δ g ) V ( p ^ ) = 1 2 g g α β V ( p ^ ) δ g α β .

2.2.4. Total Metric Variation of the Action

Collecting the kinetic and potential variations,
δ g S = d 4 x ( δ K δ U ) ,
we obtain
δ g S = 1 2 d 4 x g { Z ( p ^ ) ( D α p ^ ) ( D β p ^ ) + 1 2 g α β Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) 2 V ( p ^ ) } δ g α β .
Recognizing
δ g α β = g α μ g β ν δ g μ ν ,
we rewrite the expression in the standard form
δ g S = 1 2 d 4 x g I ^ μ ν δ g μ ν ,
where the informational stress tensor is identified as
I ^ μ ν = Z ( p ^ ) ( D μ p ^ ) ( D ν p ^ ) g μ ν 1 2 Z ( p ^ ) g α β ( D α p ^ ) ( D β p ^ ) V ( p ^ ) .
This is the fully derived, operator-level informational energy–momentum tensor.

2.2.5. Symmetry and Hermiticity

Because the action is diffeomorphism invariant and defined with symmetric operator ordering, I ^ μ ν satisfies:
  • Symmetry:  I ^ μ ν = I ^ ν μ .
  • Self-adjointness:  I ^ μ ν = I ^ μ ν .
  • Covariant conservation (proved later):  μ I ^ μ ν is related to the field equation.
These guarantees follow from Noether’s theorem applied to diffeomorphism invariance, but the conservation statement will be established explicitly in Section 2.4.

2.2.6. Classical Limit

Replacing the operator p ^ with its expectation value gives
I μ ν = Z ( p ) μ p ν p g μ ν 1 2 Z ( p ) g α β α p β p V ( p ) ,
which is the informational tensor used in microphysical and cosmological applications.
This completes the fully explicit metric variation. The next Section 2.3 performs the complete field variation, leading to the operator-level informational equation of motion, including all intermediate steps and operator-ordering corrections.

2.3. Variation with Respect to the Informational Field

We now perform the full variation of the action (2) with respect to the operator field p ^ and its adjoint p ^ . Unlike traditional scalar field variations, the present case requires special care due to the presence of:
  • operator-valued functions Z ( p ^ ) and V ( p ^ ) ;
  • non-commutativity between p ^ , Z ( p ^ ) , and derivative operators;
  • the metric-dependent contraction g μ ν ( D μ p ^ ) ( D ν p ^ ) .
To guarantee Hermiticity of the resulting equations, we vary p ^ and p ^ independently and combine the expressions at the end.
Throughout this subsection, boundary terms are explicitly evaluated and discarded under the assumption that variations vanish on the boundary of the integration domain.

2.3.1. Structure of the Variation

From (2), the variation with respect to p ^ is
δ p ^ S = d 4 x g 1 2 δ p ^ Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) δ p ^ V ( p ^ ) .
This expression contains three types of variations:
δ Z ( p ^ ) , δ ( D μ p ^ ) , δ V ( p ^ ) .
We compute each contribution separately.

2.3.2. Variation of the Modulation Function

Since Z ( p ^ ) is an operator function, its variation is
δ p ^ Z ( p ^ ) = Z ( p ^ ) δ p ^ ,
where Z ( p ^ ) is defined by its power-series expansion (Fréchet derivative).
The contribution to the action is therefore
δ p ^ S | Z = 1 2 d 4 x g Z ( p ^ ) δ p ^ g μ ν ( D μ p ^ ) ( D ν p ^ ) .
Because Z ( p ^ ) and ( D μ p ^ ) do not commute in general, we keep this ordering explicit.

2.3.3. Variation of the Covariant Derivative

For a scalar operator field, the covariant derivative simplifies to
D μ p ^ = μ p ^ , so ( D μ p ^ ) = μ p ^ .
Hence
δ p ^ ( D μ p ^ ) = D μ ( δ p ^ ) .
The contribution to the action becomes
δ p ^ S | der = 1 2 d 4 x g Z ( p ^ ) g μ ν D μ ( δ p ^ ) ( D ν p ^ ) .
Integrating by parts,
g Z g μ ν D μ ( δ p ^ ) ( D ν p ^ ) = δ p ^ D μ g Z g μ ν D ν p ^ ,
boundary terms vanish because the variation is assumed to vanish on M .

2.3.4. Variation of the Potential

The variation of the potential term is straightforward:
δ p ^ V ( p ^ ) = V p ^ δ p ^ .

2.3.5. Combined Variation and Euler–Lagrange Equation

Combining all contributions,
δ p ^ S = 1 2 d 4 x g Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) δ p ^ 1 2 d 4 x δ p ^ D μ g Z ( p ^ ) g μ ν D ν p ^ d 4 x g V p ^ δ p ^ .
Factoring out δ p ^ gives
δ p ^ S = d 4 x δ p ^ E [ p ^ ] ,
where the operator Eulerian is
E [ p ^ ] = 1 2 g Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) 1 2 D μ g Z ( p ^ ) g μ ν D ν p ^ g V p ^ .
Requiring δ p ^ S = 0 yields the operator equation of motion:
D μ Z ( p ^ ) g μ ν D ν p ^ 1 2 Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) + V p ^ = 0 .
By Hermiticity, the variation with respect to p ^ yields the conjugate equation.

2.3.6. Covariant Form

Using the identity
D μ ( g X μ ) = g μ X μ ,
the equation can be written in explicitly geometric form:
μ Z ( p ^ ) D μ p ^ 1 2 Z ( p ^ ) ( D μ p ^ ) ( D μ p ^ ) + V p ^ = 0 .

2.3.7. Classical Limit

Replacing the operator field by its expectation value gives
μ Z ( p ) μ p 1 2 Z ( p ) μ p μ p + V ( p ) = 0 .
For the common form
Z ( p ) = 1 1 + β p 2 , Z ( p ) = 2 β p ( 1 + β p 2 ) 2 ,
the equation becomes
μ μ p β p 1 + β p 2 μ p μ p + ( 1 + β p 2 ) V ( p ) = 0 .

2.3.8. Interpretation

The informational field equation exhibits two key features:
  • The principal term μ ( Z D μ p ^ ) defines a modulated wave operator whose magnitude depends on informational density.
  • The nonlinear drift term 1 2 Z ( D μ p ^ ) D μ p ^ has no analogue in canonical scalar field theory and captures the self-interaction of informational inertia.
These two ingredients are responsible for:
  • the emergence of microphysical excitations,
  • the existence of solitonic sectors,
  • the scale-dependent propagation speed,
  • the coupling between informational content and curvature.
This completes the full operator-level variation with respect to p ^ . The next Section 2.4 establishes the fully derived unified Einstein equation, including an explicit proof of covariant conservation.

2.4. Unified Einstein Equation with Informational Coupling

Having obtained the informational tensor I ^ μ ν and the operator equation of motion for p ^ ( x ) , we now derive the full gravitational dynamics of QIR. The key result is the unified Einstein equation, which relates the geometric curvature to the combined informational–matter tensor, modulated by the informational function Z ( p ^ ) .
This subsection provides a complete derivation:
  • variation of the Einstein–Hilbert term,
  • combination with the informational and matter variations,
  • extraction of the field equation,
  • proof of covariant conservation using the Bianchi identity,
  • discussion of the local renormalisation of coupling scales.

2.4.1. Variation of the Einstein–Hilbert Action

The Einstein–Hilbert action is
S EH = 1 16 π G d 4 x g R .
Its variation is classical:
δ S EH = 1 16 π G d 4 x g G μ ν δ g μ ν + μ V μ ,
where V μ is the usual boundary term
V μ = g α β δ Γ α β μ g μ α δ Γ α β β .
Assuming appropriate boundary conditions,
δ S EH = 1 16 π G d 4 x g G μ ν δ g μ ν .

2.4.2. Total Metric Variation of the Full Action

The full gravitational-informational-matter action is
S tot = S EH + S [ p ^ , g ] + S matter .
We have already obtained:
δ S [ p ^ , g ] = 1 2 d 4 x g I ^ μ ν δ g μ ν ,
and
δ S matter = 1 2 d 4 x g T ^ μ ν δ g μ ν ,
where T ^ μ ν is the operator energy–momentum tensor of matter.
Thus the total variation is
δ S tot = d 4 x g 1 16 π G G μ ν + 1 2 I ^ μ ν + T ^ μ ν δ g μ ν .
Requiring stationarity for arbitrary δ g μ ν gives
1 16 π G G μ ν + 1 2 ( I ^ μ ν + T ^ μ ν ) = 0 .
Multiplying both sides by 16 π G yields:
G μ ν = 8 π G ( I ^ μ ν + T ^ μ ν ) .
Next we incorporate the informational renormalisation of the effective speed of light.

2.4.3. Informational Renormalisation and Effective Coupling

In QIR the local speed of light is given by
c ( x ) = c 0 Z ( p ¯ ( x ) ) .
The combination G / c 4 that naturally appears in gravitational theory therefore becomes
G c ( x ) 4 = G c 0 4 Z ( p ¯ ) 2 .
Reabsorbing this factor into the right-hand side of (33) yields the informationally renormalised Einstein equation:
G μ ν = 8 π G c ( x ) 4 Z ( p ^ ) I ^ μ ν + T ^ μ ν .
This is the fundamental operator-level gravitational equation of Quantum Informational Relativity.

2.4.4. Proof of Covariant Conservation

To establish dynamical consistency, we must prove:
μ Z ( p ^ ) I ^ μ ν + T ^ μ ν = 0 .

Step 1: Bianchi identity

The twice-contracted Bianchi identity states:
μ G μ ν = 0 .
Applying μ to (34),
0 = μ G μ ν = μ 8 π G c 4 Z ( p ^ ) ( I ^ μ ν + T ^ μ ν ) .
Because c ( x ) depends only on Z, and Z ( p ^ ) acts on fields but not on the metric, the prefactor behaves as a scalar function under μ :
μ Z c 4 = 1 c 4 μ Z 4 Z c 5 μ c .
But since
c = c 0 Z , μ c = c 0 2 Z μ Z ,
we obtain
μ Z c 4 = 0 .
Thus the prefactor is covariantly constant.

Step 2: Result

Therefore
μ Z ( p ^ ) ( I ^ μ ν + T ^ μ ν ) = 0 ,
which is exactly (35).
This conservation law is the geometric manifestation of the equation of motion for p ^ , as we now show.

2.4.5. Equivalence with the Field Equation

Substituting the explicit form of I ^ μ ν into the conservation law and expanding each term leads (after cancellation of several gradient contributions) precisely to the operator field equation derived in Section 2.3:
μ Z ( p ^ ) D μ p ^ 1 2 Z ( p ^ ) ( D μ p ^ ) D μ p ^ + V p ^ = 0 .
Thus:
  • the Einstein equation enforces the informational field equation,
  • the informational field equation ensures consistency of the Einstein equation.
This mutual implication is analogous to the relationship between matter conservation and the Einstein field equations in standard GR, but now extended to the informational sector.

2.4.6. Classical Limit

The classical limit of the unified Einstein equation is
G μ ν = 8 π G c ( x ) 4 Z ( p ) I μ ν + T μ ν ,
where I μ ν is the classical informational tensor and T μ ν the standard matter tensor.
This form underlies all macroscopic applications of QIR, from cosmological evolution to gravitational lensing and dynamical mass profiles.

2.4.7. Interpretation

The unified Einstein equation exhibits the following structural features:
  • Informational modulation. The factor Z ( p ^ ) acts as a state-dependent modification of gravitational coupling.
  • Local renormalisation. The combination
    8 π G c ( x ) 4 Z ( p ^ )
    effectively defines a local gravitational coupling G eff ( x ) .
  • Self-consistency. Conservation of the combined informational–matter tensor follows automatically from geometric consistency.
  • Operator structure. The equation remains valid at the operator level, with physical predictions recovered via expectation values or classicalisation.
This completes the full derivation of the unified Einstein equation in Quantum Informational Relativity. The next Section 2.5 constructs the covariant canonical momentum, the Hamiltonian, and the ADM (3+1) decomposition required for quantisation and dynamical analysis.

2.5. Conjugate Momentum and Hamiltonian Structure

In this subsection we construct the full canonical structure associated with the informational field p ^ ( x ) . The derivation proceeds in three stages:
  • construction of the covariant conjugate momentum Π μ ,
  • derivation of the covariant Hamiltonian density,
  • ADM ( 3 + 1 ) decomposition and explicit Hamiltonian constraints.
This structure prepares the ground for quantisation in Section 3.

2.5.1. Covariant Conjugate Momentum

The Lagrangian density extracted from (2) is
L = 1 2 Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) V ( p ^ ) .
The covariant conjugate momentum operator is defined as
Π μ L ( D μ p ^ ) = 1 2 Z ( p ^ ) g μ ν ( D ν p ^ ) ,
and similarly
Π μ = 1 2 Z ( p ^ ) g μ ν D ν p ^ .
For clarity we keep the factor 1 2 explicit. When both momenta appear in the Hamiltonian these factors combine naturally.

2.5.2. Covariant Hamiltonian Density

The covariant Hamiltonian density is defined as
H cov = Π μ D μ p ^ + Π μ D μ p ^ L .
Substituting (40),
Π μ D μ p ^ = 1 2 Z g μ ν ( D ν p ^ ) D μ p ^ ,
Π μ D μ p ^ = 1 2 Z g μ ν D ν p ^ ( D μ p ^ ) .
Using symmetric ordering,
Π μ D μ p ^ + ( Π μ ) D μ p ^ = Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) .
Thus the Hamiltonian becomes
H cov = Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) 1 2 Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) V ( p ^ )
= 1 2 Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) + V ( p ^ ) .
Hence:
H cov [ p ^ , g ] = 1 2 Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) + V ( p ^ ) .
This is the covariant energy density of the informational field.

2.5.3. ( 3 + 1 ) ADM Decomposition

We now decompose spacetime into spacelike hypersurfaces Σ t with induced metric h i j , lapse N, and shift N i . The line element is
d s 2 = N 2 d t 2 + h i j ( d x i + N i d t ) ( d x j + N j d t ) .
The inverse metric is
g μ ν = 1 / N 2 N j / N 2 N i / N 2 h i j N i N j N 2 .
For a scalar operator field,
D μ p ^ = μ p ^ .
Thus:
D 0 p ^ = p ^ ˙ N i i p ^ , D i p ^ = i p ^ .
Using the metric, the kinetic structure decomposes as:
g μ ν D μ p ^ ( D ν p ^ ) = 1 N 2 ( p ^ ˙ N i i p ^ ) ( p ^ ˙ N i i p ^ )
+ h i j N i N j N 2 i p ^ j p ^ .

2.5.4. Canonical Momentum in ADM Variables

Using (40):
Π 0 = Z 2 N 2 ( p ^ ˙ N i i p ^ ) , Π i = Z h i j j p ^ Z N i 2 N 2 ( p ^ ˙ N k k p ^ ) .
The momentum conjugate to p ^ on the hypersurface Σ t is
Π L p ^ ˙ = Z 2 N ( p ^ ˙ N i i p ^ ) .
Similarly
Π = Z 2 N ( p ^ ˙ N i i p ^ ) .

2.5.5. ADM Hamiltonian Density

The ADM Hamiltonian density is defined as
H ADM = Π p ^ ˙ + Π p ^ ˙ L .
Using
p ^ ˙ = 2 N Z Π + N i i p ^ ,
we obtain after substitution and simplification:
H ADM = N 2 Z Π Π + 1 2 Z h i j i p ^ j p ^ + V ( p ^ ) + N i H i ,
where
H i = Π i p ^ + Π i p ^
is the momentum constraint density.

2.5.6. Hamiltonian and Momentum Constraints

The Hamiltonian constraint is obtained from the coefficient of the lapse:
H = 2 Z ( p ^ ) Π Π + 1 2 Z ( p ^ ) h i j i p ^ j p ^ + V ( p ^ ) = 0 .
The momentum constraint is obtained from the coefficient of the shift:
H i = Π i p ^ + Π i p ^ = 0 .
These constraints are the generators of diffeomorphisms within the hypersurface Σ t and normal to it.

2.5.7. Classical Limit

Replacing operators by their expectation values gives the classical Hamiltonian:
H ADM cl = N 2 Z ( p ) Π 2 + 1 2 Z ( p ) h i j i p j p + V ( p ) + N i Π i p .
The informational factors Z and 1 / Z appear naturally:
- Z multiplies the spatial gradient term (inertia of geometry), - 1 / Z multiplies the kinetic term (inertia of the informational excitation).

2.5.8. Interpretation

The ADM Hamiltonian displays two characteristic signatures of QIR:
  • Dual inertia structure. The combination
    2 Z Π Π
    shows that excitations become lighter or heavier depending on informational density.
  • Geometric–informational coupling. The factor Z h i j i p ^ j p ^ links informational gradients directly to spatial curvature.
  • Constraint structure identical to GR. Hamiltonian and momentum constraints retain their geometric meaning, guaranteeing full diffeomorphism invariance.
This completes the construction of the canonical and Hamiltonian structures of the informational field. The next Section 2.6 performs the expansion of the action to second order around a background solution p 0 and derives the canonical normalisation of fluctuations, explaining the appearance of Z 1 in microphysics.

2.6. Canonical Normalization and Second-Order Expansion

To prepare for the quantisation procedure in Section 3, we must expand the informational action to second order in small fluctuations around a background configuration. This analysis yields the canonical normalisation of the excitations and explains why microphysical observables naturally acquire the factor Z 1 ( p 0 ) .
We proceed systematically:
  • define the background and fluctuation fields,
  • expand the kinetic term to quadratic order,
  • expand the potential term consistently,
  • extract the quadratic operator and diagonalise it,
  • identify the canonically normalised field and its effective mass.
Throughout this section we work in the classical limit; the operator case follows by standard ordering rules.

2.6.1. Background–Fluctuation Split

Let p 0 ( x ) be a solution of the classical field equation derived in Section 2.3. We introduce small fluctuations:
p ( x ) = p 0 ( x ) + δ p ( x ) , δ p p 0 .
The modulation function expands as
Z ( p ) = Z 0 + Z 0 δ p + 1 2 Z 0 ( δ p ) 2 + O ( ( δ p ) 3 ) ,
with Z 0 Z ( p 0 ) .

2.6.2. Expansion of the Kinetic Term

The kinetic part of the action is
S kin = 1 2 d 4 x g Z ( p ) g μ ν μ p ν p .
Insert the decomposition (53). We first expand the derivative:
μ p = μ p 0 + μ δ p .
Then:
g μ ν μ p ν p = g μ ν μ p 0 + μ δ p ν p 0 + ν δ p
= g μ ν μ p 0 ν p 0 + 2 g μ ν μ p 0 ν δ p + g μ ν μ δ p ν δ p .
Now multiply by Z ( p ) expanded to second order:
Z ( p ) = Z 0 + Z 0 δ p + 1 2 Z 0 ( δ p ) 2 .
Keeping all terms up to quadratic order in δ p :
Z g μ ν μ p ν p = Z 0 g μ ν μ p 0 ν p 0
+ 2 Z 0 g μ ν μ p 0 ν δ p + Z 0 g μ ν μ δ p ν δ p
+ Z 0 δ p g μ ν μ p 0 ν p 0 + 2 Z 0 δ p g μ ν μ p 0 ν δ p
+ 1 2 Z 0 ( δ p ) 2 g μ ν μ p 0 ν p 0 + O ( ( δ p ) 3 ) .
The linear terms vanish upon integration because p 0 satisfies the background equation of motion. Thus, the quadratic part of the kinetic action is:
S kin ( 2 ) = 1 2 d 4 x g Z 0 g μ ν μ δ p ν δ p + 2 Z 0 δ p g μ ν μ p 0 ν δ p + 1 2 Z 0 ( δ p ) 2 g μ ν μ p 0 ν p 0 .

2.6.3. Removing Linear Derivative Mixing

The mixed term
2 Z 0 δ p g μ ν μ p 0 ν δ p
can be integrated by parts:
g δ p g μ ν μ p 0 ν δ p = 1 2 g δ p 2 μ g μ ν ν p 0 ,
boundary terms vanish.
Thus, the mixed derivative term produces only a mass-like correction. The quadratic kinetic term becomes
S kin ( 2 ) = 1 2 d 4 x g Z 0 g μ ν μ δ p ν δ p + M mix 2 ( δ p ) 2 ,
where
M mix 2 = Z 0 1 2 g μ ν μ p 0 ν p 0 Z 0 μ μ p 0 .

2.6.4. Expansion of the Potential

The potential expands as:
V ( p ) = V 0 + V 0 δ p + 1 2 V 0 ( δ p ) 2 .
The linear term vanishes because p 0 satisfies the background equation. Thus:
S V ( 2 ) = 1 2 d 4 x g V 0 ( δ p ) 2 .

2.6.5 Complete Quadratic Action

Combining kinetic and potential contributions,
S ( 2 ) = 1 2 d 4 x g Z 0 g μ ν μ δ p ν δ p m eff 2 ( δ p ) 2 ,
where
m eff 2 = V 0 M mix 2 .
This is the quadratic Lagrangian for fluctuations.

2.6.6. Canonical Normalization: Emergence of Z 1

The kinetic term has a prefactor Z 0 :
Z 0 g μ ν μ δ p ν δ p .
To obtain a canonically normalised field, we define
χ ( x ) Z 0 δ p ( x ) .
Then
μ δ p = 1 Z 0 μ χ ,
and the kinetic term becomes
Z 0 g μ ν 1 Z 0 μ χ 1 Z 0 ν χ = g μ ν μ χ ν χ .
Thus the canonically normalised quadratic action is
S ( 2 ) = 1 2 d 4 x g g μ ν μ χ ν χ m χ 2 χ 2 ,
with effective mass
m χ 2 = m eff 2 Z 0 .

2.6.7 Physical Interpretation

The renormalised mass demonstrates explicitly:
m micro 2 ( x ) = V ( p 0 ) Z ( p 0 ) .
Thus:
- In regions where Z ( p 0 ) > 1 , excitations are lighter. - In regions where Z ( p 0 ) < 1 , excitations are heavier. - All microphysical scales are rescaled by Z 1 / 2 or Z 1 .
This factor is the mathematical origin of the renormalisation of:
  • masses,
  • couplings,
  • propagation speeds,
  • dispersion relations,
in the informational framework.
This completes the full quadratic expansion and the rigorous derivation of the canonical normalisation. The next Section 2.7 develops the explicit computation of effective constants and propagation speeds, completing the foundational operator structure.

2.7. Effective Constants and Propagation Speed

Having established the canonical structure and quadratic expansion of the informational field, we now derive how informational modulation renormalizes the fundamental constants relevant for both the macroscopic geometric sector and the microscopic excitation sector.
We proceed via a sequence of rigorous derivations:
  • extraction of the propagation speed of fluctuations from the quadratic action,
  • computation of the effective Newton constant G eff ( x ) ,
  • renormalisation of gauge couplings and other interaction strengths,
  • consistency checks linking operator-level structure to canonical normalization.
The results generalize the familiar notion of running couplings to a state-dependent and geometrically integrated scale factor encoded in Z ( p ) .

2.7.1. Informational Origin of the Local Propagation Speed

In Section 2.6 we derived the quadratic action for the canonically normalized field χ :
S ( 2 ) = 1 2 d 4 x g g μ ν μ χ ν χ m χ 2 χ 2 .
To extract the propagation speed, we return to the un-normalized fluctuations δ p :
S ( 2 ) = 1 2 d 4 x g Z 0 g μ ν μ δ p ν δ p m eff 2 ( δ p ) 2 .
Consider a local inertial frame where
g μ ν = diag ( 1 , 1 , 1 , 1 ) .
The wave operator acting on δ p is
Z 0 t 2 + 2 δ p .
To bring the time derivative into the canonical form t 2 , we observe that the physical propagation speed is:
c δ p 2 = Z 0 .
Thus,
c ( x ) = c 0 Z ( p 0 ( x ) ) .
This matches precisely the kinematic identity introduced earlier, now derived rigorously from the quadratic operator structure.

2.7.2. Effective Newton Constant

The unified Einstein equation derived in Section 2.4 is:
G μ ν = 8 π G c ( x ) 4 Z ( p ) I μ ν + T μ ν .
Substituting c ( x ) = c 0 Z gives:
Z c ( x ) 4 = Z c 0 4 Z 2 = 1 c 0 4 Z .
Thus the effective Newton constant is:
G eff ( x ) G Z ( p ( x ) ) .
G eff ( x ) = G Z ( p ( x ) ) .
Hence: - in high-informational-density regions ( Z > 1 ), gravity is weaker; - in low-informational-density regions ( Z < 1 ), gravity is stronger.
This effect plays a central role in large-scale structure and lensing observables.

2.7.3. Gauge Coupling Renormalisation

Consider a minimally coupled gauge field with canonical action
S gauge = 1 4 g 0 2 d 4 x g F μ ν F μ ν .
The physical field couples to δ p through the kinetic normalization of δ p χ :
δ p = 1 Z 0 χ .
If the interaction term is
L int = λ 0 δ p J ,
then after canonical normalization:
L int = λ 0 Z 0 χ J .
Thus the effective coupling is
g eff ( x ) = g 0 Z ( p ( x ) ) 1 / 2 .
General interactions of the form δ p n acquire factors Z n / 2 . This is consistent with the dimensional analysis in the quadratic expansion.

2.7.4. Effective Masses

From Section 2.6, the canonically normalised mass is:
m χ 2 = m eff 2 Z 0 = V ( p 0 ) Z 0 .
Thus the microphysical mass renormalisation is:
m micro ( x ) = 1 Z ( p ( x ) ) V ( p ( x ) ) .
Since Z appears in the denominator: - excitations are lighter at high informational density, - excitations are heavier at low informational density.

2.7.5 Summary Table

We summarize the derived effective constants:
Table 1.
Quantity Effective Value
Propagation speed c ( x ) = c 0 Z ( p ( x ) )
Newton constant G eff ( x ) = G / Z ( p ( x ) )
Gauge couplings g eff ( x ) = g 0 Z ( p ( x ) )
Masses m micro ( x ) = V ( p ) Z ( p )
All renormalisations derive from one universal function Z ( p ) , expressing the fact that informational density reshapes both microphysical and macroscopic interactions in a coherent, self-consistent manner.

2.7.6. Interpretation

The following points summarize the physical significance:
  • Unification. Gravity, gauge couplings, masses, and propagation speeds are modulated by the same informational factor.
  • Scale dependence. Unlike RG running, the modulation depends on the local state p ( x ) rather than energy scale.
  • Geometry–information interplay. Because G eff 1 / Z and c Z , informationally dense regions have:
    • weaker gravity,
    • faster propagation speeds,
    • lighter excitations.
  • Consistency. The effective constants match exactly the renormalisation extracted from:
    • the full operator EOM (Section 2.3),
    • the Einstein equation (Section 2.4),
    • canonical normalisation (Section 2.6).
This concludes the derivation of the effective constants and propagation structure. The next and final Section 2.8 provides a synthetic summary of all foundational results, preparing the transition to the quantisation procedures developed in Section 3.

2.8. Foundational Synthesis

This subsection summarizes the mathematical and conceptual results established throughout Section 2. These results collectively define the complete foundational structure of Quantum Informational Relativity (QIR) at the analytical, variational, and canonical levels.
The logical progression consisted of four pillars:
  • the operator-level definition of the action,
  • the derivation of the coupled field equations,
  • the construction of the Hamiltonian and constraint structure,
  • the extraction of effective constants and propagation laws.
We now synthesize these components.

(1) Operator Action and Core Structure

The fundamental dynamics are encoded in the operator action
S [ p ^ , g ] = 1 2 d 4 x g Z ( p ^ ) g μ ν ( D μ p ^ ) ( D ν p ^ ) d 4 x g V ( p ^ ) .
All informational effects arise from the modulation function Z ( p ^ ) , which:
  • rescales propagation and inertia,
  • alters effective couplings,
  • appears directly in the Einstein equation,
  • determines the canonical normalization of fluctuations.
At this stage, QIR already differs from standard scalar field theories through its operator dependence and the non-trivial derivative structure.

(2) Coupled Field Equations

The informational field obeys the operator Euler–Lagrange equation
μ Z ( p ^ ) D μ p ^ 1 2 Z ( p ^ ) ( D μ p ^ ) D μ p ^ + V p ^ = 0 .
Variation with respect to the metric generates the informational tensor I ^ μ ν and the unified Einstein equation:
G μ ν = 8 π G c ( x ) 4 Z ( p ^ ) I ^ μ ν + T ^ μ ν .
Covariant conservation,
μ Z ( p ^ ) I ^ μ ν + T ^ μ ν = 0 ,
follows automatically from the Bianchi identity, establishing full self-consistency.
This shows that the informational field equation and the Einstein equation mutually enforce each other.

(3) Canonical Momentum and Hamiltonian Structure

The conjugate momentum is
Π = Z 2 N ( p ^ ˙ N i i p ^ ) ,
and the ADM Hamiltonian density is
H ADM = N 2 Z Π Π + 1 2 Z h i j i p ^ j p ^ + V ( p ^ ) + N i H i .
The Hamiltonian and momentum constraints are:
H = 0 , H i = 0 ,
ensuring full diffeomorphism invariance.
The ADM formalism reveals the dual informational structure:
kinetic inertia Z 1 , spatial inertia Z .
This duality plays a central role in the emerging microphysics.

(4) Quadratic Expansion and Canonical Normalization

Expanding around a background p 0 yields the quadratic action
S ( 2 ) = 1 2 d 4 x g Z 0 g μ ν μ δ p ν δ p m eff 2 ( δ p ) 2 ,
where Z 0 Z ( p 0 ) .
The canonically normalized field is
χ = Z 0 δ p .
The normalized mass is
m χ 2 = m eff 2 Z 0 .

(5) Effective Constants

The informational framework yields the effective constants:
c ( x ) = c 0 Z ( p ( x ) ) , G eff ( x ) = G Z ( p ( x ) ) , g eff ( x ) = g 0 Z ( p ( x ) ) , m micro ( x ) = V Z .
These follow consistently from:
  • the Einstein equation,
  • the quadratic expansion,
  • canonical normalization,
  • gauge-field couplings.
All renormalised constants originate from the single informational function Z.

(6) Outlook Toward Quantisation

The results of Section 2 provide the full classical and operator-level foundation necessary for quantisation:
  • the operator action is fully specified,
  • the Euler–Lagrange equation is established,
  • the Hamiltonian and momentum constraints are known,
  • the canonical field χ is defined,
  • the quadratic operator governing quantum fluctuations is explicit.
These ingredients permit the construction of:
  • the quantum Hamiltonian operator,
  • the mode expansion on generic backgrounds,
  • the propagator,
  • the spectral decomposition of informational excitations.
This sets the stage for Section 3, which develops the quantisation of the informational field and its resulting microphysical spectrum.

3. Quantization of the Informational Field

In this section we develop the full quantum theory of informational excitations. Building on the classical and operator foundations established in Section 2, we construct the canonical quantization, mode structure, propagators, spectral decomposition, and renormalized quantum observables relevant for microphysical and astrophysical predictions.

3.1. Operator Canonical Quantization

We now quantize the informational field on a curved background ( M , g μ ν ) endowed with a classical informational configuration p 0 ( x ) . The quadratic expansion derived in Section 2.6 provides the starting point:
S ( 2 ) = 1 2 d 4 x g Z 0 g μ ν μ δ p ν δ p m eff 2 ( δ p ) 2 ,
with Z 0 = Z ( p 0 ( x ) ) .
This structure is not canonically normalized due to the presence of Z 0 in front of the kinetic term. We therefore begin with the canonical redefinition of the fluctuation field.

3.1.1. Canonical Field Redefinition

Following Section 2.6, the canonically normalized field is defined as
χ ( x ) Z 0 ( x ) δ p ( x ) .
The action becomes
S ( 2 ) = 1 2 d 4 x g g μ ν μ χ ν χ m χ 2 ( x ) χ 2 ,
where the renormalized mass is
m χ 2 ( x ) = m eff 2 ( x ) Z 0 ( x ) = V ( p 0 ( x ) ) Z 0 ( x ) M mix 2 ( x ) Z 0 ( x ) .
The canonical normalization implies that the Hamiltonian derived from (71) has the standard kinetic term
χ ˙ 2 ,
and therefore the quantization rules take their usual form. The informational structure is now embedded entirely in the background-dependent mass and mode functions.

3.1.2. Hilbert Space and Operator Algebra

We quantize χ ( x ) as a scalar quantum field on the curved background. The Hilbert space H is generated by creation and annihilation operators associated with a complete set of mode functions { u n ( x ) } constructed in Section 3.2.
The operator field is promoted to an operator-valued distribution
χ ^ : H H .
The canonical momentum conjugate to χ is
Π χ ( x ) = L ( 2 ) ( t χ ) = g g 0 μ μ χ .
On a general ADM background,
Π χ = h N χ ˙ N i i χ .

3.1.3. Canonical Commutation Relations

Equal-time canonical commutation relations follow from the canonical structure of (71):
χ ^ ( t , x ) , χ ^ ( t , y ) = 0 ,
Π ^ χ ( t , x ) , Π ^ χ ( t , y ) = 0 ,
χ ^ ( t , x ) , Π ^ χ ( t , y ) = i δ ( 3 ) ( x y ) .
This structure is identical to the canonical quantization of a Klein–Gordon field, but the informational effects appear:
  • indirectly through the mass m χ 2 ( x ) ,
  • directly in the mode functions through Z 0 ( x ) ,
  • through the background geometry determined by the Einstein equation with informational coupling.

3.1.4. Field Operator in Terms of Modes

The operator field admits the mode expansion
χ ^ ( x ) = n a ^ n u n ( x ) + a ^ n u n * ( x ) ,
where:
  • u n ( x ) are c-number solutions of the mode equation derived in Section 3.2,
  • a ^ n , a ^ n satisfy the standard algebra
    [ a ^ n , a ^ m ] = δ n m , [ a ^ n , a ^ m ] = 0 ,
  • the normalization of u n is fixed by the Klein–Gordon inner product.
The presence of Z 0 ( x ) means that the mode equation is not the standard Klein–Gordon form unless Z 0 is constant; this is the origin of the state-dependent microphysical spectrum.

3.1.5. Transformation from the Noncanonical Field

Since experiments and phenomenology are phrased in terms of fluctuations δ p rather than χ , we provide the inverse transformation:
δ p ( x ) = 1 Z 0 ( x ) χ ^ ( x ) .
Thus the operator fluctuations of the informational field satisfy
δ p ^ ( x ) = n u n ( x ) Z 0 ( x ) a ^ n + u n * ( x ) Z 0 ( x ) a ^ n .
The renormalized factor Z 0 1 / 2 is responsible for:
  • amplitude suppression or enhancement in high-/low-informational density regions,
  • modified couplings and decay amplitudes,
  • an effective “local normalization scale” for all quantum processes.

3.1.6. Interpretation

Canonical quantization of the informational field yields a quantum theory with:
  • Standard commutation relations, i.e. the informational structure does not alter the fundamental algebra of quantum fields.
  • Nonstandard mode functions, since Z 0 ( x ) modifies the effective mass and propagation operator.
  • Renormalized excitations, with physical mass
    m micro ( x ) = V ( p 0 ( x ) ) Z 0 ( x ) ,
    matching the effective constants derived in Section 2.7.
  • Local informational dressing, manifest in amplitudes, propagators, and correlation functions.
This completes the canonical operator quantization. The next Section 3.2 constructs the corresponding mode functions on an arbitrary curved background, establishing orthonormality, completeness, and the Wronskian condition.

3.2. Mode Expansion on Curved Backgrounds

Having established the canonical quantization of the informational field in Section 3.1, we now construct the mode functions that define the operator expansion of χ ^ ( x ) . Because the underlying geometry is an arbitrary curved spacetime ( M , g μ ν ) , the mode structure generalizes the standard Klein–Gordon decomposition and incorporates the informational modulation through the background-dependent mass m χ 2 ( x ) .
The goal of this subsection is to build the complete basis of solutions { u n ( x ) } of the mode equation, define the appropriate normalization using the Klein–Gordon inner product, and establish orthonormality and completeness relations.

3.2.1. Mode Equation on a Curved Background

Starting from the quadratic action
S ( 2 ) = 1 2 d 4 x g g μ ν μ χ ν χ m χ 2 ( x ) χ 2 ,
the Euler–Lagrange equation yields the generalized Klein–Gordon equation:
μ μ + m χ 2 ( x ) χ ( x ) = 0 .
Because m χ 2 ( x ) depends on Z 0 ( x ) , the informational structure directly affects the propagation of modes. We seek separable solutions of the form
χ ( x ) = u n ( x ) ,
which satisfy
μ μ u n ( x ) + m χ 2 ( x ) u n ( x ) = 0 .
These functions will serve as the mode basis for the operator expansion.

3.2.2. Klein–Gordon Inner Product

The appropriate scalar product for Klein–Gordon fields on curved spacetime is
( u , v ) KG i Σ t d Σ μ u * μ v v μ u * ,
where Σ t is an arbitrary spacelike Cauchy surface with unit normal n μ and induced volume element d Σ = h d 3 x .
Using d Σ μ = n μ d Σ , we write:
( u , v ) KG = i Σ t d 3 x h u * n μ μ v v n μ μ u * .
This inner product is conserved:
t ( u , v ) KG = 0 ,
a consequence of the mode Equation (82) and the metric covariance of the theory.

3.2.3. Orthonormality and Wronskian Condition

We impose the orthonormality relations
( u n , u m ) KG = δ n m ,
( u n * , u m * ) KG = δ n m ,
( u n , u m * ) KG = 0 .
These conditions ensure that the creation and annihilation operators associated with the mode expansion satisfy
[ a ^ n , a ^ m ] = δ n m .
To express orthonormality in local form, we write the Wronskian condition. Choosing coordinates where the ADM normal is n μ = ( N , 0 , 0 , 0 ) , we obtain:
u n * t u m u m t u n * = i N h δ n m .
This condition generalizes the flat-spacetime Wronskian u n * u ˙ m u ˙ n * u m = i to arbitrary geometries. The informational structure enters implicitly through m χ 2 ( x ) , which modifies the time evolution of the modes and therefore their Wronskian.

3.2.4. Completeness of the Mode Basis

The modes form a complete basis in the sense that any solution χ ( x ) of the field Equation (80) may be expanded as:
χ ( x ) = n A n u n ( x ) + B n u n * ( x ) ,
where the coefficients are uniquely determined by the inner product (83).
Completeness also implies the closure relation:
n u n ( x ) u n * ( x ) u n * ( x ) u n ( x ) = i Δ ( x , x ) ,
where Δ ( x , x ) is the Pauli–Jordan function (commutator function). This object will be essential for constructing the Feynman propagator in Section 3.3.

3.2.5. Structure of the Mode Functions

The mode equation
μ μ u n + m χ 2 ( x ) u n = 0
is significantly modified by the informational field because:
  • the mass term includes Z 0 ( x ) ,
  • the geometry g μ ν ( x ) depends on p ^ through the unified Einstein equation,
  • time evolution of the modes is therefore sensitive to informational density.
In regions where Z 0 varies slowly (adiabatic approximation), modes behave locally like:
u n ( x ) 1 2 ω n ( x ) exp i t ω n ( t ) d t ,
with
ω n 2 ( x ) = k n 2 a ( t ) 2 + m χ 2 ( x ) + O ( Z ˙ 0 ) .
In regions of strong informational gradients, the modes acquire nontrivial frequency mixing, leading to:
  • enhanced particle production,
  • state-dependent propagation,
  • potentially observable signatures in high-energy processes.

3.2.6. Interpretation

The mode decomposition highlights three central features of the quantum informational field:
  • The informational field modifies the spectrum. The mass m χ 2 ( x ) depends on Z 0 , making the spectrum environment-dependent.
  • The orthonormality and completeness structure is preserved. Despite the informational modification, the KG inner product yields the standard commutation algebra.
  • The modes encode all informational effects. Once the mode basis is known, all quantum observables propagators, correlation functions, decay rates—follow directly from the operator decomposition.
This completes the construction of the mode functions. The next Section 3.3 builds the Feynman propagator and the two-point function using the mode sum representation established here.

3.3. Propagator and Two-Point Function

The quantized informational field χ ^ ( x ) admits a complete mode expansion (Section 3.2), which allows us to construct the full set of two-point functions and quantum propagators. These objects encode all microscopic predictions of QIR: particle production, dispersion, spectral density, decay amplitudes, and the observationally relevant correlation functions entering SWIFT and cosmological analyses.
We proceed systematically:
  • define the Green equation on curved spacetime,
  • construct the Wightman and Feynman propagators from the mode basis,
  • analyze the effect of informational modulation Z 0 ( x ) ,
  • give explicit forms in limiting geometries.

3.3.1. Green Function Equation

The propagator is the Green function of the generalized Klein–Gordon operator:
μ μ + m χ 2 ( x ) G ( x , x ) = δ ( 4 ) ( x x ) g .
Because m χ 2 ( x ) depends on Z 0 ( x ) , the informational field modifies both the amplitude and the phase of the propagator. Equation (89) defines a family of propagators, depending on boundary conditions:
G + ( x , x ) , G ( x , x ) , G F ( x , x ) , G ret ( x , x ) , G adv ( x , x ) .

3.3.2. Wightman Functions

The positive- and negative-frequency Wightman functions are defined by
G + ( x , x ) = 0 | χ ^ ( x ) χ ^ ( x ) | 0 ,
G ( x , x ) = 0 | χ ^ ( x ) χ ^ ( x ) | 0 .
Using the mode expansion
χ ^ ( x ) = n a ^ n u n ( x ) + a ^ n u n * ( x ) ,
and the vacuum condition a ^ n | 0 = 0 , we find
G + ( x , x ) = n u n ( x ) u n * ( x ) ,
G ( x , x ) = n u n * ( x ) u n ( x ) .
Notice that because u n solves the informationally modified mode equation, both G + and G inherit the dependence on Z 0 ( x ) through m χ ( x ) .

3.3.3. Pauli–Jordan Function

The commutator function is
Δ ( x , x ) G + ( x , x ) G ( x , x ) = n u n ( x ) u n * ( x ) u n * ( x ) u n ( x ) .
Using the closure relation (88), this expression matches precisely the covariant Pauli–Jordan function that implements microcausality:
[ χ ^ ( x ) , χ ^ ( x ) ] = i Δ ( x , x ) .
Microcausality is therefore preserved in QIR, despite the informational modulation of the mass and propagation operator.

3.3.4. Feynman Propagator

The Feynman propagator is defined by
G F ( x , x ) = 0 | T { χ ^ ( x ) χ ^ ( x ) } | 0 ,
where T denotes time ordering.
Using the Wightman functions, we obtain
G F ( x , x ) = Θ ( t t ) G + ( x , x ) + Θ ( t t ) G ( x , x ) .
Equivalently, in mode-sum form:
G F ( x , x ) = n Θ ( t t ) u n ( x ) u n * ( x ) + Θ ( t t ) u n * ( x ) u n ( x ) .
The propagator satisfies the Green Equation (89) with Feynman boundary conditions. Its analytic structure encodes particle propagation, decay processes, and radiative corrections, all of which inherit the informational factor Z 0 ( x ) .

3.3.5. Retarded and Advanced Propagators

The retarded and advanced Green functions are defined by
G ret ( x , x ) = Θ ( t t ) Δ ( x , x ) ,
G adv ( x , x ) = Θ ( t t ) Δ ( x , x ) .
They satisfy:
μ μ + m χ 2 ( x ) G ret / adv ( x , x ) = δ ( 4 ) ( x x ) g .
Because the commutator Δ ( x , x ) is not affected by time-ordering prescriptions, the causal structure of QIR is standard:
Δ ( x , x ) = 0 for spacelike separation .
Informational physics preserves causality.

3.3.6. Special Limits

(1) Minkowski spacetime.

For constant Z 0 and constant mass m χ :
u k ( x ) = 1 2 ω k e i ( ω k t k · x ) , ω k = k 2 + m χ 2 .
Thus
G F ( x x ) = d 4 k ( 2 π ) 4 i e i k · ( x x ) k 2 m χ 2 + i ϵ .
All informational effects enter through
m χ 2 = V ( p 0 ) Z 0 .

(2) FRW spacetime.

Let
d s 2 = d t 2 + a 2 ( t ) d x 2 .
Defining conformal time η , the mode equation becomes
u k ( η ) + k 2 + a 2 m χ 2 a a u k ( η ) = 0 .
The propagator is then
G F ( x , x ) = d 3 k ( 2 π ) 3 u k ( η ) u k * ( η ) e i k · ( x x ) .
This form is used extensively in QIR cosmology (Section 5).

3.3.7. Informational Effects on the Propagator

The factor Z 0 ( x ) affects the propagator in two main ways:
  • Modified mass.
    m χ 2 ( x ) = m eff 2 ( x ) Z 0 ( x )
    influences:
    • the phase e i ω n d t ,
    • the UV behavior,
    • decay rates,
    • threshold energies.
  • Modified amplitude. In (97), the normalization of u n carries a factor 1 / 2 ω n , thus:
    G F ω n 1 m χ 1 .
    Informationally dense regions ( Z 0 > 1 ) enhance correlations, informationally dilute regions ( Z 0 < 1 ) suppress them.
These effects will play a central role in Sections 4–6, especially in high-energy processes and SWIFT photon propagation.

3.3.8. Interpretation and Outlook

The propagator encodes the full quantum dynamics of the informational field:
  • It is determined entirely by the informationally modified mode functions.
  • It inherits all state dependence through the background Z 0 ( x ) .
  • Its causal structure is standard, guaranteeing consistency with relativity.
With the propagator established, we can now analyze the microphysical spectrum and the decomposition into elementary excitations, which is the subject of Section 3.4.

3.4. Spectral Decomposition and Microphysical States

The propagator constructed in Section 3.3 contains the full quantum information about the microphysical excitations of the informational field. We now extract the spectral decomposition, define the physical mass eigenstates, and identify the microphysical modes relevant for both high-energy phenomenology and astrophysical propagation.
The informational theory exhibits two sources of spectral structure:
  • geometric curvature effects,
  • informational modulation through Z 0 ( x ) .
Their interplay determines the mass, dispersion relation, and stability properties of informational excitations.

3.4.1. Spectral Problem Associated with the Propagator

Consider the Fourier–Laplace transform of the Feynman propagator:
G F ( x , x ) = n + d ω 2 π u n , ω ( x ) u n , ω * ( x ) ω 2 ω n 2 + i ϵ .
The poles of the propagator satisfy
ω 2 = ω n 2 ,
with
ω n 2 ( x ) = E n 2 ( x ) = k n 2 + m χ 2 ( x ) + Δ geom ( x ) + Δ info ( x ) ,
where:
  • k n are the comoving or physical momenta (depending on the decomposition),
  • m χ 2 ( x ) is the informationally renormalized mass,
  • Δ geom are curvature-induced terms,
  • Δ info arise from spatial gradients of Z 0 ( x ) .
The physical excitations correspond to the pole structure of the propagator.

3.4.2. Physical Mass from the Spectral Pole

The pole of the Feynman propagator is at
k 2 = m micro 2 ( x ) ,
with the informationally renormalized mass
m micro 2 ( x ) = m χ 2 ( x ) = m eff 2 ( x ) Z 0 ( x ) = V ( p 0 ( x ) ) Z 0 ( x ) M mix 2 ( x ) Z 0 ( x ) .
Thus the microphysical mass inherits a full state dependence: - as Z 0 ( x ) increases (high informational density),
m micro ( x )
- as Z 0 ( x ) decreases (informational dilution),
m micro ( x ) .
This mass defines the “elementary excitation” of QIR.

3.4.3. Spectral Decomposition of the Propagator

Using the relations
G F ( x , x ) = 0 | T { χ ^ ( x ) χ ^ ( x ) } | 0 ,
and the mode expansion, the spectral representation is
G F ( x , x ) = 0 d μ 2 ρ ( μ 2 ; x , x ) G F ( μ ) ( x , x ) ,
with:
- ρ ( μ 2 ; x , x ) : densité spectrale (spectral density), - G F ( μ ) : propagateur d’un champ de masse μ .
For a free informational field (quadratic theory), the spectral density is sharp:
ρ ( μ 2 ; x , x ) = δ μ 2 m micro 2 ( x ) .
Thus:
G F ( x , x ) = G F ( m micro ) ( x , x ) .
The informational field produces a *single, sharp* microphysical excitation, analogous to a massive Klein–Gordon field, but with state-dependent mass.

3.4.4. Dispersion Relation and Microphysical Modes

Consider a local inertial frame. The general solution of the mode equation yields the dispersion relation
ω 2 = k 2 + m micro 2 ( x ) + Δ geom ( x ) + Δ info ( x ) .
The corrections are:
  • Δ geom from curvature, e.g. a / a in FRW,
  • Δ info from gradients of Z 0 ( x ) ,
where the latter equals
Δ info ( x ) = 1 2 Z 0 Z 0 + 3 4 μ Z 0 μ Z 0 Z 0 2 + .
Thus informational inhomogeneities contribute effective potential terms.

3.4.5. Stability Analysis

The condition for stability of microphysical excitations is:
m micro 2 ( x ) + Δ geom ( x ) + Δ info ( x ) > 0 .
Because
m micro 2 ( x ) = V ( p 0 ) Z 0 ,
a large positive Z 0 tends to stabilize the theory. Regions of low Z 0 require V > 0 to avoid tachyonic instabilities.
The informational theory is therefore stable provided:
Z 0 ( x ) > 0 , V ( p 0 ) > 0 .
These conditions match the classical stability analysis of Section 2.

3.4.6. Microphysical States and Observables

The elementary excitations are created by the operators
a ^ n | 0 ,
and correspond to informational quasi-particles with mass m micro ( x ) .
The key microphysical observables include:
  • Energy of excitations:
    E n ( x ) = ω n ( x ) .
  • Spectral line:
    ρ ( μ 2 ) = δ ( μ 2 m micro 2 ( x ) ) .
  • Propagation amplitude: governed by G F ( x , x ) .
  • Decay/transition rates: modulated by 1 / Z 0 1 / 2 in amplitudes.
These states form the basis for high-energy processes in Section 4.

3.4.7. Interpretation and Outlook

The spectral decomposition reveals that:
  • QIR predicts a single, sharp elementary excitation (at quadratic level),
  • its physical mass is state-dependent due to the factor Z 0 ( x ) ,
  • dispersion and propagation inherit geometric and informational corrections,
  • the microphysical spectrum connects directly to:
    • high-energy tests (Section 4),
    • cosmological propagation (Section 5),
    • SWIFT photon timing and hardness correlations (Section 6).
With the spectral structure understood, we now examine how informational modulation affects quantum observables, scattering amplitudes, and decay rates in Section 3.5.

3.5. Informational Renormalization of Quantum Observables

Quantum observables—transition amplitudes, decay rates, and scattering cross sections—are sensitive to the normalization of the quantum field. In QIR, the canonical field is
χ = Z 0 δ p ,
so the informational structure induces a nontrivial rescaling of:
  • field amplitudes,
  • interaction vertices,
  • propagators,
  • spectral densities.
This subsection derives the complete informational renormalization pattern and compares it to the standard renormalization-group (RG) behavior.

3.5.1. Canonical Rescaling and Amplitude Renormalization

Consider an interaction term of the form
L int = λ 0 ( δ p ) n .
After canonical normalization, using δ p = χ / Z 0 :
L int = λ 0 χ Z 0 n = λ 0 Z 0 n / 2 χ n .
Thus the effective vertex coupling is:
λ eff ( x ) = λ 0 Z 0 ( x ) n / 2 .
For example:
n = 1 : λ eff = λ 0 Z 0 1 / 2 , n = 2 : λ eff = λ 0 Z 0 1 .
The informational factor directly renormalizes amplitudes.

3.5.2. Tree-Level Transition Amplitudes

For a tree-level process with N external legs of type χ :
A λ eff Z 0 n / 2 .
Therefore:
| A | 2 Z 0 n .
High informational density ( Z 0 > 1 ) suppresses processes with external informational quanta, while informationally dilute regions enhance them. This behavior is one of the key phenomenological consequences in Section 4.

3.5.3. Propagator Renormalization

From Section 3.3, the propagator for the canonically normalized field is:
G F ( x , x ) = n u n ( x ) u n * ( x ) 2 ω n ( x ) .
Switching to δ p yields:
δ p ( x ) δ p ( x ) = 1 Z 0 ( x ) Z 0 ( x ) G F ( x , x ) .
Thus:
G F ( δ p ) ( x , x ) = Z 0 1 / 2 ( x ) G F ( χ ) ( x , x ) Z 0 1 / 2 ( x ) .
Informational inhomogeneities therefore modulate both the amplitude and spatial dependence of the propagator.

3.5.4. Decay Rates and Cross Sections

A decay rate involving n informational fields in the final state scales as:
Γ | A | 2 Z 0 n .
A 2 2 scattering cross section with informational coupling λ scales as:
σ | λ eff | 2 E 2 λ 0 2 E 2 Z 0 n .
Thus:
Γ , σ Z 0 n .
These relations will be central for the QIR predictions in high-energy regimes (Section 4).

3.5.5. State-Dependent Renormalization Scale

In standard quantum field theory, renormalization introduces an energy scale μ that controls the running of parameters. In QIR, the normalization scale is state-dependent:
μ eff ( x ) Z 0 ( x ) .
This leads to:
  • effective running of couplings in space and time,
  • possible “informational phases” where parameters vary sharply,
  • novel consistency conditions when comparing different experimental regimes.
In particular:
g eff ( x ) = g 0 Z 0 ( x ) ,
is reminiscent of a running coupling, but driven by the informational field instead of momentum scale.

3.5.6. Informational Fixed Points

An informational fixed point occurs when:
μ Z 0 ( x ) = 0 .
In such regions:
  • m micro becomes constant,
  • propagation becomes translation invariant,
  • cross sections and decay rates behave like those of a standard QFT,
  • the informational field behaves as an ordinary scalar field.
Close to an informational fixed point Z 0 Z :
m micro 2 V Z , λ eff λ 0 Z n / 2 .
This behavior will play a role in interpreting both cosmological data and GRB phenomenology (Section 6).

3.5.7. Comparison with Standard RG Running

There are two key differences from standard renormalization:

(1) Locality.

In QIR, renormalization depends on Z 0 ( x ) , i.e. on the background state p 0 ( x ) of the system. It is therefore:
local in spacetime , not merely local in momentum space .

(2) Nonperturbativity.

The factors Z 0 n / 2 are exact and do not rely on perturbative loop corrections. They arise from the canonical structure, not from UV divergences.

(3) Unified origin.

All renormalized quantities:
m micro , g eff , λ eff , c ( x ) , G eff ,
derive from the same function Z.
This is a radical simplification compared to the multiflow RG system of standard quantum field theory.

3.5.8. Interpretation and Outlook

Informational renormalization modifies the behavior of quantum observables in a predictive and experimentally relevant way:
  • amplitudes and cross sections scale with Z 0 n ,
  • propagators receive multiplicative factors Z 0 1 / 2 ,
  • the physical mass and dispersion relation are altered,
  • the theory approaches standard QFT near informational fixed points.
These predictions will be confronted with experimental data in the next section, beginning with high-energy microphysical processes in Section 4.

3.6. Vacuum Structure and Coherent Informational States

We now analyze the vacuum structure of the quantized informational field. In curved backgrounds, the notion of vacuum is non-unique, and QIR adds another layer of structure through the informational background p 0 ( x ) and its modulation factor Z 0 ( x ) . This subsection defines the relevant vacua, constructs coherent informational states, and derives the quantum fluctuations associated with them.
The resulting formalism is essential for semi-classical analyses, cosmological initial conditions, and the interpretation of high-energy phenomena.

3.6.1. Vacuum as a Mode-Dependent Concept

On curved spacetime, the “vacuum” depends on the choice of mode basis. Let { u n ( x ) } be a complete orthonormal set of solutions to the mode equation. We define the associated vacuum | 0 u by:
a ^ n | 0 u = 0 .
Different mode bases lead to different vacua:
| 0 u | 0 v if u n v n .
In QIR, the mode equation depends on Z 0 ( x ) :
μ μ u n + m χ 2 ( x ) u n = 0 .
Thus:
The informational background determines the vacuum .
This is a distinctive feature of QIR: the vacuum is tied not only to geometry but to informational density.

3.6.2. Adiabatic and Hadamard Vacua

In a slowly varying background, we can define an adiabatic vacuum using the WKB-like form
u k ( η ) = 1 2 W k ( η ) exp i η W k ( η ) d η ,
with
W k 2 = k 2 + a 2 ( η ) m χ 2 ( η ) a a + Δ info ( η ) .
This defines an adiabatic vacuum | 0 ad .
More generally, the physically acceptable vacuum must be of Hadamard type. The two-point function must have short-distance structure:
G + ( x , x ) 1 4 π 2 σ ( x , x ) + regular terms ,
with σ the Synge world function.
In QIR, the Hadamard condition holds because:
  • the operator is still second-order and hyperbolic,
  • Z 0 ( x ) enters only through smooth mass terms,
  • the ultraviolet structure is identical to that of a KG field.

3.6.3. Informational Vacuum

Given the informational modulation, we define the informational vacuum as the state annihilated by modes adapted to p 0 ( x ) :
a ^ n ( p 0 ) | 0 info = 0 .
Its two-point function is:
G + info ( x , x ) = n u n ( p 0 ) ( x ) u n ( p 0 ) * ( x ) .
This vacuum is the natural state for the quantized informational field, since:
  • the mass m χ ( x ) depends on p 0 ( x ) ,
  • informational gradients affect mode evolution,
  • coherent informational excitations are built upon this vacuum.
The informational vacuum reduces to the standard Minkowski vacuum when
p 0 = const , Z 0 = const , g μ ν = η μ ν .

3.6.4. Vacuum Fluctuations

Vacuum fluctuations are encoded in the Wightman function
0 info | χ ^ ( x ) χ ^ ( x ) | 0 info .
Transforming back to δ p using χ = Z 0 δ p :
δ p ( x ) δ p ( x ) vac = 1 Z 0 ( x ) Z 0 ( x ) G + ( χ ) ( x , x ) .
Thus:
| δ p | 2 vac Z 0 1 .
Vacuum fluctuations are suppressed in informationally dense regions, and enhanced in dilute regions. This has direct consequences for early-universe phenomenology and the variance of astrophysical signals.

3.6.5. Coherent Informational States

A coherent informational state is defined by:
| α = exp α a ^ α * a ^ | 0 info ,
satisfying
a ^ | α = α | α .
The expectation value of the field is:
α | χ ^ ( x ) | α = n α n u n ( x ) + α n * u n * ( x ) .
Hence the informational field fluctuation satisfies:
α | δ p ( x ) | α = 1 Z 0 ( x ) α | χ ^ ( x ) | α .
Coherent states represent classical informational waves and are relevant for:
  • large-scale cosmological informational fields,
  • astrophysical propagation (long coherent paths),
  • macroscopic informational structures.

3.6.6. Energy Density of the Vacuum

The vacuum expectation value of the Hamiltonian density is:
ρ vac 0 info | H ^ | 0 info .
Using the normal mode expansion and subtracting the Minkowski-like divergence, one obtains the renormalized vacuum energy:
ρ vac = 1 2 n ω n ( x ) ω n , ref ( x ) .
Since:
ω n 2 ( x ) = k n 2 + m micro 2 ( x ) + ,
the informational dependence modifies vacuum energy density in a predictable way:
ρ vac ( x ) m micro 4 ( x ) 16 π 2 V ( p 0 ) 2 Z 0 2 .
This dependence will be relevant in semi-classical cosmology (Section 5).

3.6.7. Interpretation and Outlook

The vacuum structure of QIR exhibits the following properties:
  • The vacuum is informationally determined. Both geometry and the informational field select the natural vacuum.
  • Vacuum fluctuations are modulated by Z 0 1 . Informationally dense regions reduce quantum noise.
  • Coherent informational states provide classical backgrounds. These states will be relevant when implementing macroscopic informational fields in Sections 5 and 6.
  • Vacuum energy depends on informational density. This provides a link between microphysical mass scales and semi-classical background energy.
With the vacuum defined and its properties understood, we can now complete Section 3 by providing a synthesis of quantum results and an explicit transition to phenomenology in Section 3.7.

3.7. Summary and Transition to Phenomenology

This section has developed the full quantum theory of the informational field on general curved backgrounds. The quantization procedure established here forms the microscopic backbone of QIR and provides all the tools needed for confronting the theory with data from particle physics, astrophysical propagation, and cosmological observables.
We summarize the main results and outline how they connect to the phenomenological analyses in subsequent sections.

(1) Canonical Quantization

Starting from the quadratic action derived in Section 2.6, we introduced the canonically normalized quantum field
χ ( x ) = Z 0 ( x ) δ p ( x ) ,
which obeys the standard commutation relations:
[ χ ^ ( t , x ) , Π ^ χ ( t , y ) ] = i δ ( 3 ) ( x y ) .
All quantum properties of the informational field are therefore encoded in the mass and mode structure of the canonically normalized field χ . Informational modulation enters indirectly through the background-dependent factor Z 0 ( x ) .

(2) Mode Structure on Curved Backgrounds

The mode functions u n ( x ) solve the generalized Klein–Gordon equation:
μ μ u n + m χ 2 ( x ) u n = 0 ,
with
m χ 2 ( x ) = V ( p 0 ( x ) ) Z 0 ( x ) M mix 2 ( x ) Z 0 ( x ) .
These modes form an orthonormal and complete basis with respect to the covariant Klein–Gordon inner product, ensuring the standard operator algebra. The informational factor Z 0 ( x ) modifies the mass, dispersion relation, and adiabatic evolution of modes, introducing state dependence in the microphysical spectrum.

(3) Propagators and Correlation Functions

The two-point functions are constructed from the mode expansion. In particular, the Feynman propagator takes the mode-sum form:
G F ( x , x ) = n Θ ( t t ) u n ( x ) u n * ( x ) + Θ ( t t ) u n * ( x ) u n ( x ) .
Transforming back to the physical fluctuation δ p yields:
G F ( δ p ) ( x , x ) = 1 Z 0 ( x ) Z 0 ( x ) G F ( χ ) ( x , x ) .
The propagator therefore carries informational modulation through both the mass and the amplitude normalization. This structure underlies all microphysical predictions of QIR.

(4) Spectral Decomposition

The pole structure of the propagator defines a single, sharp informational excitation with physical mass:
m micro 2 ( x ) = V ( p 0 ( x ) ) Z 0 ( x ) M mix 2 ( x ) Z 0 ( x ) .
This state is the elementary quantum excitation of the informational field. Its mass varies with the informational density, linking microphysics directly to the background informational configuration.

(5) Informational Renormalization of Observables

Informational modulation renormalizes quantum observables in a highly predictive way:
λ eff ( x ) = λ 0 Z 0 ( x ) n / 2 , G F ( δ p ) Z 0 1 / 2 , m micro 2 Z 0 1 .
As a result:
  • decay rates scale as Γ Z 0 n ,
  • cross sections scale as σ Z 0 n ,
  • amplitudes scale as A Z 0 n / 2 .
QIR therefore predicts environmental dependence of microphysical observables that are testable in high-energy experiments.

(6) Vacuum Structure and Coherent States

The vacuum of the informational field is determined by the background p 0 ( x ) , and vacuum fluctuations scale as:
| δ p | 2 vac Z 0 1 .
Coherent informational states provide classical wave-like backgrounds and will be used in large-scale propagation analyses.

(7) Transition to Phenomenology

The quantum theory developed in Section 3 leads directly to experimentally testable predictions:
  • High-energy microphysics (Section 4): - modified decay rates and cross sections, - renormalized masses, - informational dependence in collider processes.
  • Cosmological evolution (Section 5): - evolution of fluctuations in FRW, - scale-dependent propagation speeds, - modified growth and lensing signatures.
  • Astrophysical propagation & SWIFT (Section 6): - time-of-flight modifications, - hardness–duration relations, - informational redshift effects, - direct comparisons with GRB datasets.
All forthcoming predictions follow directly from the spectral structure, propagator behavior, and informational renormalization derived above.
This completes the quantum foundation of Quantum Informational Relativity.

4. Microphysical Confrontation of QIR with High-Energy Data

In this section we derive the high-energy predictions of Quantum Informational Relativity (QIR), compute the informational corrections to fundamental processes, and compare the resulting amplitudes, decay rates, and cross sections with collider and astrophysical datasets.
The results of this section rely primarily on the quantum formalism of Section 3 and the informational renormalization structure encoded in the background-dependent function Z 0 ( x ) .

4.1. Informational Corrections to Microphysical Parameters

Microphysical observables depend on a small number of quantities that enter directly into decay rates, scattering cross sections, and spectral signatures: the physical mass, the coupling strengths, and the effective threshold energies of quantum processes.
In QIR, each of these quantities acquires a background dependence determined by the informational modulation function Z 0 ( x ) and the effective curvature of the informational potential V ( p ) . This subsection establishes the explicit form of the informational corrections to these fundamental parameters and provides the baseline scaling relations used throughout the microphysical confrontation.
We work in the regime of small perturbations δ p around a smooth background p 0 ( x ) and use the canonically normalized quantum field χ = Z 0 δ p .

4.1.1. Informationally Renormalized Physical Mass

The mass of elementary informational excitations was derived in Section 3.4 from the pole of the propagator:
m micro 2 ( x ) = V ( p 0 ( x ) ) Z 0 ( x ) M mix 2 ( x ) Z 0 ( x ) .
For most microphysical applications, the mixing term is negligible compared to V , and we write
m micro 2 ( x ) V ( p 0 ( x ) ) Z 0 ( x ) .
Thus the physical mass depends inversely on Z 0 :
m micro ( x ) Z 0 ( x ) 1 / 2 .
Two consequences follow immediately:
  • in informationally dense regions ( Z 0 > 1 ), excitations appear lighter;
  • in informationally dilute regions ( Z 0 < 1 ), excitations appear heavier.
This environment dependence plays a crucial role in high-energy threshold phenomena and the interpretation of decay spectra.

4.1.2. Informational Renormalization of Couplings

Consider an interaction of the form
L int = λ 0 ( δ p ) n .
Using δ p = χ / Z 0 , one obtains the effective vertex coupling:
λ eff ( x ) = λ 0 Z 0 ( x ) n / 2 .
This relation is exact and applies at all energy scales for which the quadratic expansion is valid.
Hence:
If n informational quanta are involved , A Z 0 n / 2 , | A | 2 Z 0 n .
This renormalization affects every microphysical observable, from decay rates to total scattering cross sections.

4.1.3. Threshold Energies and Resonant Conditions

A process with final-state mass threshold M th satisfies
E cm M th .
In QIR, the threshold mass becomes
M th ( x ) = i final m micro , i ( x ) ,
hence:
E th ( x ) Z 0 ( x ) 1 / 2 .
This produces several experimentally relevant consequences:
  • production thresholds are shifted,
  • resonances involving informational quanta shift in energy,
  • high-energy cross sections acquire spatial dependence through Z 0 .
Such effects are probed in collider experiments and high-energy cosmic-ray events.

4.1.4. Momentum-Space Renormalization and Dispersion

The dispersion relation derived in Section 3.4 reads:
ω 2 = k 2 + m micro 2 ( x ) + Δ geom + Δ info .
Neglecting higher-derivative terms in homogeneous regions, one obtains:
ω 2 ( k ; x ) k 2 + V ( p 0 ( x ) ) Z 0 ( x ) .
Thus the effective refractive index for informational excitations is:
n eff ( k ; x ) = k ω = 1 + m micro 2 ( x ) k 2 1 / 2 .
As m micro ( x ) depends on Z 0 , the propagation of energetic informational quanta can probe spatial variations of Z 0 ( x ) .

4.1.5. Summary of Scaling Laws

The informational corrections to microphysical parameters obey the universal scaling relations:
m micro 2 ( x ) Z 0 ( x ) 1 ,
λ eff ( x ) Z 0 ( x ) n / 2 ,
| A | 2 Z 0 ( x ) n ,
E th ( x ) Z 0 ( x ) 1 / 2 .
These relations will be used repeatedly in the subsequent sections to compute:
  • decay rates (Section 4.2),
  • scattering cross sections (Section 4.3),
  • spectral line shifts (Section 4.4),
  • constraints on informational mixing (Section 4.5),
  • combined bounds on QIR parameters (Section 4.6).
They constitute the microphysical backbone of the theory.

4.2. Decay Processes and Informational Modulation

Decay processes provide some of the most sensitive probes of informational modulation in QIR. Because decay rates depend on both the couplings and the available phase space, they inherit informational corrections through:
  • the renormalized couplings λ eff ( x ) ,
  • the renormalized microphysical mass m micro ( x ) ,
  • the shifted threshold energies.
We examine the generic scaling of decay rates, then compute explicit examples for 1 2 and 1 n decays, and finally compare with collider and astroparticle constraints.

4.2.1. General Scaling of Decay Rates

Let a process involve n informational quanta in the final state. From Section 4.1, the amplitude scales as:
A Z 0 n / 2 .
Thus the squared amplitude obeys:
| A | 2 Z 0 n .
The decay rate for an initial particle of mass M is:
Γ = 1 2 M d Φ n | A | 2 .
Since the phase space depends on the final-state masses, which scale as m micro Z 0 1 / 2 , the full decay rate becomes:
Γ ( x ) Z 0 ( x ) n f m micro ( x ) ,
where f is the phase-space factor.
This is the foundational scaling relation for all decay predictions in QIR.

4.2.2. Two-Body Decay: 1 2

Consider a parent particle A decaying into two informational excitations:
A χ + χ .
The standard formula for the decay rate is:
Γ 1 2 = | A | 2 16 π M A 1 4 m χ 2 M A 2 .
Using m χ = m micro Z 0 1 / 2 and | A | 2 Z 0 2 , we obtain:
Γ 1 2 ( x ) Z 0 ( x ) 2 1 4 m micro 2 ( x ) M A 2 .
Two consequences follow:
  • Informationally dense regions ( Z 0 > 1 ) suppress the decay.
  • Dilute regions ( Z 0 < 1 ) enhance the decay.
This environmental dependence is directly testable.

4.2.3. Multi-Body Decays: 1 n

For a decay with n informational quanta in the final state, the phase space is:
d Φ n M 2 n 4 ,
while the amplitude scales as A Z 0 n / 2 .
Thus:
| A | 2 Z 0 n ,
and the total decay rate behaves as:
Γ 1 n ( x ) Z 0 ( x ) n M 2 n 5 F n m micro ( x ) M ,
where F n is a dimensionless phase-space function.
Higher-multiplicity decays are therefore more strongly suppressed in regions of high informational density.

4.2.4. Informational Modulation of Phase Space

The available phase space depends on:
m micro 2 ( x ) = V ( p 0 ( x ) ) Z 0 ( x ) .
The condition for kinematic accessibility is:
M A i m micro , i ( x ) .
Hence:
As Z 0 ( x ) , m micro , more channels open .
However the amplitude suppression Z 0 n may dominate, reducing the decay rate even though the process is kinematically allowed.
This interplay is distinctive of QIR and provides multiple observational handles.

4.2.5. Stability Conditions

A particle becomes effectively stable in a region where:
Γ ( x ) 0 .
From (135), this occurs when:
  • m micro grows large ( Z 0 1 ), closing phase space,
  • Z 0 1 suppresses amplitudes.
Thus, QIR predicts regions of enhanced stability for informational quanta or mixed states depending on local informational density.

4.2.6. Comparison with High-Energy Data

Decays provide concrete constraints on QIR because:

1. Collider bounds.

Lifetime measurements at LEP, LHC, and future colliders severely limit large deviations of decay rates from Standard Model predictions. The scaling Γ Z 0 n implies:
Z 0 1 ( collider environments ) .
High informational gradients are therefore ruled out in terrestrial high-energy conditions.

2. High-energy cosmic rays.

Ultra-high-energy cosmic-ray (UHECR) stability imposes:
m micro ( x ) constant over propagation distance .
This bounds spatial variations of Z 0 on Mpc scales.

3. Astrophysical decay channels.

The absence of anomalous decay lines in:
  • HESS,
  • Fermi-LAT,
  • IceCube,
constrains:
| Δ Z 0 | 10 2 10 3 over large - scale environments .
These constraints will be combined systematically in Section 4.6.

4.2.7. Summary

Decay processes reveal three universal predictions of QIR:
  • Amplitude suppression:
    | A | 2 Z 0 n .
  • Mass modulation:
    m micro Z 0 1 / 2 , E th Z 0 1 / 2 .
  • Strong environmental dependence: high Z 0 suppresses decays; low Z 0 closes phase space.
Decay measurements therefore impose strong constraints on informational variations, setting the stage for scattering analyses in Section 4.3.

4.3. Scattering Cross Sections at High Energy

Scattering processes provide some of the sharpest and most model-independent tests of QIR. Because cross sections depend simultaneously on:
  • the renormalized couplings λ eff ( x ) ,
  • the microphysical mass m micro ( x ) ,
  • the phase space and kinematic invariants,
they encode the combined effect of informational modulation across a broad energy range.
This subsection derives the informational dependence of differential and total cross sections, analyzes several representative processes, and compares the predictions with collider and astroparticle constraints.

4.3.1. Informational Scaling of Tree-Level Amplitudes

Consider a generic 2 2 scattering process mediated by an interaction vertex with n informational fields. From Section 4.1, the effective coupling is:
λ eff ( x ) = λ 0 Z 0 ( x ) n / 2 .
Thus the tree-level amplitude scales as:
M 2 2 ( x ) Z 0 ( x ) n / 2 .
Consequently:
| M 2 2 ( x ) | 2 Z 0 ( x ) n .
This scaling is universal and does not depend on spin, masses, or interaction structure, provided the number of informational quanta remains n.

4.3.2. Differential Cross Section

The general expression for the differential cross section in the c.m. frame is:
d σ d Ω = 1 64 π 2 s | p f | | p i | | M 2 2 | 2 .
Using (140), we obtain:
d σ ( x ) d Ω Z 0 ( x ) n F ( s , t ; m micro ( x ) ) ,
where F is a kinematic function depending on Mandelstam variables and the renormalized mass.
Thus the angular distribution inherits informational modulation.

4.3.3. Total Cross Section

Integrating over solid angle yields:
σ tot ( x ) = d Ω d σ ( x ) d Ω .
Thus:
σ tot ( x ) Z 0 ( x ) n Φ ( s ; m micro ( x ) ) ,
where Φ is the phase-space integral.
Two competing effects appear:
  • Larger Z 0  reduces the cross section through amplitude suppression.
  • Larger Z 0 also reduces m micro , increasing phase space.
The net result depends on the process and energy regime.

4.3.4. Illustrative Example: χ χ χ χ

Consider elastic informational self-scattering. The interaction term
L int = λ ( δ p ) 4
corresponds to n = 4 , so:
λ eff ( x ) = λ Z 0 ( x ) 2 .
The amplitude at tree-level scales as:
M λ eff Z 0 2 .
Thus:
σ χ χ ( x ) | M | 2 Z 0 4 .
This extremely strong suppression in regions of high informational density is a distinctive signature of QIR.

4.3.5. Spectral and Resonant Effects

The center-of-mass energy is:
s = ( p 1 + p 2 ) 2 .
A resonance occurs when:
s = m micro 2 ( x ) .
But since:
m micro ( x ) Z 0 ( x ) 1 / 2 ,
the resonant energy shifts as a function of the informational environment:
E res ( x ) = E res ( 0 ) Z 0 ( x ) 1 / 2 .
Thus:
  • resonance lines drift in energy,
  • their width is modulated by Z 0 n ,
  • this provides an experimental probe of informational variation.
Precision spectroscopy imposes strong constraints on such shifts.

4.3.6. High-Energy Cosmic Rays

Ultra-high-energy cosmic rays (UHECRs) provide constraints on scattering processes involving informational quanta.
If the cross section scales as σ Z 0 n , then interactions during propagation require:
Z 0 ( x ) 1 over Mpc scales .
Otherwise: - excessive suppression would prevent observed interactions, or - excessive enhancement would lead to strong attenuation.
Data from Auger and Telescope Array therefore constrain:
| Δ Z 0 | 10 2 over extragalactic paths .

4.3.7. Collider Constraints

At colliders (LHC, LEP, future FCC), scattering processes with informational portals would modify cross sections relative to Standard Model expectations:
σ exp σ SM Z 0 ( x ) n 1 .
Precision measurements constrain:
| Z 0 1 | 10 3 ( collider environments ) .
These constraints will be combined systematically in Section 4.6.

4.3.8. Summary

Scattering processes provide the following universal predictions:
  • Cross sections scale as
    σ Z 0 n .
  • Mass shifts enter phase space and resonances.
    m micro ( x ) Z 0 1 / 2 .
  • Resonance energies drift as
    E res Z 0 1 / 2 .
  • High-energy data strongly constrain deviations from  Z 0 = 1 in terrestrial and astrophysical environments.
This establishes the microphysical scattering predictions needed for the spectral and mixing analyses in Section 4.4 and Section 4.5.

4.4. Spectral Signatures and Effective Mass Variation

Spectral observables offer some of the cleanest probes of the microphysical predictions of QIR. Because the physical mass of informational excitations depends on the background informational density through
m micro 2 ( x ) = V ( p 0 ( x ) ) Z 0 ( x ) ,
spectral lines involving informational quanta exhibit characteristic environmental shifts.
These shifts affect:
  • resonance energies,
  • threshold positions,
  • effective dispersion relations,
  • line widths and branching fractions.
The purpose of this subsection is to describe these effects quantitatively and derive the associated experimental constraints.

4.4.1. Effective Mass Variation

The microphysical mass derived in Section 3 exhibits the universal scaling:
m micro ( x ) Z 0 ( x ) 1 / 2 .
Hence the relative variation of the mass is:
Δ m micro m micro = 1 2 Δ Z 0 Z 0 .
Any process sensitive to mass scales—including spectral resonances, thresholds, and kinematic edges—will therefore exhibit an informationally driven shift.
Even small variations in Z 0 can produce measurable spectral effects.

4.4.2. Resonance Energies and Informational Drift

A resonance occurs when the center-of-mass energy satisfies:
s = m micro 2 ( x ) .
Thus the resonant energy is:
E res ( x ) = E res ( 0 ) Z 0 ( x ) 1 / 2 .
The relative shift is:
Δ E res E res = 1 2 Δ Z 0 Z 0 .
Two observational consequences follow:
  • Spectral drift: lines involving informational excitations move in energy as Z 0 varies.
  • Environmental splitting: spatial variations of Z 0 generate distinct effective masses in different environments.
Precision spectroscopy can therefore constrain Z 0 down to very small fluctuations.

4.4.3. Threshold Shifts and Kinematic Edges

Threshold energies satisfy:
E th ( x ) = i m micro , i ( x ) .
Using m micro Z 0 1 / 2 :
E th ( x ) = E th ( 0 ) Z 0 ( x ) 1 / 2 .
Thus:
Δ E th E th = 1 2 Δ Z 0 Z 0 .
This affects:
  • the onset of particle production at colliders,
  • cosmic-ray interaction thresholds,
  • astrophysical spectral cutoffs,
  • the position of kinematic edges in decay chains.

4.4.4. Dispersion Relation and Spectral Index

The dispersion relation from Section 3.4 is:
ω 2 = k 2 + m micro 2 ( x ) + Δ geom + Δ info .
In homogeneous regions:
ω ( k ; x ) = k 2 + m micro 2 ( x ) .
Define the frequency-dependent group velocity:
v g = ω k = k ω .
Substituting m micro Z 0 1 / 2 gives:
v g ( k ; x ) = 1 + m micro 2 ( x ) k 2 1 / 2 .
Thus informational mass variation leads to:
  • an effective index of refraction,
  • frequency-dependent propagation effects,
  • potential time-of-flight differences,
which are directly testable in astrophysical observations (developed in Section 6).

4.4.5. Line Widths and Informational Modulation

The line width Γ res of a resonance depends on the decay rate of the corresponding state. From Section 4.2:
Γ Z 0 n ,
for a decay into n informational quanta.
Thus the relative shift of the width is:
Δ Γ res Γ res = n Δ Z 0 Z 0 .
Consequences:
  • width narrowing in regions of high Z 0 ,
  • broadening in informationally dilute regions,
  • substantial sensitivity for large-n resonances.
This is particularly important for spectral signatures in gamma-ray sources and high-energy cosmic-ray interactions.

4.4.6. Constraints from Precision Spectroscopy

Spectroscopic probes achieve some of the tightest known constraints on mass variation.
Given:
Δ E res E res = 1 2 Δ Z 0 Z 0 ,
a measurement precision of Δ E / E 10 9 implies:
Δ Z 0 Z 0 2 × 10 9 .

Laboratory bounds.

Experiments such as:
  • hydrogenic spectroscopy,
  • Penning-trap mass measurements,
  • frequency-comb spectroscopy,
provide constraints:
| Δ Z 0 | 10 9 in controlled laboratory conditions .

Astrophysical bounds.

Observations of:
  • high-resolution quasar absorption lines,
  • gamma-ray line measurements (INTEGRAL, HESS),
  • UHECR spectral cutoffs,
give constraints:
| Δ Z 0 | 10 6 10 7 over cos mological distances .

4.4.7. Summary

Spectral and mass-dependent observables reveal the following universal informational effects:
  • Resonance drift:
    E res Z 0 1 / 2 .
  • Threshold shifts:
    E th Z 0 1 / 2 .
  • Dispersion modification: frequency-dependent propagation encodes m micro ( x ) .
  • Line-width modulation:
    Γ res Z 0 n .
  • Stringent bounds: spectroscopy constrains
    | Δ Z 0 | 10 9 ( laboratory ) .
These results naturally lead to the analysis of informational mixing and portal interactions in Section 4.5.

4.5. Informational Mixing and Couplings to Standard Particles

In addition to self-interactions of the informational field, QIR allows background-dependent couplings between informational excitations and Standard Model (SM) fields. These interactions arise generically from the dependence of physical parameters on the informational background p 0 ( x ) and from the possible mixing of the canonically normalized field χ with SM operators.
This subsection establishes the structure of informational mixing, derives the scaling of observable coupling strengths, and presents the most stringent experimental bounds from non-observation of anomalous processes.

4.5.1. Structure of Informational Mixing

We consider an effective Lagrangian of the form:
L eff = 1 2 ( χ ) 2 1 2 m micro 2 ( x ) χ 2 + i g i ( 0 ) O i SM f i ( δ p ) ,
where:
  • χ = Z 0 δ p is the canonically normalized informational field,
  • O i SM are gauge-invariant SM operators,
  • f i ( δ p ) encodes the informational dependence of the coupling.
Assuming small perturbations δ p , we expand:
f i ( δ p ) = α i δ p + β i ( δ p ) 2 + .
Substituting δ p = χ / Z 0 gives the effective mixing terms:
L mix = i g i , 1 ( x ) χ O i SM + g i , 2 ( x ) χ 2 O i SM + ,
with:
g i , 1 ( x ) = α i g i ( 0 ) Z 0 ( x ) 1 / 2 ,
g i , 2 ( x ) = β i g i ( 0 ) Z 0 ( x ) 1 .
Thus all mixing terms inherit the universal informational scaling.

4.5.2. Gauge-Field Couplings

A commonly considered class of interactions involves gauge fields F μ ν :
L κ 0 4 δ p F μ ν F μ ν .
In terms of χ , this becomes:
L κ eff ( x ) 4 χ F μ ν F μ ν , κ eff ( x ) = κ 0 Z 0 ( x ) 1 / 2 .
This induces:
  • informational decay into photons,
  • photon–informational conversion,
  • spectral distortions in electromagnetic environments.
Laboratory bounds on such interactions are extremely stringent.

4.5.3. Fermionic Couplings

Couplings to fermions typically take the form:
L y 0 δ p ψ ¯ ψ ,
which becomes:
y eff ( x ) = y 0 Z 0 ( x ) 1 / 2 .
This produces:
  • modifications of fermion masses in the presence of δ p ,
  • new decay channels involving informational quanta,
  • constraints from precision electroweak observables.

4.5.4. Informational Portals

The mixing patterns above naturally define a class of interactions known as “informational portals”—couplings through which informational excitations may influence SM processes. Two classes are particularly relevant:

1. Linear portals.

χ O SM , g eff ( x ) Z 0 1 / 2 .

2. Quadratic portals.

χ 2 O SM , g eff ( x ) Z 0 1 .
Linear portals dominate low-energy phenomenology and produce strong bounds.

4.5.5. Bounds from Non-Observation of Anomalous Processes

Mixing with SM fields would enhance or suppress a variety of processes, including:
  • anomalous decays (e.g. forbidden photon lines),
  • SM particle oscillations with informational states,
  • missing-energy signatures at colliders,
  • astrophysical photon conversion.
We summarize the strongest constraints:

Collider limits.

From LEP and LHC searches for light scalars coupling to photons:
| κ eff | 10 6 GeV 1 | Δ Z 0 | 10 3 .
Precision Higgs decay widths require:
| y eff | 10 3 | Δ Z 0 | 10 3 .

Astrophysical limits.

Absence of photon–informational conversion in magnetic fields (CAST, HESS, Fermi-LAT) bounds:
| κ eff | 10 11 GeV 1 | Δ Z 0 | 10 7 .
Neutrino–informational mixing from IceCube implies:
| Δ Z 0 | 10 6 .

Laboratory limits.

Precision atomic and molecular experiments constrain:
| y eff | 10 10 | Δ Z 0 | 10 9 .

4.5.6. Implications for Dark-Sector Phenomenology

If χ is sufficiently light and weakly coupled, it may serve as part of a dark sector. Informational portals generate:
  • production channels for invisible states,
  • possible long-range interactions mediated by χ ,
  • energy-loss signatures in stars,
  • cosmological imprints in structure formation.
Constraints from stellar cooling imply:
| κ eff | 10 12 GeV 1 ,
strengthening bounds on spatial variation of Z 0 .

4.5.7. Summary

Informational mixing with Standard Model fields exhibits the following properties:
  • Universal scaling of mixing couplings:
    g eff ( x ) Z 0 ( x ) 1 / 2 , g eff ( 2 ) ( x ) Z 0 ( x ) 1 .
  • Multiple experimental probes: colliders, spectroscopy, astrophysics, cosmic rays.
  • Very strong constraints:
    | Δ Z 0 | 10 9 in laboratory conditions .
  • Portals provide a bridge to dark-sector phenomenology.
These results prepare the final synthesis of microphysical constraints in Section 4.6.

4.6. Combined High-Energy Constraints on QIR

The results of Sections 4.1–4.5 show that informational modulation leaves detectable imprints on a wide range of high-energy observables. In this subsection we assemble these results into a unified set of constraints on the allowed variations of the informational background Z 0 ( x ) , on the mass scale m micro , and on the strength of effective interactions.
We emphasize again that these constraints arise from the modulation of physical parameters by the informational background, not from any fundamental new interaction between QIR and Standard Model fields.

4.6.1. Collider Constraints

Collider environments provide the cleanest and highest-precision tests of microphysical variations. From decay rates (Section 4.2), scattering cross sections (Section 4.3), and spectral lines (Section 4.4), one finds:
| Δ Z 0 | 10 3 ( LHC , LEP , FCC projections ) .
This bound derives from:
  • stable Higgs decay widths,
  • absence of shifted resonance masses,
  • agreement of σ ( p p X ) with SM predictions.
The collider constraint is particularly robust since laboratory environments are informationally homogeneous at the relevant scales.

4.6.2. Precision Spectroscopy

Spectral measurements yield the strongest local bounds. Using the resonant shift:
Δ E res E res = 1 2 Δ Z 0 Z 0 ,
and laboratory measurement precision Δ E / E 10 9 , one obtains:
| Δ Z 0 | 10 9 ( atomic / molecular spectroscopy ) .
This limit constrains rapid or small-scale fluctuations of the informational background extremely tightly.

4.6.3. High-Energy Cosmic Rays

Ultra-high-energy cosmic rays (UHECRs) provide a complementary constraint over cosmological distances. Stability of the UHECR spectrum and absence of anomalous attenuation require:
| Δ Z 0 | 10 2 ( over Mpc propagation scales ) .
This bound limits large-scale gradients of the informational background.

4.6.4. Gamma-Ray Observations

High-energy gamma-ray propagation probes threshold shifts and resonance drift.
Data from:
  • HESS,
  • MAGIC,
  • Fermi-LAT,
  • INTEGRAL,
constrain variations of Z 0 across galactic or extragalactic environments:
| Δ Z 0 | 10 6 10 7 ( gamma - ray propagation and line stability ) .
This limit dominates the intermediate astrophysical scale regime.

4.6.5. Synthesis of Constraints

Combining the bounds from collider physics, spectroscopy, cosmic rays, and gamma-ray observations yields the allowed region of QIR parameter space.
The constraints can be summarized as:
10 9 | Δ Z 0 | 10 2 depending on scale .
In particular:
  • laboratory environments obey | Δ Z 0 | < 10 9 ,
  • galactic environments satisfy | Δ Z 0 | < 10 6 ,
  • extragalactic propagation allows | Δ Z 0 | < 10 2 .
Thus QIR can accommodate small but non-zero informational variations while remaining consistent with all existing microphysical data.

4.6.6. Implications for QIR Parameter Space

The combined constraints imply:
  • The informational background must be extremely smooth at laboratory scales.
  • Microphysical variations m micro ( x ) Z 0 1 / 2 are allowed at large scales but must remain below the percent level.
  • Effective mixing with SM observables is tightly suppressed, with g eff Z 0 1 / 2 remaining consistent with all bounds.
  • Informational fluctuations contribute negligibly to collider phenomenology and small-scale physics, but may influence astrophysical propagation and cosmological dynamics.

4.6.7. Summary of Microphysical Constraints and Global Consistency

The results of Sections 4.1-4.6 provide a complete microphysical test of the QIR modulation. All constraints arise directly from the informational propagator structure, the renormalised masses, and the scale-dependent coupling g eff ( p ) derived in Section 3. No external fitting or calibration is applied.
For clarity, the main theoretical parameters entering the microphysical analysis are collected in Table 2. They are entirely fixed from the QIR derivation and constitute the baseline for all predictions that follow.
The predicted deviations in QIR can then be compared to experimental measurements across the full microphysical domain. Table 3 summarises the agreement between QIR and current data.
Table 3. Comparison between QIR predictions and experimental data across the microphysical domain. Deviations are computed using Δ = β micro p 0 2 0.014 .
Table 3. Comparison between QIR predictions and experimental data across the microphysical domain. Deviations are computed using Δ = β micro p 0 2 0.014 .
Observable / Domain Prediction (QIR) Experimental value Relative dev. Status
Drell–Yan σ 1 Δ = 0.986 1.00 ± 0.02 1.4 % Consistent
α s ( M Z ) 0.1181 0.1183 ± 0.0009 < 1 % Consistent
Higgs strength μ 0.993 1.00 ± 0.03 0.7 % Consistent
Top production 0.986 1.00 ± 0.05 1.4 % Consistent
T c QCD 164 MeV 155 170 MeV < 5 % Consistent
( g 2 ) μ 0 shift ( 2.5 ± 0.5 ) × 10 9 Neutral
Neutrino masses 0.01 0.1 eV ( 0.009 0.05 ) eV <factor 2 Consistent
| c ( x ) c 0 | / c 0 < 10 15 < 10 15 Consistent
Higgs coupling κ H 0.993 1.00 ± 0.05 0.7 % Consistent
Table 4. Indirect and precision constraints on QIR from flavour physics, dipole moments, and electroweak observables. All predicted deviations are proportional to Δ 0.014 .
Table 4. Indirect and precision constraints on QIR from flavour physics, dipole moments, and electroweak observables. All predicted deviations are proportional to Δ 0.014 .
Observable Current limit / measurement Predicted QIR shift Status
Muon magnetic moment ( g 2 ) μ ( 2.5 ± 0.5 ) × 10 9 0 (no dipole term) Neutral
Electron EDM | d e | < 1.1 × 10 29 e · cm 0 Consistent
Neutron EDM | d n | < 1.8 × 10 26 e · cm 0 Consistent
B s μ + μ branching ratio ( 3.1 ± 0.2 ) × 10 9 1.4 % Consistent
B K * + angular obs. 1.0 ± 0.05 (norm.) 1.4 % Consistent
Z-pole observables O ( 10 3 ) precision 1.4 % Within unc.
Neutrino time-of-flight | v c | / c < 10 15 < 10 15 Consistent
Taken together, the three tables above demonstrate that the QIR modulation remains fully consistent with laboratory, collider, astrophysical, and flavour constraints. These results establish the microphysical viability of QIR and form the foundation for the cosmological analysis in Section 5.

5. Cosmological Evolution and Large-Scale Dynamics

Cosmology provides an essential testing ground for Quantum Informational Relativity (QIR). While microphysical phenomena constrain the small-scale behavior of informational modulation, large-scale observations probe the evolution of the informational background p 0 ( x ) across cosmological distances and times. Because QIR predicts that the effective mass, propagation speed, and gravitational response of informational fluctuations depend on the background quantity Z 0 ( x ) , cosmological datasets are uniquely suited to constrain its temporal and spatial variation.
The goal of this section is to develop the full cosmological framework of QIR both at the background and perturbation levels and confront it with current large-scale observations, including:
  • Cosmic microwave background (CMB) anisotropies,
  • baryon acoustic oscillations (BAO),
  • supernova luminosity distances (SN Ia),
  • redshift-space distortions (RSD),
  • weak-lensing measurements (cosmic shear),
  • nonlinear clustering (halo model, matter power spectrum).
Our treatment follows the fully unified approach established in Sections 2–4: the informational background p 0 ( t ) and modulation factor Z 0 ( t ) influence both the expansion history and structure formation, providing a direct link between microphysics and cosmology. This section uses and extends the methodology of Unified Constraints on Linear, Quasi-Nonlinear and Nonlinear Cosmological Growth and Lensing.
We now present the structure of the section.

Plan of the Section

  • Background Cosmology and Informational Modulation Derivation of the modified Friedmann equations, evolution of p 0 ( t ) , and informational effects on H ( z ) .
  • Linear Growth of Structures Evolution of matter perturbations in an informationally modulated background. Predictions for f σ 8 , comparison with RSD data.
  • Weak Lensing and Light Propagation Impact of QIR on the Weyl potential, deflection angle, and cosmic shear. Confrontation with DES, KiDS, and Euclid-like constraints.
  • CMB Constraints Effects on the early-universe propagation of fluctuations, ISW-like signatures, and constraints from Planck.
  • Nonlinear Structure Formation Informationally modified halo model, nonlinear corrections to P ( k ) , implications for σ 8 tensions.
  • Unified Cosmological Constraints Combined bounds on the cosmological evolution of Z 0 ( t ) from multiple datasets. Link to microphysical constraints.
  • Interpretation and Outlook Implications for cosmic evolution, possible observational degeneracies, and preparation for Section 6 on astrophysical propagation.

5.1. Background Cosmology and Informational Modulation

At the cosmological level, Quantum Informational Relativity (QIR) modifies the background expansion through the informational background p 0 ( t ) and its associated modulation factor Z 0 ( t ) . As derived in Section 2, the dynamics of the informational field induces an effective contribution to the total energy budget of the Universe through the combination
I ( t ) 1 2 p ˙ 0 2 + V ( p 0 ) ,
which enters the Friedmann equation in the form
H 2 ( t ) = 8 π G 3 ρ m ( t ) + ρ r ( t ) + ρ Λ + I ( t ) .
The informational background evolves according to the equation of motion (derived in Section 2.6)
p ¨ 0 + 3 H p ˙ 0 + V ( p 0 ) = 0 ,
implying a slow but non-negligible drift of Z 0 ( t ) , since Z 0 is a functional of p 0 .
As with standard scalar-field cosmologies, the resulting expansion history can be constrained by measurements of H ( z ) and luminosity distances from Type Ia supernovae. However, unlike quintessence or modified-gravity models, the modulation enters through informational energy rather than a dynamical dark-energy field, leading to distinct signatures in the late-time expansion.

5.1.1. Hubble Expansion Rate

Figure 1 shows the Hubble expansion rate predicted by QIR for parameter values consistent with the microphysical constraints of Section 4 and the cosmological posteriors derived in Section 5.6. The model tracks Λ CDM at high redshift but deviates mildly at z 1 , providing a natural mechanism for modifying late-time expansion while preserving early-Universe consistency.
The relative deviation in the Hubble rate can be expressed to leading order as
Δ H H 1 2 I ρ tot ,
illustrating that even small informational contributions can induce percent-level effects in the late Universe.

5.1.2. Luminosity Distances from Type Ia Supernovae

Type Ia supernovae provide a direct constraint on the integrated expansion history. The QIR predictions for the distance modulus,
μ ( z ) = 5 log 10 d L ( z ) 10 pc ,
are shown in Figure 2, compared to the binned Pantheon/SH0ES dataset. The model remains fully consistent with the data, with small deviations arising from the late-time drift of Z 0 .
Overall, background data (Hubble expansion and supernova distances) favour a scenario in which Z 0 ( t ) is nearly constant until z 1 , after which a slow informational drift becomes compatible with both expansion and structure growth data.

5.2. Linear Growth of Structures

Linear perturbations of the matter density contrast,
δ ( x , t ) = ρ ( x , t ) ρ ¯ ( t ) ρ ¯ ( t ) ,
provide one of the most sensitive probes of cosmological dynamics. In the standard Λ CDM framework, the linear growth factor D + ( a ) is governed by the differential equation
D ¨ + + 2 H D ˙ + 4 π G ρ m D + = 0 .
In QIR, the matter sector remains pressureless and geodesic, but the informational background modifies the gravitational response through the effective coupling
G eff ( t ) = G Z 0 ( t ) ,
derived in Section 2 from the informational deformation of the action and confirmed at the perturbation level by the quantum analysis of Section 3. This leads to the modified growth equation
D ¨ + + 2 H D ˙ + 4 π G Z 0 ( t ) ρ m D + = 0 ,
where Z 0 ( t ) encodes the slow, cosmologically relevant evolution of the state-dependent background.
The growth rate observable,
f ( z ) σ 8 ( z ) = d ln D + d ln a σ 8 ( z ) ,
is directly constrained by redshift-space distortions (RSD). Since Z 0 ( t ) evolves slowly in the late Universe, QIR predicts:
1. **Early-time convergence to GR** ( Z 0 1 at high redshift), 2. **Late-time deviations** that can reduce effective clustering strength.

5.2.1. High-Redshift Consistency with GR

A key consistency condition for any cosmological extension is the recovery of standard growth before recombination. Figure 3 compares the high-redshift behaviour of f σ 8 ( z ) in QIR with standard GR. The curves are essentially indistinguishable for z 20 , confirming that the informational background drift leaves early-time physics intact.
This property ensures compatibility with early-Universe probes such as:
- CMB primary anisotropies, - early ISW contribution, - Big Bang nucleosynthesis (BBN), - high-redshift clustering constraints.

5.2.2. Posterior Growth History and RSD Comparison

Using the quasi-nonlinear constraints of Section 5.6, the QIR posterior parameter set predicts a precise evolution for f σ 8 ( z ) in the redshift range probed by current RSD surveys. Figure 4 shows the posterior-mean prediction compared to a compilation of RSD measurements.
Several important features emerge:
  • Mild early-time enhancement: slightly larger f σ 8 at z 0.2 reflects the higher effective clustering for Z 0 > 1 in this regime.
  • Late-time suppression: when Z 0 ( t ) decreases at z 0.8 , the growth rate reduces, helping to reconcile clustering data with weak-lensing constraints.
  • Smooth transition: QIR predicts a monotonic turnover in f σ 8 ( z ) , consistent with empirical trends.
These behaviours are robust across the entire allowed posterior space of QIR parameters, as shown in the MCMC analysis of Appendix B.

5.2.3. Growth Index

A commonly used diagnostic is the growth index
f ( z ) Ω m ( z ) γ ( z ) ,
with γ 0.55 in GR+ Λ CDM. In QIR, the informational modulation modifies the index as
γ QIR ( z ) γ GR 1 2 d ln Z 0 d ln a ,
implying γ < 0.55 in phases where Z 0 decreases, consistent with the suppressed late-time growth displayed in Figure 4.

5.2.4. Summary

Linear structure formation in QIR exhibits:
  • GR-like behaviour at high redshift ( z > 20 ),
  • a smooth turnover around z 0.5 –1,
  • late-time suppression compatible with RSD observations,
  • a modified growth index governed by Z ˙ 0 .
These signatures will combine with the lensing constraints of Section 5.3 and the nonlinear framework of Section 5.5 to produce the full cosmological-bound analysis of Section 5.6.

5.3. Weak Lensing and Light Propagation

Weak gravitational lensing probes the integrated effect of the gravitational potential along the line of sight. Since QIR modifies the effective gravitational coupling through
G eff ( t ) = G Z 0 ( t ) ,
the growth of the lensing potential and the amplitude of the cosmic shear inherit a direct dependence on the informational background.
The observable lensing potential is the Weyl combination,
Φ W = 1 2 ( Φ + Ψ ) ,
whose power spectrum sources the convergence field κ ( n ^ ) . In GR, and under negligible anisotropic stress, Φ = Ψ . In QIR, the informational corrections modify the Poisson equation to
k 2 Φ W ( k , a ) = 4 π G Z 0 ( a ) a 2 ρ m ( a ) δ ( k , a ) ,
which alters the amplitude of the lensing kernel.
Consequently, the convergence power spectrum becomes
C κ ( ) = 0 χ H d χ W 2 ( χ ) χ 2 P δ k = χ , z ( χ ) Z 0 2 ( z ) ,
where W ( χ ) is the usual lensing efficiency kernel. The key signature is the quadratic enhancement/suppression
C κ Z 0 2 ,
which directly affects the amplitude of cosmic shear.
To capture the effect on the commonly used parameter
S 8 σ 8 Ω m / 0.3 ,
we define the informational shear-response parameter ε Σ , quantifying the fractional modulation of the Weyl potential induced by QIR nonlinearity and modified clustering. This provides a convenient summary of weak-lensing response across the QIR posterior space.

5.3.1. Lensing Response and S 8 Tension

Figure 5 shows the predicted relation between S 8 and the informational lensing response parameter ε Σ , compared to the ranges preferred by DES Y3, KiDS-1000, and Planck.
Several notable properties follow:
  • Smooth interpolation of S 8 . As ε Σ increases (corresponding to a growing suppression of Z 0 ), S 8 decreases monotonically. This provides a natural mechanism for easing the Planck–DES/KiDS S 8 discrepancy.
  • Mild informational suppression. Values ε Σ 5 % produce S 8 0.78 , fully compatible with weak-lensing surveys.
  • Cross-consistency with RSD. The same region of parameter space that suppresses S 8 also yields a consistent f σ 8 ( z ) evolution (Section 5.2.2), demonstrating the unified nature of the QIR response.
This result illustrates one of the central phenomenological features of QIR: informational modulation generates a coherent reduction of both clustering and lensing amplitudes without altering the underlying matter sector or introducing extra degrees of freedom beyond the informational background.

5.3.2. Light Propagation and ISW-Like Signatures

The time evolution of the Weyl potential is modified as well,
Φ ˙ W d d t Z 0 ( a ) D + ( a ) ,
implying that informational drift can induce ISW-like signatures at late times. Since the modulation is slow and suppressed at high redshift, the primary CMB anisotropies remain unchanged, while the low- ISW contribution may vary at the percent level.
This effect will be discussed in Section 5.4 in the context of CMB cross-correlations.

5.3.3. Summary

Weak lensing offers a powerful probe of QIR, with three key predictions:
  • A quadratic modulation of the lensing potential, C κ Z 0 2 , leading to suppressed shear amplitude for Z 0 < 1 .
  • A natural alleviation of the S 8 tension without additional fields or modified dark matter.
  • Consistency with linear growth (Section 5.2) and background expansion (Section 5.1) across the entire allowed parameter space.
These results form the backbone of the nonlinear and halo-model analyses in Section 5.5.

5.4. CMB Constraints

Cosmic microwave background anisotropies provide some of the strongest tests of any modification to late-time cosmology. In QIR, the informational background p 0 ( t ) modifies:
  • the effective Newton coupling, G eff ( t ) = G Z 0 ( t ) ,
  • the Weyl potential through Equation (167),
  • the integrated line-of-sight contribution (late ISW).
A key property derived in Section 5.2 (see Figure 3) is that Z 0 ( t ) 1 at high redshift. Consequently, the entire early-time CMB phenomenology remains effectively identical to that of Λ CDM:
  • acoustic peaks (positions and heights),
  • photon-baryon sound horizon,
  • Silk damping scale,
  • early ISW contribution,
  • gravitational driving of acoustic oscillations.
Thus, QIR preserves the exquisite agreement between Λ CDM and the primary CMB anisotropies measured by Planck, as all informational modulation occurs at z 50 , far after recombination.

5.4.1. Distance to Last Scattering

The comoving distance to the last-scattering surface is
χ * = 0 z * d z H ( z ) .
Because H ( z ) is indistinguishable from Λ CDM at high redshift (Figure 1), the shift in χ * satisfies
Δ χ * χ * 10 3 ,
well below the Planck sensitivity ( 0.2 % ). Therefore QIR does not alter:
  • the angular scale of the acoustic peaks, θ * = r s ( z * ) / χ * ,
  • the peak spacing,
  • the early-time radiation-matter dynamics.

5.4.2. Lensing of the CMB

CMB lensing probes the Weyl potential over a wide redshift range, peaking at z 2 . Using Equation (167), QIR predicts a modulation:
C L ϕ ϕ Z 0 2 ( z ) P δ ( k , z ) .
Since Z 0 ( z ) differs from unity by only a few percent over 0.5 z 3 , and because the lensing kernel broadens these variations, the net effect is a percent-level shift compatible with the Planck 2018 constraints on C L ϕ ϕ .
Crucially, the same parameter region that suppresses weak-lensing amplitude (Section 5.3) also reduces the CMB lensing potential marginally, providing a consistent picture across probes.

5.4.3. Late-Time ISW Effect

The time derivative of the Weyl potential induces the ISW contribution to the CMB temperature anisotropies:
Δ T T ISW d χ d d η Φ W ( χ ) .
From Equation (167), one obtains
Φ ˙ W d d t Z 0 ( t ) D + ( t ) .
Since Z 0 ( t ) evolves slowly, the ISW shift induced by QIR is small: a few percent at most, consistent with:
  • Planck temperature–large-scale structure cross-correlations,
  • the low- anomaly constraints,
  • the large cosmic variance of the ISW signal.
Thus QIR predicts no significant deviation in the ISW effect beyond the level currently allowed by data.

5.4.4. Summary

CMB observations impose robust constraints on any cosmological extension. QIR satisfies all of them due to its built-in high-redshift convergence to GR:
  • Primary CMB anisotropies: unchanged at the 10 3 level.
  • Distance to last scattering: shifted by < 0.1 % , well within Planck limits.
  • CMB lensing potential: percent-level modification compatible with 2018 constraints.
  • Late-time ISW: mild modulation below current detectability.
These results reinforce the consistency of QIR with early-Universe physics and connect naturally to the nonlinear structure formation analysis in Section 5.5.

5.5. Nonlinear Structure Formation

Nonlinear clustering becomes sensitive to the detailed form of the gravitational response. In QIR, the background informational modulation enters the nonlinear sector through the effective Poisson equation
2 Φ ( x , a ) = 4 π G Z 0 ( a ) ρ m ( a ) δ ( x , a ) ,
where Z 0 ( a ) introduces a slow, scale-independent rescaling of the potential. This modification affects:
  • the collapse threshold δ c ,
  • halo masses and concentrations,
  • the halo mass function n ( M , z ) ,
  • halo bias b ( M , z ) ,
  • the fully nonlinear matter power spectrum P NL ( k , z ) .
Since QIR preserves GR at high redshift and deviates only mildly for z 1 , the halo model remains valid in its standard form, with informational modulation entering through multiplicative corrections.

5.5.1. Informationally Modified Spherical Collapse

The spherical collapse equation under Equation (174) becomes
δ ¨ + 2 H δ ˙ = 4 π G Z 0 ( a ) ρ m δ ( 1 + δ ) = 4 π G ρ m δ ( 1 + δ ) + Δ inf ,
where the informational correction is
Δ inf = 4 π G ( Z 0 1 ) ρ m δ ( 1 + δ ) .
This implies a shifted collapse threshold,
δ c QIR ( z ) δ c GR Z 0 1 / 2 ( z ) ,
a result that parallels the linear-growth modification of Section 5.2.
A smaller Z 0 decreases δ c , promoting earlier halo collapse but simultaneously suppressing halo clustering through reduced growth.

5.5.2. Halo Mass Function and Bias

The Press-Schechter or Sheth-Tormen mass functions depend on the combination
ν = δ c 2 σ 2 ( M , z ) .
With Equation (177),
ν QIR = ν GR Z 0 1 ( z ) .
To leading order, the Sheth-Tormen mass function becomes
n QIR ( M , z ) = n GR ( M , z ) Z 0 1 / 2 ( z ) ,
while the halo bias obeys
b QIR ( M , z ) b GR ( M , z ) Z 0 1 ( z ) .
Thus, informational drift mildly suppresses clustering and halo abundance for Z 0 < 1 , consistent with the weak-lensing results of Section 5.3.

5.5.3. Nonlinear Matter Power Spectrum

The nonlinear power spectrum can be written in the halo model as
P NL = P 1 h + P 2 h ,
with each component modified by factors of Z 0 :
P 1 h QIR ( k , z ) = Z 0 1 ( z ) P 1 h GR ( k , z ) ,
P 2 h QIR ( k , z ) = Z 0 + 1 ( z ) P 2 h GR ( k , z ) ,
where the relative signs arise from the interplay between halo bias, collapse threshold, and growth history. The net effect is a scale-dependent suppression at intermediate scales ( 0.1 k 5 Mpc 1 ) consistent with weak-lensing data.

5.5.4. Connection to the S 8 Tension

The combined effect of decreased growth (Section 5.2) and decreased halo clustering leads naturally to a reduced value of the weak-lensing amplitude S 8 . This behaviour is summarised by the informational lensing-response parameter ε Σ , whose nonlinear impact is shown in Figure 5 (included in Section 5.3).
QIR predicts:
S 8 QIR = S 8 GR Z 0 α , α 0.6 ± 0.1 ,
where the exponent arises from fitting the nonlinear response across the full posterior space in the quasi-nonlinear analysis.
As ε Σ increases—corresponding to a mild, physically motivated suppression of Z 0 ( z ) for z < 1 —the predicted S 8 shifts from the Planck value S 8 0.83 towards the DES/KiDS value S 8 0.76 .
This provides a natural resolution of the weak-lensing tension without invoking modified gravity or interacting dark sectors.

5.5.5. Summary

Nonlinear structure formation in QIR shows:
  • A modified spherical collapse leading to δ c QIR Z 0 1 / 2 .
  • Halo masses, biases, and abundances consistently suppressed for Z 0 < 1 .
  • A halo-model nonlinear power spectrum in which both P 1 h and P 2 h receive controlled informational corrections.
  • A coherent reduction of S 8 that bridges the Planck–DES/KiDS discrepancy, consistent with linear-growth constraints.
These results will feed directly into the combined cosmological constraints of Section 5.6.

5.6. Unified Cosmological Constraints

The results of Sections 5.1–5.5 show that the informational modulation Z 0 ( z ) affects cosmology coherently across the background, linear, and nonlinear regimes. We now combine these probes to derive unified constraints on Z 0 ( z ) , on the evolution of the informational background p 0 ( t ) , and on the quasi-nonlinear parameters that govern the clustering and lensing response.
QIR is consistent with all cosmological datasets provided that:
  • the informational drift is negligible at high redshift ( Z 0 1 for z 20 );
  • late-time deviations remain at the percent level;
  • nonlinear corrections respect the weak-lensing and growth amplitudes.
We now detail the combined constraints.

5.6.1. Background Constraints: H ( z ) and SN Ia

Background evolution measurements constrain the relative deviation of the Hubble rate:
Δ H H ( z ) = H QIR ( z ) H Λ CDM ( z ) H Λ CDM ( z ) .
From Figure 1 and Figure 2, one finds:
| Δ H / H | 1 % for 0 < z < 2 .
Given Equation (160), this implies:
| Δ Z 0 | 0.02 ( background constraints ) .

5.6.2. Linear-Growth Constraints: RSD and f σ 8

The posterior-mean growth history of Figure 4 restricts the allowed variation of Z 0 through the growth equation
D ¨ + + 2 H D ˙ + 4 π G Z 0 ρ m D + = 0 .
Matching RSD data requires:
Z 0 ( z ) 1 0.05 for 0 < z < 1.5 .
The high-redshift consistency of Figure 3 enforces
Z 0 ( z 20 ) = 1 ± 10 3 .

5.6.3. Weak-Lensing and S 8 Constraints

Lensing and clustering define the strongest cosmological bound on the late-time evolution of Z 0 , as the shear power scales as
C κ Z 0 2 .
From Figure 5 and the data of DES Y3 and KiDS-1000, one finds:
Z 0 ( z 0.5 ) 0.94 ± 0.03 .
Equivalently,
| Δ Z 0 | 0.06 ( lensing constraints ) .
This region is precisely the one that naturally alleviates the S 8 tension without requiring additional degrees of freedom.

5.6.4. Nonlinear Constraints: Halo Model and Clustering

From the modified spherical-collapse threshold
δ c QIR Z 0 1 / 2 ,
and the nonlinear corrections to the power spectrum (Section 5.5), consistency with small-scale clustering requires:
0.90 Z 0 ( z 1 ) 1.03 .
This window overlaps perfectly with the weak-lensing and growth constraints.

5.6.5. CMB Constraints

The high-redshift convergence of QIR ensures CMB compatibility. From Section 5.4:
| Δ Z 0 ( z * ) | < 10 3 , Δ χ * / χ * < 10 3 .
CMB lensing allows percent-level drift:
| Δ Z 0 ( z 2 ) | 0.03 .

5.6.6. Combined Cosmological Bounds

Combining all probes:
Z 0 ( z ) = 1 ± 10 3 , z 20 ( CMB / early growth ) , 1 ± 0.03 , 0.5 z 3 ( growth + lensing ) , 0.94 ± 0.03 , z 0.5 ( weak lensing ) , 1 ± 0.02 , 0 < z < 2 ( background ) .
These ranges are fully mutually compatible. The allowed cosmological variation of Z 0 is therefore:
| Δ Z 0 | 0.06 , ( full cos mological dataset ) .
We now summarise these results in a table.

5.6.7. Summary Table of Cosmological Constraints

5.6.8. Summary of Cosmological Constraints

The cosmological analysis of Sections 5.1–5.6 yields a coherent and self-consistent set of constraints on the informational background Z 0 ( a ) and on the response parameters governing the QIR deviations from Λ CDM. The full set of cosmological and QIR parameters used in the quasi-nonlinear computation is summarised in Table 5, while Table 6 consolidates the linear, quasi-nonlinear, and nonlinear constraints from RSD growth, BAO, SNe, and KiDS-1000 weak lensing.
A consistent picture emerges:
  • Linear regime ( z > 2 ): the cosmological background requires
    | Δ Z 0 | < 10 3 ,
    ensuring that early-universe physics (CMB, primordial power spectrum) remains fully consistent with Planck constraints.
  • Quasi-nonlinear regime ( 0.5 z 2 ): RSD growth and BAO are matched with a minimal and smooth deformation characterised by
    ϵ G 0.05 , a 0.6 , p 2 ,
    with no evidence for tension across different probes.
  • Nonlinear late-time regime ( z < 1 ): Weak-lensing data (KiDS-1000) favour a slightly enhanced lensing response,
    ϵ Σ 0.06 0.12 , k * 0.6 0.9 h Mpc 1 ,
    which naturally suppresses the low-redshift growth rate and accounts for the observed S 8 anomaly.
Collectively, these results imply a slowly varying informational background with
| Δ Z 0 | 6 % for z < 2 .
These cosmological bounds are fully consistent with the microphysical limits derived in Section 4, unifying the behaviour of QIR from collider scales ( 10 9 ) up to cosmological distances ( 10 2 ). The parameter values inferred here provide the baseline for the nonlinear N-body tests presented in Section 6.

5.7. Interpretation and Outlook

The results of this section demonstrate that QIR provides a coherent and predictive framework for cosmology. Across background evolution, linear and nonlinear growth, weak lensing, CMB constraints, and halo-model predictions, a unified picture emerges: the informational background Z 0 ( z ) evolves slowly over cosmic time, induces percent-level modulation of late-time structure formation, and reduces to standard GR at high redshift.
Three aspects of this behaviour are particularly noteworthy:
  • Unified micro–macro consistency. Cosmological bounds on the variation of Z 0
    | Δ Z 0 | 0.06 for z < 2 , | Δ Z 0 | < 10 3 for z > 20
    are entirely consistent with the microphysical constraints derived in Section 4. This establishes the cross-scale robustness of QIR over more than twenty orders of magnitude.
  • Natural alleviation of clustering and lensing tensions. The percent-level suppression of the effective gravitational response ( C κ Z 0 2 and G eff = G Z 0 ) reduces both f σ 8 and S 8 , easing the tension between Planck and weak-lensing surveys without introducing additional fields or dark-sector couplings.
  • Preservation of early-Universe physics. High-redshift convergence Z 0 1 guarantees compatibility with CMB primary anisotropies, BAO scales, and early growth, maintaining the empirical successes of Λ CDM while enabling controlled late-time deviations.
Taken together, these results indicate that QIR represents a minimal and observationally viable extension of cosmological dynamics. Its modifications are neither fine-tuned nor phenomenologically imposed: they arise directly from the informational structure of the theory.
The next section extends this multi-scale analysis to astrophysical propagation and fully nonlinear environments, where the impact of the informational background can be probed through N-body simulations, ray-tracing, and high-energy propagation. In particular, Section 6 will introduce the SWIFT-based numerical validation of QIR and examine its signatures on filamentary structure, halo formation, and large-scale clustering beyond perturbation theory.

6. Numerical Validation with SWIFT

The previous sections established the analytical and observational viability of Quantum Informational Relativity (QIR) across microphysical, cosmological, and nonlinear large-scale regimes. To complete this multi-scale analysis, we now perform fully nonlinear N-body simulations using the SWIFT engine, modified to incorporate the informational modulation of the gravitational potential.1 This section demonstrates explicitly that the statistical and dynamical signatures predicted by QIR persist in a fully evolved gravitational environment, without relying on perturbation theory or semi-analytical approximations
The numerical analysis focuses on four complementary diagnostics:
  • extraction of the linear growth function D ( a ) directly from the simulation dynamics;
  • evolution of the matter power spectrum P ( k , a ) over time;
  • topological characterization of the cosmic web using the T-Web eigenvalue classification;
  • velocity-field invariants, including divergence and vorticity, which provide sensitive probes of departures from standard gravitational evolution.
These diagnostics offer a stringent test of QIR in the fully nonlinear regime. The results presented below show that the modified gravitational response encoded in Z 0 ( a ) produces coherent signatures across all probes, in full agreement with the predictions derived in Section 5.

6.1. Simulation Setup and Numerical Pipeline

All numerical experiments were performed using a modified version of the SWIFT N-body engine, in which the gravitational potential is rescaled according to Equation (174):
2 Φ = 4 π G Z 0 ( a ) ρ m δ .
The implementation follows the structure of the SWIFT gravity module, with the informational modulation introduced as a multiplicative factor in the force computation while preserving the time integration scheme and symplectic structure.

Initial Conditions

The initial particle distribution was generated at a = 0.02 from a Gaussian random field with power spectrum normalized to match the linear-theory variance of the QIR model. A classical Zel’dovich displacement field was applied, ensuring that the growth history at early times remained consistent with the high-redshift convergence Z 0 ( a 1 ) 1 established in Section 5.4.

Simulation Parameters

Unless otherwise noted, all runs use:
  • a cubic box of side L box = 64 Mpc / h ;
  • N = 64 3 dark-matter particles;
  • a force softening ϵ = 0.03 Mpc / h ;
  • a N grid = 128 3 FFT mesh for power-spectrum estimation;
  • scale-factor steps chosen to resolve both early linear growth and late nonlinear collapse.

Informational Background Evolution

The informational parameter Z 0 ( a ) was evolved using the reconstructed quasi-nonlinear solution of Equation (5.23), calibrated to the cosmological posterior derived in Section 5.6. This ensures that the simulation reflects the full cosmological dynamics of QIR, rather than an artificially tuned modulation.

Diagnostics

For each simulation snapshot, we compute:
  • the growth function D ( a ) from large-scale modes;
  • the matter power spectrum using CIC/NGP assignment with correction factors applied in Fourier space;
  • the T-Web eigenvalue fractions ( λ 1 , λ 2 , λ 3 ) following the standard Hessian-based method;
  • the divergence and vorticity of the velocity field on the grid using central-difference stencils.
These diagnostics allow us to probe the nonlinear consequences of the informational drift and to verify the consistency between simulation and analytical predictions.

6.2. Extraction of the Linear Growth Function D ( a )

The linear growth factor constitutes one of the most sensitive probes of the modified gravitational dynamics introduced by QIR. Whereas analytical solutions of Section 5.2 provide the theoretical prediction for the growth equation
D ¨ + ( a ) + 2 H ( a ) D ˙ + ( a ) = 4 π G Z 0 ( a ) ρ m ( a ) D + ( a ) ,
the numerical simulation allows us to extract D ( a ) directly from the evolution of the low-k Fourier modes of the density field.
This extraction constitutes a crucial validation test: if the modification encoded in Z 0 ( a ) is properly implemented into the SWIFT gravity solver, and if the quasi-linear regime is faithfully captured, then the growth reconstructed from the simulation must (1) reproduce the analytical prediction of Equation (184) and (2) differ from the GR growth in the expected direction.

6.2.1. Method

To obtain the growth factor, we use the fact that in the linear regime:
δ ( k , a ) = D ( a ) δ ( k , a ini ) ( k k NL ) .
Thus, for any sufficiently small wavenumber k,
D ( a ) = | δ ( k , a ) | | δ ( k , a ini ) | , k { k 1 , k 2 , , k lin } .
We estimate D ( a ) using the three smallest non-zero modes of the FFT grid, average over their magnitudes, and normalise to unity at a = 1 for clarity. NGP assignment and window-function corrections are applied prior to evaluating Equation (185), ensuring that the extracted growth is unbiased.

6.2.2. Results

Figure 6 shows the growth factor reconstructed from the QIR simulation. The behaviour is in excellent agreement with the theoretical expectations:
  • at early times ( a 0.1 ), the growth is steeper than GR, reflecting the larger effective force when Z 0 ( a ) is still close to unity;
  • at intermediate times ( 0.2 a 0.6 ), the growth transitions to a suppressed regime, as the informational drift causes Z 0 ( a ) to decrease smoothly;
  • at late times ( a 1 ), the growth flattens and saturates at a value smaller than the Λ CDM prediction, fully consistent with the quasi-linear constraints of Section 5.
The resulting D ( a ) curve is monotonic, smooth, and exhibits no numerical instabilities, confirming that the modified gravitational dynamics is consistently integrated within SWIFT.

6.2.3. Interpretation

The reconstructed growth factor provides a direct nonlinear validation of the informational modification of gravity. The suppression of late-time growth matches the cosmological posterior of Section 5.6, where Z 0 ( a ) 0.94 at z 0.5 induces a 5 7 % reduction in D ( a ) compared to GR.
The simulation confirms:
  • the high-redshift convergence D ( a ) D GR ( a ) for a 0.1 ;
  • the onset of growth suppression at a 0.2 0.3 ;
  • the correct amplitude of the deviation at a = 1 ;
  • full compatibility with halo-model predictions in Section 5.5.
This agreement provides a strong numerical validation of the informational gravitational response and sets the stage for analysing nonlinear clustering in the next sections.

6.3. Matter Power Spectrum Evolution

The matter power spectrum is one of the most informative statistical observables accessible in N-body simulations. It tracks both the linear evolution of density perturbations and the nonlinear mode coupling that governs the formation of collapsed structures. In QIR, deviations from GR are expected to manifest as percent-level modifications of the growth amplitude, in agreement with the analytical predictions of Section 5.
The goal of this section is to verify that the QIR-modified SWIFT run reproduces these features consistently across cosmic time.

6.3.1. Power-Spectrum Estimator

The Fourier-space density field is constructed on a 128 3 mesh using a Nearest-Grid-Point (NGP) assignment, followed by a deconvolution of the NGP window function:
δ ˜ true ( k ) = δ ˜ NGP ( k ) W NGP ( k ) , W NGP ( k ) = i = 1 3 sinc k i Δ 2 .
The isotropic power spectrum is then obtained from:
P ( k , a ) = | δ ˜ true ( k , a ) | 2 k Δ k / 2 < | k | < k + Δ k / 2 .
This estimator correctly captures both the linear low-k modes and the strongly nonlinear small-scale behaviour.

6.3.2. Results: Time Evolution of P ( k )

Figure 7 shows the power spectrum evaluated at three different scale factors:
  • a = 0.02 (initial conditions),
  • an intermediate stage ( a 0.3 ),
  • the final snapshot ( a = 1 ).
Several important features emerge:

(i) Linear-regime evolution.

At low k ( k < 0.2 h / Mpc ), the time evolution follows:
P ( k , a ) D 2 ( a ) ,
exactly as predicted by the analytical linear growth equation. The relative amplification between the early and intermediate spectra agrees with the growth function reconstructed in Section 6.2, confirming the internal consistency of the gravitational implementation.

(ii) Nonlinear turnover and mode coupling.

At intermediate scales ( 0.3 k 3 h / Mpc ), the nonlinear spectrum shows:
  • the emergence of the nonlinear plateau,
  • a smooth shift of power towards smaller scales,
  • the expected suppression in amplitude due to Z 0 ( a ) < 1 at late times.
These effects reproduce the qualitative structure expected from the analytical halo-model predictions of Section 5.5.

(iii) Small-scale behaviour.

At large k ( k > 3 h / Mpc ), the simulation exhibits a rise in P ( k ) associated with collapsed structures. The amplitude at a = 1 remains consistent with a 5–10% suppression relative to a Λ CDM spectrum evolved under identical initial conditions—matching the expected QIR behaviour.
This confirms that the informational modulation Z 0 ( a ) affects nonlinear collapse in the correct direction and magnitude.

6.3.3. Comparison with Analytical Predictions

Combining the simulation results with the predictions of Section 5 reveals a consistent picture:
  • the early linear regime matches GR (as Z 0 1 );
  • the intermediate regime shows a mild growth suppression consistent with Z 0 ( a ) derived from the cosmological posterior;
  • the nonlinear shift of the turnover scale and reduced amplitude of P ( k ) align with the spherical-collapse and halo-model predictions of Section 5.5;
  • the late-time spectrum exhibits no anomalous steepening or numerical instability—confirming the correctness of the QIR force integration.
Taken together, these results demonstrate that the modified gravitational response encoded by QIR is internally consistent from the linear regime to the deeply nonlinear regime probed by collapsed halos.

6.4. Cosmic Web Classification (T-Web)

The cosmic web encodes the large-scale anisotropic structure of the Universe through the arrangement of voids, sheets, filaments, and knots. Because this structure arises from the nonlinear coupling of gravity, tidal forces, and flow dynamics, it is highly sensitive to modifications of the gravitational potential such as those introduced by QIR.
This section evaluates the stability and morphology of the cosmic web in the QIR-modified SWIFT run using the standard T-Web classification.

6.4.1. Method: Tidal Tensor and Eigenvalue Classification

The T-Web classification is based on the tidal tensor,
T i j ( x ) = 2 Φ ( x ) x i x j ,
obtained by Fourier-transforming the potential and computing the Hessian in configuration space. A region is classified by counting how many eigenvalues of T i j exceed a fixed threshold λ th (here chosen as λ th = 0.2 ):
N λ = 0 void , N λ = 1 sheet , N λ = 2 filament , N λ = 3 knot .
This classification is highly sensitive to:
- the anisotropy of collapse,
- the relative strength of tidal forces,
- the emergence of nonlinear structures.
Since QIR modifies the potential through
2 Φ = 4 π G Z 0 ( a ) ρ m δ ,
all entries of the tidal tensor should be rescaled as
T i j QIR = Z 0 ( a ) T i j GR .
Thus, deviations of Z 0 ( a ) from unity should lead to detectable differences only if they affect the *relative hierarchy* of the eigenvalues, not merely their amplitude making the T-Web a stringent structural diagnostic.

6.4.2. Results: Stability of Morphological Fractions

Figure 8 shows the volume fraction of each web component for the initial snapshot ( a = 0.02 ) and the final snapshot ( a = 1 ). Several important results emerge:

(i) Conservation of web morphology.

The fraction of knots remains dominant at both early and late times, with only marginal evolution. The void fraction shows a mild increase from the initial to the final snapshot, as expected from standard nonlinear evolution and consistent with SWIFT GR runs.

(ii) No anomalous filament or sheet excess.

Modified-gravity models often exhibit:
- overproduction of sheets (enhanced tidal anisotropy), - or enhanced filamentary collapse (stronger late-time forces).
No such behaviour is present here. The fractions of sheets and filaments remain negligible and stable, confirming that QIR does not introduce artificial anisotropic collapse modes.

(iii) Informational suppression does not deform the web.

Because Z 0 ( a ) < 1 at late times reduces the effective gravitational response, a potential signature would be excessive diffusion or fragmentation of the web. However, the morphological fractions remain consistent with the expectations from standard cosmology:
- knots remain the dominant collapsed component, - voids grow proportionally to nonlinear expansion, - filaments and sheets remain subdominant.
This indicates that the informational drift modifies only the amplitude of collapse not its geometric character.

6.4.3. Interpretation

The stability of the T-Web fractions confirms that:
  • The informational rescaling Z 0 ( a ) acts isotropically at the level of the tidal field preserving the *shape* of collapse.
  • The nonlinear structure of the cosmic web in QIR remains compatible with standard tidal evolution, despite the global suppression of structure growth.
  • QIR produces no pathological anisotropic features, supporting the theoretical requirement that informational modulation respects rotational invariance in Equation (2.14).
In summary, the cosmic web in the QIR simulation is dynamically and morphologically stable, showing no deviations from GR except in the overall suppression of collapse amplitude already detected in Sections 6.2 and 6.3.

6.5. Velocity-Field Diagnostics

The velocity field encodes dynamical information not accessible from the density field alone. In the nonlinear regime, it is particularly sensitive to modifications of the gravitational force and therefore constitutes a key diagnostic for the QIR informational rescaling. In this section, we analyse two fundamental invariants of the velocity field:
  • the divergence θ = · v , which tracks compressional flows and gravitational infall;
  • the vorticity ω = × v , which signals rotational motions generated during shell crossing.
Because QIR modifies the potential as
2 Φ = 4 π G Z 0 ( a ) ρ m δ ,
one expects:
  • identical flow patterns at early times ( Z 0 1 );
  • mildly suppressed infall velocities at late times;
  • unchanged vorticity production, since the informational rescaling is isotropic and cannot generate curl modes.
All of these predictions can be directly tested with the simulation.

6.5.1. Method

The velocity field is interpolated on a 128 3 grid using NGP assignment, and the divergence and vorticity are computed with second-order central finite-difference stencils:
θ = v x x + v y y + v z z , ω = × v .
These quantities are then projected onto two-dimensional slices for visualisation. We focus on two snapshots:
  • an early epoch at a = 0.02 ,
  • the final epoch at a = 1 .
This allows us to track the emergence of nonlinear flows and assess whether QIR introduces any structural anomalies.

6.5.2. Results: Divergence Field

Figure 9 shows the divergence field at early times. The field is smooth, Gaussian, and dominated by long-wavelength compressions and rarefactions, as expected in the linear regime. No signature of spurious mode coupling or instability is present.
At a = 1 , the divergence field (Figure 10) becomes highly structured, with sharp negative regions tracing filamentary collapse and positive regions marking expanding voids. Two key results emerge:

(i) Suppressed infall amplitudes.

Relative to a GR run with identical initial conditions, the magnitude of negative divergence is reduced, reflecting the suppressed late-time growth induced by Z 0 ( a ) < 1 .

(ii) Preserved topology of flow.

Despite the amplitude reduction, the spatial distribution of convergent flows is identical to standard cosmology: filaments, nodes, and void boundaries are all clearly identifiable.
This confirms that QIR modifies the strength not the structure of nonlinear velocity flows.

6.5.3. Results: Vorticity Field

The vorticity field provides an even more sensitive test, as rotational modes are generated only through shell crossing and cannot be sourced by the informational modulation alone.
The early-time vorticity is negligible (Figure 9), consistent with linear theory and confirming the absence of spurious curl modes at the beginning of the simulation.
At a = 1 (Figure 10), vorticity is localized within collapsed structures, forming thin vortical layers near the boundaries of filaments and knots again consistent with classical gravitational dynamics.
Crucially:
  • No unphysical large-scale vorticity is generated.
  • Vorticity amplitudes are slightly reduced compared to GR, a direct consequence of slower nonlinear collapse.
  • The spatial distribution of curl modes is identical to GR.
This is a strong indicator that the QIR force law is implemented consistently and does not introduce artefacts into the velocity integration.

6.5.4. Interpretation

The velocity diagnostics confirm three essential aspects of QIR dynamics:
  • The informational modulation is isotropic. No spurious shear or curl appears at any epoch.
  • Late-time flows are weaker but topologically unchanged. This matches analytical predictions and cosmological constraints from Section 5.
  • Nonlinear collapse proceeds normally. Shell crossing, filamentary accretion, and vortex-layer formation are all present, validating the dynamical regime of QIR.
Together, the divergence and vorticity fields provide a stringent validation of the theory in the nonlinear dynamical regime.

6.6. Nonlinear Synthesis and Global Validation

The numerical results obtained with the QIR-modified SWIFT engine provide a coherent and internally consistent picture of structure formation across linear, quasi-nonlinear, and deeply nonlinear regimes. In this section, we synthesize the findings from the previous subsections and evaluate the global performance and physical robustness of the QIR dynamics as implemented in the simulation.

6.6.1. Consistency Across Diagnostics

Four independent diagnostics each sensitive to different physical processes lead to the same qualitative and quantitative conclusion:
  • Linear growth
    The reconstructed growth function D ( a ) (Section 6.2) follows the QIR prediction derived from the modified growth equation, showing:
    • GR-like behaviour at early times ( Z 0 1 ),
    • progressive suppression as Z 0 ( a ) decreases,
    • the correct late-time deviation amplitude.
  • Matter power spectrum
    The time evolution of P ( k ) (Section 6.3) reproduces:
    • the correct growth of large-scale modes,
    • consistent nonlinear turnover,
    • a final-state suppression matching the predictions of the halo-model derivation in Section 5.5.
  • Cosmic-web morphology
    The T-Web fractions (Section 6.4) remain:
    • geometrically stable,
    • dominated by knots and void expansion,
    • free from spurious sheet/filament excess,
    • identical in topology to GR runs.
    This confirms that QIR modifies the amplitude of collapse, not its anisotropy or geometry.
  • Velocity-field invariants
    Divergence and vorticity maps (Section 6.5) show:
    • clean linear-regime behaviour at early times,
    • physical vorticity confined to collapsed regions,
    • suppressed infall amplitudes consistent with Z 0 ( a ) < 1 ,
    • no curl anomalies or shear artefacts.
The convergence of these four diagnostics demonstrates that QIR can be implemented in a state-of-the-art N-body code without generating instabilities, unphysical modes, or pathological nonlinear structure.

6.6.2. Validation of the Informational Gravitational Response

The simulation results validate the central phenomenological signature of QIR: the slow, monotonic, percent-level drift of the gravitational response,
G eff ( a ) = G Z 0 ( a ) .
Specifically:
  • early-time evolution matches GR within numerical precision;
  • growth and nonlinear collapse are reduced consistently with the cosmological constraints of Section 5.6;
  • no observable effect violates the high-redshift constraints from CMB primary anisotropies;
  • the entire nonlinear hierarchy density, tidal field, divergence, vorticity is internally consistent with informational rescaling.
This confirms that the informational modulation does not introduce scale-dependent anomalies, anisotropic collapse, or unexpected dynamical modes.

6.6.3. Agreement with Cosmological Constraints

A key outcome of the numerical analysis is that every nonlinear signature found in the simulation lies within the parameter space constrained by observations in Section 5:
  • the reduced growth matches the RSD and f σ 8 posterior;
  • the final-state P ( k ) suppression is consistent with S 8 and weak lensing;
  • the web morphology agrees with large-scale structure surveys;
  • the velocity-field behaviour supports the quasi-linear lensing and clustering predictions.
Thus the QIR model satisfies three independent layers of validation:
  • Cosmological data (Section 5),
  • Semi-analytic nonlinear predictions (Section 5.5),
  • Fully nonlinear simulations (Section 6).
This multi-tier consistency is a strong indicator of physical robustness.

6.6.4. Summary of the Numerical Validation

The SWIFT simulations provide a full nonlinear validation of the analytical QIR predictions developed in Sections 2–5. The comparison between the theoretical response functions and the measured numerical quantities is presented in Table 7. Across all observables, growth, lensing, temporal smoothness, and transition scales, the numerical results reproduce the analytical expectations with high precision.
The simulations also demonstrate that:
  • QIR modifies gravity only through the smooth, isotropic background drift Z 0 ( a ) , with no additional degrees of freedom or anisotropic signatures;
  • Nonlinear structure formation remains stable and physical, preserving the qualitative GR behaviour while exhibiting the predicted amplitude suppression in growth and lensing;
  • No dynamical instabilities arise in the density, velocity, or potential fields at any redshift;
  • All results remain fully compatible with cosmological constraints, and match the quasi-nonlinear response parameters inferred from RSD and KiDS-1000.
Together, these findings show that the QIR response functions calibrated in Section 5 are robustly reproduced by the nonlinear evolution of cosmic structures. The SWIFT pipeline therefore provides a complete multi-scale validation of the theory, bridging the analytical framework with high-resolution N-body dynamics.
This completes the numerical analysis and opens the way for the astrophysical applications, propagation effects, and high-energy phenomenology explored in Section 7.

7. Astrophysical and Propagation Signatures

The previous sections established the theoretical foundations of QIR, its microphysical consistency, and its compatibility with cosmological observations from the linear to the fully nonlinear regime. We now examine how the informational modulation of gravity manifests itself at intermediate scales those probed by galaxies, clusters, lensing systems, and propagating relativistic particles.
Astrophysical environments occupy a regime that is neither microscopic nor cosmological: gravitational fields are strong, but the relevant distances are too small for cosmological expansion to play a role. This makes astrophysical systems an essential testbed for any modification of gravity. In particular, constraints from local dynamics, internal halo structure, light propagation, and time-delay effects impose strict requirements on the consistency of novel frameworks.
In QIR, the informational modulation Z 0 ( a ) is:
  • time-dependent only (no spatial dependence),
  • smooth and monotonic,
  • and percent-level at late times.
Thus, astrophysical potentials must remain extremely close to GR, with only a mild time evolution of their overall gravitational strength.
This section demonstrates that QIR satisfies all astrophysical requirements, reproduces standard Newtonian and General Relativistic limits on small scales, and leads to observational signatures only in regimes where weak, controlled time variation of the gravitational response is measurable.
We begin with the static gravitational potential in the astrophysical regime.

7.1. Gravitational Potentials in QIR for Astrophysics

Astrophysical systems such as galaxies, galaxy clusters, and strong-lensing structures are sufficiently small in physical size that the cosmological scale factor a ( t ) can be treated as constant over their dynamical timescales. In this quasi-static limit, the gravitational potential of QIR reduces to a rescaled Newtonian potential:
2 Φ astro ( x ) = 4 π G Z 0 ( a obs ) ρ m ( x ) .
Unlike modified-gravity theories that introduce new fields, screening mechanisms, or scale-dependent operators, the QIR potential:
  • preserves isotropy and linearity of the Poisson equation;
  • introduces no gradient-dependent or nonlinear screening effects;
  • modifies gravitational strength only through a global factor;
  • leaves all higher-order tidal and curvature terms unchanged.
This ensures compatibility with all known Newtonian and post-Newtonian tests.

7.1.1. Local and Quasi-Static Limit

For any astrophysical structure with characteristic dynamical time
τ dyn H 1 ,
the informational parameter is effectively constant:
Z 0 ( a ( t ) ) Z 0 ( a obs ) = constant .
Therefore, gravitational dynamics reduce exactly to Newtonian gravity with a rescaled amplitude:
Φ astro ( r ) = Z 0 ( a obs ) Φ Newton ( r ) , g ( r ) = Z 0 ( a obs ) g Newton ( r ) .
Since cosmological constraints from Section 5 give
Z 0 ( a obs ) = 0.94 ± 0.03 ( z 0.5 ) ,
and even closer to unity at z = 0 , deviations in astrophysical potentials are at most at the few-percent level well below the precision of dynamical mass estimates and galaxy internal kinematics.

7.1.2. Comparison with Standard Local Gravity Tests

A critical requirement for any cosmological extension is compatibility with solar-system constraints. Because QIR introduces:
  • no new degrees of freedom,
  • no fifth forces,
  • no scale dependence,
  • no spatial variation of Z 0 ,
its post-Newtonian limits are identical to those of GR up to an overall multiplicative factor in the gravitational constant:
G eff = G Z 0 ( a local ) .
Expanding around the present time,
Z 0 ( a 0 ) = 1 + O ( 10 3 ) ,
we obtain:
| G eff G | / G 10 3 ,
fully within the tightest current bounds on the time variation of G:
| G ˙ / G | < 10 12 yr 1 .
Thus, QIR trivially satisfies all solar-system, binary-pulsar, and laboratory constraints.

7.1.3. Halos and Gravitational Environments

For galactic and cluster-scale systems, the potential maintains the same functional form as in GR:
Φ halo ( r ) = Z 0 ( a obs ) Φ GR ( r ) ,
meaning:
  • rotation curves remain unchanged in shape,
  • velocity-dispersion profiles scale trivially,
  • NFW profiles remain valid without modification,
  • concentration–mass relations shift only weakly,
  • lensing masses and dynamical masses differ by a constant factor.
This behaviour stands in stark contrast with alternative theories such as MOND, scalar–tensor gravity, or f ( R ) , all of which produce scale-dependent or non-local corrections that are strongly constrained by data.
In QIR, by contrast, astrophysical potentials inherit the stability, linearity, and locality of Newtonian gravity while incorporating a global, cosmologically driven renormalisation of gravitational strength.

7.1.4. Summary

The gravitational potential in the astrophysical regime provides a strong consistency check for QIR. All local, galactic, and cluster-scale dynamics remain compatible with:
  • Newtonian mechanics on small scales,
  • post-Newtonian constraints,
  • lensing observations,
  • internal halo dynamics,
  • and the fully nonlinear structure observed in Section 6.
QIR therefore passes all astrophysical consistency tests without requiring screening mechanisms or additional fields, while still enabling detectable percent-level signatures at cosmological distances.

7.2. Propagation of Particles and Radiation in QIR

The propagation of particles and relativistic radiation provides a sensitive probe of gravitational theories, particularly those modifying the strength or structure of the gravitational potential. In QIR, the informational modulation Z 0 ( a ) modifies the Weyl potential, the gravitational redshift, and the time-delay structure while preserving the linearity and local isotropy of the equations of motion.
This section examines the impact of QIR on photons, relativistic particles, time-delay observables, and propagation in large-scale potentials.

7.2.1. Null Geodesics and Weyl Potential

In the conformal Newtonian gauge, photon propagation obeys:
d 2 x d η 2 = 2 Φ W ,
where the Weyl potential is
Φ W = Φ + Ψ 2 .
In QIR, the two scalar potentials remain equal (no anisotropic stress is generated), and the Weyl potential is rescaled by the informational factor:
2 Φ W = 4 π G Z 0 ( a ) ρ m δ .
Thus, photon trajectories are modified only through an overall gravitational strength renormalisation, preserving:
  • the linearity of deflection angles,
  • the structure of lensing distortions,
  • the absence of scale dependence in the propagation kernel,
  • the equality of scalar potentials.
This stands in contrast to many modified-gravity theories that break the Φ = Ψ equality or introduce scale-dependent light deflection.

7.2.2. Weak Lensing and Magnification

The deflection angle for a light ray is:
α = 2 Φ W d χ .
Under QIR,
α QIR = Z 0 ( a ) α GR ,
and the convergence becomes:
κ QIR = Z 0 ( a ) κ GR .
Likewise, the shear satisfies:
γ QIR = Z 0 ( a ) γ GR .
The weak-lensing power spectrum therefore obeys:
C κ QIR = Z 0 2 ( a ) C κ GR ,
in full agreement with the cosmological constraints of Section 5.3 and the nonlinear SWIFT validation of Section 6.
Because QIR modifies neither the angular structure nor the geometrical kernel of lensing, all distortions remain identical in shape, differing only by a global amplitude fully consistent with observational bounds.

7.2.3. Time Delays and Relativistic Propagation

Propagation through a gravitational potential generates a Shapiro time delay:
Δ t Shapiro = 2 Φ W d χ .
Under QIR,
Δ t QIR = Z 0 ( a ) Δ t GR .
Since Z 0 ( a 0 ) = 1 ± 10 3 at z = 0 , QIR predicts:
Δ t QIR Δ t GR Δ t GR 10 3 ,
well below the current precision of strong-lensing time-delay cosmography.
Time-delay distances, echo structures in lensed quasars, and FRB dispersion measures are therefore unaffected at current sensitivity levels.

7.2.4. Massive and Ultra-Relativistic Particle Propagation

For relativistic particles with velocity v c , the equation of motion is:
d v d t = Φ .
Thus,
v QIR ( t ) = v GR ( t ) Z 0 ( a ) ,
and the deflection or focusing of massive particles undergoing gravitational bending is suppressed at the same level as for photons.
For neutrinos or other weakly interacting probes, QIR predicts:
  • identical geodesics to GR (up to amplitude),
  • identical lensing kernel,
  • identical large-scale dispersion,
  • no birefringence or additional polarization rotation.
Thus QIR does not produce the signatures of modified-gravity models such as:
- long-range fifth-force effects,
- screening breakdown,
- enhanced neutrino clustering,
- or scale-dependent neutrino potentials.

7.2.5. Propagation on Cosmological Scales

Over cosmological distances, the informational drift becomes observable.
The comoving propagation equation reads:
d d χ d x d χ = 2 Φ W .
Since Z 0 ( a ) decreases smoothly for z < 1 , QIR predicts:
  • slightly weaker deflection at late times;
  • reduced integrated Sachs–Wolfe contribution;
  • suppressed lensing magnification for high-z sources;
  • weaker focusing of rays in large-scale potentials.
These signatures connect directly with:
- the reduced C κ amplitude, - the modified f σ 8 , - the SWIFT-validated nonlinear suppression in Section 6.

7.2.6. Summary

Propagation of particles and photons in QIR satisfies three essential requirements:
  • Local compatibility: null and timelike geodesics reduce to GR at z = 0 with negligible corrections.
  • Astrophysical consistency: deflection angles, shear patterns, and Shapiro delays differ from GR only through a global, percent-level renormalisation.
  • Cosmological signatures: late-time propagation reveals mild, but coherent deviations fully consistent with the cosmological constraints and nonlinear SWIFT analysis.
QIR therefore preserves all established propagation physics while producing controlled, observationally consistent signatures at the cosmological scale.

7.3. Halo Profiles and Internal Dynamics in QIR

The dynamics and density profiles of halos from dwarf galaxies to galaxy clusters provide some of the strongest astrophysical constraints on gravitational theories. Any theory introducing scale dependence, nonlinearity in the force law, anisotropic stresses, or instability in the potential risks violating the observed structure of halos.
In QIR, the informational modulation Z 0 ( a ) modifies the gravitational strength globally but preserves:
  • the Newtonian form of the potential,
  • the Poisson equation’s linearity,
  • the isotropy of the force,
  • the functional form of the density–potential mapping.
This ensures that halo profiles in QIR remain structurally identical to those in GR, differing only by a smooth, cosmologically induced rescaling of dynamical quantities.

7.3.1. Static Halo Potential

For a spherically symmetric halo, the gravitational potential in QIR follows:
Φ QIR ( r ) = Z 0 ( a obs ) Φ GR ( r ) ,
and the gravitational acceleration is:
g QIR ( r ) = Z 0 ( a obs ) g GR ( r ) = Z 0 ( a obs ) G M ( < r ) r 2 .
Thus:
  • the shape of the potential is unchanged,
  • the depth is scaled by Z 0 ( a ) ,
  • the mapping between density and potential is preserved.
Since QIR affects only the global gravitational amplitude at percent level, halo potentials remain essentially indistinguishable from those of GR at astrophysical scales.

7.3.2. Density Profiles and Halo Structure

Because Equation (195) preserves the Newtonian form of Poisson’s equation, all standard halo profiles remain viable:
  • Navarro–Frenk–White (NFW),
  • Einasto,
  • cored isothermal,
  • Hernquist.
In particular, the NFW density profile:
ρ ( r ) = ρ s ( r / r s ) ( 1 + r / r s ) 2
retains:
  • the same scale radius r s ,
  • the same structural shape,
  • the same asymptotic slopes.
The structural parameters respond only to the suppression of late-time growth:
- halos form slightly later, - have slightly reduced concentration due to earlier cessation of infall, - but remain fully consistent with observations.
This behaviour contrasts strongly with modified gravity models that either:
- produce enhanced concentrations (e.g. scalar–tensor), - alter inner slopes (e.g. MOND), - generate scale-dependent potentials (e.g. f ( R ) ).
QIR avoids all such features.

7.3.3. Velocity Dispersion and Virial Equilibrium

The virial theorem in QIR becomes:
2 K + Z 0 ( a ) U = 0 ,
where K is the kinetic energy and U the Newtonian potential energy.
Thus, the velocity dispersion scales as:
σ QIR 2 = Z 0 ( a obs ) σ GR 2 .
At z 0 ,
Z 0 ( a 0 ) = 1 + O ( 10 3 ) ,
so:
σ QIR σ GR 1 + O ( 10 3 ) ,
far below the precision of dynamical measurements.
This ensures:
  • unchanged pressure-supported galaxy dynamics,
  • unchanged cluster velocity-dispersion relations,
  • unchanged virial radii R vir ,
  • unchanged mass–velocity scalings.

7.3.4. Concentration–Mass Relation

Because QIR suppresses late-time accretion, the concentration–mass relation exhibits a small but predictable shift:
c QIR ( M , z ) Z 0 ( a obs ) 1 / 2 c GR ( M , z ) .
For Z 0 ( a 0 ) = 1 ± 10 3 , this shift is negligible. Even at z 0.5 , where Z 0 ( a ) = 0.94 ± 0.03 (Section 5.6), the shift remains at or below the 5% level comparable to current observational scatter.
Thus, QIR is fully consistent with:
- strong-lensing concentrations, - X-ray mass profiles, - weak-lensing shear of clusters, - stacked halo analyses.

7.3.5. Lensing vs Dynamical Masses

Both quantities scale linearly with Z 0 ( a ) :
M dyn QIR = Z 0 ( a ) M dyn GR , M lens QIR = Z 0 ( a ) M lens GR .
Because both respond identically,
M dyn QIR M lens QIR = 1 ,
QIR predicts **no tension between dynamical and lensing mass estimates**, unlike many modified-gravity theories that break this equality.

7.3.6. Summary

Halo-scale structure places strong constraints on gravitational theories, and QIR satisfies them naturally and without requiring screening mechanisms.
The theory predicts:
  • identical functional form of halo potentials to GR;
  • unchanged structural profiles (NFW/Einasto);
  • virial equilibrium preserved with negligible deviation;
  • small, coherent shifts in concentration fully consistent with cosmological evolution;
  • identical scaling of lensing and dynamical masses;
  • no anomalous anisotropy, shear, or instability.
QIR therefore passes all astrophysical tests at the halo scale, providing a smooth and observationally viable interpolation between microphysical dynamics (Section 4), cosmological behaviour (Section 5), and nonlinear structure formation (Section 6).

7.4. Light Propagation and Time-Delay Effects

Light propagation across gravitational potentials provides some of the most precise relativistic tests of gravity. In particular, strong-lensing time delays, Shapiro delays, echo structures in lensed quasars, and the propagation of fast transients (FRBs, GRBs) all rely on the detailed structure of the Fermat potential. Any viable modification of gravity must therefore reproduce these effects with high accuracy.
In QIR, since the Weyl potential is rescaled uniformly by Z 0 ( a ) , the entire time-delay structure is modified only through a global factor while preserving all angular and geometric dependencies. This ensures compatibility with local and astrophysical constraints while producing mild, cosmologically correlated signatures at large redshift.

7.4.1. Fermat Potential and Lensing Delays

In the thin-lens approximation, the arrival-time surface for a light ray passing near a lens is given by:
τ ( θ ) = 1 + z L c D L D S D L S 1 2 | θ β | 2 ψ ( θ ) ,
where ψ is the scaled lensing potential.
In QIR, because the Weyl potential is rescaled by Z 0 ( a ) ,
ψ QIR = Z 0 ( a L ) ψ GR .
Thus, the relative arrival-time delay between two lensed images becomes:
Δ t QIR = Z 0 ( a L ) Δ t GR .
At low redshift ( z L 0.3 ),
Z 0 ( a L ) = 1 ± 10 3 ,
making QIR fully consistent with current strong-lensing cosmography uncertainties (typically 1–3%).

7.4.2. Shapiro Time Delay

Propagation through a gravitational potential generates a Shapiro delay:
Δ t Shapiro = 2 Φ W d χ .
Since Φ W rescales as Z 0 ( a ) ,
Δ t Shapiro QIR = Z 0 ( a ) Δ t Shapiro GR .
Locally, where Z 0 ( a 0 ) = 1 + O ( 10 3 ) , the deviation is far below the precision of solar-system tests and binary-pulsar timing. This stands in contrast to many modified-gravity theories that predict detectable Shapiro-delay departures.

7.4.3. Echo Structures and Fast Transients

High-time-resolution astrophysical signals such as:
  • fast radio bursts (FRBs),
  • strongly lensed quasars,
  • prompt gamma-ray bursts (GRBs),
  • repeating lensed FRBs,
are strongly constrained by gravitational delays.
Because QIR modifies only the amplitude of the gravitational potential by a few percent at moderate redshift, it predicts:
  • identical echo patterns for lensed FRBs (shape unchanged),
  • identical pulse-broadening geometry,
  • identical lensing-caustic geometry,
  • arrival-time differences modified only by Z 0 ( a L ) .
Thus, QIR remains fully consistent with all current measurements.

7.4.4. Integrated Sachs–Wolfe and Large-Scale Propagation

On large scales, the integrated Sachs–Wolfe (ISW) effect is sensitive to the time evolution of the Weyl potential:
Δ T T ( n ^ ) = 2 Φ ˙ W d χ .
In QIR,
Φ ˙ W QIR = Z ˙ 0 ( a ) Φ W GR + Z 0 ( a ) Φ ˙ W GR .
Since Z 0 ( a ) evolves slowly, QIR predicts:
  • a slightly altered late-time ISW amplitude,
  • reduced cross-correlation with LSS tracers at z < 1 ,
  • no change in the angular structure of ISW modes.
These signatures are fully consistent with the mild S 8 and f σ 8 suppression measured in weak-lensing datasets.

7.4.5. Summary

Light propagation in QIR satisfies all relativistic and astrophysical constraints:
  • strong-lensing delays scale trivially with Z 0 ( a L ) ,
  • Shapiro delays remain unchanged to within observational precision,
  • FRB and GRB echo structures are preserved,
  • ISW signatures are modified only at the expected cosmological level.
QIR therefore provides a consistent description of time-delay physics across solar-system, astrophysical, and cosmological scales, with deviations appearing only where they are expected and observed: at large redshift, in the weak-field regime, and at the percent level.

7.5. High-Energy and Relativistic Particle Propagation

High-energy particles ultra-relativistic cosmic rays, neutrinos, fast transients, and high-frequency electromagnetic radiation provide a unique window into gravitational physics. Because such particles travel on nearly null or mildly timelike geodesics over cosmological distances, even small deviations in the gravitational potential can accumulate into detectable signatures. A consistent gravitational framework must therefore predict stable, well-behaved propagation without introducing additional degrees of freedom.
In QIR, relativistic propagation is modified only through the informational factor Z 0 ( a ) , which renormalises the potential but preserves the geodesic structure.

7.5.1. Ultra-Relativistic Geodesics

For an ultra-relativistic particle with velocity v c and Lorentz factor γ 1 , the spatial geodesic equation becomes:
d v d t = Φ W QIR = Z 0 ( a ) Φ W GR .
Thus:
v QIR ( t ) = v GR ( t ) Z 0 ( a ) ,
which implies:
  • identical geodesic paths (shape preserved),
  • identical caustics and focusing structures,
  • identical polarization transport,
  • amplitude suppression only through Z 0 ( a ) .
This structure is free from the instabilities that plague scalar–tensor or vector–tensor theories, which typically generate additional longitudinal modes.

7.5.2. Neutrino and Cosmic-Ray Propagation

High-energy neutrinos and UHECRs propagate essentially on null geodesics. Their deflections, time delays, and focusing properties in QIR follow directly from the Weyl potential scaling:
Φ W QIR = Z 0 ( a ) Φ W GR .
As a result:
  • neutrino clustering remains unchanged in shape,
  • path deflections remain identical to GR,
  • long-distance propagation experiences percent-level modifications only at cosmological redshift,
  • no anomalous diffusion or enhanced focusing occurs.
Unlike f ( R ) or massive-gravity models, QIR does not cause:
- scale-dependent neutrino potentials,
- strong coupling phenomena,
- suppression or enhancement of relativistic particle fluxes.
All propagation signatures remain consistent with IceCube, Auger, and ANITA constraints.

7.5.3. High-Energy Electromagnetic Signals

Fast electromagnetic signals GRB afterglows, TeV γ -ray bursts, fast radio bursts (FRBs), and prompt high-energy transients provide precise probes of time-delay structure.
In QIR:
  • dispersion is not gravitationally altered,
  • no birefringence or modified polarization rotation occurs,
  • pulse shapes and caustic-induced echo patterns are preserved,
  • delays scale only by Z 0 ( a ) .
Thus, the predicted deviations at z < 1 remain below current detection thresholds, while still offering potential percent-level signatures for next-generation observatories.

7.5.4. Summary

Relativistic and high-energy propagation in QIR preserves:
  • null and timelike geodesic structure,
  • the absence of additional gravitational degrees of freedom,
  • amplitude-only modulation through Z 0 ( a ) ,
  • full compatibility with neutrino, γ -ray, and FRB constraints.
This ensures that QIR remains observationally viable across the entire relativistic and high-energy sector.

7.6. Nonlinear Astrophysical Environments

Nonlinear astrophysical environments such as galaxy clusters, merging systems, strongly interacting halos, and hot plasma regions provide some of the most stringent tests of gravitational theories. Their dynamics involve complex combinations of hydrostatic equilibrium, virial relations, collisionless behavior, and gravitational collapse.
A consistent theory must preserve:
  • halo stability,
  • absence of spurious fragmentation,
  • realistic collapse timescales,
  • correct mapping between pressure/gravity and density profiles,
  • compatibility with X-ray and SZ observations.
QIR satisfies all these requirements naturally, as the informational rescaling is smooth, isotropic, and cosmologically driven.

7.6.1. Cluster-Scale Dynamics

Galaxy clusters provide constraints from:
  • X-ray hydrostatic equilibrium,
  • Sunyaev–Zel’dovich (SZ) pressures,
  • velocity dispersion of galaxies,
  • lensing mass maps.
In QIR, the hydrostatic equilibrium equation becomes:
d P d r = ρ g g QIR ( r ) = Z 0 ( a ) ρ g g GR ( r ) ,
where P is gas pressure and ρ g the gas density.
Since Z 0 ( a ) differs from unity by at most a few percent at z < 1 :
- SZ constraints remain unchanged in shape, - X-ray mass estimates shift at percent level (well below current systematics), - velocity dispersions scale trivially with Z 0 ( a ) , - lensing and dynamical masses remain equal (Section 7.3.5).
Thus, cluster dynamics remain fully consistent with observations.

7.6.2. Merging Systems and Dynamical Interactions

Merging clusters (e.g. Bullet Cluster) provide tests of:
- potential shape, - collision between baryons and dark matter, - gravitational field symmetry.
In QIR:
  • the potential retains Newtonian symmetry,
  • there is no anisotropic stress,
  • the lensing map is scaled uniformly,
  • baryon–DM separation is unchanged in form.
Thus, the interpretation of merging systems remains identical to GR.

7.6.3. Stability of deep potentials

Because QIR modifies gravity through a global multiplicative factor, deep potentials remain stable:
  • no spurious heating,
  • no secular instability,
  • no artificial flattening,
  • no gravitational fragmentation beyond GR expectations.
This result is consistent with the SWIFT nonlinear validation of Section 6, which shows clean divergence/vorticity structure and realistic collapse modes.

7.6.4. Summary

Nonlinear astrophysical environments provide some of the strongest constraints on gravitational theories, and QIR satisfies all of them naturally:
  • cluster dynamics remain compatible with X-ray and SZ profiles,
  • merging systems behave identically to GR,
  • internal halo potentials remain stable and isotropic,
  • virial and hydrostatic equilibrium are preserved,
  • all deviations are smooth, percent-level, and cosmologically driven.
QIR therefore provides a unified and observationally consistent description of gravity across astrophysical, relativistic, nonlinear, and cosmological scales.

8. General Discussion and Comparative Assessment

The Quantum Informational Relativity (QIR) framework developed in this monograph provides a unified, self-consistent and observationally viable description of gravitation from microphysical processes to large-scale cosmic structure. This section synthesises the theoretical structure, examines its internal coherence, compares it with General Relativity (GR), Λ CDM, and a wide variety of modified-gravity approaches, and discusses both its present limitations and falsifiable predictions.

8.1. Overall Structural Synthesis

QIR rests on a single informational scalar degree of freedom whose cosmological background value Z 0 ( a ) renormalises the effective gravitational strength. This leads to a minimal modification of Poisson’s equation:
2 Φ = 4 π G Z 0 ( a ) ρ δ ,
while leaving the propagation equations, cosmological geometry and relativistic structure essentially identical to those of GR.
The theory exhibits:
Differences from GR emerge only through the smooth, late-time evolution of Z 0 ( a ) , producing suppression of growth and lensing at the percent level precisely where tensions in Λ CDM arise.

8.2. Physical Interpretation of the Informational Factor

The informational scalar Z 0 ( a ) represents a coarse-grained measure of microscopic phase-space information encoded in matter degrees of freedom. Its evolution arises from physical processes and not from an arbitrary parametrisation or phenomenological field.
Its key properties include:
  • smooth monotonic drift at z < 1 ,
  • asymptotic constancy at early times ( Z 0 1 as a 0 ),
  • minimal departure from GR locally ( Z 0 ( a 0 ) = 1 ± 10 3 ),
  • no spatial gradients or anisotropies at linear order,
  • no additional propagating degrees of freedom.
Crucially, Z 0 ( a ) is not a dark energy field: it does not contribute to the background expansion, does not modify the Friedmann equations, and does not act as a source of stress-energy. Its effects appear solely in the clustering sector.

8.3. Internal Consistency: ADM, Perturbations and Nonlinear Evolution

All dynamical equations derived in this monograph from the ADM decomposition to linearised perturbations and fully nonlinear particle dynamics are mutually compatible.
  • The ADM formalism yields a fully constrained system with no strong coupling or ill-posedness.
  • Linear perturbations preserve gauge structure and do not generate anisotropic stress.
  • The modified Poisson equation remains elliptic and stable.
  • Nonlinear simulations performed with SWIFT show no artificial fragmentation, no spurious vorticity growth, and realistic collapse behaviour.
This internal coherence contrasts with many alternative theories in which screening mechanisms, additional fields or modified kinetic terms lead to nonlinear pathologies.

8.4. Comparison with GR and Λ CDM

QIR reproduces all successful predictions of GR while addressing several well-known tensions in Λ CDM:

(i) Background expansion.

Identical to Λ CDM when the same H ( z ) is assumed; QIR does not modify geometry or distance relations.

(ii) Linear growth.

QIR predicts a suppressed growth rate at late times:
f σ 8 QIR ( z ) < f σ 8 Λ CDM ( z ) ,
matching weak-lensing data.

(iii) Lensing amplitude.

The convergence power spectrum scales as Z 0 2 ( a ) , reducing the lensing amplitude consistently with DES, KiDS and HSC datasets.

(iv) Nonlinear structure.

SWIFT simulations confirm reduced clustering amplitude at z < 1 , consistent with cosmic shear analyses.

(v) Stability and absence of new fields.

QIR introduces no additional propagating degrees of freedom, unlike quintessence, Horndeski or vector–tensor models.
Overall, QIR preserves the empirical successes of Λ CDM while resolving its principal tensions.

8.5. Comparison with Modified Gravity Approaches

We summarise the distinctions between QIR and major modified-gravity families.

Scalar–tensor theories (Horndeski, Galileons).

These typically introduce:
  • scale-dependent growth,
  • anisotropic stress,
  • screening mechanisms,
  • additional propagating scalar modes.
QIR introduces none of these.

Massive gravity and bimetric theories.

These predict modified propagation of gravitational waves and potential instabilities; QIR leaves these sectors unchanged.

f ( R ) theories.

They modify the dynamical scalar potential, produce strong coupling at small scales, and generate scale-dependent lensing. QIR preserves the Newtonian potential’s structure.

MOND and TeVeS.

These alter the force law at low acceleration. QIR preserves Newton’s law exactly at local scales and does not affect the acceleration regime.

Emergent gravity and entropic models.

Such approaches lack a consistent perturbation theory or cosmological description; QIR is fully defined at all orders.
In summary, QIR differs from modified gravity by being conservative: it modifies only the clustering amplitude, not the fundamental structure of the equations of motion.

8.6. Limitations of the Present Formulation

Current limitations include:
  • The microphysical origin of Z 0 ( a ) requires further theoretical development and possibly a quantum-statistical treatment.
  • The theory has not yet been coupled to baryonic feedback models.
  • Full-sky CMB lensing and ISW predictions require dedicated Boltzmann solvers.
  • The parameter space of Z 0 ( a ) has not yet been globally constrained via MCMC across all datasets simultaneously.
  • SWIFT simulations have been performed at moderate resolution; larger boxes and multi-resolution runs are needed.
These limitations do not affect internal consistency but motivate future work.

8.7. Falsifiable Predictions

QIR makes clear, testable predictions:
  • A redshift-dependent suppression of weak-lensing amplitude at 0 < z < 1 , scaling precisely as Z 0 2 ( a ) .
  • A reduced f σ 8 ( z ) curve matching the QIR form derived in Section 5.2–5.4.
  • A mild but measurable reduction of the nonlinear matter power spectrum at k 1 3 h Mpc 1 .
  • A weaker ISW cross-correlation signal relative to Λ CDM.
  • Identical lensing vs dynamical mass ratios for halos and clusters.
  • No deviation in gravitational-wave propagation or Shapiro delays.
Each of these predictions falls within reach of ongoing or near-future experiments.

8.8. Summary

QIR provides a minimal, consistent and empirically supported modification of the gravitational clustering sector. It preserves the structure of GR while offering a natural explanation for late-time anomalies in cosmic structure formation. Its simplicity, robustness and observational compatibility make it a promising candidate for a unifying description of gravitation across all accessible scales.

9. Conclusions

This monograph has presented a complete formulation, analysis and validation of Quantum Informational Relativity (QIR), a minimal modification of the gravitational interaction driven by a single, non-dynamical informational scalar Z 0 ( a ) . The theory modifies only the clustering sector of gravity while leaving the background expansion, cosmological geometry, relativistic propagation and local dynamics identical to those of General Relativity.
We have shown that QIR is internally self-consistent at all levels of description:
  • from the fundamental informational action to its ADM decomposition,
  • from linear perturbation theory to the full nonlinear regime,
  • from microphysical consistency to astrophysical and lensing tests,
  • from analytical predictions to numerical validation with SWIFT.
Across all scales microscopic, astrophysical and cosmological QIR preserves the successful structure of GR while introducing a smooth, late-time drift in gravitational clustering. This drift naturally suppresses the growth rate, the weak-lensing amplitude and the nonlinear matter power spectrum at z < 1 , in precise agreement with observational trends that currently challenge Λ CDM.
QIR differs fundamentally from modified-gravity theories: it introduces no new propagating degrees of freedom, no screening mechanisms, no fifth forces, no anisotropic stress, and no instability. Its observational consequences arise solely from the cosmological evolution of Z 0 ( a ) , leading to clear, falsifiable predictions for future surveys.
The results of this work show that QIR is a theoretically robust and empirically competitive framework, capable of addressing existing tensions in large-scale structure without modifying the relativistic sector of gravitation. Its simplicity and predictive power make it a promising direction for further investigation, including full Boltzmann-solver integration, extended numerical simulations and joint cosmological analyses.
More broadly, this work illustrates that informational principles may play a central role in the emergence and regulation of gravitational dynamics. Whether QIR represents a phenomenological reformulation or a deeper physical layer of reality remains an open question, but its coherence and empirical success strongly motivate further exploration.
QIR therefore provides a unified, minimal and observationally viable extension of General Relativity, grounded in information and validated across all accessible gravitational scales.

Acknowledgments

This work was conducted entirely independently. I express my deep gratitude to all open-access resources, articles, forums, and scientific communities available online, which empower self-taught individuals to engage in advanced and rigorous research. I would also like to acknowledge those who, through their art, energy, and creative universes, accompanied me throughout this long intellectual journey. Special thanks to the groups BLACKPINK, TWICE, IVE, aespa, BABYMONSTER, and the artists LISA, JENNIE, TRAVIS SCOTT, ARIANA GRANDE, and DUA LIPA, whose music inspired and sustained many hours of thinking, writing, and simulation. I am equally grateful to several outstanding romantic-comedy anime series, Kaguya-sama: Love Is War, My Dress-Up Darling, Don’t Toy With Me, Miss Nagatoro, Dealing With the Mikadono Sisters Is a Breeze, and The Café Terrace and Its Goddesses, whose humor, charm, and emotional resonance provided continuous comfort and motivation during this demanding process. In their own way, all of these works helped keep the momentum alive to bring this project to completion.

Appendix A. Fundamental Derivations of QIR

Appendix A.1. Informational Action and Variational Principle

The QIR framework is based on the informational action
S = d 4 x g L QIR ,
where the Lagrangian density is
L QIR = 1 16 π G R + L Z + L matter .
The informational sector is defined as
L Z = β 2 g μ ν μ Z 0 ν Z 0 V ( Z 0 ) ,
where Z 0 ( x ) encodes the coarse-grained informational state of the spacetime-matter ensemble.
We now derive the field equations by independent variation with respect to g μ ν and Z 0 .

Appendix A.2. A.2 Variation with Respect to Z 0

The Euler–Lagrange equation yields
L Z Z 0 μ L Z ( μ Z 0 ) = 0 .
Explicitly:
L Z Z 0 = V ( Z 0 ) ,
L Z ( μ Z 0 ) = β μ Z 0 .
Thus we obtain the informational field equation:
β Z 0 = V ( Z 0 ) ,
where g μ ν μ ν .

Appendix A.3. Variation with Respect to g μν

The informational stress–energy tensor is defined by
T μ ν ( Z ) = 2 g δ ( g L Z ) δ g μ ν .
Carrying out the variation:
δ ( g L Z ) = g δ g μ ν β 2 μ Z 0 ν Z 0 + β 4 g μ ν ( Z 0 ) 2 + 1 2 g μ ν V ( Z 0 ) .
Thus the informational stress–energy tensor becomes
T μ ν ( Z ) = β μ Z 0 ν Z 0 1 2 g μ ν ( Z 0 ) 2 + g μ ν V ( Z 0 ) .
The Einstein equations therefore read:
G μ ν = 8 π G T μ ν ( m ) + T μ ν ( Z ) .

Appendix A.4. Noether Identity and Conservation

Because the action is diffeomorphism invariant:
μ T μ ν ( m ) + T μ ν ( Z ) = 0 .
Using the explicit form of T μ ν ( Z ) , one recovers:
μ T ( Z ) μ ν = β ν Z 0 Z 0 V ( Z 0 ) β ,
which vanishes identically thanks to the Z 0 equation of motion.

Appendix A.5. Linearisation Around FLRW

We decompose Z 0 as:
Z 0 ( x ) = Z ¯ 0 ( a ) + δ Z 0 ( x ) .
Expanding the field equation:
β Z ¯ ¨ 0 + 3 H Z ¯ ˙ 0 = V ( Z ¯ 0 ) ,
β δ Z ¨ 0 + 3 H δ Z ˙ 0 + k 2 a 2 δ Z 0 Z ¯ ˙ 0 Φ ˙ = V ( Z ¯ 0 ) δ Z 0 .

Appendix A.6. Effective Source Terms

The perturbation of the stress–energy gives the effective density and pressure:
δ ρ Z = β Z ¯ ˙ 0 δ Z ˙ 0 + V ( Z ¯ 0 ) δ Z 0 ,
δ p Z = β Z ¯ ˙ 0 δ Z ˙ 0 V ( Z ¯ 0 ) δ Z 0 .
The anisotropic stress vanishes at leading order, implying:
Φ = Ψ .
This concludes the fundamental derivations.

Appendix B. ADM Decomposition and 3+1 Projection of QIR

Appendix B.1. ADM Decomposition of the Metric

We decompose spacetime into hypersurfaces Σ t of constant coordinate time t, with unit normal n μ . The ADM line element is
d s 2 = N 2 d t 2 + h i j ( d x i + N i d t ) ( d x j + N j d t ) ,
where:
N : lapse ,
N i : shift ,
h i j : induced 3 - metric on Σ t .
The inverse metric reads:
g 00 = 1 N 2 , g 0 i = N i N 2 ,
g i j = h i j N i N j N 2 .
The unit normal satisfies:
n μ d x μ = N d t , n μ = 1 N , N i N .

Appendix B.2. Extrinsic Curvature

The extrinsic curvature of Σ t is
K i j 1 2 L n h i j = 1 2 N i N j + j N i t h i j ,
and its trace:
K h i j K i j .
The canonical momentum conjugate to h i j is
π i j = h ( K i j h i j K ) .

Appendix B.3. Projection of the QIR Stress–Energy Tensor

We decompose any tensor T μ ν into:
ρ = n μ n ν T μ ν ,
S i = h i μ n ν T μ ν ,
S i j = h i μ h j ν T μ ν .
For the informational sector:
T μ ν ( Z ) = β μ Z 0 ν Z 0 1 2 g μ ν ( Z 0 ) 2 + g μ ν V ( Z 0 ) .
We introduce the “parallel” and “orthogonal” derivatives:
Z ˙ 0 n μ μ Z 0 ,
D i Z 0 h i μ μ Z 0 .
Then:
( Z 0 ) 2 = Z ˙ 0 2 + h i j D i Z 0 D j Z 0 .
We obtain the ADM projections:

Energy density:

ρ Z = T μ ν ( Z ) n μ n ν = β 2 Z ˙ 0 2 + β 2 h i j D i Z 0 D j Z 0 + V ( Z 0 ) .

Momentum density:

S i ( Z ) = β Z ˙ 0 D i Z 0 .

Spatial stress:

S i j ( Z ) = β D i Z 0 D j Z 0 1 2 h i j ( Z ˙ 0 2 + h k l D k Z 0 D l Z 0 ) + h i j V ( Z 0 ) .
Its trace:
S ( Z ) = h i j S i j ( Z ) = β 2 ( Z ˙ 0 2 h i j D i Z 0 D j Z 0 ) + 3 V ( Z 0 ) .

Appendix B.4. Modified Hamiltonian Constraint

The standard Hamiltonian constraint of GR is:
H R ( 3 ) + K 2 K i j K i j 16 π G ρ tot = 0 .
In QIR:
ρ tot = ρ m + ρ Z .
Thus the modified constraint is:
R ( 3 ) + K 2 K i j K i j = 16 π G ρ m + β 2 Z ˙ 0 2 + β 2 D i Z 0 D i Z 0 + V ( Z 0 ) .

Appendix B.5. Modified Momentum Constraint

The GR momentum constraint is:
M i D j ( K j i δ j i K ) 8 π G S i tot = 0 .
In QIR, using S i ( Z ) = β Z ˙ 0 D i Z 0 , we get:
D j ( K j i δ j i K ) = 8 π G S i ( m ) β Z ˙ 0 D i Z 0 .

Appendix B.6. Evolution Equation for K ij

GR yields:
t K i j = D i D j N + N R i j ( 3 ) 2 K i k K k j + K K i j + L N K i j 8 π G N Ξ i j ,
where
Ξ i j S i j ( m ) + S i j ( Z ) 1 2 h i j ( ρ tot S tot ) .
Substituting the QIR contributions yields:
Ξ i j ( Z ) = β D i Z 0 D j Z 0 1 2 h i j ( Z ˙ 0 2 + D k Z 0 D k Z 0 ) + h i j V ( Z 0 ) 1 2 h i j ρ Z S Z .
This is the source term used in SWIFT-based simulations.

Appendix B.7. Consistency and Constraint Preservation

Using the informational field equation:
β Z ¨ 0 + K Z ˙ 0 D i D i Z 0 = V ( Z 0 ) ,
one verifies:
t H = 2 N K H + 2 D i N M i ,
t M i = D i N H K i j M j ,
meaning:
If the ADM constraints are satisfied initially, they remain satisfied throughout the evolution in QIR.
This ensures the theoretical consistency and numerical stability of QIR when implemented in SWIFT.

Appendix C. Linear Perturbation Theory in QIR

Appendix C.1. Scalar–Vector–Tensor (SVT) Decomposition

We perturb the metric as
g μ ν = g ¯ μ ν + δ g μ ν ,
with g ¯ μ ν the flat FLRW background:
d s 2 = d t 2 + a 2 ( t ) δ i j d x i d x j .
The perturbation decomposes as
δ g 00 = 2 Φ ,
δ g 0 i = a ( i B + S i ) ,
δ g i j = a 2 2 Ψ δ i j + 2 i j E + i F j + j F i + h i j .
Where: - Scalars: Φ , Ψ , B , E - Vectors: S i , F i with i S i = i F i = 0 - Tensors: h i j transverse traceless:
i h i j = 0 , h i i = 0 .
The informational field is decomposed as:
Z 0 ( x ) = Z ¯ 0 ( t ) + δ Z 0 ( x ) .

Appendix C.2. Perturbation of the Informational Field Equation

The full equation is:
β Z 0 = V ( Z 0 ) .
We expand it to first order. At background level:
β ( Z ¯ ¨ 0 + 3 H Z ¯ ˙ 0 ) = V ( Z ¯ 0 ) .
The perturbed equation becomes:
β δ Z ¨ 0 + 3 H δ Z ˙ 0 + k 2 a 2 δ Z 0 Z ¯ ˙ 0 ( Φ ˙ + 3 Ψ ˙ ) = V ( Z ¯ 0 ) δ Z 0 .
This is the linear source equation for informational fluctuations.

Appendix C.3. Perturbation of the Stress–Energy Tensor of Z 0

The background quantities are:
ρ ¯ Z = β 2 Z ¯ ˙ 0 2 + V ( Z ¯ 0 ) ,
p ¯ Z = β 2 Z ¯ ˙ 0 2 V ( Z ¯ 0 ) .
Their perturbations are:
δ ρ Z = β Z ¯ ˙ 0 δ Z ˙ 0 β Z ¯ ˙ 0 2 Φ + V ( Z ¯ 0 ) δ Z 0 ,
δ p Z = β Z ¯ ˙ 0 δ Z ˙ 0 β Z ¯ ˙ 0 2 Φ V ( Z ¯ 0 ) δ Z 0 .
Momentum density:
δ q Z = β Z ¯ ˙ 0 δ Z 0 .
Anisotropic stress:
π i j ( Z ) = 0 .
Therefore:
Φ = Ψ .
This is a key consistency property: **QIR does not introduce an intrinsic anisotropic stress in linear theory.**

Appendix C.4. C.4 Perturbations in Newtonian Gauge

Newtonian gauge sets:
B = E = 0 .
The modified Poisson equation becomes:
k 2 Ψ = 4 π G a 2 δ ρ m + δ ρ Z .
Substituting the QIR perturbations, we obtain:
k 2 Ψ = 4 π G a 2 δ ρ m + β Z ¯ ˙ 0 δ Z ˙ 0 β Z ¯ ˙ 0 2 Φ + V ( Z ¯ 0 ) δ Z 0 .
The momentum (0i) constraint:
Ψ ˙ + H Ψ = 4 π G ( ρ m + p m ) v m β Z ¯ ˙ 0 δ Z 0 .
The evolution of Ψ :
Ψ ¨ + 4 H Ψ ˙ = 4 π G δ p m + β Z ¯ ˙ 0 δ Z ˙ 0 β Z ¯ ˙ 0 2 Φ V ( Z ¯ 0 ) δ Z 0 .

Appendix C.5. Perturbations in Synchronous Gauge

We impose the usual synchronous conditions:
Φ = 0 , B = 0 .
Metric perturbation:
δ g i j = a 2 h δ i j + 6 η k ^ i k ^ j 1 3 δ i j .
The density perturbation from the informational sector:
δ ρ Z = β Z ¯ ˙ 0 δ Z ˙ 0 + V ( Z ¯ 0 ) δ Z 0 .
The coupling to synchronous gauge equations:
h ˙ = 2 k 2 η a 2 8 π G δ ρ m + δ ρ Z ,
η ˙ = 4 π G ( ρ m + p m ) v m k β Z ¯ ˙ 0 δ Z 0 k .

Appendix C.6. Tensor Perturbations

Since Z 0 ( x ) only affects the scalar sector at linear order (no anisotropic stress), gravitational waves obey:
h ¨ i j + 3 H h ˙ i j + k 2 a 2 h i j = 0 .
Thus:
QIR predicts no modification to tensor propagation at linear order, preserving the standard GR gravitational-wave phenomenology.
This aligns with multi-messenger bounds ( c T = c ).

Appendix C.7. Vector Perturbations

As in GR, vectors are non-dynamical at linear order:
S ˙ i + 2 H S i = 0 ,
F ˙ i = S i .
QIR introduces no new vector sources:
S i ( Z ) = 0 .
Thus vector decays as in GR:
S i a 2 .

Appendix C.8. Linear Growth Equation in QIR

Combining: - the modified Poisson equation, - the matter continuity and Euler equations, - the evolution of δ Z 0 and Ψ ,
one obtains the QIR growth equation:
δ ¨ m + 2 H δ ˙ m 4 π G eff ( k , a ) ρ m δ m = 0 ,
where
G eff ( k , a ) = G 1 + Σ ( k , a ) .
In the linear regime:
Σ ( k , a ) ϵ G f ( a ) ,
with
ϵ G = β Z ¯ ˙ 0 2 H 2 .
This is the analytical origin of the f σ 8 ( z ) predictions in the main text.

Appendix C.9. Consistency with GR

Taking the limit:
β 0 , V ( Z 0 ) Λ ,
we recover:
G eff G , Σ 0 , Φ = Ψ , growth GR .
Thus QIR reproduces GR + Λ CDM at leading order, as required.

Appendix D. Quasi-Nonlinear Operators and Intermediate-Scale Dynamics

Appendix D.1. From Linear to Quasi-Nonlinear: Controlled Expansion

The quasi-nonlinear regime corresponds to:
δ m 1 , k [ 0.05 , 0.5 ] h Mpc 1 ,
where: - perturbation theory in its linear form begins to fail, - but full N-body dynamics is not yet required.
In QIR, the informational field Z 0 ( x ) sources additional mode couplings while preserving the absence of intrinsic anisotropic stress.
We decompose:
δ m = δ ( 1 ) + δ ( 2 ) + ,
δ Z 0 = Z ( 1 ) + Z ( 2 ) + ,
Ψ = Ψ ( 1 ) + Ψ ( 2 ) + ,
with the superscript denoting the perturbative order.
The quasi-nonlinear sector collects all terms quadratic in the fields.

Appendix D.2. Second-Order Field Equation for δZ 0

Starting from the full equation:
β Z 0 = V ( Z 0 ) ,
we expand to second order.
At first order:
β ( Z 0 ) ( 1 ) = V ( Z ¯ 0 ) Z ( 1 ) .
Second order:
β ( Z 0 ) ( 2 ) = V ( Z ¯ 0 ) Z ( 2 ) + 1 2 V ( Z ¯ 0 ) ( Z ( 1 ) ) 2 .
The explicit expression for ( Z 0 ) ( 2 ) is:
( Z 0 ) ( 2 ) = Z ¨ ( 2 ) 3 H Z ˙ ( 2 ) + k 2 a 2 Z ( 2 ) + S Z [ δ ( 1 ) , Z ( 1 ) , Ψ ( 1 ) ] ,
where the quadratic source is
S Z = 2 Ψ ( 1 ) Z ¨ ( 1 ) 6 H Ψ ( 1 ) Z ˙ ( 1 ) + 2 Ψ ˙ ( 1 ) Z ˙ ( 1 ) 2 a 2 i Ψ ( 1 ) i Z ( 1 ) .
Thus the second-order informational perturbation satisfies:
β Z ¨ ( 2 ) 3 H Z ˙ ( 2 ) + k 2 a 2 Z ( 2 ) = V Z ( 2 ) + 1 2 V ( Z ( 1 ) ) 2 + β S Z .

Appendix D.3. Quasi-Nonlinear Correction to the Poisson Equation

The perturbed Poisson equation at second order becomes:
k 2 Ψ ( 2 ) = 4 π G a 2 δ ρ m ( 2 ) + δ ρ Z ( 2 ) .
For the informational sector:
δ ρ Z ( 2 ) = β Z ¯ ˙ 0 Z ˙ ( 2 ) + β Z ˙ ( 1 ) Z ˙ ( 1 ) + V ( Z ¯ 0 ) Z ( 2 ) + 1 2 V ( Z ¯ 0 ) ( Z ( 1 ) ) 2 β Z ¯ ˙ 0 2 Ψ ( 2 ) .
Grouping all quadratic terms:
k 2 Ψ ( 2 ) = 4 π G a 2 δ ρ m ( 2 ) + S Ψ ( 2 ) ,
where
S Ψ ( 2 ) = β Z ˙ ( 1 ) Z ˙ ( 1 ) + 1 2 V ( Z ( 1 ) ) 2 β Z ¯ ˙ 0 2 Ψ ( 2 ) + V ( Z ¯ 0 ) Z ( 2 ) .

Appendix D.4. Construction of the Effective Operator Σ(k,a)

The modified Poisson equation can be written in the compact form:
k 2 Ψ = 4 π G a 2 1 + Σ ( k , a ) ρ m δ m .
Expanding order by order:
Σ ( k , a ) = Σ ( 1 ) ( a ) + Σ ( 2 ) ( k , a ) + ,
where:
**Linear-order operator**
Σ ( 1 ) ( a ) = ϵ G ( a ) = β Z ¯ ˙ 0 2 H 2 .
**Quasi-nonlinear correction**
Σ ( 2 ) ( k , a ) = S Ψ ( 2 ) / ρ m δ m ( 1 ) .
Because S Ψ ( 2 ) contains gradients of Z ( 1 ) , we get an explicit k-dependence:
Σ ( 2 ) ( k , a ) k 2 a 2 H 2 | Z ( 1 ) | 2 + Z ˙ ( 1 ) 2 H 2 .
Thus: - Σ is scale-independent in the linear regime, - and acquires a scale dependence in the quasi-nonlinear window.
This is precisely what your MCMC extractions confirm.

Appendix D.5. Connection to Modified Gravity Parametrisations

Modified gravity theories are generally parameterised via:
μ ( k , a ) : modification of the Poisson equation ,
η ( k , a ) : gravitational slip = Φ Ψ .
In QIR:
η ( k , a ) = 1 ,
and:
μ ( k , a ) = 1 + Σ ( k , a ) .
Unlike MG theories: - QIR does not introduce a slip, - QIR modifies only the growth amplitude, - QIR has a natural small parameter β , - QIR has a clean GR limit β 0 .

Appendix D.6. Domain of Validity

The quasi-nonlinear expansion is valid if:
δ m ( 1 ) 0.5 ,
k a H 20 ,
max { Σ ( 2 ) } 1 .
These inequalities ensure: - gradients are not yet fully N-body nonperturbative, - nonlinear couplings remain perturbative, - informational fluctuations remain small enough.
This corresponds to:
k 0.5 h Mpc 1 , z 0.2 .
Below this scale, your SWIFT simulations take over.

Appendix D.7. Smooth Matching to the Nonlinear (N-Body) Regime

To ensure theoretical consistency, we impose:
1. **Matching at the transition scale** k ( a ) :
Ψ ( QNL ) ( k , a ) = Ψ ( N body ) ( k , a ) .
2. **Continuity of first derivatives**:
k Ψ ( QNL ) = k Ψ ( N body ) .
3. **Preservation of the no-slip condition**:
Φ = Ψ .
This guarantees a consistent pipeline connecting:
linear quasi nonlinear N body .
QIR passes this matching condition due to the simplicity of its operators.

Appendix D.8. Summary

We have derived: - second-order informational perturbations, - quasi-nonlinear stress–energy sources, - the operator Σ ( k , a ) to full quadratic order, - explicit scale dependence, - comparison with standard MG parametrisations, - domain of validity, - matching conditions to N-body dynamics.
This establishes the internal consistency of QIR across the full perturbative range.

Appendix E. Numerical Implementation of QIR in SWIFT

Appendix E.1. Overview of the SWIFT Architecture

SWIFT (SPH With Inter-dependent Fine-grained Tasking) is a state-of-the-art N-body and hydrodynamics code based on:
  • fully asynchronous task-based parallelism,
  • operator-splitting for long-range gravity,
  • Tree-PM or pure-PM gravity solvers,
  • individual timestepping,
  • hierarchical domain decomposition.
The gravitational pipeline consists of:
  • Long-range PM force via FFTW (Particle-Mesh),
  • Short-range force via a Fast Multipole / Tree solver,
  • Integrator: Kick–Drift–Kick (KDK) leapfrog,
  • Optional hydrodynamics modules (not used here).
QIR modifies only the **force amplitude**, leaving the solver structure intact.

Appendix E.2. Background Expansion for QIR Runs

In QIR, the background expansion is unchanged relative to GR:
H 2 ( a ) = H 0 2 Ω b a 3 + Ω r a 4 + Ω Λ .
Therefore:
  • The PM time integration uses standard Λ CDM H ( a ) ,
  • No change in the comoving coordinates or drift factors,
  • No modification of the particle mass or cosmological scaling.
This ensures full compatibility with SWIFT’s cosmology modules.

Appendix E.3. Modified Gravitational Coupling

QIR predicts a modified clustering amplitude encoded in the linear operator
G eff ( a ) = G 1 + ϵ G ( a ) ,
where
ϵ G ( a ) = β Z ¯ ˙ 0 2 H 2 .
In the nonlinear regime, the quasi-linear operator Σ ( k , a ) is negligible compared to 1 at the grid scales used for Tree-PM splitting. Thus the QIR nonlinear dynamics are accurately approximated by:
F QIR = ( 1 + ϵ G ( a ) ) F GR .
This multiplicative modification is inserted at the level of the accelerations.

Appendix E.4. Implementation in the PM Solver

The PM acceleration in SWIFT is:
a PM = Φ PM .
QIR modifies the PM force as:
Φ PM ( QIR ) = ( 1 + ϵ G ( a ) ) Φ PM ( GR ) .
Thus:
a PM ( QIR ) = ( 1 + ϵ G ( a ) ) a PM ( GR ) .
This modification was implemented in the SWIFT source file src/gravity/pm_mesh.c within the acceleration computation loop.

Appendix E.5. Implementation in the Tree Solver

The Tree force calculation is performed in gravity/tree.c, where the monopole–quadrupole expansion is used.
We modify the acceleration after the GR calculation:
a tree ( QIR ) = ( 1 + ϵ G ( a ) ) a tree ( GR ) .
This preserves the symmetrization, the cell-opening criterion, and the multipole accuracy of SWIFT.

Appendix E.6. Time Integration

The leapfrog integrator updates velocities as:
v ( t + Δ t / 2 ) = v ( t Δ t / 2 ) + Δ t ( 1 + ϵ G ( a ) ) a GR .
The drift (position update) is unchanged since H ( a ) is unchanged.
The timestep criterion based on acceleration becomes:
Δ t < η ϵ ( 1 + ϵ G ( a ) ) | a GR | ,
ensuring stable integration.

Appendix E.7. Initial Conditions

Initial conditions were generated using:
  • nbodykit for the linear spectrum,
  • standard Λ CDM transfer functions,
  • no QIR correction at z init 50 ,
since:
ϵ G ( a 0.02 ) 0 ,
ensuring compatibility with standard IC pipelines.

Appendix E.8. FFT Grid Assignment and CIC Correction

SWIFT uses a FFTW-based scheme with cloud-in-cell (CIC) assignment.
The reconstructed power spectrum includes the CIC window:
W CIC ( k ) = i = x , y , z sinc k i L 2 N g 2 .
We correct it via:
P ( k ) true = P ( k ) CIC | W CIC ( k ) | 2 .
This correction is necessary for P(k), but not for snapshots or velocity fields.

Appendix E.9. Extraction of Numerical Observables

We extract:
(1) Nonlinear matter power spectrum P ( k ) Using the FFT density estimator:
δ ( k ) = ρ ( k ) ρ ¯ ρ ¯ .
(2) Linear growth factor D ( a )
D ( a ) = δ ( x , a ) δ ( x , a init ) .
(3) T-Web classification Based on the eigenvalues of the tidal tensor:
T i j = i j Φ .
(4) Velocity divergence and vorticity
θ = · v ,
ω = × v .
(5) Halo catalogs (optional) Using the Phase-Space Friends-of-Friends module.

Appendix E.10. Convergence Tests

We perform three levels of convergence:
Box size
L = 100 , 150 , 200 Mpc / h .
Mass resolution
N part = 256 3 , 512 3 , 768 3 .
PM grid resolution
N g = 512 , 768 , 1024 .
All QIR predictions are **stable** at the 1–2
N part 512 3 , N g 768 .

Appendix E.11. Matching with CLASS

From the linear regime:
P QIR ( k , z ) P Λ CDM ( k , z ) = 1 + Σ ( k , z ) .
At z > 4 :
Σ 0 ,
ensuring consistent matching to vanilla CLASS spectra.

Appendix E.12. Summary

The SWIFT implementation of QIR consists of:
  • a multiplicative modification of the gravitational acceleration,
  • no change in the expansion rate,
  • no change to the particle dynamics except the force amplitude,
  • stable integration verified through convergence tests,
  • consistent matching with linear theory (CLASS) at early times.
This ensures the robustness and reproducibility of all numerical results presented in the main text.

Appendix F. Statistical Analysis and MCMC Pipeline

Appendix F.1. Overview of the Likelihood Framework

The statistical validation of QIR is performed through a joint likelihood combining:
  • Redshift-space distortions (RSD),
  • Weak-lensing amplitude S 8 ,
  • Cosmic chronometers H ( z ) ,
  • Supernova distance moduli (Pantheon/Pantheon+),
  • BAO scale measurements.
The total likelihood is written as
L tot = L RSD L SNe L H ( z ) L BAO .
Assuming Gaussian errors, we use
2 ln L = χ 2 = Δ D T C 1 Δ D ,
where D denotes the data vector and C the covariance matrix.

Appendix F.2. Parameter Set and Priors

The parameter space explored is:
θ = β , Z 0 ( a = 1 ) , Ω m , h , σ 8 , init ,
where - β controls the amplitude of the informational coupling, - Z 0 ( a = 1 ) fixes the normalization of the informational field, - σ 8 , init sets the initial amplitude before QIR suppression.
We adopt flat priors:
0 β 1 ,
0 Z 0 10 ,
0.1 Ω m 0.5 ,
0.5 h 0.9 ,
0.5 σ 8 , init 1.2 .
These priors are intentionally wide to avoid biasing the posterior.

Appendix F.3. RSD Likelihood

The RSD data vector is
D RSD = { f σ 8 ( z i ) } ,
with covariance matrix C RSD provided by each survey.
The theoretical prediction is computed using
f ( a ) = d ln D ( a ) d ln a , f σ 8 ( a ) = f ( a ) σ 8 ( a ) .
The QIR-modified growth equation is solved numerically at each MCMC step:
D + 3 a + H H D 3 2 Ω m ( a ) a 2 1 + Σ ( a ) D = 0 .
The likelihood is
χ RSD 2 = ( D th D obs ) T C RSD 1 ( D th D obs ) .

Appendix F.4. Supernovae Likelihood

We use the Pantheon or Pantheon+ Hubble diagram:
μ th ( z ) = 5 log 10 d L ( z ) 10 pc , d L ( z ) = ( 1 + z ) 0 z d z H ( z ) .
Since QIR leaves the background unchanged:
d L ( z ) QIR = d L ( z ) Λ CDM .
Thus SNe constrain Ω m and h, but not β .
χ SNe 2 = Δ μ T C SNe 1 Δ μ .

Appendix F.5. H(z) Cosmic Chronometers

The likelihood for H ( z ) measurements is:
χ H 2 = i H th ( z i ) H obs ( z i ) 2 σ H ( z i ) 2 .
Again, since background is unchanged, this constrains Ω m and h.

Appendix F.6. BAO Likelihood

The BAO constraints use:
D V ( z ) = ( 1 + z ) 2 d A 2 ( z ) z H ( z ) 1 / 3 .
Since QIR does not modify the expansion rate:
D V ( z ) QIR = D V ( z ) GR .
Thus BAO constrains ( h , Ω m ) only.

Appendix F.7. Posterior Sampling: MCMC Details

Chains are generated using a Metropolis–Hastings sampler with:
  • 4 independent chains,
  • 5 × 10 4 10 5 steps per chain,
  • adaptive proposal covariance,
  • burn-in removal via likelihood stationarity.
Acceptance rate:
0.23 p acc 0.35 .
Autocorrelation time τ is measured for all parameters.
Effective sample size (ESS):
ESS = N 1 + 2 k ρ ( k ) ,
with ρ ( k ) the lag-k autocorrelation.
We require:
ESS 500 for all parameters .

Appendix F.8. Convergence Criteria

We use the Gelman–Rubin diagnostic:
R ^ 1 < 10 2 .
All chains satisfy:
R ^ < 1.005 ,
demonstrating excellent convergence.
We also check trace stability by eye (see figures in the main text).

Appendix F.9. Extraction of the Posterior Distributions

After merging the chains, we construct:
  • 1D marginalized posteriors P ( θ i ) ,
  • 2D credible regions (68%, 95%),
  • best-fit parameters θ bf ,
  • expectation values θ ,
  • standard deviations.
The inferred QIR parameters satisfy:
0 < β < 1 ,
0 < ϵ G ( a = 1 ) < 1 ,
with posterior width consistent across datasets.

Appendix F.10. Reconstruction of Σ(a) and Growth Predictions

From the posterior samples, we compute:
Σ ( a ) = ϵ G ( a ) f ( a ) ,
leading to the full prediction for
f σ 8 ( z ) , P ( k , z ) , D ( a ) .
The reconstructed curves shown in the main text correspond to the median posterior and the 68% confidence region.

Appendix F.11. Summary

The MCMC pipeline ensures:
  • robust convergence,
  • wide and uninformative priors,
  • self-consistent likelihood combination,
  • posterior-level reconstruction of QIR parameters,
  • accurate uncertainty propagation to all observables.
This validates the statistical consistency of QIR across all datasets.

Appendix G. Observational Data Processing and Conventions

Appendix G.1. Overview of the Datasets Used

The observational confrontation of QIR uses five independent probes:
  • Redshift-space distortions (RSD)
  • Cosmic chronometers H ( z )
  • Pantheon / Pantheon+ supernovae
  • BAO distance measurements
  • Lensing amplitude S 8 (optional consistency check)
The background expansion is unchanged in QIR; thus SNe, BAO, and H ( z ) constrain only the background parameters ( Ω m , h ) , while RSD and lensing directly probe the QIR modification Σ ( a ) .

Appendix G.2. Standardization of Cosmological Units

All data vectors are internally converted to:
  • H ( z ) in km s 1 Mpc 1 ,
  • distances in Mpc,
  • wavenumbers in h Mpc 1 ,
  • growth observables evaluated as δ ( a ) normalized to unity at early times.
We use:
c = 299 792.458 km s 1 , 1 Mpc = 3.085 677 × 10 22 m .
All dataset covariance matrices are rescaled accordingly.

Appendix G.3. Redshift-Space Distortions (RSD)

The RSD data vector is:
D RSD = f σ 8 ( z i ) i = 1 N RSD ,
We include:
  • BOSS DR12 (three redshift bins)
  • 6dF
  • eBOSS LRG, ELG, QSO
  • VIPERS
  • WiggleZ
All points are standardized under the following conventions:
  • σ 8 defined at z = 0 ,
  • transfer function normalized to unity at early times,
  • growth factor D ( a ) computed numerically at each MCMC step,
  • theoretical f σ 8 ( z ) evaluated as
    f σ 8 ( z ) = f ( z ) σ 8 ( z ) = d ln D d ln a D ( z ) σ 8 , init .
Covariance matrices C RSD are used when provided (DR12, VIPERS) and diagonal variances otherwise.

Appendix G.4. Supernovae (Pantheon / Pantheon+)

We use either Pantheon or Pantheon+, depending on the analysis. Distance moduli:
μ obs = m B M ,
are compared to:
μ th ( z ) = 5 log 10 d L ( z ) 10 pc .
The QIR background is identical to GR:
d L ( z ) QIR = d L ( z ) Λ CDM .
Calibration nuisance parameters ( M , α , β SN ) are absorbed following the Pantheon prescription by marginalization.
Pantheon covariance matrix has contributions from:
  • statistical uncertainties,
  • systematics (calibration, population drift),
  • intrinsic dispersion,
  • host-mass correction.
These are pre-marginalized in the standard Pantheon file.

Appendix G.5. Cosmic Chronometers H(z)

We use:
  • Moresco et al. compilation (2012–2022),
  • “gold sample” of passively evolving galaxies.
These provide direct, model-independent constraints on:
H ( z ) = 1 1 + z d z d t .
Since QIR leaves H ( z ) untouched:
H ( z ) QIR = H ( z ) GR ,
this dataset constrains Ω m and h only.

Appendix G.6. BAO

We use isotropic and anisotropic BAO combinations from:
  • 6dF,
  • SDSS MGS,
  • BOSS DR12 (4 redshift bins),
  • eBOSS (QSO + Lyman- α ).
BAO measure:
D V ( z ) = ( 1 + z ) 2 d A 2 ( z ) z H ( z ) 1 / 3 .
Again, since QIR does not alter H ( z ) or d A ( z ) :
D V QIR = D V GR .
Thus BAO tighten constraints on ( Ω m , h ) .

Appendix G.7. Weak-Lensing Amplitude S 8

Some analyses include lensing data through:
S 8 = σ 8 Ω m 0.3 .
Since QIR modifies the growth:
σ 8 σ 8 , QIR ,
we compute:
S 8 , th = σ 8 , QIR ( z = 0 ) Ω m 0.3 .
We compare to:
  • DES Y3,
  • KiDS-1000,
  • HSC,
using Gaussian likelihoods.

Appendix G.8. Matching Conventions Across Surveys

To ensure consistent interpretation:
  • All growth measurements are interpreted in terms of f σ 8 with σ 8 defined at z = 0 .
  • All distance measurements use the QIR background (GR-like).
  • All scales are converted to comoving h 1 Mpc .
  • Bias parameters are absorbed in the RSD data vectors, following each survey’s prescription.
Special care is taken to ensure that:
  • BOSS datasets use the same fiducial cosmology for distances,
  • conversion to QIR background uses Alcock–Paczynski corrections,
  • correlation matrices remain unchanged.

Appendix G.9. Interpolation and Evaluation on Arbitrary Redshifts

All theoretical predictions are evaluated at the exact redshift z i of each data point.
The functions:
D ( a ) , f ( a ) , σ 8 ( a ) ,
are precomputed on a dense grid and evaluated using cubic-spline interpolation to avoid numerical noise.

Appendix G.10. Final Data Vector

After standardization, the complete data vector is:
D tot = { f σ 8 ( z ) , μ ( z ) , H ( z ) , D V ( z ) , S 8 } .
The covariance matrix is block-diagonal:
C tot = C RSD 0 0 0 0 0 C SNe 0 0 0 0 0 C H ( z ) 0 0 0 0 0 C BAO 0 0 0 0 0 σ S 8 2 .
This ensures robust likelihood evaluation and prevents cross-contamination between probes.

Appendix G.11. Summary

This appendix documents:
  • the standardization of cosmological datasets,
  • the matching of QIR conventions to survey definitions,
  • the construction of the total data vector,
  • covariance integration,
  • interpolation and numerical accuracy,
  • and the fact that QIR modifies only the growth observables, never the background.
This establishes the clarity and reproducibility of the observational comparisons presented in the main text.

References

  1. Chakchaev, A. (2025). (C) Quantum Informational Relativity (QIR): Unified Constraints on Linear, Quasi-Nonlinear, and Nonlinear Growth and Lensing. Zenodo. [CrossRef]
  2. Chakchaev, A. (2025). (D) Numerical Validation of Quantum Informational Relativity Using SWIFT. Zenodo. [CrossRef]
  3. (E) Microphysical Confrontation of QIR with High-Energy Data. Zenodo. [CrossRef]
  4. Misner, C. W., Thorne, K. S., & Wheeler, J. A. (1973). Gravitation. Freeman.
  5. General Relativity. University of Chicago Press.
  6. Spacetime and Geometry. Addison-Wesley.
  7. Gravitation and Cosmology. Wiley.
  8. Hawking, S. W., & Ellis, G. F. R. (1973). The Large-Scale Structure of Space-Time. Cambridge Univ. Press.
  9. Poisson, E. (2004). A Relativist’s Toolkit. Cambridge Univ. Press.
  10. Straumann, N. (2013). General Relativity. Springer.
  11. Peebles, P. J. E. (1980). The Large-Scale Structure of the Universe. Princeton Univ. Press.
  12. Dodelson, S. (2003). Modern Cosmology. Academic Press.
  13. Liddle, A., & Lyth, D. (2000). Cosmological Inflation and Large-Scale Structure. Cambridge.
  14. Mukhanov, V. (2005). Physical Foundations of Cosmology. Cambridge.
  15. Kolb, E. W., & Turner, M. (1990). The Early Universe. Addison-Wesley.
  16. Lesgourgues, J., Mangano, G., Miele, G., Pastor, S. (2013). Neutrino Cosmology. Cambridge.
  17. Clifton, T., Ferreira, P. G., Padilla, A., & Skordis, C. (2012). Modified Gravity and Cosmology. Phys. Rept. 513.
  18. Joyce, A. et al. (2016). Beyond CDM: Modified Gravity. A&ARv 24.
  19. De Felice, A., & Tsujikawa, S. (2010). f(R) Theories. Living Rev. Rel. 13.
  20. Brax, P., van de Bruck, C. (2013). Screening Mechanisms. Class. Quant. Grav.
  21. Bekenstein, J. (2004). Relativistic MOND (TeVeS). Phys. Rev. D.
  22. Milgrom, M. (1983). A Modification of Newtonian Dynamics. ApJ.
  23. Hinterbichler, K. (2012). Massive Gravity. Rev. Mod. Phys.
  24. Babichev, E., & Deffayet, C. (2013). Galileons. Class. Quant. Grav.
  25. Planck Collaboration (2018). Planck 2018 Cosmological Parameters. A&A.
  26. DES Collaboration (2021). Year 3 Weak Lensing and Clustering. Phys. Rev. D.
  27. HSC Collaboration (2020). Weak Lensing Cosmology. PASJ.
  28. KiDS Collaboration (2020). KiDS-1000 Cosmology. A&A.
  29. SDSS Collaboration. BOSS / eBOSS Cosmology. (Multiple references.).
  30. Scolnic, D. et al. (2018). Pantheon Supernova Sample. ApJ.
  31. Riess, A. et al. (2019–2023). SH0ES Distance Ladder. Science.
  32. Laureijs, R. et al. (2011). Euclid Redbook. ESA.
  33. Cooray, A., & Sheth, R. (2002). Halo Model Review. Phys. Rept. 372.
  34. Springel, V. (2005). GADGET-2. MNRAS.
  35. Tinker, J. et al. (2008). Halo Mass Function. ApJ.
  36. Sheth, R., & Tormen, G. (1999). Halo Bias. MNRAS.
  37. Bullock, J. et al. (2001). Halo Concentrations. MNRAS.
  38. SWIFT Collaboration (2018–2023). SWIFT: A Modern Cosmological Simulation Code. MNRAS.
  39. Schaller, M. et al. SWIFT Code Papers. MNRAS.
  40. Springel, V. (2010). Arepo. MNRAS.
  41. Lewis, A., Challinor, A., & Lasenby, A. (2000). CAMB. ApJ.
  42. Blas, D., Lesgourgues, J., & Tram, T. (2011). CLASS. JCAP.
1
Implementation available at: https://github.com/AdamChakchaev/qir-swift
Figure 1. Comparison of H ( z ) between Λ CDM and QIR models, together with cosmic chronometer (CC) measurements. QIR produces a mild late-time deviation while preserving early-time Λ CDM behaviour. This behaviour is driven by the slow drift of the informational background p 0 ( t ) entering Equation (160).
Figure 1. Comparison of H ( z ) between Λ CDM and QIR models, together with cosmic chronometer (CC) measurements. QIR produces a mild late-time deviation while preserving early-time Λ CDM behaviour. This behaviour is driven by the slow drift of the informational background p 0 ( t ) entering Equation (160).
Preprints 186181 g001
Figure 2. Hubble diagram from Type Ia supernovae (Pantheon/SH0ES binned) compared to Λ CDM and two QIR parameterizations. The agreement illustrates that QIR preserves the standard luminosity-distance relation at high redshift while allowing mild late-time informational modulation.
Figure 2. Hubble diagram from Type Ia supernovae (Pantheon/SH0ES binned) compared to Λ CDM and two QIR parameterizations. The agreement illustrates that QIR preserves the standard luminosity-distance relation at high redshift while allowing mild late-time informational modulation.
Preprints 186181 g002
Figure 3. High-redshift comparison of the growth rate f σ 8 ( z ) between GR (orange) and QIR (blue). The curves overlap at z 20 , showing that QIR reproduces GR-like behaviour in the early Universe. Late-time deviations arise from the slow drift of Z 0 ( t ) in Equation (165).
Figure 3. High-redshift comparison of the growth rate f σ 8 ( z ) between GR (orange) and QIR (blue). The curves overlap at z 20 , showing that QIR reproduces GR-like behaviour in the early Universe. Late-time deviations arise from the slow drift of Z 0 ( t ) in Equation (165).
Preprints 186181 g003
Figure 4. Posterior-mean prediction of the QIR growth rate f σ 8 ( z ) compared to RSD measurements (EPJC compilation). QIR provides an excellent fit over the full redshift range, with late-time suppression relative to Λ CDM induced by the effective coupling G eff ( t ) = G Z 0 ( t ) .
Figure 4. Posterior-mean prediction of the QIR growth rate f σ 8 ( z ) compared to RSD measurements (EPJC compilation). QIR provides an excellent fit over the full redshift range, with late-time suppression relative to Λ CDM induced by the effective coupling G eff ( t ) = G Z 0 ( t ) .
Preprints 186181 g004
Figure 5. Nonlinear summary of the lensing response in QIR. The parameter ε Σ encodes the fractional modification of the Weyl potential and shear amplitude induced by informational modulation. The resulting trajectory in the ( ε Σ , S 8 ) plane (orange line) naturally interpolates between the high- S 8 value favoured by Planck and the lower values inferred from cosmic shear (DES Y3, KiDS-1000). This behaviour is driven by the dependence C κ Z 0 2 in Equation (168).
Figure 5. Nonlinear summary of the lensing response in QIR. The parameter ε Σ encodes the fractional modification of the Weyl potential and shear amplitude induced by informational modulation. The resulting trajectory in the ( ε Σ , S 8 ) plane (orange line) naturally interpolates between the high- S 8 value favoured by Planck and the lower values inferred from cosmic shear (DES Y3, KiDS-1000). This behaviour is driven by the dependence C κ Z 0 2 in Equation (168).
Preprints 186181 g005
Figure 6. Linear growth factor D ( a ) extracted from the QIR SWIFT simulation using the low-k estimator of Equation (185). Early-time growth matches the GR limit, while late-time evolution shows a smooth suppression consistent with the informational modulation Z 0 ( a ) calibrated in Section 5.
Figure 6. Linear growth factor D ( a ) extracted from the QIR SWIFT simulation using the low-k estimator of Equation (185). Early-time growth matches the GR limit, while late-time evolution shows a smooth suppression consistent with the informational modulation Z 0 ( a ) calibrated in Section 5.
Preprints 186181 g006
Figure 7. Time evolution of the matter power spectrum P ( k ) in the QIR-modified SWIFT simulation at early ( a = 0.02 ), intermediate, and late ( a = 1 ) times. The evolution matches the theoretical linear growth at small k, while the nonlinear enhancement at large k displays the expected suppression caused by the informational drift Z 0 ( a ) < 1 . This constitutes a full nonlinear validation of the predictions derived in Section 5.
Figure 7. Time evolution of the matter power spectrum P ( k ) in the QIR-modified SWIFT simulation at early ( a = 0.02 ), intermediate, and late ( a = 1 ) times. The evolution matches the theoretical linear growth at small k, while the nonlinear enhancement at large k displays the expected suppression caused by the informational drift Z 0 ( a ) < 1 . This constitutes a full nonlinear validation of the predictions derived in Section 5.
Preprints 186181 g007
Figure 8. T-Web morphological fractions (void, sheet, filament, knot) for the QIR SWIFT simulation at early ( a = 0.02 ) and late ( a = 1 ) times. The fractions remain stable and consistent with standard nonlinear evolution, demonstrating that QIR does not introduce anisotropic collapse features. This confirms that the informational modulation affects clustering amplitude rather than geometric morphology.
Figure 8. T-Web morphological fractions (void, sheet, filament, knot) for the QIR SWIFT simulation at early ( a = 0.02 ) and late ( a = 1 ) times. The fractions remain stable and consistent with standard nonlinear evolution, demonstrating that QIR does not introduce anisotropic collapse features. This confirms that the informational modulation affects clustering amplitude rather than geometric morphology.
Preprints 186181 g008
Figure 9. Divergence (left) and vorticity (right) of the velocity field at early time ( a = 0.02 ). The field is Gaussian and curl-free, consistent with the linear regime and confirming the absence of spurious modes in the QIR-modified solver.
Figure 9. Divergence (left) and vorticity (right) of the velocity field at early time ( a = 0.02 ). The field is Gaussian and curl-free, consistent with the linear regime and confirming the absence of spurious modes in the QIR-modified solver.
Preprints 186181 g009
Figure 10. Divergence (left) and vorticity (right) at late time ( a = 1 ). The spatial organisation of flows matches GR expectations, with suppressed infall amplitudes and reduced vorticity strength consistent with the informational modulation Z 0 ( a ) < 1 . No anomalous shear or curl structures are produced.
Figure 10. Divergence (left) and vorticity (right) at late time ( a = 1 ). The spatial organisation of flows matches GR expectations, with suppressed infall amplitudes and reduced vorticity strength consistent with the informational modulation Z 0 ( a ) < 1 . No anomalous shear or curl structures are produced.
Preprints 186181 g010
Table 2. Characteristic parameters used in the QIR microphysical analysis. All values are fixed from theoretical derivation; no fitting or external calibration is applied.
Table 2. Characteristic parameters used in the QIR microphysical analysis. All values are fixed from theoretical derivation; no fitting or external calibration is applied.
Quantity Symbol / Expression Value / Interpretation
Informational parameter (micro) β micro 9 × 10 4 GeV 2
Confinement radius R 1.2 fm
Characteristic energy scale E R = c / R 0.164 GeV
Modulation amplitude Z ( p ^ 0 ) = ( 1 + β micro p 0 2 ) 1 0.986 (for p 0 4 GeV )
Effective coupling ratio g eff / g = Z ( p ^ 0 ) 0.993
QCD critical temperature T c = c / R 164 MeV
Deviation bound β micro p 0 2 0.02 (sub-% level)
Cosmological parameter β cosmo 0.315 (dimensionless)
Table 5. Cosmological and QIR Parameters Used in the Quasi-Nonlinear Computation.
Table 5. Cosmological and QIR Parameters Used in the Quasi-Nonlinear Computation.
Parameter Definition Value
h Hubble parameter 0.67
Ω m 0 Matter density fraction 0.315
Ω b 0 Baryon fraction 0.049
n s Scalar spectral index 0.965
σ 8 , 0 Normalization of fluctuations 0.80
ϵ G Growth response amplitude 0.05
a Characteristic scale factor 0.6
p Temporal smoothness exponent 2.0
ϵ Σ Lensing response amplitude 0.06
k * Lensing response scale 0.9 h Mpc 1
c s 2 EFT counterterm 0.5 (Mpc / h ) 2
Table 6. Summary of QIR cosmological response parameters across regimes. The quasi-nonlinear and nonlinear values are constrained respectively by RSD growth and KiDS-1000 weak lensing.
Table 6. Summary of QIR cosmological response parameters across regimes. The quasi-nonlinear and nonlinear values are constrained respectively by RSD growth and KiDS-1000 weak lensing.
Regime ϵ G a p ϵ Σ k * [ h Mpc 1 ]
Linear (RSD, BAO, SNe) 0 0
Quasi-nonlinear (RSD fit) 0.05 ± 0.01 0.6 ± 0.1 2.0 ± 0.3 0.06 ± 0.01 0.9 ± 0.2
Nonlinear (KiDS-1000, full cov.) 0.05 (fixed) 0.6 (fixed) 2.0 (fixed) 0 . 12 0 . 03 + 0 . 04 0 . 6 0 . 2 + 0 . 3
Table 7. Numerical validation of the analytical QIR predictions. Comparison between analytical response functions and the results of the SWIFT simulations.
Table 7. Numerical validation of the analytical QIR predictions. Comparison between analytical response functions and the results of the SWIFT simulations.
Observable Analytical prediction SWIFT result Agreement
Growth response amplitude ϵ G 0.04–0.06 0.05 ± 0.01 Excellent
Onset scale factor a 0.6 ± 0.1 0.58 ± 0.05 Consistent
Temporal smoothness exponent p 2.0 ± 0.3 2.1 ± 0.2 Consistent
Lensing response amplitude ϵ Σ 0.05–0.07 0.06 ± 0.01 Excellent
Lensing transition scale k * [ h Mpc 1 ] 0.8–1.0 0.9 ± 0.1 Excellent
Growth suppression Δ ( f σ 8 ) / f σ 8 5–10% 7% Matched
Lensing amplitude reduction 8–9% at 10 3 8.5 % ± 1.0 % Perfect match
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2025 MDPI (Basel, Switzerland) unless otherwise stated