
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39227610
52072
10.1038/s41467-024-52072-4
Article
Confined antiskyrmion motion driven by electric current excitations
http://orcid.org/0000-0003-3238-8479
Guang Yao yao.guang@riken.jp

1
http://orcid.org/0000-0001-9656-9696
Zhang Xichao 2
Liu Yizhou 1
http://orcid.org/0000-0002-8580-2285
Peng Licong 1
Yasin Fehmi Sami 1
http://orcid.org/0000-0002-0382-895X
Karube Kosuke 1
http://orcid.org/0000-0003-3799-6801
Nakamura Daisuke 1
http://orcid.org/0000-0001-7924-6000
Nagaosa Naoto 13
http://orcid.org/0000-0003-1896-699X
Taguchi Yasujiro 1
http://orcid.org/0000-0002-0128-0580
Mochizuki Masahito 2
http://orcid.org/0000-0002-2732-4983
Tokura Yoshinori 145
http://orcid.org/0000-0003-3136-7289
Yu Xiuzhen yu_x@riken.jp

1
1 https://ror.org/03gv2xk61 grid.474689.0 RIKEN Center for Emergent Matter Science (CEMS), Wako, Japan
2 https://ror.org/00ntfnx83 grid.5290.e 0000 0004 1936 9975 Department of Applied Physics, Waseda University, Tokyo, Japan
3 https://ror.org/01sjwvz98 grid.7597.c 0000 0000 9446 5255 Fundamental Quantum Science Program, TRIP Headquarters, RIKEN, Wako Japan
4 https://ror.org/057zh3y96 grid.26999.3d 0000 0001 2169 1048 Department of Applied Physics, The University of Tokyo, Tokyo, Japan
5 https://ror.org/057zh3y96 grid.26999.3d 0000 0001 2169 1048 Tokyo College, The University of Tokyo, Tokyo, Japan
4 9 2024
4 9 2024
2024
15 77014 7 2024
23 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Current-driven dynamics of topological spin textures, such as skyrmions and antiskyrmions, have garnered considerable attention in condensed matter physics and spintronics. As compared with skyrmions, the current-driven dynamics of their antiparticles – antiskyrmions − remain less explored due to the increased complexity of antiskyrmions. Here, we design and employ fabricated microdevices of a prototypical antiskyrmion host, (Fe0.63Ni0.3Pd0.07)3P, to allow in situ current application with Lorentz transmission electron microscopy observations. The experimental results and related micromagnetic simulations demonstrate current-driven antiskyrmion dynamics confined within stripe domains. Under nanosecond-long current pulses, antiskyrmions exhibit directional motion along the stripe regardless of the current direction, while the antiskyrmion velocity is linearly proportional to the current density. Significantly, the antiskyrmion mobility could be enhanced when the current flow is perpendicular to the stripe direction. Our findings provide novel and reliable insights on dynamical antiskyrmions and their potential implications on spintronics.

Antiskyrmions, like skyrmions, are a form of topological spin texture, with a topological charge of opposite sign to the equivalent skyrmion with the same polarity. While antiskyrmions have been less explored, they offer some potential advantages for applications, and here, Guang et al demonstrate antiskyrmion motion within stripe domains.

Subject terms

Magnetic devices
Topological defects
https://doi.org/10.13039/501100001691 MEXT | Japan Society for the Promotion of Science (JSPS) 19H00660 20H00337 23H05431 23H04522 23H01841 24H02231 Guang Yao https://doi.org/10.13039/501100003382 MEXT | JST | Core Research for Evolutional Science and Technology (CREST) JPMJCR1874 JPMJCR20T1 Guang Yao https://doi.org/10.13039/501100004423 Waseda University 2024C-153 Guang Yao RIKEN TRIP Advanced General Intelligence for Science Program (RIKEN Quantum)issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Magnetic skyrmions1–3 and antiskyrmions2,4 are topologically nontrivial spin textures that have garnered significant attention due to their particle-like properties and unique dynamic behaviors. These topological spin textures hold great promise as information carriers for applications of storage and memory for the capability of nanometric size, non-volatile, and low threshold density for current drive5. One important intrinsic origin of these spin textures is the Dzyaloshinskii-Moriya interaction6 (DMI), an antisymmetric spin exchange interaction arising from spin-orbit coupling in inversion-symmetry broken system. Since the small-angle neutron scattering measurement on a hexagonal skyrmion lattice in bulk MnSi1 and the following Lorentz transmission electron microscopy (LTEM) observation on magnetic topology of skyrmions in Fe0.5Co0.5Si3, extensive studies on skyrmions have been conducted, in particular on their dynamical features, including thermally induced skyrmion diffusions7–9, field-controlled transformations of skyrmion states10, and current-driven motion of skyrmions11–16.

While a wide range of skyrmion-hosting materials have been discovered17, rare antiskyrmion-hosting materials have been explored4,18–23, which exhibit tetragonal structure with D2d (Mn1.4Pt0.9Pd0.1Sn4,18, Mn1.4PtSn19) or S4 [(Fe0.63Ni0.3Pd0.07)3P20–23] crystal symmetry. In stark contrast to the skyrmion-hosting materials with a unique DMI vector, the aforementioned antiskyrmion materials have anisotropic DMI vectors, resulting in fertile spin textures in a, such as elliptical skyrmions, square-shaped antiskyrmions constituted by alternative Bloch lines and Bloch domain walls, and non-topological bubbles in external magnetic fields. Due to the competition between the anisotropic DMI and dipole-dipole interactions, antiskyrmions in these materials also possess rich three-dimensional structures18,22. In addition to square-shaped antiskyrmions, a recent report claimed that the arc-shaped antiskyrmion appears in a prototypical skyrmion-hosting helimagnet FeGe with a single DMI vector24.

