@essrt/physics 25.0.1 → 26.0.1

package/package.json

@@ -1,9 +1,9 @@
 {
   "name": "@essrt/physics",
-  "version": "25.0.1",
+  "version": "26.0.1",
   "description": "ESSRT (Expanding Solar System Resonance Theory) physics core — 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": "v27.0",
+    "modelVersion": "v27.1",
     "_note": "The model identity this package version ships (IP-unified-architecture §10: package semver tracks the API axis; the model version + CONSTANTS_HASH/COEFFICIENTS_HASH in src/constants/ carry the scientific identity). The version-pinning gate enforces this pairing against public/input/model-version.json."
   },
   "type": "module",
@@ -39,6 +39,7 @@
     "./math/angles": { "types": "./types/math/angles.d.cts", "default": "./src/math/angles.cjs" },
     "./cardinal/one-source-structure": { "types": "./types/cardinal/one-source-structure.d.cts", "default": "./src/cardinal/one-source-structure.cjs" },
     "./moon/deep-ecc-channel": { "types": "./types/moon/deep-ecc-channel.d.cts", "default": "./src/moon/deep-ecc-channel.cjs" },
+    "./moon/ecc-sensitivities": { "types": "./types/moon/ecc-sensitivities.d.cts", "default": "./src/moon/ecc-sensitivities.cjs" },
     "./moon/deep-modes-artifact": { "types": "./types/moon/deep-modes-artifact.d.cts", "default": "./src/moon/deep-modes-artifact.cjs" },
     "./earth/deep-orbital-history": { "types": "./types/earth/deep-orbital-history.d.cts", "default": "./src/earth/deep-orbital-history.cjs" },
     "./earth/sidereal-year-channel": { "types": "./types/earth/sidereal-year-channel.d.cts", "default": "./src/earth/sidereal-year-channel.cjs" },

package/src/chain-cycles/index.cjs

@@ -28,7 +28,9 @@
  * NUMERICS ARE LOAD-BEARING: fixed 10-yr grid over ±250 kyr, per-cell
  * 3-point Simpson at build, LINEAR interpolation on read; beyond the dense
  * grid a COARSE cumulative table (1-kyr cells, the same per-cell Simpson,
- * anchored at the dense edge) grown lazily in 250-kyr chunks in the direction
+ * anchored at the dense edge; read with a 3-point Simpson over the partial
+ * cell, so the count's RATE is continuous — plan 23 follow-up, see the tier
+ * note below) grown lazily in 250-kyr chunks in the direction
  * asked for — the owner's playback at 1,000 yr/s past year 250,000 had fallen
  * to 190 ms of JS per frame: the former fallback, an adaptive Simpson of up to
  * 1,024 samples per call cached by exact year pair, re-integrated every one
