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.


Far-field Lorenz-Mie scattering in an absorbing host medium. II: Improved...

Mishchenko, M., 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.

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.

PDF of Publication: 
Download from publisher's website.
Research Program: 
Radiation Science Program (RSP)