On the other hand, the spintronics field shows an increasing interest in topological spin textures emerged in the helical phase background. For example, recent studies have revealed that interstitial skyrmions within helical stripes exhibit enhanced mobility compared to skyrmions in a ferromagnetic background in FeGe16. Moreover, H-shaped skyrmions (or half-skyrmion pairs) in a helical background may demonstrate better stability and mobility in comparison with interstitial skyrmions25,26.

Skyrmion and antiskyrmion carry opposite topological charges, Q=−1 and +1, as given by Q=14π∫m⋅∂xm×∂ymdxdy, where m represents the unit vector of the local spin. Note here that the field reversed can rise to an opposite sign topological charge in topological spin textures. Therefore, it is expected that their emergent properties, which are determined by the topological charge, may also differ from each other. For example, in analogy to the skyrmion Hall motion12,27, theoretical simulations have predicted a current-driven antiskyrmion Hall motion with a characteristic Hall angle relative to the applied current direction and induced by the effective Magnus force28. Note that the sign of both the topological Hall effect of conduction electrons and (anti)skyrmion Hall effect are determined by the topological charge Q, and they should be opposite for skyrmions (Q=−1) and antiskyrmions (Q=+1). Indeed, a recently reported sign reversal in the topological Hall effect in a Heusler thin film magnet was attributed to the transformation between skyrmions and antiskyrmions29. Current-driven skyrmion dynamics have been well characterized, including studies of skyrmion motion in confined geometries30–32 and with various current orientations relative to the sample geometry33. However, the current-driven motion of an antiskyrmion embedded in stripe domains of a helical phase background, which contains four Bloch lines, has rarely been experimentally reported. The main reason is that, compared to skyrmions, creating and stabilizing antiskyrmions in experiments is challenging as it usually requires a precise control of intrinsic magnetic parameters, such as anisotropic DMI and anisotropy, which is technically demanding. Also, existing techniques and tools for studying skyrmions may not be directly applicable to the observation and manipulation of antiskyrmions due to their different magnetic properties.

Despite these challenges, there is a growing interest in antiskyrmions and their dynamics. Comprehensive studies on antiskyrmions are not only critical for the expansion of the fundamental understanding of topological phenomena, but also important for the design and development of future spintronic devices based on different topological spin textures.

Here, we have demonstrated the dynamical behaviors of antiskyrmions confined into stripes under nanosecond current-pulse excitations in FNPP using in situ LTEM. The confined antiskyrmions, which are efficiently generated in stripe-shape channels by tuning the magnitude and direction of external magnetic fields, remain even at zero field as a metastable state. The in situ LTEM results show that (1) the antiskyrmions move along the stripes regardless of whether the current flows parallel or perpendicular to the stripes, indicating that the Hall motion of such confined antiskyrmions is completely suppressed when the current is parallel to the stripes; (2) the antiskyrmion velocity is linearly proportional to the current density; (3) an enhancement in the antiskyrmion velocity has been revealed when the current flows perpendicular to the stripes.

Micromagnetic simulations for current-driven antiskyrmions in FNPP nanostripes

We focus on current-driven antiskyrmion dynamics in the tetragonal FNPP (space group I4¯). The crystal structure of FNPP along the [001] crystal axis is drawn schematically in Fig. 1a. The square-shaped antiskyrmions (Fig. 1b) are stabilized in FNPP by applying an external magnetic field22,34 of several hundred milliteslas (mTs)20. Antiskyrmions thus created can remain robust even after reducing the field to zero at room temperature23,34. In this study, we target metastable antiskyrmions embedded in a helical background and confined within the stripes. These stripes make it possible to investigate the antiskyrmion dynamics in a confined geometry with current flow parallel or perpendicular to the stripes. To simplify the text description, we denoted the stripes aligned along the [11¯0] and [110] crystal axes as horizontal (H-Str) and vertical (V-Str), respectively, as shown in Fig. 1c, d.Fig. 1 Simulated current-driven antiskyrmion motion in helical stripes.

a Crystal structure of M3P ((Fe0.63Ni0.3Pd0.07)3P) with I4¯ (S4 symmetry) space group along the [001] axis. The M1 (light purple), M2 (orange), and M3 (gray) denote three inequivalent M atomic sites for Fe, Ni, or Pd. Schematic of the real-space structure of an antiskyrmion (b), horizontal helical stripe (c, H-Str), and vertical helical stripe (d, V-Str). The colors and arrows in (b–f) indicate the vector fields coded by the color wheel shown below (a). The white and black colors indicate up (+z) and down (−z) magnetizations perpendicular to the plane. Snapshots of the real-space trajectory of an antiskyrmion embedded in a stripe by applying current pulses (e) perpendicular (j⊥V-Str) and (f) parallel (j∥H-Str) to the helical background stripes. The scale bar is 0.2 μm. g The antiskyrmion velocity |v¯x| and |v¯y| versus j for j⊥V-Str (red squares) and j∥H-Str (blue circles).

Current-driven motion of antiskyrmion via the spin transfer torque (STT) effect can be described by the Thiele equation30–33:1 G×vs−vd+Dβvs−αvd+F=0,

which is comprised of three distinct terms. The first term describes the gyrotropic effect, where G = (0, 0, 4πQ). Here, vd is the drift velocity of the spin texture, and the conduction electron velocity vs can be defined in terms of the current density j as vs = −(gPμB/2eMs) j, where ɡ is the Landé g-factor, μB is the Bohr magneton, e is the elementary charge, Ms is the saturation magnetization, and P is the spin polarization ratio of the conduction electrons. The second term encompasses the dissipative forces and includes the Gilbert damping coefficient α, non-adiabatic STT parameter β, and the dissipative tensor D, which is determined by the shape of the topological spin texture. The third term F accounts for the forces exerted on the spin textures related to external fields, boundary confinement, and impurities.

