By applying the matrix-exponential operator technique to the radiative-transfer equation in discrete form, new analytical solutions are obtained for the transmission and reflection matrices in the limiting cases x ≪ 1 and x ≫ 1, where x is the optical depth of the layer. Orthogonality of the eigenvectors of the matrix exponential apparently yields new conditions for determining Chandrasekhar’s characteristic roots. The exact law of reflection for the discrete eigenfunctions is also obtained. Finally, when used in conjunction with the doubling method, the matrix exponential should result in reductions in both computation time and loss of precision.
© 1981 Optical Society of America
P. C. Waterman, "Matrix-exponential description of radiative transfer," J. Opt. Soc. Am. 71, 410-422 (1981)