@essrt/physics 11.0.1 → 11.1.0

package/package.json

@@ -1,9 +1,9 @@
 {
   "name": "@essrt/physics",
-  "version": "11.0.1",
+  "version": "11.1.0",
   "description": "ESSRT (Expanding Solar System Resonance Theory) physics core — 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": "v16.0",
+    "modelVersion": "v16.1",
     "_note": "The model identity this package version ships (IP-unified-architecture §10: package semver tracks the API axis; the model version + CONSTANTS_HASH/COEFFICIENTS_HASH in src/constants/ carry the scientific identity). The version-pinning gate enforces this pairing against public/input/model-version.json."
   },
   "type": "module",

package/src/constants/generated.js

@@ -36,7 +36,7 @@
 export const CONSTANTS_HASH = "726b10e23d8ace09";
 
 /** Model version label — single source: public/input/model-version.json (§10 two-axis scheme). */
-export const MODEL_VERSION = "v16.0";
+export const MODEL_VERSION = "v16.1";
 
 /** Canonical preprint DOI — single source: public/input/model-version.json. */
 export const PREPRINT_DOI = "10.21203/rs.3.rs-8758810/v4";

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

@@ -0,0 +1,24 @@
+/**
+ * GENERATED — do not edit. Regenerate:
+ *   node tools/constants/generate.mjs --write
+ *
+ * Source: data/nbody-secular-series.json — the Sun planetary completion's
+ * embedded data (plan 06 I3): planetLamDotJ2000DegPerYr, the planets' J2000
+ * sidereal mean motions from the model's own ±10-Myr run (the completion's
+ * CARRIERS; the former carriers were the planet records' rounded of-date
+ * periods minus p₀ — Venus 4.1″/yr, Jupiter 1.7″/yr off the run's own
+ * motion); and eccVectors, the e·e^{iϖ} of Venus, Earth, Mars and Jupiter on
+ * a 1-kyr grid over ±50 kyr (ecliptic J2000) — the D'Alembert long-period
+ * rows' amplitude and phase of date. Earth's carrier stays the framework
+ * sidereal year (the D6 ratio-only doctrine for the run's absolute Earth
+ * rate). Two gates guard the chain: check:artifacts pins artifact ↔ engine;
+ * generate.mjs check mode pins this embed ↔ artifact. CJS for the .cjs
+ * completion consumer.
+ */
+'use strict';
+
+const SUN_COMPLETION_ARTIFACT_HASH = "5ad17b275a003f6e";
+
+const SUN_COMPLETION_ARTIFACT = Object.freeze({"planetLamDotJ2000DegPerYr":{"mercury":1494.725744748,"venus":585.177909613,"mars":191.402970848,"jupiter":30.348892629,"saturn":12.221520382,"uranus":4.285090582,"neptune":2.18458998},"eccVectors":{"t0Yr":-50000,"stepYr":1000,"bodies":{"venus":{"zQ":[0.0279902,0.0270462,0.0260531,0.0250136,0.0239295,0.0228034,0.0216389,0.0204406,0.0192142,0.0179653,0.0166999,0.0154225,0.0141391,0.0128553,0.0115784,0.0103157,0.009074,0.0078599,0.006678,0.0055331,0.0044289,0.00337,0.0023604,0.0014037,0.0005018,-0.0003439,-0.0011327,-0.0018642,-0.0025389,-0.0031565,-0.0037171,-0.0042205,-0.0046666,-0.0050566,-0.0053911,-0.0056726,-0.0059018,-0.0060797,-0.0062062,-0.0062821,-0.0063089,-0.0062889,-0.0062255,-0.0061207,-0.0059782,-0.0057996,-0.0055879,-0.0053462,-0.0050788,-0.0047909,-0.0044872,-0.0041718,-0.0038483,-0.0035199,-0.0031899,-0.0028619,-0.0025389,-0.0022246,-0.0019216,-0.0016326,-0.0013585,-0.0011005,-0.0008588,-0.0006348,-0.0004296,-0.0002445,-0.0000808,0.000062,0.0001836,0.0002854,0.0003676,0.00043,0.000472,0.0004936,0.0004954,0.000479,0.0004465,0.0003999,0.0003413,0.0002721,0.000194,0.000109,0.00002,-0.0000692,-0.0001549,-0.0002339,-0.0003033,-0.0003614,-0.000406,-0.0004346,-0.0004451,-0.0004351,-0.0004024,-0.0003454,-0.0002637,-0.0001573,-0.000027,0.0001276,0.0003064,0.0005097,0.0007375],"zP":[0.018044,0.0190048,0.0199026,0.0207345,0.0214985,0.0221921,0.0228133,0.023361,0.0238349,0.0242364,0.0245664,0.0248268,0.0250179,0.0251396,0.0251924,0.0251771,0.0250952,0.0249487,0.0247395,0.024469,0.0241382,0.0237481,0.0233005,0.022799,0.0222471,0.0216492,0.0210095,0.0203323,0.0196229,0.0188852,0.0181246,0.0173466,0.0165573,0.0157631,0.0149696,0.0141817,0.0134033,0.0126374,0.0118877,0.0111579,0.0104519,0.0097728,0.0091226,0.0085024,0.0079131,0.0073554,0.0068302,0.0063401,0.005886,0.0054692,0.0050896,0.0047468,0.0044407,0.004171,0.0039381,0.0037416,0.0035814,0.0034562,0.0033649,0.0033046,0.0032728,0.0032669,0.0032846,0.0033239,0.0033824,0.0034567,0.0035432,0.0036382,0.0037379,0.0038402,0.0039424,0.0040437,0.0041416,0.0042336,0.0043171,0.0043902,0.0044512,0.0044998,0.0045358,0.0045583,0.0045673,0.0045609,0.0045385,0.0044992,0.0044432,0.0043714,0.0042846,0.0041835,0.0040683,0.0039395,0.0037975,0.0036443,0.0034811,0.003311,0.0031368,0.0029615,0.0027871,0.002616,0.0024505,0.0022939,0.0021498]},"earth":{"zQ":[0.013227267,0.013012506,0.012834994,0.012696287,0.012597963,0.012540754,0.01252666,0.012549935,0.012604882,0.012690231,0.012800154,0.012934505,0.013090442,0.013266166,0.013451544,0.013642315,0.013829021,0.014004751,0.01416527,0.014306209,0.014427996,0.014517185,0.014571021,0.014583484,0.01454318,0.014451082,0.014306069,0.014099616,0.013836713,0.013514405,0.013124137,0.012669856,0.012147851,0.011555754,0.010903934,0.010193697,0.009426983,0.0086143,0.007749231,0.006839273,0.005893806,0.004915481,0.003921086,0.002917354,0.001909786,0.000909778,-0.000076816,-0.001048344,-0.001988304,-0.002888344,-0.003740032,-0.004530056,-0.005260301,-0.005923768,-0.006513461,-0.007030299,-0.00746818,-0.007819964,-0.008084188,-0.008260875,-0.00834821,-0.008350985,-0.008275204,-0.008118387,-0.007884356,-0.00757642,-0.007194194,-0.006743231,-0.006236619,-0.0056763,-0.005070356,-0.004430986,-0.00375356,-0.003046685,-0.002319076,-0.00157522,-0.000827409,-0.000082833,0.000652857,0.001372028,0.002074185,0.002755866,0.003406986,0.00402808,0.00460898,0.005142789,0.005637646,0.006085386,0.006487544,0.006852803,0.00717117,0.007444843,0.007677336,0.007863353,0.008009495,0.008120645,0.008198929,0.008249784,0.00827192,0.008262513,0.008227021],"zP":[0.00629199,0.006157176,0.006021018,0.005888044,0.005762291,0.005648055,0.005553525,0.00548229,0.005439151,0.005424431,0.005438065,0.005482258,0.005560686,0.005676894,0.005837022,0.006042534,0.006293455,0.006590047,0.006927276,0.007306088,0.007729922,0.008196139,0.008708393,0.009259394,0.009843797,0.010459935,0.01109414,0.011746079,0.012411664,0.013078877,0.013747935,0.014412012,0.015051841,0.015665551,0.016243589,0.016770777,0.017250892,0.017677425,0.018040396,0.018339722,0.018561665,0.018702501,0.018762478,0.018735718,0.018628747,0.018440526,0.018169036,0.017819371,0.017387029,0.01687065,0.016283928,0.015630793,0.014913266,0.014144751,0.013327617,0.012464356,0.011566053,0.010638047,0.009687504,0.008727497,0.007768976,0.006816931,0.005878676,0.004964821,0.004074941,0.003218725,0.00240902,0.001647618,0.000944991,0.000305919,-0.000270521,-0.000782068,-0.00122626,-0.001604285,-0.001905386,-0.00213232,-0.002287543,-0.002367048,-0.002382224,-0.002334107,-0.002223921,-0.002058154,-0.001834453,-0.001554207,-0.001229972,-0.000862036,-0.000458958,-0.000030798,0.000426187,0.000903095,0.001397238,0.00191096,0.002431882,0.002957636,0.003487072,0.004008594,0.004525642,0.005041098,0.005550956,0.006055939,0.006554044]},"mars":{"zQ":[0.0433563,0.03467,0.0258658,0.0170192,0.0082073,-0.0004939,-0.0090076,-0.0172605,-0.0251854,-0.0327154,-0.0397897,-0.0463513,-0.0523471,-0.0577309,-0.0624622,-0.0665068,-0.0698416,-0.0724452,-0.0743037,-0.0754094,-0.0757598,-0.0753669,-0.0742434,-0.07241,-0.0698931,-0.066724,-0.0629392,-0.058577,-0.053681,-0.0482992,-0.0424818,-0.0362828,-0.0297564,-0.0229575,-0.0159422,-0.0087638,-0.0014737,0.0058727,0.0132213,0.0205196,0.0277185,0.0347724,0.0416363,0.0482678,0.0546244,0.0606678,0.0663629,0.071674,0.0765709,0.0810267,0.085017,0.0885174,0.0915066,0.0939667,0.0958802,0.0972347,0.0980194,0.0982305,0.0978673,0.0969283,0.0954136,0.0933238,0.0906673,0.0874561,0.0837052,0.0794336,0.0746638,0.0694224,0.0637356,0.0576332,0.0511496,0.0443223,0.0371937,0.0298083,0.0222154,0.0144652,0.0066091,-0.0013021,-0.0092175,-0.017081,-0.0248374,-0.0324313,-0.0398113,-0.046924,-0.0537171,-0.060141,-0.066145,-0.0716831,-0.0767087,-0.0811772,-0.0850485,-0.0882866,-0.0908567,-0.0927255,-0.0938676,-0.0942613,-0.0938902,-0.0927442,-0.0908235,-0.0881389,-0.0847045],"zP":[0.0897466,0.0918932,0.09322,0.093721,0.0934006,0.0922683,0.090343,0.0876532,0.0842298,0.0801091,0.0753335,0.0699511,0.064016,0.0575862,0.0507249,0.0435004,0.0359792,0.0282273,0.0203139,0.012307,0.0042791,-0.0036983,-0.0115608,-0.0192412,-0.026677,-0.0338101,-0.0405856,-0.0469554,-0.0528697,-0.0582895,-0.0631778,-0.0675048,-0.0712468,-0.0743843,-0.0769046,-0.0787988,-0.0800608,-0.0806874,-0.0806826,-0.0800561,-0.0788225,-0.0769976,-0.0745995,-0.0716471,-0.0681631,-0.0641733,-0.0597044,-0.0547853,-0.0494511,-0.0437376,-0.0376793,-0.031313,-0.0246763,-0.0178092,-0.0107514,-0.0035454,0.0037629,0.0111266,0.0185008,0.0258424,0.0331059,0.0402454,0.0472133,0.0539669,0.0604642,0.0666637,0.0725249,0.0780099,0.0830837,0.0877107,0.0918571,0.0954933,0.0985927,0.1011309,0.1030866,0.104441,0.1051785,0.1052838,0.1047436,0.1035491,0.1016984,0.0991961,0.0960507,0.0922751,0.0878862,0.0829058,0.0773604,0.0712844,0.0647174,0.0577059,0.0503057,0.0425739,0.0345707,0.0263598,0.0180075,0.0095798,0.0011454,-0.0072255,-0.0154599,-0.0234891,-0.0312472]},"jupiter":{"zQ":[0.0409753,0.0433687,0.0456423,0.0477984,0.0498037,0.0516436,0.0532903,0.0546913,0.0558342,0.0567036,0.0573161,0.057663,0.0577405,0.0575449,0.0570533,0.056287,0.0552536,0.0540004,0.0525615,0.050946,0.0491766,0.0472552,0.0452242,0.0431052,0.0409485,0.0388261,0.0367588,0.0347878,0.0329308,0.0312513,0.0297929,0.0285898,0.0276805,0.027063,0.026749,0.0267075,0.0269516,0.027497,0.0283258,0.0294268,0.030747,0.0322622,0.0338998,0.0356259,0.0374175,0.0392265,0.0410199,0.0427079,0.0442559,0.0456258,0.0468048,0.0477873,0.0485481,0.0490863,0.0493524,0.04932,0.0489817,0.0483498,0.0474508,0.0462864,0.0448821,0.0432313,0.0413402,0.0392239,0.036927,0.0345012,0.0319648,0.0293507,0.026667,0.0239268,0.0211542,0.0183921,0.0157132,0.0131488,0.0107451,0.0085068,0.0064494,0.0046131,0.0030061,0.00168,0.0006473,-0.0000718,-0.0005082,-0.0006718,-0.0005465,-0.0001609,0.0004972,0.0013985,0.0025113,0.0037665,0.0050863,0.0064365,0.007773,0.0090802,0.010283,0.0113433,0.0122282,0.0128782,0.0132824,0.013423,0.0133097],"zP":[-0.0343967,-0.0331679,-0.0316586,-0.0298935,-0.027898,-0.025726,-0.0233818,-0.0208861,-0.0182549,-0.0154933,-0.0126481,-0.0097471,-0.0068584,-0.004012,-0.0012369,0.0014519,0.0040408,0.0064707,0.0087255,0.0107573,0.0125406,0.0140643,0.0153359,0.0163587,0.017093,0.0175361,0.0176648,0.0174823,0.017015,0.0163061,0.0154165,0.0143655,0.0131777,0.0118658,0.0104958,0.0091211,0.007783,0.0065585,0.0054776,0.0045665,0.0038381,0.0033371,0.0031166,0.0032035,0.0036307,0.0043725,0.0054131,0.0067327,0.0082968,0.0101046,0.0121388,0.0143914,0.0168005,0.0193307,0.0219615,0.0246359,0.02735,0.0300702,0.0327765,0.0354121,0.0379233,0.0402887,0.0424719,0.0444815,0.0462815,0.0478736,0.0492324,0.0503175,0.0511262,0.0516474,0.0519031,0.0518759,0.0515854,0.0510289,0.0501942,0.0491163,0.0478163,0.0463574,0.0447717,0.0430857,0.041321,0.0394987,0.03767,0.0358646,0.0341508,0.0325953,0.0312215,0.0300563,0.029115,0.028451,0.0280902,0.028059,0.0283637,0.0289847,0.0298923,0.0310308,0.0323937,0.0339708,0.0357433,0.0376723,0.0397061]}}},"meta":{"dumpSha256":"686d6af51a689f7bd64ac1e2d8b26c243a4beccdacd3167adfeeae078e8b990f","source":"data/nbody-secular-series.json verdict.planetLamDot (J2000 sidereal mean motions) + bodies.<X>.zQ/zP at the 1-kyr nodes (ecliptic J2000, verbatim)"}});
+
+module.exports = { SUN_COMPLETION_ARTIFACT, SUN_COMPLETION_ARTIFACT_HASH };

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

