universal-physics-tensor 0.4.5

package/dist/bridges/equations/be-28-onsager-entropy-production.d.ts

@@ -0,0 +1,175 @@
+/**
+ * Bridge Equation 28 — Maximum Entropy Production Principle
+ * (Onsager linear-response entropy-production scalar reformulation).
+ *
+ * **Original (un-encodable) form.** The Dewar 2003/2005 variational
+ * formulation of MEPP:
+ *
+ *   δ ∫ [dS/dt - λ Σᵢ Jᵢ Xᵢ - μ(∇·v)] dt = 0
+ *
+ * is not AST-encodable. It requires three grammar extensions: (1)
+ * variational δ-operator over functionals; (2) Lagrange multipliers
+ * λ, μ for constraint enforcement; (3) discrete-index sum primitive
+ * over species. MEPP itself is also a **contested principle** —
+ * Grinstein-Linsker 2007 *J. Phys. A* 40:9717 produced counterexamples
+ * to Dewar's information-theoretic derivation, and the principle
+ * conflicts with Prigogine's *minimum* entropy production principle
+ * for near-equilibrium systems.
+ *
+ * **Encoded reformulation (Onsager linear-response entropy production).**
+ * Wave Z-G (2026-05-11): per user-confirmed design choice after
+ * surfacing the relabeling concern, replace the MEPP variational
+ * formulation with the canonical **Onsager linear-response
+ * entropy-production scalar**:
+ *
+ *   **σ = Σᵢ Jᵢ Xᵢ**     (typed scalar with dim [entropy/time])
+ *
+ * where Jᵢ are generalized fluxes (heat, particle, charge) and Xᵢ
+ * are their conjugate generalized forces (∇T/T², chemical-potential
+ * gradients, electric field, etc.). The discrete index sum is
+ * collapsed into a single typed-stub symbol `force_flux_product` —
+ * the AST has no discrete-index sum primitive, so the sum-collapsed
+ * stub absorbs the multi-species content. Dim `[entropy / time]`
+ * = `[L² M T⁻³ Θ⁻¹]` = `[W/K]` (power per kelvin).
+ *
+ * **⚠ CRITICAL WARNING — Definiendum vs. principle (honest-claude
+ * scope, REQUIRED reading).** This reformulation **does NOT encode
+ * the Maximum Entropy Production Principle (MEPP).** The original
+ * MEPP is a contested **variational principle** asserting that
+ * non-equilibrium systems organize to *maximize* entropy production
+ * subject to constraints. Due to AST grammar limitations (no
+ * variational δ, no Lagrange multiplier, no discrete-index sum),
+ * that principle cannot be represented. This equation instead
+ * encodes the foundational Onsager / Prigogine **definition** of the
+ * entropy production rate (σ) as a sum of flux-force products. This
+ * bridge is retained under the BE-28 label for historical continuity,
+ * but it represents the *definiendum* of MEPP (the quantity MEPP
+ * makes a claim about), NOT the maximization conjecture itself.
+ *
+ * Wording owed to the Gemini Pro Wave-Z cross-validation
+ * consultation (2026-05-11), which independently confirmed the
+ * "DEFENSIBLE-WITH-CAVEATS" verdict and recommended this stronger
+ * warning prefix. Both OpenAI o3 and Gemini Pro agreed: the
+ * relabeling is acceptable *only* because the user explicitly
+ * accepted the trade-off via `AskUserQuestion` and the warning is
+ * prominent. See Wave Z-G notes in `src/bridges/index.ts` for the
+ * full discussion.
+ *
+ * The reformulation is closer in spirit to a **renaming** of BE-28
+ * (from MEPP → Onsager entropy production) than to the BE-25 /
+ * BE-16 / BE-37 reformulations, which preserved their bridge labels
+ * (consciousness ↔ information; information ↔ thermodynamics;
+ * modified light propagation). MEPP's bridge label "Why nature
+ * chooses specific NESS" is NOT preserved by Onsager. Wave Z-G is
+ * a project-internal decision to "close" BE-28 at the cost of the
+ * MEPP semantic content.
+ *
+ * **AST encoding pattern.** Single typed-stub for the force-flux
+ * sum. Dim: `[L² M T⁻³ Θ⁻¹]` (entropy per time, equivalently W/K).
+ *
+ * Bracket-checks (numerical evaluator):
+ *   - σ = 1 W/K: entropy production at the rate of 1 joule of heat
+ *     per second per kelvin. Typical of a small heat engine or
+ *     dissipative process at room temperature.
+ *   - σ ≥ 0: the Second Law in differential form (entropy
+ *     production is non-negative for any irreversible process).
+ *     The numerical evaluator enforces this with a RangeError.
+ *
+ * References:
+ *   - Onsager 1931 *Phys. Rev.* 37:405 and 38:2265 (reciprocal
+ *     relations — the foundational paper).
+ *   - de Groot & Mazur 1962 *Non-Equilibrium Thermodynamics*
+ *     (canonical textbook).
+ *   - Dewar 2003 *J. Phys. A* 36:631 (MEPP information-theoretic
+ *     derivation — the original-BE-28 source, now dropped).
+ *   - Dewar 2005 *J. Phys. A* 38:L371 (MEPP follow-up).
+ *   - Grinstein & Linsker 2007 *J. Phys. A* 40:9717 (MEPP
+ *     rebuttal — counterexamples to Dewar's derivation; one of the
+ *     reasons MEPP-as-variational-principle is contested).
+ *   - Prigogine 1947 *Étude thermodynamique des phénomènes
+ *     irréversibles* (minimum-entropy-production for near-
+ *     equilibrium — the contrasting principle that conflicts with
+ *     MEPP in the linear regime).
+ *
+ * Status: speculative.
+ *
+ * Honest-claude scope notes (extended):
+ *   - Onsager linear-response is **already implicitly used** by
+ *     BE-21 (KSS η/s bound), BE-23 (SYK Planckian resistivity), and
+ *     BE-29 (Jarzynski free-energy equality). Encoding BE-28 as
+ *     "Onsager force-flux entropy production" creates some content
+ *     overlap with those bridges; the distinguishing feature is
+ *     that BE-28 encodes the σ scalar itself, not a derived
+ *     quantity (transport coefficient or fluctuation theorem).
+ *   - The MEPP variational principle is NOT captured. Anyone reading
+ *     the encoded scalar should understand that the encoding answers
+ *     "what is the entropy production rate?" but NOT "why does
+ *     nature select this rate?" — the latter is the actual MEPP
+ *     content and remains an open research question (Dewar 2003,
+ *     Grinstein-Linsker 2007 rebuttal, Endres 2017 follow-ups).
+ *   - The discrete index sum Σᵢ over species is collapsed into a
+ *     single typed-stub `force_flux_product`. A future grammar
+ *     extension with explicit discrete-index sum primitive could
+ *     unfold this stub into an explicit `J_heat·X_heat +
+ *     J_particle·X_particle + J_charge·X_charge` expression.
+ *   - σ ≥ 0 (Second Law) is enforced by the numerical evaluator but
+ *     not by the AST type system — the dimension allows negative
+ *     values; the non-negativity is a semantic property.
+ *   - The dim [entropy/time] = [L² M T⁻³ Θ⁻¹] = [W/K] is the
+ *     **total** entropy production rate (extensive quantity). The
+ *     local entropy production density σ_local with dim
+ *     [entropy/(volume · time)] = [L⁻¹ M T⁻³ Θ⁻¹] is a different
+ *     scalar; the encoding choice is the **total** rate, consistent
+ *     with the integral form ∫ J·X d³x = σ_total over the system.
+ *
+ * @see docs/specification/Part-II.md ("Bridge Equation 28")
+ * @see src/bridges/index.ts BRIDGE_EQUATIONS.find(e => e.id === 28)
+ * @module bridges/equations/be-28-onsager-entropy-production
+ */
+import type { ExprNode, DimensionValidationReport } from '../../dimensional/validator.js';
+/**
+ * Symbol for the index-collapsed force-flux product Σᵢ Jᵢ Xᵢ.
+ *
+ * The AST has no discrete-index sum primitive, so the multi-species
+ * sum is encoded as a single typed-stub with the dim of σ itself.
+ * Same idiom as BE-17 `T_torsion_squared` (typed-stub for a tensor
+ * contraction the AST cannot expand) and BE-46 `exp_factor` (typed-
+ * stub for an unexpressible transcendental subexpression).
+ *
+ * Dim `[entropy/time]` = `[W/K]`.
+ */
+export declare const BE28_FORCE_FLUX_PRODUCT: ExprNode;
+/**
+ * RHS of `σ = Σᵢ Jᵢ Xᵢ`:
+ *
+ *   force_flux_product
+ *
+ * Encoded as a single typed-stub (no internal expansion in the AST).
+ * Dim `[entropy/time]`.
+ */
+export declare const BE28_ENTROPY_PRODUCTION_RHS: ExprNode;
+/** LHS: σ (total entropy production rate). Dim `[entropy/time]`. */
+export declare const BE28_ENTROPY_PRODUCTION_LHS: ExprNode;
+export interface OnsagerEntropyInputs {
+    /** Sum of force-flux products Σᵢ Jᵢ Xᵢ in SI units [W/K]. Must be finite and ≥ 0 (Second Law). */
+    force_flux_product_W_per_K: number;
+}
+/**
+ * Evaluate the Onsager linear-response entropy production rate:
+ *
+ *   σ = Σᵢ Jᵢ Xᵢ
+ *
+ * The caller supplies the already-summed force-flux product as a
+ * scalar; the AST cannot expand the discrete-index sum, so the
+ * evaluator takes the sum as input.
+ *
+ * @returns σ in W/K (joule per second per kelvin). Must be ≥ 0
+ *   by the Second Law.
+ */
+export declare function evaluateOnsagerEntropyProduction(input: OnsagerEntropyInputs): number;
+/**
+ * Run the AST through the dimensional analyzer; LHS and RHS should
+ * both be `[entropy/time]`.
+ */
+export declare function validateBE28Dimensions(): DimensionValidationReport;
+//# sourceMappingURL=be-28-onsager-entropy-production.d.ts.map

