TY - JOUR

T1 - Uniform-over-size approximation of the internal fields for scatterers with low refractive-index contrast

AU - Inzhevatkin, Konstantin G.

AU - Yurkin, Maxim A.

N1 - Funding Information: We thank Alexander Konoshonkin for insightful discussions about geometric optics and Supercomputing center of the Novosibirsk State University for provided computational resources. The work was supported by the Russian Science Foundation (Grant No. 18-12-00052). Publisher Copyright: © 2021

PY - 2022/1

Y1 - 2022/1

N2 - We improved the Wentzel-Kramers-Brillouin (WKB) approximation for calculating the electric field inside a scatterer. The new method, named WKBr, additionally takes into account refraction (rotation) of the incident rays at the particle-medium interface. For large particles the WKBr is equivalent to the geometric optics in the limit of relative refractive index m approaching 1. Thus, the resulting accuracy is determined by (m − 1) uniformly for all particle sizes, in contrast to the original WKB, which accuracy deteriorates with increasing size. We extensively studied several versions WKBr for the case of a homogeneous sphere. These versions additionally account for Fresnel transmission and wavefront focusing for rays entering the scatterer, combination of intersecting rays, and the vanishing of the electric field in the shadow region. While keeping the same order of errors with varying (m − 1), they significantly reduce the errors for moderate values of m (from 1.1 to 1.3). Further, we tested the WKB and WKBr as an initial guess in the iterative solution of the linear system in the framework of the discrete dipole approximation (DDA). It allowed us to accelerate the DDA solution without compromising the final accuracy. The acceleration becomes more pronounced with increasing size, while the WKBr becomes clearly superior to the WKB. Moreover, the new initialization of the iterative solver extends the practical applicability of the DDA to size parameter equal to 250, 160, and 140 for m = 1.1, 1.2, and 1.3, respectively. They could not be reached previously due to the stagnation of the iterative solution. We expect similar DDA acceleration for particles of other shapes. One can easily test it for the case of the WKB, since it is already implemented in the open-source ADDA code.

AB - We improved the Wentzel-Kramers-Brillouin (WKB) approximation for calculating the electric field inside a scatterer. The new method, named WKBr, additionally takes into account refraction (rotation) of the incident rays at the particle-medium interface. For large particles the WKBr is equivalent to the geometric optics in the limit of relative refractive index m approaching 1. Thus, the resulting accuracy is determined by (m − 1) uniformly for all particle sizes, in contrast to the original WKB, which accuracy deteriorates with increasing size. We extensively studied several versions WKBr for the case of a homogeneous sphere. These versions additionally account for Fresnel transmission and wavefront focusing for rays entering the scatterer, combination of intersecting rays, and the vanishing of the electric field in the shadow region. While keeping the same order of errors with varying (m − 1), they significantly reduce the errors for moderate values of m (from 1.1 to 1.3). Further, we tested the WKB and WKBr as an initial guess in the iterative solution of the linear system in the framework of the discrete dipole approximation (DDA). It allowed us to accelerate the DDA solution without compromising the final accuracy. The acceleration becomes more pronounced with increasing size, while the WKBr becomes clearly superior to the WKB. Moreover, the new initialization of the iterative solver extends the practical applicability of the DDA to size parameter equal to 250, 160, and 140 for m = 1.1, 1.2, and 1.3, respectively. They could not be reached previously due to the stagnation of the iterative solution. We expect similar DDA acceleration for particles of other shapes. One can easily test it for the case of the WKB, since it is already implemented in the open-source ADDA code.

KW - Discrete dipole approximation

KW - Geometric optics

KW - Internal field

KW - Iterative solver

KW - Light scattering

KW - Wentzel-Kramers-Brillouin approximation

UR - http://www.scopus.com/inward/record.url?scp=85118856249&partnerID=8YFLogxK

U2 - 10.1016/j.jqsrt.2021.107965

DO - 10.1016/j.jqsrt.2021.107965

M3 - Article

AN - SCOPUS:85118856249

VL - 277

JO - Journal of Quantitative Spectroscopy and Radiative Transfer

JF - Journal of Quantitative Spectroscopy and Radiative Transfer

SN - 0022-4073

M1 - 107965

ER -