Source-linked AI summary
Retrieving the Aerosol Complex Refractive Index using PyMieScatt: A Mie Computational Package with Visualization Capabilities
Benjamin J. Sumlin, William R. Heinson, Rajan K. Chakrabarty
TL;DR
The paper addresses the difficulty of inferring a particle’s complex refractive index from measured optical properties. It presents PyMieScatt, an open-source Python package with visual and numerical inverse-Mie methods, and shows how the methods identify valid solutions while exposing ambiguity and measurement uncertainty. The package also supports broad Mie calculations for particles and size distributions.
Problem
The complex refractive index cannot be directly measured and must be inferred from scattering and absorption data, while inverse Mie equations are difficult to solve without computers.
Method
PyMieScatt retrieves complex refractive index by intersecting optical-measurement contours in n-k space and also provides a numerical survey-iteration method under spherical-particle and known-size assumptions.
Results
Polydisperse aerosol inversions were unlikely to encounter multiple-solution problems, while visual error regions revealed where solutions may exist within laboratory measurement uncertainty.
Takeaways & Limitations
PyMieScatt provides an open-source, documented toolkit for inverse refractive-index retrieval and broader optical analysis of spherical particles and ensembles.
Takeaways & Limitations
The approach assumes spherical, homogeneous particles or ensembles with known or assumed size information and known illumination wavelength.
Abstract
from arXiv · showhide
The complex refractive index m=n+ik of a particle is an intrinsic property which cannot be directly measured, it must be inferred from its extrinsic properties such as the scattering and absorption cross-sections. Bohren and Huffman called this approach describing the dragon from its tracks, since the inversion of Lorenz-Mie theory equations is intractable without the use of computers. This article describes PyMieScatt, an open-source module for Python that contains functionality for solving the inverse problem for complex m using extensive optical and physical properties as input, and calculating regions where valid solutions may exist within the error bounds of laboratory measurements. Additionally, the module has comprehensive capabilities for studying homogeneous and coated single spheres, as well as ensembles of homogeneous spheres with user-defined size distributions, making it a complete tool for studying the optical behavior of spherical particles.
2 McDonnell Center for the Space Sciences
PyMieScatt addresses the inverse Mie problem by retrieving complex refractive indices from optical measurements under constrained particle assumptions. Its visual and numerical methods expose uncertainty, multiple mathematical solutions, and the physical validity of candidate solutions.
- Motivation: PyMieScatt fills a gap in Python by providing open-source Mie calculations and inverse retrieval of complex refractive index for known particle size parameters.The package is designed for standalone calculations and custom research or educational scripts.
- Inverse methods: The survey-iteration method offers a faster numerical alternative, but its strictly numerical treatment cannot accurately represent input-measurement error regions.It estimates distinct solutions through a coarse survey and then minimizes errors from multiple initial guesses.
- Constraining solutions: For size-distribution measurements, the retrieved effective index was m = 1.576+0.029i, while overlapping uncertainty regions identify indices consistent with measurement error.The visual method shows both where valid solutions may exist and where they may not.
Supplementary Material for
The passage identifies the paper as “Retrieving the Aerosol Complex Refractive Index using PyMieScatt.”
- The paper title is “Retrieving the Aerosol Complex Refractive Index using PyMieScatt.”
- The title names PyMieScatt as the tool associated with retrieving the aerosol complex refractive index.
- The title frames aerosol complex refractive-index retrieval as the paper’s subject.
A Mie Computational Package with Visualization Capabilities
The passage lists Benjamin J. Sumlin, William R. Heinson, and Rajan K. Chakrabarty as the paper’s authors.
- Benjamin J. Sumlin is listed as an author.
- William R. Heinson is listed as an author.
- Rajan K. Chakrabarty is listed as an author with affiliation markers 1,2 and a corresponding-author marker.
1 Center for Aerosol Science and Engineering
The listed affiliation is the Center for Aerosol Science and Engineering at Washington University in St. Louis, Missouri, USA.
- The authors are affiliated with the Center for Aerosol Science and Engineering.
- The affiliation is at Washington University in St. Louis, Missouri, USA, ZIP code 63130.
2 McDonnell Center for the Space Sciences
PyMieScatt provides Python tools for computing Mie optical properties across particle types and size distributions, including inverse retrieval of refractive index. Its algorithms support visualization, coated spheres, angular scattering, and benchmarked efficiency calculations.
- Library design: PyMieScatt accepts particle diameter and wavelength as user-facing inputs while using the size parameter x=πd/λ internally.This interface matches experimental workflows where diameter and wavelength may come from separate instruments.
- Homogeneous spheres: MieQ() computes efficiencies and the asymmetry parameter from the complex refractive index and size parameter using Mie coefficients.The library derives extinction, scattering, absorption, backscatter, radiation-pressure, and asymmetry properties from the coefficient series.
- Homogeneous spheres: RayleighMieQ() and LowFrequencyMieQ() accelerate calculations for particles with x≪1 by using simplifying assumptions or only the first few coefficient terms.These functions target the Rayleigh or low-frequency regime.
- Size distributions: PyMieScatt integrates optical properties over measured or generated size distributions and converts the resulting coefficients to common laboratory units.Mie_SD() accepts binned distributions, while Mie_Lognormal() generates k-modal lognormal distributions; the 10^-6 factor converts to Mm^-1.
- Scattering outputs: The package calculates angular scattering intensities, polarization-dependent quantities, scattering-matrix elements, and θ-space or q-space outputs for particle ensembles.Scattered-field calculations use S1 and S2, with default angular output from 0–180° in 0.5° increments unless constrained.
- Inverse and coated-particle calculations: Coated-sphere functions handle core and shell diameters and refractive indices, while the survey-iteration algorithm retrieves refractive index by refining candidate neighborhoods with forward Mie calculations.The survey identifies candidate regions from scattering and absorption arrays, and iteration produces distinct solutions within neighborhoods.