Localized spin waves in isolated $k\pi$ skyrmions

The localized magnon modes of isolated $k\pi$ skyrmions on a field-polarized background are analyzed based on the Landau-Lifshitz-Gilbert equation within the terms of an atomistic classical spin model, with system parameters based on the Pd/Fe biatomic layer on Ir(111). For increasing skyrmion order $k$ a higher number of excitation modes are found, including modes with nodes in the radial eigenfunctions. It is shown that at low fields $2\pi$ and $3\pi$ skyrmions are destroyed via a burst instability connected to a breathing mode, while $1\pi$ skyrmions undergo an elliptic instability. At high fields all $k\pi$ skyrmions collapse due to the instability of a breathing mode. The effective damping parameters of the spin waves are calculated in the low Gilbert damping limit, and they are found to diverge in the case of the lowest-lying modes at the burst and collapse instabilities, but not at the elliptic instability. It is shown that the breathing modes of $k\pi$ skyrmions may become overdamped at higher Gilbert damping values.