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.9974% 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.8636% with phase-derived base eccentricities)
6. Law 6 — Saturn-Jupiter-Earth Resonance: Earth’s $H/8$ obliquity cycle equals Jupiter’s ICRF perihelion and Saturn’s ecliptic perihelion at the Fibonacci-anchor level; on the 8H lattice the identity sits at $8H/65$ — 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 (J2000 anchor; the Fibonacci-divisor structure is invariant at any epoch, the literal year count rescales at geological time — see Expanding Resonance), 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 #4) among 15 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 that defines its inclination oscillation. The anchors are set so that at the balanced year (302,645 BC), all eight planets simultaneously reach minimum inclination. Saturn is anti-phase: its cosine sign is flipped in the formula, giving it the opposite contribution to the balance sum — but its inclination still lands at the minimum alongside the others at the balanced year. This synchronised minimum recurs 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$. The Fibonacci-divisor identities are structural and invariant at any epoch; the literal $H/F_n$ year counts in the table below are J2000 snapshots that rescale at geological time (see Expanding Resonance).

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 structural balance that locks Jupiter’s ICRF perihelion period to Saturn’s ecliptic perihelion period — Fibonacci anchor at $H/8$, 8H-lattice secular at $8H/65$, with Earth’s $H/8$ obliquity one lattice step away — 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 — two-channel decomposition

The eccentricity balance reaches 99.8636% but not exactly 100%. Sensitivity analysis on the 8-planet sum rules out mis-measurement of any single planet’s mass or orbit: the required shifts are 4–6 orders of magnitude larger than DE440/JPL precision. The residual decomposes into two physically-motivated channels:

(i) Minor-body contributions (~96% of the residual). Substituting Law 4 (e_amp = K · sin(tilt) · √d / (√m · a^(3/2))) into the Law-5 weight v = √m · a^(3/2) · e / √d produces a striking cancellation: every body following Law 4 contributes a uniform v = K · sin(tilt) ≈ 1.7 × 10⁻⁶ to the eccentricity balance, independent of mass and distance (the a^(3/2) and √m factors cancel exactly). Random ± in-phase/anti-phase aggregation across N ≈ 625 such bodies gives σ ≈ 4.3 × 10⁻⁵, matching the observed residual. The relevant population is sub-200 km low-eccentricity classical-belt KBOs, where the body’s actual oscillation amplitude approximates the Law-4 intrinsic prediction. Named TNOs (Pluto, Eris, Haumea, Sedna, …) are not part of this aggregation: their actual eccentricity amplitudes are dominated by external forcing (Neptune resonance, scattering, galactic tides), not by intrinsic Law-4 dynamics. For Pluto specifically, long-term numerical integrations (Williams & Benson 1971, Malhotra & Williams 1997) give actual e_amp ≈ 0.025, whereas Law 4’s intrinsic prediction is ~0.001 — a 1:24 intrinsic-to-external split, with Law 4 correctly capturing the intrinsic baseline while Neptune’s 3:2 resonance pumps the amplitude by ~25×.

(ii) Mass uncertainty (~4% of the residual). Uranus and Neptune masses are currently constrained only by Voyager 2 flybys (relative uncertainties ~5 × 10⁻⁴); the other six are known to 10⁻⁶ – 10⁻⁸ from orbiters and ranging. A 10⁵-trial Monte Carlo perturbation gives a combined σ ≈ 1.5 × 10⁻⁶ contribution to the Law-5 gap, dominated by Uranus.

Testable predictions. (a) The Vera Rubin Observatory (LSST, 2025–2035) is expected to constrain rotation poles for ~10³ small TNOs; a statistical clustering of small classical-belt KBO obliquities near arcsin(K · √d · a^(-3/2) · √m / e) would empirically support the bidirectional reading of Law 4 beyond the eight planets. (b) A future Uranus or Neptune orbiter (e.g., NASA’s proposed Uranus Orbiter and Probe, ~2032) would shrink the mass-uncertainty channel by ~100×, leaving the minor-body signature largely intact — discriminating the two channels in the budget.

Law-4 closure as a quantitative criterion for planethood. The pattern of which bodies satisfy Law-4 closure (actual e_amp ≈ intrinsic Law-4 prediction) and which are externally dominated maps onto the IAU’s 2006 third criterion: “has cleared the neighborhood around its orbit.” Bodies that have cleared their neighborhood evolve under their own dynamics — their secular e_amp is set intrinsically by Law 4. Bodies that have not (Pluto and the plutinos, scattered-disk objects, Sedna, comets, asteroids) have their secular evolution dominated by external coupling, exactly as the IAU criterion implies. The two formalisations partition the same eight bodies. Soter (2006) proposed a mass-discriminant Λ to formalise the IAU criterion quantitatively, but Λ never became canonical. Law-4 closure provides a different quantitative formalisation that returns a clean pass/fail and matches the same 8-body partition. The agreement is not tautological: all proposed Planet Nine candidates fail Law-4 closure by 4–7 orders of magnitude independent of how the observed eccentricity is interpreted (see Planet Nine prediction), making Law-4 closure a predictive criterion rather than a post-hoc relabelling.

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

