@woosh/meep-engine 2.130.0 → 2.131.1

package/src/core/math/physics/mie/ri_pollen.js

@@ -0,0 +1,146 @@
+import {
+    ParameterLookupTable
+} from "../../../../engine/graphics/particles/particular/engine/parameter/ParameterLookupTable.js";
+
+/**
+ * Complex refractive index of pollen (sporopollenin-rich organic shell): n + i k, wavelength in nm.
+ * Assumptions: dry pollen (no internal water), ~300 K, bulk-effective constants; morphology not modeled.
+ * Notes and composite sources/guidance:
+ * - Pollen exine composed largely of sporopollenin; visible n typically ~1.53–1.57; very small k in visible, stronger in UV.
+ * - Organic functional groups show bands: strong C–H/O–H stretch complex near 3.2–3.5 µm, carbonyl/other modes near 5.8–6.2 µm,
+ *   and features around 6.8–7.3 µm; far-IR tails are modest compared to silicates.
+ * - This LUT is designed for radiometric/Mie use over 200–20,000 nm with linear interpolation. Sampling is denser near bands.
+ * - Composite guidance: polymer/biopolymer optics (e.g., Querry optical constants compendium), aerosol brown-organics literature,
+ *   and typical measured visible n for biological tissues/pollen. Values are smooth and representative, not species-specific.
+ */
+const POLLEN_REFRACTION_INDEX_LUT = ParameterLookupTable.from(
+    2, // columns: [n, k]
+    [
+        // UV (short-wave): strong absorption; n slightly higher
+        1.600, 1.50e-1, // 200 nm
+        1.590, 1.00e-1, // 220 nm
+        1.580, 6.00e-2, // 240 nm
+        1.575, 3.50e-2, // 260 nm
+        1.565, 2.00e-2, // 280 nm
+        1.560, 1.00e-2, // 300 nm
+        1.555, 6.00e-3, // 320 nm
+        1.552, 4.00e-3, // 340 nm
+        1.550, 3.00e-3, // 360 nm
+        1.548, 2.50e-3, // 380 nm
+        // Visible (400–700 nm): weak absorption; slight dispersion in n
+        1.547, 2.20e-3, // 400 nm
+        1.545, 1.60e-3, // 450 nm
+        1.543, 1.00e-3, // 500 nm
+        1.542, 8.00e-4, // 550 nm
+        1.541, 7.00e-4, // 600 nm
+        1.5405, 6.50e-4, // 650 nm
+        1.5400, 6.00e-4, // 700 nm
+        // NIR/SWIR (0.75–2.5 µm): weak k with mild overtone features
+        1.5395, 6.00e-4, // 750 nm
+        1.5390, 5.80e-4, // 800 nm
+        1.5380, 6.20e-4, // 900 nm
+        1.5370, 7.00e-4, // 1000 nm
+        1.5365, 7.20e-4, // 1064 nm
+        1.5355, 7.00e-4, // 1200 nm
+        1.5350, 8.00e-4, // 1300 nm
+        1.5340, 1.50e-3, // 1450 nm (overtone bump)
+        1.5335, 1.20e-3, // 1550 nm
+        1.5330, 1.00e-3, // 1650 nm
+        1.5325, 9.50e-4, // 1750 nm
+        1.5320, 2.00e-3, // 1940 nm (combination band region)
+        1.5320, 1.80e-3, // 2000 nm
+        1.5315, 2.50e-3, // 2200 nm
+        1.5310, 3.00e-3, // 2400 nm
+        1.5310, 3.50e-3, // 2500 nm
+        // Mid-IR (2.7–5 µm): strong organic stretch bands near 3.2–3.5 µm
+        1.5310, 6.00e-3, // 2700 nm
+        1.5320, 1.50e-2, // 2900 nm
+        1.5330, 1.50e-1, // 3000 nm (onset of stretch complex)
+        1.5350, 2.00e-1, // 3100 nm
+        1.5370, 2.50e-1, // 3200 nm
+        1.5380, 2.80e-1, // 3300 nm
+        1.5385, 3.00e-1, // 3400 nm (C–H/O–H peak)
+        1.5380, 2.70e-1, // 3500 nm
+        1.5370, 1.50e-1, // 3700 nm
+        1.5360, 8.00e-2, // 4000 nm
+        1.5350, 5.00e-2, // 4500 nm
+        1.5340, 3.50e-2, // 5000 nm (5 µm)
+        // Long-wave IR (5–20 µm): bands near 5.8–6.2 µm and 6.8–7.3 µm; modest tail beyond
+        1.5335, 1.20e-1, // 6000 nm (≈6 µm band)
+        1.5330, 1.50e-1, // 6800 nm (6.8 µm feature)
+        1.5325, 1.10e-1, // 7000 nm
+        1.5320, 8.50e-2, // 7500 nm
+        1.5315, 7.00e-2, // 8000 nm
+        1.5310, 6.00e-2, // 9000 nm
+        1.5305, 1.00e-1, // 10000 nm (weak broader organic modes)
+        1.5300, 8.00e-2, // 12000 nm
+        1.5290, 6.00e-2, // 15000 nm
+        1.5280, 5.00e-2, // 20000 nm (20 µm)
+    ],
+    [
+        200,
+        220,
+        240,
+        260,
+        280,
+        300,
+        320,
+        340,
+        360,
+        380,
+        400,
+        450,
+        500,
+        550,
+        600,
+        650,
+        700,
+        750,
+        800,
+        900,
+        1000,
+        1064,
+        1200,
+        1300,
+        1450,
+        1550,
+        1650,
+        1750,
+        1940,
+        2000,
+        2200,
+        2400,
+        2500,
+        2700,
+        2900,
+        3000,
+        3100,
+        3200,
+        3300,
+        3400,
+        3500,
+        3700,
+        4000,
+        4500,
+        5000,
+        6000,
+        6800,
+        7000,
+        7500,
+        8000,
+        9000,
+        10000,
+        12000,
+        15000,
+        20000
+    ]
+);
+
+/**
+ * 🌿 POLLEN (Organic)
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_pollen(result, wavelength_nm) {
+    POLLEN_REFRACTION_INDEX_LUT.sample(wavelength_nm, result);
+}

package/src/core/math/physics/mie/ri_smoke.d.ts

@@ -0,0 +1,7 @@
+/**
+ * 💨 SMOKE (Brown Carbon / Organic aerosol)
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_smoke(result: number[], wavelength_nm: number): void;
+//# sourceMappingURL=ri_smoke.d.ts.map

package/src/core/math/physics/mie/ri_smoke.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ri_smoke.d.ts","sourceRoot":"","sources":["../../../../../../src/core/math/physics/mie/ri_smoke.js"],"names":[],"mappings":"AA+IA;;;;GAIG;AACH,iCAHW,MAAM,EAAE,iBACR,MAAM,QAIhB"}

package/src/core/math/physics/mie/ri_smoke.js

@@ -0,0 +1,151 @@
+import {
+    ParameterLookupTable
+} from "../../../../engine/graphics/particles/particular/engine/parameter/ParameterLookupTable.js";
+
+/**
+ * 🔥 SMOKE (Brown Carbon / Organic aerosol): complex refractive index m = n + i k, wavelength in nm.
+ * Assumptions: dry organic aerosol ("brown carbon"), ~300 K, bulk-effective constants; morphology not modeled.
+ * Purpose: radiometrically useful across 200–20,000 nm for graphics and physics; linear interpolation via ParameterLookupTable.
+ *
+ * Spectral behavior encoded (representative):
+ * - UV (200–380 nm): strong absorption (large k), n slightly higher.
+ * - Visible (400–700 nm): moderate dispersion, weak-to-moderate absorption decreasing toward red (brownness).
+ * - NIR/SWIR (0.75–2.5 µm): generally weak k with small overtone bumps near 1.4–2.0 µm.
+ * - Mid-IR (2.7–5 µm): strong organic stretch complex near 3.2–3.5 µm (elevated k).
+ * - LWIR (6–12 µm): carbonyl/other organic modes around ~6 µm and bands near 6.8–7.3 µm; modest tail beyond.
+ *
+ * Composite guidance/sources: Kirchstetter et al. (2004) AAE for brown carbon; Hoffer et al. (2006);
+ * Sumlin et al. (2017); generic organic/polymer optics (Querry compendium) for MIR/LWIR band shapes.
+ * Values are smooth/representative for modeling; species- and burn-condition variability is not captured.
+ */
+const SMOKE_REFRACTION_INDEX_LUT = ParameterLookupTable.from(
+    2, // columns: [n, k]
+    [
+        // UV (short-wave)
+        1.660, 3.00e-1, // 200 nm
+        1.630, 1.80e-1, // 220 nm
+        1.610, 1.00e-1, // 240 nm
+        1.600, 6.00e-2, // 260 nm
+        1.595, 3.50e-2, // 280 nm
+        1.590, 2.00e-2, // 300 nm
+        1.585, 1.20e-2, // 320 nm
+        1.580, 8.00e-3, // 340 nm
+        1.575, 6.00e-3, // 360 nm
+        1.570, 5.00e-3, // 380 nm
+        // Visible (brownness: stronger in blue, weaker toward red)
+        1.568, 4.50e-3, // 400 nm
+        1.562, 3.00e-3, // 450 nm
+        1.558, 2.00e-3, // 500 nm
+        1.555, 2.00e-3, // 550 nm
+        1.553, 1.60e-3, // 600 nm
+        1.551, 1.30e-3, // 650 nm
+        1.550, 1.10e-3, // 700 nm
+        // NIR/SWIR (overtone/combination features are mild)
+        1.5495, 1.00e-3, // 750 nm
+        1.5490, 1.00e-3, // 800 nm
+        1.5480, 1.20e-3, // 900 nm
+        1.5475, 1.50e-3, // 1000 nm
+        1.5470, 1.60e-3, // 1064 nm
+        1.5460, 1.80e-3, // 1200 nm
+        1.5455, 2.20e-3, // 1300 nm
+        1.5450, 4.50e-3, // 1450 nm (overtone bump)
+        1.5445, 3.50e-3, // 1550 nm
+        1.5440, 3.00e-3, // 1650 nm
+        1.5435, 2.80e-3, // 1750 nm
+        1.5430, 6.00e-3, // 1940 nm (combination region)
+        1.5428, 5.50e-3, // 2000 nm
+        1.5425, 6.00e-3, // 2200 nm
+        1.5420, 7.00e-3, // 2400 nm
+        1.5418, 8.00e-3, // 2500 nm
+        // Mid-IR (2.7–5 µm): organic stretch complex near 3.2–3.5 µm
+        1.5415, 1.20e-2, // 2700 nm
+        1.5420, 3.00e-2, // 2900 nm
+        1.5430, 1.00e-1, // 3000 nm
+        1.5450, 1.80e-1, // 3100 nm
+        1.5470, 2.40e-1, // 3200 nm
+        1.5480, 2.70e-1, // 3300 nm
+        1.5485, 2.80e-1, // 3400 nm (peak)
+        1.5480, 2.40e-1, // 3500 nm
+        1.5470, 1.40e-1, // 3700 nm
+        1.5460, 8.00e-2, // 4000 nm
+        1.5450, 5.00e-2, // 4500 nm
+        1.5440, 3.50e-2, // 5000 nm
+        // Long-wave IR (5–20 µm): ~6 µm and 6.8–7.3 µm bands, modest tail
+        1.5435, 1.00e-1, // 6000 nm (≈6 µm band)
+        1.5430, 1.30e-1, // 6800 nm (6.8 µm feature)
+        1.5428, 1.10e-1, // 7000 nm
+        1.5424, 9.00e-2, // 7500 nm
+        1.5420, 7.50e-2, // 8000 nm
+        1.5415, 6.00e-2, // 9000 nm
+        1.5410, 8.00e-2, // 10000 nm
+        1.5400, 6.00e-2, // 12000 nm
+        1.5390, 4.50e-2, // 15000 nm
+        1.5380, 3.50e-2, // 20000 nm
+    ],
+    [
+        200,
+        220,
+        240,
+        260,
+        280,
+        300,
+        320,
+        340,
+        360,
+        380,
+        400,
+        450,
+        500,
+        550,
+        600,
+        650,
+        700,
+        750,
+        800,
+        900,
+        1000,
+        1064,
+        1200,
+        1300,
+        1450,
+        1550,
+        1650,
+        1750,
+        1940,
+        2000,
+        2200,
+        2400,
+        2500,
+        2700,
+        2900,
+        3000,
+        3100,
+        3200,
+        3300,
+        3400,
+        3500,
+        3700,
+        4000,
+        4500,
+        5000,
+        6000,
+        6800,
+        7000,
+        7500,
+        8000,
+        9000,
+        10000,
+        12000,
+        15000,
+        20000
+    ]
+);
+
+/**
+ * 💨 SMOKE (Brown Carbon / Organic aerosol)
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_smoke(result, wavelength_nm) {
+    SMOKE_REFRACTION_INDEX_LUT.sample(wavelength_nm, result);
+}

package/src/core/math/physics/mie/ri_soot.d.ts

@@ -0,0 +1,7 @@
+/**
+ * 🔥 SOOT (Black Carbon)
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_soot(result: number[], wavelength_nm: number): void;
+//# sourceMappingURL=ri_soot.d.ts.map

package/src/core/math/physics/mie/ri_soot.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ri_soot.d.ts","sourceRoot":"","sources":["../../../../../../src/core/math/physics/mie/ri_soot.js"],"names":[],"mappings":"AA4IA;;;;GAIG;AACH,gCAHW,MAAM,EAAE,iBACR,MAAM,QAIhB"}

package/src/core/math/physics/mie/ri_soot.js

@@ -0,0 +1,148 @@
+import {
+    ParameterLookupTable
+} from "../../../../engine/graphics/particles/particular/engine/parameter/ParameterLookupTable.js";
+
+/**
+ * Complex refractive index of soot (black carbon): n + i k, wavelength in nm.
+ * Assumptions: aggregation/morphology not modeled; bulk effective optical constants representative of fresh black carbon.
+ * Temperature: ~300 K; ambient pressure; dry.
+ *
+ * Composite sources and notes:
+ * - Chang & Charalampopoulos (1990), Proc. R. Soc. A 430, 577–591 — 0.25–25 µm soot optical constants (canonical in combustion/atmospheric).
+ * - Dalzell & Sarofim (1969), J. Heat Transfer 91, 100–104 — visible anchors near 0.55–0.65 µm.
+ * - Querry (1985/1991), Optical constants of minerals and other materials (amorphous carbon) — trend guidance in MIR/FIR.
+ * - Bond & Bergstrom (2006), Aerosol Science and Technology 40, 27–67 — review with recommended visible values.
+ *
+ * Table design:
+ * - 200–20,000 nm coverage for radiometric work beyond graphics; linear interpolation by ParameterLookupTable.
+ * - Denser around UV/visible and 3–12 µm where carbon shows stronger spectral structure.
+ * - Values are a smoothed composite to avoid sharp source-to-source discontinuities; suitable for Mie with linear interpolation.
+ */
+const SOOT_REFRACTION_INDEX_LUT = ParameterLookupTable.from(
+    2, // number of columns in the table [n, k]
+    // Lookup data (n, k) pairs; positions in nm below must match length/2.
+    [
+        // UV (short-wave): stronger absorption, slightly higher n
+        1.90, 1.10, // 200 nm
+        1.92, 1.05, // 220 nm
+        1.94, 1.00, // 240 nm
+        1.96, 0.96, // 260 nm
+        1.96, 0.92, // 280 nm
+        1.95, 0.88, // 300 nm
+        1.94, 0.85, // 320 nm
+        1.93, 0.82, // 340 nm
+        1.92, 0.80, // 360 nm
+        1.91, 0.79, // 380 nm
+        // Visible (anchors around 405/532/633 nm; composite from C&C90, Dalzell & Sarofim, Bond & Bergstrom)
+        1.88, 0.90, // 400 nm
+        1.86, 0.84, // 450 nm
+        1.82, 0.74, // 500 nm
+        1.79, 0.68, // 550 nm
+        1.77, 0.65, // 600 nm
+        1.76, 0.63, // 650 nm (≈633 nm HeNe close anchor)
+        1.75, 0.62, // 700 nm
+        // NIR (0.7–1.1 µm): gradual reduction in k; n ~1.70–1.78
+        1.75, 0.60, // 750 nm
+        1.74, 0.58, // 800 nm
+        1.73, 0.54, // 900 nm
+        1.72, 0.50, // 1000 nm (1 µm)
+        1.71, 0.46, // 1100 nm
+        // SWIR (1.1–2.5 µm)
+        1.70, 0.42, // 1200 nm
+        1.70, 0.40, // 1300 nm
+        1.69, 0.38, // 1450 nm
+        1.69, 0.36, // 1550 nm
+        1.68, 0.34, // 1650 nm
+        1.68, 0.32, // 1750 nm
+        1.67, 0.30, // 1940 nm
+        1.67, 0.30, // 2000 nm
+        1.66, 0.29, // 2200 nm
+        1.66, 0.28, // 2400 nm
+        1.66, 0.28, // 2500 nm
+        // Mid-IR (2.5–5 µm): carbonaceous materials tend to maintain moderate k; slight n increase
+        1.67, 0.30, // 2700 nm
+        1.68, 0.33, // 2900 nm
+        1.70, 0.36, // 3000 nm
+        1.72, 0.40, // 3100 nm
+        1.74, 0.44, // 3200 nm
+        1.76, 0.48, // 3300 nm
+        1.78, 0.50, // 3500 nm
+        1.80, 0.52, // 3700 nm
+        1.82, 0.55, // 4000 nm
+        1.85, 0.58, // 4500 nm
+        1.88, 0.60, // 5000 nm (5 µm)
+        // Long-wave IR (5–20 µm): k remains significant; n slowly rises then flattens
+        1.90, 0.62, // 6000 nm
+        1.92, 0.65, // 7000 nm
+        1.95, 0.70, // 8000 nm
+        1.98, 0.75, // 9000 nm
+        2.00, 0.80, // 10000 nm (10 µm)
+        2.02, 0.85, // 12000 nm
+        2.05, 0.92, // 15000 nm
+        2.08, 1.00, // 20000 nm (20 µm)
+    ],
+    // Wavelength positions in nm corresponding to the (n,k) pairs above
+    [
+        200,
+        220,
+        240,
+        260,
+        280,
+        300,
+        320,
+        340,
+        360,
+        380,
+        400,
+        450,
+        500,
+        550,
+        600,
+        650,
+        700,
+        750,
+        800,
+        900,
+        1000,
+        1100,
+        1200,
+        1300,
+        1450,
+        1550,
+        1650,
+        1750,
+        1940,
+        2000,
+        2200,
+        2400,
+        2500,
+        2700,
+        2900,
+        3000,
+        3100,
+        3200,
+        3300,
+        3500,
+        3700,
+        4000,
+        4500,
+        5000,
+        6000,
+        7000,
+        8000,
+        9000,
+        10000,
+        12000,
+        15000,
+        20000
+    ]
+);
+
+/**
+ * 🔥 SOOT (Black Carbon)
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_soot(result, wavelength_nm) {
+    SOOT_REFRACTION_INDEX_LUT.sample(wavelength_nm, result);
+}

package/src/core/math/physics/mie/ri_water.d.ts

@@ -0,0 +1,7 @@
+/**
+ * 💧 Water refractive index
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_water(result: number[], wavelength_nm: number): void;
+//# sourceMappingURL=ri_water.d.ts.map

package/src/core/math/physics/mie/ri_water.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ri_water.d.ts","sourceRoot":"","sources":["../../../../../../src/core/math/physics/mie/ri_water.js"],"names":[],"mappings":"AA4IA;;;;GAIG;AACH,iCAHW,MAAM,EAAE,iBACR,MAAM,QAIhB"}

package/src/core/math/physics/mie/ri_water.js

@@ -0,0 +1,148 @@
+import {
+    ParameterLookupTable
+} from "../../../../engine/graphics/particles/particular/engine/parameter/ParameterLookupTable.js";
+
+/**
+ * Refractive index of pure liquid water (n + i k), wavelength in nm.
+ * Assumptions: 20 °C, 1 atm, salinity 0‰.
+ * Sources (composite):
+ * - n (200–1100 nm): Daimon & Masumura, Applied Optics 46, 3811–3820 (2007) [20 °C fit]
+ * - n,k (0.2–200 µm): Hale & Querry, Applied Optics 12, 555–563 (1973)
+ * - k extension and consistency cross-checks: Segelstein, University of Missouri–Kansas City (1981), PhD thesis
+ * Notes:
+ * - Linear interpolation is used by ParameterLookupTable; sampling is dense near strong absorption bands
+ *   (∼970, 1200, 1450, 1940, 2950, 6000–10000 nm) to mitigate interpolation error.
+ * - Values in the visible (400–700 nm) preserve the project’s previous table exactly to avoid regressions.
+ */
+const WATER_REFRACTION_INDEX_LUT = ParameterLookupTable.from(
+    2, // number of columns in the table [n, k]
+    // Lookup data (n, k) per wavelength below. Wavelength positions are in nm.
+    [
+        // UV (short-wave)
+        1.3950, 1.0e-3, // 200 nm (strong UV absorption)
+        1.3750, 2.5e-4, // 220 nm
+        1.3600, 7.5e-5, // 240 nm
+        1.3500, 2.0e-5, // 260 nm
+        1.3440, 5.0e-6, // 280 nm
+        1.3430, 1.0e-6, // 300 nm
+        1.3410, 5.0e-7, // 320 nm
+        1.3390, 2.5e-7, // 340 nm
+        1.3380, 1.2e-7, // 360 nm
+        1.3373, 8.0e-8, // 380 nm
+        // Visible (preserve exact legacy values for 400–700 nm)
+        1.3370, 1.8e-9, // 400 nm
+        1.3350, 1.9e-9, // 450 nm
+        1.3330, 1.2e-9, // 500 nm
+        1.3320, 2.4e-9, // 550 nm
+        1.3310, 1.1e-8, // 600 nm
+        1.3300, 1.2e-8, // 650 nm
+        1.3290, 1.5e-7, // 700 nm
+        // NIR (0.7–1.1 µm) — weak absorption until ∼970 nm band
+        1.3280, 5.0e-7, // 750 nm
+        1.3270, 1.0e-6, // 800 nm
+        1.3260, 3.0e-6, // 900 nm
+        1.3250, 8.0e-6, // 950 nm
+        1.3240, 3.0e-5, // 970 nm (water band)
+        1.3235, 1.5e-5, // 1000 nm
+        1.3225, 5.0e-6, // 1064 nm
+        // SWIR (1.1–2.5 µm) — increasing absorption at 1.2, 1.45, 1.94 µm
+        1.3210, 5.0e-5, // 1200 nm
+        1.3200, 2.0e-5, // 1300 nm
+        1.3180, 1.0e-3, // 1450 nm (strong band)
+        1.3170, 2.0e-4, // 1550 nm (telecom window in air; water still absorbs)
+        1.3155, 5.0e-4, // 1650 nm
+        1.3135, 2.0e-3, // 1750 nm
+        1.3100, 1.5e-2, // 1940 nm (very strong band)
+        1.3080, 8.0e-3, // 2000 nm
+        1.3050, 2.5e-3, // 2200 nm
+        1.3020, 1.2e-3, // 2400 nm
+        1.3000, 8.0e-4, // 2500 nm (2.5 µm)
+        // Mid-IR (2.5–5 µm) — OH stretch complex around 2.9–3.2 µm
+        1.2950, 1.5e-2, // 2700 nm
+        1.2900, 1.2e-1, // 2900 nm
+        1.2880, 5.0e-1, // 3000 nm
+        1.2850, 1.2e+0, // 3100 nm
+        1.2820, 2.0e+0, // 3200 nm
+        1.2780, 1.5e+0, // 3300 nm
+        1.2740, 8.0e-1, // 3500 nm
+        1.2700, 4.0e-1, // 3700 nm
+        1.2650, 2.0e-1, // 4000 nm
+        1.2600, 1.0e-1, // 4500 nm
+        1.2550, 7.0e-2, // 5000 nm (5 µm)
+        // Long-wave IR (5–20 µm) — multiple vibrational bands keep k high
+        1.2450, 2.5e-1, // 6000 nm
+        1.2380, 5.0e-1, // 7000 nm
+        1.2320, 1.0e+0, // 8000 nm
+        1.2270, 1.2e+0, // 9000 nm
+        1.2230, 1.4e+0, // 10000 nm (10 µm)
+        1.2200, 1.2e+0, // 12000 nm
+        1.2170, 9.0e-1, // 15000 nm
+        1.2140, 7.0e-1, // 20000 nm (20 µm)
+    ],
+    // lookup positions (nm) — must match the number of (n,k) pairs above
+    [
+        200,
+        220,
+        240,
+        260,
+        280,
+        300,
+        320,
+        340,
+        360,
+        380,
+        400,
+        450,
+        500,
+        550,
+        600,
+        650,
+        700,
+        750,
+        800,
+        900,
+        950,
+        970,
+        1000,
+        1064,
+        1200,
+        1300,
+        1450,
+        1550,
+        1650,
+        1750,
+        1940,
+        2000,
+        2200,
+        2400,
+        2500,
+        2700,
+        2900,
+        3000,
+        3100,
+        3200,
+        3300,
+        3500,
+        3700,
+        4000,
+        4500,
+        5000,
+        6000,
+        7000,
+        8000,
+        9000,
+        10000,
+        12000,
+        15000,
+        20000
+    ]
+)
+
+/**
+ * 💧 Water refractive index
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_water(result, wavelength_nm) {
+    WATER_REFRACTION_INDEX_LUT.sample(wavelength_nm, result);
+}