package/dist/bridges/equations/be-28-onsager-entropy-production.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"be-28-onsager-entropy-production.d.ts","sourceRoot":"","sources":["../../../src/bridges/equations/be-28-onsager-entropy-production.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+HG;AAEH,OAAO,KAAK,EAAE,QAAQ,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAqB1F;;;;;;;;;;GAUG;AACH,eAAO,MAAM,uBAAuB,EAAE,QAGrC,CAAC;AAEF;;;;;;;GAOG;AACH,eAAO,MAAM,2BAA2B,EAAE,QAAkC,CAAC;AAE7E,oEAAoE;AACpE,eAAO,MAAM,2BAA2B,EAAE,QAAqC,CAAC;AAIhF,MAAM,WAAW,oBAAoB;IACnC,kGAAkG;IAClG,0BAA0B,EAAE,MAAM,CAAC;CACpC;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,gCAAgC,CAAC,KAAK,EAAE,oBAAoB,GAAG,MAAM,CAepF;AAID;;;GAGG;AACH,wBAAgB,sBAAsB,IAAI,yBAAyB,CASlE"}

package/dist/bridges/equations/be-28-onsager-entropy-production.js

@@ -0,0 +1,203 @@
+/**
+ * Bridge Equation 28 — Maximum Entropy Production Principle
+ * (Onsager linear-response entropy-production scalar reformulation).
+ *
+ * **Original (un-encodable) form.** The Dewar 2003/2005 variational
+ * formulation of MEPP:
+ *
+ *   δ ∫ [dS/dt - λ Σᵢ Jᵢ Xᵢ - μ(∇·v)] dt = 0
+ *
+ * is not AST-encodable. It requires three grammar extensions: (1)
+ * variational δ-operator over functionals; (2) Lagrange multipliers
+ * λ, μ for constraint enforcement; (3) discrete-index sum primitive
+ * over species. MEPP itself is also a **contested principle** —
+ * Grinstein-Linsker 2007 *J. Phys. A* 40:9717 produced counterexamples
+ * to Dewar's information-theoretic derivation, and the principle
+ * conflicts with Prigogine's *minimum* entropy production principle
+ * for near-equilibrium systems.
+ *
+ * **Encoded reformulation (Onsager linear-response entropy production).**
+ * Wave Z-G (2026-05-11): per user-confirmed design choice after
+ * surfacing the relabeling concern, replace the MEPP variational
+ * formulation with the canonical **Onsager linear-response
+ * entropy-production scalar**:
+ *
+ *   **σ = Σᵢ Jᵢ Xᵢ**     (typed scalar with dim [entropy/time])
+ *
+ * where Jᵢ are generalized fluxes (heat, particle, charge) and Xᵢ
+ * are their conjugate generalized forces (∇T/T², chemical-potential
+ * gradients, electric field, etc.). The discrete index sum is
+ * collapsed into a single typed-stub symbol `force_flux_product` —
+ * the AST has no discrete-index sum primitive, so the sum-collapsed
+ * stub absorbs the multi-species content. Dim `[entropy / time]`
+ * = `[L² M T⁻³ Θ⁻¹]` = `[W/K]` (power per kelvin).
+ *
+ * **⚠ CRITICAL WARNING — Definiendum vs. principle (honest-claude
+ * scope, REQUIRED reading).** This reformulation **does NOT encode
+ * the Maximum Entropy Production Principle (MEPP).** The original
+ * MEPP is a contested **variational principle** asserting that
+ * non-equilibrium systems organize to *maximize* entropy production
+ * subject to constraints. Due to AST grammar limitations (no
+ * variational δ, no Lagrange multiplier, no discrete-index sum),
+ * that principle cannot be represented. This equation instead
+ * encodes the foundational Onsager / Prigogine **definition** of the
+ * entropy production rate (σ) as a sum of flux-force products. This
+ * bridge is retained under the BE-28 label for historical continuity,
+ * but it represents the *definiendum* of MEPP (the quantity MEPP
+ * makes a claim about), NOT the maximization conjecture itself.
+ *
+ * Wording owed to the Gemini Pro Wave-Z cross-validation
+ * consultation (2026-05-11), which independently confirmed the
+ * "DEFENSIBLE-WITH-CAVEATS" verdict and recommended this stronger
+ * warning prefix. Both OpenAI o3 and Gemini Pro agreed: the
+ * relabeling is acceptable *only* because the user explicitly
+ * accepted the trade-off via `AskUserQuestion` and the warning is
+ * prominent. See Wave Z-G notes in `src/bridges/index.ts` for the
+ * full discussion.
+ *
+ * The reformulation is closer in spirit to a **renaming** of BE-28
+ * (from MEPP → Onsager entropy production) than to the BE-25 /
+ * BE-16 / BE-37 reformulations, which preserved their bridge labels
+ * (consciousness ↔ information; information ↔ thermodynamics;
+ * modified light propagation). MEPP's bridge label "Why nature
+ * chooses specific NESS" is NOT preserved by Onsager. Wave Z-G is
+ * a project-internal decision to "close" BE-28 at the cost of the
+ * MEPP semantic content.
+ *
+ * **AST encoding pattern.** Single typed-stub for the force-flux
+ * sum. Dim: `[L² M T⁻³ Θ⁻¹]` (entropy per time, equivalently W/K).
+ *
+ * Bracket-checks (numerical evaluator):
+ *   - σ = 1 W/K: entropy production at the rate of 1 joule of heat
+ *     per second per kelvin. Typical of a small heat engine or
+ *     dissipative process at room temperature.
+ *   - σ ≥ 0: the Second Law in differential form (entropy
+ *     production is non-negative for any irreversible process).
+ *     The numerical evaluator enforces this with a RangeError.
+ *
+ * References:
+ *   - Onsager 1931 *Phys. Rev.* 37:405 and 38:2265 (reciprocal
+ *     relations — the foundational paper).
+ *   - de Groot & Mazur 1962 *Non-Equilibrium Thermodynamics*
+ *     (canonical textbook).
+ *   - Dewar 2003 *J. Phys. A* 36:631 (MEPP information-theoretic
+ *     derivation — the original-BE-28 source, now dropped).
+ *   - Dewar 2005 *J. Phys. A* 38:L371 (MEPP follow-up).
+ *   - Grinstein & Linsker 2007 *J. Phys. A* 40:9717 (MEPP
+ *     rebuttal — counterexamples to Dewar's derivation; one of the
+ *     reasons MEPP-as-variational-principle is contested).
+ *   - Prigogine 1947 *Étude thermodynamique des phénomènes
+ *     irréversibles* (minimum-entropy-production for near-
+ *     equilibrium — the contrasting principle that conflicts with
+ *     MEPP in the linear regime).
+ *
+ * Status: speculative.
+ *
+ * Honest-claude scope notes (extended):
+ *   - Onsager linear-response is **already implicitly used** by
+ *     BE-21 (KSS η/s bound), BE-23 (SYK Planckian resistivity), and
+ *     BE-29 (Jarzynski free-energy equality). Encoding BE-28 as
+ *     "Onsager force-flux entropy production" creates some content
+ *     overlap with those bridges; the distinguishing feature is
+ *     that BE-28 encodes the σ scalar itself, not a derived
+ *     quantity (transport coefficient or fluctuation theorem).
+ *   - The MEPP variational principle is NOT captured. Anyone reading
+ *     the encoded scalar should understand that the encoding answers
+ *     "what is the entropy production rate?" but NOT "why does
+ *     nature select this rate?" — the latter is the actual MEPP
+ *     content and remains an open research question (Dewar 2003,
+ *     Grinstein-Linsker 2007 rebuttal, Endres 2017 follow-ups).
+ *   - The discrete index sum Σᵢ over species is collapsed into a
+ *     single typed-stub `force_flux_product`. A future grammar
+ *     extension with explicit discrete-index sum primitive could
+ *     unfold this stub into an explicit `J_heat·X_heat +
+ *     J_particle·X_particle + J_charge·X_charge` expression.
+ *   - σ ≥ 0 (Second Law) is enforced by the numerical evaluator but
+ *     not by the AST type system — the dimension allows negative
+ *     values; the non-negativity is a semantic property.
+ *   - The dim [entropy/time] = [L² M T⁻³ Θ⁻¹] = [W/K] is the
+ *     **total** entropy production rate (extensive quantity). The
+ *     local entropy production density σ_local with dim
+ *     [entropy/(volume · time)] = [L⁻¹ M T⁻³ Θ⁻¹] is a different
+ *     scalar; the encoding choice is the **total** rate, consistent
+ *     with the integral form ∫ J·X d³x = σ_total over the system.
+ *
+ * @see docs/specification/Part-II.md ("Bridge Equation 28")
+ * @see src/bridges/index.ts BRIDGE_EQUATIONS.find(e => e.id === 28)
+ * @module bridges/equations/be-28-onsager-entropy-production
+ */
+import { validate, validateEquation } from '../../dimensional/validator.js';
+const sym = (name, dim) => ({ kind: 'symbol', name, dim });
+// --- Local dim constants ---
+/**
+ * [L² M T⁻³ Θ⁻¹] = [entropy/time] = [W/K] — total entropy production
+ * rate. ENTROPY = {L:2, M:1, T:-2, Θ:-1}; dividing by TIME {T:1}
+ * gives this dim.
+ */
+const ENTROPY_RATE = {
+    L: 2, M: 1, T: -3, I: 0, Theta: -1, N: 0, J: 0,
+};
+// --- Symbolic AST ---
+/**
+ * Symbol for the index-collapsed force-flux product Σᵢ Jᵢ Xᵢ.
+ *
+ * The AST has no discrete-index sum primitive, so the multi-species
+ * sum is encoded as a single typed-stub with the dim of σ itself.
+ * Same idiom as BE-17 `T_torsion_squared` (typed-stub for a tensor
+ * contraction the AST cannot expand) and BE-46 `exp_factor` (typed-
+ * stub for an unexpressible transcendental subexpression).
+ *
+ * Dim `[entropy/time]` = `[W/K]`.
+ */
+export const BE28_FORCE_FLUX_PRODUCT = sym('force_flux_product', ENTROPY_RATE);
+/**
+ * RHS of `σ = Σᵢ Jᵢ Xᵢ`:
+ *
+ *   force_flux_product
+ *
+ * Encoded as a single typed-stub (no internal expansion in the AST).
+ * Dim `[entropy/time]`.
+ */
+export const BE28_ENTROPY_PRODUCTION_RHS = BE28_FORCE_FLUX_PRODUCT;
+/** LHS: σ (total entropy production rate). Dim `[entropy/time]`. */
+export const BE28_ENTROPY_PRODUCTION_LHS = sym('sigma', ENTROPY_RATE);
+/**
+ * Evaluate the Onsager linear-response entropy production rate:
+ *
+ *   σ = Σᵢ Jᵢ Xᵢ
+ *
+ * The caller supplies the already-summed force-flux product as a
+ * scalar; the AST cannot expand the discrete-index sum, so the
+ * evaluator takes the sum as input.
+ *
+ * @returns σ in W/K (joule per second per kelvin). Must be ≥ 0
+ *   by the Second Law.
+ */
+export function evaluateOnsagerEntropyProduction(input) {
+    const { force_flux_product_W_per_K } = input;
+    if (!Number.isFinite(force_flux_product_W_per_K)) {
+        throw new RangeError(`evaluateOnsagerEntropyProduction: force_flux_product_W_per_K must be finite, got ${force_flux_product_W_per_K}`);
+    }
+    if (force_flux_product_W_per_K < 0) {
+        throw new RangeError(`evaluateOnsagerEntropyProduction: σ ≥ 0 by the Second Law; got ${force_flux_product_W_per_K}. ` +
+            `If the input represents a fluctuation that can be negative on short timescales, ` +
+            `it is not an entropy production rate — check the sign convention or the choice of conjugate forces.`);
+    }
+    return force_flux_product_W_per_K;
+}
+// --- Self-validation ---
+/**
+ * Run the AST through the dimensional analyzer; LHS and RHS should
+ * both be `[entropy/time]`.
+ */
+export function validateBE28Dimensions() {
+    const eq = validateEquation(BE28_ENTROPY_PRODUCTION_LHS, BE28_ENTROPY_PRODUCTION_RHS);
+    const lhs = validate(BE28_ENTROPY_PRODUCTION_LHS);
+    const rhs = validate(BE28_ENTROPY_PRODUCTION_RHS);
+    return {
+        ok: eq.ok,
+        lhsDim: lhs.inferredDimension,
+        rhsDim: rhs.inferredDimension,
+    };
+}
+//# sourceMappingURL=be-28-onsager-entropy-production.js.map

