PyStratify against MNPBEM

← Verification · Spheres · Plane wave

Under active development. Verification is growing as comparisons are added and rerun: the matrix, its figures and their presentation are not yet in their final form. Spotted an error, or missing a topic? Tell us.

MNPBEM computes with boundary elements, under GNU Octave. Each case computes the same quantity with PyStratify and with MNPBEM and plots both on the same axes.

Metal-like sphere

cross sections against MNPBEM's boundary elements (percent-level by design)

Open in calculator
400450500550600650700750800050100150λ [nm]C [nm²]
Cext, PyStratify Csca, PyStratify Cext, MNPBEM (2044 faces) Csca, MNPBEM (2044 faces)
CextCscaPyStratifyMNPBEM (2044 faces), every 25 nm

Source

Source typePlane wave
λ min [nm]400
λ max [nm]800
λ step [nm]2.5
R 20 nm k R 20 nm k

Sphere · Plane wave

LabelMaterialnkThickness [nm]
CMetal-likeconst0.1579173.1662220
HVacuumconst10∞

Thin nanoshell

MNPBEM's boundary elements, 11–16% low: its mesh cannot resolve a 5 nm shell

Open in calculator
400450500550600650700750800500100015002000λ [nm]Cext[nm²]
Cext, PyStratify Cext, MNPBEM (1016 faces)
CextPyStratifyMNPBEM (1016 faces), every 25 nm

Source

Source typePlane wave
λ min [nm]400
λ max [nm]800
λ step [nm]2.5
R 20 nm k R 20 nm k

Sphere · Plane wave

LabelMaterialnkThickness [nm]
CSilicaconst1.45945015
1Metal-likeconst0.1579173.166225
HVacuumconst10∞