@essrt/physics 25.0.1 → 26.0.1

package/types/deltat/deep-time-chain.d.ts

@@ -177,6 +177,19 @@ export function createDeepTimeChain(deps: {
         apsidalMeetsNodalSecondsAtAge: (tMa: number) => number | null;
         lunarLevelingSecondsAtAge: (tMa: number) => number | null;
     };
+    moonEccSensitivities: Readonly<{
+        perigeeEffective: number;
+        nodeEffective: number;
+        perigeePhysical: number;
+        nodePhysical: number;
+        kappaPerCy: number;
+        frameT2DegPerCy2: number;
+        generalPrecessionDegPerCy: number;
+        perigeeRateOfDateDegPerCy: number;
+        nodeRateOfDateDegPerCy: number;
+        perigeeRateIcrfDegPerCy: number;
+        nodeRateIcrfDegPerCy: number;
+    }>;
     meanMoonDistanceMetresAtAge: (tMa: number) => number | null;
     meanMoonDistanceAtAge: (tMa: number) => number | null;
     meanLodSecondsAtAge: (tMa: number) => number | null;
@@ -238,14 +251,6 @@ export function createDeepTimeChain(deps: {
     /** @param {number} tMa @param {number} periodJ2000Seconds @returns {number} */
     meanEarthPlanetSynodicPeriodAtAge: (tMa: number, periodJ2000Seconds: number) => number;
 }>;
-/**
- * Moon-channel eccentricity sensitivities (perigee/node), the [g/g₀]^s
- * exponents of the factored deep-time law — Meeus-effective (S_N moved
- * 1.0 → 1.018 with the v4 frame-attribution batch). ONE home since B2 (they
- * were literals in model.js, tools/lib/deep-time.js and script.js `_FW_MOON`).
- */
-export const MOON_ECC_SENSITIVITY_PERIGEE: 2.407;
-export const MOON_ECC_SENSITIVITY_NODE: 1.018;
 /**
  * The ΔT stack's switches. `resonator` gates the swing's δLOD terms; the
  * swing's ΔT term is gated by `resonatorDeltaT` when given (the Node twin

package/types/index.d.ts

@@ -2,6 +2,7 @@ export { createEpochPrimitives } from "./layer0/index.js";
 export { deriveEpochParams } from "./layer0/derive-params.js";
 export { createDerivedViews } from "./layer1/index.js";
 export { createDeepCalendar } from "./phase/deep-calendar.cjs";
+export { createDeepTimeChain } from "./deltat/deep-time-chain.js";
 export { createCardinalStructure } from "./cardinal/one-source-structure.cjs";
 export { createMoonMonthChain } from "./moon/month-chain.cjs";
 export { createChainCycleIntegrator } from "./chain-cycles/index.cjs";
@@ -61,7 +62,7 @@ export { GENERATED as DEFAULT_CONSTANTS, CONSTANTS_HASH, MODEL_VERSION, PREPRINT
 export { FITTED_COEFFICIENTS, COEFFICIENTS_HASH } from "./constants/index.js";
 export { GREGORIAN_START, GREGORIAN_START_JD, calendarToJD, jdToCalendar, jdToDecimalYear, dayFractionToHMS, julianCalendarToJD, jdToJulianCalendar } from "./time/calendar.cjs";
 export { DEG2RAD, RAD2DEG, ARCSEC_PER_TURN, ARCSEC_PER_RAD, wrap360, wrap180, wrap180Closed, wrap2pi, wrapPi } from "./math/angles.cjs";
-export { createDeepTimeChain, MOON_ECC_SENSITIVITY_PERIGEE, MOON_ECC_SENSITIVITY_NODE } from "./deltat/deep-time-chain.js";
+export { MOON_ECC_SENSITIVITY_PERIGEE, MOON_ECC_SENSITIVITY_NODE, deriveMoonEccSensitivities, MEEUS_ARGUMENT_RATES_DEG_PER_CY } from "./moon/ecc-sensitivities.cjs";
 export { computeMassChainJ2000, computeMoonPeriodsJ2000 } from "./earth/mass-chain.cjs";
 export { SECONDS_PER_DAY, DAYS_PER_JULIAN_YEAR, JULIAN_CENTURY_DAYS, SECONDS_PER_JULIAN_YEAR } from "./constants/units.cjs";
 export { computeEarthFrameOfDate, solveWheelAngleForLongitude, computeSceneWorldFrameJ2000, SCENE_EQUINOX_VIEW_AZIMUTH_DEG } from "./earth/frame-of-date.cjs";

package/types/model/deep-time-chain.d.ts

@@ -153,6 +153,19 @@ export function assembleDeepTimeChain(ctx: import("./derived-constants.js").Ctx)
             apsidalMeetsNodalSecondsAtAge: (tMa: number) => number | null;
             lunarLevelingSecondsAtAge: (tMa: number) => number | null;
         };
+        moonEccSensitivities: Readonly<{
+            perigeeEffective: number;
+            nodeEffective: number;
+            perigeePhysical: number;
+            nodePhysical: number;
+            kappaPerCy: number;
+            frameT2DegPerCy2: number;
+            generalPrecessionDegPerCy: number;
+            perigeeRateOfDateDegPerCy: number;
+            nodeRateOfDateDegPerCy: number;
+            perigeeRateIcrfDegPerCy: number;
+            nodeRateIcrfDegPerCy: number;
+        }>;
         meanMoonDistanceMetresAtAge: (tMa: number) => number | null;
         meanMoonDistanceAtAge: (tMa: number) => number | null;
         meanLodSecondsAtAge: (tMa: number) => number | null;

package/types/model/lunar-theory.d.ts

@@ -18,9 +18,9 @@ export function assembleLunarTheory(ctx: import("./time-axis.js").Ctx): {
             M0: number;
             MP0: number;
             F0: number;
-            LPR: number;
-            DR: number;
-            MR: number;
+            LPR: 481267.88123421;
+            DR: 445267.1114034;
+            MR: 35999.0502909;
             WDOT: number;
             NDOT: number;
             T2_W: number;
@@ -31,6 +31,12 @@ export function assembleLunarTheory(ctx: import("./time-axis.js").Ctx): {
             T2_LP_TIDAL: number;
             S_W: number;
             S_N: number;
+            S_W_PHYS: number;
+            S_N_PHYS: number;
+            W_ICRF: number;
+            N_ICRF: number;
+            P0: number;
+            T2F: number;
         };
         _rateCompletion: () => {
             rc: {

package/types/model/one-source-movement.d.ts

@@ -33,7 +33,7 @@ export function assembleOneSourceMovement(ctx: import("./year-day-lengths.js").C
         periOfDateDegAt: (year: number) => number;
         /** The hybrid's obliquity at a decimal year, degrees — the banked series inside its span, the α(t)-coupled ζ-tail integration beyond (the deep sampler grows ~0.1 s/Myr). @param {number} year @returns {number} */
         epsAt: (year: number) => number;
-        /** The hybrid's own domain, [yearFrom, yearTo] — the banked ζ series' span (plan 12 §8.5); null without a series. */
+        /** The hybrid's own domain, [yearFrom, yearTo] — the banked ζ series' span; null without a series. */
         obliquityHybridDomainYears: [number, number] | null;
         /** The precession of date's two factors at a decimal year (the factory's own — ONE home). @param {number} year @returns {{solarTorque: number, ellipticity: number}} */
         precessionOfDateFactorsAt: (year: number) => {

package/types/model/surface.d.ts

@@ -44,7 +44,7 @@ export function assembleSurface(ctx: import("./cardinal-instants.js").Ctx): Read
     lunisolar: Readonly<{
         /** The mean lunisolar precession period T_p(t), years — 1,296,000/the precession RATE of date (25,771.6 at J2000): THE published deep-time clock. @param {number} year @returns {number} */
         meanPeriodYearsAtYear: (year: number) => number;
-        /** THE PRECESSION RATE OF DATE ψ̇(t) = [ω/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M)³] × cos ε̄(t)/cos ε̄(0), ″/yr — the constant's scaling on the spin and the lunar torque times the cited mean obliquity's cosine (plan 12 §8 option (a); p₀ = 50.2880, the IAU 2006 input). @param {number} year @returns {number} */
+        /** THE PRECESSION RATE OF DATE ψ̇(t) = [ω/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M)³] × cos ε̄(t)/cos ε̄(0), ″/yr — the constant's scaling on the spin and the lunar torque times the cited mean obliquity's cosine (p₀ = 50.2880, the IAU 2006 input). @param {number} year @returns {number} */
         meanRateArcsecPerYrAtYear: (year: number) => number;
         /** THE PRECESSION CONSTANT OF DATE as a rate on the J2000 obliquity, ″/yr — the spin and torque scaling alone (the hybrid obliquity integration's integrand: its own ŝ·n̂ carries the cosine). @param {number} year @returns {number} */
         torqueRateArcsecPerYrAtYear: (year: number) => number;
@@ -123,9 +123,9 @@ export function assembleSurface(ctx: import("./cardinal-instants.js").Ctx): Read
     earth: Readonly<{
         perihelionLongitudeDeg: (year: number) => number;
         obliquityDeg: (year: number) => number;
-        /** The hybrid's own domain, [yearFrom, yearTo] — the banked ζ series' span, ±10 Myr about J2000 (plan 12 §8.5: the published surfaces serve `obliquityDeg` inside it only; null without a series, i.e. a bare createModel()). */
+        /** The hybrid's own domain, [yearFrom, yearTo] — the banked ζ series' span, ±10 Myr about J2000 (the published surfaces serve `obliquityDeg` inside it only; null without a series, i.e. a bare createModel()). */
         obliquityHybridDomainYears: [number, number] | null;
-        /** The cited mean obliquity of date ε̄(t) (Farhat et al. 2022's column; plan 12 §8 option (a)), degrees — a Myr-scale mean; NaN beyond 4,429 Ma. @param {number} year @returns {number} */
+        /** The cited mean obliquity of date ε̄(t) (Farhat et al. 2022's column), degrees — a Myr-scale mean; NaN beyond 4,429 Ma. @param {number} year @returns {number} */
         meanObliquityCitedDeg: (year: number) => number;
         eccentricity: (year: number) => number;
         inclinationDeg: (year: number) => number;
@@ -184,6 +184,19 @@ export function assembleSurface(ctx: import("./cardinal-instants.js").Ctx): Read
         distanceKmAtYear: (year: number) => number;
         siderealMonthDaysAtYear: (year: number) => number;
         synodicMonthDays: number;
+        eccSensitivities: Readonly<{
+            perigeeEffective: number;
+            nodeEffective: number;
+            perigeePhysical: number;
+            nodePhysical: number;
+            kappaPerCy: number;
+            frameT2DegPerCy2: number;
+            generalPrecessionDegPerCy: number;
+            perigeeRateOfDateDegPerCy: number;
+            nodeRateOfDateDegPerCy: number;
+            perigeeRateIcrfDegPerCy: number;
+            nodeRateIcrfDegPerCy: number;
+        }>;
         lonDegAtJD: (jd: number) => number;
         betaDegAtJD: (jd: number) => number;
         distanceKmAtJD: (jd: number) => number;
@@ -191,7 +204,7 @@ export function assembleSurface(ctx: import("./cardinal-instants.js").Ctx): Read
     eclipse: Readonly<{
         sunLonDegAtJD: (jd: number) => number;
         sunLonCompletedDegAtJD: (jd: number) => number;
-        /** The APPARENT Sun (completed geometric − the derived aberration constant + the leading nutation terms on the model's own arguments) — the quantity Horizons' observer ecliptic longitude is; the cardinal instants are its crossings (plan 06 R1/I1). @param {number} jdUT @returns {number} */
+        /** The APPARENT Sun (completed geometric − the derived aberration constant + the leading nutation terms on the model's own arguments) — the quantity Horizons' observer ecliptic longitude is; the cardinal instants are its crossings. @param {number} jdUT @returns {number} */
         sunApparentLonDegAtJD: (jdUT: number) => number;
         findLunarInRange: (jdStart: number, jdEnd: number) => {
             jd: number;

package/types/model.d.ts

@@ -145,6 +145,19 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         distanceKmAtYear: (year: number) => number;
         siderealMonthDaysAtYear: (year: number) => number;
         synodicMonthDays: number;
+        eccSensitivities: Readonly<{
+            perigeeEffective: number;
+            nodeEffective: number;
+            perigeePhysical: number;
+            nodePhysical: number;
+            kappaPerCy: number;
+            frameT2DegPerCy2: number;
+            generalPrecessionDegPerCy: number;
+            perigeeRateOfDateDegPerCy: number;
+            nodeRateOfDateDegPerCy: number;
+            perigeeRateIcrfDegPerCy: number;
+            nodeRateIcrfDegPerCy: number;
+        }>;
         lonDegAtJD: (jd: number) => number;
         betaDegAtJD: (jd: number) => number;
         distanceKmAtJD: (jd: number) => number;

package/types/moon/arguments.d.cts

@@ -16,6 +16,7 @@ export type MoonArgsDeg = {
  *     elpEarthFigureJ2ArcsecPerCy2: number,
  *     elpGeneralPrecessionPA_T2ArcsecPerCy2: number,
  *     eccE0: number,
+ *     moonEccSensitivities: ReturnType<typeof import('./ecc-sensitivities.cjs').deriveMoonEccSensitivities>,
  *   },
  *   fns: {
  *     eccAt: (tYr: number) => number,
@@ -56,6 +57,7 @@ export function createMoonArguments({ constants, fns }: {
         elpEarthFigureJ2ArcsecPerCy2: number;
         elpGeneralPrecessionPA_T2ArcsecPerCy2: number;
         eccE0: number;
+        moonEccSensitivities: ReturnType<typeof import("./ecc-sensitivities.cjs").deriveMoonEccSensitivities>;
     };
     fns: {
         eccAt: (tYr: number) => number;
@@ -86,9 +88,9 @@ export function createMoonArguments({ constants, fns }: {
         M0: number;
         MP0: number;
         F0: number;
-        LPR: number;
-        DR: number;
-        MR: number;
+        LPR: 481267.88123421;
+        DR: 445267.1114034;
+        MR: 35999.0502909;
         WDOT: number;
         NDOT: number;
         T2_W: number;
@@ -99,6 +101,12 @@ export function createMoonArguments({ constants, fns }: {
         T2_LP_TIDAL: number;
         S_W: number;
         S_N: number;
+        S_W_PHYS: number;
+        S_N_PHYS: number;
+        W_ICRF: number;
+        N_ICRF: number;
+        P0: number;
+        T2F: number;
     };
     _rateCompletion: () => {
         rc: {

package/types/moon/deep-ecc-channel.d.cts

@@ -17,7 +17,8 @@
  *
  * — ż_J2000 is the run's own z-rate at J2000 (script-written into the
  * artifact by generate.mjs from the banked series: "we are the source"),
- * T_q the H/12 taper the arguments' rate anchors ride. Inside the canon it
+ * T_q the arguments' secular-completion taper (SECULAR_COMPLETION_TAPER_YEARS,
+ * the one the rate anchors ride). Inside the canon it
  * restores the run's slope exactly; beyond ±T_q it is a constant vector
  * offset of ~3e-3, the size of the table's own spectrum-class scatter
  * there, never a growing polynomial. Absent the two fields (an older

package/types/moon/ecc-sensitivities.d.cts

@@ -0,0 +1,48 @@
+/** Meeus Ch. 47 mean-argument rates, degrees per Julian century (the J2000 polynomials' linear terms). */
+export const MEEUS_ARGUMENT_RATES_DEG_PER_CY: Readonly<{
+    Lp: 481267.88123421;
+    D: 445267.1114034;
+    M: 35999.0502909;
+    Mp: 477198.8675055;
+    F: 483202.0175233;
+}>;
+/** The Meeus-effective pair: Meeus-exact of-date T² by construction (the shipped record values). */
+export const MOON_ECC_SENSITIVITY_PERIGEE: 2.407;
+export const MOON_ECC_SENSITIVITY_NODE: 1.018;
+/**
+ * The physical pair and the rate split, derived from the shipped constants.
+ * @param {{
+ *   eccentricityJ2000: number,
+ *   eccentricityDotJ2000: number,
+ *   frameT2ArcsecPerCy2: number,
+ *   generalPrecessionArcsecPerYr: number,
+ * }} inputs — e₀ (the catalog secular value the Taylor checks use); ė₀ per
+ *   Julian century; the general precession's T² coefficient the of-date Meeus
+ *   polynomials embed (Lieske 1976, the ELP frame), ″/cy²; p₀, the IAU 2006 p_A
+ *   linear rate, ″/yr
+ * @returns {Readonly<{
+ *   perigeeEffective: number, nodeEffective: number,
+ *   perigeePhysical: number, nodePhysical: number,
+ *   kappaPerCy: number, frameT2DegPerCy2: number, generalPrecessionDegPerCy: number,
+ *   perigeeRateOfDateDegPerCy: number, nodeRateOfDateDegPerCy: number,
+ *   perigeeRateIcrfDegPerCy: number, nodeRateIcrfDegPerCy: number,
+ * }>}
+ */
+export function deriveMoonEccSensitivities({ eccentricityJ2000, eccentricityDotJ2000, frameT2ArcsecPerCy2, generalPrecessionArcsecPerYr }: {
+    eccentricityJ2000: number;
+    eccentricityDotJ2000: number;
+    frameT2ArcsecPerCy2: number;
+    generalPrecessionArcsecPerYr: number;
+}): Readonly<{
+    perigeeEffective: number;
+    nodeEffective: number;
+    perigeePhysical: number;
+    nodePhysical: number;
+    kappaPerCy: number;
+    frameT2DegPerCy2: number;
+    generalPrecessionDegPerCy: number;
+    perigeeRateOfDateDegPerCy: number;
+    nodeRateOfDateDegPerCy: number;
+    perigeeRateIcrfDegPerCy: number;
+    nodeRateIcrfDegPerCy: number;
+}>;

package/types/moon/month-chain.d.cts

@@ -23,13 +23,10 @@
  *   eccentricity channel's modulation; constants are the engine's derived
  *   J2000 anchors (both engines derive them from the same JSON).
  *   fns.precessionPeriodYearsAtAge is the engine's lunisolar precession period
- *   of date T_p(t), Julian years (the composed rate's period; plan 09 step
- *   5c — the tropical month's equinox turn, formerly the counter's 13/H(t),
- *   0.086 % slow: the chain's J2000 tropical month sat 5.85 ms above the IAU
- *   input, the 1/T_p form is within 0.1 ms of it).
+ *   of date T_p(t), Julian years (the composed rate's period — the tropical
+ *   month's equinox turn; within 0.1 ms of the IAU input at J2000).
  *   fns.distanceMetresAtAge is the engine's CITED distance history
- *   (deltat/moon-distance-history; plan 12 option A — the internal quartic
- *   twin this chain carried until then is deleted; null beyond the span).
+ *   (deltat/moon-distance-history; null beyond the span).
  */
 export function createMoonMonthChain({ constants, fns }: {
     constants: {