Flexible implementation of the particle shape and internal inhomogeneity in the invariant imbedding T-matrix method

We report a new implementation of the invariant imbedding T-matrix (IITM) method based on a discrete spherical grid approach for representing the particle shape and internal inhomogeneity. The new version of the IITM (referred to as the IITM-discrete) improves the flexibility of the IITM-especially for inhomogeneous particles. It is much more convenient for specifying the particle morphology in the electromagnetic wave scattering simulations. Particle shape is represented by a series of discrete spherical layers ranging from the inscribed sphere to the circumscribed sphere. Spherical layers are discretized by the centroidal Voronoi tessellation (CVT) approach. The procedure of computing the U-matrix (the only shape-dependent module in the T-matrix program) is simplified upon using the gridded particle shape and refractive index information saved in an external file. The grid resolution is a key factor that determines the numerical accuracy and computational cost. Numerical tests of IITM-discrete show its compatibility with other light scattering methods. Using IITM-discrete, we found that the internal inhomogeneity could have large impact on dust optical properties.