Fibonacci Laws of Planetary Motion

The orbits of the eight planets are not random. Their shapes (eccentricities) and tilts (inclinations) follow patterns built on Fibonacci numbers — the sequence 1, 1, 2, 3, 5, 8, 13, 21, 34, … where each number is the sum of the two before it. Fibonacci ratios converge toward the golden ratio (1.618…), which plays a central role in the stability of dynamical systems. In the solar system, Fibonacci patterns mark the orbits that survived 4.5 billion years of gravitational evolution.

Six independent laws connect planetary precession periods, inclinations, and eccentricities through Fibonacci numbers and mass-weighted quantities, predicting orbital properties for all eight planets from a single timescale: the Earth Fundamental Cycle H = 335,317 years. This value is empirically fitted — it is the unique timescale that simultaneously satisfies eight independent constraints from observed precession cycles, climate patterns, and integer day counts. The six Fibonacci-aligned laws operate within the broader Solar System Resonance Cycle (8H), where every planet’s precession periods divide 8H by integers — see Fundamental Cycles for the all-planet scaffold within which the Fibonacci patterns sit. See Mathematical Foundation for the derivation of H itself.

The Scientific Foundation

The Fibonacci Laws rest on three pillars — each independently supported by established science and peer-reviewed research.

1. Two real counter-rotating motions

The model’s foundation is not speculative. Earth undergoes two well-documented precession motions that move in opposite directions:

Both are standard astronomy with known physical causes — documented in IAU reference models (Capitaine et al. 2003 ) and planetary ephemerides. The current periods vary slightly over time; the mean periods over the full 335,317-year Earth Fundamental Cycle are the structurally significant values. The opposite directions are a consequence of the underlying physics: axial precession follows from gyroscopic torque on Earth’s equatorial bulge, while apsidal precession follows from planets pulling Earth’s perihelion forward in its orbital direction. The key observation: the apsidal-to-axial ratio of their mean periods is ~111,772 : ~25,794 = 13 : 3 — a ratio of two Fibonacci numbers.

2. Fibonacci patterns in orbital mechanics are documented

Fibonacci-related frequency ratios in the solar system are not unique to this model. They appear throughout the peer-reviewed literature:

• Molchanov (1968) in Icarus found that many planetary orbital resonances approximate Fibonacci fractions, including Venus/Earth at 8/13 = 0.615
• Aschwanden (2018) analyzed 75 orbital period ratios in the solar system (planets and moons) and found approximately 60% match Fibonacci fractions within measurement uncertainty
• Pletser (2019) in Astrophysics and Space Science confirmed that ~60% of planetary and satellite period ratios preferentially cluster near Fibonacci fractions, with these orbits associated with more regular, less inclined, and more circular configurations
• Aschwanden & Scholkmann (2017) found Fibonacci harmonic ratios in 73% of 932 exoplanet pairs — extending the pattern well beyond our solar system
• Fibonacci patterns also appear in DNA structure (34 × 21 angstroms per turn; Yamagishi & Shimabukuro 2008 ), phyllotaxis (sunflower spirals, pineapple scales), and numerous other natural systems

The Kirkwood Gaps in the asteroid belt provide dramatic visual confirmation: asteroids at simple integer resonances with Jupiter (3:1, 5:2, 2:1) have been swept away, while those at Fibonacci-related resonances (like the Hilda group at 3:2) form stable clusters.

3. KAM theory explains why

The Kolmogorov–Arnold–Moser (KAM) theorem (1954–1963) provides the rigorous mathematical explanation. In perturbed dynamical systems, orbits whose frequencies have the “most irrational” ratios are maximally stable against perturbation. The golden ratio φ ≈ 1.618 — toward which Fibonacci ratios converge — is the most irrational number in a precise mathematical sense: it is the hardest to approximate by simple fractions. This means orbits with golden-ratio-related frequencies are the last to become unstable when perturbed by gravitational interactions.

Over 4.5 billion years, the solar system has naturally selected for these maximally stable configurations. Greene (1979) confirmed computationally that the “golden” invariant torus is indeed the last to break. Morbidelli and Giorgilli (1995) showed super-exponential stability near golden-ratio frequency ratios in the asteroid belt.

1. The Six Fibonacci Laws