@@ -107,6 +107,8 @@
 
 'use strict';
 
+const { SUN_COMPLETION_ARTIFACT } = require('./sun-completion-artifact.cjs');
+
 /** J2000 mean-longitude phase anchors (deg), body order Mercury, Venus,
  *  Earth (EMB), Mars, Jupiter, Saturn — declared epoch constants
  *  (header (ii)); the RATES are injected (framework-derived). */
@@ -130,25 +132,25 @@ const ARG_D0 = 297.8501921;
  * detrended 0.033 min). Individual coefficients redistribute between the
  * near-degenerate ±M sideband families; only the composed function ships.
  *
- * PLAN 06 I2 — the LONG-PERIOD rows (the last two): the Earth–Mars–
- * Jupiter long inequality 4λ_E − 8λ_Ma + 3λ_J (carrier period 1801 yr;
- * the classical ~1783-yr term of the solar theory) and the Venus–Earth
- * term 8λ_V − 13λ_E (238 yr), DERIVED on the model's own Wisdom–Holman
- * engine (tools/explore/nbody-wh.mjs, the chain artifact's integrator) on
- * the campaign's ONE Horizons J2000 seed, Sun + eight planets, DE440
- * masses, 1PN, order 2, dt 2 d, −5100..+1100 yr: the EMB's osculating
- * mean longitude Ω+ω+M in yearly means, LSQ on [1, T, T²] + cos/sin of
- * the arguments on THESE carriers (tools/explore/i2-long-inequality.mjs).
- * Measured: 6.27″ / 1.84″; step-converged (dt 1 d: 6.272 / 1.869); the
- * integration reproduces DE441's Earth longitude to 0.68″ sd over 6000
- * yr with 0.03″ left at the 1783-yr argument. Neither could come from the
- * 200-yr D2 window (the term's in-window ramp went to the projected-out
- * secular basis; the 240-yr term folded into it). WHY they were missing:
- * the certified Sun integrates a smooth tropical year for its mean
- * longitude; against Horizons its residual carried exactly this 6″ ripple
- * (plan 06 I1 analysis) — with the same phase in the sidereal residual, a
- * perturbation of Earth's mean motion, not a frame effect. Sub-0.2″
- * candidates left out (2J−5S 0.17″, 5V−8E 0.17″: window-dependent phase).
+ * PLAN 06 I2 → I3 — the LONG-PERIOD rows live in LONG_PERIOD_TERMS below
+ * (D'Alembert form on the e-vectors of date, derived over ±20 kyr; I3). The
+ * I2 record: the Earth–Mars–Jupiter long inequality 4λ_E − 8λ_Ma + 3λ_J
+ * (the classical ~1783-yr term) and the Venus–Earth term 8λ_V − 13λ_E
+ * (≈239 yr) were first DERIVED as constant rows on the model's own
+ * Wisdom–Holman engine over −5100..+1100 yr (6.27″ / 1.84″;
+ * tools/explore/i2-long-inequality.mjs). Neither could come from the 200-yr
+ * D2 window (the term's in-window ramp went to the projected-out secular
+ * basis; the 240-yr term folded into it). WHY they were missing: the
+ * certified Sun integrates a smooth tropical year for its mean longitude;
+ * against Horizons its residual carried exactly this 6″ ripple (plan 06 I1
+ * analysis) — with the same phase in the sidereal residual, a perturbation
+ * of Earth's mean motion, not a frame effect. Sub-0.2″ candidates left out
+ * (2J−5S 0.17″, 5V−8E 0.17″: window-dependent phase). The CARRIERS moved
+ * with I3 from the planet records' rounded of-date periods minus p₀ to the
+ * model's own banked J2000 sidereal mean motions (computeCarrierRatesDegPerCy;
+ * the 70 short-period literals are KEPT — the two carrier sets share the
+ * J2000 anchors and part by ≤ 12″/yr·|k|, i.e. ≤ 0.3°·|k| at the 200-yr
+ * window's edges, the same class as the I2 of-date→sidereal move).
  * @type {Array<[number[], number[], number, number]>}
  */
 const TERMS = [
@@ -222,28 +224,128 @@ const TERMS = [
   [[0, 0, 1, -1, 0, 0], [0, 0, 1, 0, 0, 0], 0.0188, 0.0523],
   [[0, 2, -2, 0, 0, 0], [0, 0, 1, 0, 0, 0], -0.0421, -0.0340],
   [[0, 0, 1, -1, 0, 0], [0, 0, -1, 0, 0, 0], -0.0353, -0.0379],
-  // plan 06 I2 — the long-period rows (see the doc comment above)
-  [[0, 0, 4, -8, 3, 0], [0, 0, 0, 0, 0, 0], 6.2314, -0.7168],
-  [[0, 8, -13, 0, 0, 0], [0, 0, 0, 0, 0, 0], 1.8389, 0.1419],
+  // plan 06 I3 — the Venus–Earth term 8λ_V − 13λ_E (≈239 yr on the banked
+  // carriers): a CONSTANT row, derived on the model's own run over ±3100 yr
+  // (i3-long-period-dalembert.mjs `2 20000 3100 4`; the I2 row read 1.84/0.14″
+  // on the record carriers over −5100..+1100). The term is formally FIFTH order
+  // (Σk = −5) and resonant (divisor 1.5°/yr: e⁵/ν² reaches the arcsecond), so it
+  // is a CLUSTER of lines (ϖ/Ω multipliers summing to +5) beating on ~20 kyr —
+  // no single-ϖ D'Alembert form holds (the {V,E} form shifted the Sun +0.9″ at
+  // J2000 against Horizons' modern window; this row 0.74 → 0.83″ mean with the
+  // scatter 0.94 → 0.78″), and the row is a LOCAL description of the
+  // Horizons-certified era (≈1″ beyond it).
+  [[0, 8, -13, 0, 0, 0], [0, 0, 0, 0, 0, 0], 1.4825, 0.9436],
+];
+
+/**
+ * PLAN 06 I3 — the LONG INEQUALITY in D'ALEMBERT form (superseding the
+ * constant-amplitude row of I2, 6.23/−0.72″): the Earth–Mars–Jupiter long
+ * inequality 4λ_E − 8λ_Ma + 3λ_J (≈1783 yr) as a sum over bodies X of
+ *     e_X(t)·[a_X·cos(θ + ϖ_X(t)) + b_X·sin(θ + ϖ_X(t))]
+ * with e_X, ϖ_X the e-vectors OF DATE (the embedded artifact, ecliptic J2000;
+ * D'Alembert: Σk_λ + Σj_ϖ = 0 and the argument has Σk = −1, so +ϖ_X).
+ * WHY (measured 2026-10 on the model's inertial Sun against DE441 in the fixed
+ * J2000 frame, tools/explore/sun-inertial-vs-de441.cjs): the constant rows left
+ * 0.08″ at the 1783-yr period inside ±3000 yr and 3.7″ outside on both sides —
+ * the inequality is first order in the eccentricities and its composed
+ * amplitude/phase ride the e-vectors of date (the model's own N-body shows the
+ * line growing 3.2″ → 9.0″ from −18 to +14 kyr). DERIVED on the model's own
+ * Wisdom–Holman run over ±20,000 yr (tools/explore/i3-long-period-dalembert.mjs
+ * `2 20000 9100 4`: yearly means of the EMB's osculating mean longitude, a
+ * degree-4 secular detrend — diagnostic — plus these columns on the SAME
+ * embedded e-vectors; dt 2 d, step-converged against dt 1 d): the 1783-yr band
+ * left in the residual 0.05″ inside ±9100 and 0.07/0.05″ outside, the composed
+ * J2000 amplitude 6.88″ on ±20 kyr vs 6.89″ on ±9100 (window-independent; the
+ * constant row read 6.12 vs 6.89), and the ±9100-fitted rows evaluated over
+ * ±20 kyr leave 0.28/0.53″ (rms Δ 0.57″). Against Horizons' modern window the
+ * row moves the Sun by −0.36″ at J2000 with the scatter 0.94 → 0.89″ (the I2
+ * constant row was derived off-centre, −5100..+1100). The Venus–Earth term is
+ * NOT of this form — see its constant row in TERMS. The individual a_X, b_X
+ * are NOT physically separable (the three perihelia
+ * rotate only 33–80° per 9 kyr, so the columns are near-collinear — the first
+ * cut, with the WRONG sign θ − ϖ_X, "fitted" ±9 kyr with coefficients of
+ * thousands of ″/e and failed the window test); only the composed function
+ * ships, as for the ±M sideband pairs. Beyond ±50 kyr the e-vectors are held at
+ * the grid's ends (bounded, the former constant-row class).
+ * Extraction-native sign (N-body − smooth), negated in the evaluator.
+ * @type {Array<[number[], number, number, number]>} [l-multipliers, bodyIndex (0 Me … 5 S), a (″/e), b (″/e)]
+ */
+const LONG_PERIOD_TERMS = [
+  [[0, 0, 4, -8, 3, 0], 2, -785.9521, 591.0878],   // 4E−8Ma+3J · earth
+  [[0, 0, 4, -8, 3, 0], 3, -80.4721, 77.5236],     // 4E−8Ma+3J · mars
+  [[0, 0, 4, -8, 3, 0], 4, 159.5521, -303.8727],   // 4E−8Ma+3J · jupiter
 ];
+/** body index → the embedded e-vector series' key */
+const ECC_BODY_KEY = [null, 'venus', 'earth', 'mars', 'jupiter', null];
 
 const D2R = Math.PI / 180;
 
+/**
+ * The e-vector (e·cos ϖ, e·sin ϖ; ecliptic J2000) of a body at T Julian
+ * centuries TT from J2000, from the embedded 1-kyr grid (linear between nodes,
+ * held at the grid's ends). The ONE source for every runtime — an
+ * artifact-less createModel() and the browser before its series load compute
+ * the identical rows.
+ * @param {number} bodyIndex 0 Me, 1 V, 2 E, 3 Ma, 4 J, 5 S (only 1–4 are embedded)
+ * @param {number} T
+ * @returns {[number, number]}
+ */
+function eccVectorOfDateEmbedded(bodyIndex, T) {
+  const key = ECC_BODY_KEY[bodyIndex];
+  if (!key) throw new RangeError(`sun-planetary-completion: no embedded e-vector series for body index ${bodyIndex}`);
+  const grid = SUN_COMPLETION_ARTIFACT.eccVectors;
+  const b = /** @type {Record<string, {zQ: number[], zP: number[]}>} */ (grid.bodies)[key];
+  const x = Math.min(Math.max((T * 100 - grid.t0Yr) / grid.stepYr, 0), b.zQ.length - 1);
+  const i = Math.min(Math.floor(x), b.zQ.length - 2), f = x - i;
+  return [b.zQ[i] + (b.zQ[i + 1] - b.zQ[i]) * f, b.zP[i] + (b.zP[i + 1] - b.zP[i]) * f];
+}
+
+/**
+ * The completion's CARRIER rates (deg per Julian century TT) — ONE home for the
+ * model wiring and the matched-pair gate. Planets: the model's own banked J2000
+ * sidereal mean motions (the embedded artifact; the former carriers were the
+ * planet records' rounded of-date periods minus p₀ — Venus 4.1″/yr and Jupiter
+ * 1.7″/yr off the run's own motion, 83° of the Venus–Earth argument at ±9000 yr).
+ * Earth: the framework sidereal year (the D6 ratio-only doctrine for the run's
+ * absolute Earth rate). Moon elongation: sidereal month vs sidereal year.
+ * @param {{ meanSiderealYearDays: number, moonSiderealMonthDays: number }} c
+ * @returns {{ planets: number[], moonElongation: number }}
+ */
+function computeCarrierRatesDegPerCy({ meanSiderealYearDays, moonSiderealMonthDays }) {
+  const L = SUN_COMPLETION_ARTIFACT.planetLamDotJ2000DegPerYr;
+  const degPerCyOf = (/** @type {number} */ cyclesPerDay) => 360 * 36525 * cyclesPerDay;
+  return {
+    planets: [L.mercury * 100, L.venus * 100, degPerCyOf(1 / meanSiderealYearDays), L.mars * 100, L.jupiter * 100, L.saturn * 100],
+    moonElongation: degPerCyOf(1 / moonSiderealMonthDays - 1 / meanSiderealYearDays),
+  };
+}
+
 /**
  * @param {{ embWobbleArcsec: number,
- *           carrierRatesDegPerCy: { planets: number[], moonElongation: number } }} opts
+ *           carrierRatesDegPerCy: { planets: number[], moonElongation: number },
+ *           eccVectorOfDate?: (bodyIndex: number, T: number) => [number, number],
+ *           moonElongationDegAtT?: (T: number) => number }} opts
  *   - embWobbleArcsec: the DERIVED Earth-around-EMB wobble amplitude
  *     a_M·μ/AU in arcsec (μ = 1/(1+M_E/M_M));
- *   - carrierRatesDegPerCy.planets: the six FRAMEWORK mean-longitude
- *     rates (deg/Julian-century TT), body order Me,V,E,Ma,J,S — one
- *     revolution per the model's own tropical period records;
+ *   - carrierRatesDegPerCy.planets: the six mean-longitude rates (deg/
+ *     Julian-century TT), body order Me,V,E,Ma,J,S — computeCarrierRatesDegPerCy
+ *     (the model's own banked J2000 sidereal mean motions; Earth the framework
+ *     sidereal year);
  *   - carrierRatesDegPerCy.moonElongation: the framework Moon
- *     mean-elongation rate (deg/cy TT) for the EMB-wobble carrier.
- *   All computed from live constants by the model wiring so the
- *   carrier↔table matched pair tracks the constants.
+ *     mean-elongation rate (deg/cy TT) — the EMB-wobble carrier's rate when
+ *     moonElongationDegAtT is absent;
+ *   - eccVectorOfDate: the e-vectors of date for the D'Alembert rows (default:
+ *     the embedded grid — the ONE source; an override is for instruments);
+ *   - moonElongationDegAtT: the Moon's mean elongation OF DATE (deg) at T —
+ *     the framework's own argument (ṅ included). The constant-rate carrier
+ *     D₀ + D₁·T ran ~29° off it at ±9000 yr: the lunar equation (6.44″) read
+ *     as the largest short-period line left against DE441 at −9000 (2.8″ at
+ *     29.53 d; the −9000 bin's scatter 4.13 → 3.62″ on the argument of date).
+ *   Injected by the model wiring so the carrier↔table matched pair tracks the
+ *   constants and the embedded artifact.
  * @returns {{ sunPlanetaryCompletionDeg: (T: number) => number }}
  */
-function createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy }) {
+function createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy, eccVectorOfDate = eccVectorOfDateEmbedded, moonElongationDegAtT = undefined }) {
   if (!Number.isFinite(embWobbleArcsec)) {
     throw new Error('createSunPlanetaryCompletion: embWobbleArcsec must be a finite number (derived a_M·μ/AU in arcsec)');
   }
@@ -270,7 +372,15 @@ function createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy })
       for (let i = 0; i < 6; i++) th += kl[i] * l[i] + kM[i] * M[i];
       table += cA * Math.cos(th) + sA * Math.sin(th);
     }
