universal-physics-tensor 4.0.0 → 5.0.0

package/dist/bridges/be102-landauer-conductance.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"be102-landauer-conductance.js","sourceRoot":"","sources":["../../src/bridges/be102-landauer-conductance.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AACH,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,sBAAsB,CAAC;AAqBlD;;;;GAIG;AACH,MAAM,UAAU,2BAA2B,CAAC,EAAE,MAAM,EAA6B;IAC/E,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;QAC7B,MAAM,IAAI,KAAK,CAAC,oDAAoD,CAAC,CAAC;IACxE,CAAC;IACD,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,MAAM,CAAC;IAC9C,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,CAAC;AACzB,CAAC"}

package/dist/bridges/be70-einstein-relation.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"be70-einstein-relation.d.ts","sourceRoot":"","sources":["../../src/bridges/be70-einstein-relation.ts"],"names":[],"mappings":"AAgBA;;;GAGG;AACH,MAAM,WAAW,sBAAsB;IACrC,qCAAqC;IACrC,QAAQ,CAAC,YAAY,EAAE,MAAM,CAAC;IAC9B,oCAAoC;IACpC,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,6EAA6E;IAC7E,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;CACtB;AAED;;;GAGG;AACH,MAAM,WAAW,sBAAsB;IACrC,QAAQ,CAAC,YAAY,EAAE,MAAM,CAAC;IAC9B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,yBAAyB;IACzB,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B;AAED;;;;GAIG;AACH,wBAAgB,wBAAwB,CAAC,EACvC,YAAY,EACZ,GAAG,EACH,GAAG,GACJ,EAAE,sBAAsB,GAAG,sBAAsB,CAejD"}
+{"version":3,"file":"be70-einstein-relation.d.ts","sourceRoot":"","sources":["../../src/bridges/be70-einstein-relation.ts"],"names":[],"mappings":"AAiBA;;;GAGG;AACH,MAAM,WAAW,sBAAsB;IACrC,qCAAqC;IACrC,QAAQ,CAAC,YAAY,EAAE,MAAM,CAAC;IAC9B,oCAAoC;IACpC,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,6EAA6E;IAC7E,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;CACtB;AAED;;;GAGG;AACH,MAAM,WAAW,sBAAsB;IACrC,QAAQ,CAAC,YAAY,EAAE,MAAM,CAAC;IAC9B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,yBAAyB;IACzB,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B;AAED;;;;GAIG;AACH,wBAAgB,wBAAwB,CAAC,EACvC,YAAY,EACZ,GAAG,EACH,GAAG,GACJ,EAAE,sBAAsB,GAAG,sBAAsB,CAajD"}

package/dist/bridges/be70-einstein-relation.js

@@ -13,6 +13,7 @@
  * @module bridges/be70-einstein-relation
  */
 import { K_B_SI } from '../core/constants.js';
+import { assertSameCarrierSign } from './carrier-sign.js';
 /**
  * Evaluate `D = μ k_B T / q`.
  *
@@ -28,9 +29,7 @@ export function evaluateEinsteinRelation({ mu_m2_per_Vs, T_K, q_C, }) {
     if (!Number.isFinite(q_C) || q_C === 0) {
         throw new Error('evaluateEinsteinRelation: q_C must be finite and nonzero');
     }
-    if (mu_m2_per_Vs * q_C < 0) {
-        throw new Error('evaluateEinsteinRelation: mu_m2_per_Vs and q_C must have the same sign');
-    }
+    assertSameCarrierSign(mu_m2_per_Vs, q_C, 'evaluateEinsteinRelation: mu_m2_per_Vs', 'q_C');
     const D_m2_per_s = (mu_m2_per_Vs * K_B_SI * T_K) / q_C;
     return { mu_m2_per_Vs, T_K, q_C, D_m2_per_s };
 }

package/dist/bridges/be70-einstein-relation.js.map

@@ -1 +1 @@
-{"version":3,"file":"be70-einstein-relation.js","sourceRoot":"","sources":["../../src/bridges/be70-einstein-relation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,sBAAsB,CAAC;AA2B9C;;;;GAIG;AACH,MAAM,UAAU,wBAAwB,CAAC,EACvC,YAAY,EACZ,GAAG,EACH,GAAG,GACoB;IACvB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,YAAY,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,uDAAuD,CAAC,CAAC;IAC3E,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;QACvC,MAAM,IAAI,KAAK,CAAC,0DAA0D,CAAC,CAAC;IAC9E,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;QACvC,MAAM,IAAI,KAAK,CAAC,0DAA0D,CAAC,CAAC;IAC9E,CAAC;IACD,IAAI,YAAY,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC;QAC3B,MAAM,IAAI,KAAK,CAAC,wEAAwE,CAAC,CAAC;IAC5F,CAAC;IACD,MAAM,UAAU,GAAG,CAAC,YAAY,GAAG,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;IACvD,OAAO,EAAE,YAAY,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC;AAChD,CAAC"}
+{"version":3,"file":"be70-einstein-relation.js","sourceRoot":"","sources":["../../src/bridges/be70-einstein-relation.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,sBAAsB,CAAC;AAC9C,OAAO,EAAE,qBAAqB,EAAE,MAAM,mBAAmB,CAAC;AA2B1D;;;;GAIG;AACH,MAAM,UAAU,wBAAwB,CAAC,EACvC,YAAY,EACZ,GAAG,EACH,GAAG,GACoB;IACvB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,YAAY,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,KAAK,CAAC,uDAAuD,CAAC,CAAC;IAC3E,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;QACvC,MAAM,IAAI,KAAK,CAAC,0DAA0D,CAAC,CAAC;IAC9E,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;QACvC,MAAM,IAAI,KAAK,CAAC,0DAA0D,CAAC,CAAC;IAC9E,CAAC;IACD,qBAAqB,CAAC,YAAY,EAAE,GAAG,EAAE,wCAAwC,EAAE,KAAK,CAAC,CAAC;IAC1F,MAAM,UAAU,GAAG,CAAC,YAAY,GAAG,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;IACvD,OAAO,EAAE,YAAY,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC;AAChD,CAAC"}

package/dist/bridges/be88-fermi-sea.d.ts

@@ -0,0 +1,31 @@
+/**
+ * Inputs for {@link evaluateFermiSea}.
+ * @public
+ */
+export interface FermiSeaInputs {
+    /** Number density of electrons, per cubic metre. */
+    readonly n_per_m3: number;
+    /** Isotropic effective mass, kilograms. */
+    readonly m_kg: number;
+}
+/**
+ * Result of {@link evaluateFermiSea}.
+ * @public
+ */
+export interface FermiSeaResult {
+    readonly n_per_m3: number;
+    readonly m_kg: number;
+    /** Fermi wavevector, per metre. */
+    readonly k_F_per_m: number;
+    /** Fermi energy on the isotropic parabola, joules. */
+    readonly E_F_J: number;
+    /** Fermi velocity, metres per second. */
+    readonly v_F_m_per_s: number;
+}
+/**
+ * Evaluate the two-spin Fermi sphere and the isotropic parabola at k_F.
+ *
+ * @public
+ */
+export declare function evaluateFermiSea({ n_per_m3, m_kg }: FermiSeaInputs): FermiSeaResult;
+//# sourceMappingURL=be88-fermi-sea.d.ts.map

