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: Earth’s major precession periods divide $H$ by Fibonacci numbers — H/3, H/5, H/8, H/13 — a hierarchy unique to Earth
2. Law 2 — Inclination constant: $d \times \text{amp} \times \sqrt{m} = \psi$ with a single universal constant $\psi$ and pure Fibonacci divisors $d$ (all eight planets)
3. Law 3 — Inclination balance: angular-momentum-weighted inclination amplitudes of seven planets balance against Saturn alone (99.9975% with phase-derived base eccentricities)
4. Law 4 — Eccentricity amplitude constant: $e_\text{amp} = K \times \sin(\text{mean obliquity}) \times \sqrt{d} / (\sqrt{m} \times a^{1.5})$ with a single universal constant $K$ derived from Earth (all eight planets)
5. Law 5 — Eccentricity balance: mass- and distance-weighted eccentricities of seven planets balance against Saturn alone — using the same Fibonacci divisors and phase groups as Law 3 (99.8632% with phase-derived base eccentricities)
6. Law 6 — Saturn-Jupiter-Earth Resonance: Earth’s $H/8$ obliquity cycle equals both Jupiter’s ICRF perihelion and Saturn’s ecliptic perihelion — a triple identity at $H/8$ maintained by the Jupiter-Saturn mutual resonance lock

Two empirical constants, both derived from Earth, predict all oscillation amplitudes:

From these two constants, all eight planets’ inclination and eccentricity amplitudes are predicted with zero free parameters. Combined with the balance laws, the model connects orbital shapes and tilts of all eight planets through Fibonacci divisors rooted in a single timescale — the Earth Fundamental Cycle $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 Earth’s major precession periods are $H/F$ for Fibonacci $F$. 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 constant $\psi$ (derived from Earth) with pure Fibonacci divisors. All 8 planets fit within their Laplace-Lagrange secular theory bounds. The configuration is the unique mirror-symmetric solution (Config #11) among 42 surviving a five-stage selection pipeline from 7,558,272 tested — see Law 2 and Finding 2 for details.

Eccentricity amplitudes — predicted from the single constant $K$ (derived from Earth) using Fibonacci divisors, mass, distance, and axial tilt. See Law 4.

Saturn’s base eccentricity — predicted by Law 5 from the other seven planets via the eccentricity balance equation (~0.27% accuracy). The remaining seven base eccentricities are derived from the System Reset phase (n=7) — see Finding 4.

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. 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

Cycle Anchors and Balance Groups

Each planet has a per-planet cycle anchor — the ICRF perihelion longitude where MAX inclination occurs, evaluated at the balanced year (302,635 BC). Saturn is anti-phase: its cosine sign is flipped relative to all other planets. The full anti-phase alignment — Saturn at maximum while all others at minimum — occurs once per Solar System Resonance Cycle (8H).

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).

Precession Periods

Law 1: The Fibonacci Cycle Hierarchy

Earth’s major precession periods divide the Earth Fundamental Cycle $H =$ 335,317 by Fibonacci numbers — $H/3$ (inclination), $H/5$ (ecliptic), $H/8$ (obliquity), $H/13$ (axial). The Fibonacci addition rule connects them: $3 + 5 = 8$, $5 + 8 = 13$.

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{,}772 + 1/67{,}063 = 1/41{,}915$ (i.e., $3 + 5 = 8$)
• $1/67{,}063 + 1/41{,}915 = 1/25{,}794$ (i.e., $5 + 8 = 13$)
• $1/41{,}915 + 1/25{,}794 = 1/15{,}967$ (i.e., $8 + 13 = 21$)

The empirical content of Law 1 is that Earth’s major precession periods — inclination ($H/3$), ecliptic ($H/5$), obliquity ($H/8$), axial ($H/13$) — are all related to a single timescale $H$ through Fibonacci denominators. The higher-order periods ($H/21$, $H/34$) then follow as beat frequencies. The coincidence of Earth’s $H/8$ obliquity with Jupiter’s ICRF perihelion period and Saturn’s ecliptic perihelion period is the subject of Law 6. 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.3068 × 10⁻³.

Why the symbol $\psi$? A notation choice — $\psi^2$ appears in the AMD energy formula $\sqrt{a}\,\psi^2/(2d^2)$, where the Fibonacci divisor $d$ takes the role of an integer index, an algebraic form reminiscent of how integer quantum numbers appear in quantum-mechanical energy partitions — see Physical meaning of $\psi$. The resemblance is purely formal; the model makes no claim about quantum physics.

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:

