From Identity to Measurement

Worked deductions from the Primorial Unification to astronomical distances

Part I derives the mathematical framework from axioms. Part II tests it against observations.

I. The Identity

Axiom. The field equals its own rate of change: Φ = ∂_tΦ. The four terms of the identity map to the four primorial primes: P₂ → c, P₃ → g, P₅ → Φ, P₇ → ∂_tΦ.

II. The Derivation

We derive the Fibonacci recurrence from the identity in two independent ways, then confirm they agree.

Step 1 — Continuous eigenfunction.   Φ = dΦ/dt has one eigenfunction (up to scale): Φ(t) = Ce^t. The field generates its own growth.

Step 2 — Discrete self-reference.   Now pass to discrete scales. Let R = a_n/a_n−1 be the ratio between adjacent levels. The identity Φ = ∂_tΦ says the field is its own generator — in ratio language, the field at each scale decomposes into two parts:

The parameter-free combination is:

This is not an analogy. The fixed-point equation R = 1 + 1/R is Φ = ∂_tΦ in algebraic form: the field at each scale equals unity plus the field at the reciprocal scale. The self-referential structure of the identity maps directly to the self-referential fixed-point equation of the golden ratio.

Now recover the recurrence. Write R_n = a_n/a_n−1 and multiply through by a_n−1:

Step 3 — Independent confirmation by enumeration.   We can reach the same result by a separate path. The identity Φ = ∂_tΦ is linear (superposition holds), has no free parameters (all coefficients are unity), and has a non-trivial solution (e^t grows). Enumerate all linear recurrences with unit coefficients:

Order 2 is uniquely selected: it is the minimal non-trivial case, and it is the same recurrence derived from the fixed-point equation in Step 2. The two independent derivations converge on the same structure.

Step 4 — Characteristic equation.   Substituting the ansatz a_n = rⁿ:

Step 5 — The Euler bridge.   The continuous eigenvalue e and the discrete eigenvalue φ are not merely analogous — they are linked by Euler’s formula. From the characteristic equation, φ satisfies:

The denominator 5 = P₅ is the primorial prime mapped to Φ in the identity. The path from axiom to eigenvalue runs through Euler at angle π/P₅.

Corollary — π from the lattice.   The Euler bridge reads in both directions. Forward: φ = 2cos(π/5). Backward:

This is not merely an algebraic rearrangement. It reveals that the lattice contains π structurally. Each shell at distance φⁿ is a sphere — its circumference is 2πφⁿ, its surface area 4πφ²ⁿ. The ratio of circumference to diameter is:

The two constants occupy orthogonal roles: φ governs the radial structure (shell spacing), π governs the angular structure (rotational geometry). A sphere rotating through the radial lattice — as Earth does — produces a continuous wave encoding both: φ in the amplitude envelope, π in the phase. The waveform is not random; it is a quasiperiodic signal whose two fundamental frequencies are the two most irrational numbers in mathematics, locked together by the pentagon at angle π/P₅.

Step 6 — Binet’s solution.   The general solution to a_n = a_n−1 + a_n−2 is:

This is the critical bridge: the Fibonacci recurrence does not merely “involve” φ — its solutions are φⁿ (plus an exponentially vanishing correction). The lattice distances d_n = φⁿ are the dominant eigenmode of the recurrence, not an approximation.

The connection is exact. For all integers n:

Step 7 — Self-consistency.   The derivation must close: the recurrence obtained in Step 2 should independently recover the fixed-point equation it was derived from. Define R_n = a_n/a_n−1 and divide the recurrence by a_n−1:

The circle closes. Step 2 derived the recurrence from R = 1 + 1/R. Here, the recurrence independently recovers R = 1 + 1/R. The two paths are algebraic inverses of each other — one is not an assumption but a consequence that can be checked.  ◼

Step 8 — Universality.   For any positive seeds (a₀, a₁), the ratio a_n/a_n−1 → φ. The eigenvalue is independent of initial conditions — just as e^t is the universal eigenfunction of d/dt regardless of amplitude.  ◼

II.b The Selection Principle

Why φ and not some other number.

