Physical Origin — Why Fibonacci?
The Fibonacci Laws document what the solar system’s structure is. This page addresses why it exists. The answer has two parallel parts at different scales. At the planetary-spacing scale (Laws 1–5), the architecture rests on KAM theory acting on configurations frozen during formation and preserved since by conservative dynamics (§§1–3 below). At the 8H climate-lattice scale (the integer divisors 8H/N that organize Earth’s obliquity record), the secular eigenfrequencies are forced by action-angle closure of the obliquity sector (§4 below). KAM and action-angle closure do not compete — they act on different scales.
For the laws themselves, see Fibonacci Laws. For derivations, see Fibonacci Laws Derivation. For broader scientific context, see Scientific Background.
1. The Mathematical Foundation
The KAM theorem (Kolmogorov 1954, Arnold 1963, Moser 1962) proves that in perturbed dynamical systems, orbits with golden-ratio frequency ratios are the most resistant to destabilization. Fibonacci ratios (3/2, 5/3, 8/5, 13/8, …) are the best rational approximations to the golden ratio, so KAM theory favours them as the most stable discrete frequencies available to real orbits. Greene (1979) confirmed computationally that the golden-ratio torus is the last to break; Morbidelli & Giorgilli (1995) showed super-exponential stability near these ratios in the asteroid belt.
The Fibonacci patterns emerge when orbits are described in AMD-natural variables — eccentricity × √m and inclination × √m. The √m weighting is uniquely dictated by the Angular Momentum Deficit (AMD), the conserved quantity governing long-term orbital stability. At exactly α = 0.50 the Fibonacci structure achieves a spread of 0.11%; the next-best exponent gives more than 28%. See fibonacci_amd_structure.py for the AMD decomposition.
One intrinsic motion, the rest coupling-driven. Of Earth’s precession cycles, only axial precession (H/13) is physically intrinsic — it arises from the lunisolar torque on Earth’s equatorial bulge and would persist with no other planets present. Inclination precession (H/3), ecliptic precession (H/5), and the obliquity and eccentricity oscillations are all sustained by gravitational coupling with Jupiter and Saturn (see Law 6). That an intrinsic lunisolar motion and a planet-driven motion lock into a 13:3 Fibonacci ratio — the seed of the entire H hierarchy — is exactly the kind of KAM-stable configuration the dissipative disk phase would select.
2. Formation-Epoch Freezing
KAM theory explains why Fibonacci ratios are preferred — but not why the balance conditions are satisfied so precisely (99.9974% inclination, 99.8636% eccentricity). The precision comes from the solar system’s formation epoch.
Phase 1 — Protoplanetary disk (0–10 Myr). Planets form and migrate within a gas disk. Dissipative forces (gas drag, disk torques, tidal interactions) continuously push orbits toward configurations that minimize AMD. AMD measures how far orbits deviate from perfectly circular and coplanar; dissipation naturally reduces it.
Phase 2 — KAM selection. Among all AMD-minimizing configurations, those organized by Fibonacci ratios have the widest stability margins. Dissipative evolution therefore preferentially converges toward Fibonacci-organized configurations — not because Fibonacci is imposed, but because Fibonacci configurations are the deepest stability wells in the energy landscape.
Phase 3 — Disk dissipation (~3–10 Myr). When the gas disk dissipates, the dissipative mechanism shuts off and the Fibonacci configuration is frozen — like a ball settling into the deepest valley and the landscape hardening around it.
Earth’s base eccentricity as a formation anchor. Earth’s base eccentricity (e_E = 0.015386) is the model’s single irreducible eccentricity parameter — an oscillation midpoint, not the J2000 value (which is 0.01671, about 8.6% higher). The other seven base eccentricities derive from K, J2000 observations, and the System Reset anchor (n = 7) with balance-group phase offsets. Earth’s base is the only free eccentricity parameter in the model — consistent with the formation-epoch picture: the dissipative era selected a specific Fibonacci-compatible configuration with Earth’s base at this value, and conservative dynamics has preserved it since.
3. Origin Story — Planetary Spacing
| Phase | Mechanism | What it sets | Timescale |
|---|---|---|---|
| Mathematical foundation | KAM: golden ratio is the most irrational, therefore most stable | Why Fibonacci numbers specifically | Timeless |
| Dissipative selection | Disk evolution minimizes AMD while KAM selects Fibonacci wells | Specific Fibonacci divisors | 3–10 Myr |
| Long-term preservation | Conservative Hamiltonian dynamics; KAM tori protect the architecture | Why the structure persists at J2000 | 4.5 Gyr |
The analogy: a landscape with many valleys. The deepest valleys (hardest to escape) are those shaped by the golden ratio (KAM). During formation, a river (dissipative gas disk) carried material toward the lowest points. When the river dried up, everything was left sitting in the deepest valleys. Billions of years of wind and rain (gravitational perturbations) have shifted material within each valley, but nothing has been enough to push anything over a ridge. The Fibonacci Laws describe the geometry of those deepest valleys. The climate lattice — the integer divisors 8H/N recording the secular eigenfrequencies — is a parallel story at a different scale, with its own mechanism, covered next.
4. Mechanism for the Climate Lattice
KAM theory and formation-epoch freezing explain planetary-spacing stability (Laws 1–5). But the model’s climate lattice — the Solar System Resonance Cycle 8H = 2,682,536 yr and the integer divisors 8H/N that organize Earth’s obliquity and eccentricity record — is a different scale and requires a different mechanism. The empirical search identified the mechanism: action-angle closure of the obliquity sector, not KAM resonance protection.
Chirikov-KAM resonance overlap was not the mechanism
The natural first hypothesis was that the 8H lattice marks the surviving KAM-stable spectral structure — that 8H/N integers should sit in gaps between overlapping resonance widths (Chirikov’s criterion). The test was null. Across 33 spectral peaks in LA2004 eccentricity + obliquity (Welch PSD, 51-Myr solution), every peak is KAM-stable — sharp spectral peaks already imply quasi-periodic dynamics, by construction. Lattice-aligned peaks have median Chirikov K = 0.584; non-lattice peaks have median K = 0.516. The difference is not significant (Mann-Whitney p = 0.667), and the direction is opposite the KAM-overlap hypothesis prediction.
KAM stability is necessary but not sufficient: the 8H lattice is in the stable manifold, but the stable manifold contains more than the 8H lattice. Whatever selects the lattice specifically must be a more specific principle.
Action-angle closure — strongly supported in the obliquity sector
If 8H is a true closed-orbit period of the secular system, then in action-angle coordinates the obliquity trajectory must return to its starting point every 8H, and the obliquity eigenfrequencies must be integer multiples of 1 / (8H). Three direct tests on LA2004 51-Myr data confirm this — strikingly, in the obliquity sector only:
| Sub-test | Eccentricity vector (h, k) | Obliquity ε(t) |
|---|---|---|
Closure distance D(8H) vs 200 random lags | 57.5 % of random lags beat it → null | only 2.5 % of random lags beat it → ✓ |
Top 10 spectral peaks at integer divisors of 1 / (8H) | 6 / 10 within 5 % | 8/10 within 5 %, most within 0.2 % |
The obliquity result is direct: top spectral peaks have frequencies that hit integer divisors of 1 / (8H) with sub-percent precision (n = 66, 67, 68, 69, 70, 91, 94, …). The eccentricity vector does not close at 8H — but this is predicted. Earlier work (Test C-Invariant) found obliquity is 100 % on the lattice while eccentricity is only 74 %. The off-lattice 26 % is a Mercury-chaos-driven perturbation, which breaks closure for the eccentricity vector but not for obliquity (consistent with Laskar 1989/1994 inner-planet chaos restricted to eccentricity).
Why this forces integer divisibility of the lattice
If 8H is the closed-orbit period and the eigenfrequencies are commensurate at this period, then each eigenfrequency must take the form n / (8H) for some integer n. The spectrum is forced onto integer divisors of 8H — not coincidence, conservation. The climate-active subset of those integers (the L1 lattice used in the model’s climate fit) is then a selection criterion within an already-quantized set, not a free choice. This is the conservation law the lattice needs: integer divisibility falls out of action-angle periodicity.
Two mechanisms, two scales
| Scale | Mechanism | What it explains |
|---|---|---|
| Planetary spacing (Laws 1–5) | KAM theorem + formation-epoch freezing (§§1–3) | Why orbital architecture sits at Fibonacci-stable ratios |
8H climate lattice (8H/N integer divisors) | Action-angle closure of the obliquity sector | Why secular eigenfrequencies are forced onto n / (8H) |
Both mechanisms persist as the lattice evolves: the integer labels stay invariant across geological time while the absolute periods rescale — see Expanding Resonance for the time-evolution layer.
Scripts: chirikov_resonance_test.py, action_closure_test.py.
5. Hierarchy of Certainty
Not all claims carry equal evidence. An honest assessment:
| Claim | Evidence | Status |
|---|---|---|
| Fibonacci numbers appear in solar system orbital architecture | p = 1.5 × 10⁻⁴–1.0 × 10⁻⁵ (joint permutation test; 4 empirical tests in permutation null, 9 in MC nulls) → 3.62–4.26σ | Established |
| The structure acts on AMD-natural variables (√m weighting) | 0.11% spread at α = 0.5 vs >28% at any other exponent | Established |
| KAM theory provides the mathematical foundation | Proven theorem + computational confirmation | Established framework |
| The structure was set during formation and frozen at disk dissipation | N-body confirms no secular mode preserves it | Well-supported |
8H is a closed-orbit period of the obliquity sector; action-angle closure forces integer divisibility of 8H/N | Test A: only 2.5 % of random lags match the observed closure; 8/10 top obliquity peaks land within 5 % of integer divisors of 1/(8H); Chirikov-KAM null (p = 0.67) | Established (climate-lattice mechanism) |
The bottom line: The existence of non-trivial Fibonacci structure beyond chance is supported at 3.62–4.26σ. The planetary-spacing mechanism — KAM theory selecting golden-ratio-compatible orbits, frozen at formation — is well-supported by theory and simulation. The climate-lattice mechanism — action-angle closure of the obliquity sector — is directly tested with sub-percent precision (Test A). The specific values (H, ψ, K, and the value of 8H itself) are empirically determined and not yet derived from first principles. This is comparable to where Kepler’s laws stood before Newton: the patterns are real and the mechanisms are identified, but the deepest “why these specific numbers” awaits a further theoretical breakthrough.
6. References
- KAM theory: Kolmogorov 1954, Arnold 1963, Moser 1962
- Golden ratio stability: Greene 1979, Mackay 1983
- Planetary KAM: Morbidelli & Giorgilli 1995
- AMD: Laskar 1997, Laskar & Petit 2017
- Secular theory: Brouwer & van Woerkom 1950
- Chirikov resonance overlap: Chirikov 1979
- Mercury chaos in eccentricity sector: Laskar 1989, Laskar 1994
Investigation scripts on GitHub (filenames prefixed fibonacci_, chirikov_, action_closure_).
Summary
| Question | Answer |
|---|---|
| Why Fibonacci (planetary spacing)? | Golden ratio = most stable frequency ratio (KAM); Fibonacci = its discrete skeleton |
| Why so precise? | Dissipative formation drove orbits into the deepest Fibonacci stability wells |
| Why integer divisors of 8H (climate lattice)? | Action-angle closure — the obliquity sector is periodic at 8H, forcing eigenfrequencies onto n / (8H) |
| Why still true today? | KAM tori (planetary) + action-angle periodicity (climate) protect their respective lattices |
| What’s still unknown? | Why the value 8H = 2,682,536 yr specifically; why this subset of integer divisors becomes climate-active; Earth’s special role |