PyStratify against smuthi

← Verification · Spheres · Point dipole
Reference
smuthi (dipole source and a sphere, by Mie T-matrix and its own translation of the dipole field)
Kind
code: another group's code
Calculated
total rate, radial and tangential ED, of four spheres: dielectric, a metal-like 40 nm sphere 5 nm away, in water, a 1 µm n = 3.5 sphere
Agreement
converges geometrically in smuthi's l_max: 2.2e-16 (dielectric, water), 5.1e-10 (metal, l_max 60), 1.1e-6 (large, l_max 60); asserted 1e-12, 1e-12, 1e-5 (l_max 40)
A dipole 100–1000 nm from a dielectric sphere, n = 2, R = 150 nm, in air, 800 nm smuthi's T-matrix truncated at l_max = 21 nearest the sphere, 10 farthest; PyStratify at tolerance 1e-10.
20040060080010000.911.11.21.31.41.5gap [nm]Γ / Γ₀
radial, PyStratify tangential, PyStratify radial, smuthi tangential, smuthi
Hover to read values.
  • largest relative difference: 8.88×10⁻¹⁶
  • PyStratify the lines, the reference the circles
A dipole 5–60 nm from a metal-like sphere, ε = −10 + i, R = 20 nm, in air, 600 nm smuthi's T-matrix truncated at l_max = 40 nearest the sphere, 10 farthest; PyStratify at tolerance 1e-10.
10203040506010⁻¹11010²gap [nm]Γ / Γ₀
radial, PyStratify tangential, PyStratify radial, smuthi tangential, smuthi
Hover to read values.
  • largest relative difference: 1.82×10⁻⁶
  • PyStratify the lines, the reference the circles
A dipole 20–300 nm from a sphere of n = 1.5, R = 100 nm, in water (n = 1.33), 633 nm smuthi's T-matrix truncated at l_max = 40 nearest the sphere, 10 farthest; PyStratify at tolerance 1e-10.
501001502002503000.911.11.21.3gap [nm]Γ / Γ₀
radial, PyStratify tangential, PyStratify radial, smuthi tangential, smuthi
Hover to read values.
  • largest relative difference: 4.44×10⁻¹⁶
  • PyStratify the lines, the reference the circles
A dipole 50–1000 nm from a large high-index sphere, n = 3.5 + 0.01i, R = 500 nm, in air, 800 nm smuthi's T-matrix truncated at l_max = 40 nearest the sphere, 13 farthest; PyStratify at tolerance 1e-10.
20040060080010001234gap [nm]Γ / Γ₀
radial, PyStratify tangential, PyStratify radial, smuthi tangential, smuthi
Hover to read values.
  • largest relative difference: 2.55×10⁻⁵
  • PyStratify the lines, the reference the circles

Run it yourself: smuthi_spheres.py · test_smuthi_spheres.py · figures.py · SMUTHI_SPHERES.md on GitHub, PyStratify 97e52a9