package/dist/bridges/be88-fermi-sea.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"be88-fermi-sea.d.ts","sourceRoot":"","sources":["../../src/bridges/be88-fermi-sea.ts"],"names":[],"mappings":"AAkBA;;;GAGG;AACH,MAAM,WAAW,cAAc;IAC7B,oDAAoD;IACpD,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,2CAA2C;IAC3C,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;CACvB;AAED;;;GAGG;AACH,MAAM,WAAW,cAAc;IAC7B,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,mCAAmC;IACnC,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAC3B,sDAAsD;IACtD,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,yCAAyC;IACzC,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;CAC9B;AAED;;;;GAIG;AACH,wBAAgB,gBAAgB,CAAC,EAAE,QAAQ,EAAE,IAAI,EAAE,EAAE,cAAc,GAAG,cAAc,CAWnF"}

package/dist/bridges/be88-fermi-sea.js

@@ -0,0 +1,35 @@
+/**
+ * BE-88 — two-spin Fermi sphere,
+ *   k_F = (3 π² n)^{1/3},
+ *   E_F = ℏ² k_F² / (2 m*),
+ *   v_F = ℏ k_F / m*.
+ *
+ * Two spin states times the sphere (4π/3) k_F³/(2π)³ give n, so
+ * k_F³ = 3 π² n. One spin is k_F³ = 6 π² n. The isotropic parabola
+ * is a hypothesis, and so are the band, the two-spin count, and
+ * T = 0. Not a lattice band. ℏ⁻² times the second derivative of
+ * that parabola is 1/m*; this evaluator does not return that curvature.
+ *
+ * The overlay formalRef is `PhysJS.FermiSea.fermi_sea`.
+ *
+ * @module bridges/be88-fermi-sea
+ */
+import { HBAR_SI } from '../core/constants.js';
+/**
+ * Evaluate the two-spin Fermi sphere and the isotropic parabola at k_F.
+ *
+ * @public
+ */
+export function evaluateFermiSea({ n_per_m3, m_kg }) {
+    if (!Number.isFinite(n_per_m3) || n_per_m3 <= 0) {
+        throw new Error('evaluateFermiSea: n_per_m3 must be finite and > 0');
+    }
+    if (!Number.isFinite(m_kg) || m_kg <= 0) {
+        throw new Error('evaluateFermiSea: m_kg must be finite and > 0');
+    }
+    const k_F_per_m = (3 * Math.PI ** 2 * n_per_m3) ** (1 / 3);
+    const E_F_J = (HBAR_SI ** 2 * k_F_per_m ** 2) / (2 * m_kg);
+    const v_F_m_per_s = (HBAR_SI * k_F_per_m) / m_kg;
+    return { n_per_m3, m_kg, k_F_per_m, E_F_J, v_F_m_per_s };
+}
+//# sourceMappingURL=be88-fermi-sea.js.map

package/dist/bridges/be88-fermi-sea.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"be88-fermi-sea.js","sourceRoot":"","sources":["../../src/bridges/be88-fermi-sea.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AACH,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AA4B/C;;;;GAIG;AACH,MAAM,UAAU,gBAAgB,CAAC,EAAE,QAAQ,EAAE,IAAI,EAAkB;IACjE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,IAAI,CAAC,EAAE,CAAC;QAChD,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IACvE,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,IAAI,IAAI,CAAC,EAAE,CAAC;QACxC,MAAM,IAAI,KAAK,CAAC,+CAA+C,CAAC,CAAC;IACnE,CAAC;IACD,MAAM,SAAS,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,IAAI,CAAC,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3D,MAAM,KAAK,GAAG,CAAC,OAAO,IAAI,CAAC,GAAG,SAAS,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;IAC3D,MAAM,WAAW,GAAG,CAAC,OAAO,GAAG,SAAS,CAAC,GAAG,IAAI,CAAC;IACjD,OAAO,EAAE,QAAQ,EAAE,IAAI,EAAE,SAAS,EAAE,KAAK,EAAE,WAAW,EAAE,CAAC;AAC3D,CAAC"}

package/dist/bridges/be89-debye-cutoff.d.ts

