@essrt/physics 25.0.1 → 26.0.0

package/package.json

@@ -1,9 +1,9 @@
 {
   "name": "@essrt/physics",
-  "version": "25.0.1",
+  "version": "26.0.0",
   "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/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),

package/src/moon/arguments.cjs

@@ -20,8 +20,12 @@
  *    with p the IAU general precession.
  *
  * Construction (docs/66 §1, abridged): linear rates are observational J2000
- * anchors; perigee/node ride the PHASE-AWARE e_E channel integral (s_ϖ
- * 2.407 / s_Ω 1.018, Meeus-effective); Lp's planetary T² remainder rides
+ * anchors; perigee/node ride the PHASE-AWARE e_E channel integral on the
+ * PHYSICAL sensitivities (s_ϖ 2.479 / s_Ω 0.867 — derived from the
+ * Meeus-effective 2.407 / 1.018 by removing the frame's T²;
+ * moon/ecc-sensitivities.cjs, the ONE home) plus the model's own precession
+ * of date, pFix, the ONE frame term every of-date longitude carries (plan 23
+ * — so M′ and F are frame-free at every epoch); Lp's planetary T² remainder rides
  * the bounded e_E² carrier + the obliquity-line carrier (K_PL derived
  * lazily from the closed v4 budget — zero new constants); D and M are
  * identity-composed (D ≡ Lp − L_sun, M ≡ L_sun − ϖ_sun) on the injected
@@ -37,6 +41,8 @@
 // reform). ONE home — time/calendar.cjs; re-exported here for the readers
 // that import it from the arguments module.
 const { jdToDecimalYear } = require('../time/calendar.cjs');
+// The Meeus Ch. 47 argument rates — ONE home (plan 23; the sensitivities' derivation reads the same set).
+const { MEEUS_ARGUMENT_RATES_DEG_PER_CY } = require('./ecc-sensitivities.cjs');
 
 /**
  * @typedef {{Lp: number, D: number, M: number, Mp: number, F: number}} MoonArgsDeg
@@ -59,6 +65,7 @@ const SECULAR_COMPLETION_TAPER_YEARS = 27943.083333333332;
  *     elpEarthFigureJ2ArcsecPerCy2: number,
  *     elpGeneralPrecessionPA_T2ArcsecPerCy2: number,
  *     eccE0: number,
+ *     moonEccSensitivities: ReturnType<typeof import('./ecc-sensitivities.cjs').deriveMoonEccSensitivities>,
  *   },
  *   fns: {
  *     eccAt: (tYr: number) => number,
@@ -94,8 +101,9 @@ function createMoonArguments({ constants, fns }) {
     j2000JD, julianCenturyDays,
     eccentricityJ2000, eccentricityDotJ2000, eccentricityDotDotJ2000,
     elpEarthFigureJ2ArcsecPerCy2, elpGeneralPrecessionPA_T2ArcsecPerCy2,
-    eccE0,
+    eccE0, moonEccSensitivities,
   } = constants;
+  if (!moonEccSensitivities) throw new Error('moon/arguments: constants.moonEccSensitivities is required (the deep-time chain\'s derived pair — plan 23)');
   const {
     eccAt, channelIntegral, computeObliquityEarth, jdToSIyear,
     sunMeanLonDegAtJdTT, sunPerihelionLonDegAtJdTT,
@@ -108,17 +116,26 @@ function createMoonArguments({ constants, fns }) {
   const bundle = (() => {
     const LP0 = 218.3164477, D0 = 297.8501921, M0 = 357.5291092,
           MP0 = 134.9633964, F0 = 93.2720950;
-    const LPR = 481267.88123421, DR = 445267.1114034, MR = 35999.0502909,
-          MPR = 477198.8675055, FR = 483202.0175233;
+    const { Lp: LPR, D: DR, M: MR, Mp: MPR, F: FR } = MEEUS_ARGUMENT_RATES_DEG_PER_CY;   // ONE home (plan 23)
     const WDOT = LPR - MPR;   // perigee ϖ̇ of-date (+4069.0137) = ϖ̇_ICRF + p
     const NDOT = LPR - FR;    // node   Ω̇ of-date (−1934.1363) = Ω̇_ICRF + p
     // e_E channel Taylor CHECKS (documented J2000 reference values — the live
     // path is the phase-aware channelIntegral; T2_*/T3_* are NOT consumed):
     const E0 = eccentricityJ2000, EDOT0 = eccentricityDotJ2000;
     const KAPPA = 3 * E0 * EDOT0 / (1 - E0 * E0);
-    // v4 frame attribution: FRAME-EFFECTIVE exponents (physical pair after
-    // removing the IAU ṗ_A frame term: s_ϖ 2.479 / s_Ω 0.867).
-    const S_W = 2.407, S_N = 1.018;
+    // The sensitivities (plan 23; moon/ecc-sensitivities.cjs, the ONE home). The
+    // Meeus-EFFECTIVE pair S_W/S_N reproduces Meeus's of-date T² exactly — the
+    // documented check below. The live channel runs the PHYSICAL pair
+    // S_W_PHYS/S_N_PHYS (the frame's T² removed; the three-body lab's pure-gravity
+    // values) on the ICRF rates W_ICRF/N_ICRF, and the frame rides separately:
+    // the Meeus polynomial p₀·T + T2F·T² in the snapshot branch (so its T² is
+    // Meeus-exact as before), the model's own precession of date (pFix) in the
+    // deep branch — ONE frame term on L′, ϖ and Ω alike.
+    const SENS = moonEccSensitivities;
+    const S_W = SENS.perigeeEffective, S_N = SENS.nodeEffective;
+    const S_W_PHYS = SENS.perigeePhysical, S_N_PHYS = SENS.nodePhysical;
+    const W_ICRF = SENS.perigeeRateIcrfDegPerCy, N_ICRF = SENS.nodeRateIcrfDegPerCy;
+    const P0 = SENS.generalPrecessionDegPerCy, T2F = SENS.frameT2DegPerCy2;
     const T2_W = S_W * WDOT * KAPPA / 2;   // −0.010318 °/cy² (Meeus ϖ −0.010320)
     const T2_N = S_N * NDOT * KAPPA / 2;   // +0.0020752 °/cy² (Meeus Ω +0.0020753)
     const EDDOT0 = eccentricityDotDotJ2000;
