@essrt/physics 1.13.0 → 1.15.0

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@essrt/physics",
-  "version": "1.13.0",
+  "version": "1.15.0",
   "description": "ESSRT (Expanding Solar System Resonance Theory) physics core \u2014 pure, injectable-constants, no I/O, no globals, no DOM. The complete model incl. fitted coefficients; every published version immutably ships one recorded model identity.",
   "essrt": {
     "modelVersion": "v11.0",

package/src/constants/generated.js

@@ -766,8 +766,8 @@ export const REFERENCE_DATA = Object.freeze({
     "sunSiderealYearDiffSeconds": 0.02,
     "frameworkSunVsJplRmsArcsec": 0.95,
     "meeusCh25SunVsJplRmsArcsec": 1.28,
-    "moonSeriesLonVsJplRmsArcsec": 2.96,
-    "moonSeriesLatVsJplRmsArcsec": 0.65,
+    "moonSeriesLonVsJplRmsArcsec": 2.84,
+    "moonSeriesLatVsJplRmsArcsec": 0.35,
     "mercuryNewtonianArcsecCy": 532,
     "mercuryObservedICRFArcsecCy": 575,
     "mercuryPark2017RateArcsecCy": 575.31,

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

@@ -235,4 +235,16 @@ function createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy })
  *  landings. */
 const PAIRED_SUN_HARMONICS_SHA256 = 'e2cf42e9770c9e0a';
 
-module.exports = { createSunPlanetaryCompletion, PAIRED_SUN_HARMONICS_SHA256 };
+/** sha256/16 of JSON.stringify([...planets, moonElongation]) — the seven
+ *  full-precision carrier rates (deg/cy TT) the TERMS table was extracted
+ *  under at N3. The model wiring recomputes the rates live from the planet
+ *  records, so a planet-period / year / month input change moves the
+ *  carriers automatically while the table stays frozen — a silent few-
+ *  arcsec stale below the api gate's ≤8″ backstop. test:model recomputes
+ *  this fingerprint from live constants (identical arithmetic to model.js)
+ *  and fails on mismatch: re-run the N3 extraction chain
+ *  (tools/explore/n2-sun-framework-carriers.mjs →
+ *  n3-carrier-swap-preview.mjs), re-embed TERMS, and update this value. */
+const PAIRED_CARRIER_RATES_SHA256 = '893e055ee12343bd';
+
+module.exports = { createSunPlanetaryCompletion, PAIRED_SUN_HARMONICS_SHA256, PAIRED_CARRIER_RATES_SHA256 };

package/src/moon/series-extension.cjs

@@ -56,6 +56,53 @@
  * (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,
+ * 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
+ * lunar orbit — the same channel whose secular piece DLT-1 derives as
+ * the Earth-figure +0.1925″/cy² T² term). Extracted from the B−A
+ * differential of epoch-phased 3-body twin integrations (J2 on vs off,
+ * identical ICs; per-run fits on each run's OWN mean angles — fitting
+ * the raw differential leaks A·δrate·t quadrature, measured as 2.4″
+ * fake rows), dt-halving converged ≤0.003″. TWO DOCUMENTED CONTROLS
+ * certify the channel and calibrate the run→instrument frame mapping:
+ * the derived pure-Ω λ term reproduces the Meeus flattening additive
+ * (+1962e-6 sin(Lp−F)) at 1.028, and the derived β sin(Lp) reproduces
+ * the −2235e-6 additive at 0.971 — both already shipped via the
+ * catalog additive family and therefore NOT re-shipped here. JPL
+ * out-of-sample (validation only, never fit): official all-phase λ
+ * 2.964 → 2.922″, β 0.653 → 0.606″; dense 2-day λ 2.378 → 2.321″,
+ * β 0.655 → 0.601″; syzygy-fleet elongation improves. Instruments:
+ * fq7-band-probe.mjs · fq7-jpl-band.mjs · fq7-j2-node.mjs ·
+ * fq7-r3-ship.mjs · fq7-r3-preview.mjs. The Ω carrier polynomial is an
+ * IAU-class instrument literal — the same DECLARED-INPUT status as the
+ * Delaunay carriers above (Ω ≡ Lp − F of the same catalog class).
+ *
+ * THE FQ-7 DUST ROUND (plan §12i) — 43 λ + 33 β derived Delaunay terms
+ * appended to LON_TERMS/LAT_TERMS below. The broadband-dust census
+ * (fq7-dust-census.mjs, secular-guarded v2: cubic detrend + T/T²
+ * quadrature absorbers on the head arguments, without which the
+ * doctrine-blocked T-modulated parameter class masquerades as fake
+ * 184-yr sideband families) ranked the post-R3 capturable content; the
+ * dense-JPL arbitration (fq7-dust-arbitration.mjs) measured EVERY
+ * trusted leader REAL at ratio 0.9–1.1; the labels bit-match MPP02
+ * main-problem catalog rows below the Meeus ~0.4″ truncation cutoff, in
+ * argument classes no earlier extraction catalog enumerated (odd kD,
+ * |kM| = 3, |kF| = 4/5). Main-problem ⇒ the D1 3-body lab derives them:
+ * one 120-yr run, λ and β columns, head fidelity 0.51%/0.07%,
+ * dt-halving drift 0.000″, lab/census ratio ∈ [0.7, 1.4] required
+ * (4 λ census rows the lab does NOT fully produce stay unshipped —
+ * derive or don't ship). Lab signs match the MPP02 catalog signs row
+ * for row (the frame-consistency control). JPL out-of-sample: official
+ * all-phase λ 2.922 → 2.840″, β 0.606 → 0.353″ — the β endpoint AT the
+ * R1-measured MPP02-vs-JPL floor (~0.33″); dense 2-day λ 2.321 →
+ * 2.219″, β 0.601 → 0.384″; syzygy fleet improves; the sign-flip
+ * control degrades every gate. Instruments: fq7-dust-census.mjs ·
+ * fq7-dust-arbitration.mjs · fq7-dust-extraction.mjs ·
+ * fq7-dust-preview.mjs.
  */
 
 'use strict';
@@ -131,6 +178,52 @@ const LON_TERMS = [
   [4, 0, -2, 2, -0.169],
   [4, -2, -2, 0, 0.161],
   [0, 0, 5, 0, 0.112],
+  // FQ-7 DUST ROUND — the argument classes no earlier catalog enumerated
+  // (odd kD, |kM| = 3, |kF| = 4); lab-derived, MPP02-catalog-certified,
+  // JPL-arbitrated (see header note):
+  [0, 0, 0, 4, 0.422],
+  [3, 0, 0, 0, 0.402],
+  [2, -3, 0, 0, 0.341],
+  [3, 1, -1, 0, 0.275],
+  [3, 0, 0, -2, -0.250],
+  [2, -3, -1, 0, 0.240],
+  [3, -1, -1, 0, -0.233],
+  [3, 1, 0, 0, 0.147],
+  [0, 3, 0, 0, -0.096],
+  [4, -1, -3, 0, 0.096],
+  [2, -1, 3, 0, 0.094],
+  [2, 0, -3, -2, -0.093],
+  [0, 0, 1, 4, 0.092],
+  [6, -1, -2, 0, 0.091],
+  [1, 1, 2, 0, 0.090],
+  [4, 0, 0, 2, -0.084],
+  [3, -1, -2, 0, -0.082],
+  [0, 0, 1, -4, -0.081],
+  [2, 0, 4, 0, 0.078],
+  [2, 0, 0, -4, -0.078],
+  [0, 1, 0, -2, 0.074],
+  [2, -1, 1, 2, -0.073],
+  [2, -1, -1, -2, -0.073],
+  [6, -1, -1, 0, 0.072],
+  [2, 1, -1, 2, 0.068],
+  [2, 0, -5, 0, 0.067],
+  [2, 1, 0, 2, 0.066],
+  [4, 0, 1, -2, -0.065],
+  [2, -2, 2, 0, 0.065],
+  [2, -2, 0, -2, 0.065],
+  [3, -1, 0, 0, 0.063],
+  [0, 2, 2, 0, -0.062],
+  [5, 0, -2, 0, -0.059],
+  [0, 1, -2, -2, 0.050],
+  [4, 1, 1, 0, -0.050],
+  [0, 1, -4, 0, -0.048],
+  [0, 1, 2, 2, 0.045],
+  [1, 0, 1, 2, 0.044],
+  [3, 0, -3, 0, -0.043],
+  [1, 0, 3, 0, -0.042],
+  [1, 2, 0, 0, -0.041],
+  [0, 3, -1, 0, -0.038],
+  [1, -2, 0, 0, 0.037],
 ];
 
 /** Latitude tail: [kD, kM, kMp, kF, sin″] — 30 derived terms ≥ 0.1″
@@ -166,6 +259,41 @@ const LAT_TERMS = [
   [1, 0, -2, -1, -0.1093],
   [2, 2, -1, 1, -0.1044],
   [1, 1, 1, 1, 0.1014],
+  // FQ-7 DUST ROUND — lab-derived, MPP02-catalog-certified,
+  // JPL-arbitrated (see header note):
+  [3, 0, 0, -1, -0.352],
+  [3, 0, -1, -1, -0.305],
+  [4, -1, -2, 1, 0.214],
+  [3, 0, -1, 1, -0.206],
+  [0, 2, -1, 1, -0.094],
+  [6, 0, -1, -1, 0.094],
+  [0, 0, 4, -1, 0.091],
+  [2, -1, 0, -3, 0.087],
+  [6, 0, -2, -1, 0.081],
+  [1, 0, -2, 1, -0.078],
+  [0, 1, -3, -1, -0.075],
+  [2, 1, -2, -1, 0.073],
+  [2, -2, 1, -1, 0.069],
+  [2, 0, -2, 3, -0.065],
+  [1, 0, 2, 1, -0.064],
+  [2, -1, -2, 1, -0.064],
+  [2, 1, 2, -1, -0.063],
+  [4, 0, 0, -3, 0.062],
+  [0, 1, 3, 1, -0.062],
+  [6, 0, -2, 1, 0.060],
+  [2, -2, 1, 1, 0.059],
+  [2, 1, 0, -3, -0.059],
+  [4, -1, 1, -1, 0.057],
+  [0, 2, 0, 1, -0.055],
+  [0, 2, 1, 1, -0.054],
+  [2, -1, -3, -1, 0.054],
+  [4, -1, -2, -1, 0.052],
+  [4, 1, -1, 1, -0.051],
+  [1, 1, -1, 1, -0.051],
+  [2, 1, -2, 1, 0.049],
+  [3, 0, -2, 1, -0.048],
+  [4, 0, 2, -1, 0.048],
+  [6, 0, -1, 1, 0.042],
 ];
 
 /** Planetary mean-longitude polynomials, degrees per Julian century TT
@@ -198,6 +326,27 @@ const PLANETARY_LON_TERMS = [
   [[0, 1, 0, -1, 0, 0, -1], -0.0077, 0.2235],
 ];
 
+/** Lunar mean-node polynomial [deg at J2000, deg/cy, deg/cy²] TT —
+ *  the Ω carrier for the J2 node family (declared input, see header). */
+const OMEGA = [125.0445479, -1934.1362891, 0.0020754];
+
+/**
+ * The J2 node family (see header): [[kD,kM,kMp,kF,kΩ], sin″] — pure
+ * sin-parity (the derivation's cos components read ≤0.005″), ADD-composed
+ * like the Delaunay tail.
+ * @type {Array<[number[], number]>}
+ */
+const NODE_LON_TERMS = [
+  [[0, 0, 1, 0, -1], -0.544],
+  [[0, 0, 1, 0, 1], 0.546],
+  [[0, 0, 0, 2, 1], 0.371],
+  [[2, 0, 0, 0, 1], 0.103],
+];
+/** @type {Array<[number[], number]>} */
+const NODE_LAT_TERMS = [
+  [[0, 0, 0, 1, -1], -0.375],
+];
+
 const D2R = Math.PI / 180;
 
 /**
@@ -219,6 +368,9 @@ function moonSeriesExtensionDeg(T) {
     const th = k[0] * l[0] + k[1] * l[1] + k[2] * l[2] + k[3] * l[3] + k[4] * l[4] + k[5] * D + k[6] * Mp;
     lon += cA * Math.cos(th) + sA * Math.sin(th);
   }
+  const Om = poly(OMEGA);
+  for (const [k, s] of NODE_LON_TERMS) lon += s * Math.sin(k[0] * D + k[1] * M + k[2] * Mp + k[3] * F + k[4] * Om);
+  for (const [k, s] of NODE_LAT_TERMS) lat += s * Math.sin(k[0] * D + k[1] * M + k[2] * Mp + k[3] * F + k[4] * Om);
   return { dLonDeg: lon / 3600, dLatDeg: lat / 3600 };
 }