@@ -0,0 +1,39 @@
+/**
+ * BE-89 — Debye cutoff of three acoustic branches,
+ *   ω_D = v_s (6 π² n)^{1/3}.
+ *
+ * Three branches filling 3n states, 3·(4π/3) k_D³/(2π)³ = 3n, give
+ * k_D³ = 6 π² n. Equating that three-branch sum to n gives
+ * k_D³ = 2 π² n. The branch count and the common speed are hypotheses.
+ *
+ * The overlay formalRef is `PhysJS.DebyeCutoff.debye_cutoff`.
+ *
+ * @module bridges/be89-debye-cutoff
+ */
+/**
+ * Inputs for {@link evaluateDebyeCutoff}.
+ * @public
+ */
+export interface DebyeCutoffInputs {
+    /** Common acoustic speed, metres per second. */
+    readonly v_m_per_s: number;
+    /** Number density of atoms, per cubic metre. */
+    readonly n_per_m3: number;
+}
+/**
+ * Result of {@link evaluateDebyeCutoff}.
+ * @public
+ */
+export interface DebyeCutoffResult {
+    readonly v_m_per_s: number;
+    readonly n_per_m3: number;
+    /** Debye frequency of one linear branch, radians per second. */
+    readonly omega_D_rad_per_s: number;
+}
+/**
+ * Evaluate ω_D = v_s (6 π² n)^{1/3}.
+ *
+ * @public
+ */
+export declare function evaluateDebyeCutoff({ v_m_per_s, n_per_m3 }: DebyeCutoffInputs): DebyeCutoffResult;
+//# sourceMappingURL=be89-debye-cutoff.d.ts.map

package/dist/bridges/be89-debye-cutoff.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"be89-debye-cutoff.d.ts","sourceRoot":"","sources":["../../src/bridges/be89-debye-cutoff.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH;;;GAGG;AACH,MAAM,WAAW,iBAAiB;IAChC,gDAAgD;IAChD,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAC3B,gDAAgD;IAChD,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;CAC3B;AAED;;;GAGG;AACH,MAAM,WAAW,iBAAiB;IAChC,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAC3B,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,gEAAgE;IAChE,QAAQ,CAAC,iBAAiB,EAAE,MAAM,CAAC;CACpC;AAED;;;;GAIG;AACH,wBAAgB,mBAAmB,CAAC,EAAE,SAAS,EAAE,QAAQ,EAAE,EAAE,iBAAiB,GAAG,iBAAiB,CASjG"}

package/dist/bridges/be89-debye-cutoff.js

@@ -0,0 +1,28 @@
+/**
+ * BE-89 — Debye cutoff of three acoustic branches,
+ *   ω_D = v_s (6 π² n)^{1/3}.
+ *
+ * Three branches filling 3n states, 3·(4π/3) k_D³/(2π)³ = 3n, give
+ * k_D³ = 6 π² n. Equating that three-branch sum to n gives
+ * k_D³ = 2 π² n. The branch count and the common speed are hypotheses.
+ *
+ * The overlay formalRef is `PhysJS.DebyeCutoff.debye_cutoff`.
+ *
+ * @module bridges/be89-debye-cutoff
+ */
+/**
+ * Evaluate ω_D = v_s (6 π² n)^{1/3}.
+ *
+ * @public
+ */
+export function evaluateDebyeCutoff({ v_m_per_s, n_per_m3 }) {
+    if (!Number.isFinite(v_m_per_s) || v_m_per_s <= 0) {
+        throw new Error('evaluateDebyeCutoff: v_m_per_s must be finite and > 0');
+    }
+    if (!Number.isFinite(n_per_m3) || n_per_m3 <= 0) {
+        throw new Error('evaluateDebyeCutoff: n_per_m3 must be finite and > 0');
+    }
+    const omega_D_rad_per_s = v_m_per_s * (6 * Math.PI ** 2 * n_per_m3) ** (1 / 3);
+    return { v_m_per_s, n_per_m3, omega_D_rad_per_s };
+}
+//# sourceMappingURL=be89-debye-cutoff.js.map

package/dist/bridges/be89-debye-cutoff.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"be89-debye-cutoff.js","sourceRoot":"","sources":["../../src/bridges/be89-debye-cutoff.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAwBH;;;;GAIG;AACH,MAAM,UAAU,mBAAmB,CAAC,EAAE,SAAS,EAAE,QAAQ,EAAqB;IAC5E,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,IAAI,SAAS,IAAI,CAAC,EAAE,CAAC;QAClD,MAAM,IAAI,KAAK,CAAC,uDAAuD,CAAC,CAAC;IAC3E,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,IAAI,CAAC,EAAE,CAAC;QAChD,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;IAC1E,CAAC;IACD,MAAM,iBAAiB,GAAG,SAAS,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,IAAI,CAAC,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC/E,OAAO,EAAE,SAAS,EAAE,QAAQ,EAAE,iBAAiB,EAAE,CAAC;AACpD,CAAC"}

package/dist/bridges/be90-debye-heat.d.ts

@@ -0,0 +1,30 @@
+/**
+ * Inputs for {@link evaluateDebyeHeat}.
+ * @public
+ */
+export interface DebyeHeatInputs {
+    /** Number of atoms. */
+    readonly N: number;
+    /** Absolute temperature, kelvin. */
+    readonly T_K: number;
+    /** Debye temperature, kelvin. */
+    readonly thetaD_K: number;
+}
+/**
+ * Result of {@link evaluateDebyeHeat}.
+ * @public
+ */
+export interface DebyeHeatResult {
+    readonly N: number;
+    readonly T_K: number;
+    readonly thetaD_K: number;
+    /** Heat capacity, joules per kelvin. */
+    readonly C_V_J_per_K: number;
+}
+/**
+ * Evaluate the low-temperature Debye heat capacity. π⁴/15 is assumed.
+ *
+ * @public
+ */
+export declare function evaluateDebyeHeat({ N, T_K, thetaD_K }: DebyeHeatInputs): DebyeHeatResult;
+//# sourceMappingURL=be90-debye-heat.d.ts.map

