@essrt/physics 22.1.0 → 23.2.0

package/src/model.js

@@ -17,7 +17,6 @@
  * src/lib/essrt.ts) operation-for-operation where FP association matters,
  * which itself mirrors tools/lib/constants.js §9.
  */
-import { deriveEpochParams } from './layer0/derive-params.js';
 import { buildPlanetChainsFromArtifactData, computePlanetElementsAtYear } from './planets/keplerian-chain.cjs';
 import { createSecularSeriesOverride } from './planets/secular-series.cjs';
 import { computeEquatorNodeOriginSFrameDeg, convertNodeSFrameToEquatorOriginDeg } from './planets/inv-plane-frame.cjs';
@@ -26,13 +25,11 @@ import { createDeepOrbitalHistory } from './earth/deep-orbital-history.cjs';
 import { CHAIN_ARTIFACT } from './planets/chain-artifact.js';
 import { createPlanetSpinChannelFromArtifacts, computeObliquityJ2000Deg } from './planets/spin-channel.cjs';
 import { PLANET_DISPLAY_NAMES } from './planets/display-names.cjs';
-import { createDeltaTCycles } from './deltat/cycles.cjs';
-import { createDeepTimeLod } from './deltat/deep-time.cjs';
-import { createMoonRecessionHistory, createSolarChannelBudget } from './deltat/recession-history.cjs';
-import { evalClimateL1OrbitalPermil, createAlphaGiaChannel } from './climate/l1-orbital.cjs';
-import { createDeepEccChannel } from './moon/deep-ecc-channel.cjs';
+import { createDeepTimeChain } from './deltat/deep-time-chain.js';
+import { computeMassChainJ2000 } from './earth/mass-chain.cjs';
+import { MOON_DISTANCE_ARTIFACT } from './deltat/moon-distance-artifact.cjs';
+import { evalClimateL1OrbitalPermil } from './climate/l1-orbital.cjs';
 import { DEEP_MODES_ARTIFACT } from './moon/deep-modes-artifact.cjs';
-import { createMoonMonthChain } from './moon/month-chain.cjs';
 import { createChainCycleIntegrator } from './chain-cycles/index.cjs';
 import { createMoonArguments, jdToDecimalYear } from './moon/arguments.cjs';
 import { createMoonSeries } from './moon/series.cjs';
@@ -51,16 +48,9 @@ import { driver2PeriodSecondsAtAge } from './planets/orbit-chain.cjs';
 // here into F.RA_DAY_OFFSET_MS; plan 09 step 5b deleted that family — it was
 // the equation of time's window term of an instrument, not a motion.)
 
