Search

Understanding the inherent optical properties (IOPs) of coccoliths and coccolithophores is important in oceanic radiative transfer simulations and remote sensing implementations. In this study, the invariant imbedding T-matrix method (II-TM) is employed to investigate the IOPs of coccoliths and coccolithophores. The Emiliania huxleyi ( Ehux ) coccolith and coccolithophore models are built based on observed biometric parameters including the eccentricity, the number of slits, and the rim width of detached coccoliths. The calcification state that specifies the amount of calcium of a single coccolith is critical in the determination of the size–volume/mass relationship (note, the volume/mass of coccoltihs at different calcification states are different although the diameters are the same). The present results show that the calcification state, namely, under-calcification, normal-calcification, or over-calcification, significantly influences the backscattering cross section and the phase matrix. Furthermore, the linear depolarization ratio of the light scattered by coccoliths is sensitive to the degree of calcification, and provides a potentially valuable parameter for interpreting oceanic remote sensing data. The phase function of an ensemble of randomly oriented coccolithophores has a similar pattern to that of individual coccoliths, but the forward scattering is dominant in the coccolithophores due to the large geometric cross sections. The linear depolarization ratio associated with coccolithophores is found to be larger than that for coccoliths as polarization is more sensitive to multiple scattering than the phase function. The simulated coccolithophore phase matrix numerical results are compared with laboratory measurements. For scattering angles larger than 100°, an increase of the phase function with respect to the scattering angle is confirmed based on the present coccolithophore model while the spherical approximation fails. [ABSTRACT FROM AUTHOR]

Abstract: The invariant imbedding T-matrix method (II-TM) is employed to compute the optical properties of randomly oriented ice crystals of various shapes including hexagonal columns, hollow columns, droxtals, bullet rosettes and aggregates. The II-TM is shown to be numerically stable and capable of obtaining the single-scattering properties of hexagonal ice crystals with size parameters up to 150. The 22° and 46° halo peaks in the phase function of compact hexagonal ice crystals begin to emerge at a size parameter of approximately 80 and tend to become insensitive to particle size as the corresponding size parameter approaches 150. Furthermore, the II-TM solutions are shown to be in agreement with their counterparts based on the discrete dipole approximation (DDA) method and the pseudo-spectral time-domain (PSTD) method. In addition, the accuracy of the improved geometric-optics method (IGOM) is examined for randomly oriented hexagonal ice crystal cases over a wide size-parameter range from the resonant to geometric-optics regimes. The II-TM is also used to study the effects of particle surface roughness and internal inclusions on the single-scattering properties of ice particles. [Copyright &y& Elsevier]

Abstract: The extended boundary condition method (EBCM) and invariant imbedding method (IIM) are two fundamentally different T-matrix methods for the solution of light scattering by nonspherical particles. The standard EBCM is very efficient but encounters a loss of precision when the particle size is large, the maximum size being sensitive to the particle aspect ratio. The IIM can be applied to particles in a relatively large size parameter range but requires extensive computational time due to the number of spherical layers in the particle volume discretization. A numerical combination of the EBCM and the IIM (hereafter, the EBCM+IIM) is proposed to overcome the aforementioned disadvantages of each method. Even though the EBCM can fail to obtain the T-matrix of a considered particle, it is valuable for decreasing the computational domain (i.e., the number of spherical layers) of the IIM by providing the initial T-matrix associated with an iterative procedure in the IIM. The EBCM+IIM is demonstrated to be more efficient than the IIM in obtaining the optical properties of large size parameter particles beyond the convergence limit of the EBCM. The numerical performance of the EBCM+IIM is illustrated through representative calculations in spheroidal and cylindrical particle cases. [Copyright &y& Elsevier]

Abstract: Three terms, “Waterman''s T-matrix method”, “extended boundary condition method (EBCM)”, and “null field method”, have been interchangeable in the literature to indicate a method based on surface integral equations to calculate the T-matrix. Unlike the previous method, the invariant imbedding method (IIM) calculates the T-matrix by the use of a volume integral equation. In addition, the standard separation of variables method (SOV) can be applied to compute the T-matrix of a sphere centered at the origin of the coordinate system and having a maximal radius such that the sphere remains inscribed within a nonspherical particle. This study explores the feasibility of a numerical combination of the IIM and the SOV, hereafter referred to as the IIM+SOV method, for computing the single-scattering properties of nonspherical dielectric particles, which are, in general, inhomogeneous. The IIM+SOV method is shown to be capable of solving light-scattering problems for large nonspherical particles where the standard EBCM fails to converge. The IIM+SOV method is flexible and applicable to inhomogeneous particles and aggregated nonspherical particles (overlapped circumscribed spheres) representing a challenge to the standard superposition T-matrix method. The IIM+SOV computational program, developed in this study, is validated against EBCM simulated spheroid and cylinder cases with excellent numerical agreement (up to four decimal places). In addition, solutions for cylinders with large aspect ratios, inhomogeneous particles, and two-particle systems are compared with results from discrete dipole approximation (DDA) computations, and comparisons with the improved geometric-optics method (IGOM) are found to be quite encouraging. [Copyright &y& Elsevier]

Abstract: The far field within the context of the Lorenz–Mie theory and the T-matrix formulation is usually expressed on the basis of the asymptotic properties of vector spherical waves. The radiation condition is taken into account by employing proper vector spherical functions as the expansion basis of the scattered field. The asymptotic behavior of the Hankel function is obtained from differential equations. The asymptotic far field can also be obtained from the Kirchhoff surface integral equation, in which the radiation condition has been implemented when it is derived from the Maxwell equations. This note is to present an explicit establishment of the relationship between the asymptotic far field and the near field in the Lorenz–Mie theory and the T-matrix formulation through the Kirchhoff surface integral. [Copyright &y& Elsevier]

Abstract: The many-body iterative T-matrix (MBIT) method is based on the T-matrix method and focuses on scatterers with large aspect ratios. In this study, the MBIT method is extended so that it can be applied to an infinite scatterer with a 1-D periodic finite structure. A semi-analytical solution of infinite scatterers is obtained. In addition, the MBIT method is theoretically extended to include scatterers without axial rotational symmetry. The direct validation of the present method is performed by comparing the analytic solution for an infinite cylinder and the numerically accurate discrete dipole approximation method whereas indirect validation is conducted by using asymptotic solutions. Furthermore, the study also demonstrates that scattering properties of infinite scatterers can be used as a substitute for finite scatterers with large aspect ratios. [Copyright &y& Elsevier]