Comparison with Modified‑Gravity Frameworks

A wide range of alternatives to General Relativity (GR) have been proposed to address persistent tensions involving galaxy rotation curves, cosmic acceleration, and structure formation. These include MOND and its relativistic extensions, scalar–tensor theories such as TeVeS, higher‑order geometric theories such as f(R) gravity, and various emergent‑gravity approaches. While these frameworks differ in detail, they share a common feature: they retain a fundamentally geometric ontology, treating gravity as a manifestation of spacetime curvature (possibly supplemented by additional fields or modified actions).

The Monad‑Field Framework adopts a different starting point. It posits a substrate ontology in which the familiar 3+1 spacetime description is interpreted as an effective continuum approximation of a deeper three‑dimensional substrate. In this view, the temporal coordinate functions as an ordering parameter associated with substrate evolution rather than a fundamental dimension. Observable gravitational phenomena arise not from curvature but from tension gradients within this substrate, governed by a constitutive relation linking substrate density to internal tension.

This section contrasts the Monad‑Field Framework with several representative modified‑gravity approaches across ontology, dynamical structure, and phenomenological implications.

Ontology of Gravity

In GR, gravity is encoded in the curvature of a four‑dimensional Lorentzian manifold. MOND and TeVeS retain this geometric structure but modify the acceleration law or introduce additional scalar and vector fields. f(R) gravity generalizes the Einstein–Hilbert action by replacing the Ricci scalar with a nonlinear function f(R), yielding higher‑order geometric field equations. Emergent‑gravity models interpret gravity as an effective entropic or information‑theoretic phenomenon but still recover GR‑like geometric dynamics at macroscopic scales.

By contrast, the Monad‑Field Framework introduces a substrate continuum in the macroscopic limit, characterized by a scalar density and an associated tension field. Gravity is interpreted as the material response of this substrate: accelerations arise from spatial gradients of tension rather than from curvature of a metric. The geometric structures used in the effective description are therefore emergent, not fundamental.

Drivers of Gravitational Dynamics

In GR, the stress–energy tensor sources curvature. MOND introduces an acceleration scale below which the effective gravitational law deviates from Newtonian expectations. TeVeS supplements the metric with additional fields whose dynamics reproduce MOND‑like behavior. f(R) gravity modifies the curvature functional, introducing effective additional degrees of freedom. Emergent‑gravity models attribute gravitational dynamics to entropic gradients.

In the Monad‑Field Framework, gravitational dynamics arise from a constitutive relation between substrate density and tension. The “source” of gravitational behavior is therefore not stress–energy in a geometric field equation but the material properties of the substrate itself. This marks a structural shift from geometric to constitutive dynamics.

Galaxy Rotation Curves

GR requires non‑baryonic dark matter to reproduce flat galaxy rotation curves. MOND is explicitly constructed to match these curves, while TeVeS provides a relativistic embedding. Some emergent‑gravity models claim to recover similar behavior through entropic corrections.

Preliminary analyses within the Monad‑Field Framework indicate that tension‑gradient profiles can reproduce observed galactic accelerations without invoking particle dark matter. Exploratory fits to systems such as NGC 1052‑DF2 suggest that the same constitutive mechanism may account for both typical and anomalously low‑DM systems, though full statistical model comparisons remain future work.

Cosmic Acceleration

In GR, late‑time cosmic acceleration is modeled via a cosmological constant Λ, which is highly fine‑tuned. f(R) gravity can mimic an effective cosmological constant, and emergent‑gravity models sometimes attribute acceleration to entropic effects.

In the Monad‑Field Framework, the global evolution of substrate tension suggests a possible mechanism for effective repulsive behavior at cosmological scales. While this provides a conceptual route toward a unified substrate‑based account of both local and large‑scale gravitational phenomena, detailed cosmological solutions and observational fits remain to be developed.

Structure Formation

Structure formation in GR relies on cold dark matter to seed perturbation growth. MOND‑like theories face challenges in reproducing cosmological structure without additional components. TeVeS improves on this but still encounters tensions with data. f(R) and emergent‑gravity models must be tuned to remain consistent with both background expansion and perturbation growth.

In the Monad‑Field Framework, density‑dependent tension gradients may provide a mechanism for gravitational clustering without cold dark matter. This offers a conceptually unified picture in which the same substrate mediates both local gravitational effects and the growth of inhomogeneities. However, explicit perturbation‑growth calculations and power‑spectrum predictions are not yet available, and this remains an open area of development.

Conceptual Position in the Modified‑Gravity Landscape

Persistent tensions involving dark matter, dark energy, and fine‑tuning have motivated exploration of alternatives to purely geometric gravity frameworks. The Monad‑Field Framework contributes to this landscape by proposing a non‑geometric substrate ontology in which gravitational phenomena emerge from constitutive tension gradients rather than curvature. It is therefore neither a deformation of Einstein’s equations nor a phenomenological modification of the acceleration law, but a distinct theoretical program grounded in continuum‑mechanical analogies and emergent‑spacetime ideas.

Summary Comparison Table

