Skip to Content
🚀 Want to help? Get in touch!
The ModelMathematical Foundation

Mathematical Foundations

This page provides the mathematical basis for the Holistic Universe Model, including how the 333,888-year Holistic-Year was derived, what constraints it satisfies, data sources, comparisons with established models, and how the model can be tested or falsified.


1. How 333,888 Years Was Derived

The Honest Starting Point

First, let’s be clear about what we know and don’t know:

  • We do NOT know why the Holistic-Year is 333,888 years from first principles
  • We DO know that 333,888 is the only value that satisfies all six constraints below simultaneously

The number was found empirically by modeling and iteration, not derived from fundamental physics. This is similar to how Kepler found his laws empirically before Newton explained them theoretically.

The Six Constraints

The length of 333,888 years is the only value that satisfies all of the following:

#ConstraintWhat It Requires
11246 AD AlignmentPerihelion must align with December solstice in 1246 AD (verified by Meeus’s formula)
2Longitude of PerihelionMust match observed progression from 90° (1246 AD) to 102.95° (2000 AD)
3Climate CyclesMust produce ~3 × 100k year pattern visible in ice core temperature records
4Eccentricity RangeMust produce eccentricity values matching observations (~0.014 to ~0.017)
5Whole Days per CycleNumber of solar days in a perihelion precession cycle must be an integer
6Mercury PrecessionMust be compatible with observed Mercury perihelion precession (~5600”/century)

The Derivation Process

  1. Start with 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. Check planetary compatibility (Mercury precession)

Only 333,888 years satisfies all constraints. Other values fail one or more tests.

Why This Number?

333,888 = 2⁵ × 3 × 13 × 269
FactorPurpose
3Gives whole inclination precession cycles (333,888 ÷ 3 = 111,296)
13Gives whole axial precession cycles (333,888 ÷ 13 = 25,683.69)
8 = 2³Gives whole obliquity cycles (333,888 ÷ 8 = 41,736)
16 = 2⁴Gives whole perihelion precession cycles (333,888 ÷ 16 = 20,868)
269Required to satisfy the solar day integer constraint

The factor 269 (a prime) is not arbitrary - it’s the value needed for the number of solar days per perihelion precession cycle to be an integer (7,621,874 days in 20,868 years).


2. The Fibonacci Observation

What We Observe

The ratio between inclination precession and axial precession is remarkably close to consecutive Fibonacci numbers:

T_incl / T_axial = 111,296 / 25,683.69 = 4.3333... = 13/3

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

What This Means (and Doesn’t Mean)

Important Distinction: The Fibonacci ratio is an observation, not an explanation. The model does not claim to know WHY this ratio exists - only that it DOES exist and produces accurate predictions.

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 Cycle Table

From the Holistic-Year, all cycles are derived by division:

CycleDivisorDuration (years)Fibonacci?
Holistic-Year1333,888F(1) = 1
Inclination Precession (ICRF)3111,296F(4) = 3
Inclination Precession (Ecliptic)566,777.6F(5) = 5
Obliquity Cycle841,736F(6) = 8
Axial Precession1325,683.69F(7) = 13
Perihelion Precession1620,868No (but 16 = 13 + 3)

Perihelion Precession Derivation

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

Earth orbits EARTH-WOBBLE-CENTER: clockwise, period = 25,684 years PERIHELION-OF-EARTH orbits Sun: counter-clockwise, period = 111,296 years Meeting frequency = 1/T_axial + 1/T_incl (opposite directions, so ADD frequencies) = 1/25,684 + 1/111,296 = 1/20,868 Therefore: They meet every 20,868 years

Note: 16 = 13 + 3, which is why 333,888 ÷ 16 gives the perihelion precession period.


3. Mean Values vs Current Values

The Key Distinction

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

ParameterModel Mean ValueCurrent ObservedDifference
Axial precession period25,683.69 yr~25,772 yr-88 yr (-0.34%)
Inclination precession (ICRF)111,296 yr~112,000 yr-704 yr (-0.63%)
Obliquity cycle41,736 yr~41,000 yr+736 yr (+1.8%)
Perihelion precession20,868 yr~21,000 yr-132 yr (-0.63%)

Is This Unfalsifiable?

A valid concern: if any discrepancy can be attributed to “not being at mean,” is the model testable?

Answer: Yes, because:

  1. The model predicts specific values at specific dates (not just means)
  2. The model predicts how values change over time (specific rates)
  3. These predictions can be compared to observations over decades

4. Calibration vs Prediction

What Was Calibrated (Inputs)