-    const D = (ARG_D0 + D1 * T) * D2R;
+    // the long-period rows — amplitude and phase on the e-vectors of date
+    for (const [kl, X, a, b] of LONG_PERIOD_TERMS) {
+      let th = 0;
+      for (let i = 0; i < 6; i++) th += kl[i] * l[i];
+      const [q, p] = eccVectorOfDate(X, T);
+      const e = Math.hypot(q, p), w = Math.atan2(p, q);
+      table += e * (a * Math.cos(th + w) + b * Math.sin(th + w));
+    }
+    const D = moonElongationDegAtT ? moonElongationDegAtT(T) * D2R : (ARG_D0 + D1 * T) * D2R;
     const arcsec = -table - embWobbleArcsec * Math.sin(D);
     return arcsec / 3600;
   }
@@ -284,19 +394,21 @@ function createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy })
 const PAIRED_SUN_HARMONICS_SHA256 = 'cbc189cea1c20292';   // eccentricity unification: Step-0 refit → N2/N3 re-derived (fidelity 0.616″, 70 terms)
 
 /** sha256/16 of JSON.stringify([...planets, moonElongation]) — the seven
- *  full-precision carrier rates (deg/cy TT) the TERMS table pairs with:
- *  the SIDEREAL carriers of plan 06 I2 (record rate − the model's J2000
- *  precession; Earth the framework sidereal year; the N3 literals kept,
- *  the two long-period rows derived on these rates). The model wiring
- *  recomputes the rates live from the planet records, so a planet-period /
- *  year / month input change moves the carriers automatically while the
- *  table stays frozen — a silent few-arcsec stale below the api gate's
- *  ≤12″ backstop. test:model recomputes this fingerprint from live
- *  constants (identical arithmetic to model.js) and fails on mismatch:
- *  re-run the extraction chain (tools/explore/n2-sun-framework-carriers.mjs
- *  → i2-sidereal-carrier-table.mjs for the short-period rows,
- *  i2-long-inequality.mjs for the long-period rows), re-embed TERMS, and
- *  update this value. */
-const PAIRED_CARRIER_RATES_SHA256 = '2d066e92bae955e4';
+ *  full-precision carrier rates (deg/cy TT) the tables pair with: since plan
+ *  06 I3 the model's own banked J2000 sidereal mean motions for the planets
+ *  (the embedded sun-completion artifact), Earth the framework sidereal year,
+ *  the Moon elongation from the sidereal month/year identity
+ *  (computeCarrierRatesDegPerCy; the N3 short-period literals kept, the
+ *  LONG_PERIOD_TERMS derived on these rates). The model wiring recomputes
+ *  the rates live, so a series re-bank, a year or a month input change moves
+ *  the carriers automatically while the tables stay frozen — a silent
+ *  few-arcsec stale below the api gate's ≤12″ backstop. test:model recomputes
+ *  this fingerprint from the live artifact + constants (identical arithmetic)
+ *  and fails on mismatch: re-run the extraction chain
+ *  (tools/explore/n2-sun-framework-carriers.mjs → i2-sidereal-carrier-table.mjs
+ *  for the short-period rows, i3-long-period-dalembert.mjs for the long-period
+ *  rows), re-embed the tables, and update this value. History: 2d066e92bae955e4
+ *  was the I2 record-based carrier set. */
+const PAIRED_CARRIER_RATES_SHA256 = '8e11456fcf8db81c';
 
