PyStratify against tmm

← Verification · Thin films · 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.

tmm computes with Byrnes's coherent and incoherent transfer matrices. Each case computes the same quantity with PyStratify and with tmm and plots both on the same axes.

Incoherent stack

PyStratify's incoherent solver against inctmm, the code it replaced

Open in calculator
5006007008009000.10.20.30.40.50.6λ [nm]R, T, A
R, PyStratify (pimtmm) T, PyStratify (pimtmm) A, PyStratify (pimtmm) R, inctmm T, inctmm A, inctmm
RTAPyStratify (pimtmm)inctmm, every 15 nm

Source

Source typePlane wave
θ from normal [°]31
λ min [nm]450
λ max [nm]900
λ step [nm]5
kpsd 12200 nmkpsd 12200 nm

Thin films · Plane wave

LabelMaterialnkThickness [nm]Coherence
↓Inconst10∞–
1Film 1const2.080.0390auto
2Thick layerconst1.470.00112000incoherent
3Film 2const1.820.02110auto
↓Outconst1.560∞–

Backside reflection

the glass's back reflection at 65° against inctmm and a phase average

Open in calculator
400500600700800900-0.100.10.20.30.40.50.6λ [nm]R, Re r_s r̄_p
Σ|rs|², PyStratify Σ|rp|², PyStratify Re Σ rs r̄p, PyStratify Rs, inctmm Rp, inctmm Re rs r̄p, phase average
Σ|rs|²Σ|rp|²Re Σ rs r̄pPyStratifyRs, inctmm, every 20 nmRp, inctmm, every 20 nmRe rs r̄p, phase average, every 20 nm

Source

Source typePlane wave
θ from normal [°]65
λ min [nm]400
λ max [nm]900
λ step [nm]10
kpsd 1.0012e+5 nmkpsd 1.0012e+5 nm

Thin films · Plane wave

LabelMaterialnkThickness [nm]Coherence
↓Inconst10∞–
1FilmoscillatorCauchy 2.0 + 0.01/λ²–120auto
2Glassconst1.520100000incoherent
↓Outconst10∞–