The six laws form a symmetric architecture anchored on Earth and extending to all eight planets. Law 1 establishes the Earth-specific timescale — Earth’s major precession periods are H divided by Fibonacci numbers, a hierarchy unique to Earth among the planets. Laws 2–5 extend the framework to all eight planets via their individual Fibonacci d-values, combining each planet’s orbital property with the square root of its mass — a weighting dictated by the Angular Momentum Deficit (see section 3 below). Laws 2 and 4 quantify amplitudes through universal constants (ψ for inclination, K for eccentricity); Laws 3 and 5 establish balance conditions where seven planets cancel against Saturn alone. Law 6 identifies the Earth–Jupiter–Saturn three-planet coupling that ties the gas-giant resonance to Earth’s spin-axis dynamics.

Law 1: The Fibonacci Cycle Hierarchy

Earth’s major precession periods divide the Earth Fundamental Cycle H 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 Earth Fundamental Cycle H = 335,317 years is the master timescale. When divided by the Fibonacci numbers 3, 5, 8, 13, …, the resulting periods correspond to observed astronomical cycles:

These periods are not independent — they follow the same addition rule as the Fibonacci sequence itself. Just as 3 + 5 = 8, the corresponding frequencies add:

1/~111,772 + 1/~67,063 = 1/~41,915

This works at every level: 1/~67,063 + 1/~41,915 = 1/~25,794, and 1/~41,915 + 1/~25,794 = 1/15,967. The beat frequency rule is a direct consequence of the Fibonacci addition, applied to the timescale hierarchy. The ratios between consecutive periods approach the golden ratio: ~111,772/~67,063 = 1.667, ~67,063/~41,915 = 1.600, converging toward 1.618…

Earth’s H/Fibonacci hierarchy is unique. Only Earth has all major precession periods at H divided by Fibonacci numbers. Jupiter’s and Saturn’s perihelion periods coincide with some of these Fibonacci values (H/5, H/8, H/21) — but those specific coincidences are the subject of Law 6. The other planets’ precession periods divide the Solar System Resonance Cycle (8H) by various integers, mostly non-Fibonacci. See Fundamental Cycles for the all-planet 8H/N table.

Law 2: The Inclination Constant

Each planet’s mass-weighted inclination amplitude, multiplied by a Fibonacci divisor, equals the same universal constant.

Inclination measures how tilted an orbit is relative to the average plane of the solar system. Each planet’s tilt oscillates slowly over hundreds of thousands of years. The mass-weighted inclination amplitude is η = amplitude × √m. When multiplied by a Fibonacci divisor d specific to each planet, the result is a single constant ψ shared by all eight planets:

d × η = ψ

where d is a pure Fibonacci number — the same sequence (3, 5, 21, 34) that appears throughout the model. The constant ψ is empirical, derived from Earth’s fitted inclination amplitude:

ψ = d_E × amp_E × √m_E = 3.3068 × 10⁻³

Why the symbol ψ? A notation choice — ψ² appears in the energy formula AMD = √a × ψ² / (2d²), 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. The resemblance is purely formal; the model makes no claim about quantum physics.

Worked example — Earth’s inclination amplitude:

Rearranging the law: amplitude = ψ / (d × √m), where mass is in solar masses (M☉ = 1.989 × 10³⁰ kg):

Earth’s orbit tilts ±0.63603° around its mean inclination of 1.48113° — oscillating between 0.845° and 2.117° over hundreds of thousands of years. The same formula, with only the divisor d and mass changed, predicts every other planet’s tilt range.

The eight divisors form a mirror-symmetric pattern across the asteroid belt (see section 4).

What ψ physically means — mass-independent energy partition: Law 2 divides each planet’s tilt by √m — heavier planets tilt less. But this mass dependence is not arbitrary: √m is the unique exponent that makes the underlying energy partition completely mass-free. The standard AMD for a planet’s inclination oscillation is m × √a × amp² / 2. Substituting amp = ψ / (d × √m) from Law 2, the m in front cancels the m in the denominator:

m × √a × [ψ / (d√m)]² / 2 = √a × ψ² / (2d²)

Each planet’s share of the eight-planet inclination oscillation energy is proportional to √a/d²:

Saturn alone carries over half the tilt oscillation energy — not because it is the most massive (Jupiter is 3× heavier), but because it has the lowest d combined with a large orbit. Earth, despite being 1000× lighter than Jupiter, carries more oscillation energy (18.2% vs 14.9%) because its d = 3 beats Jupiter’s d = 5 in the 1/d² 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%. (Shares are relative to the eight major planets, which carry 99.994% of the system’s orbital angular momentum; TNOs contribute the remainder.)