These values were used as inputs to construct the model:

InputValueSource
1246 AD alignmentPerihelion at December solsticeMeeus’s formula
Longitude of perihelion (J2000)102.95°NASA Planetary Fact Sheet
Obliquity (J2000)23.439291°IAU 2006
Eccentricity (J2000)0.01671022NASA Planetary Fact Sheet
Sidereal year31,558,149.724 sJPL Horizons

What Is Predicted (Outputs)

These values are predictions of the model, not inputs:

PredictionModel ValueComparison ValueAgreement
Obliquity at -10,00024.23°24.23° (Laskar)✓ Exact
Obliquity at +10,00022.41°22.40° (Laskar)✓ ±0.01°
Perihelion longitude 1000 AD85.3°85.4° (Meeus)✓ ±0.1°
Perihelion longitude 2500 AD111.4°111.3° (Meeus)✓ ±0.1°
Obliquity range22.15° - 24.68°22.1° - 24.5° (standard)✓ Close

The model was not tuned to match Laskar’s obliquity formula or Meeus’s perihelion values for dates other than 1246 AD. The agreement is a genuine prediction.


5. Comparison with Standard Theory

Where the Model Agrees

PhenomenonModelStandard TheoryAgreement
Axial precession rate~50.3″/yr50.2875″/yr (IAU)✓ Within 0.1%
Obliquity (J2000)23.439°23.439291°✓ Exact
Perihelion progression~17.25″/yr~17.19″/yr (Meeus)✓ Within 0.4%
Obliquity variation±1.27°±1.2° (Laskar)✓ Close

Where the Model Disagrees

PhenomenonModelStandard TheoryTestable?
Eccentricity cycle20,868 years~100k/400k yearsYes - future decades
Eccentricity range0.0139 - 0.01670.0047 - 0.0747Yes - future centuries
Long-term obliquityReturns to meanContinues changingYes - geological record
Climate driverObliquity + InclinationEccentricity (100k)Yes - ice core analysis
Historic length of solar year in daysSolar year more or less the same as todayin 1246 AD about 3 seconds longer than todayNo
Precession predictionsPrecession will have a turning point and start to increase in ~200 yearsPrecession will further decreaseYes - time
Mercury’s 43 arcsec anomalyEarth’s reference frame motion (wobble + PERIHELION-OF-EARTH)General Relativity space-time curvatureYes - anomaly should decrease to ~34”/century

Mercury’s Perihelion Anomaly: The model proposes that the famous ~43 arcsecond “anomaly” in Mercury’s perihelion precession is not caused by relativistic effects, but by Earth’s moving reference frame. The prediction: this value will decrease from ~40 to ~34 arcseconds/century as axial precession increases. See Mercury Precession for details.


6. Data Sources

Primary Sources

ConstantValueSource
J2000 Epoch2000-01-01 12:00 TTIAU Resolution B1.9 (2000)
Astronomical Unit149,597,870.700 kmIAU Resolution B2 (2012)
Earth Eccentricity (J2000)0.01671022NASA Planetary Fact Sheet 
Obliquity (J2000)23.439291°IAU 2006
Sidereal Year31,558,149.724 sJPL Horizons 
Axial Precession Rate50.2875″/yearIAU 2006 Resolution B1

Secondary Sources

Data TypeSourceReference
Obliquity formulasLaskar et al. (1993)A&A 270, 522-533
Longitude of perihelionMeeus (1998)Astronomical Algorithms, Ch. 26
Precession theoryCapitaine et al. (2003)A&A 412, 567-586
Invariable planeSouami & Souchay (2012)A&A 543, A133
Planetary ephemeridesJPL DE440/441JPL Solar System Dynamics 

7. Testable Predictions

What Would Falsify the Model?

The model would be falsified if:

ObservationWould Falsify If
Eccentricity continues decreasing linearly to ~0The 20,868-year cycle doesn’t exist
~100k climate cycle proven to be eccentricity-drivenModel’s climate mechanism is wrong
Axial precession continues downwardsThe length of solar year is not driven by the difference to the obliquity mean

Specific Predictions

Short-term (verifiable within decades)

PredictionModel ValueStandard TheoryHow to Verify
Obliquity 2050 AD23.41°~23.41°Satellite measurements
Eccentricity 2050 AD~0.01668~0.01665Ephemeris comparison
Perihelion date 2050~Jan 4.5~Jan 4-5Direct observation

Medium-term (verifiable within centuries)

