react-cheminfo 0.2.0 → 0.3.0

package/lib/orbital/core/occupancy.js

@@ -0,0 +1,85 @@
+/**
+ * How the electrons of one atom are laid out over its orbitals, and the two
+ * adjustments a *virtual* orbital needs.
+ *
+ * Split out of `atomicOrbitals.ts` so that module stays what it says it is —
+ * the join of a configuration, a screening, a radial function and a harmonic —
+ * with none of the bookkeeping in the way.
+ */
+import { groupRank } from "./screening.js";
+/** No known element occupies a shell above n = 7. */
+const HIGHEST_SHELL = 7;
+/**
+ * Slater's rules can hand a badly-screened orbital a charge at or below zero —
+ * a 7p on hydrogen, say — which has no bound hydrogen-like solution. Clamping
+ * to a small positive charge keeps such an orbital enormous and diffuse, which
+ * is the honest picture, rather than throwing.
+ */
+export const MINIMUM_CHARGE = 0.3;
+/**
+ * Spread electrons over the orbitals of a subshell by Hund's rule.
+ * @param electrons - Electrons in the subshell.
+ * @param orbitalCount - `2ℓ + 1`.
+ * @returns One count per orbital: singles first, then pairs.
+ */
+export function hundDistribution(electrons, orbitalCount) {
+    const spread = new Array(orbitalCount).fill(0);
+    const capacity = orbitalCount * 2;
+    let left = Math.min(electrons, capacity);
+    for (let pass = 0; pass < 2 && left > 0; pass++) {
+        for (let index = 0; index < orbitalCount && left > 0; index++) {
+            spread[index] = spread[index] + 1;
+            left -= 1;
+        }
+    }
+    return spread;
+}
+/**
+ * How far up the list goes when the caller does not say.
+ * @param valenceShell - Outermost occupied shell.
+ * @returns One shell past it, never below 4 and never above 7.
+ */
+export function defaultMaximumShell(valenceShell) {
+    return Math.min(HIGHEST_SHELL, Math.max(4, valenceShell + 1));
+}
+/**
+ * The configuration with one electron taken out of the outermost occupied
+ * subshell — the atom as it is while one of its electrons sits in a virtual
+ * orbital.
+ * @param configuration - The ground-state configuration to reduce.
+ * @returns The reduced configuration, empty subshells dropped.
+ */
+export function withoutOutermostElectron(configuration) {
+    let outermost = -1;
+    for (let index = 0; index < configuration.length; index++) {
+        const entry = configuration[index];
+        const best = configuration[outermost];
+        if (outermost === -1 ||
+            best === undefined ||
+            groupRank(entry) > groupRank(best)) {
+            outermost = index;
+        }
+    }
+    const reduced = [];
+    for (let index = 0; index < configuration.length; index++) {
+        const entry = configuration[index];
+        const electrons = index === outermost ? entry.electrons - 1 : entry.electrons;
+        if (electrons > 0)
+            reduced.push({ ...entry, electrons });
+    }
+    return reduced;
+}
+/**
+ * The outermost shell an atom actually occupies.
+ * @param configuration - Its configuration.
+ * @returns The largest `n` holding an electron.
+ */
+export function outermostShell(configuration) {
+    let outermost = 1;
+    for (const entry of configuration) {
+        if (entry.electrons > 0 && entry.n > outermost)
+            outermost = entry.n;
+    }
+    return outermost;
+}
+//# sourceMappingURL=occupancy.js.map

package/lib/orbital/core/occupancy.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"occupancy.js","sourceRoot":"","sources":["../../../src/orbital/core/occupancy.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAGH,OAAO,EAAE,SAAS,EAAE,MAAM,gBAAgB,CAAC;AAE3C,qDAAqD;AACrD,MAAM,aAAa,GAAG,CAAC,CAAC;AAExB;;;;;GAKG;AACH,MAAM,CAAC,MAAM,cAAc,GAAG,GAAG,CAAC;AAElC;;;;;GAKG;AACH,MAAM,UAAU,gBAAgB,CAC9B,SAAiB,EACjB,YAAoB;IAEpB,MAAM,MAAM,GAAG,IAAI,KAAK,CAAS,YAAY,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACvD,MAAM,QAAQ,GAAG,YAAY,GAAG,CAAC,CAAC;IAClC,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,QAAQ,CAAC,CAAC;IACzC,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;QAChD,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,YAAY,IAAI,IAAI,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC;YAC9D,MAAM,CAAC,KAAK,CAAC,GAAI,MAAM,CAAC,KAAK,CAAY,GAAG,CAAC,CAAC;YAC9C,IAAI,IAAI,CAAC,CAAC;QACZ,CAAC;IACH,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,mBAAmB,CAAC,YAAoB;IACtD,OAAO,IAAI,CAAC,GAAG,CAAC,aAAa,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,YAAY,GAAG,CAAC,CAAC,CAAC,CAAC;AAChE,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,wBAAwB,CACtC,aAA2C;IAE3C,IAAI,SAAS,GAAG,CAAC,CAAC,CAAC;IACnB,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,aAAa,CAAC,MAAM,EAAE,KAAK,EAAE,EAAE,CAAC;QAC1D,MAAM,KAAK,GAAG,aAAa,CAAC,KAAK,CAAsB,CAAC;QACxD,MAAM,IAAI,GAAG,aAAa,CAAC,SAAS,CAAC,CAAC;QACtC,IACE,SAAS,KAAK,CAAC,CAAC;YAChB,IAAI,KAAK,SAAS;YAClB,SAAS,CAAC,KAAK,CAAC,GAAG,SAAS,CAAC,IAAI,CAAC,EAClC,CAAC;YACD,SAAS,GAAG,KAAK,CAAC;QACpB,CAAC;IACH,CAAC;IACD,MAAM,OAAO,GAAwB,EAAE,CAAC;IACxC,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,aAAa,CAAC,MAAM,EAAE,KAAK,EAAE,EAAE,CAAC;QAC1D,MAAM,KAAK,GAAG,aAAa,CAAC,KAAK,CAAsB,CAAC;QACxD,MAAM,SAAS,GACb,KAAK,KAAK,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,SAAS,CAAC;QAC9D,IAAI,SAAS,GAAG,CAAC;YAAE,OAAO,CAAC,IAAI,CAAC,EAAE,GAAG,KAAK,EAAE,SAAS,EAAE,CAAC,CAAC;IAC3D,CAAC;IACD,OAAO,OAAO,CAAC;AACjB,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,cAAc,CAC5B,aAA2C;IAE3C,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,KAAK,MAAM,KAAK,IAAI,aAAa,EAAE,CAAC;QAClC,IAAI,KAAK,CAAC,SAAS,GAAG,CAAC,IAAI,KAAK,CAAC,CAAC,GAAG,SAAS;YAAE,SAAS,GAAG,KAAK,CAAC,CAAC,CAAC;IACtE,CAAC;IACD,OAAO,SAAS,CAAC;AACnB,CAAC"}