Jupiter’s ICRF perihelion period and Saturn’s ecliptic perihelion period lock to a single period: $8H/65 = 41{,}270$ yr — a structural balance of the solar system, not a coincidence, and the obliquity beat recorded in Earth’s climate. Earth’s own obliquity sits one $8H$-lattice step away at the Fibonacci value $H/8$ (= $8H/64$): obliquity is Earth’s axial precession ($H/13$) beating against the ecliptic, so beating against the gas giants’ actual ecliptic period $8H/39$ gives $8H/65$, while beating against Law 1’s Fibonacci anchor $H/5$ gives $H/8$. The single-integer shift in the ecliptic period ($40 \to 39$) propagates one-for-one into the obliquity beat ($64 \to 65$). The gas giants drive Earth’s spin-axis dynamics through their mutual resonance lock. At the Fibonacci-anchor level all three quantities collapse onto $H/8$, closing through $3 + 5 = 8$.

Saturn’s perihelion precesses obviously retrograde in the ecliptic frame (period $= -8H/65 \approx -H/8$). JPL’s WebGeoCalc confirms this at ~-3,400 arcsec/century from SPICE ephemeris data, consistent with the model’s prediction of -3,425 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 §12 for the full observational evidence and comparison of the two theories.

The Fibonacci anchors that place Jupiter and Saturn in this resonance are Earth’s own precession periods — ecliptic precession ($H/5$) and obliquity ($H/8$), the first terms of Law 1’s hierarchy where $3 + 5 = 8$ and $5 + 8 = 13$. Jupiter’s perihelion falls near Earth’s $H/5$ and Saturn’s near Earth’s $H/8$, which is why the anchors are Fibonacci-clean — but the planets’ actual secular periods are $8H/39$ and $8H/65$, one lattice integer off on the 8H lattice. The clean $3 + 5 = 8$ triangle is therefore Earth’s internal precession hierarchy (Law 1), not a relation among the planets’ true periods.

Why Earth, Jupiter, and Saturn hold special roles

The genuine three-body relation is the 8H-lattice secular lock: Jupiter’s ICRF perihelion and Saturn’s ecliptic perihelion coincide at $8H/65$ — the period at which they drive Earth’s obliquity (the climate-recorded $k + s_3$ beat; Earth’s own obliquity cycle is $H/8 = 8H/64$, one step off). This is why these three planets recur throughout the Fibonacci Laws: they also dominate the inclination balance (Law 3) and the eccentricity balance (Law 5), three independent structural arguments converging on the same trio. Earth’s $d = 3$ (shared only with Saturn) feeds the universal inclination constant $\psi$ directly.

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 15 survivors, only one is mirror-symmetric: Config #4 (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 #4 of 15 by eccentricity balance (99.8636%).

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.8636%), 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}}$. Verification: at $\cos\theta = +1$ the formula reduces to $e_{\text{base}} - e_{\text{amp}}$ (minimum), at $\cos\theta = -1$ to $e_{\text{base}} + e_{\text{amp}}$ (maximum), and at $\cos\theta = 0$ to $\sqrt{e_{\text{base}}^2 + e_{\text{amp}}^2}$. The asymmetric $\cos\theta$ dependence preserves the Law 5 eccentricity balance at every epoch, not only at the extrema.

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 theory (Kolmogorov 1954, Arnold 1963, Moser 1962): orbits survive if their frequency ratio satisfies a Diophantine condition — a quantitative measure of irrationality — and the golden ratio $\varphi$ (whose convergents are the Fibonacci ratios $F_{n+1}/F_n$) is the hardest irrational to approximate by rationals (Hurwitz 1891). Greene (1979) showed the last invariant torus to break as perturbation increases has frequency ratio $\varphi$; Morbidelli & Giorgilli (1995) proved golden-ratio KAM tori have superexponentially long stability times.

Over 4.6 Gyr of evolution, orbits near low-order rational resonances are disrupted (Kirkwood gaps); orbits near golden-ratio KAM tori survive. The surviving population is organised around Fibonacci ratios. The AMD-natural variables $\xi = e\sqrt{m}$ and $\eta = i\sqrt{m}$ are the action-like variables in the secular Hamiltonian, and KAM stability constrains these to Fibonacci ratios — providing the theoretical foundation for why Fibonacci numbers specifically appear. Full conceptual treatment: Physical Origin; empirical evidence: Supporting Evidence §2.

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.9974%.
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.8636%.
6. Law 6 — E–J–S resonance — Do Jupiter’s ICRF perihelion and Saturn’s ecliptic perihelion coincide at the same period (Earth’s obliquity beat)? Observed: exact (both $8H/65$).
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

