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

Ayranci, Isil
Vaillon, Rodolphe
Selçuk, Nevin
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.


Modified Superformula Contours Optimized via Genetic Algorithms for Exponentially Converging 2D Solutions of MFIE
Guler, Sadri; Onol, Can; Ergül, Özgür Salih; Sever, Emrah; Dikmen, Fatih; Tuchkin, Yury A. (2017-05-25)
An infinitely smooth parametrical representation with derivatives of all orders is used, resulting into exponentially converging solutions of magnetic field integral equation (MFIE) in 2D either for TM or TE polarized excitations. A version of superformula modified for this purpose has been subject to optimization of its parameters via genetic algorithms to provide smooth parameterization for a desired boundary in two-dimensional problems. The organization of the MFIE kernel and convergence of the solution ...
Full-Wave Computational Analysis of Optical Chiral Metamaterials
Guler, Sadri; Solak, Birol; Gür, Uğur Meriç; Ergül, Özgür Salih (2017-09-27)
We present computational analysis of optical chiral metamaterials that consist of helical metallic elements. At optical frequencies, metals are modeled as penetrable objects with plasmonic properties. A rigorous implementation based on boundary element methods and the multilevel fast multipole algorithm is used for efficient and accurate analysis of three-dimensional structures. Numerical results demonstrate interesting polarization-rotating characteristics of such arrays with helical elements, as well as t...
Accurate solutions of scattering problems involving low-contrast dielectric objects with surface integral equations
Ergül, Özgür Salih (2007-11-16)
We present the stabilization of the surface integral equationsfor accurate solutions of scattering problems involvinglow-contrast dielectric objects. Unlike volume formulations,conventional surface formulations fail to provide accurateresults for the scatteredfields when the contrast of theobject is small. Therefore, surface formulations are requiredto be stabilized by extracting the nonradiating parts of theequivalent currents. In addition to previous strategies forthe stabilization, we introduce a n...
Fast and accurate solutions of electromagnetics problems involving lossy dielectric objects with the multilevel fast multipole algorithm
Ergül, Özgür Salih (2012-03-01)
Fast and accurate solutions of electromagnetic scattering problems involving lossy dielectric objects are considered. Problems are formulated with two recently developed formulations, namely, the combined-tangential formulation (CTF) and the electric and magnetic current combined-field integral equation (JMCFIE), and solved iteratively using the multilevel fast multipole algorithm (MLFMA). Iterative solutions and accuracy of the results are investigated in detail for diverse geometries, frequencies, and con...
Adaptive Control of Necklace States in a Photonic Crystal Waveguide
Yüce, Emre; Sokolov, Sergei; Bertolotti, Jacopo; Combrie, Sylvain; Lehoucq, Gaelle; De Rossi, Alfredo; Mosk, Allard P. (2018-10-01)
Resonant cavities with high quality factor and small mode volume provide crucial enhancement of light matter interactions in nanophotonic devices that transport and process classical and quantum information. The production of functional circuits containing many such cavities remains a major challenge, as inevitable imperfections in the fabrication detune the cavities, which strongly affects functionality such as transmission. In photonic crystal waveguides, intrinsic disorder gives rise to high-Q localized ...
Citation Formats
I. Ayranci, R. Vaillon, and N. Selçuk, “Performance of discrete dipole approximation for prediction of amplitude and phase of electromagnetic scattering by particles,” JOURNAL OF QUANTITATIVE SPECTROSCOPY & RADIATIVE TRANSFER, pp. 83–101, 2007, Accessed: 00, 2020. [Online]. Available: