Fibonacci Laws of Planetary Motion — Technical Background

What this page shows

Every planet has two key orbital properties that change slowly over hundreds of thousands of years:

• Eccentricity — how elongated the orbit is (0 = perfect circle, 1 = extreme ellipse). Earth’s is currently about 0.017, meaning its orbit is nearly circular.
• Inclination — how tilted the orbit is relative to the average plane of the Solar System. All planets wobble up and down slightly over long timescales; the oscillation amplitude measures how far each planet tilts.

When each value is multiplied by $\sqrt{m}$ (where $m$ is the planet’s mass), both eccentricities and inclinations snap onto Fibonacci ratios across all eight planets. The $\sqrt{m}$ weighting is not arbitrary — it is the unique exponent that arises in the Angular Momentum Deficit (AMD), the conserved quantity governing long-term orbital stability. In AMD-natural variables ($\xi = e\sqrt{m}$ for eccentricity, $\eta = i\sqrt{m}$ for inclination), ratios between planets become simple Fibonacci fractions.

The six laws documented below are:

1. Law 1 — Fibonacci Cycle Hierarchy: dividing $H$ by successive Fibonacci numbers yields the major precession periods
2. Law 2 — Inclination constant: $d \times \text{amp} \times \sqrt{m} = \psi$ with a single universal constant $\psi$ and pure Fibonacci divisors $d$
3. Law 3 — Inclination balance: angular-momentum-weighted amplitudes cancel between two phase groups (100% balance with dual-balanced eccentricities)
4. Law 4 — Eccentricity constant: AMD partition ratios satisfy a small Fibonacci/Lucas constraint per mirror pair, predicting the four outer-planet eccentricities from the four inner-planet eccentricities (mean error 0.11%)
5. Law 5 — Eccentricity balance: an independent balance condition on eccentricities using the same divisors and phase groups (100% balance with dual-balanced eccentricities)
6. Law 6 — Saturn-Jupiter-Earth Resonance: a closed beat-frequency loop linking the three dominant precession periods

The universal constant $\psi$ can be computed purely from Fibonacci numbers and the Holistic-Year $H =$ 335,317:

From this single constant, all eight planets’ inclination amplitudes are predicted with zero free parameters. Combined with the eccentricity laws, the model connects orbital shapes and tilts of all eight planets through Fibonacci arithmetic rooted in a single timescale — the Holistic-Year $H =$ 335,317 years, which also determines each planet’s precession period (Law 1; see also Precession).

Prediction summary

Precession periods — Law 1 shows that dividing $H$ by successive Fibonacci numbers yields all major precession periods. The beat frequency rule $1/H(n) + 1/H(n+1) = 1/H(n+2)$ is an algebraic identity from the Fibonacci recurrence.