package/lib/orbital/core/palette.d.ts

@@ -0,0 +1,36 @@
+/**
+ * Colours for orbital phases and nodal surfaces.
+ *
+ * Plain data: the renderer picks a palette, this module says what the hex
+ * values are. Two phase palettes ship, and both are first-class — around 8% of
+ * male students have a red/green or red/blue-at-low-luminance deficiency, so
+ * the blue/yellow pair is a supported choice, not a fallback.
+ */
+/** Which pair of phase colours the viewer uses. */
+export type PhasePaletteId = 'textbook' | 'colourBlindSafe';
+/** The colours one phase palette provides. */
+export interface PhasePalette {
+    id: PhasePaletteId;
+    /** Name shown in the palette picker. */
+    label: string;
+    /** Lobe where the wavefunction is positive. */
+    positive: string;
+    /** Lobe where the wavefunction is negative. */
+    negative: string;
+    /** Lobe drawn with the `phase` display flag off. */
+    neutral: string;
+    /** Translucent disc marking a nodal plane. */
+    node: string;
+}
+/**
+ * The two phase palettes.
+ *
+ * `textbook` is the blue/red of printed MO diagrams and of Jmol's signed
+ * isosurfaces; `colourBlindSafe` swaps the negative lobe for the amber of the
+ * Okabe-Ito qualitative palette, which stays separable under protanopia,
+ * deuteranopia and tritanopia.
+ */
+export declare const PHASE_PALETTES: Record<PhasePaletteId, PhasePalette>;
+/** Palette used until the student picks another one. */
+export declare const DEFAULT_PHASE_PALETTE_ID: PhasePaletteId;
+//# sourceMappingURL=palette.d.ts.map

package/lib/orbital/core/palette.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"palette.d.ts","sourceRoot":"","sources":["../../../src/orbital/core/palette.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAEH,mDAAmD;AACnD,MAAM,MAAM,cAAc,GAAG,UAAU,GAAG,iBAAiB,CAAC;AAE5D,8CAA8C;AAC9C,MAAM,WAAW,YAAY;IAC3B,EAAE,EAAE,cAAc,CAAC;IACnB,wCAAwC;IACxC,KAAK,EAAE,MAAM,CAAC;IACd,+CAA+C;IAC/C,QAAQ,EAAE,MAAM,CAAC;IACjB,+CAA+C;IAC/C,QAAQ,EAAE,MAAM,CAAC;IACjB,oDAAoD;IACpD,OAAO,EAAE,MAAM,CAAC;IAChB,8CAA8C;IAC9C,IAAI,EAAE,MAAM,CAAC;CACd;AAED;;;;;;;GAOG;AACH,eAAO,MAAM,cAAc,EAAE,MAAM,CAAC,cAAc,EAAE,YAAY,CAiB/D,CAAC;AAEF,wDAAwD;AACxD,eAAO,MAAM,wBAAwB,EAAE,cAA2B,CAAC"}

package/lib/orbital/core/palette.js

@@ -0,0 +1,37 @@
+/**
+ * Colours for orbital phases and nodal surfaces.
+ *
+ * Plain data: the renderer picks a palette, this module says what the hex
+ * values are. Two phase palettes ship, and both are first-class — around 8% of
+ * male students have a red/green or red/blue-at-low-luminance deficiency, so
+ * the blue/yellow pair is a supported choice, not a fallback.
+ */
+/**
+ * The two phase palettes.
+ *
+ * `textbook` is the blue/red of printed MO diagrams and of Jmol's signed
+ * isosurfaces; `colourBlindSafe` swaps the negative lobe for the amber of the
+ * Okabe-Ito qualitative palette, which stays separable under protanopia,
+ * deuteranopia and tritanopia.
+ */
+export const PHASE_PALETTES = {
+    textbook: {
+        id: 'textbook',
+        label: 'Blue / red (textbook)',
+        positive: '#2563eb',
+        negative: '#dc2626',
+        neutral: '#9ca3af',
+        node: '#64748b',
+    },
+    colourBlindSafe: {
+        id: 'colourBlindSafe',
+        label: 'Blue / yellow (colour-blind safe)',
+        positive: '#0072b2',
+        negative: '#e69f00',
+        neutral: '#9ca3af',
+        node: '#56606b',
+    },
+};
+/** Palette used until the student picks another one. */
+export const DEFAULT_PHASE_PALETTE_ID = 'textbook';
+//# sourceMappingURL=palette.js.map

package/lib/orbital/core/palette.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"palette.js","sourceRoot":"","sources":["../../../src/orbital/core/palette.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AAoBH;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,cAAc,GAAyC;IAClE,QAAQ,EAAE;QACR,EAAE,EAAE,UAAU;QACd,KAAK,EAAE,uBAAuB;QAC9B,QAAQ,EAAE,SAAS;QACnB,QAAQ,EAAE,SAAS;QACnB,OAAO,EAAE,SAAS;QAClB,IAAI,EAAE,SAAS;KAChB;IACD,eAAe,EAAE;QACf,EAAE,EAAE,iBAAiB;QACrB,KAAK,EAAE,mCAAmC;QAC1C,QAAQ,EAAE,SAAS;QACnB,QAAQ,EAAE,SAAS;QACnB,OAAO,EAAE,SAAS;QAClB,IAAI,EAAE,SAAS;KAChB;CACF,CAAC;AAEF,wDAAwD;AACxD,MAAM,CAAC,MAAM,wBAAwB,GAAmB,UAAU,CAAC"}

package/lib/orbital/core/realHarmonics.d.ts

@@ -0,0 +1,57 @@
+/**
+ * The real spherical harmonics, ℓ = 0 to 3, written as cartesian ratios.
+ *
+ * These are the angular halves of the orbitals a chemist draws: the real
+ * combinations `p_x, p_y, p_z` rather than the complex `m = −1, 0, +1`, because
+ * only the real ones have the lobes along axes that a bond can point at.
+ *
+ * Each is normalised so that `∫|Y|² dΩ = 1` over the unit sphere, and each is
+ * given `(x, y, z)` in ångström together with their length, so sampling a grid
+ * never calls `acos` or `atan2`. The number of angular nodes is exactly ℓ, and
+ * this file is where that fact becomes visible: `s` has none, every `p` has one
+ * plane, every `d` two surfaces, every `f` three.
+ *
+ * `key` is ascii and appears in a URL — `#/elements/Fe/3dx2-y2`; `label` is the
+ * typeset subscript a student reads.
+ */
+/** Subshell letter for an angular momentum, `s p d f`. */
+export declare const SUBSHELL_LETTERS: readonly ["s", "p", "d", "f"];
+/** One real spherical harmonic. */
+export interface RealHarmonic {
+    /** Angular momentum quantum number, also the count of angular nodes. */
+    l: number;
+    /** Ascii subscript, url-safe: `x`, `z2`, `x2-y2`, `y_3x2-y2`. */
+    key: string;
+    /** Typeset subscript, e.g. `x²−y²`. Empty for `s`. */
+    label: string;
+    /**
+     * Value of the harmonic at a point.
+     * @param x - Offset from the nucleus along x, ångström.
+     * @param y - Offset from the nucleus along y, ångström.
+     * @param z - Offset from the nucleus along z, ångström.
+     * @param radius - Length of `(x, y, z)`, passed in because the caller has
+     * already computed it for the radial part.
+     * @returns The dimensionless amplitude; 0 at the nucleus for ℓ > 0.
+     */
+    evaluate: (x: number, y: number, z: number, radius: number) => number;
+}
+/**
+ * The harmonics of one shell, indexed by ℓ: `REAL_HARMONICS[2]` is the five d
+ * functions. The order inside a shell is the one the orbital list shows.
+ */
+export declare const REAL_HARMONICS: RealHarmonic[][];
+/**
+ * The harmonics of one subshell.
+ * @param l - Angular momentum quantum number, 0 to 3.
+ * @returns Its `2ℓ + 1` real harmonics, in list order.
+ * @throws {Error} When ℓ is outside the s, p, d and f shells this site draws.
+ */
+export declare function harmonicsOf(l: number): RealHarmonic[];
+/**
+ * The letter naming a subshell.
+ * @param l - Angular momentum quantum number, 0 to 3.
+ * @returns `s`, `p`, `d` or `f`.
+ * @throws {Error} When ℓ is outside them.
+ */
+export declare function subshellLetter(l: number): string;
+//# sourceMappingURL=realHarmonics.d.ts.map

package/lib/orbital/core/realHarmonics.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"realHarmonics.d.ts","sourceRoot":"","sources":["../../../src/orbital/core/realHarmonics.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,0DAA0D;AAC1D,eAAO,MAAM,gBAAgB,+BAAgC,CAAC;AAE9D,mCAAmC;AACnC,MAAM,WAAW,YAAY;IAC3B,wEAAwE;IACxE,CAAC,EAAE,MAAM,CAAC;IACV,iEAAiE;IACjE,GAAG,EAAE,MAAM,CAAC;IACZ,sDAAsD;IACtD,KAAK,EAAE,MAAM,CAAC;IACd;;;;;;;;OAQG;IACH,QAAQ,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,KAAK,MAAM,CAAC;CACvE;AAeD;;;GAGG;AACH,eAAO,MAAM,cAAc,EAAE,YAAY,EAAE,EAsF1C,CAAC;AAEF;;;;;GAKG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,MAAM,GAAG,YAAY,EAAE,CAMrD;AAED;;;;;GAKG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAMhD"}

