@essrt/physics 4.0.0 → 4.1.0

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@essrt/physics",
-  "version": "4.0.0",
+  "version": "4.1.0",
   "description": "ESSRT (Expanding Solar System Resonance Theory) physics core \u2014 pure, injectable-constants, no I/O, no globals, no DOM. The complete model incl. fitted coefficients; every published version immutably ships one recorded model identity.",
   "essrt": {
     "modelVersion": "v13.0",
@@ -32,6 +32,9 @@
     "./phase": "./src/phase/index.cjs",
     "./cardinal": "./src/cardinal/index.cjs",
     "./moon/ecc-channel": "./src/moon/ecc-channel.cjs",
+    "./moon/deep-ecc-channel": "./src/moon/deep-ecc-channel.cjs",
+    "./moon/deep-modes-artifact": "./src/moon/deep-modes-artifact.cjs",
+    "./earth/deep-orbital-history": "./src/earth/deep-orbital-history.cjs",
     "./moon/month-chain": "./src/moon/month-chain.cjs",
     "./chain-cycles": "./src/chain-cycles/index.cjs",
     "./moon/arguments": "./src/moon/arguments.cjs",

package/src/constants/generated.d.ts

@@ -741,4 +741,13 @@ export declare const REFERENCE_DATA: {
       "max": number;
     };
   };
+  readonly planetSpinObserved: {
+    "marsSpinPrecessionArcsecPerYr": number;
+    "jupiterSpinPrecessionApproxArcsecPerYr": number;
+    "saturnSpinPrecessionPresentArcsecPerYr": number;
+    "saturnPresentToLongTermFraction": number;
+    "marsObliquityJ2000Deg": number;
+    "jupiterObliquityJ2000Deg": number;
+    "saturnObliquityJ2000Deg": number;
+  };
 };

package/src/constants/generated.js

@@ -14,6 +14,7 @@
  *   moonGrailWilliams2014            target
  *   giaCoxChaoPeltier                target
  *   obliquityChapront2002            target
+ *   planetSpinObserved               target
  *   ascendingNodesSouamiSouchay      target
  *   laplaceLagrangeBounds            target
  *   jplEclipticInclinationTrends     target
@@ -913,4 +914,13 @@ export const REFERENCE_DATA = Object.freeze({
       "max": 200
     }
   },
+  planetSpinObserved: {
+    "marsSpinPrecessionArcsecPerYr": -7.606,
+    "jupiterSpinPrecessionApproxArcsecPerYr": -2.8,
+    "saturnSpinPrecessionPresentArcsecPerYr": -0.45,
+    "saturnPresentToLongTermFraction": 0.68,
+    "marsObliquityJ2000Deg": 25.19,
+    "jupiterObliquityJ2000Deg": 3.13,
+    "saturnObliquityJ2000Deg": 26.73
+  },
 });

package/src/earth/deep-orbital-history.cjs

