Performance of discrete dipole approximation for prediction of amplitude and phase of electromagnetic scattering by particles

Electromagnetic scattering provides useful signatures for nonintrusive particle characterization. Scattered wave which carries characteristic information about particles is identified completely by its intensity, polarization state and phase. Recent developments in measurement techniques have enabled measurement of phase of the scattered wave which is a source of additional information about particles. In the present study, accuracy of discrete dipole approximation (DDA) in predicting amplitude and phase of scattered wave is investigated via publicly available DDSCAT code by Draine and Flatau, which is a well-established tool for DDA and has found wide range of applications in the literature due to its flexibility. DDSCAT routine is modified to enable accurate computation of phase of complex amplitude scattering matrix (ASM) elements as well as their magnitude. DDA method was implemented by using lattice dispersion relation for dipole polarizabilities, generalized prime factor algorithm for fast-Fourier transformation and pre-conditioned bi-conjugate gradient method with stabilization for the solution of the complex linear system of equations. Accuracy of ASM elements predicted by DDA is assessed on single sphere problems with various size parameters and refractive indices by validation against Mie theory solutions. Excellent agreement between predictions and exact solutions proves the reliability of the modified DDSCAT code for prediction of amplitude and phase of scattered electromagnetic wave. Applicability conditions and requirements of the present DDA application to ensure accurate prediction of complete set of scattering parameters are mapped for single spheres, on an extensive domain of size parameters and refractive indices. A correlation is presented to estimate the magnitude and phase errors associated with given size parameter, refractive index and cubic lattice subdivision. Assessment of computational time requirements for different optical constants shows that implementation of DDA with the present specifications is unfeasible for size parameters larger than 4 when Re(m)>2 and Im(m)<0.1 at the same time, due to slow convergence rate.