package/dist/bridges/be90-debye-heat.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"be90-debye-heat.d.ts","sourceRoot":"","sources":["../../src/bridges/be90-debye-heat.ts"],"names":[],"mappings":"AAiBA;;;GAGG;AACH,MAAM,WAAW,eAAe;IAC9B,uBAAuB;IACvB,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,oCAAoC;IACpC,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,iCAAiC;IACjC,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;CAC3B;AAED;;;GAGG;AACH,MAAM,WAAW,eAAe;IAC9B,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,wCAAwC;IACxC,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;CAC9B;AAED;;;;GAIG;AACH,wBAAgB,iBAAiB,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,QAAQ,EAAE,EAAE,eAAe,GAAG,eAAe,CAYxF"}

package/dist/bridges/be90-debye-heat.js

@@ -0,0 +1,35 @@
+/**
+ * BE-90 — Debye heat capacity at low temperature,
+ *   C_V = (12 π⁴ / 5) N k_B (T / θ_D)³.
+ *
+ * The mode integral ∫₀^{ω_D} 9 N ω²/ω_D³ dω = 3 N. The energy with the
+ * Bose integral extended to infinity is the hypothesis
+ * U = 9 N k_B T (T/θ_D)³ I, and I = π⁴/15 is a hypothesis, not an
+ * evaluation of ∫ x³/(exp(x)−1) dx. Nine times π⁴/15 is 3 π⁴/5, and
+ * U = A T⁴ differentiates to the heat capacity above. The energy
+ * prefactor 3 π⁴/5 is not this heat capacity.
+ *
+ * The overlay formalRef is `PhysJS.DebyeHeat.debye_heat`.
+ *
+ * @module bridges/be90-debye-heat
+ */
+import { K_B_SI } from '../core/constants.js';
+/**
+ * Evaluate the low-temperature Debye heat capacity. π⁴/15 is assumed.
+ *
+ * @public
+ */
+export function evaluateDebyeHeat({ N, T_K, thetaD_K }) {
+    if (!Number.isFinite(N) || N <= 0) {
+        throw new Error('evaluateDebyeHeat: N must be finite and > 0');
+    }
+    if (!Number.isFinite(T_K) || T_K <= 0) {
+        throw new Error('evaluateDebyeHeat: T_K must be finite and > 0');
+    }
+    if (!Number.isFinite(thetaD_K) || thetaD_K <= 0) {
+        throw new Error('evaluateDebyeHeat: thetaD_K must be finite and > 0');
+    }
+    const C_V_J_per_K = ((12 * Math.PI ** 4) / 5) * N * K_B_SI * (T_K / thetaD_K) ** 3;
+    return { N, T_K, thetaD_K, C_V_J_per_K };
+}
+//# sourceMappingURL=be90-debye-heat.js.map

package/dist/bridges/be90-debye-heat.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"be90-debye-heat.js","sourceRoot":"","sources":["../../src/bridges/be90-debye-heat.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,sBAAsB,CAAC;AA2B9C;;;;GAIG;AACH,MAAM,UAAU,iBAAiB,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,QAAQ,EAAmB;IACrE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,6CAA6C,CAAC,CAAC;IACjE,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;QACtC,MAAM,IAAI,KAAK,CAAC,+CAA+C,CAAC,CAAC;IACnE,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,IAAI,CAAC,EAAE,CAAC;QAChD,MAAM,IAAI,KAAK,CAAC,oDAAoD,CAAC,CAAC;IACxE,CAAC;IACD,MAAM,WAAW,GAAG,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC,GAAG,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC;IACnF,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,QAAQ,EAAE,WAAW,EAAE,CAAC;AAC3C,CAAC"}

package/dist/bridges/be91-einstein-solid.d.ts

@@ -0,0 +1,30 @@
+/**
+ * Inputs for {@link evaluateEinsteinSolid}.
+ * @public
+ */
+export interface EinsteinSolidInputs {
+    /** Number of atoms. */
+    readonly N: number;
+    /** Absolute temperature, kelvin. */
+    readonly T_K: number;
+    /** Einstein temperature, kelvin. */
+    readonly thetaE_K: number;
+}
+/**
+ * Result of {@link evaluateEinsteinSolid}.
+ * @public
+ */
+export interface EinsteinSolidResult {
+    readonly N: number;
+    readonly T_K: number;
+    readonly thetaE_K: number;
+    /** Heat capacity, joules per kelvin. */
+    readonly C_V_J_per_K: number;
+}
+/**
+ * Evaluate the three-oscillator Einstein heat capacity.
+ *
+ * @public
+ */
+export declare function evaluateEinsteinSolid({ N, T_K, thetaE_K }: EinsteinSolidInputs): EinsteinSolidResult;
+//# sourceMappingURL=be91-einstein-solid.d.ts.map

package/dist/bridges/be91-einstein-solid.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"be91-einstein-solid.d.ts","sourceRoot":"","sources":["../../src/bridges/be91-einstein-solid.ts"],"names":[],"mappings":"AAiBA;;;GAGG;AACH,MAAM,WAAW,mBAAmB;IAClC,uBAAuB;IACvB,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,oCAAoC;IACpC,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,oCAAoC;IACpC,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;CAC3B;AAED;;;GAGG;AACH,MAAM,WAAW,mBAAmB;IAClC,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,wCAAwC;IACxC,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;CAC9B;AAED;;;;GAIG;AACH,wBAAgB,qBAAqB,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,QAAQ,EAAE,EAAE,mBAAmB,GAAG,mBAAmB,CAcpG"}

package/dist/bridges/be91-einstein-solid.js

