@essrt/physics 9.2.0 → 11.0.0

package/types/earth/deep-orbital-history.d.cts

@@ -13,6 +13,8 @@
  *   axialPrecessionYearsJ2000: number,
  *   obliquityJ2000Deg: number,
  *   axialPrecessionYearsAtYearFn?: (year: number) => number,
+ *   solarTorqueShareJ2000?: number,
+ *   dynamicalEllipticityRatioAtYearFn?: (year: number) => number,
  * }} deps — mode tables from the governed deep artifact (choose the ζ tier
  *   per consumer: era = its own 8-term extraction, deep = the 16-term
  *   table; NEVER a slice); anchors from the chain artifact's one home; the
@@ -25,9 +27,12 @@
  *   vs La2004 over −200 kyr vs the deep tier's 0.069° — the C-1 verdict,
  *   plan 02) and zetaModes serves only as the TAIL beyond the span (the
  *   seam at the span edge is the extraction residual, far outside every
- *   certified window).
+ *   certified window). solarTorqueShareJ2000 and
+ *   dynamicalEllipticityRatioAtYearFn: THE PRECESSION OF DATE — the two
+ *   factors on the precession constant (see alphaAtGeneral below); absent,
+ *   the constant is the secular law alone, bit-identical.
  */
-export function createDeepOrbitalHistory({ zModes, zetaModes, zetaSeries, zSeries, anchorE, anchorPeriEclipticDeg, anchorInclEclipticDeg, anchorAscNodeEclipticDeg, axialPrecessionYearsJ2000, obliquityJ2000Deg, axialPrecessionYearsAtYearFn, }: {
+export function createDeepOrbitalHistory({ zModes, zetaModes, zetaSeries, zSeries, anchorE, anchorPeriEclipticDeg, anchorInclEclipticDeg, anchorAscNodeEclipticDeg, axialPrecessionYearsJ2000, obliquityJ2000Deg, axialPrecessionYearsAtYearFn, solarTorqueShareJ2000, dynamicalEllipticityRatioAtYearFn, }: {
     zModes: ReadonlyArray<{
         omegaRadPerYr: number;
         re: number;
@@ -57,6 +62,8 @@ export function createDeepOrbitalHistory({ zModes, zetaModes, zetaSeries, zSerie
     axialPrecessionYearsJ2000: number;
     obliquityJ2000Deg: number;
     axialPrecessionYearsAtYearFn?: (year: number) => number;
+    solarTorqueShareJ2000?: number;
+    dynamicalEllipticityRatioAtYearFn?: (year: number) => number;
 }): {
     build: (t0Yr: number, t1Yr: number, stepYr: number) => {
         /** @param {number} tYr */
@@ -71,9 +78,20 @@ export function createDeepOrbitalHistory({ zModes, zetaModes, zetaSeries, zSerie
             orbitNormalY: number;
             equinoxLonJ2000Deg: number;
             equinoxLonRateDegPerYr: number;
+            generalPrecessionLonDeg: number;
+            generalPrecessionLonRateDegPerYr: number;
         };
     };
     alphaArcsecPerYr: number;
     alphaLunisolarArcsecPerYr: number;
     axialPrecessionYearsLunisolarJ2000: number;
+    /** THE PRECESSION OF DATE — the two factors this integration applies to
+     *  the secular precession constant at t years from J2000 (both 1 at
+     *  t = 0; both 1 everywhere when the options are absent). ONE home for
+     *  the published surface and the docs.
+     *  @param {number} t @returns {{solarTorque: number, ellipticity: number}} */
+    precessionOfDateFactorsAt: (t: number) => {
+        solarTorque: number;
+        ellipticity: number;
+    };
 };

package/types/earth/sidereal-channel-artifact.d.cts

@@ -1,7 +1,7 @@
 export const SIDEREAL_CHANNEL_ARTIFACT: Readonly<{
     t0Yr: -10000000;
     stepYr: 2000;
-    windowYr: 2000;
+    windowYr: 10000;
     lamDotRel: number[];
     meta: {
         dumpSha256: string;
@@ -9,4 +9,4 @@ export const SIDEREAL_CHANNEL_ARTIFACT: Readonly<{
         chaprontGate: number;
     };
 }>;
-export const SIDEREAL_CHANNEL_ARTIFACT_HASH: "528101ca3410dc9b";
+export const SIDEREAL_CHANNEL_ARTIFACT_HASH: "a6f4729c2de722af";

package/types/earth/year-lengths.d.cts

@@ -1,29 +1,23 @@
 /**
  * @param {{
- *   sampleAt: (year: number) => {e: number, periOfDateDeg: number, equinoxLonJ2000Deg: number},
+ *   sampleAt: (year: number) => {e: number, periOfDateDeg: number, equinoxLonJ2000Deg: number, generalPrecessionLonDeg: number},
  *   massLossSiderealSecondsAtYearFn: (year: number) => number,
- *   apsidalSecularDegPerYrFn?: (year: number) => number,
  * }} opts
- *   sampleAt: the one-source movement sampler, calendar-year keyed
- *   (createDeepOrbitalHistory build().at wrapped by the caller — series
- *   tier where the artifact is available, mode tail beyond).
+ *   sampleAt: the one-source movement sampler, keyed in the sampler's own
+ *   year (Julian years from J2000 TT, 2000 = J2000; createDeepOrbitalHistory
+ *   build().at wrapped by the caller — series tier where the artifact is
+ *   available, mode tail beyond).
  *   massLossSiderealSecondsAtYearFn: the caller's secular mass-loss
  *   sidereal-year law, SI seconds, IAU-anchored at J2000.
- *   apsidalSecularDegPerYrFn: the chain's SECULAR apsidal tangent of date
- *   (keplerian-chain computeApsidalSecularDegPerYr for Earth — deg/yr,
- *   J2000-ecliptic frame). When given, the anomalistic year is the
- *   mean-element construction on this smooth rate; when absent, it falls
- *   back to the cardinal structure's measured year-over-year interval
- *   (which carries the banked series' century-scale ϖ̇ wobble).
  */
-export function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn, apsidalSecularDegPerYrFn }: {
+export function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn }: {
     sampleAt: (year: number) => {
         e: number;
         periOfDateDeg: number;
         equinoxLonJ2000Deg: number;
+        generalPrecessionLonDeg: number;
     };
     massLossSiderealSecondsAtYearFn: (year: number) => number;
-    apsidalSecularDegPerYrFn?: (year: number) => number;
 }): Readonly<{
     /** Mean tropical year of date, SI seconds. @param {number} year */
     tropicalYearSecondsAtYear: (year: number) => number;
@@ -31,12 +25,16 @@ export function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn, a
     siderealYearSecondsAtYear: (year: number) => number;
     /** Anomalistic year of date, SI seconds. @param {number} year */
     anomalisticYearSecondsAtYear: (year: number) => number;
-    /** Axial precession period of date (equinox regression beat), years. @param {number} y */
+    /** Axial precession period of date (equinox regression beat), JULIAN years. @param {number} y */
     axialPrecessionYearsAtYear: (y: number) => number;
-    /** Perihelion (apsidal-vs-equinox) precession period of date, years. @param {number} y */
+    /** Perihelion (apsidal-vs-equinox) precession period of date, JULIAN years. @param {number} y */
     perihelionPrecessionYearsAtYear: (y: number) => number;
-    /** Apsidal precession period of date vs the fixed stars, years (hypersensitive: ±1 s ↔ ≈400 yr; the identifier keeps the historical "inclination" name). @param {number} y */
+    /** Apsidal precession period of date vs the fixed stars, JULIAN years (hypersensitive: ±1 s ↔ ≈400 yr; the identifier keeps the historical "inclination" name). @param {number} y */
     inclinationPrecessionYearsAtYear: (y: number) => number;
+    /** The movement's own general precession of date, ″ per Julian year (retrograde positive) — 360°/this ≡ axialPrecessionYearsAtYear to second order. @param {number} y */
+    generalPrecessionArcsecPerJulianYr: (y: number) => number;
+    /** The movement's own apsidal rate of date against the FIXED frame, ″ per Julian year (prograde positive) — the rate the Perihelion Longitudes Prec. cell shows for Earth; 360°/this ≡ inclinationPrecessionYearsAtYear exactly. @param {number} y */
+    apsidalRateFixedArcsecPerJulianYr: (y: number) => number;
     /** The per-cardinal layer (D6 rules: year lengths λ̇-corrected; offsets/spreads raw — μs-class/second-order there). */
     cardinal: Readonly<{
         /** @param {number} y @param {'VE'|'SS'|'AE'|'WS'} type */

package/types/model.d.ts

@@ -45,6 +45,10 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         ofDatePeriodYearsAtYear: (year: number) => number;
         /** f_S, the solar fraction of the J2000 precession torque. */
         solarShareJ2000: number;
+        /** THE PRECESSION OF DATE — the solar-torque factor 1 + f_S·([(1 − e²)/(1 − e₀²)]^(−3/2) − 1) at the eccentricity of date (1 at J2000; +3.2·10⁻⁴ on the million-year mean). The spin integration multiplies it into the secular rate. @param {number} year @returns {number} */
+        solarTorqueFactorOfDateAtYear: (year: number) => number;
+        /** THE PRECESSION OF DATE — J₂(t)/J₂₀, the oblateness of date from the GIA channel (ψ̇ ∝ J₂; 1 at J2000, bounded ±5·10⁻⁴; dJ₂/dt at J2000 = the observed satellite rate). @param {number} year @returns {number} */
+        ellipticityRatioOfDateAtYear: (year: number) => number;
         /** (a₀/a_M(t))³ — the lunar torque's growth on the recession history. @param {number} year @returns {number} */
         lunarTorqueFactorAtYear: (year: number) => number;
         /** f_S + (1 − f_S)(a₀/a_M)³ — the torque term the unit divides H_era by. @param {number} year @returns {number} */
@@ -57,11 +61,11 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         nodalPeriodYears: number;
         /** The obliquity beat 2π/(ψ̇(t) − |s₃|), years — the SHIPPED deep-time obliquity period (falsification leg 1; 41,224 at J2000, on the composed rate). @param {number} year @returns {number} */
         obliquityBeatYearsAtYear: (year: number) => number;
-        /** The apsidal (perihelion vs the stars) period from the engine-D chain's secular tangent, years — inside the published window only (the tangent is an extrapolation beyond the banked series: it turns negative at −5 Myr); null beyond. @param {number} year @returns {number|null} */
+        /** The apsidal (perihelion vs the stars) period of date, JULIAN years — the ONE movement's own apsidal rate (the scene's perihelion; anom/(anom − sid) of the one family, ≡ 360°/apsidalRateFixedArcsecPerJulianYr), inside the published window only; null beyond. Before 2026-09 this rode the planet chain's secular tangent, which parted from the scene's perihelion passages by seconds of anomalistic year within a few millennia (year-lengths.cjs). @param {number} year @returns {number|null} */
         apsidalPeriodYearsAtYear: (year: number) => number | null;
         /** T_aps(t) / T_p(t) — the apsidal period in of-date precession periods (4.33 at J2000, a reading; 0.84 … 9.9 across ±26 kyr, measured); null beyond the published window. @param {number} year @returns {number|null} */
         apsidalPerPrecessionAtYear: (year: number) => number | null;
-        /** T_peri(t) = 1/(1/T_p + 1/T_aps) — the perihelion-of-date period (equinox precession + inertial perihelion motion, frame arithmetic at every epoch), years, on the of-date T_p; null beyond the published window. @param {number} year @returns {number|null} */
+        /** T_peri(t) = 1/(1/T_p + 1/T_aps) — the perihelion-of-date period (equinox precession + inertial perihelion motion, frame arithmetic at every epoch), JULIAN years, on the of-date T_p — ≡ the one family's anom/(anom − trop); null beyond the published window. @param {number} year @returns {number|null} */
         periOfDatePeriodYearsAtYear: (year: number) => number | null;
         /** T_peri(t) / T_p(t) (0.812 at J2000 — the J2000 reading); null beyond the published window. @param {number} year @returns {number|null} */
         periOfDatePerPrecessionAtYear: (year: number) => number | null;
@@ -112,6 +116,8 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         axialPrecessionYearsAtYear: (y: number) => number;
         perihelionPrecessionYearsAtYear: (y: number) => number;
         inclinationPrecessionYearsAtYear: (y: number) => number;
+        generalPrecessionArcsecPerJulianYr: (y: number) => number;
+        apsidalRateFixedArcsecPerJulianYr: (y: number) => number;
         cardinal: Readonly<{
             yearLengthSeconds: (y: number, type: "VE" | "SS" | "AE" | "WS") => number;
             eocOffsetSeconds: (y: number, type: "VE" | "SS" | "AE" | "WS") => number;

package/types/moon/deep-modes-artifact.d.cts

@@ -9,6 +9,19 @@ export const DEEP_MODES_ARTIFACT: Readonly<{
         seed: string;
         masses: string;
         runCommand: string;
+        physics: {
+            lunarQuadrupole: {
+                qTilde: number;
+                aEmKm: number;
+                effectiveFactor: number;
+                provenance: string;
+            };
+            asteroids: {
+                ceres: number;
+                pallas: number;
+                vesta: number;
+            };
+        };
         dumpFile: string;
         dumpSha256: string;
         conservationMaxDE: number;
@@ -139,4 +152,4 @@ export const DEEP_MODES_ARTIFACT: Readonly<{
     anchorZDotPerYr: number[];
     slopeTaperYears: 27943.083333333332;
 }>;
-export const DEEP_MODES_ARTIFACT_HASH: "77e0baf779e85de5";
+export const DEEP_MODES_ARTIFACT_HASH: "cbdc477b7442c333";

package/types/phase/deep-calendar.d.cts

@@ -19,6 +19,10 @@ export type DeepCalendar = {
      * at −302,635 and grows quadratically).
      */
     jdToSIyear: (jd: number) => number;
+    /**
+     * the exact inverse of jdToSIyear.
+     */
+    siYearToJD: (ySI: number) => number;
     /**
      *   inverse of the phase machinery's cumulAtYear — null outside the table.
      */
@@ -60,6 +64,7 @@ export type DeepCalendar = {
  *   scene's precession rotations integrate on this axis; a fit on it agrees
  *   with the runtime by construction; the round-trip bias Y_SI − Y is −11.0 yr
  *   at −302,635 and grows quadratically).
+ * @property {(ySI: number) => number} siYearToJD the exact inverse of jdToSIyear.
  * @property {(targetCumul: number) => (number | null)} yearAtCumulIntegral
  *   inverse of the phase machinery's cumulAtYear — null outside the table.
  * @property {() => void} ensureCumulDaysTable builds the ∫ daysPerYear dt table