A symplectic FDTD algorithm for the simulations of lossy dispersive materials

Abstract A high-order symplectic finite-difference time-domain (SFDTD) algorithm, based on the matrix splitting, the symplectic integrator propagator and the auxiliary differential equation (ADE) technique, is presented. The algorithm is applied to the lossy Lorentz–Drude dispersive model. With the rigorous and artful formula derivation, the detailed formulations are provided. Moreover, the present algorithm can also be applicable to the lossy Drude and Lorentz dispersive model in a straightforward manner. An excellent agreement is achieved between the SFDTD-calculated and exact theoretical results when calculating the transmission coefficient in simulation of metal films. Focusing action of the matched left-handed materials (LHMs) slab is also achieved as the second example in the two-dimensional space. Numerical results for a more realistic structure, the simulation of periodic arrays of silver split-ring resonators (SRRs) using the Drude dispersion model, are also included, and the results agree well with those obtained by the finite element method (FEM).

[1]  Lei Wu,et al.  High-order unified symplectic FDTD scheme for the metamaterials , 2012, Comput. Phys. Commun..

[2]  Thomas Garm Pedersen,et al.  Optimization of imprintable nanostructured a-Si solar cells: FDTD study. , 2013, Optics express.

[3]  Levent Sevgi,et al.  Realistic surface modeling for a finite-difference time-domain wave propagator , 2003 .

[4]  Dennis M. Sullivan,et al.  Frequency-dependent FDTD methods using Z transforms , 1992 .

[5]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[6]  Marc Lamy de la Chapelle,et al.  Improved analytical fit of gold dispersion: Application to the modeling of extinction spectra with a finite-difference time-domain method , 2005 .

[7]  R. Ruth,et al.  Fourth-order symplectic integration , 1990 .

[8]  D. Yevick,et al.  Generalized propagation techniques-application to semiconductor rib waveguide Y-junctions , 1990, IEEE Photonics Technology Letters.

[9]  A. Kildishev,et al.  Optical Dispersion Models for Time-Domain Modeling of Metal-Dielectric Nanostructures , 2011, IEEE Transactions on Magnetics.

[10]  G. Strang,et al.  An Analysis of the Finite Element Method , 1974 .

[11]  Peng Liu,et al.  A Compact Symplectic High-Order Scheme for Time-Domain Maxwell's Equations , 2010, IEEE Antennas and Wireless Propagation Letters.

[12]  Weng Cho Chew,et al.  A comprehensive study for the plasmonic thin-film solar cell with periodic structure. , 2010, Optics express.

[13]  Yuzo Yoshikuni,et al.  A three-dimensional fourth-order finite-difference time-domain scheme using a symplectic integrator propagator , 2001 .

[14]  P. Bevelacqua,et al.  Finite-difference and pseudospectral time-domain methods applied to backward-wave metamaterials , 2004, IEEE Transactions on Antennas and Propagation.

[15]  Zhixiang Huang,et al.  Survey on Symplectic Finite-Difference Time-Domain Schemes for Maxwell's Equations , 2008, IEEE Transactions on Antennas and Propagation.

[16]  R. J. Luebbers,et al.  Piecewise linear recursive convolution for dispersive media using FDTD , 1996 .