@@ -0,0 +1,190 @@
+/**
+ * THE DEEP ORBITAL HISTORY FACTORY — the ONE home for the Stage-C hybrid
+ * mathematics (plan 02 §8, C-3): Earth's deep-time orbital quantities from
+ * the engine's own governed mode tables, and the obliquity hybrid
+ *
+ *     ds/dt = α (ŝ·n̂)(ŝ×n̂),   α = ψ̇_J2000 / cos ε₀
+ *
+ * with the orbit normal n̂(t) from an anchored ζ mode-sum and ψ̇_J2000 the
+ * injected engine-K anchor (the equinox-precession rate measured from the
+ * engine's own year lengths). ZERO fitted constants; the J2000 celestial
+ * pole sits at ecliptic longitude +90° (s₀ = (0, +sin ε₀, cos ε₀) — the
+ * −90° choice anti-correlates perfectly, the recorded convention catch).
+ *
+ * Consumers (each injecting its own certified surfaces, cross-engine
+ * matched): the Stage-C lab / verdict generator (both ζ tiers), the Node
+ * shim tools/lib/deep-orbital-history.js (the deep insolation extractor),
+ * and the browser's published ε surfaces (VFP chart, Earth-panel row) via
+ * the generated deep-modes embed. The ERA machinery (frame conversions,
+ * lunar D5, besselian, the scene's A-solve and Step 6b) stays on the
+ * certified fitted law — measured: the hybrid differs from it by 5″ rms in
+ * 1900–2100 / 19.5″ in 1600–2400, mostly hybrid-side there — the same
+ * certification split as e and ϖ.
+ *
+ * The climatic precession is the PHYSICAL of-date angle from the moving
+ * equinox γ = unit(ŝ×n̂) to the perihelion direction, measured in the orbit
+ * plane about n̂ — no linear-precession approximation at any epoch
+ * (validated: reproduces the J2000 ecliptic longitude and La2004's
+ * e·sin ϖ̃ exactly at t = 0).
+ */
+
+'use strict';
+
+/**
+ * @param {{
+ *   zModes: ReadonlyArray<{omegaRadPerYr: number, re: number, im: number}>,
+ *   zetaModes: ReadonlyArray<{omegaRadPerYr: number, re: number, im: number}>,
+ *   anchorE: number,
+ *   anchorPeriEclipticDeg: number,
+ *   anchorInclEclipticDeg: number,
+ *   anchorAscNodeEclipticDeg: number,
+ *   axialPrecessionYearsJ2000: number,
+ *   obliquityJ2000Deg: number,
+ * }} deps — mode tables from the governed deep artifact (choose the ζ tier
+ *   per consumer: era = its own 8-term extraction, deep = the 16-term
+ *   table; NEVER a slice); anchors from the chain artifact's one home; the
+ *   axial-precession years and ε₀ from the injecting engine's certified
+ *   surfaces.
+ */
+function createDeepOrbitalHistory({
+  zModes, zetaModes,
+  anchorE, anchorPeriEclipticDeg, anchorInclEclipticDeg, anchorAscNodeEclipticDeg,
+  axialPrecessionYearsJ2000, obliquityJ2000Deg,
+}) {
+  const D2R = Math.PI / 180, R2D = 180 / Math.PI;
+
+  const mkAnchored = (
+    /** @type {ReadonlyArray<{omegaRadPerYr: number, re: number, im: number}>} */ modes,
+    /** @type {[number, number]} */ anchor,
+  ) => {
+    const sum = (/** @type {number} */ t) => {
+      let re = 0, im = 0;
+      for (const m of modes) {
+        const c = Math.cos(m.omegaRadPerYr * t), s = Math.sin(m.omegaRadPerYr * t);
+        re += m.re * c - m.im * s;
+        im += m.re * s + m.im * c;
+      }
+      return [re, im];
+    };
+    const s0 = sum(0), R = [anchor[0] - s0[0], anchor[1] - s0[1]];
+    return (/** @type {number} */ t) => { const [x, y] = sum(t); return [x + R[0], y + R[1]]; };
+  };
+  const s2h = Math.sin(anchorInclEclipticDeg / 2 * D2R);
+  const zetaAt = mkAnchored(zetaModes,
+    [s2h * Math.cos(anchorAscNodeEclipticDeg * D2R), s2h * Math.sin(anchorAscNodeEclipticDeg * D2R)]);
+  const zAt = mkAnchored(zModes,
+    [anchorE * Math.cos(anchorPeriEclipticDeg * D2R), anchorE * Math.sin(anchorPeriEclipticDeg * D2R)]);
+
+  /** @param {number} t */
+  const inclNode = (t) => {
+    const [q, p] = zetaAt(t);
+    return { i: 2 * Math.asin(Math.min(1, Math.hypot(q, p))), Om: Math.atan2(p, q) };
+  };
+  /** @param {number} t */
+  const orbitNormal = (t) => {
+    const { i, Om } = inclNode(t);
+    return [Math.sin(i) * Math.sin(Om), -Math.sin(i) * Math.cos(Om), Math.cos(i)];
+  };
+
+  const psiDot = (2 * Math.PI) / axialPrecessionYearsJ2000;
+  const EPS0 = obliquityJ2000Deg * D2R;
+  const ALPHA = psiDot / Math.cos(EPS0);
+
+  const cross = (/** @type {number[]} */ a, /** @type {number[]} */ b) =>
+    [a[1] * b[2] - a[2] * b[1], a[2] * b[0] - a[0] * b[2], a[0] * b[1] - a[1] * b[0]];
+  const dot = (/** @type {number[]} */ a, /** @type {number[]} */ b) => a[0] * b[0] + a[1] * b[1] + a[2] * b[2];
+  const deriv = (/** @type {number[]} */ s, /** @type {number} */ t) => {
+    const n = orbitNormal(t);
+    const k = ALPHA * dot(s, n);
+    const c = cross(s, n);
+    return [k * c[0], k * c[1], k * c[2]];
+  };
+  const rk4 = (/** @type {number[]} */ s, /** @type {number} */ t, /** @type {number} */ h) => {
+    const k1 = deriv(s, t);
+    const k2 = deriv([s[0] + h / 2 * k1[0], s[1] + h / 2 * k1[1], s[2] + h / 2 * k1[2]], t + h / 2);
+    const k3 = deriv([s[0] + h / 2 * k2[0], s[1] + h / 2 * k2[1], s[2] + h / 2 * k2[2]], t + h / 2);
+    const k4 = deriv([s[0] + h * k3[0], s[1] + h * k3[1], s[2] + h * k3[2]], t + h);
+    const o = [0, 1, 2].map((i) => s[i] + h / 6 * (k1[i] + 2 * k2[i] + 2 * k3[i] + k4[i]));
+    const r = Math.hypot(o[0], o[1], o[2]);
+    return [o[0] / r, o[1] / r, o[2] / r];
+  };
+
+  /** One quantity bundle at time t (years from J2000) from spin axis s. */
+  function sampleAt(/** @type {number[]} */ s, /** @type {number} */ t) {
+    const n = orbitNormal(t);
+    const { i, Om } = inclNode(t);
+    const [zx, zy] = zAt(t);
+    const e = Math.hypot(zx, zy);
+    const w = Math.atan2(zy, zx) - Om;
+    const ph = [
+      Math.cos(Om) * Math.cos(w) - Math.sin(Om) * Math.sin(w) * Math.cos(i),
+      Math.sin(Om) * Math.cos(w) + Math.cos(Om) * Math.sin(w) * Math.cos(i),
+      Math.sin(w) * Math.sin(i),
+    ];
+    const g = cross(s, n);
+    const gn = Math.hypot(g[0], g[1], g[2]);
+    const gu = [g[0] / gn, g[1] / gn, g[2] / gn];
+    const cosA = dot(gu, ph);
+    const sinA = dot(n, cross(gu, ph));
+    const periOfDateDeg = ((Math.atan2(sinA, cosA) * R2D) % 360 + 360) % 360;
+    return {
+      epsDeg: Math.acos(Math.max(-1, Math.min(1, dot(s, n)))) * R2D,
+      e,
+      periOfDateDeg,
+      eSinPeri: e * Math.sin(periOfDateDeg * D2R),
+      eCosPeri: e * Math.cos(periOfDateDeg * D2R),
+      inclEclDeg: i * R2D,
+    };
+  }
+
+  /**
+   * Integrate once over [min(t0,t1,0), max(t0,t1,0)] and return a sampler
+   * on a stepYr grid (linear interpolation; angles interpolated on the
+   * wrapped difference). t in years from J2000, negative = past.
+   * @param {number} t0Yr @param {number} t1Yr @param {number} stepYr
+   */
+  function build(t0Yr, t1Yr, stepYr) {
+    const tMin = Math.min(t0Yr, t1Yr, 0), tMax = Math.max(t0Yr, t1Yr, 0);
+    const grid = new Map();
+    const S0 = [0, Math.sin(EPS0), Math.cos(EPS0)];
+    const H_STEP = 5;
+    for (const dir of [-1, +1]) {
+      let s = S0, t = 0;
+      grid.set(0, sampleAt(s, 0));
+      const end = dir < 0 ? tMin : tMax;
+      while (dir < 0 ? t > end : t < end) {
+        s = rk4(s, t, dir * H_STEP);
+        t += dir * H_STEP;
+        // Store ONLY at true step multiples: an end-of-range sample under a
+        // ROUNDED key collides with a real grid node (measured during C-3:
+        // build(+200, …, 1000) stored ε(+200) under key 0, clobbering the
+        // t = 0 sample by 94″ — the verdict generator's rms caught it).
+        // Partial edge segments are covered by at()'s a-side fallback.
+        if (Math.abs(t % stepYr) < H_STEP / 2) grid.set(Math.round(t / stepYr) * stepYr, sampleAt(s, t));
+      }
+    }
+    return {
+      /** @param {number} tYr */
+      at(tYr) {
+        const k0 = Math.floor(tYr / stepYr) * stepYr, k1 = k0 + stepYr;
+        const a = grid.get(k0), b = grid.get(k1) ?? a;
+        if (!a) throw new Error(`deep-orbital-history: ${tYr} outside the built range`);
+        const f = (tYr - k0) / stepYr;
+        const lerp = (/** @type {number} */ x, /** @type {number} */ y) => x + (y - x) * f;
+        const dAng = (/** @type {number} */ x, /** @type {number} */ y) => x + (((y - x + 540) % 360) - 180) * f;
+        return {
+          epsDeg: lerp(a.epsDeg, b.epsDeg),
+          e: lerp(a.e, b.e),
+          periOfDateDeg: ((dAng(a.periOfDateDeg, b.periOfDateDeg) % 360) + 360) % 360,
+          eSinPeri: lerp(a.eSinPeri, b.eSinPeri),
+          eCosPeri: lerp(a.eCosPeri, b.eCosPeri),
+          inclEclDeg: lerp(a.inclEclDeg, b.inclEclDeg),
+        };
+      },
+    };
+  }
+
+  return { build, alphaArcsecPerYr: ALPHA * R2D * 3600 };
+}
+
+module.exports = { createDeepOrbitalHistory };