-/**
- * Moon-channel eccentricity sensitivities (perigee/node), the [g/g₀]^s
- * exponents of the factored deep-time law. Same known-exception class as the
- * RA offsets above: single source src/script.js _FW_MOON, mirrored as
- * literals in tools/lib/deep-time.js (_ECOMP_S_W/_ECOMP_S_N). Both are
- * Meeus-effective — S_N moved 1.0 → 1.018 with the v4 frame-attribution
- * batch.
- */
-const MOON_ECC_SENSITIVITY_PERIGEE = 2.407;
-const MOON_ECC_SENSITIVITY_NODE = 1.018;
+// (The Moon-channel eccentricity sensitivities s_ϖ/s_Ω — the [g/g₀]^s exponents of
+// the factored deep-time law — have ONE home since plan 13 B2:
+// deltat/deep-time-chain.js MOON_ECC_SENSITIVITY_PERIGEE / _NODE.)
 
 /**
  * Assemble the model surfaces from a resolved constants context + fitted
@@ -108,7 +98,7 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   // comb) left with that device; no comb phase is evaluated here any more.
 
   const meanSiderealYearDays = C.yearLengthRef.siderealYear;
-  const meanSiderealYearSeconds = meanSiderealYearDays * 86400;
+  // (the sidereal year in SI seconds is the mass chain's own, below)
   // Plan 09: the J2000 day length ω₀ is the SI day (an input); the device's
   // H/13 kinematic-sidereal-days construction (86,399.99968 s) is deleted.
   const meanLengthOfDay = C.yearLengthRef.solarDay;
@@ -130,7 +120,6 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   const G_CONSTANT = C.physicalConstants.G_CONSTANT;
   const MASS_RATIO_EARTH_MOON = C.physicalConstants.MASS_RATIO_EARTH_MOON;
   const currentAUDistance = C.physicalConstants.currentAUDistance;
-  const earthMoiFactorJ2000 = C.physicalConstants.earthMoiFactorJ2000;
   const moonDistanceKm = C.moonReference.moonDistance;
   const moonSiderealMonthInput = C.moonReference.moonSiderealMonthInput;
 
@@ -138,171 +127,54 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   // eighth-integer-per-8H count read it back 7.8 ms long; deleted).
   const moonSiderealMonth = moonSiderealMonthInput;
 
-  // Mass chain: Moon Kepler → GM_EM → Earth/Moon split → Sun (§9 order)
-  const moonOrbitalShift = moonDistanceKm * (1 / (MASS_RATIO_EARTH_MOON + 1)) * (moonSiderealMonth / meanSiderealYearDays);
-  const moonDistanceCorrected = moonDistanceKm + moonOrbitalShift;
-  const GM_EARTH_MOON_SYSTEM = (4 * Math.PI * Math.PI * Math.pow(moonDistanceCorrected, 3))
-    / Math.pow(moonSiderealMonth * meanLengthOfDay, 2);
-  const GM_EARTH_ALONE = GM_EARTH_MOON_SYSTEM * (MASS_RATIO_EARTH_MOON / (MASS_RATIO_EARTH_MOON + 1));
-  const GM_MOON_ALONE = GM_EARTH_MOON_SYSTEM / (MASS_RATIO_EARTH_MOON + 1);
-  const GM_SUN_PLUS_EARTH = (4 * Math.PI * Math.PI * Math.pow(currentAUDistance, 3))
-    / Math.pow(meanSiderealYearSeconds, 2);
-  const GM_SUN = GM_SUN_PLUS_EARTH - GM_EARTH_ALONE;
-  const M_SUN = GM_SUN / G_CONSTANT;
-
-  // ── α(t): the GIA channel — the LAGGED response to the L1 ice history, with the R2 lattice pin ──
-  const CLIMATE = F.CLIMATE_FORMULA_COEFFS;
-  const CLIMATE_REGIME = CLIMATE.regimes['lr04-post-mpt'];
+  // The mass chain (Moon Kepler → GM_EM → the Earth/Moon split → the Sun, §9
+  // order) — ONE home since plan 13 B3: earth/mass-chain.cjs (the deep-time
+  // chain below computes the same scalars from the same context).
+  const mass = computeMassChainJ2000({
+    moonDistanceKm, moonSiderealMonthDays: moonSiderealMonth, meanSiderealYearDays,
+    meanLengthOfDaySeconds: meanLengthOfDay, massRatioEarthMoon: MASS_RATIO_EARTH_MOON,
+    currentAUDistanceKm: currentAUDistance, gravitationalConstantKm3KgS2: G_CONSTANT,
+  });
+  const moonOrbitalShift = mass.moonOrbitalShiftKm;
+  const moonDistanceCorrected = mass.moonDistanceCorrectedKm;
+  const GM_EARTH_MOON_SYSTEM = mass.gmEarthMoonSystemKm3S2;
+  const GM_EARTH_ALONE = mass.gmEarthAloneKm3S2;
+  const GM_MOON_ALONE = mass.gmMoonAloneKm3S2;
+  const GM_SUN_PLUS_EARTH = mass.gmSunPlusEarthKm3S2;
+  const GM_SUN = mass.gmSunKm3S2;
+  const M_SUN = mass.sunMassKg;
 
+  // ── THE DEEP-TIME CHAIN — ONE assembly (plan 13 B2) ──────────────────────
+  // The GIA α channel (the lagged response to the L1 ice history, with the R2
+  // lattice pin and the exact-argument memo), the Layer-0 bundle, the cited
+  // Earth–Moon distance history (plan 12 option A), the ΔT cycle stack (all
+  // channels ON here), |s₃|, the LOD/ΔT core and the lunar month chain are
+  // assembled once in deltat/deep-time-chain.js from this context's constants
+  // and fitted coefficients (so a counterfactual context flows through them);
+  // the Node twin and the browser destructure the same module.
+  const chain = createDeepTimeChain({
+    constants: C, fitted: F, deepModes: DEEP_MODES_ARTIFACT, moonDistanceArtifact: MOON_DISTANCE_ARTIFACT,
+    ageAnchorYear: startmodelYear,
+  });
+  const { alphaGia, earthMoiFactorAtAge, EPOCH_PARAMS, recession, dtCycles, dtCycleLodCorrectionSum, s3ArcsecPerYr, deepLod, deepEcc, moonChain } = chain;
+  const CLIMATE_REGIME = F.CLIMATE_FORMULA_COEFFS.regimes['lr04-post-mpt'];
   /** The climate formula's L1 (unlagged) — the orbital δ¹⁸O layer itself. @param {number} year @returns {number} */
   const evalClimateL1 = (year) => evalClimateL1OrbitalPermil(year, {
     l1Terms: CLIMATE_REGIME.L1,
     yStdDenormalization: CLIMATE_REGIME.denormalization.y_std,
   });