Connection to Law 3: The same structural properties that make Saturn the dominant energy carrier — d = 3 and a large orbit — also make it the sole counterweight in the inclination balance (Law 3 below). In the mass-independent energy partition above, Saturn carries 56.1%. In the mass-dependent balance of Law 3, Saturn carries exactly 50% (by definition — it equals the other seven combined). The gap is absorbed by mass: Jupiter is 3.3× heavier than Saturn, which boosts Jupiter’s balance weight relative to its energy share, redistributing the load to an exact 50/50 split.

This also reveals what determines ψ’s magnitude: it is fixed by the total inclination oscillation energy of the solar system (set at formation) and the geometric sum Σ(√a/d²) (set by Fibonacci structure and orbital spacing). Law 3 below tells us the oscillations cancel; the budget equation tells us their scale. See Fibonacci Laws Derivation — Physical meaning of ψ for the full derivation.

Law 3: The Inclination Balance

The angular-momentum-weighted inclination oscillations of seven planets balance against Saturn’s alone.

The invariable plane — the fundamental reference plane of the solar system, perpendicular to the total angular momentum vector — must remain stable. Each planet’s orbital tilt oscillates around a per-planet cycle anchor (the ICRF perihelion longitude where MAX inclination occurs, evaluated at the balanced year). For the invariable plane to stay fixed, the angular-momentum-weighted oscillations of the two groups must cancel:

In-phase group: Mercury, Venus, Earth, Mars, Jupiter, Uranus, Neptune — Anti-phase group: Saturn (alone)

How the balance works: Each planet’s structural weight is w = √(m × a) / d, where a is the semi-major axis. Because every planet shares the same ψ (Law 2), the universal constant cancels — the balance depends only on these structural weights. The sum of the seven in-phase weights must equal Saturn’s weight alone:

The balance is 99.9975% (with phase-derived base eccentricities, which enter via the angular momentum factor √(1 − e²)). Jupiter alone gets 81% of the way to matching Saturn; Uranus and Neptune contribute most of the remaining 19%. The four inner planets barely matter — their small masses and distances make their weights negligible.

Law 4: The Eccentricity Amplitude Constant

A single constant K predicts all eight eccentricity amplitudes from Fibonacci divisors, mass, distance, and axial tilt.

Each planet’s eccentricity oscillates around a base value over its eccentricity cycle. The oscillation amplitude follows a universal formula:

e_amp = K × sin(tilt) × √d / (√m × a^1.5)

where tilt is the planet’s axial tilt (obliquity) — the sine takes the angle in either degrees or radians, since astronomical software handles the conversion automatically. K = 3.4149 × 10⁻⁶, derived from Earth’s eccentricity amplitude and axial tilt. This is the eccentricity analog of ψ (Law 2) for inclination amplitudes — a single constant, zero free parameters, all eight planets predicted.

Both constants are empirical — derived from Earth — and predict all 8 planets with zero free parameters. K additionally uses semi-major axis and axial tilt, coupling the spin and orbital domains.

Law 5: The 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.

Using the same divisors and the same Saturn-vs-seven grouping, each planet receives an eccentricity weight v = √m × a³ᐟ² × e / √d, where a is the semi-major axis and e is the eccentricity. Compared to the inclination weight (w = √(m × a) / d), the eccentricity weight scales differently: a³ᐟ² instead of √a, linear e instead of no eccentricity term, and 1/√d instead of 1/d. The sum of the seven in-phase weights must equal Saturn’s weight alone:

Balance: 99.8632% (with phase-derived base eccentricities)

The load distribution differs from Law 3: Jupiter carries only 51% (vs 81% for inclination), while Uranus contributes 37% — making the eccentricity balance more evenly spread among the three outer in-phase planets.

This balance is not a structural artifact — three tests confirm it depends on the actual eccentricity values:

1. Without eccentricities, the weight formula gives only 74% balance — the observed eccentricities contribute ~26 percentage points of improvement
2. Random eccentricities substituted into the same formula give only 50–85% balance
3. The balance peaks sharply at linear eccentricity (99.8632% with phase-derived base eccentricities), dropping to 91% for eccentricity squared — linear dependence is special

This is genuinely independent of Law 3: the eccentricity weights differ from the inclination weights by over 100-fold for some planets, yet both conditions are satisfied by the same set of Fibonacci divisors and the same phase assignment. The eccentricity balance operates on linear e rather than e², distinguishing it from the quadratic Angular Momentum Deficit (AMD) — the standard conserved quantity in secular theory. What conservation law produces this linear balance remains an open question.

Because Saturn is the only planet on one side, the balance equation directly predicts Saturn’s eccentricity from the other seven: e_Saturn = (in-phase total) / (√m_Sa × a_Sa³ᐟ² / √d_Sa) = 0.01543 / 0.2872 = 0.05372 (observed J2000: 0.05386, error −0.27%). This phase-derived base value is the model’s prediction.

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.

