Mathematical Foundations

This page provides the mathematical basis for the Holistic model: how H = 335,317 years was derived, what constraints it satisfies, data sources, comparisons with established models, and how the model can be tested or falsified.

How H was derived

H was found empirically by modelling and iteration, not derived from fundamental physics. This is similar to how Kepler found his laws empirically before Newton explained them theoretically. We do not know why the Earth Fundamental Cycle is 335,317 years from first principles; we do know 335,317 is the most likely candidate that satisfies eight independent constraints simultaneously.

The eight constraints

Derivation process

1. Start with the observed 1246 AD alignment (from Meeus’s formula for longitude of perihelion).
2. Model precession rates to match observed progression to 2000 AD.
3. Find integer ratios that produce whole-number cycles.
4. Test against climate data (ice-core ~100k pattern).
5. Verify eccentricity range matches observations.
6. Verify obliquity range matches observations.
7. Verify days and years match observations.
8. Check planetary compatibility (Mercury precession).

335,317 years satisfies all constraints best. Other values fail one or more tests.

ω = 102.93735° + 1.71953°T + 0.00046°T² + ...

Why this number?

335,317 = 23 × 61 × 239

Five Fibonacci-related divisors connect H to the five Milankovitch-type cycles:

Divisors 5, 8, and 16 produce integer day counts in TOTAL_DAYS_IN_H (so each cycle completes a whole number of solar days); 3 and 13 do not. All five precession periods vary within each cycle around the mean (the current axial period is ~25,771 years vs the model’s mean of ~25,794). The number of solar days per perihelion precession cycle is exactly an integer (7,654,495 days in ~20,957 years).

The five cycles are not independent: standard orbital mechanics derives the obliquity period and perihelion precession period as beat frequencies of the others, and with all five expressed as H/n these become Fibonacci identities (13 − 5 = 8, 13 + 3 = 16). Full analysis with the H/3 / H/5 dual roles and Jupiter ecliptic period coupling: Supporting Evidence §11.

The Fibonacci observation

The ratio between inclination precession and axial precession is remarkably close to a ratio of two Fibonacci numbers:

T_incl / T_axial = ~111,772 / 25,793.62 = 4.3333... = 13/3

Both 3 and 13 are Fibonacci numbers (F₄ and F₇).

Possible interpretations: (1) coincidence — the ratio happens to be close to 13/3; (2) resonance — orbital mechanics naturally settle into stable integer ratios; (3) deeper physics — some unknown principle selects Fibonacci ratios. The model remains agnostic on the cause. What matters is that the ratio produces accurate predictions.

The Fibonacci structure has rigorous theoretical grounding in the Kolmogorov–Arnold–Moser (KAM) theorem: orbits with frequency ratios closest to the golden ratio φ are maximally stable against perturbation, because φ is the irrational number hardest to approximate by ratios of small integers, and successive Fibonacci ratios (3/2, 5/3, 8/5, 13/8, …) converge to φ. Empirical surveys (Pletser 2019 , Aschwanden & Scholkmann 2017 ) confirm Fibonacci clustering in solar-system and exoplanet period ratios. Full theoretical treatment: Physical Origin; supporting evidence: Supporting Evidence §2.

The model extends Fibonacci structure beyond period ratios: six independent Fibonacci Laws connect planetary precession periods, eccentricities, and inclination amplitudes through Fibonacci numbers and the mass-weighted quantity $e \times \sqrt{m}$.

Perihelion precession derivation

The ~20,957-year perihelion precession emerges from the meeting frequency of two counter-rotating motions:

Earth orbits its wobble center:       clockwise,        period = ~25,794 years
Earth's perihelion point orbits Sun:  counter-clockwise, period = ~111,772 years

Meeting frequency = 1/T_axial + 1/T_incl (opposite directions, so ADD frequencies)
                = 1/~25,794 + 1/~111,772
                = 1/~20,957

Note: 16 = 13 + 3, which is why H ÷ 16 gives the perihelion precession period. Full cycle table: Precession and Fibonacci Laws §Law 1.

Mean values vs current values

The model predicts mean values over the full 335,317-year cycle. Currently observed values differ because we are at a specific position in the cycle, not at the mean.

Is this unfalsifiable? A valid concern: if any discrepancy can be attributed to “not being at mean,” is the model testable? Yes — the model predicts specific values at specific dates (not just means) and predicts how those values change over time (specific rates). Those predictions can be compared to observations over decades.

Calibration vs prediction

Degrees of freedom

The model has 6 free parameters — all governing the Earth simulation; the planetary Fibonacci configuration adds no additional degrees of freedom (a unique mirror-symmetric solution emerges from exhaustive search):

Total: 6 free parameters — the planet configuration assigns a period, quantum number d, and inclination cycle anchor to each planet. Periods and cycle anchors are observationally constrained; only the d-assignment is a free choice — a discrete selection from 7,558,272 possible Fibonacci assignments, narrowed by four successive physical filters (inclination balance, eccentricity balance, Laplace–Lagrange bounds, direction match) to 15 viable configurations, of which only one is mirror-symmetric. Comparable to standard astronomical models (Laskar uses ~6 parameters for obliquity).