Substituting Earth’s values:

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.1%) 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.2% vs 14.9%) 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 phase-derived base eccentricities are the values at which this balance approaches exactness.

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 Amplitude Constant

A single universal constant $K$ predicts all eight eccentricity oscillation amplitudes from Fibonacci divisors, mass, distance, and axial tilt:

$K =$ 3.4149 × 10⁻⁶, derived from Earth’s eccentricity amplitude and axial tilt. This is the eccentricity analog of $\psi$ (Law 2) for inclination amplitudes.

Both $\psi$ and $K$ are empirical constants derived from Earth, and both predict all 8 planets with zero free parameters. $K$ additionally uses semi-major axis $a$ and axial tilt, coupling the spin and orbital domains. This coupling is bidirectional: a planet’s mean obliquity determines its eccentricity amplitude, and its eccentricity amplitude reveals its obliquity. Together, $\psi$ and $K$ link all three oscillation parameters — inclination amplitude, eccentricity amplitude, and axial tilt — through just two constants.

Worked example: Earth’s eccentricity amplitude

Earth has Fibonacci divisor $d = 3$, mean obliquity 23.41354°, and semi-major axis $a = 1.000$ AU. Step by step:

Verification — plugging $K$ back into the formula recovers the input amplitude:

The same $K$ value then predicts all other planets’ eccentricity amplitudes — using each planet’s own obliquity, $d$, mass, and semi-major axis — with zero additional free parameters.

Note on base eccentricities: Laws 2 and 4 predict oscillation amplitudes. The base (mean) eccentricities are partially constrained by Law 5, which predicts Saturn’s from the other seven. The remaining seven base eccentricities are derived from the System Reset phase (n=7) with balance-group phase offsets (90° in-phase, 270° Saturn) — see Finding 4.

Law 5: Eccentricity Balance

The mass- and distance-weighted eccentricities of seven planets balance against Saturn’s alone — using the same Fibonacci divisors and phase groups as Law 3:

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 × 10⁻³), Uranus (5.705 × 10⁻³), and Neptune (1.733 × 10⁻³), with the four inner planets contributing only 6.5 × 10⁻⁵ combined.

The ~0.14% residual

The eccentricity balance reaches 99.8632% but not exactly 100%. The residual may reflect contributions from minor bodies (dwarf planets, asteroids) not included in the 8-planet framework, or measurement uncertainties in planetary masses — particularly Uranus and Neptune, whose masses are currently known only from Voyager 2 flybys (relative uncertainties ~0.02–0.08%). A future Uranus or Neptune orbiter providing more precise mass measurements should improve the balance — a testable prediction.

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

Earth’s $H/8$ obliquity cycle equals both Jupiter’s ICRF perihelion period and Saturn’s ecliptic perihelion period — a triple identity at $H/8$ by which the gas giants gravitationally drive Earth’s spin-axis dynamics. Jupiter and Saturn maintain this configuration through their well-known mutual resonance lock.

Saturn’s perihelion precesses obviously retrograde in the ecliptic frame (period $= -H/8$). JPL’s WebGeoCalc confirms this at ~-3,400 arcsec/century from SPICE ephemeris data, consistent with the model’s prediction of -3,372 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 §13 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

Law 6 reveals why Earth, Jupiter, and Saturn hold special roles throughout the Fibonacci Laws: they are the three planets locked into a Fibonacci resonance triangle ($3 + 5 = 8$). The period denominators 3, 5, and 8 are the first three terms of Earth’s H/Fibonacci hierarchy (Law 1), tying the gas-giant resonance to Earth’s precession structure.

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 through a pipeline of successively stricter physical constraints:

Per-config optimisation. The Laplace–Lagrange bounds and direction-trend checks depend on two parameters that are not determined by the $d$-values alone: the anchor position $n \in \{0, \ldots, 7\}$ (where in the 8H Solar System Resonance Cycle the balanced year falls — this sets all cycle anchors) and the ascending node integer $N$ per planet ($\Omega$ regression period $= -8H/N$). Each of the 96 candidates passing both balance thresholds is evaluated at its own optimal $(n, N)$, making the comparison fair. Jupiter and Saturn are constrained to share the same $N$, preserving their lockstep ascending node regression.

