Note
Go to the end to download the full example code.
Sphere: Qabs vs diameter#
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
source = Gaussian(
wavelength=400 * ureg.nanometer,
polarization=0 * ureg.degree,
optical_power=1e-3 * ureg.watt,
NA=0.2 * ureg.AU,
)
scatterer = Sphere(
diameter=np.linspace(1, 800, 300) * ureg.nanometer,
property=[Material.silver],
medium_property=1 * ureg.RIU,
source=source,
)
experiment = Setup(scatterer=scatterer, source=source)
dataframe = experiment.get("Qsca")
dataframe.plot(x="scatterer:diameter", yscale="log")

<Figure size 800x500 with 1 Axes>
Total running time of the script: (0 minutes 0.353 seconds)