The derivation above establishes that φ is the eigenvalue of self-referential dynamics. But mathematics produces many eigenvalues. The physical question is: why does nature select this one?

The answer comes from three independent lines — dynamical stability, phase-space coverage, and self-referential closure — each arriving at φ by a different path.

Line 1: KAM Stability (Hamiltonian dynamics)

Consider a bounded physical system with N coupled oscillation modes at frequencies ω₁, …, ω_N. The KAM theorem (Kolmogorov 1954, Arnold 1963, Moser 1962) proves: quasi-periodic orbits survive perturbation if and only if their frequency ratios satisfy the Diophantine condition

The constant K measures how much perturbation the orbit can absorb before it breaks. The orbit with the largest K is the last to be destroyed as perturbation strength increases. It is the most stable orbit in the system.

Which frequency ratio maximizes K? The answer depends on the continued fraction expansion. For a ratio α = [a₀; a₁, a₂, …], the rational approximants p_n/q_n have denominators growing as q_n ~ ∏ a_i. The ratio hardest to approximate — the one whose q_n grow slowest — has all partial quotients equal to 1:

Theorem (Greene 1979, Mackay 1983): The last KAM torus destroyed under increasing perturbation has frequency ratio φ. This is proven for the standard map and conjectured universal for area-preserving maps. The numbers whose continued fractions end in all 1s are called noble numbers; φ is the most noble.

This is not a new claim. It is established dynamical systems theory. What is new is the application: if spacetime is a Hamiltonian system with coupled oscillation modes (it is — the Einstein field equations are a constrained Hamiltonian system), then the most dynamically stable frequency ratios in that system are noble, and the most stable of all is φ.

Line 2: Equidistribution (entropy maximization)

The Weyl equidistribution theorem says the sequence {nα mod 1} for irrational α is uniformly distributed on [0, 1). But the rate of equidistribution differs. The discrepancy D_N — the maximum deviation from uniformity after N terms — is bounded by:

Physical meaning: a system with mode ratio φ explores its phase space most uniformly in finite time. It visits every region of configuration space as quickly as possible without revisiting any region too soon. This maximizes entropy production over finite time — the system equilibrates fastest.

The Steinitz / three-distance theorem makes this concrete: for the sequence {nα mod 1}, the gaps between successive points take at most 3 distinct values. For φ, these gaps have lengths in ratio φ:1 — the gaps themselves are self-similar. No other irrational produces gaps that are simultaneously maximally uniform and self-similar.

Line 3: Self-referential closure (eigenvalue uniqueness)

Steps 1–8 derived that the eigenvalue of Φ = ∂_tΦ is φ. This third line closes the selection by proving φ is the only eigenvalue that satisfies all three conditions simultaneously:

No other real number is simultaneously the fixed point of R = 1 + 1/R, the last KAM torus, and the optimal equidistribution constant. φ is not “chosen” by the system. It is the only thing left after stability, coverage, and self-consistency eliminate everything else.

Corollary. If spacetime scales are constrained to the dominant eigenmode, a lattice with eigenvalue φ has characteristic distances at d_n = φⁿ in natural units. The natural unit is 1 light year.

Calibration. Earth’s orbit corresponds to shell n = −23. This is not asserted — it is derived: log_φ(1 AU in ly) = −22.97, rounded to the nearest integer. The algorithm is detailed in Section III.

III. The Shell Algorithm

Shell indices are not chosen to fit — they are derived by a deterministic procedure. Given any measured distance d (in light years):

There is no freedom in this assignment. The index for each object is dictated by its measured distance and the rounding operation. The residual |n* − n| measures how close the object sits to an integer shell.

Reference Constants

IV. The Ratio Theorem

The ratio between any two shell distances is a pure dimensionless number:

This formulation makes the “smuggled constants” objection structurally impossible. It is φ raised to a power, compared to a quotient of two measurements. Both measurements use the same unit, which cancels.

Ratio Test 1 — Neptune / Earth

Shells −16 and −23  ·  solar system scale

Ratio Test 2 — M87* / Proxima Centauri

Shells +37 and +3  ·  spanning local to extragalactic