package/dist/bridges/equations/be-28-onsager-entropy-production.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"be-28-onsager-entropy-production.js","sourceRoot":"","sources":["../../../src/bridges/equations/be-28-onsager-entropy-production.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+HG;AAGH,OAAO,EAAE,QAAQ,EAAE,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AAK5E,MAAM,GAAG,GAAG,CAAC,IAAY,EAAE,GAAc,EAAY,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,QAAQ,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC;AAExF,8BAA8B;AAE9B;;;;GAIG;AACH,MAAM,YAAY,GAAc;IAC9B,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;CAC/C,CAAC;AAEF,uBAAuB;AAEvB;;;;;;;;;;GAUG;AACH,MAAM,CAAC,MAAM,uBAAuB,GAAa,GAAG,CAClD,oBAAoB,EACpB,YAAY,CACb,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,2BAA2B,GAAa,uBAAuB,CAAC;AAE7E,oEAAoE;AACpE,MAAM,CAAC,MAAM,2BAA2B,GAAa,GAAG,CAAC,OAAO,EAAE,YAAY,CAAC,CAAC;AAShF;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,gCAAgC,CAAC,KAA2B;IAC1E,MAAM,EAAE,0BAA0B,EAAE,GAAG,KAAK,CAAC;IAC7C,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,0BAA0B,CAAC,EAAE,CAAC;QACjD,MAAM,IAAI,UAAU,CAClB,oFAAoF,0BAA0B,EAAE,CACjH,CAAC;IACJ,CAAC;IACD,IAAI,0BAA0B,GAAG,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,UAAU,CAClB,kEAAkE,0BAA0B,IAAI;YAC9F,kFAAkF;YAClF,qGAAqG,CACxG,CAAC;IACJ,CAAC;IACD,OAAO,0BAA0B,CAAC;AACpC,CAAC;AAED,0BAA0B;AAE1B;;;GAGG;AACH,MAAM,UAAU,sBAAsB;IACpC,MAAM,EAAE,GAAG,gBAAgB,CAAC,2BAA2B,EAAE,2BAA2B,CAAC,CAAC;IACtF,MAAM,GAAG,GAAG,QAAQ,CAAC,2BAA2B,CAAC,CAAC;IAClD,MAAM,GAAG,GAAG,QAAQ,CAAC,2BAA2B,CAAC,CAAC;IAClD,OAAO;QACL,EAAE,EAAE,EAAE,CAAC,EAAE;QACT,MAAM,EAAE,GAAG,CAAC,iBAAiB;QAC7B,MAAM,EAAE,GAAG,CAAC,iBAAiB;KAC9B,CAAC;AACJ,CAAC"}

