A recently developed FORTRAN program computing far-field optical observables for spherical particles in an absorbing medium has exhibited numerical instability arising when the product of the particle vacuum size parameter and the imaginary part of the refractive index of the host becomes sufficiently large. We offer a simple analytical explanation of this instability and propose a compact numerical algorithm for the stable computation of Lorenz–Mie coefficients based on an upward recursion formula for spherical Hankel functions of a complex argument. Extensive tests confirm an excellent accuracy of this algorithm approaching machine precision. The improved public-domain FORTRAN program is available at https:// www.giss.nasa.gov/staff/mmishchenko/Lorenz-Mie.html.
Far-field Lorenz-Mie scattering in an absorbing host medium. II: Improved stability of the numerical algorithm
Mishchenko, M.I., J.M. Dlugach, J.A. Lock, and M.A. Yurkin (2018), Far-field Lorenz-Mie scattering in an absorbing host medium. II: Improved stability of the numerical algorithm, J. Quant. Spectrosc. Radiat. Transfer, 217, 274-277, doi:10.1016/j.jqsrt.2018.05.034.
Abstract
PDF of Publication
Download from publisher's website
Research Program
Radiation Science Program (RSP)
Mission
ACE
Disclaimer: This material is being kept online for historical purposes. Though accurate at the time of publication, it is no longer being updated. The page may contain broken links or outdated information, and parts may not function in current web browsers. Visit https://espo.nasa.gov for information about our current projects.