As seen from Earth’s ecliptic frame, Saturn’s perihelion precesses retrograde (opposite to orbital motion) with a period of ~41,915 years = H/8. JPL’s WebGeoCalc confirms the ecliptic-retrograde motion at ~-3,372 arcsec/century — see Supporting Evidence §13 for how this observation distinguishes the model from the standard Great Inequality explanation. When Saturn’s retrograde motion interacts with Jupiter’s prograde motion, it creates beat frequencies that form a closed loop:

All three rows are cyclic permutations of a single Fibonacci identity: 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’s higher-level identity: combining Jupiter’s H/5 and Saturn’s H/8 contributions produces Earth’s axial precession: 1/~67,063 + 1/~41,915 = 1/~25,794 (i.e., 5 + 8 = 13), extending the triangle into the full cycle hierarchy.

The ICRF Perspective: A Fibonacci Chain

The ecliptic frame — Earth’s orbital plane — is the natural frame for solar system dynamics. Secular perturbation theory, the Laplace-Lagrange eigensystem, and angular momentum conservation all operate in this plane; the planetary precession rates emerge from mutual gravitational interactions within it. The ICRF (International Celestial Reference Frame), fixed to distant quasars, provides an inertial reference — but no physical mechanism couples the solar system’s internal precession to that external frame. The ICRF rates are therefore a kinematic consequence of the ecliptic dynamics, not an independent constraint.

So what happens when we view the same motions from the ICRF? The conversion subtracts H/13 (the general precession) — the rate at which the ecliptic reference direction itself rotates in the fixed frame. This correction is the same for all planets: it is a property of the reference frame, not the individual orbit.

In ICRF, only Earth precesses prograde. Earth’s ecliptic rate (H/16 = 61.840″/yr) is the only one that exceeds the general precession threshold. All other planets — including Jupiter — precess retrograde in the ICRF. For the three Law 6 planets:

The frame transformation produces Fibonacci subtraction identities:

• Earth: ecliptic 16 − general precession 13 = ICRF +3  →  16 − 13 = 3
• Jupiter: ecliptic 5 − general precession 13 = ICRF −8  →  5 − 13 = −8
• Saturn: ecliptic −8 − general precession 13 = ICRF −21  →  −8 − 13 = −21

These follow directly from the Fibonacci recurrence: since 5 + 8 = 13, subtracting 13 from 5 gives −8; since 8 + 13 = 21, subtracting 13 from −8 gives −21. The complete Fibonacci chain 3 → 5 → 8 → 13 → 21 is generated by a single operation: subtracting the general precession (13) from the ecliptic denominators. Since 13 is itself a Fibonacci number, the results are Fibonacci.

Earth is the sole prograde planet in the ICRF — an exception created by the Fibonacci number 13 (the general precession). Only the three Law 6 planets (Earth, Jupiter, Saturn) have pure Fibonacci ICRF denominators (3, 8, 21); other planets have non-Fibonacci ICRF denominators, consistent with Jupiter and Saturn carrying ~85% of the solar system’s angular momentum. For the full picture of Earth’s unique structural role — including the reference frame duality and the Fibonacci addition chain — see Why Earth Is Special.

2. What the Laws Predict

Precession hierarchy — from a single timescale

Law 1 predicts all major precession periods from the Earth Fundamental Cycle H = 335,317. Six Fibonacci divisions (H/3 through H/34) match observed astronomical cycles, and their beat frequency relationships are exact algebraic identities. Law 6 shows these are physically coupled through the Saturn-Jupiter-Earth resonance loop.

Inclination amplitudes — zero free parameters

The most striking result is the zero-parameter prediction of inclination amplitudes. Given only the Earth Fundamental Cycle (H = 335,317), planetary masses, and Fibonacci divisors, the model predicts how far each planet tilts during its oscillation cycle:

All eight predictions fall within the bounds of Laplace-Lagrange secular theory — the standard framework of celestial mechanics for computing long-term orbital evolution. This is the non-trivial test: the Fibonacci formula predicts every planet’s tilt range, and secular theory independently confirms each is physically allowed.

The divisor assignments (3, 5, 21, 34) are not tunable — they are specific Fibonacci numbers found by exhaustive search as the unique mirror-symmetric configuration surviving all five physical filters (see section 4).

Saturn eccentricity — predicted from Law 5

Law 5 (eccentricity balance) is one equation in eight unknowns. Since Saturn is the only anti-phase planet, the equation directly determines Saturn’s eccentricity from the other seven:

