Abstract

I present enhanced R-matrix algorithms for analysis of general multilayered diffraction gratings. The previous R-matrix algorithms are enhanced in three aspects: computational efficiency, numerical stability, and application of half R-matrix in addition to full and quarter R-matrix recursions. On the basis of the eigensolutions of rigorous coupled-wave analysis, the enhanced R-matrix algorithms deal with eigen-submatrices directly and bypass the auxiliary layer R matrix. Such exclusion of a layer matrix leads to improvements in efficiency and algorithm robustness particularly for zero or small layer thickness relative to wavelength. Application of the enhanced algorithms to grating diffraction is exploited especially for the half and quarter R-matrix recursions. Comparison of various R-matrix algorithms via a table of flop counts shows that the enhanced algorithms are more efficient apart from being well conditioned.

© 2006 Optical Society of America

Formulation and comparison of two recursive matrix algorithms for modeling layered diffraction gratings

Lifeng Li
J. Opt. Soc. Am. A 13(5) 1024-1035 (1996)

Bremmer series, R-matrix propagation algorithm, and numerical modeling of diffraction gratings

Lifeng Li
J. Opt. Soc. Am. A 11(11) 2829-2836 (1994)

Note on the S-matrix propagation algorithm

Lifeng Li
J. Opt. Soc. Am. A 20(4) 655-660 (2003)

Dispersion in photonic media and diffraction from gratings: a different modal expansion for the R-matrix propagation technique

J. Merle Elson and P. Tran
J. Opt. Soc. Am. A 12(8) 1765-1771 (1995)

Differential theory of gratings: extension to deep gratings of arbitrary profile and permittivity through the R-matrix propagation algorithm

F. Montiel and M. Neviere
J. Opt. Soc. Am. A 11(12) 3241-3250 (1994)