Of the 42 survivors, only one is mirror-symmetric: Config #11 (0.0000132% of the search space). Since Earth is locked at $d = 3$, only Scenario A (Saturn $d = 3$) can produce the Earth↔Saturn mirror pair. It ranks #11 of 42 by eccentricity balance (99.8632%).

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 inclination cycle anchor is its ICRF perihelion longitude at the System Reset (where the planet reaches its inclination extremum).
• 4 scenarios for Jupiter and Saturn: 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, cycle anchor ~21.77°
• Total: $9^5 \times 2^5 \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}$ (99.86%), dropping substantially for other powers — 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:

Law 5 predicts Saturn’s eccentricity from the other seven planets:

Why this is significant: Saturn’s eccentricity oscillates secularly between ~0.01 and ~0.09 (a factor-of-9 dynamic range). The Fibonacci divisors were originally chosen to match precession periods (Law 1) and inclination balance (Law 3); the eccentricity prediction was never optimized for. Law 4 predicts all eight eccentricity amplitudes; Law 5 additionally predicts Saturn’s base eccentricity to ~0.27% of the observed value. The remaining seven base eccentricities are phase-derived at runtime from K, J2000 observations, and the System Reset anchor (n=7) with balance-group phase offsets (90° in-phase, 270° Saturn).

Finding 5: 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.

Law 4 Details: The K Constant

The K constant formula and derivation are described in Law 4 above. Additional details:

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}}$.

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). Each planet’s phase is determined by its position in the eccentricity oscillation cycle at J2000, using the System Reset (n=7) as the common timing anchor with a balance-group phase offset (90° for in-phase planets, 270° for Saturn). At n=7, every planet simultaneously passes through its mean eccentricity — in-phase planets rising, Saturn falling — mirroring the inclination alignment:

The wobble period is the beat frequency of each planet’s axial precession and ICRF perihelion precession. The phase determines the base eccentricity (oscillation midpoint) via the law of cosines with the J2000 observed eccentricity and K-derived amplitude.

The Universal Constants $\psi$ and $K$

Both $\psi$ and $K$ are empirical constants derived from Earth’s fitted parameters:

Given $\psi$, a planet’s mass, and its Fibonacci divisor $d$, the inclination amplitude is fully determined. Given $K$, the same parameters plus distance and axial tilt determine the eccentricity amplitude.

The ratio $\psi / K \approx 968 = F_6 \times L_5^2 = 8 \times 11^2$ (0.04% error) is a striking near-integer, but whether this reflects deeper structure is an open question.

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’s 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. Base eccentricities are phase-derived at runtime from the K amplitude, J2000 observations, and the System Reset anchor (n=7) with balance-group phase offsets, encoding 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 eleven test statistics covering all six laws, explicitly accounting for the look-elsewhere effect.

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

1. Law 1 — Fibonacci denominators — The 8 inclination-cycle periods have the form $T = H \times (a/b)$. How many of the denominators $b$ fall in the Fibonacci set $\{1,2,3,5,8,13,21,34,55\}$? Observed: 7/8 (Mercury’s $b=11$ is the only non-Fibonacci).
2. Law 2 — $\psi$ full 8-planet — How constant is $d \times \text{amp} \times \sqrt{m}$ across all 8 planets using the model amplitudes? Observed: 0.0000% spread (exact by construction of $\text{amp} = \psi / (d \sqrt{m})$; permutation still meaningful).
3. Law 3 — Inclination balance — With the model’s Fibonacci $d$-values, phase groups, and phase-derived base eccentricities, how well do the structural weights $w_j$ cancel? Observed: 99.9975%.
4. Law 4 — K amplitude constant — Can a single constant $K$ predict all 8 eccentricity amplitudes via $e_\text{amp} = K \times \sin(\text{mean obliquity}) \times \sqrt{d} / (\sqrt{m} \times a^{1.5})$? Observed: 0.0000% error (exact by construction of $K$; permutation test is still valid).
5. Law 5 — Eccentricity balance — With the same $d$-values, phase groups, and phase-derived base eccentricities, how well do the eccentricity weights $v_j$ cancel? Observed: 99.8632%.
6. Law 6 — E–J–S resonance — Do the period denominators satisfy $b_E + b_J = b_S$? Observed: exact ($3 + 5 = 8$).
7. Finding 1 — Mirror symmetry — How many of the 4 inner–outer mirror pairs share the same $d$-value? Observed: 4/4 pairs match.
8. Finding 1b — d-set Fibonacci clustering — Does the set of distinct $d$-values form exactly two consecutive Fibonacci pairs? Observed: yes — $\{3, 5, 21, 34\} = (F_4,F_5) \cup (F_8,F_9)$.
9. Finding 4 — Saturn eccentricity prediction — The eccentricity balance predicts Saturn’s base eccentricity from the other 7 planets. Observed: 0.27% error.
10. Finding 6 — Solo planet identification — For each planet, compute the residual $|v_p - \sum_{j\neq p} v_j|/v_p$ and take the minimum across all 8 planets. Saturn uniquely minimizes this. Observed: 0.27% min residual.
11. Year-length beat identity — Do the sidereal/tropical and anomalistic/tropical year beats match the Fibonacci integers $H/13$ and $H/16$? Observed: 0.0048% max error (both beats match to better than $5 \times 10^{-5}$).