@@ -0,0 +1,37 @@
+/**
+ * BE-91 — Einstein solid heat capacity,
+ *   C_V = 3 N k_B (θ_E/T)² exp(θ_E/T) / (exp(θ_E/T) − 1)².
+ *
+ * Three Planck oscillators per atom, each of energy
+ * k_B θ_E / (exp(θ_E/T) − 1), differentiate to that heat capacity.
+ * The zero-point k_B θ_E/2 is constant and is not in C_V. The kernel
+ * tends to 1 as θ_E/T → 0⁺, so the high-temperature limit is 3 N k_B.
+ * One oscillator tends to N k_B. Three oscillators and the Einstein
+ * spectrum are hypotheses.
+ *
+ * The overlay formalRef is `PhysJS.EinsteinSolid.einstein_heat`.
+ *
+ * @module bridges/be91-einstein-solid
+ */
+import { K_B_SI } from '../core/constants.js';
+/**
+ * Evaluate the three-oscillator Einstein heat capacity.
+ *
+ * @public
+ */
+export function evaluateEinsteinSolid({ N, T_K, thetaE_K }) {
+    if (!Number.isFinite(N) || N <= 0) {
+        throw new Error('evaluateEinsteinSolid: N must be finite and > 0');
+    }
+    if (!Number.isFinite(T_K) || T_K <= 0) {
+        throw new Error('evaluateEinsteinSolid: T_K must be finite and > 0');
+    }
+    if (!Number.isFinite(thetaE_K) || thetaE_K <= 0) {
+        throw new Error('evaluateEinsteinSolid: thetaE_K must be finite and > 0');
+    }
+    const x = thetaE_K / T_K;
+    const ex = Math.exp(x);
+    const C_V_J_per_K = 3 * N * K_B_SI * (x * x * ex) / (ex - 1) ** 2;
+    return { N, T_K, thetaE_K, C_V_J_per_K };
+}
+//# sourceMappingURL=be91-einstein-solid.js.map

package/dist/bridges/be91-einstein-solid.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"be91-einstein-solid.js","sourceRoot":"","sources":["../../src/bridges/be91-einstein-solid.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,sBAAsB,CAAC;AA2B9C;;;;GAIG;AACH,MAAM,UAAU,qBAAqB,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,QAAQ,EAAuB;IAC7E,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QAClC,MAAM,IAAI,KAAK,CAAC,iDAAiD,CAAC,CAAC;IACrE,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;QACtC,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;IACvE,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,IAAI,CAAC,EAAE,CAAC;QAChD,MAAM,IAAI,KAAK,CAAC,wDAAwD,CAAC,CAAC;IAC5E,CAAC;IACD,MAAM,CAAC,GAAG,QAAQ,GAAG,GAAG,CAAC;IACzB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACvB,MAAM,WAAW,GAAG,CAAC,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;IAClE,OAAO,EAAE,CAAC,EAAE,GAAG,EAAE,QAAQ,EAAE,WAAW,EAAE,CAAC;AAC3C,CAAC"}

package/dist/bridges/be92-sommerfeld-heat.d.ts

@@ -0,0 +1,30 @@
+/**
+ * Inputs for {@link evaluateSommerfeldHeat}.
+ * @public
+ */
+export interface SommerfeldHeatInputs {
+    /** Electron number density, per cubic metre. */
+    readonly n_per_m3: number;
+    /** Absolute temperature, kelvin. */
+    readonly T_K: number;
+    /** Fermi energy, joules. */
+    readonly E_F_J: number;
+}
+/**
+ * Result of {@link evaluateSommerfeldHeat}.
+ * @public
+ */
+export interface SommerfeldHeatResult {
+    readonly n_per_m3: number;
+    readonly T_K: number;
+    readonly E_F_J: number;
+    /** Heat capacity per volume, joules per kelvin per cubic metre. */
+    readonly c_V_J_per_K_m3: number;
+}
+/**
+ * Evaluate the √E Sommerfeld heat capacity. The flat-density factor is not used.
+ *
+ * @public
+ */
+export declare function evaluateSommerfeldHeat({ n_per_m3, T_K, E_F_J }: SommerfeldHeatInputs): SommerfeldHeatResult;
+//# sourceMappingURL=be92-sommerfeld-heat.d.ts.map

package/dist/bridges/be92-sommerfeld-heat.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"be92-sommerfeld-heat.d.ts","sourceRoot":"","sources":["../../src/bridges/be92-sommerfeld-heat.ts"],"names":[],"mappings":"AAgBA;;;GAGG;AACH,MAAM,WAAW,oBAAoB;IACnC,gDAAgD;IAChD,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,oCAAoC;IACpC,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,4BAA4B;IAC5B,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;CACxB;AAED;;;GAGG;AACH,MAAM,WAAW,oBAAoB;IACnC,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,mEAAmE;IACnE,QAAQ,CAAC,cAAc,EAAE,MAAM,CAAC;CACjC;AAED;;;;GAIG;AACH,wBAAgB,sBAAsB,CAAC,EAAE,QAAQ,EAAE,GAAG,EAAE,KAAK,EAAE,EAAE,oBAAoB,GAAG,oBAAoB,CAY3G"}

package/dist/bridges/be92-sommerfeld-heat.js

@@ -0,0 +1,34 @@
+/**
+ * BE-92 — electronic heat capacity of a √E density of states,
+ *   c_V = (π²/2) n k_B² T / E_F.
+ *
+ * The Sommerfeld energy correction δU = (π²/6) (k_B T)² g(E_F) is a
+ * hypothesis. Its temperature derivative is c_V = (π²/3) k_B² T g(E_F).
+ * PhysJS.FermiSea.dos_factor supplies g(E_F) = (3/2) n/E_F, which
+ * produces the factor π²/2. A flat density g = n/E_F leaves π²/3.
+ * Not the Wiedemann–Franz law and not a second proof of be-61.
+ *
+ * The overlay formalRef is `PhysJS.SommerfeldHeat.electronic_heat`.
+ *
+ * @module bridges/be92-sommerfeld-heat
+ */
+import { K_B_SI } from '../core/constants.js';
+/**
+ * Evaluate the √E Sommerfeld heat capacity. The flat-density factor is not used.
+ *
+ * @public
+ */
+export function evaluateSommerfeldHeat({ n_per_m3, T_K, E_F_J }) {
+    if (!Number.isFinite(n_per_m3) || n_per_m3 <= 0) {
+        throw new Error('evaluateSommerfeldHeat: n_per_m3 must be finite and > 0');
+    }
+    if (!Number.isFinite(T_K) || T_K <= 0) {
+        throw new Error('evaluateSommerfeldHeat: T_K must be finite and > 0');
+    }
+    if (!Number.isFinite(E_F_J) || E_F_J <= 0) {
+        throw new Error('evaluateSommerfeldHeat: E_F_J must be finite and > 0');
+    }
+    const c_V_J_per_K_m3 = ((Math.PI ** 2) / 2) * n_per_m3 * K_B_SI ** 2 * T_K / E_F_J;
+    return { n_per_m3, T_K, E_F_J, c_V_J_per_K_m3 };
+}
+//# sourceMappingURL=be92-sommerfeld-heat.js.map

