Consistency checks

← 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.

Checks that need no other code: closed forms, limits and conservation laws, the same series summed in extended precision, and independent routes through PyStratify's own physics. Each case plots PyStratify against its reference on the same axes.

Dielectric cylinder

Qsca against the Bessel boundary equations solved directly, radius swept

Open in calculator
100200300400500600012345radius [nm]Qsca
E along the axis, PyStratify H along the axis, PyStratify E along the axis, Bessel equations H along the axis, Bessel equations
E along the axisH along the axisPyStratifyBessel equations, every 15 nm

Source

Source typePlane wave
θ from axis [°]90
λ [nm]600
R min [nm]20
R max [nm]600
R step [nm]5
Cross sectionkp: electric field in the wavevector–cylinder-axis plane; transverse and axial components are drawn separatelyps: electric field perpendicular to the wavevector–cylinder-axis planesR 300 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 300 nm
Cross sectionkp: electric field in the wavevector–cylinder-axis plane; transverse and axial components are drawn separatelyps: electric field perpendicular to the wavevector–cylinder-axis planesR 300 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 300 nm

Cylinder · Plane wave

LabelMaterialnkThickness [nm]
CCylinderconst1.5020–600
HAirconst10∞

Index-matched cylinder

a cylinder of its host's index is no cylinder: the field is the plane wave

Open in calculator
-200-1000100200-1-0.500.51x [nm]Re E
E along the axis, Ez, PyStratify H along the axis, Ey, PyStratify E along the axis, Ez, the plane wave H along the axis, Ey, the plane wave
E along the axis, EzH along the axis, EyPyStratifythe plane wave, every 20 nm

Source

Source typePlane wave
θ from axis [°]82.69
λ [nm]600
x min [nm]-250
x max [nm]250
x step [nm]2.5
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.50100
HHostconst1.50∞