Magnetic Skyrmions belong to the most interesting spin structures for the development of future information technology as they have been predicted to be topologically protected. To quantify their stability, we use an innovative multiscale approach to simulating spin dynamics based on the Landau-Lifshitz-Gilbert equation. The multiscale approach overcomes the micromagnetic limitations that have hindered realistic studies using conventional techniques. We first demonstrate how the stability of a Skyrmion is influenced by the refinement of the computational mesh and reveal that conventionally employed traditional micromagnetic simulations are inadequate for this task. Furthermore, we determine the stability quantitatively using our multiscale approach. As a key operation for devices, the process of annihilating a Skyrmion by exciting it with a spin polarized current pulse is analyzed, showing that Skyrmions can be reliably deleted by designing the pulse shape.Magnetic Skyrmions [1] are topological spin structures that arise in the spin pattern of ferromagnetic systems with broken inversion symmetry, such as chiral crystals [2,3] or thin magnetic films with different top and bottom interfaces [4,5]. Skyrmion lattices [6][7][8][9] constitute the ground state for some systems, while isolated Skyrmions can appear as metastable states of some magnetic nanostructures [10]. Isolated Skyrmions have been recently considered [11][12][13][14] as the building blocks for ultradense magnetic storage devices [15].Skyrmions carry a topological charge Q = ±1 defined as [16]:where A is the area of the system and m the unit magnetization vector. Since transitions that change Q are forbidden [16] in a continuum description of m, such structures are topologically protected. Nevertheless, in a real system composed of discrete magnetic moments localized on the atomic lattice sites, no strict topological protection exists [17]. Thus it is necessary to overcome a finite energy barrier to induce transformations that change Q, such as the annihilation of a Bloch line (BL) [16,[18][19][20][21][22][23][24]. The stability against external fields is indeed a key feature of Skyrmions, making them a good candidate as information carriers in next generation storage devices [25][26][27]. The fundamental prerequisites for applications are ascertaining the stability of Skyrmions, as well as reliably annihilating them. However, the computational treatment of processes involving annihilating Skyrmions is very delicate. In analytical micromagnetic theory, singularities in the exchange field tend to arise during topological transformations, making numerical simulations very susceptible to the mesh being used [28] and therefore often inaccurate. The necessity for a computational model, capable of performing quantitatively accurate simulations is therefore obvious and a key step. While more accurate atomistic simulations would overcome this problem, the computational power required to run such simulations for a sample of realistic experimental size makes th...