Ratio Test 3 — Observable Universe / Heliopause

Shells +51 and −13  ·  heliospheric to cosmological

Ratio tests have zero free parameters. φ⁷ = 29.034 is a mathematical fact. Neptune-to-Earth = 30.07 is an astronomical fact. Their agreement to 3.5% is either coincidence or structure.

V. Dataset

The eight objects were chosen by a single rule: the most precisely measured landmark at each distance decade.

The categories are chosen by physical significance, not by proximity to φ-shells. Once the category is fixed (e.g. “nearest star”), the object is determined by nature, not by the analyst.

^† An earlier draft tested Cygnus A (757 Mly, shell +42, residual 0.324, error 20.8%). M87* was substituted because it is the most precisely measured landmark at the extragalactic decade — the Event Horizon Telescope’s direct imaging constrains its distance to ±2%, whereas Cygnus A’s radio-lobe estimates carry ~15% uncertainty. Cygnus A appears in the extended catalog (Section X) as a mediocre fit. The swap is justified by measurement precision, not by better alignment.

Deduction 1 — Earth’s Orbit

Shell n = −23  ·  calibration point

This is the calibration step. The 1.3% residual comes from choosing the integer shell n = −23. A non-integer calibration (n = −23.027) would give exact agreement — the framework insists on integers.

Deduction 2 — Neptune

Shell n = −16  ·  prediction (no parameter adjustment)

Deduction 3 — Voyager 1 / Heliopause

Shell n = −13  ·  prediction

Deduction 4 — Proxima Centauri

Shell n = +3  ·  prediction  ·  nearest star to the Sun

Quarter of one percent. No parameter was adjusted between the Earth calibration and this prediction. The formula simply gives φ³ ≈ 4.236.

Deduction 5 — Sagittarius A*

Shell n = +21  ·  prediction  ·  galactic center supermassive black hole

Largest error in the set. The measured uncertainty itself is ±1.6%, and independent methods range from 25,000 to 28,000 ly. The prediction falls within the broader literature spread.

Deduction 6 — M87*

Shell n = +37  ·  prediction  ·  first black hole ever imaged

Sub-one-percent at 54 million light years. Shell +37 is 60 shells away from the Earth calibration at −23. The prediction uses no information about M87* — only φ³⁷.

Deduction 7 — GW190521

Shell n = +49  ·  prediction  ·  most massive gravitational-wave merger detected

Deduction 8 — Observable Universe Radius

Shell n = +51  ·  prediction  ·  comoving radius of the particle horizon

The edge of the observable universe. Shell +51 — seventy-four shells from the calibration point. The entire visible cosmos, predicted to 2% by one free parameter and the golden ratio.

VI. Residual Analysis

Null hypothesis. If distances are randomly distributed in log space, the continuous index n* is uniform modulo 1, so residuals |n* − round(n*)| are Uniform(0, 0.5) with expected value 0.250. The framework predicts residuals near zero.

Under the stated null model, objects sit 5.1× closer to integer φ-shells than random placement predicts. Specifically, under the null hypothesis (8 i.i.d. Uniform(0, 0.5) residuals, E = 0.25, σ_mean = 0.051), the observed mean of 0.049 gives Z = −3.94,  p < 4 × 10⁻⁵. Even the worst outlier (Sgr A*, |r| = 0.177) sits within the 0.5 maximum. Remove it and the remaining 7 average 0.030 — eight times closer than chance.

Completeness Test: All Major Planets

What happens when we include every planet, not just the two that appear above?

Result: The 8 planets together show mean residual 0.194 — closer than random (0.250) but not statistically significant (p ≈ 0.13). Venus and Jupiter are far from any shell. This is expected. The framework predicts structural scales — distance scales where matter preferentially accumulates — not that every individual planet sits on a φ-shell. The test in the residual table above asks: “Are the characteristic scales of the universe (nearest-star distance, galactic radius, cosmic horizon) structured?” — not “Is every planet at φⁿ?”

VII. Statistical Test

The Z-score in Section VI assumes residuals are independent Uniform(0, 0.5) draws. We can confirm this without distributional assumptions by direct simulation.

Monte Carlo Simulation

