A high-order accurate scheme for Maxwell's equations with a Generalized Dispersive Material (GDM) model and material interfaces

Abstract A high-order accurate scheme for solving the time-domain dispersive Maxwell's equations and material interfaces is described. Maxwell's equations are solved in second-order form for the electric field. A generalized dispersive material (GDM) model is used to represent a general class of linear dispersive materials and this model is implemented in the time-domain with the auxiliary differential equation (ADE) approach. The interior updates use our recently developed second-order and fourth-order accurate single-stage three-level space-time finite-difference schemes, and this paper extends these schemes to treat interfaces between different dispersive materials. Composite overlapping grids are used to treat complex geometry, with Cartesian grids generally covering most of the domain and local conforming grids representing curved boundaries and interfaces. Compatibility conditions derived from the interface jump conditions and governing equations are used to derive accurate numerical interface conditions that define values at ghost points. Although some compatibility conditions couple the equations for the ghost points in tangential directions due to mixed-derivatives, it is shown how to decouple the equations to avoid solving a larger system of equations for all ghost points on the interface. The stability of the interface approximations is studied with mode analysis and it is shown that the schemes retain close to a CFL-one time-step restriction. Numerical results are presented in two and three space dimensions to confirm the accuracy and stability of the schemes. The schemes are verified using exact solutions for a planar interface, a disk in two dimensions, and a solid sphere in three dimensions.

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

[2]  A. Kildishev,et al.  Designing optimal nanofocusing with a gradient hyperlens , 2017 .

[3]  Jan S. Hesthaven,et al.  High-order accurate methods in time-domain computational electromagnetics: A review , 2003 .

[4]  Jeffrey W. Banks,et al.  Upwind schemes for the wave equation in second-order form , 2012, J. Comput. Phys..

[5]  Linjie Zhou,et al.  Miniature Multilevel Optical Memristive Switch Using Phase Change Material , 2019, ACS Photonics.

[6]  W. Henshaw,et al.  Composite overlapping meshes for the solution of partial differential equations , 1990 .

[7]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[8]  Heinz-Otto Kreiss,et al.  Difference Approximations for the Second Order Wave Equation , 2002, SIAM J. Numer. Anal..

[9]  Sergey I. Bozhevolnyi,et al.  Nanofocusing of electromagnetic radiation , 2013, Nature Photonics.

[10]  William D. Henshaw,et al.  Parallel computation of three-dimensional flows using overlapping grids with adaptive mesh refinement , 2008, J. Comput. Phys..

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

[12]  Qing Huo Liu,et al.  Higher-Order Numerical Methods for Transient Wave Equations , 2003 .

[13]  Alexander V. Kildishev,et al.  A high-order accurate scheme for Maxwell's equations with a generalized dispersive material model , 2019, J. Comput. Phys..

[14]  M. Yoshita,et al.  High-precision group-delay dispersion measurements of optical fibers via fingerprint-spectral wavelength-to-time mapping , 2016 .

[15]  C. Tropea,et al.  Light Scattering from Small Particles , 2003 .

[16]  C. Balanis.  Advanced Engineering Electromagnetics , 1989 .

[17]  Vladimir M. Shalaev,et al.  Examining nanophotonics for integrated hybrid systems: a review of plasmonic interconnects and modulators using traditional and alternative materials [Invited] , 2015 .

[18]  Mario Miscuglio,et al.  All-optical nonlinear activation function for photonic neural networks [Invited] , 2018, Optical Materials Express.

[19]  M. Gorodetsky,et al.  Dissipative Kerr solitons in optical microresonators , 2015, Science.

[20]  Yuan Cheng,et al.  Near-field optics on flatland: from noble metals to van der Waals materials , 2019, Advances in Physics: X.

[21]  Wolfgang Porod,et al.  Roadmap on all-optical processing , 2019, Journal of Optics.

[22]  Jeffrey W. Banks,et al.  High-order upwind schemes for the wave equation on overlapping grids: Maxwell's equations in second-order form , 2018, J. Comput. Phys..

[23]  Phaedon Avouris,et al.  Graphene acoustic plasmon resonator for ultrasensitive infrared spectroscopy , 2019, Nature Nanotechnology.

[24]  William D. Henshaw,et al.  A High-Order Accurate Parallel Solver for Maxwell's Equations on Overlapping Grids , 2005, SIAM J. Sci. Comput..