Sphere: Qsca vs wavelength STD#

Importing the package dependencies: numpy, PyMieSim

import numpy as np
from TypedUnit import ureg

from PyMieSim.experiment.scatterer import Sphere
from PyMieSim.experiment.source import Gaussian
from PyMieSim.experiment import Setup
from PyOptik import Material

Defining the source to be employed.

source = Gaussian(
    wavelength=np.linspace(200, 1800, 300) * ureg.nanometer,
    polarization=0 * ureg.degree,
    optical_power=1e-3 * ureg.watt,
    NA=0.2 * ureg.AU,
)

Defining the ranging parameters for the scatterer distribution

scatterer = Sphere(
    diameter=np.linspace(400, 1400, 10) * ureg.nanometer,
    property=Material.silver,
    medium_property=1 * ureg.RIU,
    source=source,
)

Defining the experiment setup

experiment = Setup(scatterer=scatterer, source=source)

Measuring the properties

dataframe = experiment.get("Qsca")

Plotting the results

dataframe.plot(x="source:wavelength", std="scatterer:diameter")
sphere Qsca vs wavelength std
<Figure size 800x500 with 1 Axes>

Total running time of the script: (0 minutes 0.381 seconds)

Gallery generated by Sphinx-Gallery