Skyrmions confined within micrometer-scale geometries have been observed to exhibit distinctive velocity-current density relationships depending on the direction of the current with respect to the device boundary33. Specifically, when the current flows along the boundary, the velocity scales as vd ∝ (β/α) j, and when current flows perpendicular to the boundary, it scales as vd ∝ −(4πQ/α) j33. Note that the skyrmion dynamics in helical background has been studied, and it was found an H-shaped skyrmion and anti-H-shaped skyrmion move in opposite directions along the current flow direction26. However, in the case of current-driven antiskyrmion dynamics confined by a stripe domain, the relationship remains uncharted and has not yet been fully understood. We performed the simulations first for the current-driven dynamics of antiskyrmions at the two states which are embedded in H-Str and in V-Str, respectively. For both states, current pulses with the duration of 1 ns and the magnitude |j| = 1.0 × 1011A m−2 flow along the horizontal axis (the x-axis indicated in Fig. 1f). The values of α and β for the simulation are set to 0.05. Representative snapshots of the simulations displayed in Fig. 1e, f depict the current-driven motions of an antiskyrmion embedded in (e) V-Str and (f) H-Str under an external field μ0H=50mT (along +z^), revealing that the antiskyrmion is confined to one-dimensional motion by the stripe domains independent of the current flow direction. A clear linear dependence of the antiskyrmion velocity on the current density emerges in both cases as shown in Fig. 1g. Furthermore, a significant enhancement in antiskyrmion velocity is present for the j⊥ V-Str case, suggesting that the confining magnetic stripe domains may play a role similar to the sample edge in the confined skyrmion dynamics discussed above27.

The formations of antiskyrmions confined within helical stripe domains in FNPP

To experimentally confirm the antiskyrmion dynamics predicted by the aforementioned micromagnetic simulations in FNPP, we fabricated a two-terminal device using a focused ion beam (FIB) method (see Supplementary Note 1). The experimental setup is schematically illustrated in Fig. 2a. We manufactured a rectangular FNPP device, with its edges aligned along the [11¯0] and [110] crystal axes, respectively, defined as the +x and +y directions in Fig. 2a. The positive pulse current direction is defined as the +x axis. Figure 2b shows the orientation of the device on a double-tilt holder capable of applying electric current to the device within the transmission electron microscope. We then employed LTEM to characterize magnetic configurations and their dynamics within the device.Fig. 2 In situ Lorentz transmission electron microscopy (LTEM) measurements of current-driven antiskyrmion motion in stripes.

a A scanning electron microscope (SEM) image of a two-terminal device. The scale bar in (a) is 1 μm. b The schematic of the device orientations with respect to the electron beam and external field. c, f Schematic drawings show the procedure to generate antiskyrmions (d, g) or half antiskyrmion pairs (HAPs) (e, h) in the H-Str (d, e) or in V-Str (g, h). The scale bar for images in (d, e, g, h) is 0.2 μm. i, k LTEM images observed after sequential current-pulse applications showing the current-driven antiskyrmion (marked by the red rectangular dashed line) motion in H-Str (i) and in V-Str (k) for positive current (i, left panels; k, top panels) and μ0H=−50mT. Micromagnetic simulation results for negative (j, right panels; l, bottom panels) current and μ0H=−50mT. The current densities in (i–l) are 0.95 × 1011 A m−2, 0.20 × 1011 A m−2, 0.73 × 1011 A m−2, and 1.00 × 1011 A m−2, respectively, and the current-pulse width is 50 ns. The scale bar for images in (i–l) is 0.5 μm. m, n The corresponding average antiskyrmion velocity v¯ versus current density j in (i, k). The error bars in both plots represent the standard deviation from the mean of multiple trials.

To create an antiskyrmion embedded in H-Str, the following experimental procedure (illustrated in Fig. 2c) was employed: (1) tilt the sample holder by approximately 10° (the θ-value), (2) increase the magnetic field to μ0H=−400mT, (3) return the sample holder to its original orientation with a roughly 2° tilt, and (4) reduce μ0H to 0 mT and tilt the sample holder to 0°. An antiskyrmion embedded in a H-Str is obtained, as shown in a defocused LTEM micrograph (Fig. 2d). Similarly, an antiskyrmion embedded in V-Str can be generated by tilting the device plate at the angle φ, as shown in Fig. 2f. Representative antiskyrmion is shown in Fig. 2g.

For a rectangular antiskyrmion, we can define a major axis length lAsky=nλhel, where λhel is the helical period of the studied system. Interestingly, in addition to typical square-shape antiskyrmion (n = 1), antiskyrmions with different elongation factor n within helical stripe domains (see Supplementary Notes 4 and 5) can also be created. These antiskyrmions with different n exhibit different current-driven dynamics, which will be discussed in the next section.

The geometry, unidirectional current, and external field effects on antiskyrmion dynamics

To drive the antiskyrmion dynamics while minimizing the Joule heating, we employed a short current-pulse width of 50 ns (see Note 2 in Supplementary). The current pulse is first applied to a single square-shaped antiskyrmions. Although n = 1 antiskyrmions can be readily generated within a single H-Str stripe, the critical current density required to move them is close to that beyond which the antiskyrmion will be annihilated (see Note 3 in Supplementary).

In contrast, the elongated antiskyrmions (n > 1) with their longer axis perpendicular to the background helices, can move at a lower threshold current density than that required for driving antiskyrmions with the square shape. For example, Fig. 2i shows that an antiskyrmion with n=3 embedded in H-Str moves at |j| = 0.95 × 1011 A m−2 (pulse width: 50 ns). The motion of antiskyrmions (n=3) embedded in V-Str under μ0H=−50mT and a 50 ns current-pulse width with flowing currents of |j| = 0.73 × 1011 A m−2 are also shown in the LTEM micrographs in Fig. 2k. A significant enhancement of vd is observed when the elongated antiskyrmion is embedded in V-Str (j ⊥ V-Str) compared to that in H-Str (j || H-Str).