package/dist/bridges/be92-sommerfeld-heat.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"be92-sommerfeld-heat.js","sourceRoot":"","sources":["../../src/bridges/be92-sommerfeld-heat.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,sBAAsB,CAAC;AA2B9C;;;;GAIG;AACH,MAAM,UAAU,sBAAsB,CAAC,EAAE,QAAQ,EAAE,GAAG,EAAE,KAAK,EAAwB;IACnF,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,IAAI,CAAC,EAAE,CAAC;QAChD,MAAM,IAAI,KAAK,CAAC,yDAAyD,CAAC,CAAC;IAC7E,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,EAAE,CAAC;QACtC,MAAM,IAAI,KAAK,CAAC,oDAAoD,CAAC,CAAC;IACxE,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;QAC1C,MAAM,IAAI,KAAK,CAAC,sDAAsD,CAAC,CAAC;IAC1E,CAAC;IACD,MAAM,cAAc,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,GAAG,MAAM,IAAI,CAAC,GAAG,GAAG,GAAG,KAAK,CAAC;IACnF,OAAO,EAAE,QAAQ,EAAE,GAAG,EAAE,KAAK,EAAE,cAAc,EAAE,CAAC;AAClD,CAAC"}

package/dist/bridges/be93-curie-weiss.d.ts

@@ -0,0 +1,41 @@
+/**
+ * Inputs for {@link evaluateCurieWeiss}.
+ * @public
+ */
+export interface CurieWeissInputs {
+    /** Moment density, per cubic metre. */
+    readonly n_per_m3: number;
+    /** Landé g-factor. */
+    readonly g: number;
+    /** Spin. S = 1/2 is the equal-weight check, not a restriction. */
+    readonly spin: number;
+    /** Bohr magneton, joules per tesla. */
+    readonly muB_J_per_T: number;
+    /** Absolute temperature, kelvin. */
+    readonly T_K: number;
+    /** Weiss temperature, kelvin. Zero is the Curie law. */
+    readonly theta_K: number;
+}
+/**
+ * Result of {@link evaluateCurieWeiss}.
+ * @public
+ */
+export interface CurieWeissResult {
+    readonly n_per_m3: number;
+    readonly g: number;
+    readonly spin: number;
+    readonly muB_J_per_T: number;
+    readonly T_K: number;
+    readonly theta_K: number;
+    /** Curie constant, kelvin. */
+    readonly C_K: number;
+    /** Dimensionless susceptibility χ = C/(T−θ). */
+    readonly chi: number;
+}
+/**
+ * Evaluate the Curie–Weiss susceptibility. θ is the mean-field hypothesis.
+ *
+ * @public
+ */
+export declare function evaluateCurieWeiss(inputs: CurieWeissInputs): CurieWeissResult;
+//# sourceMappingURL=be93-curie-weiss.d.ts.map

package/dist/bridges/be93-curie-weiss.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"be93-curie-weiss.d.ts","sourceRoot":"","sources":["../../src/bridges/be93-curie-weiss.ts"],"names":[],"mappings":"AAmBA;;;GAGG;AACH,MAAM,WAAW,gBAAgB;IAC/B,uCAAuC;IACvC,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,sBAAsB;IACtB,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,kEAAkE;IAClE,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,uCAAuC;IACvC,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;IAC7B,oCAAoC;IACpC,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,wDAAwD;IACxD,QAAQ,CAAC,OAAO,EAAE,MAAM,CAAC;CAC1B;AAED;;;GAGG;AACH,MAAM,WAAW,gBAAgB;IAC/B,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,QAAQ,CAAC,CAAC,EAAE,MAAM,CAAC;IACnB,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;IAC7B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,OAAO,EAAE,MAAM,CAAC;IACzB,8BAA8B;IAC9B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,gDAAgD;IAChD,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;CACtB;AAED;;;;GAIG;AACH,wBAAgB,kBAAkB,CAAC,MAAM,EAAE,gBAAgB,GAAG,gBAAgB,CAoB7E"}

package/dist/bridges/be93-curie-weiss.js

