opcsim.mie.cscat¶

opcsim.mie.cscat(dp, wl, refr, theta1, theta2, nsteps=100, \*\*kwargs)¶

Compute the scattering cross section between two angles according to Jaenicke and Hanusch (1993).

Following the lead of Jaenicke and Hanusch (1993), we can compute the scattering cross section for a given viewing angle [\(\Theta_1\) - \(\Theta_2\)] as:

\[C_{sca}=\frac{\lambda^2}{4\pi} \int_{\Theta_1}^{\Theta_2}[i_1(\Theta) + i_2(\Theta)]sin\Theta d\Theta\]

where \(\lambda\) is the incident wavelength of light and \(i_1\) and \(i_2\) are the intensity distribution functions, calculated as:

\[i_1(\Theta)=\mid S_1(\Theta) \mid^2,\]

\[i_2(\Theta)=\mid S_2(\Theta) \mid^2\]

The integral is calculated step-wise using the numpy.trapz function.

Parameters

dp: float: The particle diameter in microns.
wl: float: The wavelength of incident light in microns.
refr: complex: The complex refractive index of the material.
theta1: float: The angle from which to begin the integration.
theta1: float: The angle from which to end the integration.
nsteps: int: The number of steps in theta to use in performing the step-wise integration.

Returns

C_{scat}: float: The scattering cross-section.