package/dist/bridges/equations/be-29-jarzynski.d.ts

@@ -0,0 +1,86 @@
+/**
+ * Bridge Equation 29 — Jarzynski free-energy difference (canonical
+ * 1997 equality, post Wave Y reformulation).
+ *
+ *   ΔF = −k_B T ln⟨exp(−W / (k_B T))⟩
+ *
+ * Reference: Jarzynski 1997 *Phys. Rev. Lett.* 78:2690 (canonical
+ * non-equilibrium free-energy equality; exact identity for any work
+ * protocol connecting two equilibrium states at temperature T).
+ *
+ * Status: speculative.
+ *
+ * Honest-claude scope notes:
+ *   - The Jarzynski equality itself is canonical; the *bridge framing*
+ *     — applying it to gravitational work in curved spacetime, the
+ *     original BE-29 framing — is the speculative element. The
+ *     gravitational-correction term `(β/(2c⁴)) ∫T^μν δg_μν √(-g) d⁴x`
+ *     is dropped here as the AST-unencodable extension.
+ *   - The AST has no `exp` primitive nor `ln` primitive; the average
+ *     `⟨exp(-βW)⟩` is encoded as a single dimensionless symbol stub
+ *     `'avg_exp_minus_betaW'`, with the exp argument β·W = W/(k_B T)
+ *     exposed as a separate ExprNode `BE29_BETAW_ARG` (verifiable
+ *     dimensionless via the dimensional analyzer).
+ *   - Both sides have dim [energy]: ΔF [J]; k_B T [J]; ⟨exp(-βW)⟩
+ *     dimensionless; ln(dimensionless) dimensionless; product [J]. ✓
+ *
+ * Numerical bracket: at quasi-static (reversible) work, W = W_rev and
+ * ⟨exp(-βW)⟩ = exp(-βW_rev), so ΔF = -k_BT · (-βW_rev) = W_rev. The
+ * Jarzynski equality reduces to the equilibrium identity in this limit.
+ *
+ * @see docs/specification/Part-II.md ("Bridge Equation 29")
+ * @see src/bridges/index.ts BRIDGE_EQUATIONS.find(e => e.id === 29)
+ * @module bridges/equations/be-29-jarzynski
+ */
+import type { ExprNode, DimensionValidationReport } from '../../dimensional/validator.js';
+/**
+ * Lemma AST: the exp-argument β·W = W/(k_B T), exposed for verification
+ * that the exponent is dimensionless.
+ *
+ * [W]/[k_B T] = [J]/[J] = [1] ✓
+ */
+export declare const BE29_BETAW_ARG: ExprNode;
+/**
+ * RHS of the Jarzynski equality: `−k_B T · ln⟨exp(-βW)⟩`.
+ *
+ * The AST has no `ln` primitive; we encode the relation as
+ *   `ΔF = -k_B T · L`
+ * where `L` is a dimensionless symbol stub representing
+ * `ln⟨exp(-W/(k_B T))⟩`. The minus sign is captured by encoding the
+ * coefficient as `-1` (a dimensionless symbol stub named `-1`).
+ *
+ * Both terms multiplicatively combine to dim [energy] = [k_B T].
+ */
+export declare const BE29_JARZYNSKI_RHS: ExprNode;
+/** LHS: ΔF has dimension [energy]. */
+export declare const BE29_JARZYNSKI_LHS: ExprNode;
+export interface JarzynskiInputs {
+    /**
+     * Sample of work values W_i (joules) from N independent realizations
+     * of the non-equilibrium protocol (forward direction).
+     */
+    W_samples_J: ReadonlyArray<number>;
+    /** Temperature T (K). Must be > 0 and finite. */
+    T_K: number;
+}
+/**
+ * Evaluate the Jarzynski free-energy difference from a sample of work
+ * values:
+ *
+ *   ΔF = −k_B T ln⟨exp(−W / (k_B T))⟩
+ *
+ * The ensemble average is approximated by the sample mean over the
+ * `W_samples_J` array. For a reliable estimate the sample size must be
+ * large enough to capture the rare-event tail; the Jarzynski estimator
+ * has well-known finite-sample bias for low-temperature / high-work-
+ * variance protocols (Park-Schulten 2003 *J. Chem. Phys.* 119:5946).
+ *
+ * @returns Free-energy difference ΔF in joules.
+ */
+export declare function evaluateJarzynski(input: JarzynskiInputs): number;
+/**
+ * Run the AST through the dimensional analyzer; both sides should be
+ * [energy].
+ */
+export declare function validateBE29Dimensions(): DimensionValidationReport;
+//# sourceMappingURL=be-29-jarzynski.d.ts.map