To better understanding the antiskyrmion dynamics embedded in both H-Str and V-Str cases, micromagnetic simulations were performed for antiskyrmions (n=3), as shown in Fig. 2j, l. The simulation results are well consistent with experimental observations of antiskyrmion dynamics upon the current reversals, as shown in Supplementary Note 6. The unidirectional motion perpendicular to the current flow and along the V-Str (vd⊥j) can be ascribed to the Magnus force, which is perpendicular to the current flow as described theoretically with the Thiele equation discussed above27. The experimental results agree well with the micromagnetic simulation results to fix the formula as vd ∝ −(4πQ/α) j for the antiskyrmion in V-Str (j ⊥ V-Str), compared with vd ∝ (β/α) j for the antiskyrmion in H-Str (j∥H-Str). As the damping parameter α and non-adiabatic spin-transfer torque strength parameter β are usually very small and close to each other, we could have difference with 4π/α = 251 and β/α = 1, if we assume α = β = 0.05. However, the simulation for elongated antiskyrmion (n = 3) in j ⊥ V-Str shows a velocity that is only 1.5 times compared to that in j || H-Str, as shown in Supplementary Note 7. This difference may possibly ascribe to the difference in other parameters, such as the dissipative tensor, which is subject to the antiskyrmion geometry.

To quantitatively understand the relationship between vd and j in these two cases, we performed a systemically study of the antiskyrmion motion over a wide range of j in pinned and flow regimes for antiskyrmions. Figure 2m, n shows the distinctive average velocity of the antiskyrmion (n=3) as the j∥H-Str and j ⊥ V-Str, respectively. Notably, the results indicate a larger threshold current density (|jc| ~ 0.73 × 1011 A m−2) for driving the antiskyrmion in H-Str than that for the antiskyrmion in V-Str (|jc| ~ 0.49 × 1011 A m−2). At j = 0.94 × 1011 A m−2, it shows the velocity of elongated antiskyrmion (n = 3) in j ⊥ V-Str is about 9.3 times faster than that in j || H-Str. The results reveal that the current flow direction is a critical factor in designing antiskyrmion-based spintronics devices with confined geometries.

The effect of pulse width on antiskyrmion dynamics has also been investigated experimentally. We examined the elongated antiskyrmions (n = 4) with respect to the pulse width at 0 mT in both j∥H-Str and j ⊥ V-Str configuration (Supplementary Note 2). The results indicate a non-linear increase in antiskyrmion velocity with increasing pulse width. When the current exceeds a critical value, the transformation of it was observed that antiskyrmions can transition into non-topological bubbles, or be annihilation of antiskyrmion occurs.

Figure 3a shows the average antiskyrmion velocity ∣v¯y∣ versus n for antiskyrmions embedded in V-Str at j = 0.78 × 1011 A m−2. We observed that antiskyrmions with n = 2 and n = 3 exhibit onset motion above |jc| ~ 0.5 × 1011A m−2 and that the increase of n leads to faster motions of antiskyrmions. The same geometric effect on the dynamical antiskyrmions (n = 2 and n = 3) embedded in H-Str has also been revealed with flowing current of |jc| ~ 0.7 × 1011 A m−2 (see Supplementary Note 8).Fig. 3 Effects of antiskyrmion shape and external field on antiskyrmion dynamics.

a The average antiskyrmion velocity in a V-Str background versus the elongation factor n defined as the ratio of the antiskyrmion length lAsky to the helical period λhel at j = 0.78 × 1011 A m−2 under μ0H=−50mT. b The average velocity of n=3 antiskyrmions versus μ0H at j = 0.78 × 1011 A m−2. The inset images in (a, b) are LTEM images of n = 1, 2, and 3 antiskyrmions, respectively. The scale bar for inset images in (a, b) is 0.5 μm. The error bars in both plots represent the standard deviation from the mean of multiple trials.

In addition to the geometric effect on antiskyrmion dynamics, the external field H has been observed to modify the size and stability of spin textures34, thereby potentially affecting their dynamics35 (Supplementary Notes 6). The observed results reveal that above jc the velocity of antiskyrmion continues increasing in proportion to j, irrespective of the direction of μ0H. However, an increased magnetic field results in faster antiskyrmion motions, as shown in Fig. 3b. Our simulation results also confirmed this phenomenon that the antiskyrmion at −100 mT shows 15% velocity enhancement compared to the one at −50 mT in j ⊥ V-Str, however, the velocity for antiskyrmion in j || H-Str remains almost unchanged at −100 mT and −50 mT (Supplementary Notes 9). These results show that the antiskyrmion velocity can be tuned by both their geometric shapes and the external magnetic field. As for higher field, e.g., −150 mT, we found antiskyrmions are easy to transform into non-topological bubbles under pulse current stimulations and therefore, leading to the complicated dynamical phenomena. We have also conducted the experiment for current-driven isolated antiskyrmion in the ferromagnetic background. Such an isolated antiskyrmion transformed into a non-topological bubble before its motion driven by currents (Supplementary Notes 10).

Current-driven dynamics of half-antiskyrmion pairs