@@ -0,0 +1,45 @@
+/**
+ * BE-93 — Curie–Weiss susceptibility,
+ *   C = μ₀ n g² μ_B² S(S+1) / (3 k_B),
+ *   χ = C / (T − θ).
+ *
+ * Linear response χ k_B T = μ₀ n (g μ_B)² ⟨S_z²⟩ with the
+ * high-temperature moment ⟨S_z²⟩ = S(S+1)/3 gives C. Equal weights
+ * on m = ±1/2 give 1/4 = S(S+1)/3 at S = 1/2. Mean field
+ * B_eff = B + λ M with θ = C λ/μ₀ gives χ = C/(T−θ). θ = 0 is C/T.
+ * A classical moment uses μ²/3. The second moment and the mean-field
+ * shift are hypotheses. Not an su(2) derivation. μ_B is an input.
+ *
+ * The overlay formalRef is `PhysJS.CurieWeiss.curie_weiss`.
+ *
+ * @module bridges/be93-curie-weiss
+ */
+import { K_B_SI } from '../core/constants.js';
+import { MU0_SI } from '../dimensional/formula-names.js';
+/**
+ * Evaluate the Curie–Weiss susceptibility. θ is the mean-field hypothesis.
+ *
+ * @public
+ */
+export function evaluateCurieWeiss(inputs) {
+    const { n_per_m3, g, spin, muB_J_per_T, T_K, theta_K } = inputs;
+    if (!Number.isFinite(n_per_m3) || n_per_m3 <= 0) {
+        throw new Error('evaluateCurieWeiss: n_per_m3 must be finite and > 0');
+    }
+    if (!Number.isFinite(g)) {
+        throw new Error('evaluateCurieWeiss: g must be finite');
+    }
+    if (!Number.isFinite(spin) || spin < 0) {
+        throw new Error('evaluateCurieWeiss: spin must be finite and ≥ 0');
+    }
+    if (!Number.isFinite(muB_J_per_T) || muB_J_per_T <= 0) {
+        throw new Error('evaluateCurieWeiss: muB_J_per_T must be finite and > 0');
+    }
+    if (!Number.isFinite(T_K) || !Number.isFinite(theta_K) || T_K === theta_K) {
+        throw new Error('evaluateCurieWeiss: T_K and theta_K must be finite and distinct');
+    }
+    const C_K = (MU0_SI * n_per_m3 * g * g * muB_J_per_T ** 2 * spin * (spin + 1)) / (3 * K_B_SI);
+    const chi = C_K / (T_K - theta_K);
+    return { ...inputs, C_K, chi };
+}
+//# sourceMappingURL=be93-curie-weiss.js.map

package/dist/bridges/be93-curie-weiss.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"be93-curie-weiss.js","sourceRoot":"","sources":["../../src/bridges/be93-curie-weiss.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,sBAAsB,CAAC;AAC9C,OAAO,EAAE,MAAM,EAAE,MAAM,iCAAiC,CAAC;AAsCzD;;;;GAIG;AACH,MAAM,UAAU,kBAAkB,CAAC,MAAwB;IACzD,MAAM,EAAE,QAAQ,EAAE,CAAC,EAAE,IAAI,EAAE,WAAW,EAAE,GAAG,EAAE,OAAO,EAAE,GAAG,MAAM,CAAC;IAChE,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,IAAI,CAAC,EAAE,CAAC;QAChD,MAAM,IAAI,KAAK,CAAC,qDAAqD,CAAC,CAAC;IACzE,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;QACxB,MAAM,IAAI,KAAK,CAAC,sCAAsC,CAAC,CAAC;IAC1D,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC;QACvC,MAAM,IAAI,KAAK,CAAC,iDAAiD,CAAC,CAAC;IACrE,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC,IAAI,WAAW,IAAI,CAAC,EAAE,CAAC;QACtD,MAAM,IAAI,KAAK,CAAC,wDAAwD,CAAC,CAAC;IAC5E,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,OAAO,CAAC,IAAI,GAAG,KAAK,OAAO,EAAE,CAAC;QAC1E,MAAM,IAAI,KAAK,CAAC,iEAAiE,CAAC,CAAC;IACrF,CAAC;IACD,MAAM,GAAG,GAAG,CAAC,MAAM,GAAG,QAAQ,GAAG,CAAC,GAAG,CAAC,GAAG,WAAW,IAAI,CAAC,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC;IAC9F,MAAM,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC,CAAC;IAClC,OAAO,EAAE,GAAG,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;AACjC,CAAC"}

package/dist/bridges/be94-pauli-paramagnetism.d.ts

@@ -0,0 +1,30 @@
+/**
+ * Inputs for {@link evaluatePauliParamagnetism}.
+ * @public
+ */
+export interface PauliParamagnetismInputs {
+    /** Electron number density, per cubic metre. */
+    readonly n_per_m3: number;
+    /** Fermi energy, joules. */
+    readonly E_F_J: number;
+    /** Bohr magneton, joules per tesla. */
+    readonly muB_J_per_T: number;
+}
+/**
+ * Result of {@link evaluatePauliParamagnetism}.
+ * @public
+ */
+export interface PauliParamagnetismResult {
+    readonly n_per_m3: number;
+    readonly E_F_J: number;
+    readonly muB_J_per_T: number;
+    /** Dimensionless Pauli susceptibility. */
+    readonly chi_P: number;
+}
+/**
+ * Evaluate the √E Pauli susceptibility. Not χ_L = −χ_P/3.
+ *
+ * @public
+ */
+export declare function evaluatePauliParamagnetism({ n_per_m3, E_F_J, muB_J_per_T, }: PauliParamagnetismInputs): PauliParamagnetismResult;
+//# sourceMappingURL=be94-pauli-paramagnetism.d.ts.map

package/dist/bridges/be94-pauli-paramagnetism.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"be94-pauli-paramagnetism.d.ts","sourceRoot":"","sources":["../../src/bridges/be94-pauli-paramagnetism.ts"],"names":[],"mappings":"AAeA;;;GAGG;AACH,MAAM,WAAW,wBAAwB;IACvC,gDAAgD;IAChD,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,4BAA4B;IAC5B,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,uCAAuC;IACvC,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;CAC9B;AAED;;;GAGG;AACH,MAAM,WAAW,wBAAwB;IACvC,QAAQ,CAAC,QAAQ,EAAE,MAAM,CAAC;IAC1B,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;IAC7B,0CAA0C;IAC1C,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;CACxB;AAED;;;;GAIG;AACH,wBAAgB,0BAA0B,CAAC,EACzC,QAAQ,EACR,KAAK,EACL,WAAW,GACZ,EAAE,wBAAwB,GAAG,wBAAwB,CAYrD"}

package/dist/bridges/be94-pauli-paramagnetism.js