package/src/index.js

@@ -121,6 +121,13 @@ export * as planetOrbitChain from './planets/orbit-chain.cjs';
 // governed artifact; the browser flag path consumes exactly these).
 export { buildPlanetChainsFromArtifactData, computePlanetElementsAtYear, computeHeliocentricEclipticFromElements, computePoissonArgRad, computeOsculatingElements, solveKeplerRad, ANCHOR_EPOCH_YEAR, ANCHOR_EPOCH_JD } from './planets/keplerian-chain.cjs';
 export { CHAIN_ARTIFACT, CHAIN_ARTIFACT_HASH } from './planets/chain-artifact.js';
+// Engine-switch Stage B, T5d-(d): the deep-time Earth-z channel (the deep
+// lunar-chain modulation rides it; era consumers stay on the H/3 channel).
+export { createDeepEccChannel } from './moon/deep-ecc-channel.cjs';
+export { DEEP_MODES_ARTIFACT, DEEP_MODES_ARTIFACT_HASH } from './moon/deep-modes-artifact.cjs';
+// Stage C-3: the ONE home of the obliquity-hybrid / deep-orbital-history
+// mathematics (published-ε surfaces + the deep insolation features).
+export { createDeepOrbitalHistory } from './earth/deep-orbital-history.cjs';
 export { createPredictivePrecession, calcPlanetPerihelionLongDeg } from './planets/predict.cjs';
 // L10 — the composition front door: one law set, N body records. Thin by
 // design; engines keep their direct call sites (see planets/model.cjs).

