Theory of isolated magnetic skyrmions: From fundamentals to room temperature applications
Spontaneous creation and annihilation dynamics and strain-limited stability of magnetic skyrmions
Frederic Rendell-Bhatti, Raymond J. Lamb, … Damien McGrouther
Deformation of the moving magnetic skyrmion lattice in MnSi under electric current flow
D. Okuyama, M. Bleuel, … T. J. Sato
Physical foundations and basic properties of magnetic skyrmions
Alexei N. Bogdanov & Christos Panagopoulos
Collective skyrmion motion under the influence of an additional interfacial spin-transfer torque
Callum R. MacKinnon, Katharina Zeissler, … Serban Lepadatu
Aggregation and collapse dynamics of skyrmions in a non-equilibrium state
Xiuzhen Yu, Daisuke Morikawa, … Yoshinori Tokura
Electric-field driven stability control of skyrmions in an ultrathin transition-metal film
Souvik Paul & Stefan Heinze
Multiple low-temperature skyrmionic states in a bulk chiral magnet
Lars J. Bannenberg, Heribert Wilhelm, … Andrey O. Leonov
A promising birthplace for skyrmions
Enrico Schierle
Wigner solids of domain wall skyrmions
Kaifeng Yang, Katsumi Nagase, … Hongwu Liu
Felix Büttner ORCID: orcid.org/0000-0002-6204-99481,
Ivan Lemesh ORCID: orcid.org/0000-0003-1130-80031 &
Geoffrey S. D. Beach ORCID: orcid.org/0000-0001-9158-74301
Scientific Reports volume 8, Article number: 4464 (2018) Cite this article
Magnetic properties and materials
Magnetic skyrmions are topological quasiparticles of great interest for data storage applications because of their small size, high stability, and ease of manipulation via electric current. However, although models exist for some limiting cases, there is no universal theory capable of accurately describing the structure and energetics of all skyrmions. The main barrier is the complexity of non-local stray field interactions, which are usually included through crude approximations. Here we present an accurate analytical framework to treat isolated skyrmions in any material, assuming only a circularly-symmetric 360° domain wall profile and a homogeneous magnetization profile in the out-of-plane direction. We establish the first rigorous criteria to distinguish stray field from DMI skyrmions, resolving a major dispute in the community. We discover new phases, such as bi-stability, a phenomenon unknown in magnetism so far. We predict materials for sub-10 nm zero field room temperature stable skyrmions suitable for applications. Finally, we derive analytical equations to describe current-driven dynamics, find a topological damping, and show how to engineer materials in which compact skyrmions can be driven at velocities >1000 m/s.
Magnetic skyrmions are spin configurations with spherical topology1,2,3,4,5, typically manifesting as circular domains with defect-free domain walls (DWs) in systems with otherwise uniform out-of-plane magnetization. Skyrmions are the smallest non-trivial structures in magnetism and they behave like particles6,7,8,9,10,11, which makes them of fundamental interest and of practical utility for high-density data storage applications12,13,14,15. Skyrmions have been investigated for decades16, but only recently has attention shifted to the detailed domain wall structure and thereby-determined topology. Two factors have driven this trend: technological advances allowing for direct imaging of the spin structure17,18,19 and the discovery that the Dzyaloshinskii-Moriya interaction (DMI) can be used to stabilize that structure. In particular, DMI can lead to a skyrmion global ground state above the Curie temperature (T c ) in a Ginzburg-Landau theory of a ferromagnet20,21. Stray fields are not included in that model, but are of critical importance for understanding skyrmions in real materials, as underlined by the fact that all isolated room temperature skyrmions so far were observed in relatively thick films with sizable saturation magnetization11,22,23,24,25,26,27,28 and hence very strong stray field interactions. To explore the full skyrmion phase diagram and understand skyrmion stability, a theory with accurate predictive power is required and such a theory must include stray field interactions.
Stray field energies are the most difficult to treat analytically due to their nonlocal nature: they involve six-dimensional integrals whose kernel locally diverges. All existing models involve approximations that are motivated only by the need to simplify, without considering the error of the approximation or the limits of validity16,20,29,30,31,32,33,34,35,36,37. Tu38, Kiselev et al.39, and Guslienko31 evaluated the full stray field integrals, but only numerically, which is prohibitively slow and technically demanding. Micromagnetic simulations are hence the only generally applicable tool to obtain quantitative predictions, but they are too slow to systematically examine skyrmion properties across a parameter space comprising four material parameters, film thickness, and magnetic field.
Here, we derive a fully analytical theory with the precision of micromagnetic simulations but orders of magnitude faster performance, providing unique access to the full skyrmion phase diagram and deep insights into the underlying physics. For instance, we can mathematically prove that there are two types of skyrmions, stray field skyrmions and DMI skyrmions, and we discuss how to experimentally distinguish between them. Our theory predicts a sharp transition separating these two skyrmion phases, and a phase pocket in which they can co-exist, leading to bi-stabilities and zero-stiffness skyrmions. These new states exist at room temperature and can have many novel applications, some of which we suggest here. In view of applications, we prove that the Co-based multilayers at the focus of most current experimental efforts are incapable of hosting room-temperature-stable sub-10nm skyrmions. However, by examining >106 material parameter combinations, we identify alternative materials suitable to host such skyrmions at room temperature without a stabilizing field. Finally, we derive simple analytical equations to describe current-driven skyrmion dynamics with accuracy comparable to micromagnetics but yielding key insights that cannot be gained from numerics. We discover a topological contribution to damping that severely reduces the mobility of small skyrmions, and propose materials that can mitigate this effect to permit sub-10 nm skyrmions to be driven by current at >1 km/s with vanishing skyrmion Hall angle.
Validation of the model
Our theory predicts the energy of isolated skyrmions in a film of arbitrary thickness and infinite in-plane extent, based on the recent experimental confirmation19,23 of an analytic and universal 360° domain wall (DW) model40 for the spin structure of all skyrmions. Although the 360° DW model is not an exact solution of the micromagnetic energy functional, it shows excellent agreement with experiments19,23 and is validated by our extensive micromagnetic simulations. Taking this model as an ansatz for the skyrmion spin structure, we derived analytic expressions for the total energy function with better than 1% accuracy over the entire parameter space (see Supplementary Information), assuming only that the structure does not vary along the out-of-plane direction. For a given set of material parameters (uniaxial anisotropy constant K u , saturation magnetization M s , exchange constant A, interface and bulk DMI strengths D i and D b , and magnetic layer thickness d) and external field H z , minimization of this function yields the equilibrium skyrmion configuration. Multilayers are included straightforwardly via the effective medium approach24,41 (for the limitations of the effective medium model, see ref.41). Our model is therefore generally applicable to skyrmions in any material, with the exception of thick multilayers with very strong stray fields and weak DMI where flux closure surface states form and where the assumption of uniform magnetization along the out-of-plane direction does not hold42. In the Supplementary Information, we provide analytical expressions for both bulk and interface DMI terms; in the discussions below, we apply our model with a focus on systems with interface DMI only, but the model itself is more general.
The spin structure is parameterized by its radius R, domain wall width Δ, domain wall angle ψ, and topological charge N (Fig. 1a). R and Δ determine the magnetization profile m z (x, y), whereas ψ specifies whether the in-plane component of the domain wall spins is radial (Néel, ψ = 0, π), azimuthal (Bloch, ψ = π/2, 3π/2), or intermediate. Although the original 360° domain wall model included an analytical function Δ(R)40, we find that stray field interactions lead to a highly-nontrivial dependence of Δ on R and they must be treated independently. For large ρ = R/Δ, skyrmions consist of an extended out-of-plane magnetized domain bounded by a narrow circular domain wall, while for \(\rho \sim 1\) the inner domain is reduced to a point-like core. We refer to these limiting cases as bubble skyrmions and compact skyrmions43, respectively, consistent with the literature, but note that many skyrmions observed recently11,22,24,44 showed intermediate values of ρ and cannot be classified distinctly.
Skyrmion profile and stabilizing energy. (a) Illustration of the 360° domain wall model. The plot shows the normalized perpendicular magnetization m z as a function of position x along the diameter of a skyrmion and defines the characteristic parameters R (radius), Δ (domain wall width), ψ (domain wall angle, inset), and N (polarity or skyrmion charge). The profile corresponds to the μ0H z = −1 T data point in (b). The negative value of the field indicates that it is oriented antiparallel to the skyrmion core. (b) Radius and domain wall width as a function of applied magnetic field μ0H z . The small solid data points are predictions of our analytical model, the solid large data points are results from micromagnetic simulations, the open points are experimental results of Romming et al.19. The insets show the relaxed spin structures obtained from micromagnetic simulations corresponding to the large solid data points. (c) Total energy, domain wall energy, and bulk energy as a function of skyrmion radius at a field of μ0H z = −2 T. At each value of R the energy has been minimized to determine Δ and ψ. Req is the equilibrium radius as plotted in (b). At the minimum, the domain wall energy has a negative slope, qualifying this skyrmion as a DMI skyrmion. (d) Decomposition of the total energy in c into individual components. Domain wall energies (yellow): DMI energy (inverted), exchange energy, effective anisotropy energy, and volume stray field energy (multiplied by 10). Bulk energies (pink): effective Zeeman energy and remaining surface stray field energies (multiplied by 10).
Figure 1b shows that our model agrees accurately with micromagnetic simulations and the experimental data of Romming et al.19. Our model yields the energy of a given skyrmion configuration in less than a millisecond on a regular personal computer and the full energy landscape in a few seconds, providing dramatic improvement over micromagnetic simulations. Moreover, it provides information not easily accessible by simulations. For example, since Δ and ψ can be minimized for any non-equilibrium R, one readily obtains the function E(R) that describes skyrmion stability and rigidity. By contrast, standard micromagnetic simulations only yield the equilibrium spin structure.
Figure 1c shows E(R) for the skyrmion in Fig. 1b at μ0H z = −2 T, giving insight into skyrmion stability. The individual energy contributions are plotted in Fig. 1d. There are two stable states: the skyrmion state at Req and the ferromagnetic ground state at R = 0. Despite their different topology, there is a path from one to the other through the singular R = 0 state. This singular state lacks a topology, hence topological quantization is lifted here. As seen in Fig. 1d, all energy terms vanish as R → 0 except for the exchange energy, leading to a finite, universal skyrmion energy at R = 0,
$$E(R=0)={E}_{0}\approx 27.3Ad.$$
This value is close to the absolute minimum 8πAd of the exchange energy of a continuous spin structure with integer topological charge3,36, and much smaller than the 38Ad predicted using a simple linear skyrmion profile33. Our identification of a topologically valid and energetically possible path to annihilation clarifies that all skyrmions can be annihilated, even in continuum models, despite the common perception that topological protection guarantees stability32,45,46. In particular, although E0 relates to topology as previously discussed by Belavin and Polyakov3, it corresponds to the skyrmion nucleation energy barrier E n . Skyrmion stability, however, is related to the annihilation energy barrier E a = E0 − E(R), which depends nontrivially on all micromagnetic parameters and has no direct topological origin.
The annihilation barrier in Fig. 1c corresponds to only \( \sim 3{k}_{B}{T}_{{\rm{RT}}}\) at TRT = 300 K, and hence the small skyrmions observed by Romming et al.19 can only exist at cryogenic temperatures. We focus in the remainder of this work on room-temperature stable skyrmions, i.e., those with E a > 50k B TRT corresponding to lifetimes >10 years, as required for applications. We note that E n and E a serve as upper bounds for the true energy barriers because skyrmions can deform in a way that is not covered by the 360° domain wall model, hence reducing the nominal energy barrier47. However, previous studies47 and our own micromagnetic simulations indicate that deformations reduce the energy barrier by <2Ad. Note that the collapse energy E0 in the atomistic simulations by Rohart et al.47 is 23Ad and 22Ad with and without deformations, respectively, significantly smaller than the continuum limit of 8πAd, which the authors attribute to lattice effects. Such effects are beyond the scope of the present continuum theory, and therefore, we limit our subsequent discussion to the validity range of the micromagnetic framework, i.e., to diameters of approximately 1 nm and larger, and subtract 2Ad from all energy barriers in an attempt to include possible deformations.
Smallest room-temperature stable size – unifying existing models
Figure 2 shows that our model reproduces the limiting cases treated previously in the literature and continuously connects them to provide otherwise inaccessible insights into the entire skyrmion phase diagram. Two common approximations exist to describe circular spin textures in out-of-plane magnetized films. The first is the wall energy model48 introduced in the 1970s to treat bubble domains. This model approximates the domain walls as infinitely thin with constant energy per unit length σDW, enabling exact solutions for the stray field energies48. A characteristic feature is that E(R) starts from zero with positive slope at R = 0 and goes through a maximum before reaching a minimum at R = Req. This minimum vanishes above a critical field and the collapse diameter is finite even at T = 0. The second model, derived by Bogdanov and Hubert34,35 for small skyrmions in ultrathin films, approximates stray field energies as effective anisotropies. The original model does not include analytic solutions, but we derive these equations in the Supplementary Information. The effective anisotropy model predicts isolated skyrmions for any finite DMI34, and they cannot be annihilated by applied fields at T = 0.
Fundamental predictions of our model. (a–c) Room temperature collapse diameter predicted by different models (always requiring 50k B T for stability, but all diagrams would look very similar for 30k B T). In all diagrams, d is the thickness of the magnetic material and the total film thickness (including non-magnetic spacer layers) is 4d. All three diagrams were derived with A = 10 pJ/m, Q = 1.4, M s = 1.4 MA/m, and H z such that R is minimum. (a) Wall energy model of bubble domains48. (b) Bogdanov and Hubert's model34,35, which considers stray fields by effective anisotropies. (c) Our model including stray field in their full non-local nature. (d) Characteristic energy plots for different values of D i at d = 31.6 nm. The yellow and pink solid lines are the domain wall and bulk energy corresponding to the red solid line, respectively. The bulk energy has been offset by 155Ad to present it in the given plot range. (e) Radial cross section of a simulated non-DMI skyrmion (green) and low DMI (stray field) skyrmion (light blue). The red and blue lines are reproduced from ref.39, where the attempt was made to distinguish proper skyrmions from stray-field stabilized bubbles based on their radial profile. The radial profile of the simulated stray field skyrmions precisely agrees with the suggested shape of a DMI skyrmion, demonstrating that radial profile is not a valid selection criterion. θ = arccos (m z ) is the out-of-plane spin angle. The inset shows a fit of the 360° domain wall model to the numerically obtained solutions of the Euler equation of a zero-M s system34. The following parameters were used for the simulations: d = 2.8 nm, A = 10 pJ/m, M s = 1.4 MA/m, Q = 1.01, D i = 0, μ0H z = −63 mT for the non-DMI skyrmion and d = 1 nm, A = 10 pJ/m, M s = 1.4 MA/m, Q = 1.01, D i = 0.8 mJ/m2, and μ0H z = −50 mT for the low DMI skyrmion (and no non-magnetic spacer layers in both cases).
The room temperature collapse diameter as a function of D i and film thickness is shown in Figs 2a,b for the wall-energy model and effective anisotropy model, respectively, using the most accurate available form of the former41. These figures correspond to a 50k B TRT stability criterion; a 30k B TRT criterion (lifetime of order seconds) leads to similar results. The wall energy model predicts room-temperature stable skyrmions for all D i and d, with a variable but size-independent ψ. These skyrmions are always relatively large. The effective anisotropy model, by contrast, predicts room-temperature stable skyrmions with size <10 nm, but only above a critical D i that is quite large in ultrathin films. With interfacial DMI, skyrmions in this model are always purely Néel (see Supplementary Information).
The predictions for our theory are shown in Fig. 2c. The most obvious feature is a sharp boundary separating \(\ll 10\,{\rm{nm}}\) and \(\gg 10\,{\rm{n}}{\rm{m}}\) skyrmions. This phase boundary qualitatively follows the stability boundary in Fig. 2b, and the predicted collapse diameter above and below this contour agrees surprisingly well with the effective anisotropy and the wall energy model, respectively. It would therefore be tempting to use Bogdanov and Hubert's model whenever it predicts room temperature stable skyrmions and the wall energy model in all other cases. This hybrid approximation would predict a collapse diameter diagram similar to the correct result in Fig. 2c, but there are important properties that this approach would omit. First, the transition from one model to the other is not always as sharp as the collapse-diameter figure suggests. At fields below the collapse field, large and small skyrmions can coexist, yielding the bi-stable region denoted in Fig. 2c. A hybrid model would be intrinsically incapable of describing this phase since it can only describe one type of skyrmion at a time. Second, the domain wall width and domain wall angle can strongly depend on the skyrmion radius, which is not predicted by either model.
Figure 2d shows how E(R) evolves with increasing DMI and how mixed states appear. Without DMI, E(R) starts with zero or positive slope at R = 0, then reaches a maximum and diverges logarithmically to −∞ if H z = 0. An applied field can compensate the divergent stray fields and stabilize a minimum at large R (not shown in the figure). With DMI, the situation changes qualitatively: the slope at R = 0 becomes negative so the minimum appears before the maximum. Initially, the minimum is too shallow and appears in condition of unstable Δ (see Supplementary Information section 6). This is the case for D i = 1.2 mJ/m2 in Fig. 2d. For larger D i , the minimum becomes deeper and shifts to larger R. At first, this new minimum is stable only at low fields (D i = 1.5 mJ/m2 in Fig. 2d) but collapses before the stray field minimum at large R vanishes. That is, the smallest skyrmions are not necessarily found at the largest possible field values. At even larger D i , the minimum at small R becomes more pronounced and can become energetically degenerate with a second minimum, as shown for D i = 2.5 mJ/m2 and D i = 3.0 mJ/m2. The important difference between these two states is that the energy barrier separating the minima is much larger than 50k B TRT for the lower DMI value (D i = 2.5 mJ/m2) and almost vanishingly small for the larger D i . Those states constitute a bi-stability and a zero stiffness skyrmion, respectively, as discussed in more detail below. At even larger DMI, the maximum between the two minima vanishes and only one stable state remains.
The spin textures in Fig. 2c are all topologically equivalent: they exhibit spherical topology and unity topological charge and hence are all magnetic skyrmions. However, there has been much disagreement in the literature as to whether different stabilization mechanisms lead to fundamentally different types of skyrmions, and, if so, how to distinguish between them. Kiselev et al.39,49 suggested to use the skyrmion profile to determine its type. Similar to Guslienko's study on stray field skyrmions in magnetic dots31, we find that the radial profile is not a robust distinguishing criterion, as shown in Fig. 2e. We find that stray fields can stabilize compact skyrmions in the absence of DMI, with sizes down to ~20 nm, despite frequent statements to the contrary22,39,50. On the other hand, stray field skyrmions with DMI below the transition line in Fig. 2c can have purely Néel domain walls, and possess either a compact or bubble-like profile (Fig. 2e). The profile can be continuously tuned from one to the other without changing the topology, either by applied field or by adjusting \(\kappa =\frac{\pi {D}_{i}}{4\sqrt{A{K}_{u}}}\) (Fig. 2e, inset), as also shown in34,35. Hence, skyrmions with DMI-stabilized Néel domain walls are not necessarily DMI-stabilized skyrmions, even if their profiles are compact, since the energy giving rise to the minimum in E(R) can still derive from stray fields.
Based on Fig. 2c and the underlying energetics, we offer here mathematically-precise and experimentally-accessible criterion to categorize skyrmions unambiguously, defining the terms "DMI skyrmions", "stray field skyrmions", and "bubble skyrmions" as labeled in Fig. 2. The mathematical definition is based on the understanding of Fig. 2d. We first note that all contributions to the total energy can be classified into domain wall energies (\(\propto \,R\) at large R) and bulk energies (all other energies), see also Figs 1c and d. Exchange, anisotropy, DMI, and volume stray field energies are domain wall energies. The Zeeman energy is a bulk energy. Surface stray fields contribute to both categories: the domain wall contribution leads to an effective reduction of the anisotropy and the bulk contribution effectively reduces the external field. We see that the minimum at small R is formed by the local energies while the minimum at large R is essentially a minimum in the bulk energies (which can be shifted due to the linear slope of the domain wall energies), see Fig. 2d. We therefore identify DMI skyrmions as minima of the domain wall energy. Stray fields may still play an important role in shaping this minimum and giving it the required depth for thermal stability, but the origin is still the minimum in the domain wall energy term that can only appear due to DMI. We include in this definition skyrmions whose domain wall energy has a negative slope at the minimum. By this definition, we include all skyrmions in materials with globally negative domain wall energy density. Physically, these states behave like minima in the domain wall energy, but they require a magnetic field to stabilize isolated skyrmions against expanding to helical states or skyrmion lattices. Finally, we call an energy minimum a "mixed state" if it is formed by minima of both domain wall and bulk energies (e.g., when both minima are at the same radius or when the energy barrier in between is so small (compared to k B T) that they form one extended stable state).
The domain wall energy minimum is the key aspect that the wall energy model neglects in its approximations. Its origin lies in "domain wall compression" wherein the reduction of Δ at small R leads to the zero radius energy of 27Ad (instead of 0 in the wall energy model). Domain wall compression leads to a steep increase of the domain wall energy (mostly exchange) at small R, which is why domain wall energy minima can be thermally stable down to much smaller radii than minima in bulk energies. This is illustrated in Fig. 1c. Even when considering that the ground state skyrmion energy already includes some compression energy, the difference between a linear interpolation of the domain wall energy according to the wall energy model and the correct result is 17Ad (27Ad instead of 10Ad). This additional 17Ad in stabilizing energy is one of the reasons why DMI skyrmions are physically different from stray field skyrmions.
Our mathematical definition is directly linked to discernibly different physical properties. First, the radius and in particular the collapse radii are different. Second, the collapse field shows distinct scaling: it is almost independent of DMI for stray field skyrmions and strongly increases with DMI for DMI skyrmions (see Fig. 3a). Third, Δ is determined by a competition between exchange and anisotropy for stray field skyrmions, but for DMI skyrmions it reflects a balance between DMI, anisotropy, and external field. That is, Δ is independent of D i for stray field skyrmions and scales linearly with D i for DMI skyrmions. Fourth, DMI skyrmions are rather insensitive to applied fields, while stray field skyrmions are quite sensitive. And fifth, ψ is variable for stray field skyrmions but not for DMI skyrmions.
Multiple minima in systems with competing DMI and stray-field energies. (a) Phase diagram of the multi-stability of skyrmions as a function of applied field and DMI strength, where white, blue, and orange indicate the regions of instability, mono-stability, and bi-stability, respectively. Stability requirement is that all energy barriers are larger than 50k B T. The red line is the collapse field of DMI skyrmions and the gray line is the collapse field of stray field skyrmions. The dashed horizontal line indicates the slice that is plotted in panel (b). The data was obtained for a magnetic layer thickness d = 31.6 nm, a total film thickness of 4d, and magnetic layer M s = 1.4 MA/m and quality factor Q = 1.01. (b) Radius as a function of applied field for skyrmions in the bi-stability region in (a) at D i = 1.4 mJ/m2. The background color indicates the level of multi-stability as defined (a). The left inset shows a simulation of the stray field and DMI skyrmion which are simultaneously present in the same simulation area at μ0H z = −254 mT. The right insets show zoom-ins on the individual skyrmions (stray field skyrmion on top, DMI skrymion at the bottom). The coloured arrows indicate the domain wall spin orientation. (c) Energy as a function of radius for a zero stiffness skyrmion. The thermal energy available to the skyrmion at room temperature is indicated by the blue shading and the dotted lines. All states between R = 2 nm and R = 11 nm are accessible by thermal excitation on a sub-nanosecond time scale. The shown data corresponds to d = 4 nm (no spacer layers), D i = 1.85 mJ/m2, Q = 1.4, and μ0H z = −100 mT. The minimum at \(R\approx 6\,{\rm{nm}}\) is an artifact within the 1% precision of our model. The model actually predicts at most two minima.
The collapse radius is a particularly useful and technologically-relevant distinguishing feature. As shown in Fig. 2c, DMI skyrmions are those with room-temperature collapse diameters <10 nm and stray field skyrmions are those collapsing at larger sizes. This observation is supported by similar diagrams as a function of all material parameters, see Supplementary Information. With a very few exceptions (at small thicknesses), the 10 nm threshold is unambiguous even with some experimental error. That is, skyrmions either collapse at much larger or much smaller diameters than 10 nm. Even though further research is required to check if this threshold would change for other classes of materials (e.g., B20 materials or ferrimagnets) and for skyrmion lattices, we conclude that the 10 nm threshold size is a practical and valid criterion to identify the type of a skyrmion experimentally, at least for transition metal ferromagnet multilayers.
Finally, we remark that all room temperature isolated skyrmions observed so far11,22,23,24,25,26,27,50,51,52,53,54,55,56 should be classified as stray field skyrmions. Isolated DMI skyrmions have been observed at cryogenic temperatures19,57,58, but not yet at room temperature. The DMI strength measured in sputtered multilayers, however, can be on the order of 2 mJ/m222,23,24,59,60,61,62,63, which we predict to be sufficient to stabilize DMI skyrmions. But the bi-stability region can make it difficult to obtain them by just shrinking stripe domains, in particular when pinning is sizable56. From our analysis is seems likely that DMI skyrmions can be observed in existing sputtered multilayer materials by using novel skyrmion generation approaches, such as spin-orbit torque nucleation51. To verify that these skyrmions are DMI skyrmions, one only needs to confirm that they can be as small as 10 nm in diameter.
Bi-stabilities and zero stiffness skyrmions
The phase diagram of skyrmions in materials with strong stray fields shows a richness that has not been explored due to the inherent limitations of existing models. A particularly fascinating example is the appearance of multiple (possibly degenerate) minima in E(R). These minima can be separated by energy barriers exceeding 50k B TRT at room temperature, indicating that both states can be simultaneously stable. Although the possible coexistence of compact and bubble-like has previously been suggested39,49, concrete model predictions of bi-stability have only been made in the case of systems stabilized by lateral constraints64. Our model shows that bistability is a more general phenomenon. As depicted in Fig. 3a, bi-stability exists in a small pocket of the phase diagram near the point where the collapse fields of stray field and DMI skyrmions coincide. The two types of skyrmions in the bi-stability region can have very different properties (Fig. 3b), confirmed by micromagnetic simulations: Their radii differ by more than one order of magnitude and their spin structure is Néel-like for the small (DMI) skyrmion and intermediate for the large (stray field) skyrmion. The different size and domain wall angle can be used to encode information or to move the skyrmions in non-collinear directions by spin orbit torques.
If the energy barrier between the multiple minima is on the order of k B T or smaller, thermal fluctuations can cause the system to dynamically oscillate between these minima with frequencies on the order of the attempt frequency (reported to be between 2.5 × 107s−1 65 and 4.5 × 10−9 s−1 47). We therefore call these states zero stiffness skyrmions. Figure 3c shows E(R) for a system in which the radius can thermally fluctuate between 2 nm and 11 nm, such that it exhibits effectively zero stiffness with respect to variations of radius within this range. Zero stiffness skyrmions exists at larger DMI than bi-stable skyrmions (see Fig. 2d). We expect that zero stiffness skyrmions have a very low resonance frequency associated with their breathing mode, which could be exploited in nonlinear skyrmion resonators66 and should have impact on their inertia11 and on skyrmion Hall angle52,53. The ultra short timescale of sizable thermal fluctuations (a factor five in radius) and the resulting randomness of spin-orbit torque driven motion could also be used to implement high-speed skyrmion randomizers and stochastic computing67.
Skyrmions for racetrack memory applications
We now consider the design of skyrmions suitable for applications, such as racetrack-type memory devices in which bit sequences are encoded by the presence and absence of skyrmions that can be shifted by current12,13,14,68. Three necessary attributes are (i) small bit sizes, (ii) long term thermal stability, and (iii) stability in zero or low applied field. Figure 4 explains why ferromagnetic films and multilayers are incapable of meeting these requirements, and identifies alternative materials that can host sub-10 nm zero-field skyrmions at room temperature.
Zero field skyrmions. (a) Diameter of zero field skyrmions as a function of anisotropy and saturation magnetization, in a 2 nm thick film with A = 10 pJ/m and D i = 2 mJ/m2. The black solid line is a plot of \(\kappa =\frac{1}{2}\) and the dashed line follows the solution of \(\sqrt{8A(Q-\mathrm{1)}{\mu }_{0}{M}_{s}^{2}}-\pi {D}_{i}={\mu }_{0}{M}_{s}^{2}d\), where the latter has been imposed as a threshold because stripe domains are expected to nucleate spontaneously if K u is smaller23. (b) Stabilizing energy (in units of k B T) for the skyrmions in (a). The boxes illustrate the range of anisotropy and saturation magnetization values reported in the literature for different material classes. (c,d) The same diagrams as in (a and b), but for an applied field of 100 mT. (e,f) Diameter of zero-M s skyrmions in a d = 5 nm (g) and a d = 10 nm (h) thick film as a function of exchange constant and DMI strength. The anisotropy constant at each point is adjusted to obtain \({E}_{a}^{{\rm{eff}}}=50{k}_{B}T\).
Zero field skyrmions reside at minima of the domain wall energy, with energy barriers \({E}_{a}^{0}\) and \({E}_{a}^{\infty }\) that prevent shrinking to zero size and infinite growth, respectively. We take the smaller of these as the effective energy barrier \({E}_{a}^{{\rm{eff}}}\) to estimate stability. Figures 4a,b show size and stability, respectively, for zero-field skyrmions as a function of K u and M s for thin films with D i = 2 mJ/m2, typical of heavy-metal/ferromagnet interfaces. Stray fields tend to reduce \({E}_{a}^{\infty }\) and hence zero field skyrmions occur preferably at low M s , but there is an inherent tradeoff in Figs 4a,b between stability and size. Moderate applied fields (Figs 4c,d) extend the stability range for large stray field skyrmions but do little to stabilize small skyrmions.
Since it is the domain wall energy whose minimum is responsible for zero-field skyrmions, stability can be enhanced by increasing film thickness since this energy scales with d. For higher M s materials, the benefit is offset by destabilizing stray fields that lead to stripe domain formation. But at low M s , sub-10 nm skyrmions can readily be achieved, as shown in Figs 4d,e for the limit M s = 0. These figures show zero-field skyrmion size versus A and D i , with K u adjusted at each point to ensure \({E}_{a}^{{\rm{eff}}}=50{k}_{B}{T}_{{\rm{RT}}}\). As a rule of thumb, κ ≈ 0.8 can be used to estimate K u from the diagram.
Figure 4b shows parameter ranges for several classes of materials, providing guidance to realize such skyrmions experimentally. The problem with thin-film ferromagnets is that the quality factor \(Q=\frac{2{K}_{u}}{{\mu }_{0}{M}_{s}^{2}}\) is usually constant when engineering multilayers or alloys, which correlates K u and M s in a way that makes the low-M s , high-K u region inaccessible. This is reflected in the diagonal lower boundary marking the accessible parameter range for such materials (red box, Fig. 4b), and in the trend for experimental film parameters (red circles) reported in the literature. This explains why room-temperature skyrmions observed in ferromagnetic films and multilayers have all been quite large. Note that all previously observed isolated zero field skyrmions in Co-based ferromagnets are only stabilized by pinning or sample boundaries23.
Ferrimagnets, by contrast, generally show little correlation between K u and M s , as seen in the representative literature values marked by triangles in Fig. 4b. In rare-earth-transition metal alloys, these parameters can be tuned independently by temperature or composition to achieve the combination of low-M s and high K u needed for stabilizing small, zero-field skyrmions. Moreover, the low M s and bulk perpendicular anisotropy in such materials allows for thicker films with strong PMA, which reduces the D i needed to stabilize small skyrmions (compare Figs 4e,f). Finally, as seen in Figs 4e,f, values for A between 4 pJ/m and 10 pJ/m are ideal, and in the rare-earth-transition-metal ferrimagnets the exchange constants are found to be just in this range. In addition to compensated ferrimagnets, natural antiferromagnets69,70,71,72,73 are quite promising for ultrasmall skyrmions, as are Co-based synthetic antiferromagnets (SAFs) that minimize stray fields, where heavy-metal interfaces can provide substantial K u and D i .
Spin-orbit torque driven motion of skyrmions
Finally, we derive analytical expressions for skyrmion dynamics that agree well with micromagnetic simulations and guide material selection to realize current-driven velocities exceeding 1000 m/s for sub-10 nm skyrmions. Figure 5a shows an intermediate skyrmion driven by damping-like spin-orbit torque (SOT) from a charge current density jHM in an adjacent heavy metal. Due to its topological charge, the skyrmion moves at an angle ξ′ with respect to jHM, a phenomenon known as the skyrmion Hall effect52,53. For Néel skyrmions ξ′ coincides with the skyrmion Hall angle ξ, but in general ξ ≠ ξ′ since the current-induced effective force F depends on domain wall angle ψ.
Spin-orbit torque driven skyrmion motion. (a–c) Simulated skyrmion motion for different ferromagnetic skyrmions, moved by a pure damping-like SOT current jHM = 1011 A/m2. (a) Intermediate skyrmion with ψ ≈ 45°, illustrating the different directions of current density j, force F, and velocity v as well as the difference between ξ and ξ′. Other parameters are M s = 1.4 MA/m, Q = 1.4, A = 10 pJ/m, D i = 0.7 mJ/m2, \({\mathscr{N}}=10\) layers (affects the vertical spin current), d = 10 nm and total film thickness (including non-magnetic spacers) of 40 nm. (b) Motion of a 2R ≈ 200 nm diameter skyrmion at μ0H z = −4.6 mT. (c) The same skyrmion as in (b), but at a field of μ0H z = −90 mT and a corresponding diameter of 2R ≈ 20 nm. (d) Mobility (velocity per current density) and skyrmion Hall angle as a function of skyrmion diameter (the same skyrmion as in (b) and (c). The size is controlled by the out-of-plane field). Material parameters are M s = 1MA/m, K u = 765kJ/m3, A = 20 pJ/m, D i = 2 mJ/m2, \({\mathscr{N}}=1\) layer, d = 1 nm without non-magnetic spacers. Simulations in (d) were performed using current densities between 1011A/m2 and 3 × 1011A/m2. (e) Mobility and skyrmion Hall angle for a SAF bilayer with fixed Ms1 of the bottom layer and variable Ms2 of the top layer. The same spin current is assumed for both layers. Other simulation parameters are Ax,y = 10 pJ/m, A z = −10 pJ/m, K u = 500 kJ/m3, D i = 2 mJ/m2, and d1 = d2 = 2.5 nm. No field is applied. (f) 3D representation of a simulated 10 nm diameter skyrmion in a SAF bilayer, moving colinear with the current at a velocity of 1000 m/s at jHM = 1012 A/m2, using Ms1 = Ms2 = 50 kA/m and other parameters as in (e). The two magnetic layers are colored black (top layer, predominantely magnetized down) and white (bottom layer, magnetized up). The skyrmion texture in the bottom layer can be seen in the reflection in mirror included in the ray-traced image. The small image shows a side view. Thermal stability of this skyrmion is 34k B T. In all simulations, α = 0.2 and \({\theta }_{{\rm{SH}}}^{{\rm{eff}}}=0.15\).
Figures 5b,c show micromagnetic simulations of SOT-driven motion for a 200 nm and 20 nm skyrmion, respectively, in a ferromagnetic Co-based multilayer. The skyrmion Hall angle and current-driven mobility μ = v/jHM are plotted versus skyrmion size in Fig. 5d. One sees that for small skyrmions ξ approaches π/2 and the mobility drops dramatically (linearly with radius) so that ultrasmall ferromagnetic skyrmions are effectively immobile. These phenomena have been noted in previous micromagnetic simulations50, but so far neither an intuitive explanation nor a materials-based solution for this technologically-critical challenge have been presented.
In the Supplementary Information we use the Thiele equation74 to derive the steady state velocity v = v(cosξ′, sinξ′) of a skyrmion driven by spin-transfer75 and spin-orbit torques based on the 360° DW model. Our treatment includes multilayers such as synthetic ferrimagnets or antiferromagnets with multiple sub-lattices i of variable d i and Ms,i, and can be extended to natural ferrimagnets by suitably scaling γ and α83. The solid lines in Fig. 5d,e show the results of our analytical model applied to ferromagnets and synthetic antiferromagnets, respectively. Although our model cannot account for internal deformations, we find excellent agreement with micromagnetic simulations for skyrmion diameters <200 nm, which distort minimally in micromagnetic simulations and are the most technologically relevant.
For damping-like SOT the angle between v and jHM is
$${\xi }^{\text{'}}=\xi -\psi +\pi \mathrm{[1}-{\rm{\Theta }}({\theta }_{{\rm{SH}}}^{{\rm{eff}}}N)]$$
$$\xi ={\rm{a}}{\rm{t}}{\rm{a}}{\rm{n}}2(\mathop{g}\limits^{ \sim },\alpha ),$$
$$\tilde{g}=-\frac{4\langle N\rangle }{{I}_{A}(\rho )},$$
$${I}_{A}(\rho )=2\rho +\frac{2}{\rho }+1.93(\rho -0.65)\exp [-\mathrm{1.48(}\rho -\mathrm{0.65)}]\mathrm{.}$$
Here α is the Gilbert damping, \({\theta }_{{\rm{SH}}}^{{\rm{eff}}}\) is the effective spin Hall angle, and I A is the exchange integral (see Supplementary Information). 〈x〉 denotes the average of x over all the layers or sublattices, weighted by dM s of each layer: 〈x〉 = (∑ i d i Ms,ix)/(∑ i d i Ms,i). The mobility is given by:
$$\mu =\sum _{i}\frac{1}{\sqrt{\mathop{g}\limits^{ \sim }{(\rho )}^{2}+{{\alpha }_{i}}^{2}}}\frac{{\rm{\Delta }}{I}_{D}(\rho )}{{I}_{A}(\rho )}\frac{|{\theta }_{{\rm{S}}{\rm{H}},i}^{{\rm{e}}{\rm{f}}{\rm{f}}}|}{{d}_{i}{M}_{s,i}}\frac{{\gamma }_{i}\hslash }{2|e|}.$$
Here, the final term represents the summation of the current-induced effective fields in each layer, with γ the gyromagnetic ratio, ħ the reduced Planck constant and e the electron charge. The first term is a damping term that includes a topological contribution proportional to 〈N〉. For compact skyrmions, \(|\tilde{g}(\rho )|\approx 0.9\gg \alpha \), i.e., the topological damping term dominates the dynamics. The ratio ΔI D (ρ)/I A (ρ) derives from the Thiele effective force, with
$${I}_{D}=\pi \rho +\frac{1}{2}\exp (-\rho )\mathrm{.}$$
For large ρ it is proportional to Δ, whereas for small ρ it is proportional to R, which explains why it is difficult to drive small skyrmions fast with SOT.
Based on these relations, there are only two ways to enhance the mobility of small, compact skyrmions: reduce the (topological) damping and enhance the current-induced effective field. The former can be accomplished using multi-sublattice materials and structures, such as ferrimagnets and SAF multilayers. The latter can be achieved using low-M s materials interfaced with strong spin Hall metals, since the current-induced effective field scales as \({\theta }_{{\rm{SH}}}^{{\rm{eff}}}/{M}_{s}\), as verified experimentally for Ta/TbCo76.
Figure 5e shows ξ and μ for the bilayer SAF structure in Fig. 5f as a function of the compensation ratio Ms1/Ms2 of the antiferromagnetically-coupled layers, with parameters such that the layers host skyrmions with diameter of 10 nm. One finds an enhancement of the mobility as the system approaches magnetic compensation due to a reduction in the topological damping. Note that the mobility takes a maximum at \({M}_{s2}/{M}_{s1} < 1\) because Ms1 is fixed in the simulations and Ms2 is varied. Thus as Ms2/Ms1 decreases from unity, the average M s also decreases, which initially leads to an increase of μ.
Prior micromagnetic simulations have suggested SAF structures to mitigate the skyrmion Hall effect by reducing the net topological charge. The corresponding reduction in topological damping is the main reason why larger skyrmion velocities are also observed70,77. We note that SOT driven domain walls in SAF structures have been shown to achieve enhanced mobility due to interlayer exchange torque78, a mechanism that is fundamentally different from what we identify here, and which is not effective for skyrmions due to the circular symmetry. Finally, Fig. 5e shows that by decreasing M s of the constituent layers to also increase the spin Hall effective field per unit current density, velocities exceeding 1000 m/s can be achieved at a reasonable current density of 1012A/m2 for 10 nm skyrmions, as required for applications.
Based on the material requirements identified above, (i) a low d i Ms,i in each layer, (ii) a large effective spin Hall angle in each layer, (iii) two (or more) antiferromagnetically coupled layers or sub-lattices with a low (ideally zero) average topological charge 〈N〉, we can now identify particular materials. We predict that one of the most promising materials is Pt/GdCo/GdFe/Ir/Ta. GdCo/GdFe bilayers can be coupled antiferromagnetically by just tuning one of the layers to be above and one to be below the compensation temperature79. The antiferromagnetic coupling is much stronger than for RKKY coupled multilayers and not sensitive to the layer thicknesses. The material has low M s and bulk perpendicular magnetic anisotropy, i.e., it can be easily grown to be 5 nm thick in a single layer. Pt and Ir lead to large additive DMI in combination with Co and Fe, respectively22,59. Finally, Pt and Ta are known to be efficient spin-orbit torque materials with opposite (additive) sign of the spin Hall angle.
In this article, we have presented the first unified theory that describes both stray field and DMI skyrmions in one coherent model. In order to make this model readily available, we have included software implementations in an online repository84.Thanks to the accurate expression of stray field interactions, we have been able to answer a number of fundamental physical questions, such as the meaning of topological stability and the proper and physically justified distinction of DMI and stray field skyrmions. Importantly, only DMI skyrmions can be <10 nm in diameter at room temperature. We have demonstrated that stray field interactions in ferromagnetic materials prevent the formation of such DMI skyrmions at at moderate applied fields and presently accessible values of D i ≤ 2 mJ/m2 and that low-M s materials are hence required to observe sub-10 nm skyrmions at room temperature. Furthermore, we find that small skyrmions in ferromagnets are much slower than in ferrimagnets or (synthetic) antiferromagnets. While this result is similar to regular domain walls, the underlying physics is different: while antiferromagnetic domain walls benefit from an exchange torque78, skyrmions instead are affected by a topological damping that can only be reduced in systems with antiferromagnetically coupled sublattices. Considering the topological damping is hence of crucial imporance for skyrmion-related materials science. It hence became clear that the future of skyrmionic lies in ferrimagnets and antiferromagnets and that transition metal ferromagnetic multilayers cannot be engineered to meet the requirements of data storage applications. Finally, we note that our model does not treat twisted states that could arise in low-DMI systems with strong stray field interactions42, nor does it treat skyrmions stabilized by more exotic phenomena such as frustrated exchange interactions80,81,82. However, the most suitable materials for applications are strong DMI systems with small, highly-mobile skyrmions, and our model accurately predicts all relevant physics for fundamental and applied studies of such skyrmions.
Skyrme, T. A unified field theory of mesons and baryons. Nuclear Physics 31, 556–569 (1962).
Article ADS MathSciNet CAS Google Scholar
Skyrme, T. H. R. A Non-Linear Theory of Strong Interactions. Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 247, 260–278 (1958).
Belavin, A. A. & Polyakov, A. M. Metastable states of two-dimensional isotropic ferromagnets. JETP Letters 22, 245 (1975).
ADS Google Scholar
Nagaosa, N. & Tokura, Y. Topological properties and dynamics of magnetic skyrmions. Nature Nanotechnology 8, 899–911 (2013).
Büttner, F. & Kläui, M. Magnetic Skyrmion Dynamics. In Skyrmions, Series in Material Science and Engineering 211–238 (CRC Press, 2016).
Tomasello, R. et al. A strategy for the design of skyrmion racetrack memories. Scientific Reports 4, 6784 (2014).
Makhfudz, I., Krüger, B. & Tchernyshyov, O. Inertia and Chiral Edge Modes of a Skyrmion Magnetic Bubble. Physical Review Letters 109, 217201 (2012).
Article ADS PubMed Google Scholar
Everschor, K., Garst, M., Duine, R. A. & Rosch, A. Current-induced rotational torques in the skyrmion lattice phase of chiral magnets. Physical Review B 84, 064401 (2011).
Iwasaki, J., Mochizuki, M. & Nagaosa, N. Current-induced skyrmion dynamics in constricted geometries. Nature Nanotechnology 8, 742–747 (2013).
Sampaio, J., Cros, V., Rohart, S., Thiaville, A. & Fert, A. Nucleation, stability and current-induced motion of isolated magnetic skyrmions in nanostructures. Nature Nanotechnology 8, 839–844 (2013).
Büttner, F. et al. Dynamics and inertia of skyrmionic spin structures. Nature Physics 11, 225–228 (2015).
Fert, A., Cros, V. & Sampaio, J. Skyrmions on the track. Nature Nanotechnology 8, 152–156 (2013).
Wiesendanger, R. Nanoscale magnetic skyrmions in metallic films and multilayers: a new twist for spintronics. Nature Reviews Materials 16044 (2016).
Rosch, A. Skyrmions: Moving with the current. Nature Nanotechnology 8, 160–161 (2013).
Kang, W., Huang, Y., Zhang, X., Zhou, Y. & Zhao, W. Skyrmion-Electronics: An Overview and Outlook. Proceedings of the IEEE 104, 2040–2061 (2016).
Malozemoff, A. P. & Slonczewski, J. C. Magnetic Domain Walls in Bubble Materials (Academic Press, New York, 1979).
Heinze, S. et al. Spontaneous atomic-scale magnetic skyrmion lattice in two dimensions. Nature Physics 7, 713–718 (2011).
Yu, X. Z. et al. Real-space observation of a two-dimensional skyrmion crystal. Nature 465, 901–904 (2010).
Romming, N., Kubetzka, A., Hanneken, C., von Bergmann, K. & Wiesendanger, R. Field-Dependent Size and Shape of Single Magnetic Skyrmions. Physical Review Letters 114, 177203 (2015).
Bogdanov, A. N. & Rößler, U. K. Chiral Symmetry Breaking in Magnetic Thin Films and Multilayers. Physical Review Letters 87, 037203 (2001).
Rößler, U. K., Bogdanov, A. N. & Pfleiderer, C. Spontaneous skyrmion ground states in magnetic metals. Nature 442, 797–801 (2006).
Moreau-Luchaire, C. et al. Additive interfacial chiral interaction in multilayers for stabilization of small individual skyrmions at room temperature. Nature Nanotechnology 11, 444–448 (2016).
Boulle, O. et al. Room-temperature chiral magnetic skyrmions in ultrathin magnetic nanostructures. Nature Nanotechnology 11, 449–454 (2016).
Woo, S. et al. Observation of room-temperature magnetic skyrmions and their current-driven dynamics in ultrathin metallic ferromagnets. Nature Materials 15, 501–506 (2016).
Yu, G. et al. Room-Temperature Creation and Spin-Orbit Torque Manipulation of Skyrmions in Thin Films with Engineered Asymmetry. Nano Letters 16, 1981 (2016).
Jiang, W. et al. Blowing magnetic skyrmion bubbles. Science 349, 283–286 (2015).
Woo, S. et al. Current-driven dynamics and inhibition of the skyrmion Hall effect of ferrimagnetic skyrmions in GdFeCo films. arXiv:1703.10310 (2017).
Tokunaga, Y. et al. A new class of chiral materials hosting magnetic skyrmions beyond room temperature. Nature Communications 6, 7638 (2015).
Kooy, C. & Enz, U. Experimental and Theoretical Study of the Domain Configuration in Thin Layers of BaFe12O19. Philips Res. Repts. 15, 7–29 (1960).
Ezawa, M. Giant Skyrmions Stabilized by Dipole-Dipole Interactions in Thin Ferromagnetic Films. Physical Review Letters 105, 197202 (2010).
Guslienko, K. Skyrmion State Stability in Magnetic Nanodots With Perpendicular Anisotropy. IEEE Magnetics Letters 6, 4000104 (2015).
Rohart, S. & Thiaville, A. Skyrmion confinement in ultrathin film nanostructures in the presence of Dzyaloshinskii-Moriya interaction. Physical Review B 88, 184422 (2013).
Bogdanov, A. N. & Yablonskii, D. A. Thermodynamically stable "vortices" in magnetically ordered crystals. The mixed state of magnets. Sov. Phys. JETP 68, 101 (1989).
Bogdanov, A. & Hubert, A. The Properties of Isolated Magnetic Vortices. physica status solidi (b) 186, 527–543 (1994).
Bogdanov, A. & Hubert, A. Thermodynamically stable magnetic vortex states in magnetic crystals. Journal of Magnetism and Magnetic Materials 138, 255–269 (1994).
Leonov, A. O. et al. The properties of isolated chiral skyrmions in thin magnetic films. New Journal of Physics 18, 065003 (2016).
Ivanov, B. A., Stephanovich, V. A. & Zhmudskii, A. A. Magnetic vortices – The microscopic analogs of magnetic bubbles. Journal of Magnetism and Magnetic Materials 88, 116–120 (1990).
Tu, Y.-O. Determination of Magnetization of Micromagnetic Wall in Bubble Domains by Direct Minimization. Journal of Applied Physics 42, 5704–5709 (1971).
Kiselev, N. S., Bogdanov, A. N., Schäfer, R. & Rößler, U. K. Chiral skyrmions in thin magnetic films: new objects for magnetic storage technologies? Journal of Physics D: Applied Physics 44, 392001 (2011).
Braun, H.-B. Fluctuations and instabilities of ferromagnetic domain-wall pairs in an external magnetic field. Physical Review B 50, 16485–16500 (1994).
Lemesh, I., Büttner, F. & Beach, G. S. D. Accurate model of the stripe domain phase of perpendicularly magnetized multilayers. Physical Review B 95, 174423 (2017).
Dovzhenko, Y. et al. Imaging the Spin Texture of a Skyrmion Under Ambient Conditions Using an Atomic-Sized Sensor. arXiv:1611.00673 (2016).
Fert, A., Reyren, N. & Cros, V. Magnetic skyrmions: advances in physics and potential applications. Nature Reviews Materials 2, 17031 (2017).
Yu, X., Tokunaga, Y., Taguchi, Y. & Tokura, Y. Variation of Topology in Magnetic Bubbles in a Colossal Magnetoresistive Manganite. Advanced Materials 29, 1603958 (2016).
Stier, M., Häusler, W., Posske, T., Gurski, G. & Thorwart, M. Skyrmion–Anti-Skyrmion Pair Creation by in-Plane Currents. Physical Review Letters 118 (2017).
Soumyanarayanan, A., Reyren, N., Fert, A. & Panagopoulos, C. Emergent phenomena induced by spin-orbit coupling at surfaces and interfaces. Nature 539, 509–517 (2016).
Rohart, S., Miltat, J. & Thiaville, A. Path to collapse for an isolated Néel skyrmion. Physical Review B 93, 214412 (2016).
Cape, J. A. & Lehman, G. W. Magnetic Domain Structures in Thin Uniaxial Plates with Perpendicular Easy Axis. Journal of Applied Physics 42, 5732–5756 (1971).
Kiselev, N. S., Bogdanov, A. N., Schäfer, R. & Rößler, U. K. Comment on "Giant Skyrmions Stabilized by Dipole-Dipole Interactions in Thin Ferromagnetic Films". Physical Review Letters 107, 179701 (2011).
Legrand, W. et al. Room-Temperature Current-Induced Generation and Motion of sub-100 nm Skyrmions. Nano Letters 17, 2703–2712 (2017).
Büttner, F. et al. Field-free deterministic ultrafast creation of magnetic skyrmions by spin–orbit torques. Nature Nanotechnology 12, 1040–1044 (2017).
Jiang, W. et al. Direct observation of the skyrmion Hall effect. Nature Physics 13, 162–169 (2016).
Litzius, K. et al. Skyrmion Hall effect revealed by direct time-resolved X-ray microscopy. Nature Physics 13, 170–175 (2017).
Hrabec, A. et al. Current-induced skyrmion generation and dynamics in symmetric bilayers. Nature Communications 8, 15765 (2017).
Montoya, S. A. et al. Tailoring magnetic energies to form dipole skyrmions and skyrmion lattices. Physical Review B 95, 024415 (2017).
Zeissler, K. et al. Pinning and hysteresis in the field dependent diameter evolution of skyrmions in Pt/Co/Ir superlattice stacks. arXiv:1706.01065 (2017).
Romming, N. et al. Writing and Deleting Single Magnetic Skyrmions. Science 341, 636–639 (2013).
Hanneken, C. et al. Electrical detection of magnetic skyrmions by tunnelling non-collinear magnetoresistance. Nature Nanotechnology 10, 1039–1042 (2015).
Soumyanarayanan, A. et al. Tunable room-temperature magnetic skyrmions in Ir/Fe/Co/Pt multilayers. Nature Materials 16, 898–904 (2017).
Emori, S. et al. Spin Hall torque magnetometry of Dzyaloshinskii domain walls. Physical Review B 90, 184427 (2014).
Hrabec, A. et al. Measuring and tailoring the Dzyaloshinskii-Moriya interaction in perpendicularly magnetized thin films. Physical Review B 90, 020402 (2014).
Jaiswal, S. et al. Investigation of the Dzyaloshinskii-Moriya interaction and room temperature skyrmions in W/CoFeB/MgO thin films and microwires. Applied Physics Letters 111, 022409 (2017).
Pollard, S. D. et al. Observation of stable Néel skyrmions in cobalt/palladium multilayers with Lorentz transmission electron microscopy. Nature Communications 8, 14761 (2017).
Article ADS PubMed PubMed Central Google Scholar
StartFragmentZelent, M. et al. Bi-stability of magnetic skyrmions in ultrathin multilayer nanodots induced by magnetostatic interaction. Physica Status Solidi Rapid Research Letters 11, 1700259 (2017).
Hagemeister, J., Romming, N., von Bergmann, K., Vedmedenko, E. Y. & Wiesendanger, R. Stability of single skyrmionic bits. Nature Communications 6, 8455 (2015).
Schwarze, T. et al. Universal helimagnon and skyrmion excitations in metallic, semiconducting and insulating chiral magnets. Nature Materials 14, 478–483 (2015).
Pinna, D. et al. Skyrmion Gas Manipulation for Probabilistic Computing. arXiv:1701.0775 (2017).
Parkin, S. S. P., Hayashi, M. & Thomas, L. Magnetic Domain-Wall Racetrack Memory. Science 320, 190–194 (2008).
Barker, J. & Tretiakov, O. A. Static and Dynamical Properties of Antiferromagnetic Skyrmions in the Presence of Applied Current and Temperature. Physical Review Letters 116, 147203 (2016).
Zhang, X., Zhou, Y. & Ezawa, M. Antiferromagnetic Skyrmion: Stability, Creation and Manipulation. Scientific Reports 6, 24795 (2016).
Hansen, P., Clausen, C., Much, G., Rosenkranz, M. & Witter, K. Magnetic and magneto-optical properties of rare-earth transition-metal alloys containing Gd, Tb, Fe, Co. Journal of Applied Physics 66, 756–767 (1989).
Hansen, P., Klahn, S., Clausen, C., Much, G. & Witter, K. Magnetic and magneto-optical properties of rare-earth transition-metal alloys containing Dy, Ho, Fe, Co. Journal of Applied Physics 69, 3194–3207 (1991).
Wu, F. et al. Epitaxial Mn2.5Ga thin films with giant perpendicular magnetic anisotropy for spintronic devices. Applied Physics Letters 94, 122503 (2009).
Thiele, A. A. Steady-State Motion of Magnetic Domains. Physical Review Letters 30, 230–233 (1973).
Li, Z. & Zhang, S. Domain-Wall Dynamics and Spin-Wave Excitations with Spin-Transfer Torques. Physical Review Letters 92, 207203 (2004).
Ueda, K., Mann, M., de Brouwer, P. W. P., Bono, D. & Beach, G. S. D. Temperature dependence of spin-orbit torques across the magnetic compensation point in a ferrimagnetic TbCo alloy film. Physical Review B 96, 064410 (2017).
Zhang, X., Zhou, Y. & Ezawa, M. Magnetic bilayer-skyrmions without skyrmion Hall effect. Nature Communications 7, 10293 (2016).
Yang, S.-H., Ryu, K.-S. & Parkin, S. Domain-wall velocities of up to 750 m s–1 driven by exchange-coupling torque in synthetic antiferromagnets. Nature Nanotechnology 10, 221–226 (2015).
Mangin, S. et al. Engineered materials for all-optical helicity-dependent magnetic switching. Nature Materials 13, 286–292 (2014).
Leonov, A. O. & Mostovoy, M. Multiply periodic states and isolated skyrmions in an anisotropic frustrated magnet. Nature Communications 6, 8275 (2015).
Lin, S.-Z. & Hayami, S. Ginzburg-Landau theory for skyrmions in inversion-symmetric magnets with competing interactions. Physical Review B 93, 064430 (2016).
Zhang, X. et al. Skyrmion dynamics in a frustrated ferromagnetic film and current-induced helicity locking-unlocking transition. Nature Communications 8, 1717 (2017).
Stanciu, C.D. et al. Ultrafast spin dynamics across compensation points in ferrimagnetic: The role of angular momentum compensation. Physical Review B 73(22) (2006).
https://github.com/felixbuettner/SkyrmionEnergy. (access date: 22 Feb. 2018).
This work was supported by the U.S. Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES) under Award #DE-SC0012371. FB thanks Alexander Stottmeister, Benjamin Krüger, and Kai Litzius for fruitful discussions and the German Science Foundation for financial support under grant number BU 3297/1-1.
Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, 02139, USA
Felix Büttner, Ivan Lemesh & Geoffrey S. D. Beach
Felix Büttner
Ivan Lemesh
Geoffrey S. D. Beach
F.B. conceived the project, derived the energy functions, generated and analyzed the data, and drafted the manuscript. I.L. derived the analytic expression for ψ. GSDB supervised the project. All authors discussed the results, the implications, and the figures and commented on the manuscript.
Correspondence to Felix Büttner.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Büttner, F., Lemesh, I. & Beach, G.S.D. Theory of isolated magnetic skyrmions: From fundamentals to room temperature applications. Sci Rep 8, 4464 (2018). https://doi.org/10.1038/s41598-018-22242-8
Sign reversal and manipulation of anomalous Hall resistivity in facing-target sputtered Pt/Mn4N bilayers
Ze-Yu Zhang
Jia-Wei Jiang
Wen-Bo Mi
Rare Metals (2023)
Ferrimagnetic spintronics
Se Kwon Kim
Hyunsoo Yang
Nature Materials (2022)
Simulating anti-skyrmions on a lattice
Juan C. Criado
Sebastian Schenk
L. A. Turnbull
Current-driven dynamics and ratchet effect of skyrmion bubbles in a ferrimagnetic insulator
Saül Vélez
Sandra Ruiz-Gómez
Pietro Gambardella
Nature Nanotechnology (2022)
History-dependent domain and skyrmion formation in 2D van der Waals magnet Fe3GeTe2
M. T. Birch
L. Powalla
G. Schütz