@essrt/physics 1.15.0 → 1.17.0

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@essrt/physics",
-  "version": "1.15.0",
+  "version": "1.17.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

@@ -764,8 +764,8 @@ export const REFERENCE_DATA = Object.freeze({
     "sunModelTrueErrorDeg": 0.003,
     "sunTropicalYearDiffSeconds": 0.1,
     "sunSiderealYearDiffSeconds": 0.02,
-    "frameworkSunVsJplRmsArcsec": 0.95,
-    "meeusCh25SunVsJplRmsArcsec": 1.28,
+    "frameworkSunVsJplRmsArcsec": 0.8,
+    "meeusCh25SunVsJplRmsArcsec": 1.22,
     "moonSeriesLonVsJplRmsArcsec": 2.84,
     "moonSeriesLatVsJplRmsArcsec": 0.35,
     "mercuryNewtonianArcsecCy": 532,

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

@@ -41,10 +41,21 @@
  *  (ii) the J2000 phase anchors ARG_L0 / PERI / D0 — epoch initial
  *       conditions ("anchored by design" class; the RATES are now
  *       framework-derived, injected);
- *  (iii) the 2lE term (+1.42″ sin 2lE) — DECLARED FITTED, not derived: a
- *       finder-annual artifact, measured chain-independent at E4/E5
- *       (re-fit noise −0.12″/−0.03″). Everything else in this file is
- *       derived. The Earth-around-EMB wobble is DERIVED: amplitude
+ *  (iii) RETIRED (FQ-7-Sun annual-channel attribution): the former
+ *       declared-fitted +1.42″ sin 2lE term. It was fitted against a
+ *       Sun-ALONE JPL comparison bridged by the leading nutation term
+ *       only, and the attribution measured it as the semiannual
+ *       nutation term −1.32″ sin(2F−2D+2Ω) ≡ −1.32″ sin 2lE in disguise
+ *       (2M residual 0.13″ under that bridge, 1.43″ ≡ the term's
+ *       negative under the fuller bridge). Nutation in longitude is a
+ *       frame rotation common to Sun AND Moon — this chain keeps both
+ *       bodies MEAN-of-date (mean obliquity, mean sidereal time; see
+ *       besselian.cjs), so a Sun-only nutation term is a frame
+ *       inconsistency for elongation. Removed: the nutation-free syzygy
+ *       fleet improves 3.877 → 3.756″ and the fuller-bridged Sun
+ *       comparison 2.344 → 2.123″ (fq7s-jpl-preview.mjs --no-2le;
+ *       fq7s-annual-channel.mjs). This file now 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
  *       negative in the subtract convention.
@@ -175,9 +186,6 @@ const TERMS = [
   [[0, 0, 1, -1, 0, 0], [0, 0, -1, 0, 0, 0], -0.0353, -0.0379],
 ];
 
-/** The declared-fitted semiannual leftover (header (iii)): +1.42″ sin 2lE. */
-const FITTED_2LE_ARCSEC = 1.42;
-
 const D2R = Math.PI / 180;
 
 /**
@@ -222,8 +230,7 @@ function createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy })
       table += cA * Math.cos(th) + sA * Math.sin(th);
     }
     const D = (ARG_D0 + D1 * T) * D2R;
-    const arcsec = -table - embWobbleArcsec * Math.sin(D)
-      + FITTED_2LE_ARCSEC * Math.sin(2 * l[2]);
+    const arcsec = -table - embWobbleArcsec * Math.sin(D);
     return arcsec / 3600;
   }
   return { sunPlanetaryCompletionDeg };

package/src/model.js

@@ -860,11 +860,28 @@ export function assembleModel(C, F) {
   // D2 completion table is unchanged (residual 2lE re-fit ≈ 0.1″ — noise).
   const sunL0Deg = C.earthOrbital.sunMeanLongitudeJ2000_deg;
   const sunTropicalRateDegPerCy = 360 * julianCenturyDays / meanSolarYearDays;
+  // FQ-7-Sun option C-small — ONE eccentricity law for Sun and Moon in the
+  // eclipse chain. Eccentricity is frame-invariant (e = |z|, z = e·e^{iϖ}),
+  // so it may carry only FIXED-FRAME lattice content: the H/3 line the Moon's
+  // E-factor already rides (moon/ecc-channel.cjs). H/16 is the OF-DATE
+  // perihelion period (13 + 3 = 16: the H/3 apsidal rotation seen from the
+  // H/13 equinox) — it belongs to ϖ_of-date, not to e. The previous form
+  // added the H/16 law's own slope (−0.84e-5/cy) on top of the H/3 line,
+  // and that slope WAS the measured Sun-side ė tension (annual channel
+  // T·sinM −3.7″/cy; JPL sides with the H/3 line: fq7s-ecc-consistency.mjs).
+  // Anchor: the H/16 law's J2000 value (the observed e, anchored by design),
+  // evaluated once — an explicit law change, not a hidden freeze; the
+  // cardinal-point path keeps the H/16 law (eccentricityAt) unchanged.
+  // Measured at the swap: JPL Sun 2.123 → 1.487″ (1900–2100), registry
+  // window 0.93 → 0.80″, syzygy fleet 3.756 → 3.719″, T·sinM −3.65 → −0.19″/cy,
+  // curvature term unchanged in size; the H/16 term at the corpus epochs
+  // detrends to 0.174 min — below the ancient corpus's discrimination.
+  const sunEccentricityJ2000 = eccentricityAt(2000);
   /** @param {number} year @returns {number} */
   const sunEccentricityAt = (year) => {
     const cos3 = -Math.cos(phaseRadians(balancedYear, year, 3));
     const cos3J2000 = -Math.cos(phaseRadians(balancedYear, 2000, 3));
-    return eccentricityAt(year) + (eccentricityBase / 2) * (cos3 - cos3J2000);
+    return sunEccentricityJ2000 + (eccentricityBase / 2) * (cos3 - cos3J2000);
   };
   const sunMeanLongitudeDegAt = (() => {
     // SW PHASE B — the CLOSED FORM on the integrated lattice phase