@galihru/mnp-mie 0.1.0

package/README.md

@@ -0,0 +1,31 @@
+# @galihru/mnp-mie
+
+Rayleigh-limit scattering formulas for small spherical particles.
+
+Implemented formulations:
+
+$$
+\alpha = 4\pi a^3\frac{\varepsilon_p-\varepsilon_m}{\varepsilon_p+2\varepsilon_m}
+$$
+
+$$
+C_{\mathrm{ext}} = k\,\mathrm{Im}(\alpha),\quad
+C_{\mathrm{sca}} = \frac{|k|^4}{6\pi}|\alpha|^2,\quad
+C_{\mathrm{abs}} = C_{\mathrm{ext}}-C_{\mathrm{sca}}
+$$
+
+## Install
+
+```bash
+npm install @galihru/mnp-mie
+```
+
+## Usage
+
+```js
+import { complex, rayleighCrossSections } from "@galihru/mnp-mie";
+
+const epsParticle = complex(-5.2, 2.1);
+const epsMedium = complex(1.77, 0.0);
+console.log(rayleighCrossSections(548.1, 50, epsParticle, epsMedium));
+```

package/package.json

@@ -0,0 +1,27 @@
+{
+  "name": "@galihru/mnp-mie",
+  "version": "0.1.0",
+  "description": "Rayleigh-limit scattering formulations for spherical particles.",
+  "type": "module",
+  "main": "./src/index.js",
+  "exports": {
+    ".": "./src/index.js"
+  },
+  "files": [
+    "src",
+    "README.md"
+  ],
+  "author": "GALIH RIDHO UTOMO <g4lihru@students.unnes.ac.id>",
+  "license": "GPL-2.0-only",
+  "repository": {
+    "type": "git",
+    "url": "https://github.com/galihru/mnpbem.git"
+  },
+  "publishConfig": {
+    "registry": "https://registry.npmjs.org/",
+    "access": "public"
+  },
+  "scripts": {
+    "test": "node --test"
+  }
+}

package/src/complex.js

@@ -0,0 +1,35 @@
+export function complex(re = 0, im = 0) {
+  return { re: Number(re), im: Number(im) };
+}
+
+export function add(a, b) {
+  return complex(a.re + b.re, a.im + b.im);
+}
+
+export function sub(a, b) {
+  return complex(a.re - b.re, a.im - b.im);
+}
+
+export function mul(a, b) {
+  return complex(a.re * b.re - a.im * b.im, a.re * b.im + a.im * b.re);
+}
+
+export function div(a, b) {
+  const denom = b.re * b.re + b.im * b.im;
+  if (denom === 0) {
+    throw new Error("Division by zero complex number");
+  }
+  return complex((a.re * b.re + a.im * b.im) / denom, (a.im * b.re - a.re * b.im) / denom);
+}
+
+export function abs(z) {
+  return Math.hypot(z.re, z.im);
+}
+
+export function imag(z) {
+  return z.im;
+}
+
+export function fromReal(x) {
+  return complex(x, 0);
+}

package/src/index.js

@@ -0,0 +1,2 @@
+export { complex } from "./rayleigh.js";
+export { rayleighCrossSections, rayleighPolarizability } from "./rayleigh.js";

package/src/rayleigh.js

@@ -0,0 +1,34 @@
+import { abs, add, complex, div, fromReal, imag, mul, sub } from "./complex.js";
+
+function mapInput(values, fn) {
+  if (Array.isArray(values)) {
+    return values.map(fn);
+  }
+  return fn(values);
+}
+
+export function rayleighPolarizability(radiusNm, epsParticle, epsMedium) {
+  const a = Number(radiusNm);
+  const top = sub(epsParticle, epsMedium);
+  const bottom = add(epsParticle, mul(fromReal(2), epsMedium));
+  return mul(fromReal(4 * Math.PI * a ** 3), div(top, bottom));
+}
+
+export function rayleighCrossSections(wavelengthNm, radiusNm, epsParticle, epsMedium) {
+  const solveOne = (wl, ep, em) => {
+    const kFactor = (2 * Math.PI * Math.sqrt(em.re)) / Number(wl);
+    const k = fromReal(kFactor);
+    const alpha = rayleighPolarizability(radiusNm, ep, em);
+    const cExt = k.re * imag(alpha);
+    const cSca = (Math.abs(k.re) ** 4 * abs(alpha) ** 2) / (6 * Math.PI);
+    const cAbs = cExt - cSca;
+    return { cExt, cSca, cAbs };
+  };
+
+  if (Array.isArray(wavelengthNm)) {
+    return wavelengthNm.map((wl, i) => solveOne(wl, epsParticle[i], epsMedium[i]));
+  }
+  return solveOne(wavelengthNm, epsParticle, epsMedium);
+}
+
+export { complex };