Generate 10,000 sets of 8 random distances drawn uniform in log space over the observed range (10⁻⁵ to 10^10.7 ly). For each set, compute the mean residual. Count how many achieve ≤ 0.049.

On Independence

The Z-score assumes the 8 residuals are independent draws. This is reasonable if the objects were selected by category (they were) and their distances are not physically linked. Earth’s orbit does not constrain the distance to M87*. The one partial dependency: Neptune’s distance is gravitationally related to the solar system’s structure, making it non-independent of Earth. If Neptune is excluded and the remaining 7 objects tested, the mean residual drops to 0.041 — the result strengthens.

VIII. Predictions & Explorer

The strongest defense against lookback bias is prediction. The framework makes testable claims about objects not in the primary dataset:

Two near shells, two far off. This is honest. The framework predicts structural scales, not that every astronomical landmark sits on one. The informative test is statistical: across a large enough pre-registered sample, do residuals cluster near zero more than chance predicts? Expanding to 20+ objects from independent catalogs (Hipparcos nearby stars, Messier objects, GWTC events) would constitute a stronger test than any 4–8 object analysis.

Shell Explorer

IX. Falsification Criteria

What would falsify the framework?

A systematic trend in the residuals (errors growing with |n|, or consistently biased in one direction) would indicate curve-fitting rather than structure. Additionally, if future high-precision measurements of Sgr A* distance (currently the largest outlier at 8.2%) move away from 24,476 ly rather than toward it, the framework loses its weakest prediction.

What about objects that don’t sit on integer shells?

Most objects don’t. The claim is not that every object in the universe occupies a perfect φⁿ — it is that characteristic structural distances (modal scales where objects preferentially form or accumulate) follow this spacing. Individual planets between shells are expected. The test is whether the nodes of the structure — the scales where things cluster — match φⁿ.

How many free parameters?

One. The assignment n = −23 for Earth’s orbit. This fixes the scale factor. Every other prediction follows from φⁿ with no adjustment. The ratio tests in Section IV have zero free parameters. Compare: ΛCDM uses six free parameters and does not predict the Earth–Sun distance at all.

How are shell numbers assigned?

By algorithm (Section III). Each object is assigned the nearest integer n such that φⁿ (in light years) approximates the measured distance. No rounding is done in the object’s favor — n is the closest integer, even when the next integer would give a smaller error.

Could this be coincidence?

The residual analysis in Section VI and Monte Carlo in Section VII provide the formal test. For 8 landmark objects, Z = −3.94,  p < 4 × 10⁻⁵. Under the stated null model, the clustering is statistically significant at better than 99.99% confidence. The all-planets test (mean |r| = 0.194, p ≈ 0.13) shows that individual objects are NOT expected to fit — the signal is in the structural landmarks.

Why φ? What is the operator?

The identity Φ = ∂_tΦ identifies the operator: ∂_t. The field is its own time derivative — the eigenmodes are solutions where the field reproduces itself under time evolution. φ is the eigenvalue because it is the unique positive fixed point of R = 1 + 1/R. In a closed self-referential system, rational eigenvalues produce resonance lock — periodic orbits that destructively interfere. φ, being maximally irrational (continued fraction [1;1,1,1,…], slowest to converge), is the only eigenvalue that survives without annihilating itself. The system does not “prefer” φ. It is the only thing left.

Pre-Registered Blind Predictions

Every test above is retrospective: objects were known before shells were computed. To move from numerology to falsifiable physics, the framework must predict before measurement. The following predictions are stated in advance of the data they test against. Each specifies the lattice, the scale range, and the expected outcome.

Protocol. Each prediction above was stated before the corresponding data was consulted. The lattice (φⁿ ly), the scale range, and the falsification threshold are all fixed. No post-hoc adjustment is permitted. If 3 or more of the 5 predictions fail their falsification threshold, the framework’s physical postulate is rejected.