Framework	Ontology	Driver of Dynamics	Galaxy Rotation Curves	Cosmic Acceleration	Structure Formation
GR + Dark Matter	4D spacetime curvature	Stress–energy tensor	Requires dark matter halos	Cosmological constant Λ	Dark matter seeds growth
MOND / TeVeS	Geometric + extra fields	Acceleration scale a₀ or scalar–vector fields	Fits without dark matter	Requires additional components	Significant tensions
f(R) Gravity	Higher‑order geometric	Curvature functional	Requires DM or tuned models	Effective Λ possible	Constrained by perturbations
Emergent Gravity	Entropic / holographic	Information gradients	GR + effective DM	Model‑dependent	Less developed
Monad‑Field	3D substrate, emergent spacetime	Density‑dependent tension gradients	Constitutive tension profiles (no DM)	Substrate‑tension evolution	Unified substrate mechanism (in development)

Monad‑Field Framework — Summary

Monad‑Field Framework — Summary (v3.3)

The Monad‑Field framework proposes a new understanding of gravitational phenomena based on a 3D relational substrate. Below is a concise summary of its ontology, effective equations, and observational phenomenology.

Part I — Ontology

Fundamental Substrate: The universe is composed of a 3D relational substrate with a local tension state S(x). This substrate is the fundamental entity, and its internal tension determines all gravitational phenomena.

Emergent Time: Time is not a dimension but an ordering parameter for sequential substrate updates. Clocks measure local state transitions.

Emergent Geometry: Spacetime, metrics, connections, and curvature are not primitive structures but emergent descriptions of the substrate’s constitutive response. The metric gᵤᵥ is interpreted as a local propagation‑impedance tensor.

Part II — Effective Continuum Equations

The framework describes the hydrodynamic limit of the 3D substrate dynamics using equations similar in form to tensor calculus but fundamentally reinterpreted.

Effective Metric:

gᵤᵥ(S) = (1 + αS) ηᵤᵥ

Here ηᵤᵥ is a coordinate fiducial. The metric encodes how excitations propagate, and S modifies the effective update rate between substrate elements.

Substrate Dynamics:

∇ᵤ∇ᵘS − V′(S) = κ|Ψ|² + A′(S) ηᵘᵛ (∂ᵤΨ*)(∂ᵛΨ)

V(S) is a saturating potential, and A(S) encodes constitutive coupling. This equation governs how the substrate’s tension state evolves.

Excitation Dynamics:

∇ᵤ∇ᵘΨ − (m² + κS) Ψ = 0

Matter excitations propagate in the emergent geometry generated by the substrate state S.

Energy–Momentum Conservation:

∇ᵤ Tᵘᵛ = 0

This follows from translational and rotational invariance of the substrate’s relational network.

Part III — Phenomenology

LIGO Post‑Merger Relaxation (Dynamic Regime)

Gravitational‑wave signals exhibit stretched‑exponential relaxation:

A(t) ∝ exp[ − (t / τ)ᵝ ]

with β ≈ 0.35. This is interpreted as the substrate relaxing after being driven into a high‑tension state.

SPARC Galaxy Rotation Curves (Static Regime)

Fits to 158 galaxies yield a characteristic substrate velocity scale V₀ and a constitutive exponent γ. The median γ ≈ 0.43, sharply peaked between 0.3 and 0.5, indicating a universal nonlinear substrate response.

The Monad‑Field static limit modifies the Newtonian baryonic expectation:

V_obs(R) ≈ V_bar(R) · (R / R₀)^(1 − γ)

NGC 1052‑DF2 — Inactive Substrate Limit

This ultra‑diffuse galaxy shows a velocity dispersion consistent with baryons alone. In the Monad‑Field interpretation, the baryonic surface density is below the substrate activation threshold.

A sigma‑clipping fit yields:

γ_DF2 = 1.00
V₀ = 1.50 km/s

This is the pure Newtonian limit: the substrate remains in its linear, uncompressed regime and contributes negligibly.

Universal Substrate Constant

The consistency of γ (from static phenomena) and β (from dynamic phenomena) across spirals, dwarfs, ultra‑diffuse galaxies, and LIGO observations:

0.3 ≲ γ, β ≲ 0.5

suggests that a single material exponent governs both static and dynamic gravitational behavior.

The absence of a “dark matter signal” in DF2 is explained by insufficient baryonic stress to activate the substrate’s nonlinear tension response — confirming a stress‑activation threshold.

Monad‑Field Framework: A Constitutive Substrate Theory of Gravitational Phenomena (v3.3)

Monad‑Field Framework: A Constitutive Substrate Theory of Gravitational Phenomena

White Paper – v3.3 (Unified Dynamic + Static Continuum with Stress‑Activation Threshold)

Abstract

The Monad‑Field framework proposes that the fundamental entity of the universe is a 3D relational substrate whose internal tension state S(x) determines all gravitational phenomena. Time is not a dimension; it emerges as an ordering parameter indexing sequential substrate updates. Geometry, curvature, and spacetime are not primitive structures but large‑scale, emergent descriptions of the substrate’s constitutive response.

The effective 3+1 continuum equations used in this paper are hydrodynamic limits of the substrate’s behavior. They remain mathematically similar to standard tensor calculus, but their interpretation is fundamentally different: the metric is a propagation‑impedance tensor, curvature is a diagnostic of tension gradients, and gravitational dynamics arise from nonlinear substrate elasticity rather than spacetime geometry.

Two independent observational regimes — LIGO post‑merger relaxation and SPARC galaxy rotation curves — reveal the same constitutive exponent (β ≈ 0.35, γ ≈ 0.43), indicating a universal substrate constant governing both dynamic and static gravitational phenomena. A third regime, the ultra‑diffuse galaxy NGC 1052‑DF2, exhibits pure Newtonian behavior (γ = 1.00, V₀ ≈ 1.5 km/s), providing direct evidence for a stress‑activation threshold below which the substrate remains inactive.

