Sommerfeld integrals over s = k_ρ/k₁ with Parratt-recursion r_s, r_p, written without pystratify

← Verification · Thin films · Point dipole
Reference
Sommerfeld integrals over s = k_ρ/k₁ with Parratt-recursion r_s, r_p, written without pystratify
Kind
method: an independent formulation inside this repo (a scattering matrix, a plane-wave superposition, a quadrature): proves the formulation
Calculated
perpendicular and parallel total rates, ED and MD, 10–300 nm above glass, silver, absorbing films; a TIR case
Agreement
1e-10; observed ≤ 5e-12
A dipole in air above 60 nm of silver on glass, 633 nm The reference integrates the Chance–Prock–Silbey integrands with Parratt's recursion for r_s and r_p, written without PyStratify; quenching rises as d⁻³ at contact and the surface plasmon of the silver sets the breakpoint of the evanescent integral.
5010015020025030035040010⁻¹11010²height above the silver [nm]Γ / Γ₀
perpendicular, PyStratify parallel, PyStratify perpendicular, Sommerfeld integral parallel, Sommerfeld integral
Hover to read values.
  • largest relative difference: 3.37×10⁻¹²
  • PyStratify the lines, the reference the circles

Run it yourself: emitter_references.py · test_planar_emission_references.py · figures.py on GitHub, PyStratify 97e52a9