Structural vs. empirical tests. 7 of the 11 tests are structural and yield $p = 1$ in the permutation null by construction. These fall into two groups.

Multiset-invariant (5 tests) — permutation cannot disturb them because the test sees only the multiset of values (or a single scalar):

• Law 1 (Fibonacci denominators — multiset invariant)
• Law 6 (E–J–S resonance — scalar integer identity)
• Finding 1 (Mirror symmetry — identity-dependent)
• Finding 1b (d-set clustering — multiset invariant)
• Year-length beat (single scalar)

Tautological (2 tests) — the model defines the test relation to hold identically:

• Law 2 ($\psi$ full 8-planet) — amplitudes are defined as $\psi/(d\sqrt m)$
• Law 4 ($K$ amplitude constant) — eccentricity amplitudes are defined from $K$

These 7 structural tests are reported in the table for transparency but excluded from the combined $p$-value. The remaining 4 tests (Laws 3, 5; Findings 4 and 6) are empirical — sensitive to permutation — and form the basis of the combined statistic. They are not statistically independent: all four use the same fundamental quantity $v_j = \sqrt{m_j}\,a_j^{3/2}\,e_j/\sqrt{d_j}$, so the combined $p$-value absorbs a correlation penalty (see below).

Three null distributions were tested:

• Permutation (exhaustive, $8! = 40\,320$ trials): same values, randomly reassigned to planets.
• Log-uniform Monte Carlo ($10^6$ 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 period denominators from $[2, 40]$, random year-length ratios.
• Uniform Monte Carlo ($10^6$ trials): same ranges, flat distribution.

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

Results (permutation: exhaustive $8!$; Monte Carlo: $10^6$ trials):

‡ Structural / multiset-invariant (Laws 1, 6; Findings 1, 1b; Year-length beat). Permutation yields $p = 1$ by construction, so these are excluded from the permutation combined statistic. Under the MC nulls (random solar systems), the tests are meaningful — random $d$-assignments, random period denominators, and random year lengths all produce distinct null distributions — so they are included in the MC joint test.

§ Tautological (Laws 2 and 4). Model amplitudes are defined as $\psi/(d\sqrt m)$ and $K$-derived, so the test relations hold identically under any null. Reported as internal-consistency checks; excluded from every combined statistic. See the warning callout above for why no permutation or MC null can test these.

Headline statistic: direct joint permutation test

The headline is a direct joint permutation test over the 4 empirical tests (Laws 3, 5; Findings 4, 6). Each test’s raw statistic is studentized against the null mean and standard deviation computed across all $8! = 40\,320$ permutations; the four studentized z-scores are summed into a single combined statistic $T = \sum_i z_i$; and the $p$-value is the fraction of null permutations with $T_\text{null} \geq T_\text{obs}$.

This is the most defensible combined statistic because:

1. Model-independent. No Gaussian, $\chi^2$, or distributional assumption is made. The joint null distribution of $T$ is computed explicitly from the 40,320 permutations.
2. Self-correcting for correlation. The four empirical tests all depend on the quantity $v_j = \sqrt{m_j}\,a_j^{3/2}\,e_j/\sqrt{d_j}$, so they are positively correlated. Because the joint null is computed from the same shared permutation that generates the correlation, no explicit correction factor is needed — the correlation is baked into the null by construction. For reference, the measured average pairwise correlation under the permutation null is $\bar r =$ 0.283 (replacing an earlier estimate of $\bar r \approx 0.5$ that proved too conservative). Under the MC nulls the measured correlation is much smaller ($\bar r \approx$ 0.032), because random $d$-assignments break the shared $v_j$ dependency.
3. No floor-clamp artifact. By construction $p \geq 1/(n+1)$, so the combined statistic never collapses to zero even for deeply significant results.

The joint test is computed under each of the three null distributions:

The MC combined statistics include five additional tests (Laws 1, 6; Findings 1, 1b; year-length beat) that are permutation-invariant but become meaningful under random solar systems. Laws 2 and 4 remain excluded under all nulls because they are tautological.

The permutation result (3.62$\sigma$) is the figure to cite. It makes no assumption about “what counts as a random planetary system” and is computed entirely from the observed planetary data. The Monte Carlo results (4.75$\sigma$ and 4.75$\sigma$) test the more general question — “would random systems reproduce this?” — and give a stronger answer (roughly 1 in 500,000 vs 1 in 10,000 for permutation). The two MC nulls agree to within $0.01\sigma$, indicating insensitivity to the choice between log-uniform and uniform sampling.

Comparison with supporting approximations. Two older combining methods are reported for transparency but are approximations of the joint test:

• Stouffer’s $Z$ with measured $\bar r$ (permutation: 1.8 × 10⁻⁵, $\approx$ 4.1$\sigma$). Converts each per-test $p$-value to a one-tailed $z$ and sums with a variance inflation factor $1 + (k-1)\bar r$. Uses the measured correlation 0.283 (previously hard-coded at 0.5, which was too conservative).
• Fisher’s method (permutation: 8.7 × 10⁻⁸). Combines log-$p$-values into a $\chi^2$ statistic. Highly sensitive to extreme small $p$-values and to the floor-clamp at $p = 1/n_{\text{trials}}$; reported as a cross-check only.

Both give stronger apparent results than the direct joint test, but both rely on Gaussian / $\chi^2$ distributional assumptions that the direct test avoids. The gap is a useful indicator of how much the distributional approximations are contributing versus the raw data.

Jackknife: robustness to dropping any single planet

A leave-one-out jackknife re-runs the direct joint permutation test (now over $7! = 5\,040$ permutations) with each planet removed from the system in turn. This quantifies how much of the combined signal rides on any one body.

Dropping Jupiter collapses the signal to $\approx 0.5\sigma$, and dropping Uranus collapses it to $\approx 0.9\sigma$. This is consistent with — and required by — the model’s structural claims: the $v_j$ balance is global (each side of the Fibonacci $d$-ladder must contribute), so removing any high-$v$ planet destroys the balance. The pattern matches what the six laws predict: a cooperative structure in which all eight planets participate, with the outer planets (Jupiter, Saturn, Uranus, Neptune) carrying the largest $v$-weights.

Reaching the full 3.62$\sigma$ requires all eight planets. No individual planet alone drives the result.

Bottom line

All 4 empirical tests reach $p < 0.01$ under the permutation null, and all 9 MC-combinable tests reach $p < 0.05$ under both Monte Carlo nulls. The direct joint combined $p$-value is robust across the three null distributions, comfortably exceeds the conventional $3\sigma$ “evidence” threshold, and lies just below the particle-physics $5\sigma$ “discovery” threshold ($p \approx 3 \times 10^{-7}$). Cited as a single figure, the model’s statistical signal is 3.62$\sigma$ under the most conservative null and 4.75$\sigma$–4.75$\sigma$ under the more general Monte Carlo nulls.

Exoplanet Context: A Question for Future Work

Compact multi-planet exoplanet systems raise the question of whether Fibonacci period-ratio patterns appear more broadly than the Solar System. The two best-characterized candidates are:

• TRAPPIST-1 (7 planets, TTV-measured): 5 of 6 adjacent period ratios sit near Fibonacci fractions (8:5, 5:3, 3:2, 3:2, 3:2), giving 83% Fibonacci content. Consistent with the Solar System’s period-ratio clustering.
• Kepler-90 (8 planets): 5 of 7 adjacent period ratios (71%) are near Fibonacci fractions, though only two planets have measured masses.

These observations are consistent with prior exoplanet findings that Fibonacci fractions cluster preferentially in orbital period distributions (Pletser 2019, Aschwanden 2018).

Why these are not used as confirmation of the Fibonacci Laws

Neither system can test the deeper Fibonacci Laws presented above:

• TRAPPIST-1 is a mean-motion resonance chain. Fibonacci-like period ratios (mainly 3:2) are partly a by-product of the resonance structure rather than an independent realization of the same mechanism. The system’s eccentricities span only ~0.002–0.01 — a factor-of-4 dynamic range compared with the Solar System’s 141$\times$ range — leaving Laws 3 and 5 (the inclination and eccentricity balance conditions) essentially untestable. Transit inclinations measure viewing geometry, not dynamical amplitudes, so Law 2 is inaccessible.
• Kepler-90 has only two TTV-measured masses and eccentricities, which precludes any mass-weighted ($\xi$, $\eta$) test.
• TTV mass and eccentricity uncertainties (typically 5–8% relative) are larger than any fine structural match that could be claimed from derived quantities.

The period-ratio observations in these systems are therefore suggestive rather than confirmatory. They motivate future work once more TTV-characterized exoplanet systems become available, but they are not used as support for any of the eleven significance tests above, which rest entirely on Solar System data.

Relation to Prior Work and What’s Novel

Existing Fibonacci research in planetary science

Molchanov (1968) — Integer resonances in orbital frequencies

Molchanov proposed that planetary orbital frequencies satisfy simultaneous linear equations with small integer coefficients. His framework used general small integers — not specifically Fibonacci numbers — and dealt exclusively with orbital frequencies. He did not incorporate planetary masses, eccentricities, or inclinations.

Key difference: Molchanov found approximate integer relations among frequencies. Our laws use specifically Fibonacci numbers, applied to mass-weighted eccentricities and inclinations.

Aschwanden (2018) — Harmonic resonances in orbital spacing

Aschwanden identified five dominant harmonic ratios governing planet and moon orbit spacing, achieving ~2.5% accuracy on semi-major axis predictions. The analysis is purely kinematic — planetary mass, eccentricity, and inclination do not appear.

Key difference: Aschwanden analyzed distance/period ratios between consecutive pairs. Our laws analyze mass-weighted values of individual planets.

Pletser (2019) — Fibonacci prevalence in period ratios

Pletser found ~60% alignment of period ratios with Fibonacci fractions, with the tendency increasing for minor planets with smaller eccentricities and inclinations. He used eccentricity and inclination only as selection filters.

Key difference: Pletser did not analyze whether eccentricities or inclinations themselves form Fibonacci ratios. Mass-weighting is entirely absent. None of our six laws appear in his analysis.

Comparison table

For a detailed discussion of what builds on established theory versus what is genuinely novel (9 items), see Relation to Existing Physics.

Assessment

The balance conditions (Laws 3 and 5) combine known conservation principles with a novel Fibonacci structure that modulates the planetary weights. The conservation laws guarantee that inclination and eccentricity oscillations balance around the invariable plane — but they do not predict that integer Fibonacci divisors should preserve that balance to such high precision. Laws 2 and 4 (the amplitude constants $\psi$ and $K$) are the most genuinely novel claims — no existing theory predicts that Fibonacci divisors should produce universal constants for either inclination or eccentricity amplitudes across all eight planets.

The key unresolved question is why Fibonacci numbers work: do they encode something about the secular eigenmode structure (real physics), or is the Fibonacci restriction a coincidence made possible by having enough number choices? The mirror symmetry and the triple-constraint uniqueness of Config #11 argue against coincidence, but a theoretical derivation from first principles — or a successful prediction for an independent system such as exoplanetary or satellite systems — would be needed to settle the question definitively.

Vector Balance and Eigenmode Analysis

The scalar balance laws (Law 3, Law 5) determine the unique d-value configuration. A separate question is whether the angular momentum perturbation vectors cancel at all times, not just on average.

The vector balance condition

Each planet’s orbital tilt creates a 2D angular momentum perturbation:

where $L_j = m_j\sqrt{a_j(1-e_j^2)}$ is the angular momentum, $i_j$ is the inclination, and $\Omega_j$ is the ascending node longitude. For the invariable plane to be stable, the sum of all 8 vectors must equal zero at every instant:

Single-mode vs multi-mode

Single-mode model: Each planet’s ascending node precesses at one constant rate (the model’s 8H/N period from ascendingNodeCyclesIn8H). Different planets have different rates, so the vectors rotate at different speeds. The cancellation geometry breaks over time, and the vector balance can drop to ~72% at some epochs.

Multi-mode model: Each planet’s ascending node position is reconstructed as the sum of 7 eigenmode oscillations (matching the secular eigenfrequencies $s_1, s_2, \ldots, s_8$, with $s_5 = 0$ excluded as the invariable plane itself). The eigenvector amplitudes are solved by least-squares from the J2000 state and rates, with an angular momentum constraint enforced per mode.