Part I – Ontology

1.1 Fundamental substrate

The universe consists of a 3D relational substrate with local tension state S(x). Continuity of S is an effective approximation; the underlying substrate may be discrete or graph‑based. The equilibrium state S = 0 corresponds to isotropic propagation at speed c. The substrate has a finite tension capacity S_max, preventing singularities: as S → S_max, the medium saturates rather than diverging.

1.2 Emergent time

Time is an ordering parameter for substrate updates. A sequence of states {S₀, S₁, S₂, …} (or a continuous evolution label τ) defines “before” and “after”. Clocks measure local state transitions, not motion along a fundamental time axis. Physical time emerges from propagation latency and relaxation rates of the substrate.

1.3 Emergent geometry

Spacetime (3+1) is a large‑scale effective description of substrate dynamics. Geometry is not a primitive cause; it is a constitutive response of the medium.

Metric g_μν – a functional of S and its gradients, encoding how excitations propagate.
Connection Γ – derived from g_μν, describing how changes are transmitted.
Curvature R – a diagnostic of substrate inhomogeneity, not a fundamental field.

Differential geometry is used as computational bookkeeping for the substrate’s constitutive behavior, not as an ontological commitment.

Part II – Effective Continuum Equations

The following equations describe the hydrodynamic limit of the 3D substrate dynamics. The time coordinate t is a coarse‑grained evolution parameter indexing successive substrate states, not a fundamental dimension. All geometric objects are explicit functionals of S and its derivatives.

2.1 Effective metric

g_μν(S) = (1 + α S) η_μν

Here η_μν is a coordinate fiducial with no physical degrees of freedom; it represents the equilibrium propagation state of the substrate. The metric g_μν is interpreted as a local propagation impedance tensor: as S increases, the effective update rate between neighboring substrate elements changes, perceived macroscopically as gravitational time dilation and curvature.

2.2 Substrate dynamics

∇_μ∇^μS − V′(S) = κ |Ψ|² + A′(S) η^μν(∂_μΨ*)(∂_νΨ)

with a saturating potential, for example V(S) = ½ β S² + ¼ γ S⁴ + …, and A(S) = 1 + α S. The constitutive coupling α(S) can be taken to approach zero as S → S_max, ensuring that tension‑gradient forces vanish at saturation and preventing singular cores.

2.3 Excitation dynamics

∇_μ∇^μΨ − (m² + κ S) Ψ = 0

Matter excitations propagate in the emergent geometry generated by the substrate state S.

2.4 Energy–momentum conservation

∇_μ T^μν = 0

Here T^μν includes both substrate tension energy and excitation energy. Conservation follows from translational and rotational invariance of the substrate’s relational network (Noether structure), not from a fundamental geometric identity.

Part III – Phenomenology

3.1 Dynamic regime – LIGO post‑merger relaxation

Post‑merger gravitational‑wave signals exhibit stretched‑exponential relaxation:

A(t) ∝ exp [ − (t / τ)^β ], β ≈ 0.35

This is characteristic of a nonlinear, memory‑bearing medium. In the Monad‑Field interpretation, β is a dynamic constitutive exponent describing how the substrate relaxes after being driven into a high‑tension state by a merger.

3.2 Static regime – SPARC rotation curves

To probe the static behavior of the substrate, the Monad‑Field static‑limit rotation law is fit to 158 galaxies from the SPARC catalog. For each galaxy, two parameters are extracted:

V₀ – characteristic substrate velocity scale (local refresh rate).
γ – constitutive exponent controlling how substrate tension responds to baryonic loading.

Across the sample:

Median V₀ = 27.3 ± 26.0 km/s
Median γ = 0.434 ± 0.446

The distribution of γ is sharply peaked between 0.3 and 0.5, indicating a universal nonlinear substrate response across diverse galactic environments.

3.3 Static Monad‑Field rotation law

The Newtonian baryonic expectation is:

V_bar(R) = √[ V_gas²(R) + V_disk²(R) + V_bul²(R) ].

The Monad‑Field static limit modifies this via a scale‑dependent constitutive factor:

V_obs(R) ≈ V_bar(R) · (R / R₀)^{1 − γ},

where R₀ is a characteristic baryonic scale. The fitted γ for each galaxy is the exponent that maps the Newtonian deficit into the substrate’s tension‑gradient response.

3.4 NGC 3198 – baryon‑heavy spiral

NGC 3198 is a classic flat‑rotation‑curve spiral. Fitting the Monad‑Field static‑limit model yields:

V₀ = 38.27 ± 1.75 km/s
γ = 0.331 ± 0.014

The Monad‑Field total curve (baryonic + substrate) reproduces the observed flat outer profile without invoking a dark‑matter halo. The exponent γ ≈ 0.33 is statistically consistent with the LIGO relaxation exponent β ≈ 0.35, suggesting that both static and dynamic phenomena are governed by the same underlying constitutive law.

3.5 DDO 154 – substrate‑dominated dwarf

DDO 154 is a gas‑dominated dwarf galaxy where the baryonic contribution to the rotation curve is minimal. The fit yields:

