Dimension-variable invariant imbedding (DVIIM) T-matrix computational method for the light scattering simulation of atmospheric nonspherical particles

The invariant imbedding (IIM) T-matrix method has shown great potential in light scattering field. However, the T-matrix need to be calculated through the matrix recurrence formula derived from the Helmholtz equation, thus its computational efficiency is much lower than Extended Boundary Condition Method (EBCM). To alleviate this problem, the Dimension-Variable Invariant Imbedding (DVIIM) T-matrix method is presented in this paper. Compared with the traditional IIM T-matrix model, the dimensions of the T-matrix and relevant matrices are gradually increasing as the iteration performed step by step, thus the unnecessary operations of large matrices can be avoided in early iterations. To optimally determine the dimension of these matrices in each iterative calculation, the spheroid-equivalent scheme (SES) is also proposed. The effectiveness of the DVIIM T-matrix method is validated from the modeling accuracy and calculation efficiency. The simulation results show that compared with traditional T-matrix method, its modeling efficiency can be improved notably, especially for the particles with large size and aspect ratio, where for the spheroid with a aspect ratio of 0.5, the computational time is cut down by 25%. Though the dimension of the T matrix is cut down in the early iterations, the computational precision of DVIIM T-matrix model is not decreased notably, and a good agreement is achieved between the calculation results of DVIIM T-matrix method, IIM T-matrix method and other well-validated models (like EBCM and DDACSAT), where the relative errors of the integral scattering parameters (e.g., extinction, absorption, scattering cross sections) are generally less than 1%.