package/lib/orbital/core/realHarmonics.js

@@ -0,0 +1,144 @@
+/**
+ * The real spherical harmonics, ℓ = 0 to 3, written as cartesian ratios.
+ *
+ * These are the angular halves of the orbitals a chemist draws: the real
+ * combinations `p_x, p_y, p_z` rather than the complex `m = −1, 0, +1`, because
+ * only the real ones have the lobes along axes that a bond can point at.
+ *
+ * Each is normalised so that `∫|Y|² dΩ = 1` over the unit sphere, and each is
+ * given `(x, y, z)` in ångström together with their length, so sampling a grid
+ * never calls `acos` or `atan2`. The number of angular nodes is exactly ℓ, and
+ * this file is where that fact becomes visible: `s` has none, every `p` has one
+ * plane, every `d` two surfaces, every `f` three.
+ *
+ * `key` is ascii and appears in a URL — `#/elements/Fe/3dx2-y2`; `label` is the
+ * typeset subscript a student reads.
+ */
+/** Subshell letter for an angular momentum, `s p d f`. */
+export const SUBSHELL_LETTERS = ['s', 'p', 'd', 'f'];
+const FOUR_PI = 4 * Math.PI;
+const S = Math.sqrt(1 / FOUR_PI);
+const P = Math.sqrt(3 / FOUR_PI);
+const D_Z2 = Math.sqrt(5 / (16 * Math.PI));
+const D_PLANAR = Math.sqrt(15 / FOUR_PI);
+const D_X2Y2 = Math.sqrt(15 / (16 * Math.PI));
+const F_Z3 = Math.sqrt(7 / (16 * Math.PI));
+const F_XZ2 = Math.sqrt(21 / (32 * Math.PI));
+const F_ZX2Y2 = Math.sqrt(105 / (16 * Math.PI));
+const F_XYZ = Math.sqrt(105 / FOUR_PI);
+const F_CUBIC = Math.sqrt(35 / (32 * Math.PI));
+/**
+ * The harmonics of one shell, indexed by ℓ: `REAL_HARMONICS[2]` is the five d
+ * functions. The order inside a shell is the one the orbital list shows.
+ */
+export const REAL_HARMONICS = [
+    [{ l: 0, key: '', label: '', evaluate: () => S }],
+    [
+        { l: 1, key: 'x', label: 'x', evaluate: (x, y, z, r) => (P * x) / r },
+        { l: 1, key: 'y', label: 'y', evaluate: (x, y, z, r) => (P * y) / r },
+        { l: 1, key: 'z', label: 'z', evaluate: (x, y, z, r) => (P * z) / r },
+    ],
+    [
+        {
+            l: 2,
+            key: 'z2',
+            label: 'z²',
+            evaluate: (x, y, z, r) => (D_Z2 * (3 * z * z - r * r)) / (r * r),
+        },
+        {
+            l: 2,
+            key: 'xz',
+            label: 'xz',
+            evaluate: (x, y, z, r) => (D_PLANAR * x * z) / (r * r),
+        },
+        {
+            l: 2,
+            key: 'yz',
+            label: 'yz',
+            evaluate: (x, y, z, r) => (D_PLANAR * y * z) / (r * r),
+        },
+        {
+            l: 2,
+            key: 'xy',
+            label: 'xy',
+            evaluate: (x, y, z, r) => (D_PLANAR * x * y) / (r * r),
+        },
+        {
+            l: 2,
+            key: 'x2-y2',
+            label: 'x²−y²',
+            evaluate: (x, y, z, r) => (D_X2Y2 * (x * x - y * y)) / (r * r),
+        },
+    ],
+    [
+        {
+            l: 3,
+            key: 'z3',
+            label: 'z³',
+            evaluate: (x, y, z, r) => (F_Z3 * z * (5 * z * z - 3 * r * r)) / (r * r * r),
+        },
+        {
+            l: 3,
+            key: 'xz2',
+            label: 'xz²',
+            evaluate: (x, y, z, r) => (F_XZ2 * x * (5 * z * z - r * r)) / (r * r * r),
+        },
+        {
+            l: 3,
+            key: 'yz2',
+            label: 'yz²',
+            evaluate: (x, y, z, r) => (F_XZ2 * y * (5 * z * z - r * r)) / (r * r * r),
+        },
+        {
+            l: 3,
+            key: 'xyz',
+            label: 'xyz',
+            evaluate: (x, y, z, r) => (F_XYZ * x * y * z) / (r * r * r),
+        },
+        {
+            l: 3,
+            key: 'z_x2-y2',
+            label: 'z(x²−y²)',
+            evaluate: (x, y, z, r) => (F_ZX2Y2 * z * (x * x - y * y)) / (r * r * r),
+        },
+        {
+            l: 3,
+            key: 'x_x2-3y2',
+            label: 'x(x²−3y²)',
+            evaluate: (x, y, z, r) => (F_CUBIC * x * (x * x - 3 * y * y)) / (r * r * r),
+        },
+        {
+            l: 3,
+            key: 'y_3x2-y2',
+            label: 'y(3x²−y²)',
+            evaluate: (x, y, z, r) => (F_CUBIC * y * (3 * x * x - y * y)) / (r * r * r),
+        },
+    ],
+];
+/**
+ * The harmonics of one subshell.
+ * @param l - Angular momentum quantum number, 0 to 3.
+ * @returns Its `2ℓ + 1` real harmonics, in list order.
+ * @throws {Error} When ℓ is outside the s, p, d and f shells this site draws.
+ */
+export function harmonicsOf(l) {
+    const shell = REAL_HARMONICS[l];
+    if (shell === undefined) {
+        throw new RangeError(`no real harmonics for angular momentum ${l}`);
+    }
+    return shell;
+}
+/**
+ * The letter naming a subshell.
+ * @param l - Angular momentum quantum number, 0 to 3.
+ * @returns `s`, `p`, `d` or `f`.
+ * @throws {Error} When ℓ is outside them.
+ */
+export function subshellLetter(l) {
+    const letter = SUBSHELL_LETTERS[l];
+    if (letter === undefined) {
+        throw new RangeError(`no subshell letter for angular momentum ${l}`);
+    }
+    return letter;
+}
+//# sourceMappingURL=realHarmonics.js.map

