universal-physics-tensor 4.0.0 → 6.0.0

package/dist/bridges/be77-hagen-poiseuille.d.ts

@@ -14,7 +14,7 @@
  */
 /**
  * Inputs for {@link evaluateHagenPoiseuille}.
- * @public
+ * @internal
  */
 export interface HagenPoiseuilleInputs {
     /** Pipe radius, metres. */
@@ -28,7 +28,7 @@ export interface HagenPoiseuilleInputs {
 }
 /**
  * Result of {@link evaluateHagenPoiseuille}.
- * @public
+ * @internal
  */
 export interface HagenPoiseuilleResult {
     readonly R_m: number;
@@ -44,7 +44,7 @@ export interface HagenPoiseuilleResult {
  * The 8 is the integral of the no-slip parabola. It is not the Fanning
  * factor 16, and it is not a square-duct eigenvalue.
  *
- * @public
+ * @internal
  */
 export declare function evaluateHagenPoiseuille({ R_m, deltaP_Pa, mu_Pa_s, L_m, }: HagenPoiseuilleInputs): HagenPoiseuilleResult;
 //# sourceMappingURL=be77-hagen-poiseuille.d.ts.map

package/dist/bridges/be77-hagen-poiseuille.js

@@ -18,7 +18,7 @@
  * The 8 is the integral of the no-slip parabola. It is not the Fanning
  * factor 16, and it is not a square-duct eigenvalue.
  *
- * @public
+ * @internal
  */
 export function evaluateHagenPoiseuille({ R_m, deltaP_Pa, mu_Pa_s, L_m, }) {
     if (!Number.isFinite(R_m) || R_m <= 0) {

package/dist/bridges/be78-euler-buckling.d.ts

@@ -14,7 +14,7 @@
  */
 /**
  * Inputs for {@link evaluateEulerBuckling}.
- * @public
+ * @internal
  */
 export interface EulerBucklingInputs {
     /** Young's modulus, pascal. */
@@ -26,7 +26,7 @@ export interface EulerBucklingInputs {
 }
 /**
  * Result of {@link evaluateEulerBuckling}.
- * @public
+ * @internal
  */
 export interface EulerBucklingResult {
     readonly E_Pa: number;
@@ -38,7 +38,7 @@ export interface EulerBucklingResult {
 /**
  * Evaluate `P_cr = π² E I / L²` for pinned ends.
  *
- * @public
+ * @internal
  */
 export declare function evaluateEulerBuckling({ E_Pa, I_m4, L_m }: EulerBucklingInputs): EulerBucklingResult;
 //# sourceMappingURL=be78-euler-buckling.d.ts.map

package/dist/bridges/be78-euler-buckling.js

@@ -15,7 +15,7 @@
 /**
  * Evaluate `P_cr = π² E I / L²` for pinned ends.
  *
- * @public
+ * @internal
  */
 export function evaluateEulerBuckling({ E_Pa, I_m4, L_m }) {
     if (!Number.isFinite(E_Pa) || E_Pa <= 0) {

package/dist/bridges/be79-pull-in.d.ts

@@ -1,6 +1,6 @@
 /**
  * Inputs for {@link evaluatePullIn}.
- * @public
+ * @internal
  */
 export interface PullInInputs {
     /** Linear spring stiffness, newtons per metre. */
@@ -12,7 +12,7 @@ export interface PullInInputs {
 }
 /**
  * Result of {@link evaluatePullIn}.
- * @public
+ * @internal
  */
 export interface PullInResult {
     readonly k_N_per_m: number;
@@ -27,7 +27,7 @@ export interface PullInResult {
  * `ε0` is the vacuum permittivity. The gap at the fold is `2 g0/3`,
  * not `g0/2`. Fringing is not in the theorem.
  *
- * @public
+ * @internal
  */
 export declare function evaluatePullIn({ k_N_per_m, g0_m, A_m2 }: PullInInputs): PullInResult;
 //# sourceMappingURL=be79-pull-in.d.ts.map

package/dist/bridges/be79-pull-in.js

@@ -17,7 +17,7 @@ import { EPS0_SI } from '../dimensional/formula-names.js';
  * `ε0` is the vacuum permittivity. The gap at the fold is `2 g0/3`,
  * not `g0/2`. Fringing is not in the theorem.
  *
- * @public
+ * @internal
  */
 export function evaluatePullIn({ k_N_per_m, g0_m, A_m2 }) {
     if (!Number.isFinite(k_N_per_m) || k_N_per_m <= 0) {

package/dist/bridges/be80-mott-gurney.d.ts

@@ -13,7 +13,7 @@
  */
 /**
  * Inputs for {@link evaluateMottGurney}.
- * @public
+ * @internal
  */
 export interface MottGurneyInputs {
     /**
@@ -30,7 +30,7 @@ export interface MottGurneyInputs {
 }
 /**
  * Result of {@link evaluateMottGurney}.
- * @public
+ * @internal
  */
 export interface MottGurneyResult {
     readonly eps: number;
@@ -46,7 +46,7 @@ export interface MottGurneyResult {
  * Drift, Poisson, and the injecting contact E(0) = 0 are hypotheses.
  * This is not the Child–Langmuir current.
  *
- * @public
+ * @internal
  */
 export declare function evaluateMottGurney({ eps, mu_m2_per_Vs, V_volts, d_m, }: MottGurneyInputs): MottGurneyResult;
 //# sourceMappingURL=be80-mott-gurney.d.ts.map

package/dist/bridges/be80-mott-gurney.js

@@ -17,7 +17,7 @@
  * Drift, Poisson, and the injecting contact E(0) = 0 are hypotheses.
  * This is not the Child–Langmuir current.
  *
- * @public
+ * @internal
  */
 export function evaluateMottGurney({ eps, mu_m2_per_Vs, V_volts, d_m, }) {
     if (!Number.isFinite(eps) || eps <= 0) {

package/dist/bridges/be81-child-langmuir.d.ts

@@ -1,6 +1,6 @@
 /**
  * Inputs for {@link evaluateChildLangmuir}.
- * @public
+ * @internal
  */
 export interface ChildLangmuirInputs {
     /** Particle mass, kilograms. */
@@ -12,7 +12,7 @@ export interface ChildLangmuirInputs {
 }
 /**
  * Result of {@link evaluateChildLangmuir}.
- * @public
+ * @internal
  */
 export interface ChildLangmuirResult {
     readonly m_kg: number;
@@ -27,7 +27,7 @@ export interface ChildLangmuirResult {
  * `e` is `E_SI`, the elementary charge, not an input and not Euler's
  * number. `ε0` is the vacuum permittivity.
  *
- * @public
+ * @internal
  */
 export declare function evaluateChildLangmuir({ m_kg, V_volts, d_m }: ChildLangmuirInputs): ChildLangmuirResult;
 //# sourceMappingURL=be81-child-langmuir.d.ts.map

package/dist/bridges/be81-child-langmuir.js

@@ -20,7 +20,7 @@ import { EPS0_SI } from '../dimensional/formula-names.js';
  * `e` is `E_SI`, the elementary charge, not an input and not Euler's
  * number. `ε0` is the vacuum permittivity.
  *
- * @public
+ * @internal
  */
 export function evaluateChildLangmuir({ m_kg, V_volts, d_m }) {
     if (!Number.isFinite(m_kg) || m_kg <= 0) {

package/dist/bridges/be82-shockley-diode.d.ts

@@ -1,6 +1,6 @@
 /**
  * Inputs for {@link evaluateShockleyDiode}.
- * @public
+ * @internal
  */
 export interface ShockleyDiodeInputs {
     /** Saturation current, amperes. */
@@ -12,7 +12,7 @@ export interface ShockleyDiodeInputs {
 }
 /**
  * Result of {@link evaluateShockleyDiode}.
- * @public
+ * @internal
  */
 export interface ShockleyDiodeResult {
     readonly I_s_A: number;
@@ -26,7 +26,7 @@ export interface ShockleyDiodeResult {
  *
  * `e` is the elementary charge. Euler's number is `exp`, not a bare `e`.
  *
- * @public
+ * @internal
  */
 export declare function evaluateShockleyDiode({ I_s_A, V_volts, T_K }: ShockleyDiodeInputs): ShockleyDiodeResult;
 //# sourceMappingURL=be82-shockley-diode.d.ts.map

package/dist/bridges/be82-shockley-diode.js

@@ -19,7 +19,7 @@ import { E_SI, K_B_SI } from '../core/constants.js';
  *
  * `e` is the elementary charge. Euler's number is `exp`, not a bare `e`.
  *
- * @public
+ * @internal
  */
 export function evaluateShockleyDiode({ I_s_A, V_volts, T_K }) {
     if (!Number.isFinite(I_s_A)) {

package/dist/bridges/be83-thomson.d.ts

@@ -13,7 +13,7 @@
  */
 /**
  * Inputs for {@link evaluateThomsonCoefficient}.
- * @public
+ * @internal
  */
 export interface ThomsonCoefficientInputs {
     /** Absolute temperature, kelvin. */
@@ -26,7 +26,7 @@ export interface ThomsonCoefficientInputs {
 }
 /**
  * Result of {@link evaluateThomsonCoefficient}.
- * @public
+ * @internal
  */
 export interface ThomsonCoefficientResult {
     readonly T_K: number;
@@ -40,7 +40,7 @@ export interface ThomsonCoefficientResult {
  * The functional Kelvin relation and the Thomson split are hypotheses.
  * This is not Π = S T.
  *
- * @public
+ * @internal
  */
 export declare function evaluateThomsonCoefficient({ T_K, dS_dT_V_per_K2, }: ThomsonCoefficientInputs): ThomsonCoefficientResult;
 //# sourceMappingURL=be83-thomson.d.ts.map

package/dist/bridges/be83-thomson.js

@@ -17,7 +17,7 @@
  * The functional Kelvin relation and the Thomson split are hypotheses.
  * This is not Π = S T.
  *
- * @public
+ * @internal
  */
 export function evaluateThomsonCoefficient({ T_K, dS_dT_V_per_K2, }) {
     if (!Number.isFinite(T_K)) {

package/dist/bridges/be84-four-point.d.ts

@@ -15,7 +15,7 @@
  */
 /**
  * Inputs for {@link evaluateFourPointSheet}.
- * @public
+ * @internal
  */
 export interface FourPointSheetInputs {
     /** Inner-pair voltage, volts. */
@@ -25,7 +25,7 @@ export interface FourPointSheetInputs {
 }
 /**
  * Result of {@link evaluateFourPointSheet}.
- * @public
+ * @internal
  */
 export interface FourPointSheetResult {
     readonly V_volts: number;
@@ -39,7 +39,7 @@ export interface FourPointSheetResult {
  * The radial 1/r field and linear superposition are hypotheses. The
  * spacing s cancels and is not an input. A sink at 4s is not this factor.
  *
- * @public
+ * @internal
  */
 export declare function evaluateFourPointSheet({ V_volts, I_A }: FourPointSheetInputs): FourPointSheetResult;
 //# sourceMappingURL=be84-four-point.d.ts.map

package/dist/bridges/be84-four-point.js

@@ -19,7 +19,7 @@
  * The radial 1/r field and linear superposition are hypotheses. The
  * spacing s cancels and is not an input. A sink at 4s is not this factor.
  *
- * @public
+ * @internal
  */
 export function evaluateFourPointSheet({ V_volts, I_A }) {
     if (!Number.isFinite(V_volts)) {

package/dist/bridges/be85-shot-noise.d.ts

@@ -1,6 +1,6 @@
 /**
  * Inputs for {@link evaluateShotNoise}.
- * @public
+ * @internal
  */
 export interface ShotNoiseInputs {
     /**
@@ -11,7 +11,7 @@ export interface ShotNoiseInputs {
 }
 /**
  * Result of {@link evaluateShotNoise}.
- * @public
+ * @internal
  */
 export interface ShotNoiseResult {
     readonly I_A: number;
@@ -24,7 +24,7 @@ export interface ShotNoiseResult {
  * `e` is the elementary charge. The two-sided convention e I is not
  * this result.
  *
- * @public
+ * @internal
  */
 export declare function evaluateShotNoise({ I_A }: ShotNoiseInputs): ShotNoiseResult;
 //# sourceMappingURL=be85-shot-noise.d.ts.map

package/dist/bridges/be85-shot-noise.js

@@ -20,7 +20,7 @@ import { E_SI } from '../core/constants.js';
  * `e` is the elementary charge. The two-sided convention e I is not
  * this result.
  *
- * @public
+ * @internal
  */
 export function evaluateShotNoise({ I_A }) {
     if (!Number.isFinite(I_A)) {

package/dist/bridges/be86-reynolds-analogy.d.ts

@@ -15,7 +15,7 @@
  */
 /**
  * Inputs for {@link evaluateReynoldsAnalogy}.
- * @public
+ * @internal
  */
 export interface ReynoldsAnalogyInputs {
     /**
@@ -26,7 +26,7 @@ export interface ReynoldsAnalogyInputs {
 }
 /**
  * Result of {@link evaluateReynoldsAnalogy}.
- * @public
+ * @internal
  */
 export interface ReynoldsAnalogyResult {
     readonly C_f: number;
@@ -39,7 +39,7 @@ export interface ReynoldsAnalogyResult {
  * Prandtl number and the wall slopes are hypotheses, not inputs. The
  * function does not accept Pr ≠ 1 and still return this equality.
  *
- * @public
+ * @internal
  */
 export declare function evaluateReynoldsAnalogy({ C_f }: ReynoldsAnalogyInputs): ReynoldsAnalogyResult;
 //# sourceMappingURL=be86-reynolds-analogy.d.ts.map

package/dist/bridges/be86-reynolds-analogy.js

@@ -19,7 +19,7 @@
  * Prandtl number and the wall slopes are hypotheses, not inputs. The
  * function does not accept Pr ≠ 1 and still return this equality.
  *
- * @public
+ * @internal
  */
 export function evaluateReynoldsAnalogy({ C_f }) {
     if (!Number.isFinite(C_f)) {

package/dist/bridges/be87-capacitor-noise.d.ts

@@ -1,6 +1,6 @@
 /**
  * Inputs for {@link evaluateCapacitorNoise}.
- * @public
+ * @internal
  */
 export interface CapacitorNoiseInputs {
     /** Absolute temperature, kelvin. */
@@ -10,7 +10,7 @@ export interface CapacitorNoiseInputs {
 }
 /**
  * Result of {@link evaluateCapacitorNoise}.
- * @public
+ * @internal
  */
 export interface CapacitorNoiseResult {
     readonly T_K: number;
@@ -21,7 +21,7 @@ export interface CapacitorNoiseResult {
 /**
  * Evaluate `⟨v²⟩ = k_B T / C` for one quadratic capacitor term.
  *
- * @public
+ * @internal
  */
 export declare function evaluateCapacitorNoise({ T_K, C_F }: CapacitorNoiseInputs): CapacitorNoiseResult;
 //# sourceMappingURL=be87-capacitor-noise.d.ts.map

package/dist/bridges/be87-capacitor-noise.js

@@ -17,7 +17,7 @@ import { K_B_SI } from '../core/constants.js';
 /**
  * Evaluate `⟨v²⟩ = k_B T / C` for one quadratic capacitor term.
  *
- * @public
+ * @internal
  */
 export function evaluateCapacitorNoise({ T_K, C_F }) {
     if (!Number.isFinite(T_K)) {

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

@@ -0,0 +1,31 @@
+/**
+ * Inputs for {@link evaluateFermiSea}.
+ * @internal
+ */
+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}.
+ * @internal
+ */
+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.
+ *
+ * @internal
+ */
+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.
+ *
+ * @internal
+ */
+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}.
+ * @internal
+ */
+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}.
+ * @internal
+ */
+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}.
+ *
+ * @internal
+ */
+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}.
+ *
+ * @internal
+ */
+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}.
+ * @internal
+ */
+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}.
+ * @internal
+ */
+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.
+ *
+ * @internal
+ */
+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.
+ *
+ * @internal
+ */
+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"}