V₀ = 23.10 ± 0.36 km/s
γ = 0.431 ± 0.012

Here the substrate provides the majority of the velocity support. The baryonic Newtonian curve fails to reach the observed velocities, while the Monad‑Field curve matches the data. This galaxy is effectively a direct map of the substrate’s tension profile.

The inner‑region deviation (R < 1 kpc), where the model overpredicts the velocity, indicates that the substrate cannot sustain arbitrarily steep central gradients. This motivates a core‑saturation law:

S_eff(R) = S(R) · (1 − e^{−R / R_c}),

where R_c is a core radius set by the substrate’s minimal resolvable scale. As R → 0, the effective gradient is suppressed, naturally producing a finite‑density core and resolving the classical cusp–core problem of cold dark matter within a single constitutive framework.

3.6 Universal substrate constant

Across extremes:

NGC 3198: γ = 0.331 ± 0.014
DDO 154: γ = 0.431 ± 0.012
SPARC median: γ ≈ 0.43
LIGO post‑merger: β ≈ 0.35

All lie in the corridor 0.3 ≲ γ, β ≲ 0.5. This is the signature of a single material exponent governing both static and dynamic gravitational phenomena. γ and β are interpreted as two manifestations of the same constitutive constant: γ controls the static tension–radius relation in galactic equilibrium, while β controls the temporal relaxation of the substrate after mergers.

3.7 NGC 1052‑DF2 – inactive substrate limit

NGC 1052‑DF2 is a diffuse dwarf galaxy whose globular clusters exhibit an unusually small velocity dispersion, consistent with the gravitational influence of baryons alone. In the Monad‑Field framework, this corresponds to a system whose baryonic surface density is below the substrate activation threshold.

Applying a sigma‑clipping Monad‑Field dispersion fit to the DF2 tracer population yields:

γ_DF2 = 1.00
V₀ = 1.50 km/s

This is the pure Newtonian limit of the constitutive law. A value of γ = 1 indicates that the substrate remains in its linear, uncompressed regime, generating no additional centripetal force. The extremely low V₀ confirms that the substrate contributes negligibly to the internal kinematics.

Taken together:

High‑stress spirals (e.g., NGC 3198) exhibit stiffened substrate response (γ ≈ 0.33).
Low‑stress dwarfs (e.g., DDO 154) exhibit softened response (γ ≈ 0.43).
Ultra‑low‑stress diffuse galaxies (DF2) exhibit no substrate activation (γ = 1.00).

The absence of a “dark matter signal” in DF2 is therefore not evidence of missing mass, but of insufficient baryonic stress to activate the substrate’s nonlinear tension response. This is a direct empirical confirmation of the stress‑activation threshold predicted by the Monad‑Field framework.

Part IV – Open Problems

Conservation law closure – derive ∇_μT^μν = 0 from a microscopic substrate action.
Microscopic update rules – specify discrete/relational dynamics from which the continuum PDEs emerge.
Emergence of Lorentz symmetry – show that the equilibrium substrate reproduces special relativity at large scales.
Constitutive propagation tensor – replace metric‑built d’Alembertian with a tensor C^μν(S, ∂S) to make geometry fully emergent computationally.
Quantization – determine whether S, Ψ, or both require quantization and what quantum phenomenology follows.
Relation to other emergent gravity programs – connect to analogue gravity, non‑metric theories, and other constitutive approaches.

Conclusion

The Monad‑Field framework replaces geometric primitives with a 3D constitutive substrate whose nonlinear tension response generates effective spacetime geometry. Time emerges from substrate change; curvature is a diagnostic of tension gradients. The same constitutive exponent appears in post‑merger gravitational‑wave relaxation and in galaxy rotation curves, across baryon‑rich spirals, baryon‑poor dwarfs, and ultra‑diffuse “no‑dark‑matter” systems.

This suggests that dark‑matter‑like galactic dynamics and post‑merger gravitational‑wave behavior are not separate mysteries but two manifestations of a single underlying material: the Monad‑Field substrate, with a stress‑dependent activation threshold that recovers Newtonian gravity in the ultra‑low‑density limit.

The substrate alone is fundamental. Geometry is its memory. Time is its rhythm.

Monad‑Field Framework — Full Constitutive Model (v3.3)

This document presents the complete constitutive equations of the Monad‑Field substrate ontology using pure Unicode math. All terms are defined, and their ontological status is explicitly stated.

1. Full Constitutive Metric

The effective metric depends on both the substrate tension S and its gradients:

gᵤᵥ = (1 + αS) ηᵤᵥ + B(S) (∂ᵤS)(∂ᵥS)

Meaning of terms:

gᵤᵥ — Effective propagation‑impedance tensor.
S — Scalar substrate tension field (fundamental).
α — Linear constitutive coupling.
B(S) — Nonlinear anisotropic constitutive coefficient.
∂ᵤS — Tension gradient; introduces direction‑dependent propagation.
ηᵤᵥ — Fiducial flat tensor (bookkeeping only).

Origin: This is the most general second‑order constitutive metric consistent with locality, isotropic equilibrium, and gradient‑driven anisotropy.

2. Full Connection (Derived from gᵤᵥ)

The Levi‑Civita connection now contains both S and ∂S terms:

Γˡᵤᵥ = ½ gˡʳ ( ∂ᵤgᵥʳ + ∂ᵥgᵤʳ − ∂ʳgᵤᵥ )