PredictionModel ValueStandard TheoryTimeframe
Obliquity 2500 AD23.30°23.30°2500
Eccentricity 2500 AD~0.01645~0.01632500
Eccentricity minimum~0.0139 (11,680 AD)~0 (27,000 AD)Diverges significantly
RA at max obliquityDecreasing from 6hFixed at 6hVerifiable by ~6,000 AD

The RA at Maximum Obliquity Prediction

The Sun’s Right Ascension at maximum obliquity (solstices) appears fixed at RA 6h (June) and 18h (December). However, the model predicts this value slowly decreases due to the interaction between axial tilt and inclination tilt.

PropertyValue
Mean RA at max obliquity~5h48m50s / ~17h48m50s
Oscillation amplitude±11 minutes
Cycle period41,736 years
Peak value~6h00m00s (reached in 1246 AD)
Predicted value by 6,000 AD~5h58m22s
Prediction graph showing RA at maximum obliquity peaked at 6h in 1246 AD and slowly decreasing to 5h58m22s by year 6,000 AD

Why this matters: No reference to changing RA values at solstices exists in current astronomical theory. If this shift becomes observable, it would be a clear validation of the model’s prediction that axial and inclination tilts interact in a specific pattern.

This may also explain long-term variations in magnetic declination .

Long-term (theoretical, not directly verifiable)

PredictionModel ValueTimeframe
Minimum obliquity22.15°~+170,000 years
Maximum obliquity24.68°~+340,000 years
Next perihelion-solstice alignment22,114 AD~20,000 years

8. Uncertainties and Limitations

Known Limitations

AspectLimitationImpact
EccentricityModel uses 20,868-year cycle; standard uses ~100k/400kLong-term predictions diverge
Delta-TEarth’s rotation rate varies unpredictablyDay length predictions uncertain
n-body effectsModel simplifies to two-body interactionsSmall perturbations not modeled
Planet orbitsNot yet fully modeled in simulationFuture work needed

Explicit Assumptions

  1. Stable solar system - The 333,888-year cycle assumes orbital stability over this timescale
  2. Two-point model - EARTH-WOBBLE-CENTER and PERIHELION-OF-EARTH 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, not theoretically derived

What the Model Does NOT Explain

  • Why the Fibonacci ratio exists
  • Why 333,888 specifically (vs some other number)
  • What causes the two precession motions
  • Whether the cycles are truly eternal or slowly changing

9. Reproducibility

Calculation Steps

To reproduce the model’s core values:

  1. Start with 1246 AD alignment - Use Meeus’s formula to verify perihelion at 90° longitude
  2. Calculate perihelion progression - From 90° (1246 AD) to 102.95° (2000 AD) = 12.95° in 754 years
  3. Derive perihelion precession rate - 360° ÷ (12.95°/754 yr) = ~20,950 years (approximate)
  4. Apply Fibonacci constraint - Find 13:3 ratio giving mean axial = 25,684 years
  5. Calculate Holistic-Year - 13 × 25,683.69 = 333,888 years
  6. Verify all divisibility constraints - 333,888 ÷ 16 all give near-integers

Verification Tools


10. References

  1. Capitaine, N., Wallace, P. T., & Chapront, J. (2003). “Expressions for IAU 2000 precession quantities”. Astronomy & Astrophysics, 412, 567-586.

  2. Laskar, J. (1993). “Orbital, precessional and insolation quantities for the Earth from -20 Myr to +10 Myr”. Astronomy & Astrophysics, 270, 522-533.

  3. Meeus, J. (1998). Astronomical Algorithms (2nd ed.). Willmann-Bell.

  4. Souami, D., & Souchay, J. (2012). “The solar system’s invariable plane”. Astronomy & Astrophysics, 543, A133.

  5. Muller, R. A., & MacDonald, G. J. (1997). “Glacial cycles and astronomical forcing”. Science, 277, 215-218.

  6. Hays, J. D., Imbrie, J., & Shackleton, N. J. (1976). “Variations in the Earth’s orbit: Pacemaker of the ice ages”. Science, 194, 1121-1132.


Summary

QuestionAnswer
Why 333,888 years?Only value satisfying all six constraints simultaneously
Why Fibonacci?Observed empirically; physical cause unknown
Is this falsifiable?Yes - specific dated predictions can be tested
Where does it differ from standard theory?Eccentricity cycle (20,868 vs 100k/400k years)
What’s calibrated vs predicted?5 inputs; obliquity/perihelion at other dates are predictions

Return to How It Works for the conceptual overview, or continue to Precession to see how these mathematical foundations manifest in observable phenomena.

Last updated on: