PyStratify against pyGDM2

← Verification · Spheres · Point dipole

Reference: pyGDM2 decay_rate
Kind: code: another group's code
Calculated: total decay rates 20 and 40 nm from the n = 2 and n = 2 + 0.1i spheres, radial and tangential
Agreement: 0.5% at 40 nm (observed ≤ 0.14%), 10% at 20 nm (observed ≤ 7.2%)

A dipole beside a dielectric sphere (n = 2, R = 40 nm) in air, 600 nm: total decay rates against pyGDM2 pyGDM2's decay_rate couples the emitter to the dipole lattice through the Green dyadic; distances from the centroid of the lattice, PyStratify for the sphere of the same volume. The lattice is coarse close to the surface (the test: 10% at 20 nm, 0.5% at 40 nm).

radial, PyStratify tangential, PyStratify radial, pyGDM2 (1405 dipoles) tangential, pyGDM2 (1405 dipoles)

Hover to read values.

largest relative difference: 0.0261435
PyStratify the lines, the reference the circles

Run it yourself: gdm_spheres.py · test_gdm_spheres.py · figures.py on GitHub, PyStratify 97e52a9