T-matrix method

The T-matrix method is a computational technique of light scattering by nonspherical particles originally formulated by P. C. Waterman (1928-2012) in 1965.^[1] The technique is also known as null field method and extended boundary technique method (EBCM). ^[2] In the method, matrix elements are obtained by matching boundary conditions for solutions of Maxwell equations.

Definition of the T-Matrix

The incident and scattered electric field are expanded into spherical vector wave functions (SVWF), which are also encountered in Mie scattering. They are the fundamental solutions of the vector Helmholtz equation and can be generated from the scalar fundamental solutions in spherical coordinates, the spherical Bessel functions of the first kind and the spherical Hankel Functions. Accordingly, there are two linearly independent sets of solutions denoted as $\mathbf{M}^1,\mathbf{N}^1$ and $\mathbf{M}^3,\mathbf{N}^3$ , respectively. They are also called regular and propagating SVWFs, respectively. With this, we can write the incident field as

$\mathbf{E}_{inc}= \sum_{n=1}^\infty \sum_{m=-n}^n a_{mn} \mathbf{M}^1_{mn}+ b_{mn} \mathbf{N}^1_{mn}.$

The scattered field is expanded into radiating SVWFs:

$\mathbf{E}_{scat}= \sum_{n=1}^\infty \sum_{m=-n}^n f_{mn} \mathbf{M}^3_{mn}+ g_{mn} \mathbf{N}^3_{mn}.$

The T-Matrix relates the expansion coefficients of the incident field to those of the scattered field.

$\begin{pmatrix} a_{mn}\\ b_{mn}\end{pmatrix} = T \begin{pmatrix} f_{mn} \\ g_{mn} \end{pmatrix}$

The T-Matrix is determined by the scatterer shape and material and for given incident field allows to calculate the scattered field.

Calculation of the T-Matrix

The standard way to actually calculate the T-Matrix method is the Null-Field Method, that relies on the Stratton-Chu equations.^[3] They basically state that the electromagnetic fields outside a given volume can be expressed as integrals over the surface enclosing the volume involving only the tangential components of the fields on the surface. If the observation point is located inside this volume, the integrals vanish.

By making use of the boundary conditions for the tangential field components on the scatterer surface $\mathbf{n} \times (\mathbf{E}_{scat} + \mathbf{E}_{inc}) = \mathbf{E}_{int}$ and $\mathbf{n} \times (\mathbf{H}_{scat} + \mathbf{H}_{inc}) = \mathbf{H}_{int}$ , where $\mathbf {n}$ is the normal vector to the scatterer surface, one can derive an integral representation of the scattered field in terms of the tangential components of the internal fields on the scatterer surface. A similar representation can be derived for the incident field.

By expanding the internal field in terms of SVWFs and exploiting their orthogonality on spherical surfaces, one arrives at an expression for the T-Matrix. Numerical codes for the evaluation of the T-Matrix can be found online .

References

↑ Waterman, P. C. "Matrix formulation of electromagnetic scattering." Proceedings of the IEEE 53, 805-812 (1965).
↑ M. I. Mishchenko, L. D. Travis, and D. W. Mackowski, T-matrix computations of light scattering by nonspherical particles: A review, J. Quant. Spectrosc. Radiat. Transfer, 55, 535-575 (1996).
↑ J.A. Stratton & L.J. Chu: Diffraction Theory of Electromagnetic Waves Phys. Rev., American Physical Society, 56 (1939)

This article is issued from Wikipedia - version of the 3/27/2015. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.