package/dist/bridges/equations/be-29-jarzynski.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"be-29-jarzynski.d.ts","sourceRoot":"","sources":["../../../src/bridges/equations/be-29-jarzynski.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiCG;AAEH,OAAO,KAAK,EAAE,QAAQ,EAAE,yBAAyB,EAAE,MAAM,gCAAgC,CAAC;AAa1F;;;;;GAKG;AACH,eAAO,MAAM,cAAc,EAAE,QAY5B,CAAC;AAEF;;;;;;;;;;GAUG;AACH,eAAO,MAAM,kBAAkB,EAAE,QAQhC,CAAC;AAEF,sCAAsC;AACtC,eAAO,MAAM,kBAAkB,EAAE,QAAiC,CAAC;AAInE,MAAM,WAAW,eAAe;IAC9B;;;OAGG;IACH,WAAW,EAAE,aAAa,CAAC,MAAM,CAAC,CAAC;IACnC,iDAAiD;IACjD,GAAG,EAAE,MAAM,CAAC;CACb;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,iBAAiB,CAAC,KAAK,EAAE,eAAe,GAAG,MAAM,CA2BhE;AAID;;;GAGG;AACH,wBAAgB,sBAAsB,IAAI,yBAAyB,CASlE"}

package/dist/bridges/equations/be-29-jarzynski.js

