STABILITY AND NUMERICAL DISPERSION OF SYMPLECTIC 4TH-ORDER TIME-DOMAIN SCHEMES FOR OPTICAL-FIELD SIMULATION

The use of a more accurate scheme is effective in reducing the require d memory resources in the explicit time-domain simulation of optical f ield propagation. A promising technique is the application of the symp lectic integrator, which can simulate the long-term evolution of a Ham iltonian system accurately, The stability condition and the numerical dispersion of schemes with fourth-order accuracy in time and space usi ng the symplectic integrator are derived for the transverse electric ( TE)-mode in two dimensions. Their stable and accurate performance is q ualitatively verified, and is also demonstrated by numerical simulatio ns of wave-converging by a perfect electric conductor wall and propaga tion along a waveguide whose refractive index difference between the c ore and cladding is more than 9%.