Because gᵤᵥ contains ∂S, the connection contains:

terms ∝ ∂S (from the conformal part)
terms ∝ ∂∂S (from the gradient‑metric part)

Interpretation: curvature is no longer just a diagnostic of S gradients, but also of how those gradients vary spatially.

3. Full Curvature Tensor

Rʳₛᵤᵥ = ∂ᵤΓʳᵥₛ − ∂ᵥΓʳᵤₛ + ΓʳᵤˡΓˡᵥₛ − ΓʳᵥˡΓˡᵤₛ

Ontology: Purely diagnostic. Curvature measures how propagation impedance varies due to tension and its gradients.

4. Full Substrate Field Equation

The substrate PDE now includes anisotropic terms from the metric:

∇ᵤ∇ᵘS − V′(S) = κ|Ψ|² + A′(S) ηᵘᵛ (∂ᵤΨ*)(∂ᵛΨ)

But because gᵤᵥ contains ∂S, the operator ∇ᵤ∇ᵘS expands to:

∇ᵤ∇ᵘS = gᵘᵛ ∂ᵤ∂ᵥS − gᵘᵛ Γˡᵤᵥ ∂ˡS

Thus the substrate equation contains:

second derivatives of S
first derivatives of S
mixed terms from ∂S ∂S

Origin: Derived from the full constitutive metric; ensures finite capacity and nonlinear elasticity.

5. Excitation (Matter) Field Equation

∇ᵤ∇ᵘΨ − (m² + κS) Ψ = 0

Because ∇ᵤ contains ∂S, matter propagation becomes direction‑dependent in regions with strong tension gradients.

6. Energy–Momentum Conservation

∇ᵤ Tᵘᵛ = 0

Origin: Noether symmetry of the substrate’s relational structure. Not geometric.

7. Weak‑Field Static Limit

When gradients are small and time‑dependence is negligible:

∇²S ≈ 4πGρ

This fixes κ and links S to the Newtonian potential.

8. Static Rotation Law (Full Quadrature)

The exact static prediction is:

V_obs²(R) = V_bar²(R) + V₀² (R / R₀)^(2γ)

Meaning:

γ — Static constitutive exponent.
V₀ — Substrate velocity scale.
R₀ — Reference radius (1 kpc).

Origin: Quadrature form derived from the substrate tension contribution to the effective potential.

9. DF2 Dispersion Law (Stress‑Activation Threshold)

σ_int(R) = (V₀ / √3) (R / R₀)^γ

DF2 result: γ = 1.00, V₀ ≈ 1.5 km/s → substrate inactive.

Interpretation: Below a critical baryonic surface density, the substrate remains linear and contributes no additional force.

10. Core‑Saturation Law

S_eff(R) = S(R) (1 − e^(−R / R_c))

R_c — Constitutive length scale; minimal resolution of the substrate.

11. Dynamic Relaxation Law (LIGO Regime)

A(t) ∝ exp[ − (t / τ)^β ]

β ≈ 0.35 — dynamic constitutive exponent.
γ ≈ 0.4 — static constitutive exponent.
Relation: Both arise from the same nonlinear substrate kernel.

import numpy as np import matplotlib.pyplot as plt from scipy.optimize import minimize # Data (same as before) gc_data = [ (39, 7.55, 1818, 7, 7), (59, 4.91, 1799, 15, 16), (71, 2.57, 1805, 8, 6), (73, 6.77, 1814, 3, 3), (77, 0.40, 1804, 6, 6), (85, 2.26, 1801, 6, 5), (91, 1.55, 1802, 10, 10), (92, 1.94, 1789, 7, 6), (98, 3.59, 1764, 14, 11), (101, 4.77, 1800, 14, 13) ] radii = np.array([d[1] for d in gc_data]) v_obs = np.array([d[2] for d in gc_data]) e_min = np.array([d[3] for d in gc_data]) e_plus = np.array([d[4] for d in gc_data]) e_sym = np.maximum(e_min, e_plus) v_sys = np.median(v_obs) v_rel = v_obs - v_sys # Negative log-likelihood function def neg_log_likelihood(params, R, v_rel, e_v): V0, gamma = params Vc = V0 * (R / 1.0)**gamma sigma_int = Vc / np.sqrt(3) total_var = sigma_int**2 + e_v**2 if np.any(total_var <= 0): return 1e15 nll = 0.5 * np.sum(np.log(total_var) + v_rel**2 / total_var) return nll # Fit function that returns gamma and V0 def fit_monad_field(R, v_rel, e_v): x0 = [8.0, 1.0] bounds = [(1.0, 30.0), (0.0, 1.5)] res = minimize(neg_log_likelihood, x0, args=(R, v_rel, e_v), bounds=bounds, method='L-BFGS-B', options={'ftol': 1e-12, 'maxiter': 2000}) return res.x[0], res.x[1], res.fun # Sigma-clipping: remove points where |v_rel| > n_sigma * std(v_rel) n_sigma = 2.0 orig_gamma, orig_V0, orig_nll = fit_monad_field(radii, v_rel, e_sym) print(f"Original fit: V0 = {orig_V0:.2f}, γ = {orig_gamma:.3f}, nll = {orig_nll:.2f}") # Iterative clipping mask = np.ones(len(radii), dtype=bool) for iteration in range(5): current_radii = radii[mask] current_v_rel = v_rel[mask] current_e = e_sym[mask] if len(current_radii) < 5: break gamma, V0, nll = fit_monad_field(current_radii, current_v_rel, current_e) # Compute residuals (absolute velocities) Vc_fit = V0 * (current_radii / 1.0)**gamma sigma_fit = Vc_fit / np.sqrt(3) resid = np.abs(current_v_rel) - sigma_fit std_resid = np.std(resid) # Find outliers outlier_mask = np.abs(current_v_rel) > (n_sigma * np.std(current_v_rel)) if not np.any(outlier_mask): break # Update mask for original data indices new_mask = mask.copy() new_mask[new_mask] = ~outlier_mask mask = new_mask print(f"Iter {iteration+1}: removed {sum(outlier_mask)} points, γ = {gamma:.3f}") # Final fit on clipped data R_clip = radii[mask] v_rel_clip = v_rel[mask] e_clip = e_sym[mask] gamma_final, V0_final, nll_final = fit_monad_field(R_clip, v_rel_clip, e_clip) print("\n=== Sigma-clipping result ===") print(f"Kept {len(R_clip)} of {len(radii)} points") print(f"Final γ = {gamma_final:.3f}") print(f"Final V0 = {V0_final:.2f} km/s") # Plot comparison plt.figure(figsize=(9,5)) plt.errorbar(radii, np.abs(v_rel), yerr=e_sym, fmt='o', color='black', alpha=0.5, capsize=3, label='All GCs (original)') if len(R_clip) < len(radii): plt.errorbar(R_clip, np.abs(v_rel_clip), yerr=e_clip, fmt='o', color='blue', capsize=3, label='Kept after clipping') R_model = np.linspace(0.4, 8, 200) sigma_model = V0_final * (R_model / 1.0)**gamma_final / np.sqrt(3) plt.plot(R_model, sigma_model, 'r-', lw=2, label=f'Monad‑Field dispersion (γ = {gamma_final:.2f})') plt.xlabel('Radius (kpc)') plt.ylabel('|v_rel| (km/s)') plt.title('NGC 1052-DF2: Sigma‑clipping test (γ ≥ 0)') plt.legend() plt.grid(alpha=0.3) plt.tight_layout() plt.show()