The d-values (from Laws 1, 2, and 3) were chosen to satisfy the inclination constraints — not tuned for eccentricity. Yet Law 5 predicts Saturn’s base eccentricity to within 0.27% of the observed value. This is a non-trivial cross-validation of the d-value choice.

Law 4 predicts all eight eccentricity amplitudes. Law 5 predicts Saturn’s base eccentricity. The remaining seven base eccentricities are derived from the System Reset phase (see What Remains Unknown below).

3. Background to Laws 2-5

Laws 2–5 quantify individual and collective orbital properties — inclination amplitudes and their balance (Laws 2–3), eccentricity amplitudes and their balance (Laws 4–5). All four laws multiply orbital properties by √m (the square root of mass). This section explains why that weighting is not a free choice, and introduces the master ratio that connects the inclination scale to the eccentricity scale.

Why the square root of mass? (Laws 2–5)

The mass-weighting in Laws 2–5 is dictated by physics. In celestial mechanics, the Angular Momentum Deficit (AMD) is the key conserved quantity governing long-term orbital stability. It measures how much an orbit deviates from a perfect circle in the reference plane.

AMD naturally splits into two independent parts:

• Eccentricity contribution: proportional to eccentricity times the square root of mass
• Inclination contribution: proportional to inclination times the square root of mass

The square root of mass is the unique exponent that makes these contributions meaningful for stability analysis. Testing all possible mass exponents from 0 to 1 confirms this: at exactly 0.50 (square root), the Fibonacci structure achieves a spread of just 0.11% across all planets. The next-best exponent gives a spread above 28% — over 250 times worse.

This connects the Fibonacci Laws to 200 years of celestial mechanics, from Lagrange and Laplace’s original secular perturbation theory (1780s) through to Laskar’s modern AMD stability criteria (1997). The Fibonacci patterns are not imposed on the orbits — they emerge naturally when orbits are described in the variables that physics itself selects.

4. Mirror Symmetry and Configuration Uniqueness

The eight Fibonacci divisors are not random — they form an exact mirror-symmetric pattern across the asteroid belt:

Each inner planet shares its Fibonacci divisor with an outer counterpart. The divisors form two consecutive Fibonacci pairs: (3, 5) for the belt-adjacent planets and (21, 34) for the outermost pairs. Earth and Saturn share divisor 3, but are the only pair with opposite balance groups (in-phase vs anti-phase) — Saturn’s unique anti-phase behavior (MAX inclination at balanced year, while all others are at MIN) is what makes it the balancing “pivot” for both Law 3 and Law 5.

This symmetry was not assumed — it was discovered

An exhaustive search tested 7,558,272 possible planet configurations — assignments of Fibonacci divisors and phase groups for 5 free planets (Earth locked at d=3; Jupiter and Saturn constrained by the resonance structure). These were filtered through a pipeline of successively stricter physical constraints:

Why “per-config optimised”? The Laplace–Lagrange bounds and direction-trend checks depend on two parameters that vary per configuration: the anchor position n (where in the 8H Solar System Resonance Cycle the balanced year falls — this sets all cycle anchors) and the ascending node integers N per planet (the ascending node regression period = −8H/N). To make the comparison fair, each of the 96 candidates is evaluated at its own optimal (n, N), not at a reference configuration’s tuning. Jupiter and Saturn are constrained to share the same N, preserving their lockstep regression over the Solar System Resonance Cycle.

Of the 42 surviving configurations, only one is mirror-symmetric: Config #11 (Me=21, Ve=34, Ma=5, Ju=5, Sa=3, Ur=21, Ne=34). It ranks #11 of 42 by eccentricity balance (99.8632%). Seven non-mirror-symmetric configurations achieve slightly higher eccentricity balance, but none share the inner–outer Fibonacci pairing across the asteroid belt.

The mirror symmetry is not an input to the model — it is its most surprising output. The asteroid belt serves as a natural mirror axis, separating the solar system into two halves whose Fibonacci structures reflect each other. All 42 configurations are available for comparison in the interactive Balance Explorer .

5. Why Fibonacci? The KAM Connection

The appearance of Fibonacci numbers in orbital mechanics is not numerology — it is a mathematical consequence of long-term stability, explained by one of the deepest results in dynamical systems theory.

The KAM theorem

In the 1950s and 1960s, three mathematicians — Kolmogorov (1954), Arnold (1963), and Moser (1962) — proved a remarkable theorem about perturbed dynamical systems. In simplified terms:

Orbits whose frequencies have “most irrational” ratios are the most resistant to perturbation.

