PyStratify against treams

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

treams computes with T-matrices. Each case computes the same quantity with PyStratify and with treams and plots both on the same axes.

Oblique cylinder

efficiencies Q (cross width over the diameter) at 60° against treams

Open in calculator
4005006007008009000.20.40.60.811.21.41.6λ [nm]Q
Qext, PyStratify Qsca, PyStratify Qext, treams Qsca, treams
QextQscaPyStratifytreams, every 20 nm

Source

Source typePlane wave
θ from axis [°]60
λ min [nm]400
λ max [nm]900
λ step [nm]10
Cross sectionk⊥Wavevector component along the cylinder axis, out of the pagekzp: electric field in the wavevector–cylinder-axis plane; transverse and axial components are drawn separatelyp⊥pzs: electric field perpendicular to the wavevector–cylinder-axis planesR 100 nm
Sidekp: electric field in the wavevector–cylinder-axis plane; transverse and axial components are drawn separatelyps: electric field perpendicular to the wavevector–cylinder-axis planesR 100 nm
Cross sectionk⊥Wavevector component along the cylinder axis, out of the pagekzp: electric field in the wavevector–cylinder-axis plane; transverse and axial components are drawn separatelyp⊥pzs: electric field perpendicular to the wavevector–cylinder-axis planesR 100 nm
Sidekp: electric field in the wavevector–cylinder-axis plane; transverse and axial components are drawn separatelyps: electric field perpendicular to the wavevector–cylinder-axis planesR 100 nm

Cylinder · Plane wave

LabelMaterialnkThickness [nm]
CCylinderconst1.50.01100
HAirconst10∞