@essrt/physics 1.10.0 → 1.12.0

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@essrt/physics",
-  "version": "1.10.0",
+  "version": "1.12.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/eclipse/sun-planetary-completion.cjs

@@ -58,6 +58,17 @@
  * recomputes it from live constants and fails on mismatch. The api
  * centerline gate (≤8″ shadow-plane) backstops gross staleness
  * independently.
+ *
+ * E4/E5 PAIRING VERDICT (20.3h Phase 0): the framework-native Sun
+ * landing REPLACED the finder Sun form (Meeus Ch. 25 → the assembled
+ * framework Sun with the f(Y)+torque drift). The re-derivation this
+ * header calls for was MEASURED INSTEAD: the residual 2lE re-fit under
+ * the new chain is noise (−0.12″ sin / −0.03″ cos vs the declared-fitted
+ * +1.42″ — tools/explore/e3b-native-sun.mjs §2), the all-phase Sun
+ * residual improved to 0.95″, and every gate stayed green — the table's
+ * planetary content is chain-independent physics and the pairing holds
+ * WITHOUT re-extraction. The +1.42″ 2lE term is thereby measured to be
+ * genuinely planetary-class, not a Meeus-EoC absorber.
  */
 
 'use strict';

package/src/model.js

@@ -837,9 +837,18 @@ export function assembleModel(C, F) {
   //          osculating e = H/16 channel + inclination coupling.
   //   L(t) = L0 + the mean tropical rate + the f(Y) drift SHAPE only,
   //          ∫(rate_SI(y) − rate_SI(2000)) dy: the rate ANCHOR stays the mean
-  //          year (eclipse-endorsed); the drift is the Step 6d year-harmonic
+  //          year (eclipse-endorsed); the drift is the Step 6c year-harmonic
   //          claim in SI/TT (year-in-days × LOD — the LOD-day part of the
-  //          raw drift is UT-vs-TT and stays ΔT's job). Trapezoid table,
+  //          raw drift is UT-vs-TT and stays ΔT's job) PLUS the derived
+  //          TORQUE term (E5): the year harmonics carry only the GEOMETRIC
+  //          equinox displacement (the tilt nodes — the exact 8:3 amplitude
+  //          signature); the classical luni-solar torque adds a precession-
+  //          RATE modulation δp = −p₀·tan ε·δε(t) on the model's own
+  //          two-component obliquity law (−A cos φ₃ + A cos φ₈). Both
+  //          lengthen the year at obliquity max, so they ADD; per-divisor
+  //          drift scale 1 + p₀·tan²ε·H/(2π·div) = 1.306 (H/8) / 1.815
+  //          (H/3) — the structure the ancient corpus blind-selected before
+  //          the derivation existed. Zero new constants. Trapezoid table,
   //          10-yr steps over −3000..3000; outside, the drift freezes at the
   //          edge (rate reverts to linear — the finder domain is the corpus
   //          era).
@@ -858,31 +867,64 @@ export function assembleModel(C, F) {
     return eccentricityAt(year) + (eccentricityBase / 2) * (cos3 - cos3J2000);
   };
   const sunMeanLongitudeDegAt = (() => {
-    const Y_LO = -3000, Y_HI = 3000, STEP = 10;
-    /** @type {Float64Array|null} */
-    let cum = null;
+    // SW PHASE B — the CLOSED FORM on the integrated lattice phase
+    // (supersedes the E5 ±3,000-yr trapezoid table; valid at EVERY epoch,
+    // no table, no domain window). drift = D_smooth + Σₖ D_k + D_torque:
+    //  · D_smooth — the deep chain's smooth SI tropical-year physics
+    //    (tropicalYearDaysBase × LOD), two-point trapezoid: the same
+    //    structure the scene twin uses, so the scene δ stays bounded.
+    //  · Each oscillatory term integrates ANALYTICALLY: the antiderivative
+    //    of cos/sin on the integrated phase is (H/2πk)·[sin/−cos] — the H
+    //    drift lives in the phase itself; the amplitude factor uses H_J2000
+    //    (ppm-class difference, negligible).
+    //  · THE REBASE (E5's corpus-preferred J2000 rate anchor) is applied
+    //    PER HARMONIC as a SIN-SATURATED ramp, (H/2πk)·sin(Δφₖ): equal to
+    //    the linear rebase for |Δφₖ| ≪ 1 (the corpus era) and bounded by
+    //    the harmonic's own period beyond — extending a linearized local
+    //    slope to deep time would be the fitted-linear-slope trap. Every
+    //    quantity is derived; zero new constants; each component and its
+    //    slope vanish at year 2000 by construction (the drift-only
+    //    property is exact, not tabulated).
+    const RATE_LIN = sunTropicalRateDegPerCy / 100;              // deg / SI yr
+    const precessionP0DegPerYr = 13 * 360 / H;
+    const tanEps = Math.tan(earthtiltMean * Math.PI / 180);
+    const A_RAD = earthInclAmplitude * Math.PI / 180;
     /** @param {number} year @returns {number} */
-    const rateSIYr = (year) => 360 * (365.25 * 86400)
-      / (tropicalYearDirectDays(year)
+    const rateSmooth = (year) => 360 * (365.25 * 86400)
+      / (tropicalYearDaysBase(year)
         * (deepLod.lodSecondsAtAge(yearToTMa(year)) ?? meanLengthOfDay));
+    /** @type {{rateSm0: number, ph0: Map<number, number>}|null} */
+    let anchor = null;
     return /** @param {number} year @returns {number} */ (year) => {
-      if (cum === null) {
-        const n = (Y_HI - Y_LO) / STEP + 1;
-        cum = new Float64Array(n);
-        const ref = rateSIYr(2000);
-        let prev = rateSIYr(Y_LO) - ref;
-        for (let i = 1; i < n; i++) {
-          const next = rateSIYr(Y_LO + i * STEP) - ref;
-          cum[i] = cum[i - 1] + 0.5 * (prev + next) * STEP;
-          prev = next;
-        }
-        const at2000 = cum[(2000 - Y_LO) / STEP];
-        for (let i = 0; i < n; i++) cum[i] -= at2000;
+      if (anchor === null) {
+        const ph0 = new Map();
+        for (const [k] of F.TROPICAL_YEAR_HARMONICS) ph0.set(k, phaseRadians(balancedYear, 2000, k));
+        for (const k of [3, 8]) if (!ph0.has(k)) ph0.set(k, phaseRadians(balancedYear, 2000, k));
+        anchor = { rateSm0: rateSmooth(2000), ph0 };
       }
-      const x = (Math.min(Math.max(year, Y_LO), Y_HI) - Y_LO) / STEP;
-      const i = Math.min(Math.floor(x), cum.length - 2);
-      const drift = cum[i] + (x - i) * (cum[i + 1] - cum[i]);
-      return sunL0Deg + sunTropicalRateDegPerCy * (year - 2000) / 100 + drift;
+      const dy = year - 2000;
+      // smooth deep physics (two-point trapezoid; exact 0 at 2000)
+      let drift = 0.5 * (rateSmooth(year) - anchor.rateSm0) * dy;
+      // year-harmonic ripple: rate ≈ −(RATE_LIN/T̄)·h(y); antiderivative +
+      // sin-saturated per-harmonic rebase
+      const HK = H / (2 * Math.PI);
+      for (const [k, sK, cK] of F.TROPICAL_YEAR_HARMONICS) {
+        const ph = phaseRadians(balancedYear, year, k);
+        const ph0 = /** @type {number} */ (anchor.ph0.get(k));
+        const scale = -(RATE_LIN / meanSolarYearDays) * (HK / k);
+        const anti = (-sK) * (Math.cos(ph) - Math.cos(ph0)) + cK * (Math.sin(ph) - Math.sin(ph0));
+        const h0 = sK * Math.sin(ph0) + cK * Math.cos(ph0);
+        drift += scale * (anti - h0 * Math.sin(ph - ph0));
+      }
+      // torque (E5): rate = −p₀·tanε·δε, δε = A(−cos φ₃ + cos φ₈); same
+      // antiderivative + per-component sin-saturated rebase
+      for (const [k, sgn] of [[3, -1], [8, 1]]) {
+        const ph = phaseRadians(balancedYear, year, k);
+        const ph0 = /** @type {number} */ (anchor.ph0.get(k));
+        const scale = -precessionP0DegPerYr * tanEps * A_RAD * sgn * (HK / k);
+        drift += scale * ((Math.sin(ph) - Math.sin(ph0)) - Math.cos(ph0) * Math.sin(ph - ph0));
+      }
+      return sunL0Deg + sunTropicalRateDegPerCy * dy / 100 + drift;
     };
   })();
 

package/src/moon/series-extension.cjs

@@ -62,7 +62,14 @@ const DELAUNAY = {
 };
 
 /** Longitude tail: [kD, kM, kMp, kF, sin″] — 47 derived terms ≥ 0.1″
- *  (RMS content 2.02″). */
+ *  (RMS content 2.02″), plus the 20.3h ROUND-2 block at the end: five
+ *  dt-halving-converged terms beyond the D2 catalog's order-7 bound
+ *  (m20h-extraction-r2.mjs, the same D1 3-body lab). The lead term
+ *  [6,0,−2,0] extracted at 0.572″ vs the MPP02-comparison target 0.57″
+ *  (m20h-alias-breaker.mjs) — the integrator predicted the independent
+ *  reference again. 4-gate out-of-sample (m20h-gate-preview.mjs):
+ *  all-phase λ 3.03 → 2.96″, β unchanged, fleet improves, centerlines
+ *  2.71 → 2.69″ with the 2001 max 7.3 → 6.7″. */
 const LON_TERMS = [
   [2, 0, 1, 2, -0.9900],
   [2, 0, -4, 0, 0.9465],
@@ -111,6 +118,12 @@ const LON_TERMS = [
   [1, 0, -3, 0, -0.1252],
   [2, 0, 2, 2, -0.1236],
   [1, -1, 1, 0, -0.1225],
+  // 20.3h ROUND 2 — beyond the order-7 catalog bound (see header note):
+  [6, 0, -2, 0, 0.572],
+  [6, 0, -3, 0, 0.294],
+  [4, 0, -2, 2, -0.169],
+  [4, -2, -2, 0, 0.161],
+  [0, 0, 5, 0, 0.112],
 ];
 
 /** Latitude tail: [kD, kM, kMp, kF, sin″] — 30 derived terms ≥ 0.1″