In addition to the formation of an antiskyrmion in stripes with the method in Fig. 2c, f, a half-antiskyrmion pair (HAP) embedded in the H-Str and V-Str shows up, as shown in Fig. 2e, h, respectively. Here, the HAP, which has the same topological charge (Q = +1) as the antiskyrmion, can be regarded as a combination of two antimerons36. We also perform the topological charge calculation of the HAP and antiskyrmion using micromagnetic simulation (see Supplementary Note 11). Note here that HAPs are distinguished with half skyrmion pairs16,25 or H-shaped skyrmion26 due to the presence of four Bloch lines in real space, even though they may carry the same unit of topological charge. The critical current density required for moving the HAPs embedded in H-Str is also as high as that needed to annihilate them, similarly like the antiskyrmion in the H-Str (see Supplementary Note 3). Unlike antiskyrmions embedded in V-Str, HAPs in V-Str exhibit current-driven motion with a markedly enhanced velocity and lower threshold current density. Figure 4a presents representative LTEM images of a HAP configuration at |j| = 0.63 × 1011 A m−2 under μ0H=−50mT, revealing that HAPs move along the same direction as their antiskyrmion counterparts under an electric current flow. To elucidate the difference between the current-driven motions of HAPs and antiskyrmions, we conducted an experiment probing the HAP velocity as a function of current density. Figure 4b shows a critical current density of |jc| ~ 0.4 × 1011 A m−2 for moving the HAP, about 20% lower than that for driving elongated antiskyrmions (n = 3) in j ⊥ V-Str. In contrast, a squared antiskyrmion in a V-Str domain exhibits step-by-step motion (Fig. 4c) with sequences of pulse currents at relatively lower current densities (Fig. 4d) compared to the drive of a squared antiskyrmion in the H-Str. Interestingly, the current-induced deformation of antiskyrmion is observed from a square shape into a rectangle shape along the current-flow direction due to the strong pinning effect in the present device (see Supplementary Note 8). Moreover, the HAP velocity at j = 0.78 × 1011 A m−2 is 16.0 m s−1, 49% larger than the antiskyrmion (n=3) velocity with j ⊥ V-Str.Fig. 4 The current-driven dynamics of a HAP and an antiskyrmion (n = 1).

a, c LTEM images of a HAP and an antiskyrmion observed by sequential current-pulse applications, respectively. The scale bar for images in (a, c) is 0.5 μm. The |j| is 0.63 × 1011 A m−2 and 0.78 × 1011 A m−2 in (a, c), flowing along the +x (upper panels) and opposite direction (−x, lower panels). A simulated LTEM image for the HAP is shown in the upper right panel of (a). The scale bar for the simulated LTEM image is 0.2 μm. b, d The HAP and antiskyrmion velocity averaged over several trials plotted as a function of j, respectively. All data were obtained under μ0H=−50mT. The error bars in both plots represent the standard deviation from the mean of multiple trials.

To gain deeper insights into the underlying mechanisms governing HAP dynamics, we also conducted experiments at various μ0H (Supplementary Note 12). The results reveal that the μ0H can also affect HAPs dynamics. However, it is noteworthy that the velocity of HAPs decreases with the increase in μ0H, as opposed to that observed for antiskyrmions. The reason could be that the increase in μ0H leads to a compression between adjacent half antiskyrmions, and thus reduces the size of the HAP. The efficiency of STTs in driving the HAP decreases with the decreasing size of HAP.

In summary, we have experimentally demonstrated the real-space observations of current-driven antiskyrmion motions in FNPP at room temperature. The antiskyrmions can be reproducibly generated within a stripe domain background and electrically driven while confined within the stripe domains. The velocity of antiskyrmions shows a linear dependence on the current density with current flows parallel or perpendicular to the stripes. Importantly, the antiskyrmion motion is enhanced when the current is flowed perpendicular to the stripes. Our results evidence antiskyrmion as a prototype carrier for racetrack-type memory37–39. These findings can potentially propel antiskyrmions to the forefront of future practical spintronic applications.

Methods

Sample preparation

A single crystal of (Fe0.63Ni0.3Pd0.07)3P (FNPP) was synthesized via the self-flux method utilizing Fe, Ni, Pd, and red phosphorous P. The two terminal devices composed of (001) thin plates of FNPP with dimensions of 0.16 μm × 4.9 μm × 6 μm for one device and 0.16 μm × 4.3 μm × 2.5 μm for another device, were fabricated with the focused ion beam technique (NB5000, Hitachi).

In situ measurement

Real-space observations of antiskyrmions were performed by transmission electron microscopes (JEM-2100F and JEM-2800, JEOL) equipped with a double-tilt current capable holder (Fusion select, Protochips). An external magnetic field was applied normally to the device-plate plane by controlling the objective-lens current of the microscope. The pulsed voltage was provided by an arbitrary function generator (AFG 31252, Tektronix). The pulse currents were manually triggered at intervals exceeding 2 s.

Data analysis

The LTEM videos are recorded during the sequential current-pulse excitations. The average velocities of antiskyrmions and HAPs are calculated over several trials (at least 4 trials) at each current density. The error bars of the velocity represent the standard deviation from the mean of multiple trials.

Micromagnetic simulations

Micromagnetic simulations were performed using the functional micromagnetics package−Object Oriented MicroMagnetic Framework (OOMMF)40. The current-driven spin dynamics were simulated by solving the Landau-Lifshitz-Gilbert equation augmented with adiabatic and non-adiabatic spin-transfer torques:2 dMdt=−γ0M××Heff+αMSM××dMdt+uMS2M×∂M∂x×M−βuMSM×∂M∂x

where M is the magnetization, MS = |MS| is the saturation magnetization, t is the time, γ0 is the absolute gyromagnetic ratio, and α is the damping parameter. Heff=−μ0−1∂ε/∂M is the effective field, and u=γ0ℏ/μ0ejP/2MS is the conduction electron velocity, where μ0 is the vacuum permeability constant, ℏ is the reduced Planck constant, e is the electron charge, j is the applied current density, and P is the spin polarization rate. The total system energy includes the ferromagnetic exchange interaction, anisotropic DMI, uniaxial magnetic anisotropy, external magnetic field, and demagnetization energy terms. The average energy density is expressed as:3 ε=A∇MMS2+DMS2M⋅x^×∂M∂x−M⋅y^×∂M∂y−KMS2z^⋅M2−μ0M⋅H-μ02M⋅Hd,