-module.exports = { createSunPlanetaryCompletion, PAIRED_SUN_HARMONICS_SHA256, PAIRED_CARRIER_RATES_SHA256 };
+module.exports = { createSunPlanetaryCompletion, computeCarrierRatesDegPerCy, eccVectorOfDateEmbedded, PAIRED_SUN_HARMONICS_SHA256, PAIRED_CARRIER_RATES_SHA256 };

package/src/model.js

@@ -36,7 +36,7 @@ import { createMoonMonthChain } from './moon/month-chain.cjs';
 import { createChainCycleIntegrator } from './chain-cycles/index.cjs';
 import { createMoonArguments, jdToDecimalYear } from './moon/arguments.cjs';
 import { createMoonSeries } from './moon/series.cjs';
-import { createSunPlanetaryCompletion } from './eclipse/sun-planetary-completion.cjs';
+import { createSunPlanetaryCompletion, computeCarrierRatesDegPerCy } from './eclipse/sun-planetary-completion.cjs';
 import { createEclipseFinders } from './eclipse/finders.cjs';
 import { createBesselian } from './eclipse/besselian.cjs';
 import { driver2PeriodSecondsAtAge } from './planets/orbit-chain.cjs';
@@ -1276,20 +1276,21 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   // (the composed table reproduces the N3 one there: JPL 1900–2100
   // all-phase sd unchanged to 0.01″). The parity gate's fingerprint mirrors
   // this arithmetic (test/create-model-parity.test.mjs).