package/src/model.js

@@ -27,6 +27,8 @@ import { createDeepTimeLod } from './deltat/deep-time.cjs';
 import { createMoonRecessionHistory, createSolarChannelBudget } from './deltat/recession-history.cjs';
 import { evalClimateL1OrbitalPermil } from './climate/l1-orbital.cjs';
 import { createMoonEccChannel } from './moon/ecc-channel.cjs';
+import { createDeepEccChannel } from './moon/deep-ecc-channel.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';
@@ -333,9 +335,13 @@ export function assembleModel(C, F, laws = {}) {
   // measured, the H/16 law's present ė (−0.84e-5/cy) is 5× below the
   // observed −4.20e-5 while this law reads −4.31e-5 (Phase-0 record:
   // tools/explore/fq7s-h3-law-candidate.mjs; JPL Sun 1.49″, registry 0.80″,
-  // syzygy 3.72″). Consumers: the eclipse Sun (equation of centre), the
-  // besselian Sun distance, the Moon channel (E-factor, perigee/node T²),
-  // and the cardinal-point braid — one eccentricity everywhere.
+  // syzygy 3.72″). Consumers (engine-switch decision (ii) restatement,
+  // plan 02 §8): the eclipse Sun (equation of centre), the besselian Sun
+  // distance and the cardinal-point braid — the CLOCK-side machinery. The
+  // LUNAR chain (E-factor, perigee/node modulation/T², argument Δe² and
+  // of-date rates) rides the ONE deep e (deep-ecc-channel.cjs — the
+  // engine's own ±10-Myr z-vector); this H/3 line is its epoch-local
+  // tangent (E18), agreeing within 4.2e-5 across the historical era.
   // ONE implementation for all three runtimes: moon/ecc-channel.cjs (the
   // Node engine's deep-time.js and the browser's script.js instantiate the
   // same channel with the same inputs). Its phase counter runs from J2000
@@ -349,6 +355,14 @@ export function assembleModel(C, F, laws = {}) {
     inclinationCycleAnchorDeg: C.earthOrbital.earthInclinationCycleAnchor,
     eccentricityJ2000: C.earthOrbital.earthEccentricityJ2000,
   });
+  // Engine-switch decision (ii) (plan 02 §8): the ONE deep e — the engine's
+  // own ±10-Myr mode table, anchored form — feeds the ENTIRE lunar chain
+  // (modulation, cycle counts, arguments eccAt/channelIntegral, E-factor).
+  // The Sun/clock machinery (eclipse Sun EoC, besselian Sun distance,
+  // cardinal braid) stays on the H/3 channel above — a certification
+  // split, not a physics one (the two agree within 4.2e-5 in-era). The
+  // laws hook deliberately does NOT reach this channel.
+  const deepEcc = createDeepEccChannel(DEEP_MODES_ARTIFACT);
   /** @param {number} year @returns {number} */
   const eccentricityAt = laws.eccentricityAt ?? ((year) => moonEcc.eccAt(year - 2000));
   /** de/dyear of the one law — the cardinal braid's equation-of-centre
@@ -682,7 +696,13 @@ export function assembleModel(C, F, laws = {}) {
       meanLodSecondsAtAge: /** @param {number} tMa */ (tMa) => deepLod.lodSecondsAtAge(tMa),
       meanSiderealYearSecondsAtAge: /** @param {number} tMa */ (tMa) => deepLod.siderealYearSecondsAtAge(tMa),
       meanHAtAge: /** @param {number} tMa */ (tMa) => deepLod.hAtAge(tMa),
-      modulation: /** @param {number} tMa @param {number} s */ (tMa, s) => moonEcc.modulation(tMa, s),
+      // 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,
     },
   });