package/lib/orbital/core/realHarmonics.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"realHarmonics.js","sourceRoot":"","sources":["../../../src/orbital/core/realHarmonics.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,0DAA0D;AAC1D,MAAM,CAAC,MAAM,gBAAgB,GAAG,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAU,CAAC;AAsB9D,MAAM,OAAO,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;AAE5B,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC;AACjC,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC;AACjC,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;AAC3C,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,OAAO,CAAC,CAAC;AACzC,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;AAC9C,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;AAC3C,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;AAC7C,MAAM,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;AAChD,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,OAAO,CAAC,CAAC;AACvC,MAAM,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;AAE/C;;;GAGG;AACH,MAAM,CAAC,MAAM,cAAc,GAAqB;IAC9C,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC;IACjD;QACE,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE;QACrE,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE;QACrE,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE;KACtE;IACD;QACE;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,IAAI;YACT,KAAK,EAAE,IAAI;YACX,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;SACjE;QACD;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,IAAI;YACT,KAAK,EAAE,IAAI;YACX,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;SACvD;QACD;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,IAAI;YACT,KAAK,EAAE,IAAI;YACX,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;SACvD;QACD;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,IAAI;YACT,KAAK,EAAE,IAAI;YACX,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;SACvD;QACD;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,OAAO;YACZ,KAAK,EAAE,OAAO;YACd,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;SAC/D;KACF;IACD;QACE;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,IAAI;YACT,KAAK,EAAE,IAAI;YACX,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CACvB,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;SACrD;QACD;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,KAAK;YACV,KAAK,EAAE,KAAK;YACZ,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;SAC1E;QACD;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,KAAK;YACV,KAAK,EAAE,KAAK;YACZ,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;SAC1E;QACD;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,KAAK;YACV,KAAK,EAAE,KAAK;YACZ,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;SAC5D;QACD;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,SAAS;YACd,KAAK,EAAE,UAAU;YACjB,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;SACxE;QACD;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,UAAU;YACf,KAAK,EAAE,WAAW;YAClB,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CACvB,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;SACpD;QACD;YACE,CAAC,EAAE,CAAC;YACJ,GAAG,EAAE,UAAU;YACf,KAAK,EAAE,WAAW;YAClB,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CACvB,CAAC,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;SACpD;KACF;CACF,CAAC;AAEF;;;;;GAKG;AACH,MAAM,UAAU,WAAW,CAAC,CAAS;IACnC,MAAM,KAAK,GAAG,cAAc,CAAC,CAAC,CAAC,CAAC;IAChC,IAAI,KAAK,KAAK,SAAS,EAAE,CAAC;QACxB,MAAM,IAAI,UAAU,CAAC,0CAA0C,CAAC,EAAE,CAAC,CAAC;IACtE,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,cAAc,CAAC,CAAS;IACtC,MAAM,MAAM,GAAG,gBAAgB,CAAC,CAAC,CAAC,CAAC;IACnC,IAAI,MAAM,KAAK,SAAS,EAAE,CAAC;QACzB,MAAM,IAAI,UAAU,CAAC,2CAA2C,CAAC,EAAE,CAAC,CAAC;IACvE,CAAC;IACD,OAAO,MAAM,CAAC;AAChB,CAAC"}

package/lib/orbital/core/sample.d.ts

@@ -0,0 +1,52 @@
+/**
+ * The one job the viewer asks for, shaped so it can cross a worker boundary.
+ *
+ * The request carries an atomic number and an orbital id, never an
+ * `AtomicOrbital`: that object holds its spherical harmonic as a
+ * function*, which the structured clone algorithm cannot copy. Rebuilding the
+ * list from two plain values costs microseconds and keeps the boundary honest.
+ *
+ * A site that would rather not spend ~25 ms of its main thread per orbital
+ * implements this contract in a worker and hands `AtomicOrbitalViewer` the
+ * `sample` prop; one that does not gets the in-process version for free.
+ */
+import type { OrbitalGrid } from './grid.ts';
+/** Which orbital to sample, and how finely. */
+export interface AtomicSampleRequest {
+    /** Proton count, 1 to 118. */
+    atomicNumber: number;
+    /** Which orbital of that element, e.g. `3dz2`. */
+    orbitalId: string;
+    /**
+     * Samples along each edge of the cube. The cost is the cube of this.
+     * @default 56
+     */
+    resolution?: number;
+}
+/** What one sampling job produces. */
+export interface AtomicSampleResult {
+    /** The signed field, ready for an isosurface. */
+    grid: OrbitalGrid;
+    /** Radii of the orbital's radial nodes, ångström, ascending. */
+    nodeRadii: number[];
+}
+/**
+ * Sample one atomic orbital onto a grid.
+ *
+ * Synchronous and free of React, molstar and the DOM, so it runs unchanged on
+ * the main thread, in a worker, and under vitest.
+ * @param request - See {@link AtomicSampleRequest}.
+ * @returns The field, and where its radial nodes are.
+ * @throws {Error} When the element or the orbital id names nothing.
+ */
+export declare function runAtomicSample(request: AtomicSampleRequest): AtomicSampleResult;
+/** How a caller supplies the sampling — in process, or through a worker. */
+export type AtomicSampler = (request: AtomicSampleRequest) => Promise<AtomicSampleResult>;
+/**
+ * The default sampler: {@link runAtomicSample}, yielded to the event loop so
+ * the spinner has painted before the main thread is taken.
+ * @param request - See {@link AtomicSampleRequest}.
+ * @returns The field, and where its radial nodes are.
+ */
+export declare const sampleInProcess: AtomicSampler;
+//# sourceMappingURL=sample.d.ts.map

