Your first simulation#

Use PyMieSim.Simulation for one source, one scatterer, and an optional detector. This is the simplest way to calculate scattering properties.

from PyMieSim import (
    Gaussian,
    PolarizationState,
    Simulation,
    Sphere,
    ureg,
)

source = Gaussian(
    wavelength=750 * ureg.nanometer,
    polarization=PolarizationState(angle=0 * ureg.degree),
    optical_power=1e-3 * ureg.watt,
    numerical_aperture=0.2,
)

simulation = Simulation(
    scatterer=Sphere(
        diameter=200 * ureg.nanometer,
        material=1.5,
        medium=1.0,
    ),
    source=source,
)

result = simulation.run("Qsca", "Qext")
print(result)

Next: Parameter sweeps.