@@ -130,7 +147,7 @@ function createMoonArguments({ constants, fns }) {
     const T2_LP_TIDAL = (-25.86 / 3600) / 2;
     const T2_LP_PLANETARY = -0.0015786 - T2_LP_TIDAL;   // +7.247″/cy²
     const T2_LP = T2_LP_TIDAL + T2_LP_PLANETARY;
-    return { LP0, D0, M0, MP0, F0, LPR, DR, MR, WDOT, NDOT, T2_W, T2_N, T3_W, T3_N, T2_LP, T2_LP_TIDAL, S_W, S_N };
+    return { LP0, D0, M0, MP0, F0, LPR, DR, MR, WDOT, NDOT, T2_W, T2_N, T3_W, T3_N, T2_LP, T2_LP_TIDAL, S_W, S_N, S_W_PHYS, S_N_PHYS, W_ICRF, N_ICRF, P0, T2F };
   })();
 
   // ── Bounded Lp carriers (K_PL / C_OBL derived lazily — zero new constants)
@@ -235,6 +252,23 @@ function createMoonArguments({ constants, fns }) {
   }
   const obliquityCarrier = makeDeepCarrier(obliquityIntegrand, obliquityCarrierDense);
 
+  /** The frame share's REMAINDER for ϖ and Ω (plan 23). The physical exponents gave
+   *  up the frame's T² the Meeus-effective pair carried (T2F, Lieske's); pFix
+   *  supplies the model's own in-window acceleration (ṗΔ_eff/2, self-measured);
+   *  the difference rides the obliquity carrier's integrand, exactly as L′'s
+   *  figure+frame part does — so every of-date longitude's T² stays Meeus-exact
+   *  and ONE frame term (pFix + this remainder) rides L′, ϖ and Ω alike. Bounded
+   *  (the obliquity cycle), ~0.1° at most. Without the (d′) pair (pdSlope 0) the
+   *  whole Lieske share rides here and pFix is 0 — the same frame, one carrier.
+   *  @param {number} T @returns {number} */
+  function frameRemainder(T) {
+    if (T === 0) return 0;
+    const pdHalf = rateCompInit().pdSlopeDegPerCy2 / 2;
+    const t2Obl = (elpEarthFigureJ2ArcsecPerCy2 + elpGeneralPrecessionPA_T2ArcsecPerCy2) / 3600 - pdHalf;   // L′'s carrier normalization (oblInit)
+    const t2Rem = bundle.T2F - pdHalf;
+    return t2Obl === 0 ? 0 : obliquityCarrier(T) * (t2Rem / t2Obl);
+  }
+
   // ── Secular-ë completion carrier (ṅ campaign, plan §12i item 0 (c)) ──────
   // The certified snapshot's Lp T³/T⁴ tail is DERIVED physics (v4 D3,
   // tools/explore/v4-d3-tails.js): the Adams–Laplace channel's curvature