package/lib/orbital/core/sample.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"sample.d.ts","sourceRoot":"","sources":["../../../src/orbital/core/sample.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AASH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,WAAW,CAAC;AAG7C,+CAA+C;AAC/C,MAAM,WAAW,mBAAmB;IAClC,8BAA8B;IAC9B,YAAY,EAAE,MAAM,CAAC;IACrB,kDAAkD;IAClD,SAAS,EAAE,MAAM,CAAC;IAClB;;;OAGG;IACH,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB;AAED,sCAAsC;AACtC,MAAM,WAAW,kBAAkB;IACjC,iDAAiD;IACjD,IAAI,EAAE,WAAW,CAAC;IAClB,gEAAgE;IAChE,SAAS,EAAE,MAAM,EAAE,CAAC;CACrB;AAED;;;;;;;;GAQG;AACH,wBAAgB,eAAe,CAC7B,OAAO,EAAE,mBAAmB,GAC3B,kBAAkB,CAiBpB;AAED,4EAA4E;AAC5E,MAAM,MAAM,aAAa,GAAG,CAC1B,OAAO,EAAE,mBAAmB,KACzB,OAAO,CAAC,kBAAkB,CAAC,CAAC;AAEjC;;;;;GAKG;AACH,eAAO,MAAM,eAAe,EAAE,aAK7B,CAAC"}

package/lib/orbital/core/sample.js

@@ -0,0 +1,54 @@
+/**
+ * The one job the viewer asks for, shaped so it can cross a worker boundary.
+ *
+ * The request carries an atomic number and an orbital id, never an
+ * `AtomicOrbital`: that object holds its spherical harmonic as a
+ * function*, which the structured clone algorithm cannot copy. Rebuilding the
+ * list from two plain values costs microseconds and keeps the boundary honest.
+ *
+ * A site that would rather not spend ~25 ms of its main thread per orbital
+ * implements this contract in a worker and hands `AtomicOrbitalViewer` the
+ * `sample` prop; one that does not gets the in-process version for free.
+ */
+import { sampleAtomicOrbital } from "./atomicGrid.js";
+import { atomicOrbitalsOf, findAtomicOrbital, hydrogenicParametersOf, } from "./atomicOrbitals.js";
+import { radialNodeRadii } from "./hydrogenic.js";
+/**
+ * Sample one atomic orbital onto a grid.
+ *
+ * Synchronous and free of React, molstar and the DOM, so it runs unchanged on
+ * the main thread, in a worker, and under vitest.
+ * @param request - See {@link AtomicSampleRequest}.
+ * @returns The field, and where its radial nodes are.
+ * @throws {Error} When the element or the orbital id names nothing.
+ */
+export function runAtomicSample(request) {
+    const { atomicNumber, orbitalId, resolution } = request;
+    const orbitals = atomicOrbitalsOf(atomicNumber);
+    const orbital = findAtomicOrbital(orbitals, orbitalId);
+    if (orbital === null) {
+        throw new Error(`Element ${atomicNumber} has no orbital ${orbitalId}. Known: ${orbitals
+            .map((entry) => entry.id)
+            .join(', ')}.`);
+    }
+    const options = {};
+    if (resolution !== undefined)
+        options.resolution = resolution;
+    return {
+        grid: sampleAtomicOrbital(orbital, options),
+        nodeRadii: radialNodeRadii(hydrogenicParametersOf(orbital)),
+    };
+}
+/**
+ * The default sampler: {@link runAtomicSample}, yielded to the event loop so
+ * the spinner has painted before the main thread is taken.
+ * @param request - See {@link AtomicSampleRequest}.
+ * @returns The field, and where its radial nodes are.
+ */
+export const sampleInProcess = async (request) => {
+    await new Promise((resolve) => {
+        setTimeout(resolve, 0);
+    });
+    return runAtomicSample(request);
+};
+//# sourceMappingURL=sample.js.map

