Numerical instability is usually observed when the propagation matrix method is used to calculate the reflectance and transmittance spectra for the thick one-dimensional inhomogeneous optical structures and media. To remove this numerical instability we applied two procedures, the normalization and the singular-value decomposition, for the propagation matrix and the matrix involved in calculating the matrix of reflection coefficients, respectively. Examples of a cholesteric liquid crystal and a helical structure of ferroelectric liquid crystals with a twist defect show that the modified propagation matrix method is able to accurately calculate the reflectance spectra for thick structures.
© 2008 Optical Society of America
Original Manuscript: May 27, 2008
Manuscript Accepted: June 25, 2008
Published: August 22, 2008
Zhao Lu, "Accurate calculation of reflectance spectra for thick one-dimensional inhomogeneous optical structures and media: stable propagation matrix method," Opt. Lett. 33, 1948-1950 (2008)