-  // ONE home (climate/l1-orbital.cjs): k derived from the Cox–Chao rate, τ the
-  // record-measured relaxation time (plan 06 D7). Same construction in
-  // src/script.js, tools/lib/deep-time.js and the website's essrt.ts.
-  const alphaGia = createAlphaGiaChannel({
-    l1Terms: CLIMATE_REGIME.L1,
-    yStdDenormalization: CLIMATE_REGIME.denormalization.y_std,
-    relaxationKyr: C.deepTime.alphaGiaRelaxationKyr,
-    alphaGiaRateJ2000PerYr: C.deepTime.alphaGiaRateJ2000PerYr,
-    alphaJ2000: earthMoiFactorJ2000,
-    // the channel's second observable, J₂(t)/J₂₀ — scaled on the OBSERVED
-    // dJ₂/dt; consumed by the spin integration below (the precession of date)
-    j2RateJ2000PerYr: C.deepTime.j2RateJ2000PerYr,
-    j2J2000: C.physicalConstants.earthJ2,
-  });
-
-  // Exact-argument memo (mirrors src/script.js + tools/lib/deep-time.js —
-  // the round-5 perf campaign): same tMa → bit-identical value; bounded by
-  // wholesale clear. (Plan 09 step 7a: the lattice-α pin branch left with
-  // the ∫dt/H_era phase table, the one table this model built under it.)
-  /** @type {Map<number, number>} */
-  const alphaMemo = new Map();
-  /** @param {number} tMa @returns {number} */
-  const earthMoiFactorAtAge = (tMa) => {
-    let v = alphaMemo.get(tMa);
-    if (v === undefined) {
-      v = alphaGia.alphaAt(2000 - tMa * 1e6);
-      if (alphaMemo.size >= 8192) alphaMemo.clear();
-      alphaMemo.set(tMa, v);
-    }
-    return v;
-  };
-
-  // ── Layer 0 + deep-time core ──────────────────────────────────────────────
-  const EPOCH_PARAMS = deriveEpochParams({
-    solarLuminosityW: C.physicalConstants.solarLuminosityW,
-    solarWindKgPerS: C.physicalConstants.solarWindMassLossKgPerS,
-    speedOfLightKmPerS: C.physicalConstants.speedOfLight,
-    alpha1PerMa: C.deepTime.alpha1PerMa,
-    alpha3PerMa3: C.deepTime.alpha3PerMa3,
-    alpha4PerMa4: C.deepTime.alpha4PerMa4,
-    meanSiderealYearSeconds,
-    lodJ2000Seconds: meanLengthOfDay,   // plan 09: ω₀ = the SI day (was the H/13 kinematic day)
-    sunMassKg: M_SUN,
-    gmEarthAloneKm3S2: GM_EARTH_ALONE,
-    gmMoonAloneKm3S2: GM_MOON_ALONE,
-    gravitationalConstantKm3KgS2: G_CONSTANT,
-    earthMoiFactorJ2000,
-    earthDiameterKm: C.bodyDiametersKm.earth,
-    moonDistanceKm,
-    moonOrbitalEccentricity: C.moonReference.moonOrbitalEccentricityBase,
-    gmEarthMoonSystemKm3S2: GM_EARTH_MOON_SYSTEM,
-    // the solar precession torque's inputs (plan 06 Phase 3: the unit H(t) = 13·T_p,composed)
-    gmSunKm3S2: GM_SUN,
-    astronomicalUnitKm: currentAUDistance,
-    earthOrbitalEccentricityJ2000: C.earthOrbital.earthEccentricityJ2000,
-    moonEclipticInclinationJ2000Deg: C.moonReference.moonEclipticInclinationJ2000,
-    // p₀'s period — the tropical year of date's equinox term (plan 09 step 5d)
-    axialPrecessionJ2000Years: certifiedAxialPrecessionJ2000Years(),
-  });
-
-  // ── Driver 1½: regime-aware recession history + the solar channels ────────
-  // The quartic stays bit-identical ≤ jointMa (the Wells/Wu-gated era);
-  // beyond, the fitted staircase spline to the Roche crossing, and the
-  // ocean-leak/thermal-pump channels make L_EM time-dependent.
-  const REGIME = C.deepTime.recessionRegime;
-  const recession = createMoonRecessionHistory({
-    aMoonNowMetres: EPOCH_PARAMS.moonDistanceNowM,
-    alpha1PerMa: EPOCH_PARAMS.alpha1PerMa,
-    alpha3PerMa3: EPOCH_PARAMS.alpha3PerMa3,
-    alpha4PerMa4: EPOCH_PARAMS.alpha4PerMa4,
-    regime: {
-      jointMa: REGIME.jointMa,
-      knotAgesMa: REGIME.knotAgesMa,
-      knotDistancesKm: REGIME.knotDistancesKm,
-      genesisMa: REGIME.genesisMa,
-      rocheLimitKm: REGIME.rocheLimitKm,
-    },
-  });
   const moonDistanceMetresAtAge = recession.distanceMetresAtAge;
