Sphere: Qsca vs diameter

Importing the package dependencies: numpy, PyMieSim

import numpy as np

from PyMieSim.experiment.scatterer import Sphere, CoreShell
from PyMieSim.experiment.source import Gaussian
from PyMieSim.experiment.detector import Photodiode

from PyMieSim.experiment import Setup
from PyMieSim import units
from PyMieSim.units import nanometer, degree, watt, AU, RIU
from PyOptik import Material
Material.print_available()

+-----------------------+-----------------------+
| Sellmeier Materials   | Tabulated Materials   |
+=======================+=======================+
| lithium_niobate       | calcium               |
+-----------------------+-----------------------+
| BK7                   | copper                |
+-----------------------+-----------------------+
| air                   | sodium                |
+-----------------------+-----------------------+
| silicon               | aluminium             |
+-----------------------+-----------------------+
| fused_silica          | iron                  |
+-----------------------+-----------------------+
| SF5                   | gold                  |
+-----------------------+-----------------------+
| ZBLAN                 | silver                |
+-----------------------+-----------------------+
| flint                 | nickel                |
+-----------------------+-----------------------+
| germanium             | zinc                  |
+-----------------------+-----------------------+
| soda_lime_glass       |                       |
+-----------------------+-----------------------+
| water                 |                       |
+-----------------------+-----------------------+
| argon                 |                       |
+-----------------------+-----------------------+
| polystyren            |                       |
+-----------------------+-----------------------+
| BAK1                  |                       |
+-----------------------+-----------------------+
| crown                 |                       |
+-----------------------+-----------------------+

Defining the source to be employed.

source_1 = Gaussian(
    wavelength=[405, 810] * nanometer,
    polarization=0 * degree,
    optical_power=1e-3 * watt,
    NA=0.2 * AU
)

scatterer_1 = Sphere(
    diameter=np.linspace(10, 1000, 150) * nanometer,
    medium_property=[1.33] * RIU,
    property=Material.polystyren,
    source=source_1
)

detector_1 = Photodiode(
    NA=[0.9 / 1.33, 1.3 / 1.33, 0.2, 0.1] * units.AU,
    gamma_offset=[0] * units.degree,
    phi_offset=90 * units.degree,

)

Defining the experiment setup

experiment_1 = Setup(scatterer=scatterer_1, source=source_1, detector=detector_1)

dataframe_1 = experiment_1.get('coupling', scale_unit=True, drop_unique_level=True)

# dataframe_1._plot(x='scatterer:diameter', std='detector:NA', show=True)
dataframe_1.plot(x='scatterer:diameter', show=True)









# plt.show()

/opt/hostedtoolcache/Python/3.11.11/x64/lib/python3.11/site-packages/PyOptik/material/base_class.py:60: UserWarning: Wavelength range goes from 405.0 nm to 810.0 nm which is outside the allowable range of 436.79999999999995 nm to 1.052 µm µm. [Material: polystyren]
  warnings.warn(

<Axes: xlabel='scatterer:diameter', ylabel='microwatt'>