What makes a number “most irrational”? The golden ratio (1.618…) holds this distinction in a precise mathematical sense — it is the number hardest to approximate by simple fractions. Since Fibonacci ratios (3/2, 5/3, 8/5, 13/8, …) converge to the golden ratio, orbits with Fibonacci-related frequencies are the last to become unstable when perturbed by other planets.

Natural selection of orbits

The solar system is 4.5 billion years old. Over this vast timescale:

1. The early solar system contained many more objects on various orbits
2. Orbits with simple-fraction frequency ratios (like 2:1 or 3:1) experienced resonances — repeated gravitational kicks that destabilized them
3. Orbits with golden-ratio-related frequencies avoided these resonances
4. What survives today are the maximally stable configurations — those organized by Fibonacci numbers

Greene (1979) showed computationally that the “golden” invariant torus — the orbit with the most irrational frequency ratio — is indeed the last to break under perturbation. This has been confirmed in numerous studies of Hamiltonian systems.

From theory to observation

The KAM theorem predicts that stable orbits should have frequencies near (not exactly at) golden-ratio relationships. This is precisely what the Fibonacci Laws show: the secular eigenfrequencies of the planets are close to Fibonacci ratios, but not exactly on them — as expected for a real, slightly perturbed system. The Kirkwood gaps in the asteroid belt provide dramatic visual confirmation: asteroids at simple resonances with Jupiter (3:1, 5:2, 2:1) have been swept away, while those between the gaps survive.

6. Relation to Existing Physics

The Fibonacci Laws combine established physics with genuinely new predictions. Understanding the boundary is important.

What builds on established theory

Law 1 (Fibonacci Cycle Hierarchy) builds on the known precession periods of the planets. The individual periods are established results of celestial mechanics — what is new is recognizing them as a single Fibonacci-divided hierarchy.

Law 3 (Inclination Balance) is rooted in angular momentum conservation. The invariable plane is defined as the plane perpendicular to the total angular momentum vector — inclination oscillations must balance around it. The novel contribution is that dividing each planet’s angular momentum by a Fibonacci divisor preserves this balance to 99.9975% (with phase-derived base eccentricities entering via the angular momentum factor √(1 − e²)).

Phase angles are per-planet values (ICRF perihelion longitude at the balanced year). Saturn’s anti-phase behavior (cosine sign flipped relative to all other planets) is consistent with its known retrograde precession in secular theory. The full anti-phase alignment occurs once per Solar System Resonance Cycle (8H).

Law 5 (Eccentricity Balance) connects to Angular Momentum Deficit (AMD) conservation, a known conserved quantity. However, the linear dependence on eccentricity (rather than quadratic, as in AMD) and the 1/√d scaling distinguish it from standard formulations.

Law 6 (Saturn-Jupiter-Earth Resonance) describes beat frequency relationships between known precession periods. The individual periods are well-established — the closed-loop structure connecting them through Fibonacci addition is the novel observation.

What appears genuinely new

1. Fibonacci division of a single timescale — No known physical theory predicts that all major precession periods should equal H/F(n) for a single H and successive Fibonacci numbers F(n). The existence of a master timescale H is without precedent.

2. Fibonacci structure in eccentricities and inclinations — Previous work (Molchanov, Pletser) found Fibonacci patterns only in orbital periods. The Fibonacci Laws extend this to two additional orbital elements — eccentricity and inclination — through the mass-weighted variables ξ = e × √m and η = i × √m.

3. Fibonacci quantization of oscillation amplitudes — No known physical theory predicts that Fibonacci divisors should produce universal constants for either inclination amplitudes (Law 2, constant ψ) or eccentricity amplitudes (Law 4, constant K). Both are empirical, derived from Earth, and predict all 8 planets.

4. Two universal amplitude constants ψ and K — Both connect Fibonacci divisors to oscillation amplitudes across all eight planets with zero free parameters. No theoretical framework explains why they exist.

5. Mirror symmetry across the asteroid belt — No known law predicts that inner and outer planets should pair with identical Fibonacci divisors in distance order.

6. Simultaneous satisfaction of six independent constraints — Pure Fibonacci divisors satisfy Laplace-Lagrange bounds, inclination balance, eccentricity balance, eccentricity-inclination partition, timescale hierarchy, and resonance closure at the same time. Out of 7,558,272 configurations tested, only 42 survive all five physical filters — and only one of those is mirror-symmetric.