Compact multi-planet exoplanet systems raise the question of whether Fibonacci patterns appear more broadly than the Solar System. TRAPPIST-1 (7 planets, TTV-measured) shows 5 of 6 adjacent period ratios near Fibonacci fractions (8:5, 5:3, 3:2, 3:2, 3:2) — 83% Fibonacci content; Kepler-90 (8 planets) shows 5 of 7 (71%). Consistent with broader exoplanet findings (Pletser 2019, Aschwanden 2018).

Neither system can test the deeper Fibonacci Laws: TRAPPIST-1 is a mean-motion resonance chain (so Fibonacci-like period ratios are partly a by-product of the resonance structure), with eccentricities spanning only ~0.002–0.01 (vs the Solar System’s 141× range), so Laws 3 and 5 are essentially untestable, and transit inclinations measure viewing geometry rather than dynamical amplitudes (so Law 2 is inaccessible). Kepler-90 has only two TTV-measured masses and eccentricities. TTV uncertainties (typically 5–8%) are larger than any fine structural match. The period-ratio observations are therefore suggestive rather than confirmatory and 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 #4 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. Each planet’s orbital tilt creates a 2D perturbation $\vec{P}_j = L_j \sin(i_j) (\sin\Omega_j, \cos\Omega_j)$ with $L_j = m_j\sqrt{a_j(1-e_j^2)}$; for the invariable plane to be stable, $\sum_{j=1}^{8} \vec{P}_j(t) = 0$ at every instant.

In the model’s single-mode approximation (one constant 8H/N rate per planet), the vectors rotate at different speeds, and the cancellation geometry breaks over time — vector balance can drop to ~72% at some epochs. In a multi-mode reconstruction (7 eigenmodes per planet matching $s_1, s_2, \ldots, s_8$ excluding $s_5 = 0$), vector balance is maintained at 100% — but this is a mathematical property of the solver (56 free parameters, only 23 constraints; 33 degrees of freedom), not a unique validation of any specific frequency set. The invariable plane is stable by definition; any eigenmode decomposition that reproduces the J2000 state will automatically maintain that.

Full ascending-node periods (8H/N integer table, JPL J2000-fixed-frame fits, Jupiter–Saturn lockstep at N=36, cumulative residual ~5.8″/cy) are canonical at Invariable Plane §Ascending node periods. The model’s structural advantage there: all 7 periods derive from a single constant $H$, while Laskar’s are 7 independent measurements with no known relationship to each other.

Open Questions

For the full list of open questions, see What Remains Unknown.

References

Fibonacci and orbital resonance

1. Molchanov, A.M. (1968). “The resonant structure of the Solar System.” Icarus, 8(1-3), 203-215. ScienceDirect 
2. Backus, G.E. (1969). “Critique of ‘The Resonant Structure of the Solar System’ by A.M. Molchanov.” Icarus, 11, 88-92.
3. Molchanov, A.M. (1969). “Resonances in complex systems: A reply to critiques.” Icarus, 11(1), 95-103.
4. Aschwanden, M.J. (2018). “Self-organizing systems in planetary physics: Harmonic resonances of planet and moon orbits.” New Astronomy, 58, 107-123. arXiv 
5. Aschwanden, M.J., & Scholkmann, F. (2017). “Exoplanet Predictions Based on Harmonic Orbit Resonances.” Galaxies, 5(4), 56. MDPI 
6. Pletser, V. (2019). “Prevalence of Fibonacci numbers in orbital period ratios in solar planetary and satellite systems and in exoplanetary systems.” Astrophysics and Space Science, 364, 158. arXiv 

KAM theory and orbital stability

7. Kolmogorov, A.N. (1954). “On the conservation of conditionally periodic motions under small perturbation of the Hamiltonian.” Doklady Akad. Nauk SSSR, 98, 527-530.
8. Arnold, V.I. (1963). “Proof of a theorem of A.N. Kolmogorov on the invariance of quasi-periodic motions under small perturbations of the Hamiltonian.” Russian Math. Surveys, 18(5), 9-36.
9. Greene, J.M. (1979). “A method for determining a stochastic transition.” J. Math. Phys., 20, 1183-1201.
10. Morbidelli, A., & Giorgilli, A. (1995). “Superexponential stability of KAM tori.” J. Stat. Phys., 78, 1607-1617.
11. Laskar, J. (1989). “A numerical experiment on the chaotic behaviour of the Solar System.” Nature, 338, 237-238.
12. Laskar, J. (1997). “Large scale chaos and the spacing of inner planets.” Astronomy & Astrophysics, 317, L75-L78. — Introduces the Angular Momentum Deficit (AMD) framework underlying the model’s mass-weighted variables.
13. Celletti, A., & Chierchia, L. (2007). “KAM stability and celestial mechanics.” Memoirs of the AMS, 187.