Monad‑Field Anisotropic Lensing Signature

In the full Monad‑Field constitutive model, the effective metric contains a gradient term g_μν ⊃ B(S) (∂_μS)(∂_νS). This makes light propagation anisotropic in regions where the substrate tension S has strong spatial gradients. As a result, the weak‑lensing shear pattern is no longer purely tangential, but acquires an orientation‑dependent component aligned with ∇S.

Numerical Demonstration

We model an elliptical baryonic mass distribution (cluster‑like) on a 2D grid, solve the Poisson equation ∇²S = ρ for the substrate field S, and compute its gradients ∂_xS and ∂_yS. From these we construct an anisotropic lensing potential

ψaniso ∝ (∂xS)² − (∂yS)² .

The shear components are then obtained from ψ_aniso via second derivatives:

γ₁ = (∂_x² − ∂_y²) ψ / 2
γ₂ = ∂_x∂_y ψ

The resulting maps show a clear quadrupolar shear pattern whose orientation follows the major axis of the baryonic mass distribution, i.e. the direction of the dominant substrate tension gradients.

Radialization and Anisotropy

In standard General Relativity, weak‑lensing shear is purely tangential around the mass. In the Monad‑Field model, the anisotropic term B(S) (∂S)(∂S) causes the shear along the major axis to tilt toward the radial direction (“radialization”), breaking the circular symmetry of the tangential pattern. This can be quantified by an anisotropy ratio

A = γmajor / γminor ,

where γ_major and γ_minor are the shear magnitudes at fixed radius along the major and minor axes. GR predicts A ≈ 1 (after circular averaging), while the Monad‑Field anisotropic metric generically predicts A > 1.

Stress‑Activation and Shear Nulls

The anisotropic signal is strongest where |∇S| is large and vanishes in low‑gradient regions, reflecting the same stress‑activation behavior seen in rotation curves and DF2. In the total shear map (Newtonian + anisotropic), localized “shear nulls” appear where the anisotropic contribution nearly cancels the tangential shear. These nulls, aligned with the baryonic major axis, are a distinctive prediction of the Monad‑Field gradient‑metric term.

Observational Test

A direct test is to stack weak‑lensing maps of galaxy clusters with known baryonic orientations and measure the anisotropy ratio A. A null result (A ≈ 1) would place strong constraints on B(S); a positive detection (A > 1), with shear aligned to the baryonic major axis and possible shear nulls, would be a unique signature of the Monad‑Field anisotropic metric g_μν ⊃ B(S) (∂_μS)(∂_νS).