@@ -446,18 +480,22 @@ function createMoonArguments({ constants, fns }) {
     if (Ntrop === null || Naps === null || Nnod === null) return null;   // tidal-lock guard
     const Tj = (jd - j2000JD) / julianCenturyDays;
     const dy = y - argsY0;
-    // pFix rides Lp ONLY: the obliquity carrier's PA reduction compensates
-    // its ṗΔ·T²/2 there, keeping total Lp T² certified. The apsidal/nodal
-    // chains' own curvature is ALREADY certified-consistent (v4 checks;
-    // pre-(d′) isolation −0.20/+0.13″/cy²), so ϖ/Ω take only the tapered
-    // rate anchors — which self-measure and absorb their Δp rate share.
-    // Deep time: Mp/F inherit Lp's bounded pFix oscillation uncancelled
-    // (physically p-free args; bounded, display-class — documented).
-    const pF = rateCompInit().pFix(dy);
+    // ONE equinox for the three of-date longitudes (plan 23): pFix = ∫(p_dyn − p_kin),
+    // the model's own precession of date beyond the chains' kinematic mean, rides
+    // L′, ϖ AND Ω alike, so the frame-free arguments M′ = L′ − ϖ and F = L′ − Ω
+    // cancel it exactly at every epoch. The apsidal/nodal chains run the PHYSICAL
+    // exponents (the frame's T² is no longer frozen inside them); the Lieske share
+    // they gave up minus pFix's in-window T² rides the obliquity carrier as
+    // frameRemainder, as L′'s figure+frame part does, so ϖ and Ω stay Meeus-exact
+    // in T²; the tapered rate anchors self-measure. (Before plan 23 pFix rode L′ only while the chains
+    // carried the frame share in the Meeus-effective exponents: in the window the
+    // two agreed to 0.006″/cy², beyond it M′ and F inherited the whole wobble —
+    // measured 2.5° at ±10 kyr, 39° at −1 Myr.)
+    const pF = rateCompInit().pFix(dy), fr = frameRemainder(Tj);
     const Lp = A.LP0 + 360 * Ntrop + planetaryCarrier(Tj) + obliquityCarrier(Tj)
       + lpSecularCompletion(Tj) + pF;
-    const w = (A.LP0 - A.MP0) + 360 * Naps;                      // perigee ϖ (of-date, advance)
-    const om = (A.LP0 - A.F0) - 360 * Nnod;                      // node Ω (of-date, regression)
+    const w = (A.LP0 - A.MP0) + 360 * Naps + pF + fr;            // perigee ϖ (of-date, advance)
+    const om = (A.LP0 - A.F0) - 360 * Nnod + pF + fr;            // node Ω (of-date, regression)
     const Lsun = sunMeanLonDegAtJdTT(jd);          // the certified Sun's mean longitude (unwrapped)
     const ws = sunPerihelionLonDegAtJdTT(jd) + 180;   // the Sun's perihelion = Earth's mean ϖ of date + 180° (wrapped; M wraps)
     void A;
@@ -503,8 +541,13 @@ function createMoonArguments({ constants, fns }) {
     const Tc = Math.max(-100, Math.min(100, T));
     const Tc2 = Tc * Tc, Tc3 = Tc2 * Tc, Tc4 = Tc3 * Tc;
     const Lp = A.LP0 + A.LPR * T + A.T2_LP * Tc2 + Tc3 / 538841 - Tc4 / 65194000;
-    const w = (A.LP0 - A.MP0) + A.WDOT * (T + channelIntegral(T, A.S_W));   // perigee ϖ (of-date)
-    const om = (A.LP0 - A.F0) + A.NDOT * (T + channelIntegral(T, A.S_N));   // node Ω (of-date)
+    // Perigee and node (of-date), plan 23: the ICRF rate on the PHYSICAL channel
+    // exponent plus the Meeus polynomial's frame, p₀·T + T2F·T² (clamped with the
+    // other tails) — the T² is Meeus-exact as before, the frame share no longer
+    // rides the eccentricity channel's wiggle. The deep branch carries the model's
+    // own precession of date (pFix) in its place.
+    const w = (A.LP0 - A.MP0) + A.W_ICRF * (T + channelIntegral(T, A.S_W_PHYS)) + A.P0 * T + A.T2F * Tc2;   // perigee ϖ (of-date)
+    const om = (A.LP0 - A.F0) + A.N_ICRF * (T + channelIntegral(T, A.S_N_PHYS)) + A.P0 * T + A.T2F * Tc2;   // node Ω (of-date)
     const Lsun = sunMeanLonDegAtJdTT(jdTT);               // the certified Sun's mean longitude of date
     const ws = sunPerihelionLonDegAtJdTT(jdTT) + 180;     // the Sun's perihelion of date
     return {
@@ -520,12 +563,13 @@ function createMoonArguments({ constants, fns }) {
     const T2 = T * T, T3 = T2 * T, T4 = T3 * T;
     /** @param {number} x */
     const wrap = (x) => ((x % 360) + 360) % 360;
+    const R = MEEUS_ARGUMENT_RATES_DEG_PER_CY;   // the rates' ONE home (plan 23)
     return {
-      Lp: wrap(218.3164477 + 481267.88123421 * T - 0.0015786 * T2 + T3 / 538841 - T4 / 65194000),
-      D:  wrap(297.8501921 + 445267.1114034 * T - 0.0018819 * T2 + T3 / 545868 - T4 / 113065000),
-      M:  wrap(357.5291092 +  35999.0502909 * T - 0.0001536 * T2 + T3 / 24490000),
-      Mp: wrap(134.9633964 + 477198.8675055 * T + 0.0087414 * T2 + T3 / 69699 - T4 / 14712000),
-      F:  wrap( 93.2720950 + 483202.0175233 * T - 0.0036539 * T2 - T3 / 3526000 + T4 / 863310000),
+      Lp: wrap(218.3164477 + R.Lp * T - 0.0015786 * T2 + T3 / 538841 - T4 / 65194000),
+      D:  wrap(297.8501921 + R.D * T - 0.0018819 * T2 + T3 / 545868 - T4 / 113065000),
+      M:  wrap(357.5291092 + R.M * T - 0.0001536 * T2 + T3 / 24490000),
+      Mp: wrap(134.9633964 + R.Mp * T + 0.0087414 * T2 + T3 / 69699 - T4 / 14712000),
+      F:  wrap( 93.2720950 + R.F * T - 0.0036539 * T2 - T3 / 3526000 + T4 / 863310000),
     };
   }
 

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

@@ -12,30 +12,24 @@
  * — 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.
+ * DOMAIN (decision (ii)): 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). Earth's eccentricity of date is the
+ * N-body series on every other path as well (plan 21: one e-vector of date
+ * on every surface).
  *
- * 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.
+ * MEASURED at adoption (banked in the plan record): the attractor-mean
+ * modulation shifts were +2,526 ppm (perigee channel) and +1,067 ppm (node)
+ * against the single-line law it replaced — 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; 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):
+ * LOAD-BEARING conventions (the certified numerics):
  *  - channelIntegral uses composite Simpson with step ≤ ~4,000 yr
  *    (N ≥ 2, even);
  *  - modulation(tMa, s) takes age in Myr, positive = past, ≡ 1 at J2000;
@@ -63,7 +57,8 @@
  *
  * — ż_J2000 is the run's own z-rate at J2000 (script-written into the
  * artifact by generate.mjs from the banked series: "we are the source"),