@@ -0,0 +1,132 @@
+/**
+ * Bridge Equation 29 — Jarzynski free-energy difference (canonical
+ * 1997 equality, post Wave Y reformulation).
+ *
+ *   ΔF = −k_B T ln⟨exp(−W / (k_B T))⟩
+ *
+ * Reference: Jarzynski 1997 *Phys. Rev. Lett.* 78:2690 (canonical
+ * non-equilibrium free-energy equality; exact identity for any work
+ * protocol connecting two equilibrium states at temperature T).
+ *
+ * Status: speculative.
+ *
+ * Honest-claude scope notes:
+ *   - The Jarzynski equality itself is canonical; the *bridge framing*
+ *     — applying it to gravitational work in curved spacetime, the
+ *     original BE-29 framing — is the speculative element. The
+ *     gravitational-correction term `(β/(2c⁴)) ∫T^μν δg_μν √(-g) d⁴x`
+ *     is dropped here as the AST-unencodable extension.
+ *   - The AST has no `exp` primitive nor `ln` primitive; the average
+ *     `⟨exp(-βW)⟩` is encoded as a single dimensionless symbol stub
+ *     `'avg_exp_minus_betaW'`, with the exp argument β·W = W/(k_B T)
+ *     exposed as a separate ExprNode `BE29_BETAW_ARG` (verifiable
+ *     dimensionless via the dimensional analyzer).
+ *   - Both sides have dim [energy]: ΔF [J]; k_B T [J]; ⟨exp(-βW)⟩
+ *     dimensionless; ln(dimensionless) dimensionless; product [J]. ✓
+ *
+ * Numerical bracket: at quasi-static (reversible) work, W = W_rev and
+ * ⟨exp(-βW)⟩ = exp(-βW_rev), so ΔF = -k_BT · (-βW_rev) = W_rev. The
+ * Jarzynski equality reduces to the equilibrium identity in this limit.
+ *
+ * @see docs/specification/Part-II.md ("Bridge Equation 29")
+ * @see src/bridges/index.ts BRIDGE_EQUATIONS.find(e => e.id === 29)
+ * @module bridges/equations/be-29-jarzynski
+ */
+import { validate, validateEquation } from '../../dimensional/validator.js';
+import { DIMENSIONLESS, ENERGY, TEMPERATURE, } from '../../dimensional/types.js';
+import { k_B as DIM_kB } from '../../dimensional/constants.js';
+import { PhysicalConstants } from '../../core/types.js';
+const sym = (name, dim) => ({ kind: 'symbol', name, dim });
+/**
+ * Lemma AST: the exp-argument β·W = W/(k_B T), exposed for verification
+ * that the exponent is dimensionless.
+ *
+ * [W]/[k_B T] = [J]/[J] = [1] ✓
+ */
+export const BE29_BETAW_ARG = {
+    kind: 'op', op: '/',
+    args: [
+        sym('W', ENERGY),
+        {
+            kind: 'op', op: '*',
+            args: [
+                sym('k_B', DIM_kB),
+                sym('T', TEMPERATURE),
+            ],
+        },
+    ],
+};
+/**
+ * RHS of the Jarzynski equality: `−k_B T · ln⟨exp(-βW)⟩`.
+ *
+ * The AST has no `ln` primitive; we encode the relation as
+ *   `ΔF = -k_B T · L`
+ * where `L` is a dimensionless symbol stub representing
+ * `ln⟨exp(-W/(k_B T))⟩`. The minus sign is captured by encoding the
+ * coefficient as `-1` (a dimensionless symbol stub named `-1`).
+ *
+ * Both terms multiplicatively combine to dim [energy] = [k_B T].
+ */
+export const BE29_JARZYNSKI_RHS = {
+    kind: 'op', op: '*',
+    args: [
+        sym('-1', DIMENSIONLESS),
+        sym('k_B', DIM_kB),
+        sym('T', TEMPERATURE),
+        sym('ln_avg_exp_minus_betaW', DIMENSIONLESS),
+    ],
+};
+/** LHS: ΔF has dimension [energy]. */
+export const BE29_JARZYNSKI_LHS = sym('Delta_F', ENERGY);
+/**
+ * Evaluate the Jarzynski free-energy difference from a sample of work
+ * values:
+ *
+ *   ΔF = −k_B T ln⟨exp(−W / (k_B T))⟩
+ *
+ * The ensemble average is approximated by the sample mean over the
+ * `W_samples_J` array. For a reliable estimate the sample size must be
+ * large enough to capture the rare-event tail; the Jarzynski estimator
+ * has well-known finite-sample bias for low-temperature / high-work-
+ * variance protocols (Park-Schulten 2003 *J. Chem. Phys.* 119:5946).
+ *
+ * @returns Free-energy difference ΔF in joules.
+ */
+export function evaluateJarzynski(input) {
+    const { W_samples_J, T_K } = input;
+    if (!Number.isFinite(T_K) || T_K <= 0) {
+        throw new RangeError(`evaluateJarzynski: T_K must be a finite positive number, got ${T_K}`);
+    }
+    if (!Array.isArray(W_samples_J) || W_samples_J.length === 0) {
+        throw new RangeError('evaluateJarzynski: W_samples_J must be a non-empty array');
+    }
+    for (const W of W_samples_J) {
+        if (!Number.isFinite(W)) {
+            throw new RangeError(`evaluateJarzynski: every W in W_samples_J must be finite, got ${W}`);
+        }
+    }
+    const { kB } = PhysicalConstants;
+    const beta = 1 / (kB * T_K);
+    let sum = 0;
+    for (const W of W_samples_J) {
+        sum += Math.exp(-beta * W);
+    }
+    const avg = sum / W_samples_J.length;
+    return -kB * T_K * Math.log(avg);
+}
+// --- Self-validation ---
+/**
+ * Run the AST through the dimensional analyzer; both sides should be
+ * [energy].
+ */
+export function validateBE29Dimensions() {
+    const eq = validateEquation(BE29_JARZYNSKI_LHS, BE29_JARZYNSKI_RHS);
+    const lhs = validate(BE29_JARZYNSKI_LHS);
+    const rhs = validate(BE29_JARZYNSKI_RHS);
+    return {
+        ok: eq.ok,
+        lhsDim: lhs.inferredDimension,
+        rhsDim: rhs.inferredDimension,
+    };
+}
+//# sourceMappingURL=be-29-jarzynski.js.map