-  const solarBudget = createSolarChannelBudget({
-    lTotalJ2000KgM2S: EPOCH_PARAMS.totalAngularMomentumKgM2S,
-    mMoonAloneKg: EPOCH_PARAMS.moonMassKg,
-    gmEmM3PerS2: EPOCH_PARAMS.gmEarthMoonM3S2,
-    eFactorMoon: EPOCH_PARAMS.moonEccentricityFactor,
-    beta0: REGIME.solarOceanLeakBeta0,
-    pumpStartMa: REGIME.thermalPumpStartMa,
-    pumpEndMa: REGIME.thermalPumpEndMa,
-    pumpFactor: REGIME.thermalPumpFactor,
-    jointMa: REGIME.jointMa,
-    genesisMa: REGIME.genesisMa,
-    distanceMetresAtAge: moonDistanceMetresAtAge,
-  });
-
-  const DT = F.DT_STACK;
-  const RES = F.DT_RESONATOR;
-  // Periods STATED IN YEARS (the fit files' period_yr / T0_yr): plan 06 T5
-  // measured the former 8H/n labels chance-level; layer B item 3 removed the
-  // divisor from every runtime. Twins: tools/lib/deep-time.js, src/script.js,
-  // the website's essrt.ts.
-  const dtCycles = createDeltaTCycles({
-    taperFullHalfwidthYears: C.deepTime.dtStackTaperFullHalfwidthYr,
-    taperTotalHalfwidthYears: C.deepTime.dtStackTaperTotalHalfwidthYr,
-    tropicalYearSecondsJ2000: meanTropicalYearJ2000Seconds,
-    cycles: {
-      bond: { periodYears: DT.bond.period_yr, cosCoeffSeconds: DT.bond.cos_coeff_s, sinCoeffSeconds: DT.bond.sin_coeff_s },
-      hallstatt: { periodYears: DT.hallstatt.period_yr, cosCoeffSeconds: DT.hallstatt.cos_coeff_s, sinCoeffSeconds: DT.hallstatt.sin_coeff_s },
-      jose5: { periodYears: DT.jose5.period_yr, cosCoeffSeconds: DT.jose5.cos_coeff_s, sinCoeffSeconds: DT.jose5.sin_coeff_s },
-      jose4: { periodYears: DT.jose4.period_yr, cosCoeffSeconds: DT.jose4.cos_coeff_s, sinCoeffSeconds: DT.jose4.sin_coeff_s },
-    },
-    resonator: {
-      t0Years: RES.T0_yr,
-      q: RES.Q,
-      kicks: RES.kick_epochs_year.map(/** @param {number} t @param {number} i */ (t, i) => ({
-        tYear: t,
-        cosSeconds: RES.kick_coefficients_s[i].cos,
-        sinSeconds: RES.kick_coefficients_s[i].sin,
-      })),
-      tones: RES.drive_tones.map(/** @param {{period_yr: number, phi_locked_rad: number, amp_s: number}} t */ (t) => ({ periodYears: t.period_yr, phiLockedRad: t.phi_locked_rad, ampSeconds: t.amp_s })),
-    },
-  });
-
-  /** Layer-3/4 cyclic δLOD sum (all shipped channels ON). @param {number} year @returns {number} */
-  const dtCycleLodCorrectionSum = (year) =>
-    dtCycles.cycleLodSecondsAt('bond', year)
-    + dtCycles.cycleLodSecondsAt('hallstatt', year)
-    + dtCycles.cycleLodSecondsAt('jose5', year)
-    + dtCycles.cycleLodSecondsAt('jose4', year)
-    + dtCycles.swingLodSecondsAt(year);
-
-  // Plan 06 S6 — |s₃|, the dominant nodal mode of Earth's orbit: the
-  // largest-amplitude ζ mode of the banked deep secular modes (the recipe
-  // the registry's `eclPrecYears`/`obliqCycleYears` and the browser's
-  // obliquity-beat helper use). ONE home here; the deep-time factory's
-  // ecliptic missing-motion term, the lunisolar surface and the API read it.
-  const s3ArcsecPerYr = (() => {
-    const z = DEEP_MODES_ARTIFACT.earthZeta
-      .filter((m) => Math.abs(m.omegaRadPerYr) > 1e-9)
-      .sort((a, b) => Math.hypot(b.re, b.im) - Math.hypot(a.re, a.im))[0];
-    return Math.abs((z.omegaRadPerYr * 180) / Math.PI) * 3600;
-  })();
+  /** The published surfaces' km reading — NaN beyond the cited history's span (the API's
+   *  validity span lies inside it; the deep charts stop where the input stops).
+   *  @param {number} year @returns {number} */
+  const moonDistanceKmAtYearOrNaN = (year) => { const m = moonDistanceMetresAtAge(yearToTMa(year)); return m === null ? NaN : m / 1000; };
 
   /** THE REALIZED POLE RATE OVER A WINDOW, ″/yr — the mean over [yearFrom, yearTo] of the rate the spin
    *  integration actually turns the pole at: the inertial secular rate of date × the solar-torque factor at
    *  e(t) × J₂(t)/J₂₀ × cos ε(t)/cos ε₀ (ds/dt = α (ŝ·n̂)(ŝ×n̂): the azimuthal rate about n̂ is α cos ε).
+   *  The secular rate here is the CONSTANT's scaling (lunisolarTorqueRate…), not the published rate of
+   *  date: the hybrid's own ε(t) carries the ONE cosine, so the cited mean obliquity's factor (plan 12 §8)
+   *  does not enter twice — inside the ζ-series span the hybrid IS the model's obliquity of date.
    *  Measured (2026-10): the dominant obliquity line of the model's own ε(t) over the last 5 Myr sits at
    *  this mean minus |s₃| (40,938 fitted vs 40,936 predicted), NOT at the J2000 beat (41,098) — the pole
    *  turned faster in the past (spin +0.06 ″/yr over 5 Myr, lower mean obliquity +0.05, higher mean
@@ -323,7 +195,7 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
     const kInertial = certifiedAxialPrecessionJ2000Years() / oneSourceM.axialPrecessionYearsLunisolarJ2000;
     let sum = 0, n = 0;
     for (let y = lo; y <= hi; y += stepYears) {
-      const r = deepLod.lunisolarPrecessionRateArcsecPerYrAtAge(yearToTMa(y));
+      const r = deepLod.lunisolarTorqueRateArcsecPerYrAtAge(yearToTMa(y));
       const psiDot = (r === null ? 1296000 / certifiedAxialPrecessionJ2000Years() : r) * kInertial;
       const f = oneSourceM.precessionOfDateFactorsAt(y);
       sum += psiDot * f.solarTorque * (useEllipticity ? f.ellipticity : 1) * Math.cos(oneSourceM.epsAt(y) * Math.PI / 180) / cosEps0;
@@ -332,38 +204,6 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
     return sum / n;
   };
 
-  const deepLod = createDeepTimeLod({
-    constants: {
-      lTotalEmKgm2S: EPOCH_PARAMS.totalAngularMomentumKgM2S,
-      mMoonAloneKg: EPOCH_PARAMS.moonMassKg,
-      mEarthAloneKg: EPOCH_PARAMS.earthMassKg,
-      rEarthMetres: EPOCH_PARAMS.earthRadiusM,
-      gmEmM3PerS2: EPOCH_PARAMS.gmEarthMoonM3S2,
-      eFactorMoon: EPOCH_PARAMS.moonEccentricityFactor,
-      aLockMetres: EPOCH_PARAMS.moonLockDistanceM,
-      aMoonNowMetres: EPOCH_PARAMS.moonDistanceNowM,
-      alpha1PerMa: EPOCH_PARAMS.alpha1PerMa,
-      alpha3PerMa3: EPOCH_PARAMS.alpha3PerMa3,
-      alpha4PerMa4: EPOCH_PARAMS.alpha4PerMa4,
-      lodNowH13Seconds: EPOCH_PARAMS.lodNowH13Seconds,
-      meanSiderealYearJ2000Seconds: meanSiderealYearSeconds,
-      solarMassLossFracPerYear: EPOCH_PARAMS.solarMassLossFracPerYear,
-      precessionSolarShareJ2000: EPOCH_PARAMS.precessionSolarShareJ2000,
-    },
-    moonDistanceMetresAtAge,
-    moiFactorAtAge: earthMoiFactorAtAge,
-    // S5: the composed clock's J2000 anchor — the certified year laws' beat
-    // at 2000 (ONE home below, shared with the hybrid's self-anchor); lazy.
-    precessionPeriodJ2000YearsFn: () => certifiedAxialPrecessionJ2000Years(),
-    // Plan 06 T2 item: the solar day's ecliptic missing-motion term rides the
-    // nodal period 1,296,000/|s₃| (orbital), not the spin unit's H/5.
-    nodalPeriodYearsFn: () => 1296000 / s3ArcsecPerYr,
-    cycleLodSumAt: dtCycleLodCorrectionSum,
-    swingLodAt: (year) => dtCycles.swingLodSecondsAt(year),
-    swingLodRateAt: (year) => dtCycles.swingLodRateAt(year),
-    lEmAtAgeKgm2S: solarBudget.lEmAtAgeKgm2S,
-  });
-
   // (Plan 09 step 7a: the integrated ∫dt/H_era phase table and its counter
   // `cyclesBetween` — the frozen era clock's own convention H_era =
   // H₀·LOD/LOD₀, plan 06 D8's second named counter, which every FROZEN comb
@@ -438,7 +278,7 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   // cardinal braid) reads the one-source series below (plan 06 layer B —
   // the two agree within 4.2e-5 in-era). The laws hook deliberately does
   // NOT reach this channel.
-  const deepEcc = createDeepEccChannel(DEEP_MODES_ARTIFACT);
+  // (`deepEcc` is the deep-time chain's, destructured above.)
   /** @param {number} year @returns {number} */
   // Plan 06 layer B: the PUBLISHED Earth eccentricity is the one-source series
   // (|z|, the secular orbit); the SUN's e adds the derived mean offset (see
