A new physical-geometric optics method is developed to compute the single-scattering properties of faceted particles. It incorporates a general absorption vector to accurately account for inhomogeneous wave effects, and subsequently yields the relevant analytical formulas effective and computationally efficient for absorptive scattering particles. A bundle of rays incident on a certain facet can be traced as a single beam. For a beam incident on multiple facets, a systematic beam-splitting technique based on computer graphics is used to split the original beam into several sub-beams so that each sub-beam is incident only on an individual facet. The new beam-splitting technique significantly reduces the computational burden. The present physical-geometric optics method can be generalized to arbitrary faceted particles with either convex or concave shapes and with a homogeneous or an inhomogeneous (e.g., a particle with a core) composition. The single-scattering properties of irregular convex homogeneous and inhomogeneous hexahedra are simulated and compared to their counterparts from two other methods including a numerically rigorous method.

No abstract

The optical properties of diatom chains in the ocean are studied based on a combination of the many-body iterative T-matrix (MBIT) method and an improved implementation of the ray-by-ray (RBR) geometric optics method. The MBIT, a numerically accurate method, is advantageous for scatterers with linear cells. In contrast to other popular geometric optics methods, the RBR, an approximate method, considers the interference of all outgoing rays. The two methods are verified in comparison with benchmark simulations. The simulation results of diatom chains in a wide size range can be obtained with either or both methods, and each can be applied to any optically soft particle, i.e., in the case when relative refractive index approaches unity.

No abstract

No abstract

No abstract

Integrated and differential optical properties of a single particle, such as the scattering, absorption, and extinction cross sections, single scattering albedo, asymmetry factor, and scattering phase matrix, are derived from electromagnetic scattering theory. This process depends on microphysical inputs which include particle shape, refractive index, aspect ratio, and size parameter. In this work, we use the invariant imbedding T-matrix method (IITM) to derive analytic expressions for Jacobians of these optical properties with respect to the input parameters. These IITM-derived Jacobians for spheroids, cylinders, and hexagonal prisms are validated by comparison with results calculated with the extended boundary condition method (EBCM) and further validated using finite-difference estimates. We examine the dependencies of these Jacobians as functions of the input microphysical parameters, focusing again on spheroids, cylinders, and hexagonal prisms.

The scalar radiative transfer equation in the presence of thermal radiation source is solved in detail, using the adding-doubling method; Planck functions within any given layer are assumed to possess constant, linear, or exponential parameterizations with optical thickness. The radiance profile in any zenith direction is calculated directly in terms of matrix inversions. The inputs to the model are the inherent optical properties (layer total single-scattering albedos, scattering phase functions, and optical thickness) along with temperature and altitude profiles, and the top of the atmosphere and ground surface boundary conditions. The algorithm is implemented in a state-of-the-art MATLAB program, with the cosmic microwave background as the source at the upper boundary and a Lambertian surface reflection at the lower boundary. The simulations are validated against the VLIDORT discrete ordinate radiative transfer model. Results are compared in detail for cases with linear and exponential Planck function parameterizations.