Why complete catalogs matter. The central challenge is not the mathematics of φ — that is well established. The challenge is demonstrating that the universe selects this scaling basis independently of an observer choosing interesting objects. Every retrospective test in this paper is vulnerable to the objection that structurally significant objects were chosen because they happen to sit on shells. The blind predictions above address this by specifying complete, externally-defined catalogs (Abell 1989, GWTC-3/4, Gaia DR3/DR4, BeSSeL survey) where every object is included regardless of residual. If the signal survives a catalog the predictor did not curate, the selection objection falls.

Predictions 1, 4, and 5 are testable against existing data that was not consulted during formulation. Predictions 2 and 3 require future observations. All five are needed: retroactive checks using “unseen” existing data test the model, but only genuinely future data eliminates the possibility of unconscious selection.

Dokuchaev Interior Orbits and the Incommensurability Condition

The operator identification ∂_t acquires concrete physical content through Dokuchaev’s work on stable periodic orbits inside rotating black holes (Dokuchaev 2011, arXiv:1103.6140). Inside the inner Cauchy horizon of a Kerr-Newman black hole, there exist “orbits of the third kind” — completely bound trajectories that neither escape nor fall into the central singularity. These orbits are periodic in three independent coordinates (r, θ, φ) with periods T_r, T_θ, T_φ.

The critical property: these three periods are incommensurable. All ratios T_r/T_θ/T_φ are irrational. The 3D orbit never closes. This is not incidental — it is the stability condition. If the periods were commensurate (rational ratios), the orbit would close, tracing repeated passes through the same spacetime points, producing resonant energy buildup that destabilizes the trajectory.

This is exactly the mechanism the framework identifies. The identity Φ = ∂_tΦ says the field reproduces itself under time evolution. Rational eigenvalues produce resonance lock — periodic orbits that destructively interfere. Irrational eigenvalues avoid this. And among all irrational numbers, φ is the most irrational: its continued fraction [1;1,1,1,…] converges more slowly than any other, making it the hardest to approximate by rationals.

The KAM theorem (Kolmogorov-Arnold-Moser) formalizes this. In Hamiltonian systems, quasi-periodic orbits with sufficiently irrational frequency ratios survive small perturbations. The most robust are those whose ratios are noble numbers — numbers whose continued fraction expansion ends in all 1s. φ = [1;1,1,1,…] is the most noble number. It is the last frequency ratio standing.

The black hole interior is the ultimate closed system — causally disconnected from the external universe by two horizons. No information enters or leaves the R-region between singularity and inner Cauchy horizon. This is precisely the boundary condition under which self-referential eigenvalue selection (Φ = ∂_tΦ) should operate most cleanly. The system has no external forcing; it must sustain itself from its own dynamics. The orbits that survive are those with the most irrational period ratios. The most irrational ratio is φ.

Note what this does not claim. It does not claim that Dokuchaev’s specific computed examples have φ-ratio periods. His examples are particular solutions at specific (E, L, Q) values chosen for illustration. The framework’s prediction is about the structure of the solution space — that stability domains cluster around noble-number frequency ratios — not that any single orbit has periods in golden ratio.

Quasinormal Modes: The Direct Eigenmode Test

The shell-distance catalog tests whether φ appears in spatial structure. A more direct test of the operator identification ∂_t would look at the actual dynamical eigenmodes of spacetime itself.

Black hole quasinormal modes (QNMs) are exactly this. When perturbed, a Kerr black hole rings down at complex frequencies ω_n = ω_R + iω_I that are determined entirely by the black hole’s mass and spin — no free parameters. These are the eigenvalues of ∂_t acting on perturbations of the Kerr metric. They are computed, not observed-then-fitted. The overtone spectrum (n = 0, 1, 2, …) gives a tower of frequency ratios.

Similarly, Kerr geodesic frequencies — the orbital frequency Ω_φ, radial epicyclic frequency Ω_r, and vertical epicyclic frequency Ω_θ — are the three fundamental frequencies of motion in the Kerr metric. In the eikonal limit (high multipole l), QNM frequencies reduce to these geodesic frequencies evaluated at the light ring:

These frequency ratios Ω_r/Ω_θ and Ω_r/Ω_φ vary continuously with the spin parameter a and orbital radius. At the ISCO (innermost stable circular orbit) — the boundary between stable and unstable geodesics — the ratios take specific values determined by a alone.

