@essrt/physics 13.0.0 → 15.0.0

package/package.json

@@ -1,9 +1,9 @@
 {
   "name": "@essrt/physics",
-  "version": "13.0.0",
+  "version": "15.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": "v17.0",
+    "modelVersion": "v19.0",
     "_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",
@@ -50,10 +50,7 @@
     "./moon/series-extension": { "types": "./types/moon/series-extension.d.cts", "default": "./src/moon/series-extension.cjs" },
     "./moon/apparent": { "types": "./types/moon/apparent.d.cts", "default": "./src/moon/apparent.cjs" },
     "./planets/geometry": { "types": "./types/planets/geometry.d.cts", "default": "./src/planets/geometry.cjs" },
-    "./planets/fibonacci-laws": { "types": "./types/planets/fibonacci-laws.d.cts", "default": "./src/planets/fibonacci-laws.cjs" },
     "./planets/ecc-channel": { "types": "./types/planets/ecc-channel.d.cts", "default": "./src/planets/ecc-channel.cjs" },
-    "./planets/orientation": { "types": "./types/planets/orientation.d.cts", "default": "./src/planets/orientation.cjs" },
-    "./planets/asc-node-integrator": { "types": "./types/planets/asc-node-integrator.d.cts", "default": "./src/planets/asc-node-integrator.cjs" },
     "./planets/orbit-chain": { "types": "./types/planets/orbit-chain.d.cts", "default": "./src/planets/orbit-chain.cjs" },
     "./planets/keplerian-chain": { "types": "./types/planets/keplerian-chain.d.cts", "default": "./src/planets/keplerian-chain.cjs" },
     "./planets/spin-channel": { "types": "./types/planets/spin-channel.d.cts", "default": "./src/planets/spin-channel.cjs" },

package/src/constants/generated.d.ts

@@ -2,14 +2,14 @@
 // Gives the TypeScript website full type safety at the boundary (§2g) while
 // packages/physics stays JavaScript.
 
-export declare const CONSTANTS_HASH: "726b10e23d8ace09";
+export declare const CONSTANTS_HASH: "5ce921f2efa2e845";
 
 export declare const MODEL_VERSION: string;
 
 export declare const PREPRINT_DOI: string;
 
 export declare const DEFAULT_CONSTANTS: {
-  readonly hash: "726b10e23d8ace09";
+  readonly hash: "5ce921f2efa2e845";
   readonly additionalBodies: {
     "pluto": {
       "name": string;
@@ -17,9 +17,7 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "angleCorrection": number;
       "startpos": number;
-      "perihelionEclipticFraction": [number, number];
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
       "invPlaneInclinationMean": number;
       "invPlaneInclinationAmplitude": number;
     };
@@ -29,9 +27,7 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "angleCorrection": number;
       "startpos": number;
-      "perihelionEclipticFraction": [number, number];
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
       "invPlaneInclinationMean": number;
       "invPlaneInclinationAmplitude": number;
     };
@@ -41,9 +37,7 @@ export declare const DEFAULT_CONSTANTS: {
       "type": string;
       "angleCorrection": number;
       "startpos": number;
-      "perihelionEclipticFraction": [number, number];
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
       "invPlaneInclinationMean": number;
       "invPlaneInclinationAmplitude": number;
     };
@@ -53,9 +47,7 @@ export declare const DEFAULT_CONSTANTS: {
       "orbitDistanceOverride": number;
       "angleCorrection": number;
       "startpos": number;
-      "perihelionEclipticFraction": [number, number];
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
       "invPlaneInclinationMean": number;
       "invPlaneInclinationAmplitude": number;
     };
@@ -486,13 +478,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "venus": {
@@ -500,13 +488,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "mars": {
@@ -514,13 +498,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "jupiter": {
@@ -528,13 +508,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "saturn": {
@@ -542,13 +518,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "uranus": {
@@ -556,13 +528,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
     "neptune": {
@@ -570,13 +538,9 @@ export declare const DEFAULT_CONSTANTS: {
       "eocFraction": number;
       "startpos": number;
       "angleCorrection": number;
-      "perihelionEclipticFraction": [number, number];
       "type": string;
       "mirrorPair": string;
-      "fibonacciD": number;
       "ascendingNodeInvPlane": number;
-      "inclinationCycleAnchor": number;
-      "antiPhase": boolean;
       "ascendingNodeCyclesIn8H": number;
     };
   };

package/src/constants/generated.js

@@ -33,17 +33,17 @@
  * carry it so a counterfactual is reproducible (§2d).
  * @type {string}
  */
-export const CONSTANTS_HASH = "726b10e23d8ace09";
+export const CONSTANTS_HASH = "5ce921f2efa2e845";
 
 /** Model version label — single source: public/input/model-version.json (§10 two-axis scheme). */
-export const MODEL_VERSION = "v17.0";
+export const MODEL_VERSION = "v19.0";
 
 /** Canonical preprint DOI — single source: public/input/model-version.json. */
 export const PREPRINT_DOI = "10.21203/rs.3.rs-8758810/v4";
 
 /** @type {Readonly<Record<string, unknown>>} */
 export const DEFAULT_CONSTANTS = Object.freeze({
-  hash: "726b10e23d8ace09",
+  hash: "5ce921f2efa2e845",
   additionalBodies: {
     "pluto": {
       "name": "Pluto",
@@ -51,12 +51,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "type": "I",
       "angleCorrection": 2.469281,
       "startpos": 71.555,
-      "perihelionEclipticFraction": [
-        1,
-        1
-      ],
       "ascendingNodeInvPlane": 101.06,
-      "inclinationCycleAnchor": 203.32,
       "invPlaneInclinationMean": 15.7162,
       "invPlaneInclinationAmplitude": 0.717024
     },
@@ -66,12 +61,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "type": "III",
       "angleCorrection": -1.619816,
       "startpos": 80,
-      "perihelionEclipticFraction": [
-        1,
-        1
-      ],
       "ascendingNodeInvPlane": 59.56,
-      "inclinationCycleAnchor": 23.3195,
       "invPlaneInclinationMean": 150,
       "invPlaneInclinationAmplitude": 0.1
     },
@@ -81,12 +71,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "type": "II",
       "angleCorrection": 0.047888,
       "startpos": 57.402,
-      "perihelionEclipticFraction": [
-        1,
-        1
-      ],
       "ascendingNodeInvPlane": 10.36,
-      "inclinationCycleAnchor": 203.3195,
       "invPlaneInclinationMean": 9.25,
       "invPlaneInclinationAmplitude": 0.5
     },
@@ -96,12 +81,7 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "orbitDistanceOverride": 2.76596,
       "angleCorrection": 0,
       "startpos": 0,
-      "perihelionEclipticFraction": [
-        1,
-        1
-      ],
       "ascendingNodeInvPlane": 80.89,
-      "inclinationCycleAnchor": 203.3195,
       "invPlaneInclinationMean": 0.43,
       "invPlaneInclinationAmplitude": 0.05
     }
@@ -647,16 +627,9 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "eocFraction": -0.527,
       "startpos": 83.64651142941864,
       "angleCorrection": 0.9727205912043573,
-      "perihelionEclipticFraction": [
-        8,
-        11
-      ],
       "type": "I",
       "mirrorPair": "uranus",
-      "fibonacciD": 21,
       "ascendingNodeInvPlane": 32.83,
-      "inclinationCycleAnchor": 234.52,
-      "antiPhase": false,
       "ascendingNodeCyclesIn8H": 9
     },
     "venus": {
@@ -664,16 +637,9 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "eocFraction": 0.436,
       "startpos": 249.28769121856186,
       "angleCorrection": -2.789778861448692,
-      "perihelionEclipticFraction": [
-        -8,
-        6
-      ],
       "type": "I",
       "mirrorPair": "neptune",
-      "fibonacciD": 34,
       "ascendingNodeInvPlane": 54.7,
-      "inclinationCycleAnchor": 218.64,
-      "antiPhase": false,
       "ascendingNodeCyclesIn8H": 1
     },
     "mars": {
@@ -681,16 +647,9 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "eocFraction": -0.066224,
       "startpos": 121.4634461571797,
       "angleCorrection": -2.104652936098656,
-      "perihelionEclipticFraction": [
-        8,
-        36
-      ],
       "type": "II",
       "mirrorPair": "jupiter",
-      "fibonacciD": 5,
       "ascendingNodeInvPlane": 354.87,
-      "inclinationCycleAnchor": 236.07,
-      "antiPhase": false,
       "ascendingNodeCyclesIn8H": 64
     },
     "jupiter": {
@@ -698,16 +657,9 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "eocFraction": 0.495,
       "startpos": 13.887946755371052,
       "angleCorrection": 0.9333205235053406,
-      "perihelionEclipticFraction": [
-        8,
-        39
-      ],
       "type": "III",
       "mirrorPair": "mars",
-      "fibonacciD": 5,
       "ascendingNodeInvPlane": 312.89,
-      "inclinationCycleAnchor": 287.06,
-      "antiPhase": false,
       "ascendingNodeCyclesIn8H": 36
     },
     "saturn": {
@@ -715,16 +667,9 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "eocFraction": 0.54,
       "startpos": 11.279308392461772,
       "angleCorrection": -0.1784695052193328,
-      "perihelionEclipticFraction": [
-        -8,
-        65
-      ],
       "type": "III",
       "mirrorPair": "earth",
-      "fibonacciD": 3,
       "ascendingNodeInvPlane": 118.81,
-      "inclinationCycleAnchor": 116.26,
-      "antiPhase": true,
       "ascendingNodeCyclesIn8H": 36
     },
     "uranus": {
@@ -732,16 +677,9 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "eocFraction": 0.53,
       "startpos": 44.900388945775106,
       "angleCorrection": -0.7343433510960402,
-      "perihelionEclipticFraction": [
-        1,
-        3
-      ],
       "type": "III",
       "mirrorPair": "mercury",
-      "fibonacciD": 21,
       "ascendingNodeInvPlane": 307.8,
-      "inclinationCycleAnchor": 21.33,
-      "antiPhase": false,
       "ascendingNodeCyclesIn8H": 11
     },
     "neptune": {
@@ -749,16 +687,9 @@ export const DEFAULT_CONSTANTS = Object.freeze({
       "eocFraction": 0.585,
       "startpos": 47.9552024382285,
       "angleCorrection": 2.3327665835675617,
-      "perihelionEclipticFraction": [
-        2,
-        1
-      ],
       "type": "III",
       "mirrorPair": "venus",
-      "fibonacciD": 34,
       "ascendingNodeInvPlane": 192.04,
-      "inclinationCycleAnchor": 174.04,
-      "antiPhase": false,
       "ascendingNodeCyclesIn8H": 3
     }
   },

package/src/index.js

@@ -122,12 +122,13 @@ export { createMoonApparent } from './moon/apparent.cjs';
 // sample) and the offset-circle wheel-angle solver: the scene twins' one placement.
 export { computeEarthFrameOfDate, solveWheelAngleForLongitude } from './earth/frame-of-date.cjs';
 // Phase 8.3 — the planet machinery, extracted by LAW over body records
-// (survey order: geometry → Fibonacci laws → channels → chains).
+// (survey order: geometry → channels → chains). The ψ/K amplitude laws
+// (planets/fibonacci-laws), the record-based orientation forms
+// (planets/orientation) and the ascending-node integrator
+// (planets/asc-node-integrator) are DELETED at plan 07 R5/R6: the planets'
+// elements of date have one home, the N-body chain; docs/retired-record.md.
 export { derivePlanetGeometry } from './planets/geometry.cjs';
-export * as planetFibonacciLaws from './planets/fibonacci-laws.cjs';
 export { eccentricityFromCycles, computeEccentricityIntegrated } from './planets/ecc-channel.cjs';
-export * as planetOrientation from './planets/orientation.cjs';
-export { integrateAscendingNode } from './planets/asc-node-integrator.cjs';
 export * as planetOrbitChain from './planets/orbit-chain.cjs';
 // P5/K4.6c — the engine-D Keplerian chain (pure evaluator + the embedded
 // governed artifact; the browser flag path consumes exactly these).

package/src/model.js

@@ -18,7 +18,6 @@
  * which itself mirrors tools/lib/constants.js §9.
  */
 import { deriveEpochParams } from './layer0/derive-params.js';
-import * as FL from './planets/fibonacci-laws.cjs';
 import { buildPlanetChainsFromArtifactData, computePlanetElementsAtYear } from './planets/keplerian-chain.cjs';
 import { createSecularSeriesOverride } from './planets/secular-series.cjs';
 import { computeEquatorNodeOriginSFrameDeg, convertNodeSFrameToEquatorOriginDeg } from './planets/inv-plane-frame.cjs';
@@ -26,7 +25,6 @@ import { createPhaseMachinery } from './phase/index.cjs';
 import { createYearLengths, ONE_FAMILY_WINDOW_YEARS } from './earth/year-lengths.cjs';
 import { createDeepOrbitalHistory } from './earth/deep-orbital-history.cjs';
 import { CHAIN_ARTIFACT } from './planets/chain-artifact.js';
-import { computeSecularShape } from './planets/secular-shape.cjs';
 import { createPlanetSpinChannelFromArtifacts, computeObliquityJ2000Deg } from './planets/spin-channel.cjs';
 import { createDeltaTCycles } from './deltat/cycles.cjs';
 import { createDeepTimeLod } from './deltat/deep-time.cjs';
@@ -97,14 +95,12 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   const meanAnomalisticYearDays = (meanSolarYearDays * (H / 16)) / (H / 16 - 1);
   const meanTropicalYearJ2000Seconds = meanSolarYearDays * meanLengthOfDay;
 
-  const earthInclAmplitude = C.earth.earthInvPlaneInclinationAmplitude;
-  // Plan 07 R6: `earthtiltMean` and `eccentricityAmplitude` were read here
-  // ONLY to invert the K constant from Earth's calibration. With the K law
-  // retired nothing in this module reads either — they remain on C.earth for
-  // the Earth-side devices (R9) and the registry.
-  const earthInclMean = C.earthOrbital.earthInclinationJ2000_deg
-    - earthInclAmplitude * Math.cos(((C.earthOrbital.earthPerihelionLongitudeJ2000
-      - C.earthOrbital.earthInclinationCycleAnchor) * Math.PI) / 180);
+  // Plan 07 R5/R6/R9: `earthtiltMean`, `earthInvPlaneInclinationAmplitude`
+  // and `eccentricityAmplitude` were read here only to invert the K and ψ
+  // constants from Earth's calibration and to build the H/3 inclination
+  // device. With both laws retired nothing in this module reads them — they
+  // remain on C.earth as the K device's J2000 pose inputs (the scene's wheel
+  // geometry at the chain anchor epoch) and for the registry's record keys.
   const solsticeObliquityMean = F.SOLSTICE_OBLIQUITY_MEAN_FITTED;
 
   const G_CONSTANT = C.physicalConstants.G_CONSTANT;
@@ -429,14 +425,9 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   /** de/dyear of the Sun's e — the cardinal braid's equation-of-centre
    *  derivative rides it (±0.5-yr central difference; the constant offset
    *  cancels). @param {number} year @returns {number} */
-  /** @param {number} year @returns {number} */
-  const inclinationDeg = (year) => earthInclMean
-    - earthInclAmplitude * Math.cos(phaseRadians(balancedYear, year, 3));
-  /** @param {number} year @returns {number} */
-  const ascendingNodeDeg = (year) => {
-    const period = -H / 5;
-    return (((C.earthOrbital.earthAscendingNodeInvPlane + (360.0 * (year - 2000)) / period) % 360) + 360) % 360;
-  };
+  // (Plan 07 R5/R9: Earth's H/3 inclination cosine and −H/5 linear node — the
+  // K device — stood here; `inclinationDeg`/`ascendingNodeDeg` are now the
+  // one-source construction defined with the planets' chain evaluators below.)
 
   // ── Year/day lengths ──────────────────────────────────────────────────────
   /** @param {number} year @param {number} base @param {Array<[number, number, number]>} harmonics @returns {number} */
@@ -617,6 +608,8 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
     });
     return {
       yearLengths,
+      /** The one-source sample at a decimal year (the engine year) — ε, e, ϖ of date, the Earth-orbit normal (the fields the frame of date and the invariable-plane elements read). @param {number} year */
+      sampleAt,
       /** EXPERIMENT accessors: the series' e and equinox-referenced ϖ of date. @param {number} year @returns {number} */
       eAt: (year) => sampleAt(year).e,
       /** @param {number} year @returns {number} */
@@ -716,41 +709,27 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   for (const k of PLANET_KEYS) massFraction[k] = 1 / C.physicalConstants.massRatioDE440[k];
   massFraction.earth = massEarthAlone / M_SUN;
 
-  const PSI = FL.computePsiConstant({
-    earthInvPlaneInclinationAmplitude: earthInclAmplitude,
-    massEarthAlone,
-    massSun: M_SUN,
-  });
   // (plan 07 R6: eccentricityAmplitudeK — the K constant inverted from
   // Earth's calibration — and the t2000 phase offset it fed went with the
-  // eccentricity law. systemResetN is kept where other consumers read it.)
-
-  /** @param {[number, number]|null} frac @returns {number|null} */
-  const fractionToYears = (frac) => (frac === null ? null : (H * frac[0]) / frac[1]);
+  // eccentricity law; plan 07 R5: PSI — ψ = 3·A_earth·√(m_E/m_☉) — and the
+  // inclination law it fed went the same way. systemResetN is kept where
+  // other consumers read it.)
 
-  /** @type {Record<string, Record<string, any>>} */
-  const PLANET_RECORDS = {};              // internal: the fitters read the device fields
+  // Plan 07 R1/R8: the internal PLANET_RECORDS — the device record the
+  // fitters once read (the H·num/den perihelion period from
+  // perihelionEclipticFraction, fibonacciD, antiPhase, inclinationCycleAnchor,
+  // the −8H/N node period) — had no reader left after R6/R5 and is gone with
+  // its JSON inputs. What remains is the PUBLISHED record below.
   /** @type {Record<string, Record<string, any>>} */
   const PUBLISHED_PLANET_RECORDS = {};    // what model.planets.record() serves (R8)
   for (const k of PLANET_KEYS) {
     const mp = C.planets[k];
     const ar = C.planetOrbitalElements[k];
-    const ecl = /** @type {number} */ (fractionToYears(mp.perihelionEclipticFraction));
     // Plan 06 Phase 7 commit 2: the K device's integer axial and obliquity
-    // fractions are retired. The eccentricity law's cycle period is the
-    // chain's OWN g-mode beat (the dominant mode × largest companion of the
-    // planet's eccentricity vector), its obliquity input the DERIVED J2000
-    // obliquity of the spin channel (the IAU pole against the chain's J2000
-    // plane, acute form) — both from the governed artifacts, no fractions.
-    const wobble = computeSecularShape(/** @type {any} */ (CHAIN_ARTIFACT), k).beatYears;
-    const il = FL.computeInclinationLaw({
-      fibonacciD: mp.fibonacciD,
-      massFrac: massFraction[k],
-      invPlaneInclinationJ2000: ar.invPlaneInclinationJ2000,
-      longitudePerihelion: ar.longitudePerihelion,
-      inclinationCycleAnchor: mp.inclinationCycleAnchor,
-      antiPhase: mp.antiPhase || false,
-    }, PSI);
+    // fractions are retired; the obliquity is the DERIVED J2000 obliquity of
+    // the spin channel (the IAU pole against the chain's J2000 plane, acute
+    // form) — from the governed artifacts, no fractions. (The g-mode beat the
+    // K law's period rode left with the internal record at plan 07 R1/R8.)
     const obliquityDerived = computeObliquityJ2000Deg({
       spin: C.planetSpinPhysical[k],
       anchorInclEclipticDeg: /** @type {any} */ (CHAIN_ARTIFACT).j2000AnchorElements[k].inclEclipticDeg,
@@ -758,31 +737,11 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
       obliquityJ2000Deg: C.earthOrbital.obliquityJ2000_deg,
     });
     const obliquityMean = Math.min(obliquityDerived, 180 - obliquityDerived);   // acute form of the derived J2000 obliquity
-    PLANET_RECORDS[k] = Object.freeze({
-      name: mp.name,
-      perihelionEclipticYears: ecl,
-      longitudePerihelion: ar.longitudePerihelion,
-      ascendingNodeCyclesIn8H: mp.ascendingNodeCyclesIn8H,
-      ascendingNodePeriod: -(8 * H) / mp.ascendingNodeCyclesIn8H,
-      wobblePeriod: wobble,
-      fibonacciD: mp.fibonacciD,
-      antiPhase: mp.antiPhase || false,
-      ascendingNodeInvPlane: mp.ascendingNodeInvPlane,
-      inclinationCycleAnchor: mp.inclinationCycleAnchor,
-      invPlaneInclinationJ2000: ar.invPlaneInclinationJ2000,
-      invPlaneInclinationAmplitude: il.amplitude,
-      invPlaneInclinationMean: il.mean,
-      obliquityMean,
-      orbitalEccentricityJ2000: ar.orbitalEccentricityJ2000,
-      solarYearInput: ar.solarYearInput,
-      axialTiltJ2000: ar.axialTiltJ2000,
-    });
-    // R8 (plan 07): what `model.planets.record()` PUBLISHES is not the record
-    // above. The internal record kept the ψ/K law outputs for the fitting
-    // pipeline; at R6 the K half and both fitter steps that read it are gone,
-    // so what remains above is the ψ half (until R5) and the structural inputs.
-    // The separation stands for the same reason it was made —
-    // but nothing outside needs the device, and serving it made /v1/bodies
+    // R8 (plan 07): what `model.planets.record()` PUBLISHES is the structural
+    // record, not the device. The internal record that kept the ψ/K law
+    // outputs for the fitting pipeline went piecewise: the K half at R6, the ψ
+    // half at R5, the lattice inputs at R1/R8 (nothing read it by then).
+    // Nothing outside needs the device, and serving it made /v1/bodies
     // contradict the model's own chain on the SAME quantities: Mercury's node
     // cycle 298,060 yr against the chain's 232,001, Jupiter's perihelion period
     // 68,783 against 304,456, Neptune's 670,634 against 2,136,796. The two
@@ -878,13 +837,15 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   };
   /** Ecliptic longitude of perihelion, J2000 ecliptic and equinox (the chain's ϖ). @param {string} k @param {number} year @returns {number} */
   const planetPerihelionDeg = (k, year) => planetChainElementsAt(k, year).lonPeriEclipticDeg;
-  /** Ascending node on the model's invariable plane, from the plane's ascending node on the ICRF equator (the Souami & Souchay 2012 origin; derived conversion, inv-plane-frame). @param {string} k @param {number} year @returns {number} */
-  const planetAscNodeDeg = (k, year) => {
+  /** The S&S longitude origin as an s-frame longitude (J2000-fixed, derived). @returns {number} */
+  const nodeOriginSSDeg = () => {
     if (planetNodeOriginSSDeg === null) {
       planetNodeOriginSSDeg = computeEquatorNodeOriginSFrameDeg(/** @type {any} */ (CHAIN_ARTIFACT).invariablePlane, C.earthOrbital.obliquityJ2000_deg);
     }
-    return convertNodeSFrameToEquatorOriginDeg(planetChainElementsAt(k, year).ascNodeInvPlaneDeg, planetNodeOriginSSDeg);
+    return planetNodeOriginSSDeg;
   };
+  /** Ascending node on the model's invariable plane, from the plane's ascending node on the ICRF equator (the Souami & Souchay 2012 origin; derived conversion, inv-plane-frame). @param {string} k @param {number} year @returns {number} */
+  const planetAscNodeDeg = (k, year) => convertNodeSFrameToEquatorOriginDeg(planetChainElementsAt(k, year).ascNodeInvPlaneDeg, nodeOriginSSDeg());
   /** Inclination to the model's invariable plane (exact orbit-normal rotation). @param {string} k @param {number} year @returns {number} */
   const planetInclinationDeg = (k, year) => planetChainElementsAt(k, year).inclInvPlaneDeg;
   /** Orbital eccentricity of date, from the same chain elements (R8): the
@@ -893,6 +854,52 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
    *  e-vector does not have. @param {string} k @param {number} year @returns {number} */
   const planetEccentricity = (k, year) => planetChainElementsAt(k, year).e;
 
+  // ── Earth's invariable-plane elements of date (plan 07 R5/R9) ─────────────
+  // ONE published i_inv / Ω_inv everywhere — the mirror of the browser's
+  // inclInvPlaneModel / ascNodeInvPlaneModel (identical ops): the engine's
+  // own Earth-orbit normal of date (the one-source sampler — the banked
+  // secular series inside ±10 Myr, the mode tail beyond) against the
+  // artifact's invariable plane, the node expressed in the Souami & Souchay
+  // longitude origin through the same derived conversion the planets' node
+  // takes. The H/3 cosine on the fixed apsidal carrier and the −H/5 linear
+  // node that served here until R5 were the K device: 1.578677° against
+  // 1.578422° at J2000, 1.24° against 1.91° at −100 kyr, and 112° against
+  // 94° of node there. Doc 05 banks the one-source reading at rms 0.003°
+  // vs La2010 over −500 kyr. Argument: the engine year (dynamical time,
+  // Julian from J2000 TT — the sampler's own coordinate, as the chain's).
+  const invPlaneGeometry = (() => {
+    const ip = /** @type {{ inclEclipticDeg: number, ascNodeEclipticDeg: number }} */ (/** @type {any} */ (CHAIN_ARTIFACT).invariablePlane);
+    const D2R = Math.PI / 180;
+    const si = Math.sin(ip.inclEclipticDeg * D2R), ci = Math.cos(ip.inclEclipticDeg * D2R);
+    const nInv = [si * Math.sin(ip.ascNodeEclipticDeg * D2R), -si * Math.cos(ip.ascNodeEclipticDeg * D2R), ci];
+    const xIp = [Math.cos(ip.ascNodeEclipticDeg * D2R), Math.sin(ip.ascNodeEclipticDeg * D2R), 0];
+    const yIp = [nInv[1] * xIp[2] - nInv[2] * xIp[1], nInv[2] * xIp[0] - nInv[0] * xIp[2], nInv[0] * xIp[1] - nInv[1] * xIp[0]];
+    return { ascNodeEclipticDeg: ip.ascNodeEclipticDeg, nInv, xIp, yIp };
+  })();
+  /** The engine Earth-orbit normal of date, J2000 ecliptic frame. @param {number} year @returns {number[]} */
+  const earthOrbitNormalJ2000At = (year) => {
+    const s = oneSourceM.sampleAt(year);
+    return [s.orbitNormalX, s.orbitNormalY, Math.sqrt(Math.max(0, 1 - s.orbitNormalX * s.orbitNormalX - s.orbitNormalY * s.orbitNormalY))];
+  };
+  /** Earth's inclination of date to the model's invariable plane, degrees. @param {number} year @returns {number} */
+  const inclinationDeg = (year) => {
+    const nE = earthOrbitNormalJ2000At(year), nInv = invPlaneGeometry.nInv;
+    const d = nE[0] * nInv[0] + nE[1] * nInv[1] + nE[2] * nInv[2];
+    return Math.acos(Math.min(1, Math.max(-1, d))) / (Math.PI / 180);
+  };
+  /** Earth's ascending node of date on the invariable plane, S&S origin, degrees in [0, 360). @param {number} year @returns {number} */
+  const ascendingNodeDeg = (year) => {
+    const D2R = Math.PI / 180;
+    const { ascNodeEclipticDeg, nInv, xIp, yIp } = invPlaneGeometry;
+    const nE = earthOrbitNormalJ2000At(year);
+    let N = [nInv[1] * nE[2] - nInv[2] * nE[1], nInv[2] * nE[0] - nInv[0] * nE[2], nInv[0] * nE[1] - nInv[1] * nE[0]];
+    const L = Math.hypot(N[0], N[1], N[2]);
+    N = [N[0] / L, N[1] / L, N[2] / L];
+    const sFrameDeg = ((Math.atan2(N[0] * yIp[0] + N[1] * yIp[1] + N[2] * yIp[2],
+      N[0] * xIp[0] + N[1] * xIp[1] + N[2] * xIp[2]) / D2R + ascNodeEclipticDeg) % 360 + 360) % 360;
+    return convertNodeSFrameToEquatorOriginDeg(sFrameDeg, nodeOriginSSDeg());
+  };
+
   // ── Time axis: exact JD ↔ model-year conversion ───────────────────────────
   // The model's `year` inputs live on the SI axis (the axis the fits were
   // anchored on — tools/lib `_jdToSIyear`): linear in SI 86400-s days from the