7. Saturn eccentricity prediction from Law 5 — Law 5 (eccentricity balance) is one equation in eight unknowns and uniquely determines Saturn’s base eccentricity from the other seven, to ~0.27% of the J2000 observed value. The d-values were chosen to satisfy Laws 1–3 (inclination constraints) and were never optimized for eccentricity, so this is a non-trivial cross-validation.

8. Closed resonance loop as Fibonacci identity — The Saturn-Jupiter-Earth beat frequency triangle is an algebraic consequence of 3 + 5 = 8, but no theory predicts that three planets’ precession denominators should satisfy a Fibonacci addition rule.

7. Statistical Significance

Are these patterns real, or just numerology? In 1969, Backus criticized Molchanov’s earlier work on planetary resonances for failing to prove statistical significance. The Fibonacci Laws address this head-on with a comprehensive analysis: 11 test statistics covering all 6 laws, each evaluated against three different null models — random planetary systems generated by permutation (reshuffling the real values among planets) and two types of Monte Carlo simulation (100,000 fully random systems each).

Of the 11 tests, only 4 are empirically testable: Law 3 (inclination balance), Law 5 (eccentricity balance), Finding 4 (Saturn eccentricity prediction), and Finding 6 (solo planet identification). The remaining 7 are excluded from the combined statistic for one of two reasons: they are either multiset-invariant under permutation (Laws 1 and 6, Findings 1 and 1b, the year-length beat — permutation cannot disturb them by construction) or tautological (Laws 2 and 4, where the model defines amplitudes as $\psi/(d\sqrt m)$ and $K$-derived, so the test relations hold identically). See the warning callout in the technical derivation for why no independent test of Laws 2 and 4 is possible — multi-millennial planetary oscillation amplitudes are not directly observable on human timescales, and external secular-theory amplitudes are themselves model outputs, not measurements.

The four empirical tests:

• Inclination balance (Law 3) — 99.9975% balance, p = 0.0074 (permutation), p < 10⁻⁴ (Monte Carlo)
• Eccentricity balance (Law 5) — 99.8632% balance, p = 0.00072 (permutation), p ≈ 0.0004 – 0.0005 (Monte Carlo)
• Saturn eccentricity prediction (Finding 4) — 0.27% error, p = 0.00072 (permutation), p ≈ 0.0004 – 0.0005 (Monte Carlo)
• Solo planet identification (Finding 6) — 0.27% min residual, p = 0.0040 (permutation), p ≈ 0.003 (Monte Carlo)

These four tests are not statistically independent — all four use the same fundamental quantity v_j = √m × a^(3/2) × e / √d — so their per-test p-values cannot simply be multiplied. The headline statistic is a direct joint permutation test: 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 = Σ z_i, and the p-value is the fraction of null permutations with T_null ≥ T_obs. This approach is model-independent (no distributional assumption), self-correcting for correlation (the joint null captures it directly — no variance inflation factor needed), and has no floor-clamp artifact (p ≥ 1/(n+1) by construction).

Combined p = 1.5 × 10⁻⁴ ≈ 3.6σ (direct joint permutation test)

For reference, the measured average pairwise correlation between the four tests under the permutation null is r̄ = 0.283 — smaller than the r̄ ≈ 0.5 previously assumed in older Stouffer+Brown approximations. The direct joint test does not need to correct for correlation at all because the correlation is already built into the joint null distribution.

The strongest individual results — the balance conditions and Saturn prediction — have zero look-elsewhere effect: the Fibonacci divisors, phase groups, and target planet are fixed by the model, not optimized after the fact.

Putting this in perspective: The result is robust across the three null distributions, and the MC nulls (which additionally include Laws 1, 6; F1, F1b; and the year-length beat — permutation-invariant tests that become meaningful under random solar systems) give an even stronger answer:

• 1.5 × 10⁻⁴ (≈3.62σ) — permutation null over 4 empirical tests. Most conservative.
• 1.0 × 10⁻⁶ (≈4.75σ) — log-uniform Monte Carlo over 9 tests.
• 1.0 × 10⁻⁶ (≈4.75σ) — uniform Monte Carlo over 9 tests.

All three exceed the conventional 3σ “evidence” threshold and the Monte Carlo results lie just below the particle-physics 5σ “discovery” threshold (p ≈ 3 × 10⁻⁷). The Monte Carlo results — which test against fully randomized planetary systems with random eccentricities, inclinations, and d-values — are arguably the more relevant question and give the stronger answer (~4.75σ). The permutation result (3.62σ) is the more defensive figure to cite, but the convergence of all three nulls in the 3.62–4.75σ range makes the underlying conclusion robust to the methodological choice.

