Cylinder: Goniometer

This example demonstrates how to use a goniometer setup to measure and visualize the coupling efficiency as a function of angular displacement for cylindrical scatterers using PyMieSim.

Importing the package dependencies: numpy, PyMieSim

import numpy as np
from TypedUnit import ureg

from PyMieSim.experiment.detector import Photodiode
from PyMieSim.experiment.scatterer import Cylinder
from PyMieSim.experiment.source import Gaussian
from PyMieSim.experiment import Setup
from PyOptik import Material

Defining the source

source = Gaussian(
    wavelength=1200 * ureg.nanometer,  # 1200 nm
    polarization=90 * ureg.degree,  # Polarization angle in ureg.degrees
    optical_power=1e-3 * ureg.watt,  # 1 milliureg.watt
    NA=0.2 * ureg.AU,  # Numerical Aperture
)

Defining the scatterer distribution

scatterer = Cylinder(
    diameter=2000 * ureg.nanometer,  # 2000 nm
    property=Material.BK7,  # Material of the cylinder
    medium_property=1 * ureg.RIU,  # Refractive index of the surrounding medium
    source=source,
)

Defining the detector

detector = Photodiode(
    NA=[0.5, 0.3, 0.1, 0.05] * ureg.AU,  # Array of Numerical Apertures for the detector
    phi_offset=np.linspace(-180, 180, 200)
    * ureg.degree,  # Angular displacement from -180 to 180 ureg.degrees
    gamma_offset=0 * ureg.degree,  # Gamma offset in ureg.degrees
    sampling=400 * ureg.AU,  # Number of sampling points
    polarization_filter=None,  # No polarization filter
)

Setting up the experiment

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

Measuring the coupling efficiency

dataframe = experiment.get("coupling", scale_unit=True)

Plotting the results Visualizing how the coupling efficiency varies with angular displacement.

dataframe.plot(x="detector:phi_offset")

<Figure size 800x500 with 1 Axes>