PyStratify against scattnlay

← Verification · Spheres · Plane wave
Reference
scattnlay
Kind
code: another group's code
Calculated
Q_ext, Q_sca, Q_abs, Q_back, S₁, S₂ at 0–180° for 10 spheres up to x = 1000 (double and multiple precision)
Agreement
1e-10 with 20 orders beyond Wiscombe (observed ≤ 1.8e-11); 1e-7 at PyStratify's default truncation (observed ≤ 3.7e-8)
A gold-like matryoshka (ε = −10 + i / silica / ε = −10 + i, radii in proportion 1 : 2 : 3) in water, 600 nm: against scattnlay scattnlay (Peña & Pal; Ladutenko et al.) evaluates the multilayer Mie coefficients by its own recursion; both at the default truncation. The test covers ten spheres to x = 1000, in double and multiple precision.
2040608010012001234567outer radius [nm]Q
Qext, PyStratify Qsca, PyStratify Qext, scattnlay Qsca, scattnlay
Hover to read values.
  • largest relative difference: 4.68×10⁻¹⁰
  • PyStratify the lines, the reference the circles
A lossless water-like sphere (m = 1.33) at x = 200: the scattering amplitudes against scattnlay Six decades from the forward peak to the backward glory, with the rainbow near 138°; PyStratify 20 orders beyond Wiscombe's truncation, scattnlay at its own.
02040608010012014016018010⁻¹11010²10³10⁴10⁵scattering angle [°]|S|
|S1|, PyStratify |S2|, PyStratify |S1|, scattnlay |S2|, scattnlay
Hover to read values.
  • largest difference relative to the forward peak: 6.52×10⁻¹⁴
  • PyStratify the lines, the reference the circles

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