@@ -551,6 +391,12 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
     // 108.6 ± 8.5 ″/yr at 2.46 Ga; the composed rate reads 104.5).
     const seriesArt = /** @type {any} */ (secularSeriesArtifact);
     const sb = seriesArt && seriesArt.bodies && seriesArt.bodies.earth;
+    // Plan 12 §8.5 (owner): the hybrid's own domain is the banked ζ series' span (years; null
+    // without a series). Beyond it the ζ-mode tail carries the Quaternary mean with no tidal
+    // evolution — the API refuses the obliquity there and serves the cited mean beside it.
+    const obliquityHybridDomainYears = sb
+      ? /** @type {[number, number]} */ ([2000 + seriesArt.t0Yr, 2000 + seriesArt.t0Yr + (sb.zetaQ.length - 1) * sb.stepYr])
+      : null;
     const hist = createDeepOrbitalHistory({
       zModes: DEEP_MODES_ARTIFACT.earthZ,
       zetaModes: DEEP_MODES_ARTIFACT.earthZeta,
@@ -570,7 +416,10 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
       // rendering-device residual imported as the initial condition of the
       // physics (found when the package began publishing the hybrid).
       obliquityJ2000Deg: C.earthOrbital.obliquityJ2000_deg,
-      axialPrecessionYearsAtYearFn: (yr) => deepLod.lunisolarPrecessionPeriodYearsAtAge((startmodelYear - yr) / 1e6) ?? axial0,
+      // The CONSTANT's scaling (the torque period), not the published rate of date: the
+      // integration's own ŝ·n̂ carries the obliquity's cosine, so the cited mean obliquity's
+      // factor (plan 12 §8 option (a)) must not enter twice. Bit-identical to the pre-§8 form.
+      axialPrecessionYearsAtYearFn: (yr) => deepLod.lunisolarTorquePeriodYearsAtAge((startmodelYear - yr) / 1e6) ?? axial0,
       // THE PRECESSION OF DATE (plan 06 §9 item 10): on top of the secular
       // law above, the solar torque at the eccentricity of date and the
       // oblateness of date from the GIA channel — both exactly 1 at J2000,
@@ -588,19 +437,23 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
     // and the year-length rates (±0.5-yr central differences through the
     // grid) then returned grid-segment averages instead of local rates
     // (anomalistic of date +2.63 s / +19 s, visit-order dependent).
-    const samplers = /** @type {Map<number|'deep', any>} */ (new Map());
-    let deepRangeYr = 0;
+    const samplers = /** @type {Map<number, any>} */ (new Map());
     const gridStep = (/** @type {number} */ n) => (n <= 50000 ? 100 : n <= 2000000 ? 1000 : 5000);
-    const tierSpan = (/** @type {number} */ n) => (n <= 50000 ? 50000 : n <= 2000000 ? 2000000 : Math.ceil(n * 1.25 / 5000) * 5000);
+    const tierSpan = (/** @type {number} */ n) => (n <= 50000 ? 50000 : 2000000);
     const sampleAt = (/** @type {number} */ year) => {
       const t = year - 2000, need = Math.max(20000, Math.abs(t) * 1.25);
-      const span = tierSpan(need);
-      const key = span > 2000000 ? /** @type {'deep'} */ ('deep') : span;
-      if (!samplers.has(key) || (key === 'deep' && span > deepRangeYr)) {
-        if (key === 'deep') deepRangeYr = span;
-        samplers.set(key, hist.build(span, -span, gridStep(span)));
+      if (need > 2000000) {
+        // The deep tier (5000-yr grid): walked in the queried DIRECTION only,
+        // to two nodes past the query, and EXTENDED by the factory's shared
+        // walk when a deeper epoch arrives (plan 13 A1). It was built to
+        // ±1.25·|t| in both directions and rebuilt from J2000 on every deeper
+        // probe — bit-identical values, 4× the RK4 steps of one walk.
+        const span = Math.ceil((Math.abs(t) + 2 * 5000) / 5000) * 5000;
+        return hist.build(t < 0 ? -span : span, 0, 5000).at(t);
       }
-      return samplers.get(key).at(t);
+      const span = tierSpan(need);
+      if (!samplers.has(span)) samplers.set(span, hist.build(span, -span, gridStep(span)));
+      return samplers.get(span).at(t);
     };
     // The anomalistic year rides THIS movement's own apsidal line (the
     // scene's perihelion) — the factory owns the construction; the planet
@@ -619,6 +472,8 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
       periOfDateDegAt: (year) => sampleAt(year).periOfDateDeg,
       /** The hybrid's obliquity at a decimal year, degrees — the banked series inside its span, the α(t)-coupled ζ-tail integration beyond (the deep sampler grows ~0.1 s/Myr). @param {number} year @returns {number} */
       epsAt: (year) => sampleAt(year).epsDeg,
+      /** The hybrid's own domain, [yearFrom, yearTo] — the banked ζ series' span (plan 12 §8.5); null without a series. */
+      obliquityHybridDomainYears,
       /** The precession of date's two factors at a decimal year (the factory's own — ONE home). @param {number} year @returns {{solarTorque: number, ellipticity: number}} */
       precessionOfDateFactorsAt: (year) => hist.precessionOfDateFactorsAt(year - 2000),
       /** The spin axis's INERTIAL precession period at J2000, years — the self-anchored
@@ -686,8 +541,9 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   const moonSiderealMonthDaysAt = (year) => {
     const tMa = yearToTMa(year);
     const lod = deepLod.lodSecondsAtAge(tMa);
-    if (lod === null) return NaN;
-    const aAppKm = moonDistanceMetresAtAge(tMa) / 1000;
+    const aAppM = moonDistanceMetresAtAge(tMa);
+    if (lod === null || aAppM === null) return NaN;
+    const aAppKm = aAppM / 1000;
     const sidYrDays = deepLod.siderealYearSecondsAtAge(tMa) / lod;
     const deltaA = aAppKm * (1 / (MASS_RATIO_EARTH_MOON + 1)) * (moonSiderealMonthInput / sidYrDays);
     const aCorrM = (aAppKm + deltaA) * 1000;
@@ -950,11 +806,7 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   // (perigee advances, node regresses), T_p the J2000 precession anchor.
   // The retired era device counted these as eighth-integers per 8H and read
   // the periods back as totalDaysInH/N (the ∓13 was H/13 cycles per H).
-  const moonApsidalOfDateJ2000Seconds = C.moonReference.moonApsidalPrecessionDaysInputICRF * 86400;
-  const moonNodalOfDateJ2000Seconds = C.moonReference.moonNodalPrecessionDaysInputICRF * 86400;
-  const moonApsidalJ2000Seconds = 1 / (1 / moonApsidalOfDateJ2000Seconds - 1 / (axialPrecessionJ2000Days * 86400));
-  const moonNodalJ2000Seconds = 1 / (1 / moonNodalOfDateJ2000Seconds + 1 / (axialPrecessionJ2000Days * 86400));
-  const moonSiderealMonthJ2000Seconds = moonSiderealMonth * 86400;
+  const { moonApsidalJ2000Seconds } = chain;
   // (Plan 09 step 5c: the per-H of-date cycle counts and the month chain's H²
   // scaling of the of-date periods are deleted; the of-date perigee/node are
   // the chain's count compositions.)
@@ -964,43 +816,9 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   const moonAnomalisticMonthDays = 1 / (1 / moonSiderealMonth - 86400 / moonApsidalJ2000Seconds);
   const moonSynodicMonthDays = moonSiderealMonth * meanSiderealYearDays / (meanSiderealYearDays - moonSiderealMonth);
 
-  // moonEcc — the model's ONE eccentricity law — is created above (with
-  // eccentricityAt); the Moon channel's E-factor e(t)/e(J2000) and its
-  // perigee/node T² channel ride it, e(J2000) being the observed value exactly.
-
-  // Layer-2 month/precession chain (Brouwer-Clemence m² scaling × the
-  // e_E-line modulation)
-  const moonChain = createMoonMonthChain({
-    constants: {
-      aMoonNowMetres: EPOCH_PARAMS.moonDistanceNowM,
-      alpha1PerMa: EPOCH_PARAMS.alpha1PerMa,
-      alpha3PerMa3: EPOCH_PARAMS.alpha3PerMa3,
-      alpha4PerMa4: EPOCH_PARAMS.alpha4PerMa4,
-      gmEarthMoonM3PerS2: EPOCH_PARAMS.gmEarthMoonM3S2,
-      massRatioEarthMoon: MASS_RATIO_EARTH_MOON,
-      moonSiderealMonthInputDays: moonSiderealMonthInput,
-      meanSiderealYearJ2000Seconds: meanSiderealYearSeconds,
-      moonApsidalJ2000Seconds,
-      moonNodalJ2000Seconds,
-      moonSiderealMonthJ2000Seconds,
-      sPerigee: MOON_ECC_SENSITIVITY_PERIGEE,
-      sNode: MOON_ECC_SENSITIVITY_NODE,
-    },
-    fns: {
-      meanLodSecondsAtAge: /** @param {number} tMa */ (tMa) => deepLod.lodSecondsAtAge(tMa),
-      meanSiderealYearSecondsAtAge: /** @param {number} tMa */ (tMa) => deepLod.siderealYearSecondsAtAge(tMa),
-      // plan 09 step 5c: the tropical month's equinox turn is one per T_p(t)
-      precessionPeriodYearsAtAge: /** @param {number} tMa */ (tMa) => deepLod.lunisolarPrecessionPeriodYearsAtAge(tMa),
-      // Decision (ii) (owner-approved MODEL CHANGE, plan 02 §8): the whole
-      // lunar chain — this modulation, the integrated cycle counts it
-      // drives, the arguments' eccAt/channelIntegral and the E-factor —
-      // reads the ONE deep e. The era shift in the ancient eclipse stack is
-      // accepted and rebaselined with its explanation (deep e is
-      // La2004-corroborated in-era).
-      modulation: /** @param {number} tMa @param {number} s */ (tMa, s) => deepEcc.modulation(tMa, s),
-      distanceMetresAtAge: moonDistanceMetresAtAge,
-    },
-  });
+  // The Layer-2 month/precession chain (Brouwer-Clemence m² scaling × the
+  // e_E-line modulation on the ONE deep e — decision (ii), plan 02 §8) is
+  // the deep-time chain's `moonChain` above.
 
   // Chain-cycle integrator. S5/S12 conventions: age anchor = startmodelYear
   // (the scene's t_Ma convention), grid anchor C(2000) = 0 — grid anchor ≠
