@essrt/physics 1.3.0 → 1.3.1

package/README.md

@@ -53,7 +53,7 @@ The canonical assembly, wired once inside the package:
 import { createModel, DEFAULT_CONSTANTS } from '@essrt/physics';
 
 const model = createModel();
-model.identity.modelVersion;          // 'v10.0' + both content hashes
+model.identity.modelVersion;          // the shipped model identity (e.g. 'v11.0') + both content hashes
 model.epoch.hAtYear(2000 - 380e6);    // 306189 — Devonian H (the Wells 1963 match)
 model.earth.obliquityDeg(2000);       // 23.4393
 model.cardinal.jd(2000, 'SS');        // 2451716.575 — June solstice 2000 (USNO)

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@essrt/physics",
-  "version": "1.3.0",
+  "version": "1.3.1",
   "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/deltat/deep-time.cjs

@@ -113,7 +113,15 @@ function createDeepTimeLod(deps) {
   function siderealYearSecondsAtAge(t_Ma) {
     if (t_Ma === 0) return K.meanSiderealYearJ2000Seconds;
     const mass_loss_fraction = K.solarMassLossFracPerYear * t_Ma * 1e6;
-    return K.meanSiderealYearJ2000Seconds * (1 - 2 * mass_loss_fraction);
+    // Exact Kepler under the model's mass-history convention (a·M = const,
+    // a(t) = a₀·(1−Δm)): T ∝ M⁻² ⇒ T(t) = T₀·(1−Δm)² — the SAME Driver-2
+    // law the planet chains use (orbit-chain driver2PeriodSecondsAtAge).
+    // The previous (1 − 2Δm) was this law's first-order Taylor — the one
+    // form in the model inconsistent with it, and the entire source of the
+    // day-count invariant's 0.17 ppm residual: with the product form,
+    // H·(days/yr)·(AU₀/AU)² = TOTAL_DAYS_IN_H holds EXACTLY
+    // (tools/explore/deep-time-sensitivity.js §5; docs/99 §near-invariant).
+    return K.meanSiderealYearJ2000Seconds * (1 - mass_loss_fraction) * (1 - mass_loss_fraction);
   }
 
   /** @param {number} t_Ma @returns {number} tropical = sidereal · (1 − 13/H(t)) */

package/src/layer0/index.js

@@ -118,7 +118,12 @@ export const createEpochPrimitives = ({ params: p, alphaAtAgeMa }) => {
    */
   const siderealYearSecondsCore = (t) => {
     if (t === 0) return p.siderealYearJ2000Seconds;
-    return p.siderealYearJ2000Seconds * (1 - 2 * p.solarMassLossFracPerYear * t * 1e6);
+    // Exact Kepler product form (T ∝ M⁻² under a·M = const, a(t) = a₀·(1−Δm))
+    // — the same Driver-2 law as the planet chains and deltat/deep-time.cjs.
+    // MUST stay the bit-exact twin of siderealYearSecondsAtAge there (the
+    // transparency gate holds the pair at 84/84 round-trip).
+    const dm = p.solarMassLossFracPerYear * t * 1e6;
+    return p.siderealYearJ2000Seconds * (1 - dm) * (1 - dm);
   };
 
   /**