# ============================================================ # Monad-Field Anisotropic Lensing - Final Version # Improvements: # 1. Direct shear from ∂ᵢS ∂ⱼS (no extra potential) # 2. Quiver plot of shear orientation # 3. Comparison with Newtonian (tangential) shear + residual # 4. Stress‑activation threshold (B(S) as sigmoid of |∇S|²) # ============================================================ import numpy as np import matplotlib.pyplot as plt from scipy.ndimage import gaussian_filter from scipy.fft import fft2, ifft2, fftfreq # ------------------------------------------------------------ # 1. Setup: grid and elliptical baryonic mass distribution # ------------------------------------------------------------ N = 256 extent = 5.0 # half‑size in Mpc x = np.linspace(-extent, extent, N) y = np.linspace(-extent, extent, N) X, Y = np.meshgrid(x, y) dx = x[1] - x[0] # Elliptical cluster (axis ratio 0.6, orientation 30 deg) axis_ratio = 0.6 theta = 30.0 * np.pi / 180.0 cos_t, sin_t = np.cos(theta), np.sin(theta) X_rot = X * cos_t + Y * sin_t Y_rot = -X * sin_t + Y * cos_t R_ell = np.sqrt((X_rot/1.5)**2 + (Y_rot/0.9)**2) rho_baryon = 1.0 / (1.0 + R_ell**2)**1.5 # smooth Plummer # ------------------------------------------------------------ # 2. Compute substrate field S from Poisson equation ∇²S = 4πGρ # (using FFT; G=1 for simplicity) # ------------------------------------------------------------ rho_fft = fft2(rho_baryon) kx = 2 * np.pi * fftfreq(N, d=dx) ky = 2 * np.pi * fftfreq(N, d=dx) k2 = kx[:, None]**2 + ky[None, :]**2 k2[0,0] = 1.0 S_fft = -rho_fft / k2 S = np.real(ifft2(S_fft)) S -= S.min() S = gaussian_filter(S, sigma=1.0) # ------------------------------------------------------------ # 3. Gradients of S # ------------------------------------------------------------ dSdx = np.gradient(S, dx, axis=1) dSdy = np.gradient(S, dx, axis=0) grad_mag_sq = dSdx**2 + dSdy**2 # ------------------------------------------------------------ # 4. Activation threshold: B(S) as sigmoid of |∇S|² # (substrate stiffens only above a critical gradient) # ------------------------------------------------------------ # Parameters from DF2: low gradient -> B ≈ 0; high gradient -> B → B0 threshold = 0.5 # (gradient in (Mpc⁻¹)² – tune for visual effect) B0 = 0.3 # maximum coupling strength B_S = B0 * np.tanh(grad_mag_sq / threshold) # ------------------------------------------------------------ # 5. Anisotropic shear directly from ∂ᵢS ∂ⱼS (no extra potential) # In weak lensing, the shear field γ = (γ1, γ2) is proportional to # the anisotropic part of the metric perturbation. # Here we set γ₁ ∝ (∂ₓS)² - (∂ᵧS)² , γ₂ ∝ 2 ∂ₓS ∂ᵧS # ------------------------------------------------------------ gamma1_aniso = B_S * (dSdx**2 - dSdy**2) gamma2_aniso = B_S * (2 * dSdx * dSdy) # ------------------------------------------------------------ # 6. Newtonian (GR) shear from the baryonic mass (tangential) # Using the standard lensing potential φ_N from ρ_baryon # ------------------------------------------------------------ phi_fft = -rho_fft / k2 # φ_N = -∇⁻²ρ phi_N = np.real(ifft2(phi_fft)) phi_N = gaussian_filter(phi_N, sigma=1.0) # Second derivatives for convergence κ and shear phi_xx = np.real(ifft2(-kx[:,None]**2 * phi_fft)) phi_yy = np.real(ifft2(-ky[None,:]**2 * phi_fft)) phi_xy = np.real(ifft2(-kx[:,None]*ky[None,:] * phi_fft)) kappa_N = 0.5 * (phi_xx + phi_yy) # convergence (not used directly) gamma1_N = 0.5 * (phi_xx - phi_yy) gamma2_N = phi_xy # ------------------------------------------------------------ # 7. Total shear = Newtonian + Anisotropic (Monad-Field) # ------------------------------------------------------------ gamma1_tot = gamma1_N + gamma1_aniso gamma2_tot = gamma2_N + gamma2_aniso # ------------------------------------------------------------ # 8. Plotting # ------------------------------------------------------------ fig, axes = plt.subplots(2, 3, figsize=(15, 10)) # (a) Baryonic mass im = axes[0,0].imshow(rho_baryon, extent=[-extent,extent,-extent,extent], origin='lower') axes[0,0].set_title('Baryonic mass (elliptical)') plt.colorbar(im, ax=axes[0,0]) # (b) Substrate tension S im = axes[0,1].imshow(S, extent=[-extent,extent,-extent,extent], origin='lower') axes[0,1].set_title('Substrate field S') plt.colorbar(im, ax=axes[0,1]) # (c) Activation factor B(S) im = axes[0,2].imshow(B_S, extent=[-extent,extent,-extent,extent], origin='lower') axes[0,2].set_title('Activation B(S) (sigmoid of |∇S|²)') plt.colorbar(im, ax=axes[0,2]) # (d) Anisotropic shear (γ₁, γ₂) as colour, with quiver overlay # Downsample for quiver (skip every 16 pixels) step = 16 skip = (slice(None, None, step), slice(None, None, step)) # Magnitude and angle of the anisotropic shear gamma_aniso_mag = np.sqrt(gamma1_aniso**2 + gamma2_aniso**2) phi_aniso = 0.5 * np.arctan2(gamma2_aniso, gamma1_aniso) # factor 0.5 for tensor # Quiver vectors: (cosθ, sinθ) scaled by magnitude U = gamma_aniso_mag[skip] * np.cos(phi_aniso[skip]) V = gamma_aniso_mag[skip] * np.sin(phi_aniso[skip]) ax = axes[1,0] im = ax.imshow(gamma_aniso_mag, extent=[-extent,extent,-extent,extent], origin='lower', cmap='viridis') ax.quiver(X[skip], Y[skip], U, V, headwidth=0, headlength=0, alpha=0.6, color='red') ax.set_title('Anisotropic shear (Monad-Field only)\nwith orientation sticks') plt.colorbar(im, ax=ax) # (e) Total shear (Newtonian + anisotropic) gamma_tot_mag = np.sqrt(gamma1_tot**2 + gamma2_tot**2) phi_tot = 0.5 * np.arctan2(gamma2_tot, gamma1_tot) U_tot = gamma_tot_mag[skip] * np.cos(phi_tot[skip]) V_tot = gamma_tot_mag[skip] * np.sin(phi_tot[skip]) ax = axes[1,1] im = ax.imshow(gamma_tot_mag, extent=[-extent,extent,-extent,extent], origin='lower', cmap='plasma') ax.quiver(X[skip], Y[skip], U_tot, V_tot, headwidth=0, headlength=0, alpha=0.6, color='white') ax.set_title('Total shear (Newtonian + Monad)') plt.colorbar(im, ax=ax) # (f) Residual: Monad-only part (already shown) – but we can show difference in orientation # Alternative: the difference between total and Newtonian orientation ax = axes[1,2] im = ax.imshow(gamma1_aniso, extent=[-extent,extent,-extent,extent], origin='lower') ax.set_title('γ₁ (anisotropic component)') plt.colorbar(im, ax=ax) plt.tight_layout() plt.savefig('monad_field_lensing_final.png', dpi=150) plt.show() # ------------------------------------------------------------ # 9. Interpretation output # ------------------------------------------------------------ print("\n=== Monad-Field Lensing Signature ===") print("The anisotropic shear pattern (red sticks) aligns with the major") print("axis of the baryonic mass, while the Newtonian tangential shear") print("would be circular (not shown separately). The total shear (white sticks)") print("combines both, creating a distinctive quadrupole asymmetry.") print("\nA null detection in stacked lensing of clusters with known") print("baryonic orientation would constrain the coupling B(S).") print("A positive detection would be a unique signature of the") print("anisotropic term B(S) ∂_μS ∂_νS in the Monad-Field metric.")

