The skyrmion racetrack is a promising concept for future information technology. There, binary bits are carried by nanoscale spin swirls–skyrmions–driven along magnetic strips. Stability of the skyrmions is a critical issue for realising this technology. Here we demonstrate that the racetrack skyrmion lifetime can be calculated from first principles as a function of temperature, magnetic field and track width. Our method combines harmonic transition state theory extended to include Goldstone modes, with an atomistic spin Hamiltonian parametrized from density functional theory calculations. We demonstrate that two annihilation mechanisms contribute to the skyrmion stability: At low external magnetic field, escape through the track boundary prevails, but a crossover field exists, above which the collapse in the interior becomes dominant. Considering a Pd/Fe bilayer on an Ir(111) substrate as a well-established model system, the calculated skyrmion lifetime is found to be consistent with reported experimental measurements. Our simulations also show that the Arrhenius pre-exponential factor of escape depends only weakly on the external magnetic field, whereas the pre-exponential factor for collapse is strongly field dependent. Our results open the door for predictive simulations, free from empirical parameters, to aid the design of skyrmion-based information technology.

The mechanism and activation energy for the annihilation of a magnetic skyrmion is studied by finding the minimum energy path for the transition in a system described by a Heisenberg-type Hamiltonian extended to include dipole-dipole, Dzyaloshinskii-Moriya, and anisotropy interactions so as to represent a Co monolayer on a Pt(111) surface. The annihilation mechanism involves isotropic shrinking of the skyrmion and slow increase of the energy until the transition state is reached after which the energy drops abruptly as the ferromagnetic final state forms. The maximum energy along the minimum energy path, which gives an estimate of the activation energy within the harmonic approximation of transition state theory, is found to be in excellent agreement with direct Langevin dynamics simulations at relatively high temperature carried out by Rohart et al. [ Phys. Rev. B 93, 214412 (2016)]. The dipole-dipole interaction, the computationally most demanding term in the Hamiltonian, is found to be important but its effect on the stability of the skyrmion and shape of the transition path can be mimicked accurately by reducing the anisotropy constant in the Hamiltonian.

A rate theory for thermally activated transitions in spin systems is presented. It is based on a transition-state approximation derived from Landau-Lifshitz equations of motion and quadratic expansion of the energy surface at minima and first order saddle points. While the flux out of the initial state vanishes at first order saddle points, the integrated flux over the hyperplanar transition state is nonzero and gives a rate estimate in good agreement with direct dynamical simulations of test systems over a range in damping constant. The preexponential factor obtained for transitions in model systems representing nanoclusters with 3 to 139 transition metal adatoms is on the order of 1011 to 1013s$-${}1, similar to that of atomic rearrangements.

Various transitions that a magnetic Skyrmion can undergo are found in calculations using a method for climbing up the energy surface and converging onto first order saddle points. In addition to collapse and escape through a boundary, the method identifies a transition where the Skyrmion divides and forms two Skyrmions. The activation energy for this duplication process can be similar to that of collapse and escape. A tilting of the external magnetic field for a certain time interval is found to induce the duplication process in a dynamical simulation. Such a process could turn out to be an important avenue for the creation of Skyrmions in future magnetic devices.

Theoretical calculations of thermal spin transitions in nanoscale clusters on a surface are presented. The mechanisms for magnetization reversal are identified and the activation energy and pre-exponential factor for the rate are evaluated using a recently developed harmonic transition state theory and a Heisenberg-type Hamiltonian. A maximum is found in the pre-exponential factor as a function of cluster size corresponding to a crossover from a uniform rotation mechanism to temporary domain wall formation. As the islands grow, the energy barrier increases up to a limit where the domain wall is fully established. For larger islands, the minimum energy path becomes flat resulting in a significant recrossing correction to the transition state theory estimate of the rate. The parameters of the Hamiltonian are chosen to mimic Fe clusters on a W(110) surface, a system that has previously been studied extensively in the laboratory and the calculated results are found to be in close agreement with the reported measurements.