Based on the Padé approximation and multistep method, we propose an implicit high-order unconditionally stable complex envelope algorithm to solve the time-dependent Maxwell's equations. Unconditional numerical stability can be achieved simultaneously with a high-order accuracy in time. As we adopt the complex envelope Maxwell's equations, numerical dispersion and dissipation are very small even at comparatively large time steps. To verify the capability of our algorithm, we compare the results of the proposed method with the exact solutions.
ASJC Scopus subject areas
- Atomic and Molecular Physics, and Optics