Vol. 14

Latest Volume
All Volumes
All Issues

Efficient 4x4 Propagation Matrix Method Using a Fourth-Order Symplectic Integrator for the Optics of One-Dimensional Continuous Inhomogeneous Materials

By Zhao Lu
Progress In Electromagnetics Research Letters, Vol. 14, 1-9, 2010


Understanding the propagation of light in continuous inhomogeneous materials is important to design optical structures and devices. To have accurately numerical calculations Berreman's 4×4 propagation matrix method is generally used, and layer approximation, i.e., the whole one-dimensional continuous inhomogeneous material is divided into many small homogeneous layers, is assumed. However, this layer approximation is only correct up to the second-order of the layer thickness. To efficiently solve Berreman's first-order differential equation, a simple fourth-order symplectic integrator is presented. The efficiency of the fourth-order symplectic integrator was studied for a cholesteric liquid crystal. Numerical results of reflectance spectra show that the fourth-order symplectic integrator is highly efficient in contrast to the extensively used fast 4×4 propagation matrix.


Zhao Lu, "Efficient 4x4 Propagation Matrix Method Using a Fourth-Order Symplectic Integrator for the Optics of One-Dimensional Continuous Inhomogeneous Materials," Progress In Electromagnetics Research Letters, Vol. 14, 1-9, 2010.


    1. Berreman, D. W., "Optics in stratified and anisotropic media: 4×4 matrix formulation," J. Opt. Soc. Am., Vol. 62, 502-510, 1972.

    2. Abdulhalim, I., "Analytic propagation matrix method for linear optics of arbitrary biaxial layered media," J. Opt. A: Pure Appl. Opt., Vol. 1, 646-653, 1999.

    3. Lu, Z., "Accurate and efficient calculation of light propagation in one-dimensional inhomogeneous anisotropic media through extrapolation ," J. Opt. Soc. Am. A, Vol. 24, 236-242, 2007.

    4. Wohler, H., G. Hass, M. Fritsch, and D. A. Mlynski, "Faster 4×4 matrix method for uniaxial inhomogeneous media," J. Opt. Soc. Am. A, Vol. 5, 1554-1557, 1988.

    5. Hawkeye, M. M. and M. J. Brett, "Narrow bandpass optical filters fabricated with one-dimensionally periodic inhomogeneous thin films," J. Appl. Phys., Vol. 100, 044322, 2006.

    6. Brett, M. J., M, and M. Hawkeye, "New materials at a glance," Science, Vol. 319, 1192-1193, 2008.

    7. Suzuki, M., Quantum Monte Carlo Methods in Condensed Matter Physics, World Scientific Pub Co Inc, 1994.

    8. Shankar, R., Principles of Quantum Mechanics, 2nd Ed., Springer, 1994.

    9. Chin, S. A., "Symplectic integrators from composite operator factorizations ," Phys. Lett. A, Vol. 226, 344-348, 1997.

    10. Chin, S. A. and D. W. Kidwell, "Higher-order force gradient symplectic algorithms," Phys. Rev. E, Vol. 62, 8746-8752, 2000.

    11. Chin, S. A. and C. R. Chen, "Gradient symplectic algorithms for solving the SchrÄodinger equation with time-dependent potentials," J. Chem. Phys., Vol. 117, 1409-1415, 2002.

    12. Nehring, J., "Light propagation and reflection by absorbing cholesteric layers," J. Chem. Phys., Vol. 75, 4326-4337, 1981.

    13. Lu, Z., "Accurate calculation of reflectance spectra for thick one-dimensional inhomogeneous optical structures and media: Stable propagation matrix method," Opt. Lett., Vol. 33, 1948-1950, 2008.