@essrt/physics 5.0.0 → 6.0.1

package/types/deltat/deep-time.d.cts

@@ -62,6 +62,25 @@ export type DeepTimeLodDeps = {
      * - gated analytic swing rate
      */
     swingLodRateAt: (year: number) => number;
+    /**
+     * - the model's DERIVED
+     * J2000 axial precession period: the certified year-length laws' beat at
+     * 2000, sid/(sid − trop) ≈ 25,771.4 yr (plan 06 S5 — one J2000 precession
+     * reading; the same anchor the hybrid obliquity self-anchors on). Read
+     * LAZILY on the first composed-rate use, never at construction: every
+     * runtime's year laws read THIS factory's bases.
+     */
+    precessionPeriodJ2000YearsFn: () => number;
+    /**
+     * - the period of Earth's orbit
+     * plane's nodal regression on the invariable plane, 1,296,000/|s₃| with s₃
+     * the dominant ζ mode of the banked deep secular modes (≈ 68,751 yr). An
+     * ORBITAL quantity (μ-tier: it does not follow Earth's spin), the reference
+     * the solar day's "ecliptic missing motion" term rides (plan 06 T2 item —
+     * formerly the spin unit's H/5, 67,063 yr, 2.5 % off and scaling with the
+     * spin at deep time). Lazy, like the precession anchor.
+     */
+    nodalPeriodYearsFn: () => number;
     /**
      * - OPTIONAL time-dependent
      * Earth-Moon angular momentum (the Driver-1½ solar channels,
@@ -99,6 +118,19 @@ export type DeepTimeLodDeps = {
  * @property {(year: number) => number} cycleLodSumAt - gated δLOD sum (incl. swing)
  * @property {(year: number) => number} swingLodAt - gated swing δLOD alone
  * @property {(year: number) => number} swingLodRateAt - gated analytic swing rate
+ * @property {() => number} precessionPeriodJ2000YearsFn - the model's DERIVED
+ *   J2000 axial precession period: the certified year-length laws' beat at
+ *   2000, sid/(sid − trop) ≈ 25,771.4 yr (plan 06 S5 — one J2000 precession
+ *   reading; the same anchor the hybrid obliquity self-anchors on). Read
+ *   LAZILY on the first composed-rate use, never at construction: every
+ *   runtime's year laws read THIS factory's bases.
+ * @property {() => number} nodalPeriodYearsFn - the period of Earth's orbit
+ *   plane's nodal regression on the invariable plane, 1,296,000/|s₃| with s₃
+ *   the dominant ζ mode of the banked deep secular modes (≈ 68,751 yr). An
+ *   ORBITAL quantity (μ-tier: it does not follow Earth's spin), the reference
+ *   the solar day's "ecliptic missing motion" term rides (plan 06 T2 item —
+ *   formerly the spin unit's H/5, 67,063 yr, 2.5 % off and scaling with the
+ *   spin at deep time). Lazy, like the precession anchor.
  * @property {(tMa: number) => number} [lEmAtAgeKgm2S] - OPTIONAL time-dependent
  *   Earth-Moon angular momentum (the Driver-1½ solar channels,
  *   recession-history.cjs). Absent → the J2000 constant, which the budget
@@ -111,11 +143,13 @@ export function createDeepTimeLod(deps: DeepTimeLodDeps): {
     lodHoursAtAge: (t_Ma: number) => number | null;
     hAtAge: (t_Ma: number) => number | null;
     lunisolarPrecessionRateArcsecPerYrAtAge: (t_Ma: number) => number | null;
+    lunisolarPrecessionPeriodYearsAtAge: (t_Ma: number) => number | null;
     lunarTorqueFactorAtAge: (tMa: number) => (number | null);
     precessionTorqueTermAtAge: (tMa: number) => (number | null);
     siderealYearSecondsAtAge: (t_Ma: number) => number;
     tropicalYearSecondsAtAge: (t_Ma: number) => number;
     tropicalYearDaysAtAge: (t_Ma: number) => number | null;
+    eclipticLodCorrectionSecondsAtAge: (t_Ma: number) => number | null;
     yearInDaysAtAge: (t_Ma: number) => number | null;
     deltaTRawSecondsAtAge: (t_Ma: number) => number;
     lodSecondsWithCorrectionsAtAge: (t_Ma: number) => number | null;

package/types/earth/precession-composed.d.cts

@@ -17,15 +17,18 @@ export function computeSolarTorqueShare(c: {
 }): number;
 /**
  * @param {{
- *   p0ArcsecPerYr: number,
+ *   p0ArcsecPerYr: number | (() => number),
  *   solarShare: number,
  *   lodSecondsAtAge: (tMa: number) => (number | null),
  *   lodJ2000Seconds: number,
  *   moonDistanceMetresAtAge: (tMa: number) => (number | null),
  *   moonDistanceJ2000Metres: number,
  *   yearToTMa: (year: number) => number,
- * }} deps - p0 = 1,296,000/(H/13) (the model's J2000 rate); lodJ2000Seconds
- *   the SAME day basis lodSecondsAtAge(0) returns.
+ * }} deps - p0 = the model's derived J2000 rate, 1,296,000/axial0 (S5; a
+ *   number or a LAZY getter — the callers' certified year laws read the
+ *   deep-time factory this is built inside, so the anchor is resolved on
+ *   the first RATE use, never at construction; the torque term needs no
+ *   anchor); lodJ2000Seconds the SAME day basis lodSecondsAtAge(0) returns.
  * @returns {{
  *   composedRateArcsecPerYrAtAge: (tMa: number) => (number | null),
  *   composedPeriodYearsAtAge: (tMa: number) => (number | null),
@@ -39,7 +42,7 @@ export function computeSolarTorqueShare(c: {
  * }}
  */
 export function createComposedPrecession(deps: {
-    p0ArcsecPerYr: number;
+    p0ArcsecPerYr: number | (() => number);
     solarShare: number;
     lodSecondsAtAge: (tMa: number) => (number | null);
     lodJ2000Seconds: number;

package/types/model.d.ts

@@ -37,11 +37,11 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         tropicalYearSecondsAtYear: (year: number) => number;
     }>;
     lunisolar: Readonly<{
-        /** The mean lunisolar precession period T_p(t), years — the composed torque rate's period. @param {number} year @returns {number} */
+        /** The mean lunisolar precession period T_p(t), years — the composed torque rate's period on the derived J2000 anchor (25,771.4 at J2000). @param {number} year @returns {number} */
         meanPeriodYearsAtYear: (year: number) => number;
-        /** The composed rate ψ̇(t) = [ω/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M)³], ″/yr. @param {number} year @returns {number} */
+        /** The composed rate ψ̇(t) = [ω/ω₀]·p₀·[f_S + (1 − f_S)(a₀/a_M)³], ″/yr (p₀ = 50.2883, derived). @param {number} year @returns {number} */
         meanRateArcsecPerYrAtYear: (year: number) => number;
-        /** The of-date period — the one-family beat inside ±2 Myr, the mean beyond. @param {number} year @returns {number} */
+        /** The of-date period — the one-family beat inside ±2 Myr, the composed mean beyond. @param {number} year @returns {number} */
         ofDatePeriodYearsAtYear: (year: number) => number;
         /** f_S, the solar fraction of the J2000 precession torque. */
         solarShareJ2000: number;
@@ -49,13 +49,19 @@ export function assembleModel(C: Readonly<Record<string, any>>, F: Readonly<Reco
         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} */
         torqueTermAtYear: (year: number) => number;
-        /** The hybrid's precession constant α = p₀ / cos ε₀, ″/yr (ε₀ the J2000 obliquity input). */
+        /** The hybrid's precession constant α = p₀ / cos ε₀, ″/yr (p₀ the derived J2000 rate, ε₀ the J2000 obliquity input; 54.81). */
         torqueConstantJ2000ArcsecPerYr: number;
+        /** |s₃|, the dominant nodal mode of Earth's orbit — the largest-amplitude ζ mode of the banked deep secular modes, ″/yr (18.85; the obliquity beat's partner). ONE home for the registry, the browser and the API. */
+        nodalModeS3ArcsecPerYr: number;
+        /** The ecliptic (nodal) precession period 1,296,000/|s₃|, years (68,751) — an orbital quantity, fixed at every epoch at the two-body level. */
+        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} */
         apsidalPeriodYearsAtYear: (year: number) => number | null;
-        /** T_aps(t) / T_p(t) — the apsidal period in mean precession periods (4.33 at J2000, a fitted reading; 0.84 … 9.5 across ±26 kyr, measured); null beyond the published window. @param {number} year @returns {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; 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), years, on the of-date T_p; 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;