@@ -757,11 +777,12 @@ export function assembleModel(C, F, laws = {}) {
       eccentricityDotDotJ2000: C.earthOrbital.earthEccentricityDotDotJ2000,
       elpEarthFigureJ2ArcsecPerCy2: C.moonMeeus.elpW1T2Decomposition_arcsecPerCy2.earthFigureJ2,
       elpGeneralPrecessionPA_T2ArcsecPerCy2: C.moonMeeus.elpW1T2Decomposition_arcsecPerCy2.generalPrecessionPA_T2_Lieske1976,
-      eccE0: moonEcc.e0,
+      // Decision (ii): the lunar chain reads ONE deep e end to end
+      eccE0: deepEcc.e0,
     },
     fns: {
-      eccAt: /** @param {number} tYr */ (tYr) => moonEcc.eccAt(tYr),
-      channelIntegral: /** @param {number} T @param {number} s */ (T, s) => moonEcc.channelIntegral(T, s),
+      eccAt: /** @param {number} tYr */ (tYr) => deepEcc.eccAt(tYr),
+      channelIntegral: /** @param {number} T @param {number} s */ (T, s) => deepEcc.channelIntegral(T, s),
       computeObliquityEarth: obliquitySnapshotDeg,
       jdToSIyear: yearFromJD,
       tropicalOrbitsBetween: mcTropical,
@@ -809,10 +830,10 @@ export function assembleModel(C, F, laws = {}) {
     return Number.isFinite(dT) ? dT : 0;
   };
 
-  // Bounded Meeus E-factor from the H/3 line (framework-native branch only;
-  // the pure-Meeus polynomial A/B branch stays engine-local)
+  // Bounded Meeus E-factor from the ONE deep e (decision (ii); the
+  // pure-Meeus polynomial A/B branch stays engine-local)
   /** @param {number} dDays @returns {number} */
-  const fwEFactor = (dDays) => moonEcc.eFactorAt(dDays / C.foundational.inputmeanlengthsolaryearindays);
+  const fwEFactor = (dDays) => deepEcc.eFactorAt(dDays / C.foundational.inputmeanlengthsolaryearindays);
 
   // D2 derived additional-argument rates (deg/cy, J2000 8H-lattice months;
   // record: tools/explore/derive-a1a2a3.js)

package/src/moon/deep-ecc-channel.cjs

@@ -0,0 +1,132 @@
+/**
+ * The DEEP Earth-eccentricity channel — engine-switch Stage B, the T5d-(d)
+ * decision (ii) (plan 02 §8, owner-approved MODEL CHANGE, pre-registered
+ * acceptance in the same record).
+ *
+ * Evaluates Earth's e(t) from the engine's own ±10-Myr NAFF mode table
+ * (data/nbody-deep-secular-modes.json, embedded verbatim by generate.mjs as
+ * deep-modes-artifact.cjs) in the ANCHORED form the era chain also uses:
+ *
+ *     z(t) = Σ (re + i·im)·e^{iωt} + (z_J2000 − Σ(0)),   e = |z|
+ *
+ * — the constant remainder pins the unresolved short/ultra-long content at
+ * the JPL J2000 seed anchor, so e(J2000) is exact by construction.
+ *
+ * DOMAIN (decision (ii), proper physics — no mixed device): the ENTIRE
+ * lunar chain reads THIS channel — eccAt (the Δe² argument terms),
+ * channelIntegral (perigee/node of-date rates), modulation (the
+ * month/precession chain and its integrated cycle counts), eFactorAt (the
+ * Meeus E-factor). The Sun/clock machinery (eclipse Sun equation of
+ * centre, besselian Sun distance, cardinal braid) stays on the H/3 law
+ * (moon/ecc-channel.cjs) — a CERTIFICATION split, not a physics one: the
+ * two e's agree within 4.2e-5 wherever the domains overlap, and the H/3
+ * line keeps its E18 role as the epoch-local tangent inside the fitted
+ * clock stack.
+ *
+ * MEASURED at adoption (banked in the plan record): vs the H/3 law the
+ * attractor-mean modulation shifts are +2,526 ppm (perigee channel,
+ * s = 2.407) and +1,067 ppm (node, s = 1.018) — month/LOD shifts of
+ * ~10–25 ppm class, all 41 paleo anchors in bands; the ancient eclipse
+ * stack's cycle counts shift (χ² IMPROVES — recorded with its
+ * explanation: this e is La2004-corroborated in-era, 3× closer than the
+ * H/3 line; theory-corroborated, never tuned). Beyond ~5 Myr the table is
+ * spectrum-class (phases decohere) — the quasi-periodic ATTRACTOR (mean e
+ * and the 405-kyr-class beat structure) is the claim, never a pointwise
+ * ephemeris.
+ *
+ * LOAD-BEARING conventions (carried VERBATIM from the era channel —
+ * moon/ecc-channel.cjs — so the certified numerics match):
+ *  - channelIntegral uses composite Simpson with step ≤ ~4,000 yr
+ *    (N ≥ 2, even);
+ *  - modulation(tMa, s) takes age in Myr, positive = past, ≡ 1 at J2000;
+ *  - eFactorAt(tYr) = eccAt(tYr) / e0, ≡ 1 at J2000.
+ */
+
+'use strict';
+
+/**
+ * @param {{
+ *   earthZ: ReadonlyArray<{omegaRadPerYr: number, re: number, im: number}>,
+ *   anchorE: number,
+ *   anchorPeriEclipticDeg: number,
+ * }} artifact — the embedded deep-modes artifact (earth z-modes + the J2000
+ *   anchor pair joined from the chain artifact's one home).
+ */
+function createDeepEccChannel({ earthZ, anchorE, anchorPeriEclipticDeg }) {
+  const D2R = Math.PI / 180;
+  const zA = [anchorE * Math.cos(anchorPeriEclipticDeg * D2R), anchorE * Math.sin(anchorPeriEclipticDeg * D2R)];
+  /** @param {number} tYr @returns {[number, number]} */
+  const modeSum = (tYr) => {
+    let re = 0, im = 0;
+    for (const m of earthZ) {
+      const c = Math.cos(m.omegaRadPerYr * tYr), s = Math.sin(m.omegaRadPerYr * tYr);
+      re += m.re * c - m.im * s;
+      im += m.re * s + m.im * c;
+    }
+    return [re, im];
+  };
+  const s0 = modeSum(0);
+  const R = [zA[0] - s0[0], zA[1] - s0[1]];   // the constant remainder
+
+  /** Earth's e at tYr years from J2000 (negative = past).
+   *  @param {number} tYr @returns {number} */
+  function eccAt(tYr) {
+    const [x, y] = modeSum(tYr);
+    return Math.hypot(x + R[0], y + R[1]);
+  }
+
+  /** de/dyear — analytic from the mode sum: d|z|/dt = (z·ż)/|z|.
+   *  @param {number} tYr @returns {number} */
+  function eccRateAt(tYr) {
+    const [x0, y0] = modeSum(tYr);
+    const x = x0 + R[0], y = y0 + R[1];
+    let dx = 0, dy = 0;
+    for (const m of earthZ) {
+      const c = Math.cos(m.omegaRadPerYr * tYr), s = Math.sin(m.omegaRadPerYr * tYr);
+      dx += m.omegaRadPerYr * (-m.re * s - m.im * c);
+      dy += m.omegaRadPerYr * (m.re * c - m.im * s);
+    }
+    const r = Math.hypot(x, y);
+    return r === 0 ? 0 : (x * dx + y * dy) / r;
+  }
+
+  const e0 = eccAt(0);                          // ≡ anchorE by construction
+  const g0 = Math.pow(1 - e0 * e0, -1.5);
+
+  /** [g(e(t))/g₀]^s — the perigee/node rate modulation, ≡ 1 at J2000.
+   *  tMa is age in Myr (positive = past, the deep-time chain convention).
+   *  @param {number} tMa @param {number} s @returns {number} */
+  function modulation(tMa, s) {
+    if (tMa === 0) return 1;
+    const e = eccAt(-tMa * 1e6);
+    return Math.pow(Math.pow(1 - e * e, -1.5) / g0, s);
+  }
+
+  /** ∫₀ᵀ [(g(e(t′))/g₀)^s − 1] dt′ in Julian centuries — the phase-aware
+   *  channel-rate integral. Composite Simpson, step ≤ ~4 kyr (the era
+   *  channel's certified discretization, verbatim).
+   *  @param {number} T @param {number} s @returns {number} */
+  function channelIntegral(T, s) {
+    if (T === 0) return 0;
+    /** @param {number} t */
+    const f = (t) => {
+      const e = eccAt(t * 100);   // t in cy → years
+      return Math.pow(Math.pow(1 - e * e, -1.5) / g0, s) - 1;
+    };
+    const N = Math.max(2, 2 * Math.ceil(Math.abs(T) * 100 / 8000));
+    const h = T / N;
+    let sum = f(0) + f(T);
+    for (let i = 1; i < N; i++) sum += f(i * h) * (i % 2 ? 4 : 2);
+    return sum * h / 3;
+  }
+
+  /** Bounded Meeus E-factor: E ≡ e(t)/e(J2000), ≡ 1 at J2000.
+   *  @param {number} tYr @returns {number} */
+  function eFactorAt(tYr) {
+    return eccAt(tYr) / e0;
+  }
+
+  return { eccAt, eccRateAt, modulation, channelIntegral, eFactorAt, e0, g0 };
+}
+
+module.exports = { createDeepEccChannel };

package/src/moon/deep-modes-artifact.cjs

@@ -0,0 +1,19 @@
+/**
+ * GENERATED — do not edit. Regenerate:
+ *   node tools/constants/generate.mjs --write
+ *
+ * Source: data/nbody-deep-secular-modes.json (the governed deep-time
+ * Earth-z mode table, engine-switch Stage B), earth z-modes + verdict
+ * emitted VERBATIM, plus the J2000 anchor pair (e, ϖ) joined from
+ * data/nbody-secular-frequencies.json j2000AnchorElements.earth — the
+ * anchor's ONE home. Two gates guard the chain: check:artifacts pins
+ * artifact ↔ engine; generate.mjs check mode pins this embed ↔ artifact.
+ * CJS because the lunar-chain consumers are CommonJS modules.
+ */
+'use strict';
+
+const DEEP_MODES_ARTIFACT_HASH = "b74be2d9b00a3460";
+
+const DEEP_MODES_ARTIFACT = Object.freeze({"meta":{"integrator":"wh","order":2,"dtDays":2,"spanYears":20000000,"sampleDays":20000,"gr":true,"seed":"JPL Horizons J2000 heliocentric state vectors (tools/explore/j2000-state.mjs, the one home)","masses":"DE440 mass ratios (astro-reference physicalConstants)","runCommand":"node tools/explore/lattice-long-window-test.mjs years=20000000 integrator=wh dt=2 order=2 gr=1 frame=both sample=20000","dumpFile":"tools/explore/lattice-long-window-ecliptic-20000000-gr.local.json (untracked, 337 MB)","dumpSha256":"c6f37b93f8f53c5b72d9e2935085145820d384ed61ecdf69f3831895d5f0a516","conservationMaxDE":1.3523274883097007e-8,"decimation":4,"naffTerms":18,"naffZetaTerms":16,"naffZetaTermsEra":8,"frame":"ecliptic-J2000 (z-modes for the lunar-chain e; ζ-modes for the Stage-C obliquity hybrid — the invariable-frame ζ table remains unbanked, no consumer)","laskarRef":"Laskar, J. et al. (2004), A&A 428, 261–285, Table 3 (theory reference labels)"},"verdict":{"strongestBeatPeriodKyr":409.398962153882,"strongestBeatArcsecPerYr":3.165616232101899,"la2004BeatPeriodKyr":405.6972922210049,"rockMetronomeKyr":405.6,"companion124Kyr":123.63439229002236,"companion95Kyr":94.95801992141148,"earthLeadingModesArcsecPerYr":[4.257339363798321,7.422955595900221,17.905475677851488,17.35600814229943],"topBeatLines":[{"periodKyr":409.398962153882,"productAmp":0.0003189871275662614},{"periodKyr":94.95801992141148,"productAmp":0.00023593726318903166},{"periodKyr":123.63439229002236,"productAmp":0.0002021269449318684},{"periodKyr":98.94135212634195,"productAmp":0.0001820123713850866},{"periodKyr":130.4734867701679,"productAmp":0.00015592960633097004},{"periodKyr":2358.6470831218417,"productAmp":0.00011533256795834209},{"periodKyr":974.251698194696,"productAmp":0.00007391281237249426},{"periodKyr":706.1267630803258,"productAmp":0.00006332094707822819}],"note":"g5 matches La2004 to 0.0002 ″/yr; the beat gap vs 405.6 is g2 (chaotic-diffusion window sensitivity + omitted separate-Moon/asteroid terms) — the anatomy is the plan’s Stage-B record."},"earthZ":[{"omegaRadPerYr":0.000020640163686955943,"re":0.01663782532735386,"im":0.009773828910268845},{"omegaRadPerYr":0.00003598750427161015,"re":-0.0155511555977695,"im":-0.005606883657279738},{"omegaRadPerYr":0.00008680819575396443,"re":-0.011017310146969267,"im":0.005302931248387234},{"omegaRadPerYr":0.00008414430196835267,"re":0.005902460235142885,"im":-0.0073575521441138925},{"omegaRadPerYr":0.000027089406032463197,"re":0.0003820712762735637,"im":0.00381132243985138},{"omegaRadPerYr":0.00008287265957255809,"re":0.0018740756883102302,"im":0.0007570047078217853},{"omegaRadPerYr":0.00008556055947193001,"re":0.00065046714460136,"im":0.001636722269192129},{"omegaRadPerYr":0.00003343872861356042,"re":-0.0015906877043739744,"im":0.0008131867314299489},{"omegaRadPerYr":0.00008733249835491452,"re":0.0019695845162428804,"im":0.0007638996859228755},{"omegaRadPerYr":0.00008624160326238742,"re":-0.0007912287558015328,"im":0.0008454083667608703},{"omegaRadPerYr":0.00013693637603446235,"re":-0.00091728504412424,"im":0.0011947455784994275},{"omegaRadPerYr":0.00008149221104720623,"re":-0.0002919323858022318,"im":0.0013378769276789224},{"omegaRadPerYr":0.000035340366616489944,"re":-0.0014017009332458236,"im":-0.0005418283702703829},{"omegaRadPerYr":0.00003411910644328693,"re":0.0009644062512081837,"im":-0.0009278208477004621},{"omegaRadPerYr":0.000026443547751027787,"re":-0.0003484265086069058,"im":-0.0019049634588757702},{"omegaRadPerYr":0.00008946536777652535,"re":-0.0010363626781434278,"im":-0.00021237189770655706},{"omegaRadPerYr":0.000027681653760243513,"re":0.00026309559235420505,"im":0.0013827244513428105},{"omegaRadPerYr":0.00003665339440166644,"re":0.001341581161662452,"im":0.00035305205249495507}],"earthZeta":[{"omegaRadPerYr":-2.96732779276946e-14,"re":-0.004153530797421736,"im":0.013136287139943213},{"omegaRadPerYr":-0.00009139185627510659,"re":-0.0039341668744113665,"im":-0.007077440344808186},{"omegaRadPerYr":-0.000027083329194200705,"re":0.0041703974742109,"im":0.0016342020564862083},{"omegaRadPerYr":-0.00008606223132085272,"re":0.0015441443378585364,"im":-0.0036871518361985415},{"omegaRadPerYr":-0.00003413935334637193,"re":-0.002956817542124509,"im":0.0024359440436811054},{"omegaRadPerYr":-0.00008872604497912732,"re":-0.00013216158559522037,"im":-0.002797304033559708},{"omegaRadPerYr":-0.00003326020694677573,"re":0.0008627085184454156,"im":-0.0008468691172706413},{"omegaRadPerYr":-0.00003470927143831441,"re":0.0013226438289662852,"im":-0.0011121945465339342},{"omegaRadPerYr":-0.00009406049953471362,"re":0.001297329508193755,"im":0.0008926331980509192},{"omegaRadPerYr":-0.00002648259433928881,"re":0.001695647505347929,"im":0.0006707872920285527},{"omegaRadPerYr":-0.00012773804995598164,"re":-0.0008156188441322624,"im":0.0010583300732716742},{"omegaRadPerYr":-0.00003267570575865712,"re":0.00028878599128559593,"im":-0.0004498471691632459},{"omegaRadPerYr":-0.000033725843136033775,"re":0.0012152922715564112,"im":-0.0007872735440399818},{"omegaRadPerYr":-0.000014507526525427572,"re":-0.00064166716018201,"im":0.0005786259131500737},{"omegaRadPerYr":-0.00009265988888533535,"re":-0.0005654435501795496,"im":0.00014913577755089046},{"omegaRadPerYr":-0.000003353112962211696,"re":0.0006060116994915085,"im":0.00025189133351366343}],"earthZetaEra":[{"omegaRadPerYr":-2.96732779276946e-14,"re":-0.0041400895842324354,"im":0.01314865914923023},{"omegaRadPerYr":-0.00009139185627510659,"re":-0.0038874970989988566,"im":-0.0070425256250290546},{"omegaRadPerYr":-0.000027083329194200705,"re":0.004075700185805294,"im":0.0016137774467260707},{"omegaRadPerYr":-0.00008606223132085272,"re":0.0015261120191598105,"im":-0.0036992031256430353},{"omegaRadPerYr":-0.00003413935334637193,"re":-0.0031492966403078486,"im":0.002570649693650614},{"omegaRadPerYr":-0.00008872604497912732,"re":-0.00014686941549255626,"im":-0.0027917154781614112},{"omegaRadPerYr":-0.00003326020694677573,"re":0.0005693304183459907,"im":-0.0006693278314130071},{"omegaRadPerYr":-0.00003470927143831441,"re":0.0012956912361530007,"im":-0.0010932966486614268}],"anchorE":0.016702435651957018,"anchorPeriEclipticDeg":462.9179213996577,"anchorInclEclipticDeg":0.00010345820355474172,"anchorAscNodeEclipticDeg":140.29217988414328});
+
+module.exports = { DEEP_MODES_ARTIFACT, DEEP_MODES_ARTIFACT_HASH };

package/src/moon/ecc-channel.cjs

@@ -8,14 +8,17 @@
  * max. The observed J2000 e (−0.86%), ė (+1.7%) and the sign of ë are
  * PREDICTIONS of this line, not inputs.
  *
- * ONE law for the eclipse chain (FQ-7-Sun option C-small, doc 66 §5): the
- * eclipse Sun's equation of centre rides this same H/3 line (anchored at
- * its J2000 value, model.js sunEccentricityAt) — eccentricity is
- * frame-invariant and may carry only fixed-frame lattice periods; the
- * H/16 perihelion cycle is the H/3 rotation seen from the H/13 equinox
- * (13 + 3 = 16) and belongs to ϖ, not e (doc 108). Only the
- * cardinal-point/deep-time path still keeps the H/16 law (`cardinal`'s
- * injected eccentricityAt) — a recorded, not yet resolved, split.
+ * DOMAIN (engine-switch decision (ii) restatement, plan 02 §8): this H/3
+ * line is the CLOCK-side e — the eclipse Sun's equation of centre (FQ-7-Sun
+ * option C-small, doc 66 §5), the besselian Sun distance, the cardinal
+ * braid and the scene machinery ride it. The LUNAR chain (E-factor,
+ * modulation, channelIntegral, argument Δe²) rides the ONE deep e
+ * (deep-ecc-channel.cjs — the engine's own ±10-Myr z-vector); the two
+ * agree within 4.2e-5 across the historical era, and this line keeps its
+ * E18 role as the epoch-local tangent. Eccentricity remains
+ * frame-invariant (fixed-frame periods only); the H/16 perihelion cycle
+ * (the H/3 rotation seen from the H/13 equinox, 13 + 3 = 16) belongs to
+ * ϖ, not e (doc 108).
  *
  * Extracted VERBATIM from src/script.js (_FW_ECC, _fwEarthEcc,
  * _eCompModulation, _fwChannelIntegral, the framework branch of _fwEFactor),