package/dist/bridges/equations/be-29-jarzynski.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"be-29-jarzynski.js","sourceRoot":"","sources":["../../../src/bridges/equations/be-29-jarzynski.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiCG;AAGH,OAAO,EAAE,QAAQ,EAAE,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AAC5E,OAAO,EAEL,aAAa,EACb,MAAM,EACN,WAAW,GACZ,MAAM,4BAA4B,CAAC;AACpC,OAAO,EAAE,GAAG,IAAI,MAAM,EAAE,MAAM,gCAAgC,CAAC;AAC/D,OAAO,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAExD,MAAM,GAAG,GAAG,CAAC,IAAY,EAAE,GAAc,EAAY,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,QAAQ,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC;AAExF;;;;;GAKG;AACH,MAAM,CAAC,MAAM,cAAc,GAAa;IACtC,IAAI,EAAE,IAAI,EAAE,EAAE,EAAE,GAAG;IACnB,IAAI,EAAE;QACJ,GAAG,CAAC,GAAG,EAAE,MAAM,CAAC;QAChB;YACE,IAAI,EAAE,IAAI,EAAE,EAAE,EAAE,GAAG;YACnB,IAAI,EAAE;gBACJ,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC;gBAClB,GAAG,CAAC,GAAG,EAAE,WAAW,CAAC;aACtB;SACF;KACF;CACF,CAAC;AAEF;;;;;;;;;;GAUG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAa;IAC1C,IAAI,EAAE,IAAI,EAAE,EAAE,EAAE,GAAG;IACnB,IAAI,EAAE;QACJ,GAAG,CAAC,IAAI,EAAE,aAAa,CAAC;QACxB,GAAG,CAAC,KAAK,EAAE,MAAM,CAAC;QAClB,GAAG,CAAC,GAAG,EAAE,WAAW,CAAC;QACrB,GAAG,CAAC,wBAAwB,EAAE,aAAa,CAAC;KAC7C;CACF,CAAC;AAEF,sCAAsC;AACtC,MAAM,CAAC,MAAM,kBAAkB,GAAa,GAAG,CAAC,SAAS,EAAE,MAAM,CAAC,CAAC;AAcnE;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,iBAAiB,CAAC,KAAsB;IACtD,MAAM,EAAE,WAAW,EAAE,GAAG,EAAE,GAAG,KAAK,CAAC;IACnC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;QACtC,MAAM,IAAI,UAAU,CAClB,gEAAgE,GAAG,EAAE,CACtE,CAAC;IACJ,CAAC;IACD,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,WAAW,CAAC,IAAI,WAAW,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC5D,MAAM,IAAI,UAAU,CAClB,0DAA0D,CAC3D,CAAC;IACJ,CAAC;IACD,KAAK,MAAM,CAAC,IAAI,WAAW,EAAE,CAAC;QAC5B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;YACxB,MAAM,IAAI,UAAU,CAClB,iEAAiE,CAAC,EAAE,CACrE,CAAC;QACJ,CAAC;IACH,CAAC;IACD,MAAM,EAAE,EAAE,EAAE,GAAG,iBAAiB,CAAC;IACjC,MAAM,IAAI,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;IAC5B,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,MAAM,CAAC,IAAI,WAAW,EAAE,CAAC;QAC5B,GAAG,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;IAC7B,CAAC;IACD,MAAM,GAAG,GAAG,GAAG,GAAG,WAAW,CAAC,MAAM,CAAC;IACrC,OAAO,CAAC,EAAE,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;AACnC,CAAC;AAED,0BAA0B;AAE1B;;;GAGG;AACH,MAAM,UAAU,sBAAsB;IACpC,MAAM,EAAE,GAAG,gBAAgB,CAAC,kBAAkB,EAAE,kBAAkB,CAAC,CAAC;IACpE,MAAM,GAAG,GAAG,QAAQ,CAAC,kBAAkB,CAAC,CAAC;IACzC,MAAM,GAAG,GAAG,QAAQ,CAAC,kBAAkB,CAAC,CAAC;IACzC,OAAO;QACL,EAAE,EAAE,EAAE,CAAC,EAAE;QACT,MAAM,EAAE,GAAG,CAAC,iBAAiB;QAC7B,MAAM,EAAE,GAAG,CAAC,iBAAiB;KAC9B,CAAC;AACJ,CAAC"}