Jackknife robustness. A leave-one-out permutation test (5,040 trials per planet drop) reveals that no single planet drives the result: dropping any of Mercury, Venus, Earth, Mars, Saturn, or Neptune leaves a signal of 1.5–2.8σ. However, dropping Jupiter (→ 0.5σ) or Uranus (→ 0.9σ) collapses the combined statistic — both planets carry large v-weights and their absence destroys the balance. The full 3.62σ requires all eight planets. See the jackknife table in the technical derivation.

8. 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. Two best-characterized candidates:

• TRAPPIST-1: 5 of 6 adjacent period ratios near Fibonacci fractions (83%)
• Kepler-90: 5 of 7 adjacent period ratios near Fibonacci fractions (71%)

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

However, neither system can test the Fibonacci Laws presented in this document. TRAPPIST-1 is a mean-motion resonance chain, so Fibonacci-like period ratios are partly a by-product of resonance rather than an independent realization of the same mechanism. Its eccentricities span only ~0.002–0.01 (a factor-of-4 dynamic range compared with the solar system’s 141×), leaving Laws 3 and 5 essentially untestable. Kepler-90 has only two mass measurements. TTV mass and eccentricity uncertainties (5–8% relative) are larger than any fine structural match that could be claimed.

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

9. Three Fibonacci Levels

The Fibonacci structure operates at three distinct levels — each independently verifiable:

Level 1: Fibonacci d-values (Laws 2–5)

Each planet’s inclination amplitude is governed by a Fibonacci divisor d. These are pure Fibonacci numbers (3, 5, 21, 34) forming mirror-symmetric pairs across the asteroid belt. They determine the scalar balance — the genuine physical constraint that selects the unique configuration.

Level 2: ICRF perihelion periods are H/Fibonacci (Law 1)

All major precession periods divide the Earth Fundamental Cycle by Fibonacci numbers: H/3, H/5, H/8, H/13. This extends to all planets: each ICRF perihelion period is an H-based fraction. The Solar System Resonance Cycle (8H = 2,682,536 years) is the super-period where all planetary cycles complete integer revolutions simultaneously.

Level 3: Ascending node periods are 8H/N

The ascending node regression periods are integer divisors of the Solar System Resonance Cycle (8H/N), with Jupiter and Saturn locked to a shared N=36. Across all 7 fitted planets, the 8H/N integers reproduce JPL’s J2000-fixed-frame ascending-node trends with a cumulative residual error of ~5.8″/century (≈0.8″/century per planet; see Invariable Plane: Ascending Node Periods).

System Reset

The System Reset (≈ 2,649,854 BC) is the epoch within each Solar System Resonance Cycle when all seven fitted planets simultaneously reach their inclination extremes — in-phase planets at minimum, Saturn at maximum. At intermediate balanced years (spaced one Earth Fundamental Cycle apart), only a subset of planets are at their extremes. The System Reset occurs once per 8H and marks the start of the current octave. See Fundamental Cycles: System Reset for details.

Scalar balance vs vector balance

An important distinction: the scalar balance (Laws 3 and 5) is the genuine physical constraint. It determines the Fibonacci d-values and cannot be achieved by arbitrary configurations. The vector balance — whether the 2D angular momentum perturbation vectors cancel at all times — is guaranteed by eigenmode structure for ANY set of frequencies, not just specific ones. This is because the eigenmode solver has 33 spare degrees of freedom (56 parameters minus 23 constraints). The scalar balance is what makes the model’s predictions unique; the vector balance is a mathematical property of the solver, not a validation of specific ascending node rates.

10. What Remains Unknown

The Fibonacci Laws predict orbital properties with remarkable precision, but several questions remain open:

• What determines the base eccentricities? Law 4 predicts all eight eccentricity amplitudes. The base eccentricities (mean values around which each planet oscillates) are derived from the System Reset phase — the epoch (n=7) when every planet passes through its mean eccentricity, with in-phase planets rising (phase 90°) and Saturn falling (phase 270°). This mirrors the inclination alignment at the same epoch. The eccentricity balance (Law 5) emerges naturally at 99.8632%, not forced.

• What produces the eccentricity balance? The inclination balance follows from angular momentum conservation — a well-understood physical principle. What conservation law or perturbation mechanism produces the eccentricity balance? The linear (rather than quadratic) dependence on eccentricity distinguishes it from known conserved quantities like AMD.

• Are the laws universal? The Fibonacci Laws presented here rest entirely on solar system data. Whether any form of Fibonacci structure appears in other stable planetary systems remains an open question that future TTV-characterized exoplanet observations may address.

Summary

← Timekeeping | Physical Origin →