The Computational Optimization of the Invariant Imbedding T Matrix Method for the Particles with N-Fold Symmetry

Abstract

The invariant imbedding T-matrix (IIM T-matrix) model is regarded as one of the most promising models for calculating the scattering parameters of non-spherical particles. However, the IIM T-matrix model needs to be iterated along the radial direction when calculating the T-matrix, which involves complex calculations such as matrix inversion and multiplication. Therefore, how to improve its computational efficiency is an important problem to be solved. Focused on particles with N-fold symmetric geometry, this paper deduced the symmetry in the calculation process of the IIM T-matrix model, derived the block iteration scheme of the T-matrix, and contracted the IIM T-matrix program for particles with N-fold symmetric geometry. Discrete Dipole Approximation (DDA) and Geometrical Optics Approximation (IGOA) were employed to verify the accuracy of the improved IIM T-matrix model. The results show that the six phase matrix elements (P11, P12/P11, P22/P11, P33/P11, P34/P11 and P44/P11) calculated by our model are in good agreement with other models. The computational efficiency of the improved IIM T-matrix model was further investigated. As demonstrated by the results, the computational efficiency for the particles with N-fold symmetry improved by nearly 70% with the improvement of the symmetry of U matrix and T matrix. In conclusion, the improved model can remarkably reduce the calculation time while maintaining high accuracy.