- * T_q the H/12 taper the arguments' rate anchors ride. Inside the canon it
+ * T_q the arguments' secular-completion taper (SECULAR_COMPLETION_TAPER_YEARS,
+ * the one the rate anchors ride). Inside the canon it
  * restores the run's slope exactly; beyond ±T_q it is a constant vector
  * offset of ~3e-3, the size of the table's own spectrum-class scatter
  * there, never a growing polynomial. Absent the two fields (an older

package/src/moon/ecc-sensitivities.cjs

@@ -0,0 +1,87 @@
+/**
+ * The lunar perigee/node eccentricity sensitivities — ONE home (plan 23).
+ *
+ * The Sun's mean perturbation on the lunar perigee and node scales as
+ * g(e_E) = (1 − e_E²)^(−3/2); the of-date rates ride the channel
+ * rate(t) = rate₀ · [g(e_E(t))/g₀]^s. The pair the record fixes is the
+ * MEEUS-EFFECTIVE one: s_ϖ = 2.407 / s_Ω = 1.018 reproduce Meeus Ch. 47's
+ * of-date T² coefficients exactly (doc 66 §1, the (ė, s) degeneracy). Each of
+ * those T² coefficients contains the equinox's own acceleration — the general
+ * precession's T² term, a frame bookkeeping share that has nothing to do with
+ * Earth's eccentricity. Removing it gives the PHYSICAL pair the three-body
+ * laboratory reproduces from pure gravity (2.479 / 0.867; doc 66 §1):
+ *
+ *     s_phys = (s_eff · Ẋ_of-date · κ/2 − T2_frame) / (Ẋ_ICRF · κ/2),
+ *     κ = 3 e₀ ė₀ / (1 − e₀²),   Ẋ_ICRF = Ẋ_of-date − p₀
+ *
+ * — derived here from the injected constants, never typed. The runtime (plan
+ * 23) carries s_phys in the channel and ONE frame term, the model's own
+ * precession of date, on every of-date lunar longitude alike, so the
+ * frame-free arguments M′ = L′ − ϖ and F = L′ − Ω cancel it exactly at every
+ * epoch. Measured before the change (plan 23 §0): with the frame share frozen
+ * inside s_eff on the eccentricity channel, M′ and F inherited L′'s frame term
+ * uncancelled — 2.5° at ±10 kyr, 39° at −1 Myr — and the perigee and node
+ * carried a 3 % share of the channel's own wiggle that is not theirs; in the
+ * historical window the change is ≤ 21″ of ϖ/Ω at ±3000 (seconds of eclipse
+ * time).
+ *
+ * The Meeus Ch. 47 argument RATES live here too (the ONE set the arguments
+ * module, its pure-Meeus reference and this derivation share).
+ */
+
+'use strict';
+
+/** Meeus Ch. 47 mean-argument rates, degrees per Julian century (the J2000 polynomials' linear terms). */
+const MEEUS_ARGUMENT_RATES_DEG_PER_CY = Object.freeze({
+  Lp: 481267.88123421,   // the Moon's mean longitude L′
+  D: 445267.1114034,     // the mean elongation
+  M: 35999.0502909,      // the Sun's mean anomaly
+  Mp: 477198.8675055,    // the Moon's mean anomaly M′
+  F: 483202.0175233,     // the argument of latitude
+});
+
+/** The Meeus-effective pair: Meeus-exact of-date T² by construction (the shipped record values). */
+const MOON_ECC_SENSITIVITY_PERIGEE = 2.407;
+const MOON_ECC_SENSITIVITY_NODE = 1.018;
+
+/**
+ * The physical pair and the rate split, derived from the shipped constants.
+ * @param {{
+ *   eccentricityJ2000: number,
+ *   eccentricityDotJ2000: number,
+ *   frameT2ArcsecPerCy2: number,
+ *   generalPrecessionArcsecPerYr: number,
+ * }} inputs — e₀ (the catalog secular value the Taylor checks use); ė₀ per
+ *   Julian century; the general precession's T² coefficient the of-date Meeus
+ *   polynomials embed (Lieske 1976, the ELP frame), ″/cy²; p₀, the IAU 2006 p_A
+ *   linear rate, ″/yr
+ * @returns {Readonly<{
+ *   perigeeEffective: number, nodeEffective: number,
+ *   perigeePhysical: number, nodePhysical: number,
+ *   kappaPerCy: number, frameT2DegPerCy2: number, generalPrecessionDegPerCy: number,
+ *   perigeeRateOfDateDegPerCy: number, nodeRateOfDateDegPerCy: number,
+ *   perigeeRateIcrfDegPerCy: number, nodeRateIcrfDegPerCy: number,
+ * }>}
+ */
+function deriveMoonEccSensitivities({ eccentricityJ2000, eccentricityDotJ2000, frameT2ArcsecPerCy2, generalPrecessionArcsecPerYr }) {
+  const R = MEEUS_ARGUMENT_RATES_DEG_PER_CY;
+  const e0 = eccentricityJ2000, ed = eccentricityDotJ2000;
+  const kappaPerCy = 3 * e0 * ed / (1 - e0 * e0);
+  const generalPrecessionDegPerCy = generalPrecessionArcsecPerYr * 100 / 3600;
+  const frameT2DegPerCy2 = frameT2ArcsecPerCy2 / 3600;
+  const perigeeRateOfDateDegPerCy = R.Lp - R.Mp;   // ϖ̇ of date (+4069.0137) = ϖ̇_ICRF + p₀
+  const nodeRateOfDateDegPerCy = R.Lp - R.F;       // Ω̇ of date (−1934.1363) = Ω̇_ICRF + p₀
+  const perigeeRateIcrfDegPerCy = perigeeRateOfDateDegPerCy - generalPrecessionDegPerCy;
+  const nodeRateIcrfDegPerCy = nodeRateOfDateDegPerCy - generalPrecessionDegPerCy;
+  const half = kappaPerCy / 2;
+  const perigeePhysical = (MOON_ECC_SENSITIVITY_PERIGEE * perigeeRateOfDateDegPerCy * half - frameT2DegPerCy2) / (perigeeRateIcrfDegPerCy * half);
+  const nodePhysical = (MOON_ECC_SENSITIVITY_NODE * nodeRateOfDateDegPerCy * half - frameT2DegPerCy2) / (nodeRateIcrfDegPerCy * half);
+  return Object.freeze({
+    perigeeEffective: MOON_ECC_SENSITIVITY_PERIGEE, nodeEffective: MOON_ECC_SENSITIVITY_NODE,
+    perigeePhysical, nodePhysical,
+    kappaPerCy, frameT2DegPerCy2, generalPrecessionDegPerCy,
+    perigeeRateOfDateDegPerCy, nodeRateOfDateDegPerCy, perigeeRateIcrfDegPerCy, nodeRateIcrfDegPerCy,
+  });
+}
+
+module.exports = { MEEUS_ARGUMENT_RATES_DEG_PER_CY, MOON_ECC_SENSITIVITY_PERIGEE, MOON_ECC_SENSITIVITY_NODE, deriveMoonEccSensitivities };

package/src/moon/month-chain.cjs

@@ -33,10 +33,7 @@
  *    pair is 2.479/0.867 after removing the IAU frame term — doc 66 §1).
  *  - The of-date perigee/node periods are COMPOSITIONS of the inertial chain
  *    with the precession of date (the engines' `…OfDateCyclesBetween` count
- *    identities N_trop − N_anom and N_drac − N_trop). The former 'OfDate'
- *    family here — N·(H(t)/H₀)² counts per anchor unit and period ∝ T_yr/H —
- *    was the era device's counter law (doc 99's "Lunar Precession
- *    Invariant"); it left at plan 09 step 5c with the per-H counts.
+ *    identities N_trop − N_anom and N_drac − N_trop).
  *  - apsidalMeetsNodal / lunarLeveling existed ONLY in the browser before
  *    this move (the Node scene aliased apsidal-of-date on a cancellation
  *    argument valid for the paired scene nodes only) — the Node engine gains
