Finite-difference time-domain simulation of strong-field ionization: Perfectly matched layer approach
A Finite-Difference Time-Domain (FDTD) scheme with Perfectly Matched Layers (PMLs) is considered for solving the time-dependent Schrödinger equation, and simulate the ionization of an electron initially bound to a one-dimensional $δ$-potential, when applying a strong time-oscillating electric field. The performance of PMLs based on different absorption functions are compared, where we find slowly growing functions to be preferable. PMLs are shown to be able to reduce the computational domain, and thus the required numerical resources, by several orders of magnitude. This is demonstrated by testing the proposed method against an FDTD approach without PMLs and a very large computational domain. We further show that PMLs outperform the well known Exterior Complex Scaling (ECS) technique for short-range potentials when implemented in FDTD, though ECS remains superior for long-range potentials. The accuracy of the method is furthermore demonstrated by comparing with known numerical and analytical results for the $δ$-potential.