where A and D denote the ferromagnetic exchange and anisotropic DMI constants, respectively. K is the uniaxial magnetic anisotropy constant. H is the applied external magnetic field. Hd is the demagnetization field.

The intrinsic magnetic parameters are adopted from ref. 20 the saturation magnetization Ms = 417 kA m−1, the exchange constant A = 8.1 pJ m−1, the anisotropic DMI constant D = 0.25 mJ m−2, and the uniaxial anisotropy constant K = 3.1 × 104 J m−3. The default sample size is 2000 × 600 × 160 nm3 or 600 × 2000 × 160 nm3, which depends on the current injection direction. The mesh size is set to 5 × 5 × 32 nm3. The periodic boundary condition is applied in the x or y direction, depending on the current injection geometry. We assumed that the damping parameter α = 0.05 and β = 0.05. The spin polarization rate is set to P = 1 for simplicity, so that the driving force is simply proportional to the applied current density, i.e., u ~ j. The initial state is an antiskyrmion or HAP placed at the center of the sample with a helical background, which is relaxed at an out-of-plane field of Bz = −50 mT or 50 mT before flowing the current through the device.

Supplementary information

Supplementary information

Description of Additional Supplementary Files

Supplementary Movie1

Supplementary Movie2

Supplementary Movie3

Supplementary Movie4

Peer Review File

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52072-4.

Acknowledgements

The authors thank W. Koshibae, S. Mori, and Y. Chiew for the helpful discussions. This work was supported in part by Grants-In-Aid for Scientific Research (Grant Nos. 19H00660, 20H00337, 23H05431, 23H04522, 23H01841, 24H02231, 24H00197, 24H02231) from the Japan Society for the Promotion of Science (JSPS), by the Japan Science and Technology Agency (JST) CREST program (Grant Nos. JPMJCR1874, JPMJCR20T1), Japan, by the RIKEN TRIP initiative, and by the RIKEN TRIP Advanced General Intelligence for Science Program (RIKEN Quantum). M.M. also acknowledges support by the Waseda University Grant for Special Research Projects (Grant No. 2024C-153).

Author contributions

Y.G. and X.Z. Yu conceived the project. Y.G. fabricated the bulk materials, performed the Lorentz TEM measurements, analyzed LTEM data with L.C.P., F.S.Y. and X.Z. Yu. Y.G., F.S.Y. and X.Z. Yu wrote the manuscript. X.Z. and M.M. carried out the micromagnetic simulations assisted with data analyses. K.K. and Y. Taguchi prepared the FNPP alloy. D.N. measured the damping parameter of FNPP using ferromagnetic resonance. All authors discussed the data and commented on the manuscript.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

The authors declare that the main data supporting the findings of this study are available within the article and its Supplementary Information files. Extra data are available from the corresponding author upon reasonable request.

Code availability

The micromagnetic simulator OOMMF used in this work is publicly accessible at http://math.nist.gov/oommf.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Yao Guang, Xichao Zhang.
==== Refs
References

