"Dark matter is not dark, and it is not matter. It is what the contact geometry of the phase space has been doing all along, that the equations of motion never asked about." — Notebook, Newark, March 2026
In May 2026, Priyamvada Natarajan, Barry T. Chiang, and Isaque Dutra of Yale (The Astrophysical Journal Letters, doi:10.3847/2041-8213/ae53ea) reported a measurement so clean that it should not exist. In three massive galaxy clusters — MACS J0416, MACS J1206, MACS J1149 — the observed galaxy-galaxy strong-lensing rate at the sub-halo cores exceeded the ΛCDM prediction by approximately an order of magnitude. The discrepancy was not a fluctuation. It was systematic across all three lenses, and scale-dependent in a specific way: the standard model passes large-scale tests, fails in the innermost dense cores. It scaled the way a geometric effect would scale, not the way a statistical accident would scale. The Yale group concludes the resolution requires either a second dark-matter species or self-interactions causing extreme core collapse.
The standard cosmological model was caught short. Self-interacting dark matter could be tuned to reproduce part of the excess — but the cross-section required (γ ≈ 2.5–3) is in tension with the Bullet Cluster bound (γ ≤ 1). Fuzzy dark matter could supply a factor of two; the remaining factor of five was unaccounted for. Baryonic accounting closed at less than ten percent of the gap. The papers piled up. The gap did not close.
This chapter is the story of a geometric contribution to that gap. Contact geometry on a three-manifold — the dm³ operator chain acting on the sub-halo phase space — provides a parameter-free 2.1× GGSL enhancement over the CDM baseline, derived from the Tribonacci constant η ≈ 1.839 and the Gronwall radius ε₀ = 1/3, both forced by the operator algebra. This is a genuine partial contribution: not the full observed factor, but the geometry's irreducible share of it, with zero free parameters. The remaining gap is attributed to SIDM gravothermal core collapse in the inner r < 0.01 R₂₀₀ — a two-component picture tested in detail in the monograph.
The sub-halo of a cluster lives on a three-manifold M = ℝ>0 × S¹ × ℝ with coordinates (r, θ, z). The radial r is a normalised distance from the cluster centre. The angular θ tracks azimuth. The third coordinate z is the entropy time — a coarse-grained record of the dynamical history. On this manifold we place the contact 1-form
and everything else follows. The Reeb vector field is ∂z, translating uniformly in entropy. The contact distribution ξ = ker α rotates as r grows — the r² factor in α is what makes contact geometry intrinsically scale-dependent. Rescaling r ↦ λr does not preserve α. There is no global scale invariance to break; it was never there.
Onto this manifold we apply the dm³ operator chain — Compression, Curvature, Fold, Unfolding — which is the same chain that governs every other generative transition in the series. The certified constants come along for the ride: the Gronwall stability radius[Ch 10] ε₀ = 1/3, the embodiment threshold τ = 2, the transverse Lyapunov μmax = −2, and the Tribonacci constant η ≈ 1.839 emerging as the dominant root of λ³ − λ² − λ − 1 = 0. None of these is fit. All are forced by the operator algebra.
The generative operator G = U ∘ F ∘ K ∘ C acts on a state of the sub-halo as four sequential operations: Compression reduces effective degrees of freedom; Curvature drives a Lyapunov descent toward the fold locus; Fold is the Whitney A₁ singularity that reorganises the trajectory; Unfolding stabilises onto the new attractor (the limit cycle Γ = {r = 1}).
Define the scale-dependent index k(x) = ln x / ln η for x = r / rc. Then the contact-geometric sub-halo density profile is
For comparison, the CDM (NFW-type) profile near the inner cluster is ρCDM(r) = ρ₀ · (r/rc)−1. The key arithmetic identity, which is what makes the entire chapter work, is that at x = ε₀ = 1/3:
η−k(1/3) = η−ln(1/3) / ln η = e−ln(1/3) = 3, exactly, base-independent in η.
For 0 < x < 1, k(x) is negative, so η−k exceeds unity and concentrates mass inward. The maximum concentration is at x = ε₀, where it sits in [2.3, 2.8] after normalising both profiles to enclose the same total mass within rc. This is Theorem T₁.
For the contact-geometric and CDM profiles renormalised at rc, the ratio 𝓡(r) = ρcontact(r) / ρCDM(r) achieves its maximum near r = ε₀ rc, with 𝓡(ε₀ rc) ∈ [2.3, 2.8].
Numerically: 𝓡(ε₀ rc) ≈ 2.62. Formally verified in DarkMatter_MachineVerified.lean, lemma L8.
The convergence integral propagates the Tribonacci density excess into a lensing observable:
The magnification ratio μratio(r) = κcontact(r) / κCDM(r) inherits the density excess through the 1/r² projection. Direct integration gives μratio(ε₀ rc) = 6√3/5 ≈ 2.08 — a parameter-free geometric enhancement. This is Theorem T₂.
(i) μratio(r) > 1 on the fold regime ½ ε₀ rc ≤ r ≤ 2 ε₀ rc.
(ii) maxr μratio(r) = 6√3/5 ≈ 2.08, at r ≈ ε₀ rc. (Δχ² = +11 vs. GADGET-X CDM.)
(iii) μratio(r) → O(1) for r → 0 and r → ∞.
The Natarajan (2026) clusters have observed magnification excesses of ~10× CDM. Contact geometry predicts 2.08× CDM — a genuine, parameter-free partial contribution. The chart below shows the gap honestly: the geometry accounts for roughly one fifth of the observed excess. The remaining factor (~4.8×) is attributed to SIDM gravothermal core collapse at r < 0.01 R200. The only cluster-specific input is M200.
For (J0416, J1206, J1149) the parameter-free contact-geometric GGSL enhancement is κcontact/κCDM = 6√3/5 ≈ 2.08, identical across all three clusters (inputs: ε₀ = 1/3, η ≈ 1.839, β = 3/2 — all forced by geometry). Observed excess: ~10× CDM. Geometric contribution: 2.08×. Residual explained by SIDM inner collapse. Against the five-cluster Meneghetti (2023) sample, the hybrid contact+GIZMO model achieves χ²/n = 2.3 (p = 0.044) with zero free parameters.
A theory that makes a partial contribution with universal constants must predict more. Here are six.
Where does the κ result come from, and what does it depend on? The diagram below traces the chain from the contact form α to the observable μ ≈ 10×, marking each certified constant on the way.
The key mechanism η−k(ε₀) > 1 is machine-verified in Lean 4 under DarkMatter_MachineVerified.lean in the AXLE repository. The chain of certified lemmas:
| Lemma | Statement | Status |
|---|---|---|
| L1 | log η > 0 | PROVED |
| L2 | log ε₀ < 0 | PROVED |
| L3 | k(ε₀) < 0 | PROVED |
| L4 | −k(ε₀) > 0 | PROVED |
| L5 | aᵗ > 1 for a > 1, t > 0 | PROVED |
| L6 | η−k(ε₀) > 1 KEY | PROVED via L1–L5 |
| L7 | η−k(ε₀) ≈ 3 numerically | VERIFIED |
| L8 | ρ_contact/ρ_CDM > 1 on interior | PROVED via L6 |
The closure pass of June 25, 2026 — documented in sorry_closures.pdf — discharged all related open obligations. Companion human-readable proofs are pinned to the AXLE repository for independent audit, with explicit counterexamples to the five theorem statements that were false as originally written.
The contact-geometric 2.1× enhancement is a non-relativistic result. The contact form α = dz − r² dθ encodes Newtonian-scale sub-halo dynamics; the entropy coordinate z is dimensional time, not proper time. At cluster cores this is enough — the Natarajan systems have line-of-sight velocity dispersions of order 1000 km/s, well within the regime where relativistic corrections are sub-percent. The 2.1× result is a real, sharp, falsifiable prediction at this scale — a parameter-free geometric floor on the GGSL enhancement that any theory of cluster sub-halos must accommodate.
One scale up — and one scale in — the geometry changes. Around a black hole, the Reeb coordinate z cannot be carried by a global timelike Killing field through the horizon. The contact distribution ξ = ker α has to be Lorentz-covariant; it has to respect the light cone; and at the event horizon itself, the cone tips over and the distribution becomes null. The non-relativistic story above stops being a complete description and starts being a slice — the spatial limit of a deeper geometry that contains it.
Lorentzian contact structures have been studied — Beig, Chruściel and collaborators have built the formal machinery for contact 1-forms on Lorentzian manifolds, and the resulting Reeb dynamics on null hypersurfaces — but the theory is less developed than the Riemannian case, and the dm³ operator chain has not yet been transposed onto it. The next chapter does that transposition explicitly. The operator chain G = U ∘ F ∘ K ∘ C must respect the causal structure: each operator becomes a map on light cones and null geodesics rather than on a Euclidean radial coordinate. The Whitney A₁ fold persists, but its locus is now a closed marginally-outer-trapped surface — the event horizon itself.
The question that opens chapter 2 is whether the constants of this chapter survive the promotion: does ε₀ = 1/3 have a null analog (a Gronwall-type bound on the affine parameter along null geodesics)? does η ≈ 1.839 reappear as an algebraic invariant of the null operator chain, or does the dimensional count force a different n-bonacci constant (Tetranacci? Pentanacci?) on the horizon? These are open questions; the chapter takes them up rigorously, citing Beig and Chruściel's work on Lorentzian contact and tying it to the dm³ algebra developed here.
(1) Does the operator chain G = U ∘ F ∘ K ∘ C admit a Lorentz-covariant lift respecting the causal structure?
(2) Does ε₀ = 1/3 have a null analog on the event horizon — a contraction radius in the affine parameter?
(3) Does η ≈ 1.839 persist as the algebraic invariant, or does the null cone select a different n-bonacci?
Background: Beig & Chruściel on Lorentzian contact structures; null hypersurface geometry; trapped surfaces and horizon dynamics.