@@ -0,0 +1,33 @@
+/**
+ * BE-94 — Pauli paramagnetism for a √E density of states,
+ *   χ_P = μ₀ μ_B² (3 n) / (2 E_F).
+ *
+ * The Zeeman imbalance M = μ_B² g(E_F) B is a hypothesis, and
+ * χ_P = μ₀ M/B is μ₀ μ_B² g(E_F). dos_factor supplies
+ * g(E_F) = (3/2) n/E_F. A flat density leaves the factor 1,
+ * χ_P = μ₀ μ_B² n/E_F. Not Landau diamagnetism. μ_B is an input.
+ *
+ * The overlay formalRef is `PhysJS.PauliParamagnetism.pauli`.
+ *
+ * @module bridges/be94-pauli-paramagnetism
+ */
+import { MU0_SI } from '../dimensional/formula-names.js';
+/**
+ * Evaluate the √E Pauli susceptibility. Not χ_L = −χ_P/3.
+ *
+ * @public
+ */
+export function evaluatePauliParamagnetism({ n_per_m3, E_F_J, muB_J_per_T, }) {
+    if (!Number.isFinite(n_per_m3) || n_per_m3 <= 0) {
+        throw new Error('evaluatePauliParamagnetism: n_per_m3 must be finite and > 0');
+    }
+    if (!Number.isFinite(E_F_J) || E_F_J <= 0) {
+        throw new Error('evaluatePauliParamagnetism: E_F_J must be finite and > 0');
+    }
+    if (!Number.isFinite(muB_J_per_T) || muB_J_per_T <= 0) {
+        throw new Error('evaluatePauliParamagnetism: muB_J_per_T must be finite and > 0');
+    }
+    const chi_P = (MU0_SI * muB_J_per_T ** 2 * 3 * n_per_m3) / (2 * E_F_J);
+    return { n_per_m3, E_F_J, muB_J_per_T, chi_P };
+}
+//# sourceMappingURL=be94-pauli-paramagnetism.js.map

package/dist/bridges/be94-pauli-paramagnetism.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"be94-pauli-paramagnetism.js","sourceRoot":"","sources":["../../src/bridges/be94-pauli-paramagnetism.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,iCAAiC,CAAC;AA2BzD;;;;GAIG;AACH,MAAM,UAAU,0BAA0B,CAAC,EACzC,QAAQ,EACR,KAAK,EACL,WAAW,GACc;IACzB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,QAAQ,IAAI,CAAC,EAAE,CAAC;QAChD,MAAM,IAAI,KAAK,CAAC,6DAA6D,CAAC,CAAC;IACjF,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;QAC1C,MAAM,IAAI,KAAK,CAAC,0DAA0D,CAAC,CAAC;IAC9E,CAAC;IACD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAC,IAAI,WAAW,IAAI,CAAC,EAAE,CAAC;QACtD,MAAM,IAAI,KAAK,CAAC,gEAAgE,CAAC,CAAC;IACpF,CAAC;IACD,MAAM,KAAK,GAAG,CAAC,MAAM,GAAG,WAAW,IAAI,CAAC,GAAG,CAAC,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;IACvE,OAAO,EAAE,QAAQ,EAAE,KAAK,EAAE,WAAW,EAAE,KAAK,EAAE,CAAC;AACjD,CAAC"}

package/dist/bridges/be95-ginzburg-landau.d.ts

@@ -0,0 +1,46 @@
+/**
+ * BE-95 — sign of one Ginzburg–Landau trial wall,
+ *   (1/κ² − 2) times the gradient integral.
+ *
+ * In the normalization with gradient coefficient 1/κ², quartic
+ * (1/2)(1−f²)², and field B², the density at κ² = 1/2 is
+ * (√2 f' − a f)² + (B + (1−f²)/√2)² minus √2 times the derivative
+ * of a(1−f²). Vanishing squares and equal endpoints make that wall
+ * integral zero. A trial profile with that critical integral has
+ * energy (1/κ² − 2) times the gradient integral: negative when
+ * κ > 1/√2, zero at κ = 1/√2, and positive when κ < 1/√2. The
+ * positive side is this trial, not every minimizer. The GL density
+ * and the profile are hypotheses. This evaluator returns the
+ * dimensionless factor, not a search over minimizers.
+ *
+ * The overlay formalRef is `PhysJS.GinzburgLandau.type_boundary`.
+ *
+ * @module bridges/be95-ginzburg-landau
+ */
+/**
+ * Inputs for {@link evaluateGinzburgLandau}.
+ * @public
+ */
+export interface GinzburgLandauInputs {
+    /** Ginzburg–Landau parameter κ. */
+    readonly kappa: number;
+}
+/**
+ * Result of {@link evaluateGinzburgLandau}.
+ * @public
+ */
+export interface GinzburgLandauResult {
+    readonly kappa: number;
+    /**
+     * Factor (1/κ² − 2) in front of the positive gradient integral
+     * of the trial wall.
+     */
+    readonly trial_factor: number;
+}
+/**
+ * Evaluate the trial-wall factor (1/κ² − 2). Not the energy of every minimizer.
+ *
+ * @public
+ */
+export declare function evaluateGinzburgLandau({ kappa }: GinzburgLandauInputs): GinzburgLandauResult;
+//# sourceMappingURL=be95-ginzburg-landau.d.ts.map

package/dist/bridges/be95-ginzburg-landau.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"be95-ginzburg-landau.d.ts","sourceRoot":"","sources":["../../src/bridges/be95-ginzburg-landau.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;GAkBG;AAEH;;;GAGG;AACH,MAAM,WAAW,oBAAoB;IACnC,mCAAmC;IACnC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;CACxB;AAED;;;GAGG;AACH,MAAM,WAAW,oBAAoB;IACnC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB;;;OAGG;IACH,QAAQ,CAAC,YAAY,EAAE,MAAM,CAAC;CAC/B;AAED;;;;GAIG;AACH,wBAAgB,sBAAsB,CAAC,EAAE,KAAK,EAAE,EAAE,oBAAoB,GAAG,oBAAoB,CAM5F"}