Appendix A – Microscopic Motivation for the Emergent Metric

The effective continuum metric

gμν(S) = (1 + αS) ημν

is not an arbitrary ansatz. It emerges naturally from a simple discrete substrate model in which the local propagation delay (update time) between neighbouring relational elements depends linearly on the substrate tension S. This appendix provides a minimal derivation in 1+1 dimensions; the extension to 3+1 dimensions follows by isotropy.

A.1 Discrete 1D Substrate

Consider a one‑dimensional chain of relational nodes with spacing ε (the fundamental lattice scale). Each node i carries a scalar state S_i representing the local substrate tension. Time is not a pre‑existing dimension; it is an emergent ordering parameter labelling sequential state updates. The physical time t is the cumulative number of update steps scaled by a fundamental clock rate τ₀.

A signal propagates from node i to node i+1 in a time

Δti→i+1 = τ₀ (1 + α Si)

This rule encodes that higher substrate tension slows down the transmission of information, similar to a denser optical medium.

A.2 Continuum Limit

Let ε → 0 and N → ∞ with macroscopic length L = Nε fixed. The position coordinate is x = iε. The total time for a signal to travel from x = 0 to x = L is

t(L) = ∫₀ᴸ (τ₀ / ε) (1 + α S(x)) dx

Thus the local propagation speed becomes

v(x) = dx/dt = c₀ / (1 + α S(x))

where c₀ = ε / τ₀ is the speed in the undisturbed substrate (S = 0).

A.3 Emergent Invariant Interval (1+1 Dimensions)

For a light‑like signal, the invariant interval satisfies ds² = 0. In a medium with position‑dependent speed, the line element is

ds² = −c(x)² dt² + dx² = −c₀² (1 + αS)² dt² + dx²

In the weak‑field limit (αS ≪ 1), expanding to linear order gives

ds² ≈ (1 + 2αS) ( −c₀² dt² + dx² )

Defining the effective metric g_μν = (1 + αS) η_μν (absorbing constants into α), we recover the postulated conformal form. The extension to 3+1 dimensions follows from spatial isotropy of the substrate’s equilibrium state.

A.4 Anisotropic Extension

The simple isotropic rule above does not produce the anisotropic term B(S) ∂_μS ∂_νS used in the full constitutive metric. That term arises naturally when the propagation delay depends not only on the local S but also on the gradient of S. For example, a discrete rule of the form

Δti→i+1 = τ₀ ( 1 + α Si + β (Si+1 − Si)² )

leads, in the continuum limit, to an additional contribution proportional to (∂S)² in the effective metric. This shows that the full metric ansatz is a natural extension of the basic tension‑dependent delay model.

A.5 Relationship to the Continuum Field Equations

The discrete update rule does not by itself yield the nonlinear equation for S. However, once the metric is expressed as a functional of S, the usual principles of effective field theory (or a variational principle for the substrate) can be invoked to postulate the wave equation for S.

The derivation demonstrates that the conformal factor (1 + αS) is not an ad‑hoc assumption but a direct continuum manifestation of a simple relational rule: higher substrate tension slows down the update rate between neighbouring elements.