In this article, a two-dimensional (2D) unconditionally stable finite-difference time-domain (FDTD) approach is proposed for graphene electromagnetic (EM) device simulation. The weighted Laguerre polynomials (WLPs) are utilized to resolve stability concerns, and graphene is modelled as a thin conductive layer incorporating the surface boundary condition (SBC) in WLP-FDTD scheme. The transmittance of EM signal propagating through two graphene layers is calculated for 0-10 THz to verify the effectiveness of the proposed method. The simulation results agree excellently with the results calculated from the analytical and other numerical models. The proposed SBC-WLP-FDTD method provides an alternative numerical approach to simulate graphene-like materials with improved computing efficiency.
2. Yu, X. and C. D. Sarris, "A perfectly matched layer for subcell FDTD and applications to the modeling of graphene structures," IEEE Antennas Wireless Propag. Lett., Vol. 11, 1080-1083, 2012.
3. Ahmed, I., E. H. Khoo, and E. Li, "Efficient modeling and simulation of graphene devices with the LOD-FDTD method," IEEE Microw. Wireless Compon. Lett., Vol. 23, No. 6, 306-308, Jun. 2013.
doi:10.1109/LMWC.2013.2258463
4. De Oliveira, R. M. S., N. R. N. M. Rodrigues, and V. Dmitriev, "FDTD formulation for graphene modeling based on piecewise linear recursive convolution and thin material sheets techniques," IEEE Antennas Wireless Propag. Lett., Vol. 14, 767-770, 2015.
doi:10.1109/LAWP.2014.2378174
5. Chen, W.-J., W. Shao, J. Quan, and S.-Y. Long, "Modeling of wave propagation in thin graphene sheets with WLP-FDTD method," Journal of Electromagnetic Waves and Applications, Vol. 30, No. 6, 780-787, Apr. 2016.
doi:10.1080/09205071.2016.1150210
6. Chen, W.-J. and S.-Y. Long, "Modeling of wave propagation in thin graphene sheets with 2-D ADE-WLP-FDTD method," 2016 IEEE International Conference on Microwave and Millimeter Wave Technology (ICMMT), Vol. 1, 434-436, Jun. 2016.
doi:10.1109/ICMMT.2016.7761799
7. Zhu, Q.-Y. and W.-J. Chen, "Modeling thin grapheme sheets with efficient 2-D WLP-FDTD method," 2017 Sixth Asia-Pacific Conference on Antennas and Propagation (APCAP), Vol. 1, Oct. 2017.
8. Nayyeri, V., M. Soleimani, and O. M. Ramahi, "Modeling graphene in the finite-difference time-domain method using a surface boundary condition," IEEE Trans. Antennas Propaga., Vol. 61, No. 8, 4176-4182, Aug. 2013.
doi:10.1109/TAP.2013.2260517
9. Hanson, G. W., "Dyadic Greens functions and guided surface waves for a surface conductivity model of graphene," J. Appl. Phys., Vol. 103, No. 6, Mar. 2008.
doi:10.1063/1.2891452
10. Chung, Y. S., T. K. Sarkar, B. H. Jung, and M. Salazar-Palma, "An unconditionally stable scheme for the finite-difference time-domain method," IEEE Trans. Microw. Theory Tech., Vol. 51, No. 3, 697-704, Mar. 2003.
doi:10.1109/TMTT.2003.808732
11. Chen, W.-J., W. Shao, J.-L. Li, and B.-Z. Wang, "Numerical dispersion analysis and key parameter selection in Laguerre-FDTD method," IEEE Microw. Wireless Compon. Lett., Vol. 23, No. 12, 629-631, Dec. 2013.
doi:10.1109/LMWC.2013.2283866