@@ -136,9 +138,19 @@ function createChainCycleIntegrator({
   // Per chain and per direction, a cumulative table in COARSE_STEP cells that
   // starts at the dense edge with the edge's own value and grows by COARSE_CHUNK
   // years whenever a year beyond its reach is asked for. The cells are built
-  // with the same 3-point Simpson as the dense grid; the read is the same
-  // linear interpolation. The period laws vary on Myr scales at deep time, so
-  // 1-kyr cells hold the count to far below a cycle — the old fallback sampled
+  // with the same 3-point Simpson as the dense grid; the read is a 3-point
+  // Simpson over the partial cell from the cell's start (three integrand
+  // evaluations per read), NOT the dense grid's linear interpolation: a
+  // linear read makes the count's rate piecewise constant per cell, and the
+  // of-date node and perigee are SMALL differences of two large counts
+  // (N_drac − N_trop: 13.4 − 13.37 cycles/yr), so each table's per-cell slope
+  // step — 1e-6 of its own rate, the month's tidal drift over 1 kyr — landed
+  // fully on the difference: a ±0.3-d sawtooth per cell on the node period
+  // beyond ±250 kyr (measured on the Node Regression panel, plan 23
+  // follow-up). On the dense 10-yr grid the same step is 100× smaller and
+  // below every reader, so its linear read — the byte-identical in-window
+  // path — stays. The period laws vary on Myr scales at deep time, so 1-kyr
+  // cells hold the count to far below a cycle — the old fallback sampled
   // ~1 per kyr too, but re-did the whole span on every call.
   const COARSE_STEP = 1000, COARSE_CHUNK = 250000, COARSE_LIMIT_YEARS = 600e6;
   /** @typedef {{ cum: number[], ok: boolean }} CoarseSide */
@@ -178,9 +190,13 @@ function createChainCycleIntegrator({
     if (distance > COARSE_LIMIT_YEARS) return null;
     if (side.cum.length < Math.ceil(distance / COARSE_STEP) + 1) growCoarse(periodFnSeconds, side, edgeYear, dir, distance);
     if (!side.ok || side.cum.length < Math.ceil(distance / COARSE_STEP) + 1) return null;
-    const idxF = distance / COARSE_STEP;
-    const i = Math.floor(idxF);
-    return side.cum[i] + (idxF - i) * (side.cum[i + 1] - side.cum[i]);
+    const i = Math.floor(distance / COARSE_STEP);
+    // the partial cell from its start y0 to y: a 3-point Simpson, signed by (y − y0)
+    const y0 = edgeYear + dir * i * COARSE_STEP;
+    if (y === y0) return side.cum[i];
+    const f0 = integrandAt(periodFnSeconds, y0), fm = integrandAt(periodFnSeconds, (y0 + y) / 2), f1 = integrandAt(periodFnSeconds, y);
+    if (f0 === null || fm === null || f1 === null) return null;
+    return side.cum[i] + (f0 + 4 * fm + f1) * ((y - y0) / 6);
   }
   /** C(y) anywhere the tables reach: dense inside the grid, coarse beyond.
    *  @param {PeriodSecondsAtAge} periodFnSeconds @param {Float64Array} tab @param {number} y @returns {number | null} */

package/src/constants/generated.js

@@ -39,7 +39,7 @@
 export const CONSTANTS_HASH = "da713fc2d846767d";
 
 /** Model version label — single source: public/input/model-version.json (§10 two-axis scheme). */
-export const MODEL_VERSION = "v27.0";
+export const MODEL_VERSION = "v27.1";
 
 /** Canonical preprint DOI — single source: public/input/model-version.json. */
 export const PREPRINT_DOI = "10.21203/rs.3.rs-8758810/v4";

package/src/deltat/deep-time-chain.js

@@ -43,14 +43,9 @@ import { createMoonMonthChain } from '../moon/month-chain.cjs';
 import { massLossScaledLinearAtAge, driver2PeriodSecondsAtAge, synodicPeriodSeconds } from '../planets/orbit-chain.cjs';
 import { computeMassChainJ2000, computeMoonPeriodsJ2000 } from '../earth/mass-chain.cjs';
 
-/**
- * Moon-channel eccentricity sensitivities (perigee/node), the [g/g₀]^s
- * exponents of the factored deep-time law — Meeus-effective (S_N moved
- * 1.0 → 1.018 with the v4 frame-attribution batch). ONE home since B2 (they
- * were literals in model.js, tools/lib/deep-time.js and script.js `_FW_MOON`).
- */
-export const MOON_ECC_SENSITIVITY_PERIGEE = 2.407;
-export const MOON_ECC_SENSITIVITY_NODE = 1.018;
+// The lunar perigee/node eccentricity sensitivities — ONE home moon/ecc-sensitivities.cjs
+// (plan 23): the Meeus-effective pair and the PHYSICAL pair derived from it there.
+import { deriveMoonEccSensitivities } from '../moon/ecc-sensitivities.cjs';
 
 /**
  * The ΔT stack's switches. `resonator` gates the swing's δLOD terms; the
@@ -127,6 +122,18 @@ export function createDeepTimeChain(deps) {
   /** p₀, ″/yr. @returns {number} */
   const precessionRateJ2000ArcsecPerYr = () => 1296000 / certifiedAxialPrecessionJ2000Years();
 
+  // The lunar perigee/node eccentricity sensitivities (plan 23): the Meeus-effective pair and the
+  // PHYSICAL pair derived from it by the T² identity (moon/ecc-sensitivities.cjs, the ONE home).
+  // The month chain and the arguments module run the physical pair; the frame share the effective
+  // pair carried rides the model's own precession of date instead (the arguments' pFix, on L′, ϖ
+  // and Ω alike), so the frame-free arguments M′ and F cancel it at every epoch.
+  const moonEccSensitivities = deriveMoonEccSensitivities({
+    eccentricityJ2000: C.earthOrbital.earthEccentricityJ2000,
+    eccentricityDotJ2000: C.earthOrbital.earthEccentricityDotJ2000,
+    frameT2ArcsecPerCy2: C.moonMeeus.elpW1T2Decomposition_arcsecPerCy2.generalPrecessionPA_T2_Lieske1976,
+    generalPrecessionArcsecPerYr: C.earthOrbital.generalPrecessionJ2000_arcsecPerYr,
+  });
+
   // ── α(t): the GIA channel (lagged L1 response), the memo and the lattice pin ──
   const CLIMATE_REGIME = F.CLIMATE_FORMULA_COEFFS.regimes['lr04-post-mpt'];
   const alphaGia = createAlphaGiaChannel({
@@ -308,7 +315,8 @@ export function createDeepTimeChain(deps) {
       moonApsidalJ2000Seconds,
       moonNodalJ2000Seconds,
       moonSiderealMonthJ2000Seconds,
-      sPerigee: MOON_ECC_SENSITIVITY_PERIGEE, sNode: MOON_ECC_SENSITIVITY_NODE,
+      // plan 23: the PHYSICAL pair in the channel; the frame share rides the arguments' pFix
+      sPerigee: moonEccSensitivities.perigeePhysical, sNode: moonEccSensitivities.nodePhysical,
     },
     fns: {
       meanLodSecondsAtAge,
@@ -405,7 +413,7 @@ export function createDeepTimeChain(deps) {
     EPOCH_PARAMS, ageAnchorYear,
     // the pieces
     alphaGia, earthMoiFactorAtAge, withLatticeAlpha, j2RatioAtYear, iEarthAtAge,
-    recession, meanObliquity, dtCycles, s3ArcsecPerYr, deepLod, l0, deepEcc, eCompModulation, moonChain,
+    recession, meanObliquity, dtCycles, s3ArcsecPerYr, deepLod, l0, deepEcc, eCompModulation, moonChain, moonEccSensitivities,
     // the family
     meanMoonDistanceMetresAtAge: recession.distanceMetresAtAge,
     meanMoonDistanceAtAge: moonChain.distanceKmAtAge,

package/src/eclipse/besselian.cjs

@@ -35,11 +35,10 @@
  *     the derived series-extension tail — deliberate for the certified
  *     finder statistics (knife-edge canon events), required here. The
  *     injected moonFullAtDaysTT is the series' sceneEvalAt — the same
- *     evaluation the browser scene ships (the extension inside it since
- *     plan 06 R3 item 1: one scene Moon).
+ *     evaluation the browser scene ships (the extension inside it: one scene
+ *     Moon).
  *
- * APPARENT PLACES (plan 06 R3 item 1 — measured, the cancellation the I2
- * Sun exposed): both bodies are evaluated mean-of-date (nutation is a common
+ * APPARENT PLACES (measured — the cancellation the I2 Sun exposed): both bodies are evaluated mean-of-date (nutation is a common
  * rotation of Sun and Moon and cancels in the elongation) but with their
  * PLANETARY ABERRATION: the Sun's derived κ = 2π a/(c T_sid √(1−e²))
  * (≈20.5″, the light-time of a source at rest in the barycentric frame) and
@@ -131,8 +130,8 @@ function createBesselian(deps) {
   }
 
   /** Elliptical Sun distance — FRAMEWORK-NATIVE: the model's own
-   *  eccentricity (Law-5-locked base + drift) and perihelion longitude
-   *  (H/16 precession) close the ellipse, with the true anomaly from the
+   *  eccentricity and perihelion longitude of date close the ellipse, with
+   *  the true anomaly from the
    *  already-computed sun longitude (v = λ − (ϖ + 180°), geocentric-perigee
    *  convention). Agrees with the Meeus Ch. 25 polynomial to ±0.001% over a
    *  year, with no fitted external coefficients, and stays honest at deep

package/src/eclipse/sun-planetary-completion.cjs

@@ -68,8 +68,7 @@
  *       inconsistency for elongation. Removed: the nutation-free syzygy
  *       fleet improves 3.877 → 3.756″ and the fuller-bridged Sun
  *       comparison 2.344 → 2.123″ (measured by the fq7s preview and
- *       annual-channel instruments, deleted at plan 09 step 7e — docs/109
- *       §22). This file now carries ZERO fitted
+ *       annual-channel instruments — docs/109 §22). This file carries ZERO fitted
  *       constants. The Earth-around-EMB wobble is DERIVED: amplitude
  *       a_M·μ/AU with μ = 1/(1+M_E/M_M), injected by the model wiring
  *       from live package constants (6.4399″ at current constants), sign
@@ -143,8 +142,8 @@ const ARG_D0 = 297.8501921;
  * D2 window (the term's in-window ramp went to the projected-out secular
  * basis; the 240-yr term folded into it). WHY they were missing: the
  * certified Sun integrates a smooth tropical year for its mean longitude;
- * against Horizons its residual carried exactly this 6″ ripple (plan 06 I1
- * analysis) — with the same phase in the sidereal residual, a perturbation
+ * against Horizons its residual carried exactly this 6″ ripple — with the
+ * same phase in the sidereal residual, a perturbation
  * of Earth's mean motion, not a frame effect. Sub-0.2″ candidates left out
  * (2J−5S 0.17″, 5V−8E 0.17″: window-dependent phase). The CARRIERS moved
  * with I3 from the planet records' rounded of-date periods minus p₀ to the
@@ -225,7 +224,7 @@ const TERMS = [
   [[0, 0, 1, -1, 0, 0], [0, 0, 1, 0, 0, 0], 0.0188, 0.0523],
   [[0, 2, -2, 0, 0, 0], [0, 0, 1, 0, 0, 0], -0.0421, -0.0340],
   [[0, 0, 1, -1, 0, 0], [0, 0, -1, 0, 0, 0], -0.0353, -0.0379],
-  // plan 06 I3 — the Venus–Earth term 8λ_V − 13λ_E (≈239 yr on the banked
+  // The Venus–Earth term 8λ_V − 13λ_E (≈239 yr on the banked
   // carriers): a CONSTANT row, derived on the model's own run over ±3100 yr
   // (i3-long-period-dalembert.mjs `2 20000 3100 4`; the I2 row read 1.84/0.14″
   // on the record carriers over −5100..+1100). The term is formally FIFTH order
@@ -394,13 +393,9 @@ function createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy, e
   return { sunPlanetaryCompletionDeg };
 }
 
-/* (Plan 09 step 5a: PAIRED_SUN_HARMONICS_SHA256 — the sha256/16 of the Step-0
- *  SUN_LONGITUDE_HARMONICS this table was derived beside, last 53530a28a60b3c4a
- *  (the R10 rotation of a comb never refit since the eccentricity unification)
- *  — is retired with the comb. The comb left every Sun path at plan 06 R4, so
- *  the residual the table was derived from is the certified Sun's; the table
- *  stands (fidelity 0.616″, 70 terms) and the carrier-rate pair below remains
- *  its enforced matched pair.) */
+/* The table was derived from the certified Sun's residual and stands (fidelity
+ *  0.616″, 70 terms); the carrier-rate pair below is its enforced matched pair.
+ *  (The retired paired-harmonics fingerprint: doc 109 §23.) */
 
 /** sha256/16 of JSON.stringify([...planets, moonElongation]) — the seven
  *  full-precision carrier rates (deg/cy TT) the tables pair with: since plan

package/src/index.js

@@ -108,7 +108,10 @@ export { DEG2RAD, RAD2DEG, ARCSEC_PER_TURN, ARCSEC_PER_RAD, wrap360, wrap180, wr
 // cited Moon-distance history, the ΔT cycle stack, |s₃|, the LOD/ΔT core, the
 // lunar month chain and the mean*AtAge family, assembled ONCE from the
 // generated constants; the Node twin and the browser destructure it.
-export { createDeepTimeChain, MOON_ECC_SENSITIVITY_PERIGEE, MOON_ECC_SENSITIVITY_NODE } from './deltat/deep-time-chain.js';
+export { createDeepTimeChain } from './deltat/deep-time-chain.js';
+// The lunar perigee/node eccentricity sensitivities — the Meeus-effective pair, the physical pair's
+// derivation and the Meeus Ch. 47 argument rates, ONE home (plan 23).
+export { MOON_ECC_SENSITIVITY_PERIGEE, MOON_ECC_SENSITIVITY_NODE, deriveMoonEccSensitivities, MEEUS_ARGUMENT_RATES_DEG_PER_CY } from './moon/ecc-sensitivities.cjs';
 // The J2000 mass chain, the Moon's J2000 periods and the unit constants — ONE home.
 export { computeMassChainJ2000, computeMoonPeriodsJ2000 } from './earth/mass-chain.cjs';
 export { SECONDS_PER_DAY, DAYS_PER_JULIAN_YEAR, JULIAN_CENTURY_DAYS, SECONDS_PER_JULIAN_YEAR } from './constants/units.cjs';

package/src/model/cardinal-instants.js

@@ -57,7 +57,7 @@ export function assembleCardinalInstants(ctx) {
   const CARDINAL_SI_YEAR_D = 365.2422;
   /** The crossing NEAREST a seed instant — the unique root within ±½ tropical
    *  year of the seed (the wrapped longitude difference is the Newton residual).
-   *  Plan 06 R4c: the seed is a JD because the caller's year may live on a
+   *  The seed is a JD because the caller's year may live on a
    *  different axis than the model year — the simulator's cardinal panel
    *  labels its year with the displayed CALENDAR date (the Julian calendar
    *  before 1582, 365.25-d years) while `year` below is the model year (the SI
@@ -79,13 +79,12 @@ export function assembleCardinalInstants(ctx) {
   /** The TRUE-UT JD at which the APPARENT Sun's longitude of date equals
    *  the cardinal target in the given year — Newton on the crossing (the
    *  Sun's rate 360°/tropical year; a handful of steps from the tropical-year
-   *  seed). Replaces the fitted cardinal model's instants (R1). Plan 06 R3
-   *  item 2: the Newton root lives on the FINDER axis (the curve-only ΔT the
+   *  seed). The Newton root lives on the FINDER axis (the curve-only ΔT the
    *  eclipse finders certify), and the returned instant is that root minus the
    *  deltaTStart bridge — true UT = TT − (deltaTStart + curve), the same
    *  bridge the besselian and the registry instruments apply. Before this the
    *  published instants were 54.55 s (0.9 min) late. Against Horizons over
-   *  ±3000 yr the remaining offset is the Sun's (plan 06 I2 measured).
+   *  ±3000 yr the remaining offset is the Sun's (measured).
    *  `year` is the MODEL year (the SI axis) — see cardinalCrossingNearJdUT.
    *  @param {number} year @param {string} type @returns {number} */
   const cardinalCrossingJdUT = (year, type) => {

package/src/model/deep-time-chain.js

@@ -12,8 +12,8 @@ import { DEEP_MODES_ARTIFACT } from '../moon/deep-modes-artifact.cjs';
 export function assembleDeepTimeChain(ctx) {
   const { C, F, late, startmodelYear } = ctx;
   // 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
+  // pin and the exact-argument memo), the Layer-0 bundle, the cited
+  // Earth–Moon distance history, 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);
@@ -39,7 +39,7 @@ export function assembleDeepTimeChain(ctx) {
    *  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)
+   *  date: the hybrid's own ε(t) carries the ONE cosine, so the cited mean obliquity's factor
    *  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

package/src/model/lunar-theory.js

@@ -23,7 +23,7 @@ export function assembleLunarTheory(ctx) {
   // tools/verify/eclipse-audit.js (finders). Deep-time and framework-native
   // are hardwired ON here: they are the shipped defaults; the A/B env
   // toggles (SG_DEEP_TIME/MOON_ARGS_PURE_MEEUS) stay an engine concern.
-  const j2000JD = J2000_JD;                       // ONE home: constants/units (plan 14 step 2)
+  const j2000JD = J2000_JD;                       // ONE home: constants/units
   const julianCenturyDays = JULIAN_CENTURY_DAYS;
 
   // J2000 Moon precession anchors (the IAU inputs themselves). The
@@ -50,7 +50,7 @@ export function assembleLunarTheory(ctx) {
   // Float64Array tables key correctly and persist.
   const chainCycles = createChainCycleIntegrator({
     ageAnchorYear: startmodelYear,
-    tropicalYearJ2000Seconds: meanTropicalYearJ2000Seconds,   // the SI tick (plan 09 step 7c′)
+    tropicalYearJ2000Seconds: meanTropicalYearJ2000Seconds,   // the SI tick
     isDeepTime: () => true,
   });
   /** @param {number} tMa @returns {number|null} */
@@ -72,7 +72,7 @@ export function assembleLunarTheory(ctx) {
   const mcAnomalistic = (a, b) => chainCycles.cyclesBetween(anomalisticMonthPeriodFn, a, b);
   /** @param {number} a @param {number} b @returns {number|null} */
   const mcJupiter = (a, b) => chainCycles.cyclesBetween(jupiterPeriodFn, a, b);
-  /** The Moon's sidereal (star-referenced) orbits — the series' A3 argument (plan 09 step 5c). @param {number} a @param {number} b @returns {number|null} */
+  /** The Moon's sidereal (star-referenced) orbits — the series' A3 argument. @param {number} a @param {number} b @returns {number|null} */
   const mcSidereal = (a, b) => chainCycles.cyclesBetween(siderealMonthPeriodFn, a, b);
   /** @param {number} a @param {number} b @returns {number|null} */
   const mcApsidalOfDate = (a, b) => {
@@ -98,6 +98,8 @@ export function assembleLunarTheory(ctx) {
       elpGeneralPrecessionPA_T2ArcsecPerCy2: C.moonMeeus.elpW1T2Decomposition_arcsecPerCy2.generalPrecessionPA_T2_Lieske1976,
       // Decision (ii): the lunar chain reads ONE deep e end to end
       eccE0: deepEcc.e0,
+      // plan 23: the perigee/node sensitivities — the chain's derived pair (ONE home)
+      moonEccSensitivities: chain.moonEccSensitivities,
     },
     fns: {
       eccAt: /** @param {number} tYr */ (tYr) => deepEcc.eccAt(tYr),
@@ -191,8 +193,7 @@ export function assembleLunarTheory(ctx) {
       // against the MODEL-START anchor (startmodelYear, the browser's
       // J2000_CALENDAR_YEAR — NOT this model's 2000.0 yearToTMa convention:
       // the 0.5-yr difference is 4.95e-11 of distance and failed the
-      // bit-exact parity gate; plan 11 step 5 dropped the −0.83-d correction
-      // from the anchor, both twins alike). Falls back to the J2000 constant
+      // bit-exact parity gate). Falls back to the J2000 constant
       // without a jd (legacy call shape) or past the domain.
       getMoonDistanceKm: /** @param {number} [jdTT] */ (jdTT) => {
         if (jdTT === undefined) return moonDistanceKm;

package/src/model/one-source-movement.js

@@ -32,8 +32,8 @@ export function assembleOneSourceMovement(ctx) {
   // path; the npm package alone has no 8 MB series), the movement runs the
   // SERIES tier exactly like the browser — killing the mode-vs-series
   // value split (measured: 6 s on the anomalistic year at J2000).
-  // The ONE-SOURCE movement: the hybrid's ε(t) (plan 06 Phase 3 S3b — THE
-  // published obliquity) and the one-family year lengths, on the same
+  // The ONE-SOURCE movement: the hybrid's ε(t) (THE published obliquity) and
+  // the one-family year lengths, on the same
   // per-tier sampler.
   const oneSourceM = (() => {
     const AEarth = /** @type {any} */ (CHAIN_ARTIFACT).j2000AnchorElements.earth;
@@ -47,7 +47,7 @@ export function assembleOneSourceMovement(ctx) {
     // 108.6 ± 8.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
+    // 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
@@ -73,13 +73,13 @@ export function assembleOneSourceMovement(ctx) {
       obliquityJ2000Deg: C.earthOrbital.obliquityJ2000_deg,
       // 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.
+      // factor must not enter twice.
       axialPrecessionYearsAtYearFn: (yr) => deepLod.lunisolarTorquePeriodYearsAtAge((startmodelYear - yr) / 1e6) ?? axial0,
       // THE PRECESSION OF DATE: 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, so the J2000 anchor is
       // untouched. The channel is
-      // read PURE (alphaGia, not the lattice-pin-aware earthMoiFactorAtAge):
+      // read PURE (alphaGia, not the pinned earthMoiFactorAtAge):
       // the samplers are cached, their values must be pure in year.
       solarTorqueShareJ2000: EPOCH_PARAMS.precessionSolarShareJ2000,
       // options.giaEllipticityOfDate === false: the orbital forcing alone (the T1 line derivation).
@@ -100,7 +100,7 @@ export function assembleOneSourceMovement(ctx) {
       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
+        // walk when a deeper epoch arrives. 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;
@@ -127,7 +127,7 @@ export function assembleOneSourceMovement(ctx) {
       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. */
+      /** The hybrid's own domain, [yearFrom, yearTo] — the banked ζ series' span; 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),
@@ -139,7 +139,7 @@ export function assembleOneSourceMovement(ctx) {
   })();
   const yearLengthsM = oneSourceM.yearLengths;
   /** THE one-family tropical year of date in SI seconds: route B inside the
-   *  one-family window, the tidal-chain year beyond (plan 06 R1 — the Sun's
+   *  one-family window, the tidal-chain year beyond (the Sun's
    *  mean longitude, the cardinal instants and the published year lengths
    *  all ride this one function). @param {number} year @returns {number} */
   const oneFamilyTropicalYearSeconds = (year) => (Math.abs(year - 2000) <= ONE_FAMILY_WINDOW_YEARS

package/src/model/planet-records.js

@@ -45,10 +45,10 @@ export function assemblePlanetRecords(ctx) {
     const gArc = chainModes.g[k] && chainModes.g[k].arcsecPerYr;
     const sArc = chainModes.s[k] && chainModes.s[k].arcsecPerYr;
     PUBLISHED_PLANET_RECORDS[k] = Object.freeze({
-      name: PLANET_DISPLAY_NAMES[k],   // presentation — planets/display-names, the deleted JSON's last leaf (plan 11 step 6)
+      name: PLANET_DISPLAY_NAMES[k],   // presentation — planets/display-names
       orbitalPeriodInputDays: ar.orbitalPeriodInputDays,
       longitudePerihelion: ar.longitudePerihelion,
-      ascendingNodeInvPlane: ar.ascendingNodeInvPlaneClosedForm,   // the closed-form verified J2000 node — astro-reference since plan 11 step 2
+      ascendingNodeInvPlane: ar.ascendingNodeInvPlaneClosedForm,   // the closed-form verified J2000 node (astro-reference)
       invPlaneInclinationJ2000: ar.invPlaneInclinationJ2000,
       orbitalEccentricityJ2000: ar.orbitalEccentricityJ2000,
       axialTiltJ2000: ar.axialTiltJ2000,

package/src/model/published-periods.js

@@ -8,7 +8,7 @@ import { ONE_FAMILY_WINDOW_YEARS } from '../earth/year-lengths.cjs';
 /** @param {import('./sun-and-eclipses.js').Ctx} ctx @returns the bindings the later slices read */
 export function assemblePublishedPeriods(ctx) {
   const { certifiedAxialPrecessionJ2000Years, deepLod, yearLengthsM, yearToTMa } = ctx;
-  // Plan 06 Phase 3 S2 — THE published of-date precession period and mean
+  // THE published of-date precession period and mean
   // tropical year: the one-family route (B) inside the banked tiers
   // (|Δyear| ≤ ONE_FAMILY_WINDOW_YEARS), the composed lunisolar period
   // beyond (S5: on the derived J2000 anchor, so the two faces agree at

package/src/model/sun-and-eclipses.js

@@ -56,7 +56,7 @@ export function assembleSunAndEclipses(ctx) {
     const cyclesTo = (year) => {
       const dy = year - 2000;
       const s = dy >= 0 ? 1 : -1, a = Math.abs(dy);
-      // WARM THE SAMPLER TO THE TARGET FIRST (plan 06 R4, measured): the
+      // WARM THE SAMPLER TO THE TARGET FIRST (measured): the
       // table below grows outward in 1-yr / 100-yr steps, and every step
       // beyond ~1.6 Myr enlarged the one-family sampler's >2-Myr tier just
       // past its span, rebuilding the whole deep integration about every 32
@@ -131,7 +131,7 @@ export function assembleSunAndEclipses(ctx) {
     * (648000 / Math.PI);
   // FQ-5 N3 — the completion's carriers are FRAMEWORK-derived from the
   // model's own planet period records (the Moon-elongation rate from the
-  // sidereal month/year identity). Plan 06 I2 — the carriers are SIDEREAL
+  // sidereal month/year identity). The carriers are SIDEREAL
   // (measured): the records' orbitalPeriodInputDays are OF-DATE periods (their
   // rates sit +1.39..1.62°/cy above the VSOP sidereal mean motions — the
   // framework precession plus input rounding), and a perturbation argument
@@ -169,7 +169,7 @@ export function assembleSunAndEclipses(ctx) {
     },
     sunLonDegAt: /** @param {number} jdUT */ (jdUT) => eclipseFinders.sunLonDegAt(jdUT),
     sunCompletionDeg: sunPlanetaryCompletionDeg,
-    // plan 06 R3 item 1: the location tier on APPARENT places — the Sun's derived κ
+    // The location tier on APPARENT places — the Sun's derived κ
     // (the same helper the cardinal instants use), the Moon's relative light-time
     // inside the besselian; the series-extension tail now lives in sceneEvalAt.
     sunAberrationDegAt: /** @param {number} year */ (year) => late.sunAberrationArcsecAt(year) / 3600,

package/src/model/surface.js

@@ -94,7 +94,7 @@ export function assembleSurface(ctx) {
       // one-family inside the window, the tidal-chain mean beyond.
       tropicalYearSecondsAtYear: publishedTropicalYearSecondsAtYear,
     }),
-    // Plan 06 Phase 3 S3 → S5 — the lunisolar precession clock as a
+    // The lunisolar precession clock as a
     // first-class surface: Earth's spin clock, ONE home for its published
     // faces, spoken in PERIODS and their RATIOS — no unit, no integer.
     // ONE J2000 reading: every face reads 25,771.6 yr at J2000 — the
@@ -108,7 +108,7 @@ export function assembleSurface(ctx) {
     lunisolar: Object.freeze({
       /** 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 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} */
+      /** 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 (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; },
@@ -212,9 +212,9 @@ export function assembleSurface(ctx) {
       // 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()). */
+      /** The hybrid's own domain, [yearFrom, yearTo] — the banked ζ series' span, ±10 Myr about J2000 (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} */
+      /** The cited mean obliquity of date ε̄(t) (Farhat et al. 2022's column), 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,
@@ -234,7 +234,7 @@ export function assembleSurface(ctx) {
       siderealDaySeconds,
       stellarDaySeconds,
     }),
-    // Plan 06 R1: the instants are the crossings of the ONE Sun (the
+    // The instants are the crossings of the ONE Sun (the
     // certified completed Sun the scene renders), UT model-JD; the RA at a
     // crossing is the target longitude by construction (λ = 0/90/180/270 ⇒
     // RA = 0/90/180/270, any obliquity); the per-type year length is the
@@ -266,6 +266,10 @@ export function assembleSurface(ctx) {
       distanceKmAtYear: /** @param {number} year @returns {number} */ (year) => moonDistanceKmAtYearOrNaN(year),
       siderealMonthDaysAtYear: moonSiderealMonthDaysAt,
       synodicMonthDays: moonSynodicMonthDays,
+      // The perigee/node eccentricity sensitivities (plan 23): the Meeus-effective
+      // pair the record fixes and the PHYSICAL pair the runtime derives from it
+      // (moon/ecc-sensitivities.cjs, the ONE home) — the registry's source.
+      eccSensitivities: chain.moonEccSensitivities,
       // The apparent-position chain (truncated Meeus Ch. 47 series on
       // framework-native arguments). JD(UT) axis — UT→TT applied internally.
       lonDegAtJD: /** @param {number} jd @returns {number} */ (jd) => moonSeries.truncatedLonDeg(jd),
@@ -288,7 +292,7 @@ export function assembleSurface(ctx) {
       sunLonCompletedDegAtJD: /** @param {number} jd @returns {number} */ (jd) =>
         eclipseFinders.sunLonDegAt(jd)
         - sunPlanetaryCompletionDeg((jdTTFromUT(jd) - j2000JD) / julianCenturyDays),
-      /** The APPARENT Sun (completed geometric − the derived aberration constant + the leading nutation terms on the model's own arguments) — the quantity Horizons' observer ecliptic longitude is; the cardinal instants are its crossings (plan 06 R1/I1). @param {number} jdUT @returns {number} */
+      /** The APPARENT Sun (completed geometric − the derived aberration constant + the leading nutation terms on the model's own arguments) — the quantity Horizons' observer ecliptic longitude is; the cardinal instants are its crossings. @param {number} jdUT @returns {number} */
       sunApparentLonDegAtJD: /** @param {number} jdUT @returns {number} */ (jdUT) => sunApparentLonDegAtJdUT(jdUT),
       findLunarInRange: /** @param {number} jdStart @param {number} jdEnd */ (jdStart, jdEnd) => eclipseFinders.findLunarEclipsesInRange(jdStart, jdEnd),
       findSolarInRange: /** @param {number} jdStart @param {number} jdEnd */ (jdStart, jdEnd) => eclipseFinders.findSolarEclipsesInRange(jdStart, jdEnd),