Inclination amplitudes — predicted from the single $\psi = 2205/(2H)$ with pure Fibonacci divisors. All 8 planets fit within their Laplace-Lagrange secular theory bounds. The configuration is the unique solution (Config #1) satisfying all four physical constraints out of 7,558,272 tested — see Law 2 and Finding 2 for details.

All eight eccentricities — Law 4 predicts the four outer-planet eccentricities from the four inner-planet eccentricities via four pair-specific Fibonacci/Lucas constraints on the AMD partition ratio $R = e_\text{base}/i_\text{mean}$, with zero free parameters beyond the inner quartet (mean prediction error 0.11%, max 0.28% at Saturn). Inner planets also form a Fibonacci ladder $\xi_V : \xi_E : \xi_{\text{Ma}} : \xi_{\text{Me}} = 1 : 5/2 : 5 : 8$, where $\xi = e \times \sqrt{m}$. See Law 4 and Finding 5 for details.

Predictive system status

Input Data

Masses and semi-major axes are derived from the Holistic computation chain ($H =$ 335,317, orbit counts, Kepler’s third law). Eccentricities are J2000 values except Earth, where the oscillation-midpoint base eccentricity (0.015386) from the 3D simulation is used — see Finding 5. Inclinations are to the invariable plane at J2000 (Souami & Souchay 2012). Inclination amplitudes are derived from the Fibonacci formula $\text{amp} = \psi / (d \times \sqrt{m})$.

Assignments

Planet Configuration

Phase Angles and Balance Groups

Each planet has a per-planet phase angle — the ICRF perihelion longitude at the balanced year (~302,635 BC). These cluster near Laplace-Lagrange eigenmodes ($\gamma_1$–$\gamma_8$) within 1–10°. Saturn is anti-phase: MAX inclination at balanced year, while all other planets are at MIN.

The group assignment is constrained by: (1) each planet’s oscillation range must fall within the Laplace-Lagrange secular theory bounds, (2) the inclination structural weights must balance (Law 3), and (3) the eccentricity weights must balance (Law 5). Phase angles cluster near LL eigenmodes ($\gamma_1$–$\gamma_8$) within 1–10°.

Precession Periods

Law 1: The Fibonacci Cycle Hierarchy

Dividing the Holistic-Year $H =$ 335,317 by successive Fibonacci numbers yields the major precession periods of the solar system.

Beat frequency rule

The beat frequencies of consecutive Fibonacci-divided periods satisfy the Fibonacci recurrence:

This is an algebraic identity, not an empirical fit. For example:

• $1/111{,}296 + 1/66{,}778 = 1/41{,}736$ (i.e., $3 + 5 = 8$)
• $1/66{,}778 + 1/41{,}736 = 1/25{,}684$ (i.e., $5 + 8 = 13$)
• $1/41{,}736 + 1/25{,}684 = 1/15{,}899$ (i.e., $8 + 13 = 21$)

The empirical content of Law 1 is that the precession periods of Earth ($H/3$), Jupiter ($H/5$), and Saturn ($-H/8$, ecliptic-retrograde) are related to a single timescale $H$ through Fibonacci denominators. The higher-order periods ($H/13$, $H/21$, $H/34$) then follow as beat frequencies. See also Precession for the physical meaning of each cycle.

Law 2: Inclination Amplitude

Each planet’s mass-weighted inclination amplitude, multiplied by its Fibonacci divisor, equals the universal $\psi$-constant:

Equivalently:

This holds for all 8 planets with a single universal $\psi$ = 3.288 × 10⁻³.

Why the symbol $\psi$? In quantum mechanics, $\psi$ denotes the wave function whose squared modulus gives the probability of finding a particle in a given state. Here, $\psi^2$ appears in the AMD energy formula $\sqrt{a}\,\psi^2/(2d^2)$, where the Fibonacci divisor $d$ plays the role of a quantum number governing each planet’s share of inclination oscillation energy — see Physical meaning of $\psi$. The analogy is structural, not physical.

Predictions and Laplace-Lagrange validation

The equation $d \times \text{amp} \times \sqrt{m} = \psi$ holds by construction. The non-trivial test is that every predicted range (mean $\pm$ amplitude) falls within the Laplace-Lagrange secular theory bounds — and all 8 do.

Worked example: Earth’s inclination amplitude

Earth has Fibonacci divisor $d = 3$ ($= F_4$). Step by step:

The mean is computed from the J2000 constraint:

Fibonacci vs IAU 2006-optimized values

The 3D simulation uses a slightly different value optimized for the IAU 2006 obliquity rate:

The 0.47% difference is within the model’s tolerance. The Fibonacci prediction represents the theoretical long-term value from the balance condition, while the IAU 2006-optimized value is calibrated to the currently observed obliquity rate.

Physical meaning of $\psi$: mass-independent AMD partition

Substituting the Law 2 formula into the standard Angular Momentum Deficit gives each planet’s inclination oscillation AMD:

Mass cancels exactly. Law 2 appears to depend on mass — heavier planets tilt less. But that apparent mass dependence is precisely what removes mass from the energy budget: $\sqrt{m}$ is the unique exponent that makes this cancellation work. Each planet’s inclination oscillation energy depends only on its orbital distance ($\sqrt{a}$) and Fibonacci quantum number ($d$), not on how massive the planet is. This is not an approximation — it is an algebraic identity (see Why $\sqrt{m}$).

The formula separates into three independent factors:

Summing over all planets and solving for $\psi$ gives a budget equation that determines $\psi$‘s magnitude:

This means $\psi$ is fixed by two things: the total inclination oscillation energy (set by formation physics) and the geometric sum $\sum \sqrt{a}/d^2$ (set by Fibonacci structure and Kepler spacing). Neither can be changed independently.

The $1/d^2$ scaling gives $d$ the role of a quantum number: each planet is assigned a Fibonacci slot that determines its share of the eight-planet inclination oscillation energy:

Saturn dominates (56%) not because of mass (Jupiter is $3\times$ heavier) but because $d = 3$ combined with a large orbit maximizes $\sqrt{a}/d^2$. Earth carries more oscillation energy than Jupiter (18% vs 15%) despite being $1000\times$ lighter — its $d = 3$ beats Jupiter’s $d = 5$ in the $1/d^2$ scaling. The Earth–Saturn pair ($d = 3$) carries 74% of the eight-planet total, and the E–J–S resonance triad (Law 6) carries 89% — consistent with their role as the dynamically dominant subsystem. (Shares are relative to the eight major planets, which carry 99.994% of the system’s orbital angular momentum; TNOs contribute the remainder.)

Law 3: Inclination Balance

The angular-momentum-weighted inclination amplitudes cancel between the two phase groups, conserving the orientation of the invariable plane:

Physical motivation

The invariable plane is the fundamental reference plane of the solar system, perpendicular to the total angular momentum vector. For this plane to remain stable, the angular-momentum-weighted inclination oscillations must cancel between the two balance groups. If they didn’t, the total angular momentum vector would wobble — violating conservation of angular momentum.

Substituting the Fibonacci formula

With $\text{amp}_j = \psi / (d_j \times \sqrt{m_j})$ and $L_j = m_j \sqrt{a_j(1 - e_j^2)}$:

Since $\psi$ is a single universal constant, it cancels from both sides:

where $w_j = \sqrt{m_j \times a_j(1 - e_j^2)} / d_j$ is the structural weight for each planet.

Structural weights

Jupiter ($d = 5$) contributes the dominant in-phase weight ($1.408 \times 10^{-2}$). Saturn ($d = 3$) alone carries the entire anti-phase contribution. The remaining six in-phase planets collectively contribute $3.29 \times 10^{-3}$ to match the Saturn–Jupiter difference. The eccentricity enters through the angular momentum factor $\sqrt{1 - e^2}$; the dual-balanced outer-planet eccentricities are the values at which this balance becomes exact.

TNO contribution

The balance considers only the 8 major planets, which carry 99.994% of the solar system’s orbital angular momentum. Trans-Neptunian Objects (TNOs) contribute the remaining ~0.006%, tilting the invariable plane by approximately 1.25″ (Li, Xia & Zhou 2019 ).

Law 4: The Eccentricity Constant

Within each mirror pair, the AMD partition ratio $R = e_\text{base} / i_\text{mean,rad}$ satisfies a small Fibonacci/Lucas constraint. The four pair constraints predict the four outer-planet base eccentricities from the four inner-planet base eccentricities to within 0.3%.

For each planet, $R$ measures how it partitions its orbital deviation between eccentricity and inclination — a high $R$ means most deviation is in orbit shape, a low $R$ means most is in orbit tilt. Inputs are base eccentricities (long-term midpoints from the eccentricity base/amplitude split) and Law-2 mean inclinations (radians).

AMD partition ratio per planet

This connects to AMD theory (Laskar 1997): $C_k \approx \Lambda_k(e^2/2 + i^2/2)$ for small $e$ and $i$. The ratio $R^2 + 1 = (e^2 + i^2)/i^2$ determines how each planet partitions its angular momentum deficit between eccentricity and inclination.

Fibonacci/Lucas pair constraints

Within each mirror pair, $R$ satisfies one Fibonacci/Lucas constraint:

Three of the four constraints are pure ratios; only the Mercury/Uranus relation is a sum-of-squares. Three of the four targets are mixed Fibonacci/Lucas ratios; only Mercury/Uranus uses pure Fibonacci. Mean constraint error: 0.126%. Max: 0.282%.

Each constraint is one equation, so the four pair relations together predict the four outer-planet eccentricities (Jupiter, Saturn, Uranus, Neptune) from the four inner-planet eccentricities (Mars, Earth, Venus, Mercury) — leaving zero free parameters beyond the inner quartet, which itself satisfies the Finding 6 ξ-ladder.

Predictions — zero free parameters

Each pair’s constraint is solved for the outer planet’s $R$ value, then converted back to eccentricity via $e_\text{pred} = R_\text{pred} \times i_\text{mean,rad}$:

• Jupiter ← Mars: $R_\text{Ju} = R_\text{Ma} \times \sqrt{144/11}$
• Saturn ← Earth: $R_\text{Sa} = R_\text{E} \times 21/4$
• Uranus ← Mercury: $R_\text{Ur} = \sqrt{11 - R^2_\text{Me}}$
• Neptune ← Venus: $R_\text{Ne} = R_\text{V} \times \sqrt{55/4}$

Mean prediction error: 0.106%. Max: 0.281% (Saturn).

The reference and observed values are base eccentricities (long-term midpoints), not J2000 snapshots. The Saturn prediction is the loosest — driven by the only linear-ratio relation among the four — yet still predicts Saturn’s base eccentricity to better than three parts in a thousand from Earth alone.

Cross-check with Law 5. Feeding the Law-4-predicted set into the eccentricity balance equation gives 99.86% balance — within 0.14 percentage points of the 99.9993% achieved by the observed values. The two laws are mutually consistent. Independence is confirmed by Monte Carlo: among 50,000 random eccentricity sets satisfying Law 5 (with Saturn $\in [0,0.5]$), zero simultaneously satisfy all four Law 4 pair constraints to within 1%.

Dynamical implication — communicating vessels. Each pair constraint implies that the two planets’ eccentricities cannot evolve independently. For Earth/Saturn the linear ratio $R_\text{Sa} = (21/4) R_\text{E}$ means Saturn’s eccentricity is locked to Earth’s: as Earth’s eccentricity decreases (currently heading toward its ${\sim}0.0139$ minimum around 11,725 AD), Saturn’s must decrease in proportion. The coupling is bidirectional: Saturn has its own perihelion precession cycle ($-H/8$ = 41,915 years, ecliptic-retrograde), which drives its own eccentricity changes that feed back into Earth’s. The physical coupling is already present in the precession formulas: Saturn is the only planet whose precession formula requires Earth’s time-varying obliquity and eccentricity as inputs (GROUP 15 terms; see Formulas). Whether the same communicating-vessel behaviour holds for the other three mirror pairs (Mars/Jupiter, Venus/Neptune, Mercury/Uranus) is a further test of this framework.

This law partially resolves Open Question 1: while no universal closed-form eccentricity formula exists, all eight eccentricities are predicted from inclinations via Fibonacci/Lucas pair constraints on the AMD partition ratio.

Law 5: Eccentricity Balance

The eccentricities satisfy an independent balance condition using the same Fibonacci divisors and balance groups:

where

Or equivalently, in terms of orbital period $T_j \propto a_j^{3/2}$:

Eccentricity balance weights

Saturn alone carries the entire anti-phase contribution. The in-phase group is dominated by Jupiter ($7.928 \times 10^{-3}$), Uranus ($5.705 \times 10^{-3}$), and Neptune ($1.734 \times 10^{-3}$), with the four inner planets contributing only $6.6 \times 10^{-5}$ combined.

Comparing the two balance conditions

The half-power difference in Fibonacci divisor scaling ($1/d$ vs $1/\sqrt{d}$) and the shift from $\sqrt{a}$ to $a^{3/2}$ reflect that the eccentricity balance operates at a different order in secular perturbation theory.

The eccentricity balance is not a structural artifact — see Finding 3 for three independence tests.

Connection to AMD theory

The established AMD formula of Laskar (1997) couples eccentricity and inclination through a single conserved quantity: $C_k = m_k \sqrt{\mu a_k}(1 - \sqrt{1-e_k^2}\cos i_k)$. For small $e$ and $i$: $C_k \approx \Lambda_k(e^2/2 + i^2/2)$.

The eccentricity balance operates on linear $e$ rather than $e^2$, suggesting it captures a first-order secular constraint distinct from the quadratic AMD conservation. The physical mechanism producing this linear balance remains an open question.

Law 6: The Saturn-Jupiter-Earth Resonance

Saturn’s ecliptic-retrograde precession creates a closed beat-frequency loop with Jupiter and Earth — the physical mechanism linking the Fibonacci timescale (Law 1) to orbital structure (Laws 2–5).

Saturn is the only planet whose perihelion precesses retrograde in the ecliptic frame (period $= -H/8$). JPL’s WebGeoCalc confirms the ecliptic-retrograde motion at ~−3,400 arcsec/century from SPICE ephemeris data, consistent with the model’s geocentric prediction of −3,385 arcsec/century. Standard celestial mechanics attributes the observed ecliptic-retrograde motion to a transient phase of the ~900-year Great Inequality (Laplace 1784), predicting the rate should eventually reverse. The model instead treats Saturn’s ecliptic-retrograde precession as a permanent feature — a testable distinction. See Supporting Evidence §14 for the full observational evidence and comparison of the two theories.

When Saturn’s retrograde motion interacts with Jupiter’s prograde motion ($H/5$), the resulting beat frequencies form a closed loop:

All three rows are cyclic permutations of $3 + 5 = 8$. Each planet’s period is the beat frequency of the other two:

               Saturn (H/8)
              ╱            ╲
          8−5=3          8−3=5
            ╱                ╲
    Earth (H/3) ──3+5=8── Jupiter (H/5)

This also connects to Law 1: combining Jupiter ($H/5$) and Saturn ($H/8$) produces axial precession: $F_5/H + F_6/H = F_7/H$ = 1/25,794 (i.e., $5 + 8 = 13$).

Why this is a Fibonacci identity

The closed loop is an algebraic consequence of the Fibonacci addition rule applied to the period denominators:

Since $F_n/H + F_{n+1}/H = F_{n+2}/H$, the beat frequency rule is just the Fibonacci recurrence in disguise. The empirical content is that Earth, Jupiter, and Saturn have period denominators $b_E = 3$, $b_J = 5$, $b_S = 8$ — three consecutive Fibonacci numbers with Saturn as the sum.

Connection to the ψ-constant

The same three planets and their period denominators determine the inclination constant:

Law 6 thus reveals why Earth, Jupiter, and Saturn hold special roles throughout the Fibonacci Laws: they are the three planets locked into a Fibonacci resonance triangle, and their period denominators drive both the timescale hierarchy (Law 1) and the inclination constant (Law 2).

Findings

Findings are empirical observations and consequences that follow from the six laws.

Finding 1: Mirror Symmetry

As shown in the Planet Configuration table, each inner planet shares its $d$ with its outer counterpart across the asteroid belt: Mars ↔ Jupiter ($F_5$), Earth ↔ Saturn ($F_4$), Venus ↔ Neptune ($F_9$), Mercury ↔ Uranus ($F_8$). Earth–Saturn is the only pair with opposite balance groups (in-phase vs anti-phase). The divisors form two consecutive Fibonacci pairs: $(3, 5)$ for the belt-adjacent planets and $(21, 34)$ for the outer planets.

Finding 2: Configuration Uniqueness

The exhaustive search tested 7,558,272 configurations against four independent physical constraints:

Each filter independently eliminates most of the search space. The most restrictive three-filter intersection (Saturn-solo + balance + Laplace-Lagrange) yields only 5 configurations. Adding mirror symmetry leaves exactly one: Config #1 (0.0000132% of the search space).

Since Earth is locked at $d = 3$, only Scenario A (Saturn $d = 3$) can satisfy the Earth↔Saturn mirror symmetry. No other planet configuration produces the required mirror pairing while simultaneously satisfying all physical constraints.

The search space covers:

• Fibonacci divisors: $d \in \{1, 2, 3, 5, 8, 13, 21, 34, 55\}$ for Mercury, Venus, Mars, Uranus, Neptune
• Balance group: in-phase or anti-phase for each planet (2 options each). Each planet’s phase angle is its own ICRF perihelion longitude at the balanced year, which determines whether it reaches MIN (in-phase) or MAX (anti-phase) inclination at that moment.
• 4 scenarios for Jupiter and Saturn (fixed per scenario): A: Ju=5, Sa=3 — B: Ju=8, Sa=5 — C: Ju=13, Sa=8 — D: Ju=21, Sa=13
• Earth: locked at $d = 3$, in-phase group, phase angle ~21.77° (Earth’s own ICRF perihelion at balanced year)
• Total: $9 \times 2 \times 9 \times 2 \times 9 \times 2 \times 9 \times 2 \times 9 \times 2 \times 4 = 7{,}558{,}272$ configurations

Finding 3: Eccentricity Balance Independence

The eccentricity balance (Law 5) is genuinely independent from the inclination balance (Law 3). Three tests confirm it is not a structural artifact:

1. Coefficient test: The weight formula without eccentricities ($\sqrt{m} \times a^{3/2} / \sqrt{d}$) gives only 74% balance — the actual eccentricity values contribute 26 percentage points of improvement
2. Random test: Substituting random eccentricities into the same weight formula gives 50–85% balance across 1,000 trials
3. Power test: The balance peaks sharply at $e^{1.0}$ (100% with dual-balanced eccentricities), dropping to 98.5% for $e^{0.9}$ and 98.4% for $e^{1.1}$, and to 91% for $e^2$ — linear eccentricity is special

Additionally, the ratio $v_j/w_j$ varies by a factor of ~150 across planets, confirming the two balance conditions are not proportional.

Finding 4: Saturn Eccentricity Prediction and Law Convergence

Since Saturn is the sole anti-phase planet, the eccentricity balance directly predicts its eccentricity from the other seven:

Two structurally independent laws predict Saturn’s eccentricity from the same model inputs:

Law convergence: Laws 4 and 5 agree on Saturn’s eccentricity to within 0.28%. Law 4 derives it from the linear ratio $R_\text{Sa}/R_\text{E} = 21/4$ ($F_8/L_3$) using only the Earth-Saturn pair. Law 5 derives it from the global balance equation involving all eight planets simultaneously.

Why this is significant: Saturn’s eccentricity oscillates secularly between ~0.01 and ~0.09 (a factor-of-9 dynamic range). Two structurally different Fibonacci/Lucas constraints — one using only the Earth-Saturn pair, the other using all eight planets — converge on the same value within 0.28%. The Fibonacci divisors were originally chosen to match precession periods (Law 1) and inclination balance (Law 3); the eccentricity predictions were never optimized for.

Overconstrained system: The 4 pair constraints (Law 4) plus Law 5 provide 5 equations for 4 unknowns (the four outer-planet eccentricities, given the four inner ones). The system is overconstrained by one equation — Saturn’s eccentricity is predicted by both Law 4 (from Earth) and Law 5 (from the global balance) independently. That two equations drawing on different subsets of planetary data converge confirms that the eccentricity balance is not an independent free parameter but an emergent consequence of the Fibonacci/Lucas pair structure.

Finding 5: Inner Planet Eccentricity Ladder

The mass-weighted eccentricities of the four inner planets form a Fibonacci ratio sequence:

where $\xi = e \times \sqrt{m}$. Consecutive ratios are $5/2$, $2$, $8/5$ — all converging toward the golden ratio $\varphi \approx 1.618$.

Earth’s base eccentricity (0.015386), independently determined from the 3D simulation, anchors the ladder at position $5/2$:

Finding 6: No Universal Eccentricity Base Constant

A comprehensive search across ~30 physical and Fibonacci parameters confirms that no combination $\xi \times f(\text{observables})$ yields a planet-universal constant for base eccentricities analogous to the inclination $\psi$.

The search tested all power-law combinations of: semi-major axis, orbital period, mass, inclination precession period, perihelion precession period, AMD, Hill radius, orbital angular momentum, secular eigenfrequencies ($g$- and $s$-modes), Kozai integral, mean inclination, eccentricity half-angle $h(e)$, Fibonacci divisors ($d$, $b$), and cross-parameter ratios — across 1-, 2-, 3-, and 4-parameter combinations with fractional exponents.

Results (excluding tautologies that restate the inclination law):

The fundamental reason: the inclination divisors $d$ span only 11× (from 3 to 34), allowing a single $\psi$ to work. The eccentricity multipliers $k$ span 141× (from 1 to 141) with no smooth dependence on any present-day observable. The eccentricity ladder is a formation constraint imprinted during the dissipative epoch, not derivable from current orbital mechanics.

Finding 7: Universal Eccentricity Amplitude Constant K

While no universal constant exists for base eccentricities, the eccentricity amplitude (the oscillation range around the base) does follow a universal formula:

where $K = 3.4149 \times 10^{-6}$, derived from Earth’s mean obliquity (23.41354°) and eccentricity amplitude (0.001356). This is the eccentricity analog of $\psi$ for inclination amplitudes.

Eccentricity at any time:

where $\theta = 360° \times (t - t_0) / T_{\text{cycle}}$ and $h_1 = \sqrt{e_{\text{base}}^2 + e_{\text{amp}}^2} - e_{\text{base}}$.

Candidate relations to $\psi$ (numerical coincidences, not proven identities):

The $\psi^{11/5}$ relation uses $L_5 = 11$ (5th Lucas number) and $F_5 = 5$ (5th Fibonacci number) — pure Fibonacci/Lucas structure, but 11 is not itself a Fibonacci number. $K$ remains an empirical constant pending further investigation.

Inner/outer regime distinction. The K formula applies differently across the asteroid belt. For inner planets (Venus, Earth, Mars), the tilt-driven mechanism dominates: inverting the formula predicts obliquities to within 0.03°. For outer planets (Jupiter–Neptune), the predicted amplitudes are 3–149× smaller than the observed J2000-to-base differences — their eccentricity oscillations are dominated by Laplace–Lagrange secular eigenmode exchange, not the tilt mechanism. K determines the tilt-coupled component; secular eigenmodes provide the remainder.

Balance preservation. The Law 5 weight change from eccentricity oscillation simplifies to $\delta v = K \times \sin(\text{tilt})$ — mass and distance cancel completely. This ensures the eccentricity balance is preserved at every epoch.

J2000 eccentricity phase angles (for the $e(t)$ formula):

Inner planet phases are analytically derived (Earth: $\phi = \omega + 90° = 102.95° + 90°$). Outer planet phases are set by secular eigenmode structure.

Deriving $\psi$ from the Holistic Year

The inclination constant $\psi$ was initially a calibration parameter — measured from the data but not derived from first principles. A systematic search reveals it can be expressed exactly in terms of $H$ and Fibonacci numbers:

This predicts $\psi$ = 3.288 × 10⁻³. Combined with the Fibonacci divisors and the balance condition, this single value determines all eight inclination amplitudes — and every predicted range falls within the Laplace-Lagrange secular theory bounds.

Physical interpretation: the E–J–S period connection

The Fibonacci indices 5 and 8 are the period denominators of Jupiter and Saturn (from Law 1: $H/5$ and $H/8$). The denominator factor 2 = $F_3 = F_{b_E}$ is Earth’s contribution:

where $b_E = 3$, $b_J = 5$, $b_S = 8$ are the period denominators of Earth, Jupiter, and Saturn, satisfying the Fibonacci addition $3 + 5 = 8$.

Saturn’s $F_b$ appears squared, reflecting its special role as the “sum planet” ($b_S = b_E + b_J$). These same three planets dominate the inclination balance condition — Jupiter and Saturn are the dominant contributors to the in-phase and anti-phase groups respectively.

Equivalently, $\psi \times T_{V/N} = F_5 \times F_8^2 = 2205$, where $T_{V/N} = 2H$ is the period shared by Venus and Neptune. This connects the amplitude structure (how much planets tilt) to the period structure (how fast they oscillate).

The Master Ratio $R = 311$

Since $\psi$ is fixed by $H$ and Fibonacci numbers, the eccentricity scale is determined by a single ratio:

Venus’s base eccentricity ($e_V$ = 0.006757) is the value that makes $R$ exactly 311 — the Fibonacci primitive root prime nearest to the formation-determined scale. This simultaneously achieves 100% eccentricity balance (Law 5) and places Venus within 0.3% of its J2000 observed eccentricity. Since 311 is prime and not expressible as a Fibonacci product, the two structures remain connected but not unified. A Fibonacci primitive root prime has the Fibonacci sequence mod $p$ with the maximal possible period (π(311) = 310), ensuring maximum compatibility with the Fibonacci architecture. The eccentricity ladder is a formation constraint ($p < 10^{-5}$ from Monte Carlo), making $R$ the model’s single irreducible parameter.

For the full analysis of why 311 — including the two-layer explanation, the TRAPPIST-1 connection, and the FPR prime properties — see Physical Origin §3.

Physical Foundations

Why $\sqrt{m}$: the Angular Momentum Deficit connection

The quantity $\xi = e \times \sqrt{m}$ (or $\eta = i \times \sqrt{m}$) has a direct connection to the Angular Momentum Deficit (AMD), the standard dynamical quantity in celestial mechanics:

where $\xi_i = e_i \sqrt{m_i}$ and $\eta_i = i_i \sqrt{m_i}$ are exactly the mass-weighted quantities used in the Fibonacci laws. The $\sqrt{m}$ exponent is uniquely determined by the AMD decomposition — no other mass power produces this sum-of-squares form.

Empirical confirmation: Testing the $\psi$-constant ($d \times i \times m^\alpha$) across Venus, Earth, and Neptune for different mass exponents:

The exponent $\alpha = 0.50$ is optimal by a factor of $\sim 250\times$ over the next-best alternative.

A deeper consequence: substituting Law 2 into the AMD formula causes mass to cancel entirely, revealing that the Fibonacci divisor $d$ acts as a quantum number governing mass-independent energy partition — see Physical meaning of $\psi$ under Law 2.

Connection to KAM theory

The appearance of Fibonacci numbers in mass-weighted orbital parameters has a theoretical explanation through KAM (Kolmogorov-Arnold-Moser) theory:

1. KAM theorem (Kolmogorov 1954, Arnold 1963, Moser 1962): Orbits survive if their frequency ratio satisfies a Diophantine condition — a quantitative measure of irrationality.

2. Golden ratio optimality (Hurwitz 1891): The golden ratio $\varphi = [1; 1, 1, 1, \ldots]$ has the slowest-converging continued fraction, making it the hardest irrational number to approximate by rationals. Fibonacci ratios $F_{n+1}/F_n$ are its convergents.

3. Greene-Mackay criterion (Greene 1979): The last invariant torus to break as perturbation increases has frequency ratio equal to $\varphi$. Morbidelli & Giorgilli (1995) proved golden-ratio KAM tori have superexponentially long stability times.

4. Natural selection over 4.6 Gyr: Orbits near low-order rational resonances are disrupted (Kirkwood gaps). Orbits near golden-ratio KAM tori survive. The surviving population is organized around Fibonacci ratios.

5. AMD variables: The quantities $\xi = e\sqrt{m}$ and $\eta = i\sqrt{m}$ are the natural action-like variables in the secular Hamiltonian. KAM stability conditions constrain these to Fibonacci ratios.

This chain — KAM stability → Diophantine condition → golden ratio → Fibonacci convergents → Fibonacci structure in $\xi$ and $\eta$ — provides the theoretical foundation for why Fibonacci numbers specifically appear.

Time independence

The Fibonacci laws use oscillation amplitudes (secular eigenmode properties), not instantaneous orbital elements. The $\psi$-constants and balance conditions use amplitudes $\eta = i_\text{amp} \times \sqrt{m}$ — properties of secular eigenmodes that are exactly time-independent by construction. The eccentricity ladder acts on base eccentricities that encode a formation-epoch constraint, preserved in the eigenmodes since the dissipative era.

Statistical Significance

Molchanov’s 1968 work was criticized by Backus (1969) for not proving statistical significance. To address this, a comprehensive significance analysis was performed using three independent null models and fourteen test statistics covering all six laws, explicitly accounting for the look-elsewhere effect.

Methodology. Fourteen test statistics were computed for the real Solar System and compared against random planetary systems:

1. Pairwise Fibonacci count — How many of the 56 pairwise $\xi$-ratios and $\eta$-ratios match Fibonacci ratios within 5%? Observed: 18/56.
2. Eccentricity ladder (Finding 5) — For the best 4-planet subset, how many $\xi$-ratios match Fibonacci ratios within 3%? Observed: 2 matches (mean error 3.35%).
3. $\psi$-constant (3-planet) — How small is the relative spread of $d \times i \times \sqrt{m}$ for the best 3-planet subset with Fibonacci weights? Observed: 16.97%.
4. Cross-parameter ratio — Does $(\Lambda_\text{ecc} + \Lambda_\text{amp}) / \Lambda_\text{mean}$ match an extended Fibonacci ratio? Observed: 2.16% error.
5. Inclination balance (Law 3) — With the model’s Fibonacci $d$-values, phase groups, and base eccentricities, how well do the structural weights $w_j$ cancel? Observed: 100%.
6. Eccentricity balance (Law 5) — With the same $d$-values, phase groups, and base eccentricities, how well do the eccentricity weights $v_j$ cancel? Observed: 100%.
7. Saturn eccentricity prediction (Finding 4) — The eccentricity balance predicts Saturn’s eccentricity from the other 7 planets. Observed: 0.00% error (base eccentricities satisfy balance by construction).
8. $\psi$-constant, full 8-planet (Law 2) — How small is the spread of $d \times i_\text{J2000} \times \sqrt{m}$ across all 8 planets using J2000 inclinations? Observed: 234.7% spread.
9. Precession hierarchy (Law 1) — How many pairwise ratios of precession periods $H/F_n$ match Fibonacci ratios? Observed: 10 pairs.
10. Pair constraints (Law 4) — For each of the four mirror pairs, does the pair-specific $R$ combination (3 ratios + 1 sum-of-squares, $R = e_\text{base}/i_\text{mean,rad}$) match a Fibonacci/Lucas target? Observed: 4/4 pairs match within 0.3%.
11. E–J–S resonance (Law 6) — Do the period denominators satisfy $b_E + b_J = b_S$? Observed: exact ($3 + 5 = 8$).
12. Mirror symmetry (Finding 1) — How many of the 4 inner–outer mirror pairs share the same $d$-value? Observed: 4/4 pairs match.
13. K amplitude constant (Finding 6) — Can a single constant $K$ predict all 8 eccentricity amplitudes via $e_\text{amp} = K \times \sin(\text{tilt}) \times \sqrt{d} / (\sqrt{m} \times a^{1.5})$? Observed: 0.00% max error (amplitudes derived from $K$; permutation test is still valid).
14. Eccentricity Balance Scale (Finding 7) — Can each planet’s base eccentricity be predicted from the other 7 using the balance scale weights? Observed: 0.02% RMS error.

Tests 5–7 use the model’s fixed Fibonacci $d$-assignments and phase groups — zero look-elsewhere effect. Tests 1–4 involve implicit optimization (searching over subsets, weight assignments, or Fibonacci ratios), but this is accounted for automatically because each random trial computes the same optimized statistics. Tests 10–12 depend on planet identity (mirror-pair assignments, period fractions, and $d$-assignments), making them insensitive to the permutation test but testable via Monte Carlo.

Three null distributions were tested:

• Permutation (exhaustive, $8! = 40{,}320$ trials): same eccentricity and amplitude values, randomly reassigned to planets. For Tests 5–7: model $d$-values are kept, eccentricities are shuffled among planets. Tests 11–12 are identity-dependent and always match in this null ($p = 1$).
• Log-uniform Monte Carlo ($100{,}000$ trials): random eccentricities from $[0.005, 0.25]$, amplitudes from $[0.01°, 3.0°]$, means from $[0.1°, 10°]$ (log-uniform). Random $d$-values from $\{1, 2, 3, 5, 8, 13, 21, 34, 55\}$, random precession periods, random period fractions from $\{1, 2, 3, 5, 8, 13, 21, 34\}$.
• Uniform Monte Carlo ($100{,}000$ trials): same ranges, flat distribution. Same null model for Tests 5–12.

Masses are fixed at solar system values throughout (conservative).

Results (permutation: exhaustive $8!$; Monte Carlo: $100{,}000$ trials):