Time-domain simulation of electromagnetic field using a symplectic integrator

A new high-order-accurate method for time-domain simulation of a electromagnetic held is proposed. This method is analogous to the one used in classical mechanics for continuous system analysis. The symplectic integrator in classical mechanics is utilized to calculate the time dependence of the held. Using the same grid spacing and time step, it is demonstrated that the proposed method is more accurate than the conventional finite-difference time-domain (FDTD) method for a two-dimensional problem. This method is also applicable to full-vectorial Maxwell's equations in three dimensions.