package/lib/orbital/core/sample.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"sample.js","sourceRoot":"","sources":["../../../src/orbital/core/sample.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAGH,OAAO,EAAE,mBAAmB,EAAE,MAAM,iBAAiB,CAAC;AACtD,OAAO,EACL,gBAAgB,EAChB,iBAAiB,EACjB,sBAAsB,GACvB,MAAM,qBAAqB,CAAC;AAE7B,OAAO,EAAE,eAAe,EAAE,MAAM,iBAAiB,CAAC;AAuBlD;;;;;;;;GAQG;AACH,MAAM,UAAU,eAAe,CAC7B,OAA4B;IAE5B,MAAM,EAAE,YAAY,EAAE,SAAS,EAAE,UAAU,EAAE,GAAG,OAAO,CAAC;IACxD,MAAM,QAAQ,GAAG,gBAAgB,CAAC,YAAY,CAAC,CAAC;IAChD,MAAM,OAAO,GAAG,iBAAiB,CAAC,QAAQ,EAAE,SAAS,CAAC,CAAC;IACvD,IAAI,OAAO,KAAK,IAAI,EAAE,CAAC;QACrB,MAAM,IAAI,KAAK,CACb,WAAW,YAAY,mBAAmB,SAAS,YAAY,QAAQ;aACpE,GAAG,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC;aACxB,IAAI,CAAC,IAAI,CAAC,GAAG,CACjB,CAAC;IACJ,CAAC;IACD,MAAM,OAAO,GAAsB,EAAE,CAAC;IACtC,IAAI,UAAU,KAAK,SAAS;QAAE,OAAO,CAAC,UAAU,GAAG,UAAU,CAAC;IAC9D,OAAO;QACL,IAAI,EAAE,mBAAmB,CAAC,OAAO,EAAE,OAAO,CAAC;QAC3C,SAAS,EAAE,eAAe,CAAC,sBAAsB,CAAC,OAAO,CAAC,CAAC;KAC5D,CAAC;AACJ,CAAC;AAOD;;;;;GAKG;AACH,MAAM,CAAC,MAAM,eAAe,GAAkB,KAAK,EAAE,OAAO,EAAE,EAAE;IAC9D,MAAM,IAAI,OAAO,CAAC,CAAC,OAAO,EAAE,EAAE;QAC5B,UAAU,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;IACzB,CAAC,CAAC,CAAC;IACH,OAAO,eAAe,CAAC,OAAO,CAAC,CAAC;AAClC,CAAC,CAAC"}

package/lib/orbital/core/screening.d.ts

@@ -0,0 +1,46 @@
+/**
+ * Slater's rules: how much of the nucleus an electron actually feels.
+ *
+ * A hydrogen-like orbital needs one number, the effective nuclear charge, and
+ * taking the full `Z` would draw carbon's 2p six times too small. Slater's rules
+ * estimate the shielding `S` from the other electrons, and `Z_eff = Z − S` is
+ * what the radial function is built with.
+ *
+ * They are an empirical fit from 1930, not a calculation — but they are the ones
+ * every general-chemistry course teaches, they reproduce the periodic trends
+ * that matter (Z_eff rises across a period, barely moves down a group), and a
+ * student can check the arithmetic here against the one they did by hand.
+ *
+ * The electrons are grouped as Slater grouped them, with s and p sharing a
+ * group and each d and f standing alone:
+ *
+ * ```
+ * [1s] [2s2p] [3s3p] [3d] [4s4p] [4d] [4f] [5s5p] [5d] [5f] [6s6p] [6d] [7s7p]
+ * ```
+ */
+import type { Subshell, SubshellOccupancy } from './electronConfiguration.ts';
+/** The shielding an orbital feels, and what is left of the nuclear charge. */
+export interface Screening {
+    /** Sum of the shielding contributions, `S`. */
+    shielding: number;
+    /** `Z − S`, the charge the hydrogen-like radial function is built with. */
+    effectiveCharge: number;
+}
+/**
+ * Apply Slater's rules to one subshell of one atom.
+ * @param atomicNumber - Proton count, the `Z` the shielding is subtracted from.
+ * @param configuration - Every occupied subshell of the atom.
+ * @param subshell - The subshell whose electron is being screened.
+ * @returns Its shielding and effective nuclear charge.
+ */
+export declare function slaterScreening(atomicNumber: number, configuration: readonly SubshellOccupancy[], subshell: Subshell): Screening;
+/**
+ * Where a subshell sits in Slater's left-to-right ordering of groups.
+ *
+ * s and p of the same shell share a group, so they share a rank; d and f each
+ * get their own, placed after the s/p group of the same `n`.
+ * @param subshell - The subshell to rank.
+ * @returns A comparable rank; equal ranks mean the same Slater group.
+ */
+export declare function groupRank(subshell: Subshell): number;
+//# sourceMappingURL=screening.d.ts.map

package/lib/orbital/core/screening.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"screening.d.ts","sourceRoot":"","sources":["../../../src/orbital/core/screening.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;GAmBG;AAEH,OAAO,KAAK,EAAE,QAAQ,EAAE,iBAAiB,EAAE,MAAM,4BAA4B,CAAC;AAE9E,8EAA8E;AAC9E,MAAM,WAAW,SAAS;IACxB,+CAA+C;IAC/C,SAAS,EAAE,MAAM,CAAC;IAClB,2EAA2E;IAC3E,eAAe,EAAE,MAAM,CAAC;CACzB;AAED;;;;;;GAMG;AACH,wBAAgB,eAAe,CAC7B,YAAY,EAAE,MAAM,EACpB,aAAa,EAAE,SAAS,iBAAiB,EAAE,EAC3C,QAAQ,EAAE,QAAQ,GACjB,SAAS,CA+BX;AAED;;;;;;;GAOG;AACH,wBAAgB,SAAS,CAAC,QAAQ,EAAE,QAAQ,GAAG,MAAM,CAEpD"}

package/lib/orbital/core/screening.js