The framework predicts:

This test is fundamentally different from the shell-distance catalog. Shell distances test whether φ organizes spatial scales — a claim vulnerable to observer selection. QNM frequency ratios test whether φ organizes temporal eigenmodes — the actual solutions to ∂_t acting on spacetime. There is no observer selection: the QNM spectrum is a property of the metric, computable from first principles. If the signal is there, it was always there. If it is not, the ∂_t operator interpretation is wrong.

X. Extended Catalog

The primary dataset (Section V) uses 8 landmarks chosen by physical significance. To test whether the lattice signal persists across a broader, less curated sample, we compute shell indices for 120 objects spanning 90 orders of magnitude — from Earth’s radius to the observable universe.

Selection rule: include every well-measured object at each scale decade whose distance is known to better than ±15%. No object is excluded for having a large residual. The catalog is deliberately non-selective — it includes known poor fits (Jupiter, Sirius, Pluto) alongside strong ones.

Notation

Earth & Solar Geometry

Solar System

Nearby Stars

Exoplanet Host Stars

Nebulae

Milky Way & Satellite Galaxies

Galaxies & Clusters

Large-Scale Structure & Cosmological

Catalog Statistics

The extended catalog is deliberately weaker than the primary dataset — it includes individual objects (planets, dwarf galaxies) rather than only structural landmarks. The signal attenuates as expected: the primary 8 landmarks give Z = −3.94; the full 120-object catalog gives Z = −1.75. The lattice is a structural phenomenon, not a universal law governing every individual distance — but even the non-selective sample shows excess clustering at integer shells.

Notable Discoveries on Lattice Lines

Objects not in the primary dataset that fall within 3% of an integer shell:

JWST-Era & O4 Gravitational Wave Additions

Objects discovered or precisely measured since 2022, added without selection on residual:

Distance note. Distances in this catalog use the values most commonly cited in the astronomical literature at the time of compilation (July 2026). Gaia DR3 (2022) refined several nearby measurements — in most cases by <5%, well within the lattice’s shell spacing. Where Gaia corrections are large (e.g. Pleiades moved from 392 → 446 ly; Ring Nebula from 2,283 → 2,570 ly), the catalog uses the Gaia value. Some earlier sources may report different distances for the same object; the lattice test is robust to ±10% distance uncertainty because the shell spacing is multiplicative (φ ≈ 1.618, a 62% gap between shells).

Terrestrial Scales: Ancient Sites & Planetary Geometry

The lattice extends below the astronomical catalog into planetary-scale distances. Great-circle distances between ancient monumental sites — selected for cultural significance, not for lattice fit — can be tested against the same φⁿ ly grid.

The distances below are WGS-84 great-circle calculations between site coordinates. Conversion: 1 ly = 9.461×10¹² km.

Six of six distances tested fall within |Δ| < 0.12 of integer shells. But the remarkable result is not the count — it is which distances land closest: Giza–Easter Island (the longest pair in the ancient-sites set, joining two civilisations separated by the Pacific) sits 0.12% from φ×10,000 km, landing on shell n = −42 with |Δ| = 0.048.

φ-Ratios Between Site Distances

The Ratio Theorem (Section IV) predicts that distances on adjacent shells have ratio φ. Several ancient-site pairs exhibit this:

The Machu Picchu–Karnak / Angkor Wat–Karnak ratio recovers φ to five significant figures. These three sites form a triangle on the globe whose longest two edges stand in the golden ratio.

Statistical caveat. Eighteen ancient sites yield 153 pairwise great-circle distances and 11,628 possible ratios of pairs. At this combinatorial depth, some near-φ hits are expected by chance — roughly 20% of random distances land within |Δ| < 0.10 of an integer shell. What is not expected is that the tightest shell hit (|Δ| = 0.048) belongs to the single most culturally prominent pair (Giza–Easter Island), or that the tightest ratio hit (0.008% from φ) connects three independent civilisations (Khmer, Inca, Egyptian). A rigorous test would require pre-registration of which pairs to examine.

φⁿ ly  —  from Earth’s radius to the edge of the observable universe.