@essrt/physics 1.10.0 → 1.11.0

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@essrt/physics",
-  "version": "1.10.0",
+  "version": "1.11.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).
@@ -861,10 +870,19 @@ export function assembleModel(C, F) {
     const Y_LO = -3000, Y_HI = 3000, STEP = 10;
     /** @type {Float64Array|null} */
     let cum = null;
+    const precessionP0DegPerYr = 13 * 360 / H;
+    const tanEps = Math.tan(earthtiltMean * Math.PI / 180);
+    /** E5 — the two-component obliquity excursion, radians, integrated phase.
+     * @param {number} year @returns {number} */
+    const obliquityExcursionRad = (year) => (
+      -earthInclAmplitude * Math.cos(phaseRadians(balancedYear, year, 3))
+      + earthInclAmplitude * Math.cos(phaseRadians(balancedYear, year, 8))
+    ) * Math.PI / 180;
     /** @param {number} year @returns {number} */
     const rateSIYr = (year) => 360 * (365.25 * 86400)
       / (tropicalYearDirectDays(year)
-        * (deepLod.lodSecondsAtAge(yearToTMa(year)) ?? meanLengthOfDay));
+        * (deepLod.lodSecondsAtAge(yearToTMa(year)) ?? meanLengthOfDay))
+      - precessionP0DegPerYr * tanEps * obliquityExcursionRad(year);
     return /** @param {number} year @returns {number} */ (year) => {
       if (cum === null) {
         const n = (Y_HI - Y_LO) / STEP + 1;

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″