@@ -0,0 +1,72 @@
+/**
+ * Slater's rules: how much of the nucleus an electron actually feels.
+ *
+ * A hydrogen-like orbital needs one number, the effective nuclear charge, and
+ * taking the full `Z` would draw carbon's 2p six times too small. Slater's rules
+ * estimate the shielding `S` from the other electrons, and `Z_eff = Z − S` is
+ * what the radial function is built with.
+ *
+ * They are an empirical fit from 1930, not a calculation — but they are the ones
+ * every general-chemistry course teaches, they reproduce the periodic trends
+ * that matter (Z_eff rises across a period, barely moves down a group), and a
+ * student can check the arithmetic here against the one they did by hand.
+ *
+ * The electrons are grouped as Slater grouped them, with s and p sharing a
+ * group and each d and f standing alone:
+ *
+ * ```
+ * [1s] [2s2p] [3s3p] [3d] [4s4p] [4d] [4f] [5s5p] [5d] [5f] [6s6p] [6d] [7s7p]
+ * ```
+ */
+/**
+ * Apply Slater's rules to one subshell of one atom.
+ * @param atomicNumber - Proton count, the `Z` the shielding is subtracted from.
+ * @param configuration - Every occupied subshell of the atom.
+ * @param subshell - The subshell whose electron is being screened.
+ * @returns Its shielding and effective nuclear charge.
+ */
+export function slaterScreening(atomicNumber, configuration, subshell) {
+    const rank = groupRank(subshell);
+    const isCore = subshell.n === 1;
+    const isDiffuse = subshell.l >= 2;
+    let shielding = 0;
+    for (const occupied of configuration) {
+        const otherRank = groupRank(occupied);
+        if (otherRank > rank)
+            continue;
+        if (otherRank === rank) {
+            // An electron does not shield itself, so one is left out — but only from
+            // its *own* subshell. 2s and 2p share a group, and every one of carbon's
+            // two 2s electrons shields its 2p electron in full.
+            const others = Math.max(0, occupied.electrons - (sameSubshell(occupied, subshell) ? 1 : 0));
+            shielding += others * (isCore ? 0.3 : 0.35);
+            continue;
+        }
+        if (isDiffuse) {
+            // A d or f electron is screened completely by everything inside it.
+            shielding += occupied.electrons;
+            continue;
+        }
+        shielding +=
+            occupied.electrons * (occupied.n === subshell.n - 1 ? 0.85 : 1);
+    }
+    return {
+        shielding,
+        effectiveCharge: atomicNumber - shielding,
+    };
+}
+/**
+ * Where a subshell sits in Slater's left-to-right ordering of groups.
+ *
+ * s and p of the same shell share a group, so they share a rank; d and f each
+ * get their own, placed after the s/p group of the same `n`.
+ * @param subshell - The subshell to rank.
+ * @returns A comparable rank; equal ranks mean the same Slater group.
+ */
+export function groupRank(subshell) {
+    return subshell.n * 10 + (subshell.l <= 1 ? 0 : subshell.l);
+}
+function sameSubshell(first, second) {
+    return first.n === second.n && first.l === second.l;
+}
+//# sourceMappingURL=screening.js.map

package/lib/orbital/core/screening.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"screening.js","sourceRoot":"","sources":["../../../src/orbital/core/screening.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;GAmBG;AAYH;;;;;;GAMG;AACH,MAAM,UAAU,eAAe,CAC7B,YAAoB,EACpB,aAA2C,EAC3C,QAAkB;IAElB,MAAM,IAAI,GAAG,SAAS,CAAC,QAAQ,CAAC,CAAC;IACjC,MAAM,MAAM,GAAG,QAAQ,CAAC,CAAC,KAAK,CAAC,CAAC;IAChC,MAAM,SAAS,GAAG,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC;IAClC,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,KAAK,MAAM,QAAQ,IAAI,aAAa,EAAE,CAAC;QACrC,MAAM,SAAS,GAAG,SAAS,CAAC,QAAQ,CAAC,CAAC;QACtC,IAAI,SAAS,GAAG,IAAI;YAAE,SAAS;QAC/B,IAAI,SAAS,KAAK,IAAI,EAAE,CAAC;YACvB,yEAAyE;YACzE,yEAAyE;YACzE,oDAAoD;YACpD,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CACrB,CAAC,EACD,QAAQ,CAAC,SAAS,GAAG,CAAC,YAAY,CAAC,QAAQ,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAChE,CAAC;YACF,SAAS,IAAI,MAAM,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;YAC5C,SAAS;QACX,CAAC;QACD,IAAI,SAAS,EAAE,CAAC;YACd,oEAAoE;YACpE,SAAS,IAAI,QAAQ,CAAC,SAAS,CAAC;YAChC,SAAS;QACX,CAAC;QACD,SAAS;YACP,QAAQ,CAAC,SAAS,GAAG,CAAC,QAAQ,CAAC,CAAC,KAAK,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACpE,CAAC;IACD,OAAO;QACL,SAAS;QACT,eAAe,EAAE,YAAY,GAAG,SAAS;KAC1C,CAAC;AACJ,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,SAAS,CAAC,QAAkB;IAC1C,OAAO,QAAQ,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,QAAQ,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;AAC9D,CAAC;AAED,SAAS,YAAY,CAAC,KAAe,EAAE,MAAgB;IACrD,OAAO,KAAK,CAAC,CAAC,KAAK,MAAM,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,MAAM,CAAC,CAAC,CAAC;AACtD,CAAC"}