What was calibrated (inputs)

These values were directly used to determine the model’s parameters:

What is predicted (NOT used in calibration)

These values are genuine predictions — they were NOT used to construct the model:

*Eccentricity and inclination J2000 values were NOT used to find H. The model predicts them from the structure.

What is constrained (validation, not proof)

The circularity concern

If the eight constraints in §1 were all used to find H, then matching them doesn’t validate the model — it just confirms the fitting worked. Not all constraints are equal:

1. True inputs (1246 AD alignment, J2000 values): the model was tuned to match these.
2. Structural constraints (Fibonacci divisors, integer days): these constrain the search but don’t guarantee any specific value.
3. Post-hoc validations (Mercury, planetary inclinations): these were checked AFTER H was determined.

Honest assessment: the model’s predictive power should be judged by how well it predicts values at dates other than J2000 and 1246 AD, and whether its long-term predictions (eccentricity cycle, precession reversal) prove correct — NOT by how well it matches the data used to construct it.

Independent historical verification: the 3D simulation has been validated against 623 independently recorded astronomical events spanning approximately 2000 BC to 4000 AD — solstice and equinox dates, perihelion passages, and eclipse timings. Accuracy varies by epoch: ±1 day for ancient observations, ±1 hour for medieval records, ±1 minute for modern measurements. These historical observations were not used to calibrate the model — they serve as independent evidence that the model’s orbital mechanics produce correct results across millennia. Full dataset: verification data reference .

Comparison with standard theory

Where the model agrees

All values are J2000 instantaneous values — model output vs established astronomical references at the same epoch.

Where the model disagrees

The Mercury anomaly disagreement is canonical at Mercury Precession. The Saturn ecliptic-retrograde permanent-feature claim is canonical at Supporting Evidence §12.

Data sources

Primary sources

Secondary sources

Testable predictions

The model produces 25 specific, testable predictions organised by timeframe (near-term decades through deep-time Gyr-scale). Full details: Predictions.

The model would be falsified if:

*The BepiColombo falsifier assumes the analysis pipeline reports the raw measured perihelion advance rather than a GR-inclusive ephemeris fit total — see Mercury Precession Test methodology.

See Predictions: Verification Pathways for all 25 predictions with specific values, near-term tests (RA shift, BepiColombo, Planet Nine and small-KBO obliquity via LSST), and verification pathways.

Uncertainties and limitations

Explicit assumptions

1. Stable solar system — the 335,317-year cycle assumes orbital stability over this timescale.
2. Two-point model — the wobble center and Earth’s perihelion point are mathematical constructs.
3. Mean values exist — the model assumes precession rates oscillate around fixed means.
4. Fibonacci ratio is real — the 3:13 ratio is empirically observed; KAM theory provides plausibility but not a first-principles derivation of this specific ratio.

What the model does NOT explain

• Why 13:3 specifically (vs other Fibonacci pairs like 8:3, 21:5) — KAM theory predicts Fibonacci stability in general, not this particular ratio.
• Why 335,317 specifically (vs some other number) at the J2000 epoch (the cycles slowly rescale at geological time per Expanding Resonance — the integer labels are invariant; only the absolute periods drift).
• What causes the two precession motions.

Reproducibility

The model provides formulas to calculate obliquity, eccentricity, inclination, longitude of perihelion, and day/year lengths at any year — see Formulas Reference for ready-to-use expressions.

Example comparison — obliquity values calculated using the model’s formula vs established sources:

Verify against external data

1. JPL Horizons (ssd.jpl.nasa.gov/horizons ) — query Earth’s orbital elements for any date; compare perihelion dates, eccentricity, longitude of perihelion.
2. Laskar’s Tables (A&A 1993) — compare obliquity values for ±10,000 years; model matches within ±0.2° for this range.
3. Meeus’s Formulas (Astronomical Algorithms, 1998) — verify 1246 AD perihelion-solstice alignment; compare longitude of perihelion progression.
4. 3D Simulation (3d.holisticuniverse.com ) — visual verification of all precession movements; adjust date and observe changes in real-time.

How the model was derived (summary)

1. Anchor point: Meeus’s formula places perihelion at 90° longitude (December solstice) in 1246 AD.
2. Observed progression: longitude moved from 90° (1246 AD) to 102.947° (2000 AD) = 12.947° in 754 years.
3. Fibonacci constraint: the 13:3 ratio between axial and inclination precession was identified empirically.
4. Earth Fundamental Cycle: 13 × 25,793.62 years ≈ 335,317 years.
5. Balanced Year: 1246 AD − (14.5 × ~20,957) ≈ 302,635 BC.

All other values (obliquity range, eccentricity range, day/year lengths) follow from these foundational parameters.

← Mercury Precession | Supporting Evidence →