PyStratify against treams: achiral spheres and multilayers (x ≤ 15)

← Verification · Spheres · Plane wave

Reference: treams, achiral spheres and multilayers (x ≤ 15)
Kind: code: another group's code
Calculated: Q_ext, Q_sca, Q_abs of 8 cases incl. nanoshells, a matryoshka and a 10-layer graded sphere
Agreement: 1e-10 with 20 orders beyond Wiscombe; observed ≤ 1.1e-14

A nanoshell (10 nm silica core, 10 nm shell of ε = −10 + i) treams (Beutel et al.) builds the T-matrix of a layered sphere from its own vector spherical waves; the fixed metal index keeps the shell's plasmon at one wavelength, so the comparison is of the two codes, not of a material.

Q_ext, PyStratify Q_sca, PyStratify Q_ext, treams Q_sca, treams

Hover to read values.

largest relative difference: 8.1×10⁻¹⁵
PyStratify the lines, the reference the circles

Run it yourself: mie_codes.py · test_mie_codes.py · figures.py on GitHub, PyStratify 97e52a9