@@ -70,13 +67,10 @@
  *   eccentricity channel's modulation; constants are the engine's derived
  *   J2000 anchors (both engines derive them from the same JSON).
  *   fns.precessionPeriodYearsAtAge is the engine's lunisolar precession period
- *   of date T_p(t), Julian years (the composed rate's period; plan 09 step
- *   5c — the tropical month's equinox turn, formerly the counter's 13/H(t),
- *   0.086 % slow: the chain's J2000 tropical month sat 5.85 ms above the IAU
- *   input, the 1/T_p form is within 0.1 ms of it).
+ *   of date T_p(t), Julian years (the composed rate's period — the tropical
+ *   month's equinox turn; within 0.1 ms of the IAU input at J2000).
  *   fns.distanceMetresAtAge is the engine's CITED distance history
- *   (deltat/moon-distance-history; plan 12 option A — the internal quartic
- *   twin this chain carried until then is deleted; null beyond the span).
+ *   (deltat/moon-distance-history; null beyond the span).
  */
 function createMoonMonthChain({ constants, fns }) {
   const {
@@ -136,8 +130,7 @@ function createMoonMonthChain({ constants, fns }) {
 
   /** Tropical month in seconds (equinox-referenced): the sidereal month less
    *  the equinox's own motion over one month — one turn per the precession
-   *  period of date T_p(t). (Plan 09 step 5c: formerly 13 turns per the
-   *  counter H(t), 0.086 % slow.) @param {number} tMa @returns {number | null} */
+   *  period of date T_p(t). @param {number} tMa @returns {number | null} */
   function tropicalMonthSecondsAtAge(tMa) {
     const tSm = siderealMonthSecondsAtAge(tMa);
     if (tSm === null) return null;
@@ -146,10 +139,8 @@ function createMoonMonthChain({ constants, fns }) {
     return tSm * (1 - tSm / (/** @type {number} */ (tP) * tYr));
   }
 
-  // (Plan 09 step 5c: the H² of-date family — apsidal/nodal "cycles per H"
-  // N·(H(t)/H₀)² and the of-date periods ∝ T_yr/H — is deleted with the per-H
-  // counts. The of-date perigee and node of date are the engines' count
-  // compositions of this chain's inertial periods with the tropical month.)
+  // The of-date perigee and node are the engines' count compositions of this
+  // chain's inertial periods with the tropical month.
 
   /** Perigee precession period in seconds — Brouwer–Clemence m² scaling ×
    *  e_E-channel modulation (the factored law: period = invariant mean /

package/src/moon/series-extension.cjs

@@ -52,14 +52,14 @@
  * A1-class correlated-subsample effect (the tracked events sit −1.6″
  * below the +1.4″ fleet mean and sample the near-annual arguments at
  * correlated phases; the era-split refutes anchor double-counting and
- * the modern era improves most, 3.57 → 3.34). The api gate (8″ then; 12″
- * since plan 06 R3 item 1 — the physical ΔT-convention floor) holds
+ * the modern era improves most, 3.57 → 3.34). The api gate (12″ — the
+ * physical ΔT-convention floor) holds
  * (max 6.2″). Coefficients are stored with the run→real argument
  * mapping already applied (the extraction's heliocentric-planet vs
  * geocentric-sun π-offset is folded into the signs).
  *
  * THE J2 NODE FAMILY (FQ-7 Round 3, plan §12i) — five derived terms on
- * Ω-family arguments (NOT expressible on the integer Delaunay lattice,
+ * Ω-family arguments (NOT expressible as integer combinations of the Delaunay arguments,
  * which is why every earlier extraction missed them; the 200-yr band
  * probe named them, the dense 2-day JPL sample confirmed them real, and
  * the channel isolation identified them as EARTH'S J2 acting on the

package/src/moon/series.cjs

@@ -11,8 +11,8 @@
  *
  * Full scene form: 60 longitude + 60 latitude terms with E/E² scaling on
  * M-bearing terms, A1 (Venus, Meeus-observed rate — the 18V−16E−M′
- * near-resonance makes it hypersensitive; no lattice identity), A2/A3
- * (D2-derived lattice rates, CHAIN-INTEGRATED under deep time through
+ * near-resonance makes it hypersensitive; no chain identity), A2/A3
+ * (D2-derived rates, CHAIN-INTEGRATED under deep time through
  * their identified content), all six latitude corrections, the EoC-half
  * subtraction (the off-centre orbit geometry already provides half the
  * equation of centre), and the two-term ellipse distance.
@@ -26,7 +26,7 @@
 
 'use strict';
 
-// Plan 06 R3 item 1: the DERIVED series-extension tail (the D2 planetary and
+// The DERIVED series-extension tail (the D2 planetary and
 // node terms, moon/series-extension.cjs) is part of the ONE scene Moon —
 // formerly added by the besselian tier alone, so the rendered Moon and the
 // eclipse-location Moon were two different Moons. The truncated finder forms
@@ -98,8 +98,6 @@ function createMoonSeries({ constants, fns }) {
       const nAps = apsidalOfDateCyclesBetween(yA0, yA);
       // A3 is the STAR-referenced lunar longitude (Meeus' rate 481,266.484 °/cy
       // ≡ 360°/sidereal month): the chain's sidereal-orbit count directly.
-      // Plan 09 step 5c — it was the tropical count minus the counter's 13/H
-      // precession turns (0.086 % slow on the precession).
       const nSid = siderealOrbitsBetween(yA0, yA);
       const nJ = jupiterOrbitsBetween(yA0, yA);
       if (nT !== null && nAps !== null && nSid !== null && nJ !== null) {
@@ -175,9 +173,8 @@ function createMoonSeries({ constants, fns }) {
    *  Returns everything the scene blocks write:
    *  thetaAddRad — the hierarchy θ increment (EoC-half-subtracted Σl);
    *  lonDeg — full ecliptic longitude (GEOMETRIC, mean equinox of date)
-   *  incl. the derived series-extension tail. (The fitted mean-longitude
-   *  anchor moonMeeusLpCorrection was retired to 0 at plan 06 R3 item 1 and
-   *  its parameter deleted at plan 11 step 2: measured against Horizons' apparent
+   *  incl. the derived series-extension tail. (No fitted mean-longitude
+   *  anchor — doc 109 §23: measured against Horizons' apparent
    *  Moon over 1970–2049 the series needs no anchor, −1.0″ ± 1.5″ flat in
    *  elongation; the former +32.75″ was the eclipse tier's mean Sun
    *  compensated — its aberration κ, the I2 long inequality and the trend-ΔT

package/types/deltat/deep-time-chain.d.ts

@@ -177,6 +177,19 @@ export function createDeepTimeChain(deps: {
         apsidalMeetsNodalSecondsAtAge: (tMa: number) => number | null;
         lunarLevelingSecondsAtAge: (tMa: number) => number | null;
     };
+    moonEccSensitivities: Readonly<{
+        perigeeEffective: number;
+        nodeEffective: number;
+        perigeePhysical: number;
+        nodePhysical: number;
+        kappaPerCy: number;
+        frameT2DegPerCy2: number;
+        generalPrecessionDegPerCy: number;
+        perigeeRateOfDateDegPerCy: number;
+        nodeRateOfDateDegPerCy: number;
+        perigeeRateIcrfDegPerCy: number;
+        nodeRateIcrfDegPerCy: number;
+    }>;
     meanMoonDistanceMetresAtAge: (tMa: number) => number | null;
     meanMoonDistanceAtAge: (tMa: number) => number | null;
     meanLodSecondsAtAge: (tMa: number) => number | null;
@@ -238,14 +251,6 @@ export function createDeepTimeChain(deps: {
     /** @param {number} tMa @param {number} periodJ2000Seconds @returns {number} */
     meanEarthPlanetSynodicPeriodAtAge: (tMa: number, periodJ2000Seconds: number) => number;
 }>;
-/**
- * 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 ΔT stack's switches. `resonator` gates the swing's δLOD terms; the
  * swing's ΔT term is gated by `resonatorDeltaT` when given (the Node twin

package/types/index.d.ts

@@ -2,6 +2,7 @@ export { createEpochPrimitives } from "./layer0/index.js";
 export { deriveEpochParams } from "./layer0/derive-params.js";
 export { createDerivedViews } from "./layer1/index.js";
 export { createDeepCalendar } from "./phase/deep-calendar.cjs";
+export { createDeepTimeChain } from "./deltat/deep-time-chain.js";
 export { createCardinalStructure } from "./cardinal/one-source-structure.cjs";
 export { createMoonMonthChain } from "./moon/month-chain.cjs";
 export { createChainCycleIntegrator } from "./chain-cycles/index.cjs";
@@ -61,7 +62,7 @@ export { GENERATED as DEFAULT_CONSTANTS, CONSTANTS_HASH, MODEL_VERSION, PREPRINT
 export { FITTED_COEFFICIENTS, COEFFICIENTS_HASH } from "./constants/index.js";
 export { GREGORIAN_START, GREGORIAN_START_JD, calendarToJD, jdToCalendar, jdToDecimalYear, dayFractionToHMS, julianCalendarToJD, jdToJulianCalendar } from "./time/calendar.cjs";
 export { DEG2RAD, RAD2DEG, ARCSEC_PER_TURN, ARCSEC_PER_RAD, wrap360, wrap180, wrap180Closed, wrap2pi, wrapPi } from "./math/angles.cjs";
-export { createDeepTimeChain, MOON_ECC_SENSITIVITY_PERIGEE, MOON_ECC_SENSITIVITY_NODE } from "./deltat/deep-time-chain.js";
+export { MOON_ECC_SENSITIVITY_PERIGEE, MOON_ECC_SENSITIVITY_NODE, deriveMoonEccSensitivities, MEEUS_ARGUMENT_RATES_DEG_PER_CY } from "./moon/ecc-sensitivities.cjs";
 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";
 export { computeEarthFrameOfDate, solveWheelAngleForLongitude, computeSceneWorldFrameJ2000, SCENE_EQUINOX_VIEW_AZIMUTH_DEG } from "./earth/frame-of-date.cjs";

package/types/model/deep-time-chain.d.ts

@@ -153,6 +153,19 @@ export function assembleDeepTimeChain(ctx: import("./derived-constants.js").Ctx)
             apsidalMeetsNodalSecondsAtAge: (tMa: number) => number | null;
             lunarLevelingSecondsAtAge: (tMa: number) => number | null;
         };
+        moonEccSensitivities: Readonly<{
+            perigeeEffective: number;
+            nodeEffective: number;
+            perigeePhysical: number;
+            nodePhysical: number;
+            kappaPerCy: number;
+            frameT2DegPerCy2: number;
+            generalPrecessionDegPerCy: number;
+            perigeeRateOfDateDegPerCy: number;
+            nodeRateOfDateDegPerCy: number;
+            perigeeRateIcrfDegPerCy: number;
+            nodeRateIcrfDegPerCy: number;
+        }>;
         meanMoonDistanceMetresAtAge: (tMa: number) => number | null;
         meanMoonDistanceAtAge: (tMa: number) => number | null;
         meanLodSecondsAtAge: (tMa: number) => number | null;

package/types/model/lunar-theory.d.ts

@@ -18,9 +18,9 @@ export function assembleLunarTheory(ctx: import("./time-axis.js").Ctx): {
             M0: number;
             MP0: number;
             F0: number;
-            LPR: number;
-            DR: number;
-            MR: number;
+            LPR: 481267.88123421;
+            DR: 445267.1114034;
+            MR: 35999.0502909;
             WDOT: number;
             NDOT: number;
             T2_W: number;
@@ -31,6 +31,12 @@ export function assembleLunarTheory(ctx: import("./time-axis.js").Ctx): {
             T2_LP_TIDAL: number;
             S_W: number;
             S_N: number;
+            S_W_PHYS: number;
+            S_N_PHYS: number;
+            W_ICRF: number;
+            N_ICRF: number;
+            P0: number;
+            T2F: number;
         };
         _rateCompletion: () => {
             rc: {

package/types/model/one-source-movement.d.ts

@@ -33,7 +33,7 @@ export function assembleOneSourceMovement(ctx: import("./year-day-lengths.js").C
         periOfDateDegAt: (year: number) => number;
         /** 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: number) => number;
-        /** 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: [number, number] | null;
         /** 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: number) => {

package/types/model/surface.d.ts

@@ -44,7 +44,7 @@ export function assembleSurface(ctx: import("./cardinal-instants.js").Ctx): Read
     lunisolar: Readonly<{
         /** 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: number) => 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} */
+        /** 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: number) => number;
         /** 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: number) => number;
@@ -123,9 +123,9 @@ export function assembleSurface(ctx: import("./cardinal-instants.js").Ctx): Read
     earth: Readonly<{
         perihelionLongitudeDeg: (year: number) => number;
         obliquityDeg: (year: number) => number;
-        /** 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: [number, number] | null;
-        /** 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: (year: number) => number;
         eccentricity: (year: number) => number;
         inclinationDeg: (year: number) => number;
@@ -184,6 +184,19 @@ export function assembleSurface(ctx: import("./cardinal-instants.js").Ctx): Read
         distanceKmAtYear: (year: number) => number;
         siderealMonthDaysAtYear: (year: number) => number;
         synodicMonthDays: number;
+        eccSensitivities: Readonly<{
+            perigeeEffective: number;
+            nodeEffective: number;
+            perigeePhysical: number;
+            nodePhysical: number;
+            kappaPerCy: number;
+            frameT2DegPerCy2: number;
+            generalPrecessionDegPerCy: number;
+            perigeeRateOfDateDegPerCy: number;
+            nodeRateOfDateDegPerCy: number;
+            perigeeRateIcrfDegPerCy: number;
+            nodeRateIcrfDegPerCy: number;
+        }>;
         lonDegAtJD: (jd: number) => number;
         betaDegAtJD: (jd: number) => number;
         distanceKmAtJD: (jd: number) => number;
@@ -191,7 +204,7 @@ export function assembleSurface(ctx: import("./cardinal-instants.js").Ctx): Read
     eclipse: Readonly<{
         sunLonDegAtJD: (jd: number) => number;
         sunLonCompletedDegAtJD: (jd: number) => 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 (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: (jdUT: number) => number;
         findLunarInRange: (jdStart: number, jdEnd: number) => {
             jd: number;

package/types/model.d.ts

@@ -145,6 +145,19 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         distanceKmAtYear: (year: number) => number;
         siderealMonthDaysAtYear: (year: number) => number;
         synodicMonthDays: number;
+        eccSensitivities: Readonly<{
+            perigeeEffective: number;
+            nodeEffective: number;
+            perigeePhysical: number;
+            nodePhysical: number;
+            kappaPerCy: number;
+            frameT2DegPerCy2: number;
+            generalPrecessionDegPerCy: number;
+            perigeeRateOfDateDegPerCy: number;
+            nodeRateOfDateDegPerCy: number;
+            perigeeRateIcrfDegPerCy: number;
+            nodeRateIcrfDegPerCy: number;
+        }>;
         lonDegAtJD: (jd: number) => number;
         betaDegAtJD: (jd: number) => number;
         distanceKmAtJD: (jd: number) => number;

package/types/moon/arguments.d.cts

@@ -16,6 +16,7 @@ export type MoonArgsDeg = {
  *     elpEarthFigureJ2ArcsecPerCy2: number,
  *     elpGeneralPrecessionPA_T2ArcsecPerCy2: number,
  *     eccE0: number,
+ *     moonEccSensitivities: ReturnType<typeof import('./ecc-sensitivities.cjs').deriveMoonEccSensitivities>,
  *   },
  *   fns: {
  *     eccAt: (tYr: number) => number,
@@ -56,6 +57,7 @@ export function createMoonArguments({ constants, fns }: {
         elpEarthFigureJ2ArcsecPerCy2: number;
         elpGeneralPrecessionPA_T2ArcsecPerCy2: number;
         eccE0: number;
+        moonEccSensitivities: ReturnType<typeof import("./ecc-sensitivities.cjs").deriveMoonEccSensitivities>;
     };
     fns: {
         eccAt: (tYr: number) => number;
@@ -86,9 +88,9 @@ export function createMoonArguments({ constants, fns }: {
         M0: number;
         MP0: number;
         F0: number;
-        LPR: number;
-        DR: number;
-        MR: number;
+        LPR: 481267.88123421;
+        DR: 445267.1114034;
+        MR: 35999.0502909;
         WDOT: number;
         NDOT: number;
         T2_W: number;
@@ -99,6 +101,12 @@ export function createMoonArguments({ constants, fns }: {
         T2_LP_TIDAL: number;
         S_W: number;
         S_N: number;
+        S_W_PHYS: number;
+        S_N_PHYS: number;
+        W_ICRF: number;
+        N_ICRF: number;
+        P0: number;
+        T2F: number;
     };
     _rateCompletion: () => {
         rc: {

package/types/moon/deep-ecc-channel.d.cts

@@ -17,7 +17,8 @@
  *
  * — ż_J2000 is the run's own z-rate at J2000 (script-written into the
  * artifact by generate.mjs from the banked series: "we are the source"),
- * T_q the H/12 taper the arguments' rate anchors ride. Inside the canon it
+ * T_q the arguments' secular-completion taper (SECULAR_COMPLETION_TAPER_YEARS,
+ * the one the rate anchors ride). Inside the canon it
  * restores the run's slope exactly; beyond ±T_q it is a constant vector
  * offset of ~3e-3, the size of the table's own spectrum-class scatter
  * there, never a growing polynomial. Absent the two fields (an older

package/types/moon/ecc-sensitivities.d.cts

@@ -0,0 +1,48 @@
+/** Meeus Ch. 47 mean-argument rates, degrees per Julian century (the J2000 polynomials' linear terms). */
+export const MEEUS_ARGUMENT_RATES_DEG_PER_CY: Readonly<{
+    Lp: 481267.88123421;
+    D: 445267.1114034;
+    M: 35999.0502909;
+    Mp: 477198.8675055;
+    F: 483202.0175233;
+}>;
+/** The Meeus-effective pair: Meeus-exact of-date T² by construction (the shipped record values). */
+export const MOON_ECC_SENSITIVITY_PERIGEE: 2.407;
+export const MOON_ECC_SENSITIVITY_NODE: 1.018;
+/**
+ * The physical pair and the rate split, derived from the shipped constants.
+ * @param {{
+ *   eccentricityJ2000: number,
+ *   eccentricityDotJ2000: number,
+ *   frameT2ArcsecPerCy2: number,
+ *   generalPrecessionArcsecPerYr: number,
+ * }} inputs — e₀ (the catalog secular value the Taylor checks use); ė₀ per
+ *   Julian century; the general precession's T² coefficient the of-date Meeus
+ *   polynomials embed (Lieske 1976, the ELP frame), ″/cy²; p₀, the IAU 2006 p_A
+ *   linear rate, ″/yr
+ * @returns {Readonly<{
+ *   perigeeEffective: number, nodeEffective: number,
+ *   perigeePhysical: number, nodePhysical: number,
+ *   kappaPerCy: number, frameT2DegPerCy2: number, generalPrecessionDegPerCy: number,
+ *   perigeeRateOfDateDegPerCy: number, nodeRateOfDateDegPerCy: number,
+ *   perigeeRateIcrfDegPerCy: number, nodeRateIcrfDegPerCy: number,
+ * }>}
+ */
+export function deriveMoonEccSensitivities({ eccentricityJ2000, eccentricityDotJ2000, frameT2ArcsecPerCy2, generalPrecessionArcsecPerYr }: {
+    eccentricityJ2000: number;
+    eccentricityDotJ2000: number;
+    frameT2ArcsecPerCy2: number;
+    generalPrecessionArcsecPerYr: number;
+}): Readonly<{
+    perigeeEffective: number;
+    nodeEffective: number;
+    perigeePhysical: number;
+    nodePhysical: number;
+    kappaPerCy: number;
+    frameT2DegPerCy2: number;
+    generalPrecessionDegPerCy: number;
+    perigeeRateOfDateDegPerCy: number;
+    nodeRateOfDateDegPerCy: number;
+    perigeeRateIcrfDegPerCy: number;
+    nodeRateIcrfDegPerCy: number;
+}>;

package/types/moon/month-chain.d.cts

@@ -23,13 +23,10 @@
  *   eccentricity channel's modulation; constants are the engine's derived
  *   J2000 anchors (both engines derive them from the same JSON).
  *   fns.precessionPeriodYearsAtAge is the engine's lunisolar precession period
- *   of date T_p(t), Julian years (the composed rate's period; plan 09 step
- *   5c — the tropical month's equinox turn, formerly the counter's 13/H(t),
- *   0.086 % slow: the chain's J2000 tropical month sat 5.85 ms above the IAU
- *   input, the 1/T_p form is within 0.1 ms of it).
+ *   of date T_p(t), Julian years (the composed rate's period — the tropical
+ *   month's equinox turn; within 0.1 ms of the IAU input at J2000).
  *   fns.distanceMetresAtAge is the engine's CITED distance history
- *   (deltat/moon-distance-history; plan 12 option A — the internal quartic
- *   twin this chain carried until then is deleted; null beyond the span).
+ *   (deltat/moon-distance-history; null beyond the span).
  */
 export function createMoonMonthChain({ constants, fns }: {
     constants: {