Journal
JOURNAL OF QUANTITATIVE SPECTROSCOPY & RADIATIVE TRANSFER
Volume 123, Issue -, Pages 153-168Publisher
PERGAMON-ELSEVIER SCIENCE LTD
DOI: 10.1016/j.jqsrt.2013.01.023
Keywords
Spheroids; T-matrix; Extended boundary-condition method; Light scattering
Categories
Funding
- Royal Society of New Zealand
Ask authors/readers for more resources
We propose, describe, and demonstrate a new numerically stable implementation of the extended boundary-condition method (EBCM) to compute the T-matrix for electromagnetic scattering by spheroidal particles. Our approach relies on the fact that for many of the EBCM integrals in the special case of spheroids, a leading part of the integrand integrates exactly to zero, which causes catastrophic loss of precision in numerical computations. This feature was in fact first pointed out by Waterman in the context of acoustic scattering and electromagnetic scattering by infinite cylinders. We have recently studied it in detail in the case of electromagnetic scattering by particles. Based on this study, the principle of our new implementation is therefore to compute all the integrands without the problematic part to avoid the primary cause of loss of precision. Particular attention is also given to choosing the algorithms that minimise loss of precision in every step of the method, without compromising on speed. We show that the resulting implementation can efficiently compute in double precision arithmetic the T-matrix and therefore optical properties of spheroidal particles to a high precision, often down to a remarkable accuracy (10(-10) relative error), over a wide range of parameters that are typically considered problematic. We discuss examples such as high-aspect ratio metallic nanorods and large size parameter (approximate to 35) dielectric particles, which had been previously modelled only using quadruple-precision arithmetic codes. (C) 2013 Elsevier Ltd. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available