A PML utilizing k-vector information as applied to the whistler mode in a magnetized plasma

The perfectly matched layer (PML) in its most current form has shortcomings in certain cases involving propagation in anisotropic media, such as the whistler mode in a magnetized plasma, for which the vector component normal to the PML of the group velocity vector and the k-vector are anti-parallel with each other. We present a new type of PML that utilizes information on the k-vector direction by applying relevant spatial derivatives to the PML update equations. Derived expressions for the numerical reflection coefficient are used to quantify the performance of the PML for any incident angle. Numerical reflection calculations as well as time domain calculations are performed for whistler mode propagation in a magnetized plasma. We demonstrate the stability of the new PML and find that for the reflection coefficient Gamma, values of up to -40 dB can be realized