A fast directional algorithm for high-frequency electromagnetic scattering

This paper is concerned with the fast solution of high-frequency electromagnetic scattering problems using the boundary integral formulation. We extend the O(N log N) directional multilevel algorithm previously proposed for the acoustic scattering case to the vector electromagnetic case. We also detail how to incorporate the curl operator of the magnetic field integral equation into the algorithm. When combined with a standard iterative method, this results in an almost linear complexity solver for the combined field integral equations. In addition, the butterfly algorithm is utilized to compute the far field pattern and radar cross section with O(N log N) complexity.     

Characterization of the light scattering by ensembles of randomly distributed soot aggregates

Soot particles formed in combustion processes commonly exist in the form of ensembles of randomly distributed aggregates of small, nearly spherical monomers. In this paper, these randomly distributed aggregates are numerically generated using a combination of the cluster–cluster aggregation algorithm with the Monte Carlo method. Moreover, an efficient and accurate numerical method is presented for characterizing the light scattering by these complex soot particles illuminated by plane wave and Gaussian beam. This method exploits the unique features of the hybrid finite element-boundary integral method and, more importantly, the unique features of soot aggregates. It is designed in such a manner that it first decomposes the original problem into many sub-regions, where each primary particle is regarded as a sub-region, and then it employs the edge-based finite element method to deal with each sub-region. The sub-regions communicate through the near-field Green’s function. To reduce computational burdens, an iterative domain decomposition method in combination with parallel conjugate gradient method is adopted to solve the coupling system of equations. As an illustration, we present some of our preliminary numerical results. The results are expected to provide useful insights into the optical properties of soot particles formed in combustion processes.     

Investigation of resonant properties of metal core–shell nanoparticles using T-matrix calculations

Using T-matrix method, plasmon resonance properties of metal core–shell nanoparticles are systematically investigated. It is shown that dielectric/metal core–shell structure may be excited stronger at resonance than metal/dielectric and hetero-metal ones, but the resonance states are extremely sensible to the layers thickness. For three-layer nanospheres, resonance properties will be dominated by a sub-10 nm silver outermost shell, while only weakened by a silica one. Finally, tiny eccentric distance between the centers of core and shell in eccentric two-layer nanoparticles may fundamentally change the resonance mode of nanoparticles, and results in higher local electrical field enhancement than concentric nanospheres.     

Observations and calculations of two-dimensional angular optical scattering (TAOS) patterns of a single levitated cluster of two and four microspheres

We use single bi-sphere particles levitated in an electrodynamic balance to record two-dimensional angular scattering patterns at different angles of the coordinate system of the aggregate relative to the incident laser beam. Due to Brownian motion the particle covers the whole set of possible angles with time and allows to select patterns with high symmetry for analysis. These are qualitatively compared to numerical calculations. A small cluster of four spheres shows complex scattering patterns, comparison with computations suggest a low compactness for these clusters. An experimental procedure is proposed for studying restructuring effects occurring in mixed particles upon evaporation.     

Incomplete active polarimetry: Measurement of the block-diagonal scattering matrix

The problem of measurement of block-diagonal scattering matrix as being part of the so-called incomplete or partial Mueller polarimetry is addressed. It has been shown that non-zero elements of the scattering matrix with block-diagonal structure can be measured using only two elliptical polarizations of input light. Optimal pairs of these polarizations are derived. Utilization of the optimal input polarizations results in halving measurement time and decreasing measurement errors by 30%.     

Preconditioning techniques for iterative solvers in the Discrete Sources Method

Different preconditioning techniques for the iterative method MinRes as solver for the Discrete Sources Method (DSM) are presented. This semi-analytical method is used for light scattering computations by particles in the Mie scattering regime. Its numerical schema includes a linear least-squares problem commonly solved using the QR decomposition method. This could be the subject of numerical difficulties and instabilities for very large particles or particles with extreme geometry. In these cases, we showed that iterative methods with preconditioning techniques can provide a satisfying solution.In our previous paper, we studied four different iterative solvers (RGMRES, BiCGStab, BiCGStab(l), and MinRes) considering the performance and the accuracy of a solution. Here, we study several preconditioning techniques for the MinRes method for a variety of oblate and prolate spheroidal particles of different size and geometrical aspect ratio. Using preconditioning techniques we highly accelerated the iterative process especially for particles with a higher aspect ratio.     

Cation disorder in NaW2O6+δ·nH2−zO post-ion exchange with K, Rb, Sr, and Cs

The structure of the defect pyrochlore NaW₂O_6+δ·nH_2−zO after ion exchange with K, Rb, Sr or Cs for Na has been investigated using thermal analysis, solid-state nuclear magnetic resonance, laboratory X-ray and neutron diffraction methods. Neutron diffraction studies show that both the A-type cations (Na⁺, K⁺, Rb⁺, and/or Cs⁺) and the water molecules reside within the channels that form in the 111 direction of the W₂O₆ framework and that these strongly interact. The analytical results suggest that the water and A-type cations compete for space in the tunnels within the W₂O₆ pyrochlore framework, with the total number of water molecules and cations being approximately constant in the six samples investigated. The interplay between the cations and water explains the non-linear dependence of the a lattice parameter on the choice of cation. It appears that the ion-exchange capacity of the material will be controlled by the amount of water initially present in the sample.     

Analysis of the radiative transfer equation with highly assymetric phase function

This paper considers a scalar radiative transfer problem with high scattering anisotropy. Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal, which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2+ is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.