The system of coupled-wave equations for electromagnetic calculations of lamellar gratings is transformed to a new operator–vector form. The numerical procedure is based on truncation of the transformed system and proves to be stable, to be free of ill conditioning, and to preserve the power-conservation requirement for a lossless dielectric with high accuracy. To execute the procedure a very compact matlab-based program is developed, and numerical simulations for thick intrinsic silicon gratings are performed. Zero-reflectance phenomena at normal incidence for both TE and TM polarizations are studied. The ratios of the grating dimensions to be wavelengths at which these anomalies occur are found numerically. It is shown that by keeping the period- and slot-width-to-wavelength ratios constant and by increasing the slot depth one can repeat the anomalies. An antiblazing property at oblique incidence is also considered. The connection with recent directional polarized-emission experiments on intrinsic silicon gratings is discussed.
© 1994 Optical Society of America
Original Manuscript: June 14, 1993
Revised Manuscript: September 22, 1993
Published: July 20, 1994
S. Hava, M. Auslender, and D. Rabinovich, "Operator approach to electromagnetic coupled-wave calculations of lamellar gratings: infrared optical properties of intrinsic silicon gratings," Appl. Opt. 33, 4807-4813 (1994)