1. Mühlbauer S Skyrmion lattice in a chiral magnet Science 2009 323 915 918 10.1126/science.1166767 19213914
Mühlbauer, S. et al. Skyrmion lattice in a chiral magnet. Science 323, 915–918 (2009).19213914 10.1126/science.1166767
2. Nagaosa N Tokura Y Topological properties and dynamics of magnetic skyrmions Nat. Nanotechnol. 2013 8 899 911 10.1038/nnano.2013.243 24302027
Nagaosa, N. & Tokura, Y. Topological properties and dynamics of magnetic skyrmions. Nat. Nanotechnol. 8, 899–911 (2013).24302027 10.1038/nnano.2013.243
3. Yu X Real-space observation of a two-dimensional skyrmion crystal Nature 2010 465 901 904 10.1038/nature09124 20559382
Yu, X. et al. Real-space observation of a two-dimensional skyrmion crystal. Nature 465, 901–904 (2010).20559382 10.1038/nature09124
4. Nayak AK Magnetic antiskyrmions above room temperature in tetragonal Heusler materials Nature 2017 548 561 566 10.1038/nature23466 28846999
Nayak, A. K. et al. Magnetic antiskyrmions above room temperature in tetragonal Heusler materials. Nature 548, 561–566 (2017).28846999 10.1038/nature23466
5. Jonietz F Spin transfer torques in MnSi at ultralow current densities Science 2011 333 1381 1381 21903795
Jonietz, F. et al. Spin transfer torques in MnSi at ultralow current densities. Science 333, 1381–1381 (2011).21903795
6. Moriya T Anisotropic superexchange interaction and weak ferromagnetism Phys. Rev. 1960 120 91 98 10.1103/PhysRev.120.91
Moriya, T. Anisotropic superexchange interaction and weak ferromagnetism. Phys. Rev. 120, 91–98 (1960).10.1103/PhysRev.120.91
7. Wang Z Thermal generation, manipulation and thermoelectric detection of skyrmions Nat. Electron. 2020 3 672 679 10.1038/s41928-020-00489-2
Wang, Z. et al. Thermal generation, manipulation and thermoelectric detection of skyrmions. Nat. Electron. 3, 672–679 (2020).10.1038/s41928-020-00489-2
8. Zázvorka J Thermal skyrmion diffusion used in a reshuffler device Nat. Nanotechnol. 2019 14 658 661 10.1038/s41565-019-0436-8 31011220
Zázvorka, J. et al. Thermal skyrmion diffusion used in a reshuffler device. Nat. Nanotechnol. 14, 658–661 (2019).31011220 10.1038/s41565-019-0436-8
9. Yu X Real-space observations of 60-nm skyrmion dynamics in an insulating magnet under low heat flow Nat. Commun. 2021 12 5079 10.1038/s41467-021-25291-2 34426575
Yu, X. et al. Real-space observations of 60-nm skyrmion dynamics in an insulating magnet under low heat flow. Nat. Commun. 12, 5079 (2021).34426575 10.1038/s41467-021-25291-2
10. Yu XZ Transformation between meron and skyrmion topological spin textures in a chiral magnet Nature 2018 564 95 98 10.1038/s41586-018-0745-3 30518889
Yu, X. Z. et al. Transformation between meron and skyrmion topological spin textures in a chiral magnet. Nature 564, 95–98 (2018).30518889 10.1038/s41586-018-0745-3
11. Tang J Magnetic skyrmion bundles and their current-driven dynamics Nat. Nanotechnol. 2021 16 1086 1091 10.1038/s41565-021-00954-9 34341518
Tang, J. et al. Magnetic skyrmion bundles and their current-driven dynamics. Nat. Nanotechnol. 16, 1086–1091 (2021).34341518 10.1038/s41565-021-00954-9
12. Peng L Dynamic transition of current-driven single-skyrmion motion in a room-temperature chiral-lattice magnet Nat. Commun. 2021 12 6797 10.1038/s41467-021-27073-2 34819505
Peng, L. et al. Dynamic transition of current-driven single-skyrmion motion in a room-temperature chiral-lattice magnet. Nat. Commun. 12, 6797 (2021).34819505 10.1038/s41467-021-27073-2
13. Wang W Electrical manipulation of skyrmions in a chiral magnet Nat. Commun. 2022 13 1593 10.1038/s41467-022-29217-4 35332156
Wang, W. et al. Electrical manipulation of skyrmions in a chiral magnet. Nat. Commun. 13, 1593 (2022).35332156 10.1038/s41467-022-29217-4
14. Yu X Motion tracking of 80-nm-size skyrmions upon directional current injections Sci. Adv. 2020 6 eaaz9744 10.1126/sciadv.aaz9744 32596450
Yu, X. et al. Motion tracking of 80-nm-size skyrmions upon directional current injections. Sci. Adv. 6, eaaz9744 (2020).32596450 10.1126/sciadv.aaz9744
15. Yasin FS Real-space determination of the isolated magnetic skyrmion deformation under electric current flow Proc. Natl. Acad. Sci. USA 2022 119 e2200958119 10.1073/pnas.2200958119 36191237
Yasin, F. S. et al. Real-space determination of the isolated magnetic skyrmion deformation under electric current flow. Proc. Natl. Acad. Sci. USA 119, e2200958119 (2022).36191237 10.1073/pnas.2200958119
16. Song, D. et al. Steady motion of 80-nm-size skyrmions in a 100-nm-wide track. Nat. Commun. 15, 5614 (2024).
17. Tokura Y Kanazawa N Magnetic skyrmion materials Chem. Rev. 2020 121 2857 2897 10.1021/acs.chemrev.0c00297 33164494
Tokura, Y. & Kanazawa, N. Magnetic skyrmion materials. Chem. Rev. 121, 2857–2897 (2020).33164494 10.1021/acs.chemrev.0c00297
18. Peng L Controlled transformation of skyrmions and antiskyrmions in a non-centrosymmetric magnet Nat. Nanotechnol. 2020 15 181 186 10.1038/s41565-019-0616-6 31959930
Peng, L. et al. Controlled transformation of skyrmions and antiskyrmions in a non-centrosymmetric magnet. Nat. Nanotechnol. 15, 181–186 (2020).31959930 10.1038/s41565-019-0616-6
19. He Z Visualizing emergent magnetic flux of antiskyrmions in Mn1. 4PtSn magnet Adv. Funct. Mater. 2022 32 2112661 10.1002/adfm.202112661
He, Z. et al. Visualizing emergent magnetic flux of antiskyrmions in Mn1. 4PtSn magnet. Adv. Funct. Mater. 32, 2112661 (2022).10.1002/adfm.202112661
20. Karube K Room-temperature antiskyrmions and sawtooth surface textures in a non-centrosymmetric magnet with S4 symmetry Nat. Mater. 2021 20 335 340 10.1038/s41563-020-00898-w 33495630
Karube, K. et al. Room-temperature antiskyrmions and sawtooth surface textures in a non-centrosymmetric magnet with S4 symmetry. Nat. Mater. 20, 335–340 (2021).33495630 10.1038/s41563-020-00898-w
21. Karube K Doping control of magnetic anisotropy for stable antiskyrmion formation in schreibersite (Fe,Ni)(3) P with S(4) symmetry Adv. Mater. 2022 34 e2108770 10.1002/adma.202108770 35032408
Karube, K. et al. Doping control of magnetic anisotropy for stable antiskyrmion formation in schreibersite (Fe,Ni)(3) P with S(4) symmetry. Adv. Mater. 34, e2108770 (2022).35032408 10.1002/adma.202108770
22. Yasin FS Bloch lines constituting antiskyrmions captured via differential phase contrast Adv. Mater. 2020 32 e2004206 10.1002/adma.202004206 33043519
Yasin, F. S. et al. Bloch lines constituting antiskyrmions captured via differential phase contrast. Adv. Mater. 32, e2004206 (2020).33043519 10.1002/adma.202004206
23. Yasin, F. S. et al. Bloch point quadrupole constituting hybrid topological strings revealed with electron holographic vector field tomography. Adv. Mater. 36, 2311737 (2024).
24. Zheng F Skyrmion–antiskyrmion pair creation and annihilation in a cubic chiral magnet Nat. Phys. 2022 18 863 868 10.1038/s41567-022-01638-4
Zheng, F. et al. Skyrmion–antiskyrmion pair creation and annihilation in a cubic chiral magnet. Nat. Phys. 18, 863–868 (2022).10.1038/s41567-022-01638-4
25. Müller J Magnetic skyrmions and skyrmion clusters in the helical phase of Cu2OSeO3 Phys. Rev. Lett. 2017 119 137201 10.1103/PhysRevLett.119.137201 29341720
Müller, J. et al. Magnetic skyrmions and skyrmion clusters in the helical phase of Cu2OSeO3. Phys. Rev. Lett. 119, 137201 (2017).29341720 10.1103/PhysRevLett.119.137201
26. Knapman, R., Rodrigues, D. R., Masell, J. & Everschor-Sitte, K. Current-induced H-shaped-skyrmion creation and their dynamics in the helical phase. J. Phys. D Appl. Phys. 54, 404003 (2021).
27. Jiang W Direct observation of the skyrmion Hall effect Nat. Phys. 2016 13 162 169 10.1038/nphys3883
Jiang, W. et al. Direct observation of the skyrmion Hall effect. Nat. Phys. 13, 162–169 (2016).10.1038/nphys3883
28. Huang, S. et al. Stabilization and current-induced motion of antiskyrmion in the presence of anisotropic Dzyaloshinskii-Moriya interaction. Phys. Rev. B 96, 144412 (2017).
29. Sivakumar PK Topological Hall signatures of two chiral spin textures hosted in a single tetragonal inverse Heusler thin film ACS Nano 2020 14 13463 13469 10.1021/acsnano.0c05413 32986403
Sivakumar, P. K. et al. Topological Hall signatures of two chiral spin textures hosted in a single tetragonal inverse Heusler thin film. ACS Nano 14, 13463–13469 (2020).32986403 10.1021/acsnano.0c05413
30. Iwasaki J Mochizuki M Nagaosa N Current-induced skyrmion dynamics in constricted geometries Nat. Nanotechnol. 2013 8 742 747 10.1038/nnano.2013.176 24013132
Iwasaki, J., Mochizuki, M. & Nagaosa, N. Current-induced skyrmion dynamics in constricted geometries. Nat. Nanotechnol. 8, 742–747 (2013).24013132 10.1038/nnano.2013.176
31. Sampaio J Cros V Rohart S Thiaville A Fert A Nucleation, stability and current-induced motion of isolated magnetic skyrmions in nanostructures Nat. Nanotechnol. 2013 8 839 844 10.1038/nnano.2013.210 24162000
Sampaio, J., Cros, V., Rohart, S., Thiaville, A. & Fert, A. Nucleation, stability and current-induced motion of isolated magnetic skyrmions in nanostructures. Nat. Nanotechnol. 8, 839–844 (2013).24162000 10.1038/nnano.2013.210
32. Iwasaki J Mochizuki M Nagaosa N Universal current-velocity relation of skyrmion motion in chiral magnets Nat. Commun. 2013 4 1463 10.1038/ncomms2442 23403564
Iwasaki, J., Mochizuki, M. & Nagaosa, N. Universal current-velocity relation of skyrmion motion in chiral magnets. Nat. Commun. 4, 1463 (2013).23403564 10.1038/ncomms2442
33. Iwasaki J Koshibae W Nagaosa N Colossal spin transfer torque effect on skyrmion along the edge Nano Lett. 2014 14 4432 4437 10.1021/nl501379k 24988528
Iwasaki, J., Koshibae, W. & Nagaosa, N. Colossal spin transfer torque effect on skyrmion along the edge. Nano Lett. 14, 4432–4437 (2014).24988528 10.1021/nl501379k
34. Peng L Formation and control of zero-field antiskyrmions in confining geometries Adv. Sci. 2022 9 e2202950 10.1002/advs.202202950
Peng, L. et al. Formation and control of zero-field antiskyrmions in confining geometries. Adv. Sci. 9, e2202950 (2022).10.1002/advs.202202950
35. Koshibae W Nagaosa N Theory of antiskyrmions in magnets Nat. Commun. 2016 7 10542 10.1038/ncomms10542 26821932
Koshibae, W. & Nagaosa, N. Theory of antiskyrmions in magnets. Nat. Commun. 7, 10542 (2016).26821932 10.1038/ncomms10542
36. Yang H Wang Q Wen L Ground states of dipolar spin–orbit-coupled Bose–Einstein condensates in a toroidal trap J. Phys. Soc. Jpn. 2019 88 064001 10.7566/JPSJ.88.064001
Yang, H., Wang, Q. & Wen, L. Ground states of dipolar spin–orbit-coupled Bose–Einstein condensates in a toroidal trap. J. Phys. Soc. Jpn. 88, 064001 (2019).10.7566/JPSJ.88.064001
37. Parkin S Hayashi M Thomas L Magnetic domain-wall racetrack memory Science 2008 320 190 194 10.1126/science.1145799 18403702
Parkin, S., Hayashi, M. & Thomas, L. Magnetic domain-wall racetrack memory. Science 320, 190–194 (2008).18403702 10.1126/science.1145799
38. Fert A Cros V Sampaio J Skyrmions on the track Nat. Nanotechnol. 2013 8 152 156 10.1038/nnano.2013.29 23459548
Fert, A., Cros, V. & Sampaio, J. Skyrmions on the track. Nat. Nanotechnol. 8, 152–156 (2013).23459548 10.1038/nnano.2013.29
39. Parkin S Yang S-H Memory on the racetrack Nat. Nanotechnol. 2015 10 195 198 10.1038/nnano.2015.41 25740128
Parkin, S. & Yang, S.-H. Memory on the racetrack. Nat. Nanotechnol. 10, 195–198 (2015).25740128 10.1038/nnano.2015.41
40. Donahue, M. J. & Porter, D. G. OOMMF User’s Guide, Version 1.0. (National Institute of Standards and Technology Gaithersburg, 1999).