@@ -1568,7 +1386,7 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
       // + the calibrated ΔT cycle stack. 86,400.001780 s at J2000.
       lodRealSecondsAtYear: /** @param {number} year @returns {number} */ (year) => lodRealSeconds(year),
       alphaAtYear: /** @param {number} year @returns {number} */ (year) => earthMoiFactorAtAge(yearToTMa(year)),
-      moonDistanceKmAtYear: /** @param {number} year @returns {number} */ (year) => moonDistanceMetresAtAge(yearToTMa(year)) / 1000,
+      moonDistanceKmAtYear: /** @param {number} year @returns {number} */ (year) => moonDistanceKmAtYearOrNaN(year),
       // D6: OF-DATE — mass-loss law / the banked planetary λ̇ ratio
       // (identical at J2000 where the ratio ≡ 1 by construction)
       siderealYearSecondsAtYear: /** @param {number} year @returns {number} */ (year) => yearLengthsM.siderealYearSecondsAtYear(year),
@@ -1601,10 +1419,18 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
     // ±26 kyr) and T_peri/T_p with T_peri = 1/(1/T_p + 1/T_aps) the
     // perihelion-of-date period (0.812 at J2000) — on the of-date T_p.
     lunisolar: Object.freeze({
-      /** The mean lunisolar precession period T_p(t), years — the composed torque rate's period on the derived J2000 anchor (25,771.4 at J2000). @param {number} year @returns {number} */
+      /** The mean lunisolar precession period T_p(t), years — 1,296,000/the precession RATE of date (25,771.6 at J2000): THE published deep-time clock. @param {number} year @returns {number} */
       meanPeriodYearsAtYear: (year) => { const p = deepLod.lunisolarPrecessionPeriodYearsAtAge(yearToTMa(year)); return p === null ? certifiedAxialPrecessionJ2000Years() : p; },
-      /** The composed rate ψ̇(t) = [ω/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M)³], ″/yr (p₀ = 50.2883, derived). @param {number} year @returns {number} */
+      /** THE PRECESSION RATE OF DATE ψ̇(t) = [ω/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M)³] × cos ε̄(t)/cos ε̄(0), ″/yr — the constant's scaling on the spin and the lunar torque times the cited mean obliquity's cosine (plan 12 §8 option (a); p₀ = 50.2880, the IAU 2006 input). @param {number} year @returns {number} */
       meanRateArcsecPerYrAtYear: (year) => { const r = deepLod.lunisolarPrecessionRateArcsecPerYrAtAge(yearToTMa(year)); return r === null ? 1296000 / certifiedAxialPrecessionJ2000Years() : r; },
+      /** THE PRECESSION CONSTANT OF DATE as a rate on the J2000 obliquity, ″/yr — the spin and torque scaling alone (the hybrid obliquity integration's integrand: its own ŝ·n̂ carries the cosine). @param {number} year @returns {number} */
+      torqueRateArcsecPerYrAtYear: (year) => { const r = deepLod.lunisolarTorqueRateArcsecPerYrAtAge(yearToTMa(year)); return r === null ? 1296000 / certifiedAxialPrecessionJ2000Years() : r; },
+      /** 1,296,000/the torque rate, years. @param {number} year @returns {number} */
+      torquePeriodYearsAtYear: (year) => { const p = deepLod.lunisolarTorquePeriodYearsAtAge(yearToTMa(year)); return p === null ? certifiedAxialPrecessionJ2000Years() : p; },
+      /** cos ε̄(t)/cos ε̄(0) — the rate of date's obliquity factor on the constant (1 at J2000; 1.0075 at the Devonian, 1.0245 at 2.46 Ga). @param {number} year @returns {number} */
+      precessionObliquityFactorAtYear: (year) => { const f = deepLod.precessionObliquityFactorAtAge(yearToTMa(year)); return f === null ? 1 : f; },
+      /** THE CITED MEAN OBLIQUITY OF DATE ε̄(t), degrees — Farhat et al. 2022's own column of the Earth–Moon history (a Myr-scale MEAN: 23.25° at age 0 against the J2000 instant 23.44°; 22.23° Devonian, 19.73° at 2.46 Ga); read as the shape cos ε̄(t)/cos ε̄(0) on the precession rate. NaN beyond the input's span (4,429 Ma); forward of J2000 the J2000 slope is held. @param {number} year @returns {number} */
+      meanObliquityCitedDegAtYear: (year) => { const e = chain.meanObliquityCitedDegAtAge(yearToTMa(year)); return e === null ? NaN : e; },
       /** The of-date period — the one-family beat inside ±2 Myr, the composed mean beyond. @param {number} year @returns {number} */
       ofDatePeriodYearsAtYear: publishedAxialPrecessionYearsAtYear,
       /** f_S, the solar fraction of the J2000 precession torque. */
@@ -1697,7 +1523,15 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
       // fitted 16-harmonic K comb that sat beside it as `obliquityCombDeg`
       // (a device fitted to the scene's own wheel geometry, 1.27° off
       // La2004 at −20 kyr) is deleted with the frozen era device (plan 09).
+      // Beyond the ζ series' span (±10 Myr) the hybrid carries NO tidal evolution
+      // of the mean obliquity (nothing in ds/dt = α(ŝ·n̂)(ŝ×n̂) moves ŝ·n̂'s mean):
+      // its deep mean is the Quaternary mean on the tail's mean plane (plan 12
+      // §8.3). The cited mean obliquity of date is `meanObliquityCitedDeg`.
       obliquityDeg: /** @param {number} year @returns {number} */ (year) => oneSourceM.epsAt(year),
+      /** The hybrid's own domain, [yearFrom, yearTo] — the banked ζ series' span, ±10 Myr about J2000 (plan 12 §8.5: the published surfaces serve `obliquityDeg` inside it only; null without a series, i.e. a bare createModel()). */
+      obliquityHybridDomainYears: oneSourceM.obliquityHybridDomainYears,
+      /** The cited mean obliquity of date ε̄(t) (Farhat et al. 2022's column; plan 12 §8 option (a)), degrees — a Myr-scale mean; NaN beyond 4,429 Ma. @param {number} year @returns {number} */
+      meanObliquityCitedDeg: /** @param {number} year @returns {number} */ (year) => { const e = chain.meanObliquityCitedDegAtAge(yearToTMa(year)); return e === null ? NaN : e; },
       eccentricity: eccentricityAt,
       inclinationDeg,
       ascendingNodeDeg,
@@ -1746,7 +1580,7 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
     // the single surface the panel, charts, API and website consume.
     yearLengths: yearLengthsM,
     moon: Object.freeze({
-      distanceKmAtYear: /** @param {number} year @returns {number} */ (year) => moonDistanceMetresAtAge(yearToTMa(year)) / 1000,
+      distanceKmAtYear: /** @param {number} year @returns {number} */ (year) => moonDistanceKmAtYearOrNaN(year),
       siderealMonthDaysAtYear: moonSiderealMonthDaysAt,
       synodicMonthDays: moonSynodicMonthDays,
       // The apparent-position chain (truncated Meeus Ch. 47 series on
@@ -1780,7 +1614,7 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
       // construction sites (tools/verify/eclipse-audit.js, the browser
       // _eclipse twins) spread the SAME assembly into their own
       // createEclipseFinders call instead of triplicating it (the
-      // three-runtimes rule — cf. recession-history):
+      // three-runtimes rule — cf. moon-distance-history):
       frameworkSunDeps: Object.freeze({
         sunMeanLongitudeJ2000Deg: sunL0Deg,
         tropicalRateDegPerCy: sunTropicalRateDegPerCy,

package/src/moon/arguments.cjs

@@ -38,37 +38,11 @@
 
 'use strict';
 
-/**
- * Astronomical JD → decimal year: Julian calendar before 1582-10-15,
- * Gregorian after (Meeus Ch. 7 inverse). The S3-aligned coordinate for the
- * Sun secular deviations. @param {number} jd @returns {number}
- */
-function jdToDecimalYear(jd) {
-  const J = jd + 0.5;
-  const Z = Math.floor(J);
-  const F = J - Z;
-  let A = Z;
-  if (Z >= 2299161) {
-    const alpha = Math.floor((Z - 1867216.25) / 36524.25);
-    A = Z + 1 + alpha - Math.floor(alpha / 4);
-  }
-  const B = A + 1524;
-  const C = Math.floor((B - 122.1) / 365.25);
-  const D = Math.floor(365.25 * C);
-  const E = Math.floor((B - D) / 30.6001);
-  const day = B - D - Math.floor(30.6001 * E) + F;
-  const month = (E < 14) ? E - 1 : E - 13;
-  const year = (month > 2) ? C - 4716 : C - 4715;
-  const isGregorian = (Z >= 2299161);
-  const isLeap = isGregorian
-    ? (year % 4 === 0 && (year % 100 !== 0 || year % 400 === 0))
-    : (year % 4 === 0);
-  const monthLengths = [31, isLeap ? 29 : 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31];
-  let dayOfYear = day;
-  for (let m = 0; m < month - 1; m++) dayOfYear += monthLengths[m];
-  const daysInYear = isLeap ? 366 : 365;
-  return year + (dayOfYear - 1) / daysInYear;
-}
+// Astronomical JD → decimal year (Meeus Ch. 7 inverse, the 1582-10-15
+// reform): the S3-aligned coordinate for the Sun secular deviations. ONE
+// home since plan 13 B1 — time/calendar.cjs; re-exported here for the
+// readers that import it from the arguments module.
+const { jdToDecimalYear } = require('../time/calendar.cjs');
 
 /**
  * @typedef {{Lp: number, D: number, M: number, Mp: number, F: number}} MoonArgsDeg