-  const degPerCyOf = /** @param {number} cyclesPerDay */ (cyclesPerDay) => 360 * 36525 * cyclesPerDay;
-  const carrierPrecessionDegPerCy = degPerCyOf(1 / meanSolarYearDays) - degPerCyOf(1 / meanSiderealYearDays);
-  const carrierRatesDegPerCy = {
-    planets: [
-      degPerCyOf(1 / C.planetOrbitalElements.mercury.solarYearInput) - carrierPrecessionDegPerCy,
-      degPerCyOf(1 / C.planetOrbitalElements.venus.solarYearInput) - carrierPrecessionDegPerCy,
-      degPerCyOf(1 / meanSiderealYearDays),
-      degPerCyOf(1 / C.planetOrbitalElements.mars.solarYearInput) - carrierPrecessionDegPerCy,
-      degPerCyOf(1 / C.planetOrbitalElements.jupiter.solarYearInput) - carrierPrecessionDegPerCy,
-      degPerCyOf(1 / C.planetOrbitalElements.saturn.solarYearInput) - carrierPrecessionDegPerCy,
-    ],
-    moonElongation: degPerCyOf(1 / moonSiderealMonthInput - 1 / meanSiderealYearDays),
-  };
-  const { sunPlanetaryCompletionDeg } = createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy });
+  // PLAN 06 I3 — the carriers are the model's own banked J2000 sidereal mean
+  // motions (the embedded sun-completion artifact; ONE home,
+  // computeCarrierRatesDegPerCy — the former record-based arithmetic sat
+  // 4.1″/yr off for Venus, 1.7″/yr for Jupiter: 83° of the Venus–Earth
+  // argument at ±9000 yr), Earth the framework sidereal year, the Moon
+  // elongation from the sidereal month/year identity. The long-period rows
+  // ride the embedded e-vectors of date, and the lunar equation rides the
+  // framework's own Moon elongation of date (ṅ included) — the constant-rate
+  // carrier ran ~29° off it at ±9000 yr.
+  const carrierRatesDegPerCy = computeCarrierRatesDegPerCy({ meanSiderealYearDays, moonSiderealMonthDays: moonSiderealMonthInput });
+  const { sunPlanetaryCompletionDeg } = createSunPlanetaryCompletion({
+    embWobbleArcsec,
+    carrierRatesDegPerCy,
+    moonElongationDegAtT: /** @param {number} T */ (T) => moonArgs.argsAt(j2000JD + T * julianCenturyDays).D,
+  });
   const besselian = createBesselian({
     moonFullAtDaysTT: /** @param {number} dDaysTT */ (dDaysTT) => {
       const ev = moonSeries.sceneEvalAt(dDaysTT);

package/types/eclipse/sun-completion-artifact.d.cts

@@ -0,0 +1,38 @@
+export const SUN_COMPLETION_ARTIFACT: Readonly<{
+    planetLamDotJ2000DegPerYr: {
+        mercury: number;
+        venus: number;
+        mars: number;
+        jupiter: number;
+        saturn: number;
+        uranus: number;
+        neptune: number;
+    };
+    eccVectors: {
+        t0Yr: number;
+        stepYr: number;
+        bodies: {
+            venus: {
+                zQ: number[];
+                zP: number[];
+            };
+            earth: {
+                zQ: number[];
+                zP: number[];
+            };
+            mars: {
+                zQ: number[];
+                zP: number[];
+            };
+            jupiter: {
+                zQ: number[];
+                zP: number[];
+            };
+        };
+    };
+    meta: {
+        dumpSha256: string;
+        source: string;
+    };
+}>;
+export const SUN_COMPLETION_ARTIFACT_HASH: "5ad17b275a003f6e";

package/types/eclipse/sun-planetary-completion.d.cts

@@ -1,43 +1,87 @@
 /**
  * @param {{ embWobbleArcsec: number,
- *           carrierRatesDegPerCy: { planets: number[], moonElongation: number } }} opts
+ *           carrierRatesDegPerCy: { planets: number[], moonElongation: number },
+ *           eccVectorOfDate?: (bodyIndex: number, T: number) => [number, number],
+ *           moonElongationDegAtT?: (T: number) => number }} opts
  *   - embWobbleArcsec: the DERIVED Earth-around-EMB wobble amplitude
  *     a_M·μ/AU in arcsec (μ = 1/(1+M_E/M_M));
- *   - carrierRatesDegPerCy.planets: the six FRAMEWORK mean-longitude
- *     rates (deg/Julian-century TT), body order Me,V,E,Ma,J,S — one
- *     revolution per the model's own tropical period records;
+ *   - carrierRatesDegPerCy.planets: the six mean-longitude rates (deg/
+ *     Julian-century TT), body order Me,V,E,Ma,J,S — computeCarrierRatesDegPerCy
+ *     (the model's own banked J2000 sidereal mean motions; Earth the framework
+ *     sidereal year);
  *   - carrierRatesDegPerCy.moonElongation: the framework Moon
- *     mean-elongation rate (deg/cy TT) for the EMB-wobble carrier.
- *   All computed from live constants by the model wiring so the
- *   carrier↔table matched pair tracks the constants.
+ *     mean-elongation rate (deg/cy TT) — the EMB-wobble carrier's rate when
+ *     moonElongationDegAtT is absent;
+ *   - eccVectorOfDate: the e-vectors of date for the D'Alembert rows (default:
+ *     the embedded grid — the ONE source; an override is for instruments);
+ *   - moonElongationDegAtT: the Moon's mean elongation OF DATE (deg) at T —
+ *     the framework's own argument (ṅ included). The constant-rate carrier
+ *     D₀ + D₁·T ran ~29° off it at ±9000 yr: the lunar equation (6.44″) read
+ *     as the largest short-period line left against DE441 at −9000 (2.8″ at
+ *     29.53 d; the −9000 bin's scatter 4.13 → 3.62″ on the argument of date).
+ *   Injected by the model wiring so the carrier↔table matched pair tracks the
+ *   constants and the embedded artifact.
  * @returns {{ sunPlanetaryCompletionDeg: (T: number) => number }}
  */
-export function createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy }: {
+export function createSunPlanetaryCompletion({ embWobbleArcsec, carrierRatesDegPerCy, eccVectorOfDate, moonElongationDegAtT }: {
     embWobbleArcsec: number;
     carrierRatesDegPerCy: {
         planets: number[];
         moonElongation: number;
     };
+    eccVectorOfDate?: (bodyIndex: number, T: number) => [number, number];
+    moonElongationDegAtT?: (T: number) => number;
 }): {
     sunPlanetaryCompletionDeg: (T: number) => number;
 };
+/**
+ * The completion's CARRIER rates (deg per Julian century TT) — ONE home for the
+ * model wiring and the matched-pair gate. Planets: the model's own banked J2000
+ * sidereal mean motions (the embedded artifact; the former carriers were the
+ * planet records' rounded of-date periods minus p₀ — Venus 4.1″/yr and Jupiter
+ * 1.7″/yr off the run's own motion, 83° of the Venus–Earth argument at ±9000 yr).
+ * Earth: the framework sidereal year (the D6 ratio-only doctrine for the run's
+ * absolute Earth rate). Moon elongation: sidereal month vs sidereal year.
+ * @param {{ meanSiderealYearDays: number, moonSiderealMonthDays: number }} c
+ * @returns {{ planets: number[], moonElongation: number }}
+ */
+export function computeCarrierRatesDegPerCy({ meanSiderealYearDays, moonSiderealMonthDays }: {
+    meanSiderealYearDays: number;
+    moonSiderealMonthDays: number;
+}): {
+    planets: number[];
+    moonElongation: number;
+};
+/**
+ * The e-vector (e·cos ϖ, e·sin ϖ; ecliptic J2000) of a body at T Julian
+ * centuries TT from J2000, from the embedded 1-kyr grid (linear between nodes,
+ * held at the grid's ends). The ONE source for every runtime — an
+ * artifact-less createModel() and the browser before its series load compute
+ * the identical rows.
+ * @param {number} bodyIndex 0 Me, 1 V, 2 E, 3 Ma, 4 J, 5 S (only 1–4 are embedded)
+ * @param {number} T
+ * @returns {[number, number]}
+ */
+export function eccVectorOfDateEmbedded(bodyIndex: number, T: number): [number, number];
 /** sha256/16 of JSON.stringify(FITTED_COEFFICIENTS.SUN_LONGITUDE_HARMONICS)
  *  at derivation time — the matched-pair fingerprint asserted by test:model.
  *  Unchanged from v1/v2: SUN_HARMONICS did not move in the D2 or N3
  *  landings. */
 export const PAIRED_SUN_HARMONICS_SHA256: "cbc189cea1c20292";
 /** sha256/16 of JSON.stringify([...planets, moonElongation]) — the seven
- *  full-precision carrier rates (deg/cy TT) the TERMS table pairs with:
- *  the SIDEREAL carriers of plan 06 I2 (record rate − the model's J2000
- *  precession; Earth the framework sidereal year; the N3 literals kept,
- *  the two long-period rows derived on these rates). The model wiring
- *  recomputes the rates live from the planet records, so a planet-period /
- *  year / month input change moves the carriers automatically while the
- *  table stays frozen — a silent few-arcsec stale below the api gate's
- *  ≤12″ backstop. test:model recomputes this fingerprint from live
- *  constants (identical arithmetic to model.js) and fails on mismatch:
- *  re-run the extraction chain (tools/explore/n2-sun-framework-carriers.mjs
- *  → i2-sidereal-carrier-table.mjs for the short-period rows,
- *  i2-long-inequality.mjs for the long-period rows), re-embed TERMS, and
- *  update this value. */
-export const PAIRED_CARRIER_RATES_SHA256: "2d066e92bae955e4";
+ *  full-precision carrier rates (deg/cy TT) the tables pair with: since plan
+ *  06 I3 the model's own banked J2000 sidereal mean motions for the planets
+ *  (the embedded sun-completion artifact), Earth the framework sidereal year,
+ *  the Moon elongation from the sidereal month/year identity
+ *  (computeCarrierRatesDegPerCy; the N3 short-period literals kept, the
+ *  LONG_PERIOD_TERMS derived on these rates). The model wiring recomputes
+ *  the rates live, so a series re-bank, a year or a month input change moves
+ *  the carriers automatically while the tables stay frozen — a silent
+ *  few-arcsec stale below the api gate's ≤12″ backstop. test:model recomputes
+ *  this fingerprint from the live artifact + constants (identical arithmetic)
+ *  and fails on mismatch: re-run the extraction chain
+ *  (tools/explore/n2-sun-framework-carriers.mjs → i2-sidereal-carrier-table.mjs
+ *  for the short-period rows, i3-long-period-dalembert.mjs for the long-period
+ *  rows), re-embed the tables, and update this value. History: 2d066e92bae955e4
+ *  was the I2 record-based carrier set. */
+export const PAIRED_CARRIER_RATES_SHA256: "8e11456fcf8db81c";