Stokes Parameters Computation

This example demonstrates the computation and visualization of the Stokes parameters using PyMieSim.

from PyMieSim.units import ureg

from PyMieSim.single.scatterer import Sphere
from PyMieSim.single.source import Gaussian
from PyMieSim.polarization import PolarizationState
from PyMieSim.single.setup import Setup

polarization_state = PolarizationState(angle=0 * ureg.degree)

source = Gaussian(
    wavelength=750 * ureg.nanometer,
    polarization=polarization_state,
    optical_power=1 * ureg.watt,
    numerical_aperture=0.3,
)

scatterer = Sphere(
    diameter=600 * ureg.nanometer,
    medium=1.0,
    material=1.4,
)

setup = Setup(
    scatterer=scatterer,
    source=source,
)

stokes = setup.get_representation("stokes", sampling=100)

figure = stokes.plot()