Room-temperature current-induced generation and motion of sub-100nm skyrmions

Room-temperature current-induced generation and motion of sub-100nm skyrmions


Magnetic skyrmions are nanoscale windings of the spin configuration that hold great promise for technology due to their topology-related properties and extremely reduced sizes. After the recent observation at room temperature of sub-100 nm skyrmions stabilized by interfacial chiral interaction in magnetic multilayers, several pending questions remain to be solved, notably about the means to nucleate individual compact skyrmions or the exact nature of their motion. In this study, a method leading to the formation of magnetic skyrmions in a micrometer-sized nanotrack using homogeneous current injection is evidenced. Spin-transfer-induced motion of these small electrical-current-generated skyrmions is then demonstrated and the role of the out-of-plane magnetic field in the stabilization of the moving skyrmions is also analysed. The results of these experimental observations of spin torque induced motion are compared to micromagnetic simulations reproducing a granular type, non-uniform magnetic multilayer, in order to address the particularly important role of the magnetic inhomogeneities on the current-induced motion of sub-100 nm skyrmions, for which the material grains size is comparable to the skyrmion diameter.

Skyrmion, spin–orbit torque, room temperature imaging, magnetic multilayers, current–induced motion, micromagnetics

The recent discovery of room temperature skyrmions and larger, skyrmonic bubbles in ultrathin magnetic materials Jiang et al. (2015); Chen et al. (2015); Moreau-Luchaire et al. (2016); Boulle et al. (2016); Yu et al. (2016); Woo et al. (2016) opens the way to many fundamental studies on topological spin textures, and raises the expectations of using their topology-related properties Bogdanov and Rößler (2001); Nagaosa and Tokura (2013) in room temperature devices Fert et al. (2013); Kiselev et al. (2011). Among the different types of spin textures that can be stabilized in ferromagnetic thin films, skyrmions are believed to offer the highest potential to be reduced in size down to the nanometer scale Kiselev et al. (2011); Fert et al. (2013); Sampaio et al. (2013). Their extremely small diameter Romming et al. (2013, 2015) is particularly interesting for achieving high information density in magnetic data storage and other spintronic devices. Isolated skyrmions can notably be obtained at room temperature in magnetic multilayers lacking the inversion symmetry of their stacking order Jiang et al. (2015); Chen et al. (2015); Moreau-Luchaire et al. (2016); Boulle et al. (2016); Jiang et al. (2016, 2017); Woo et al. (2016); Yu et al. (2016); Litzius et al. (2017); Hrabec et al. (); Sampaio et al. (2013); Fert et al. (2013), in which the interfacial Dzyaloshinskii-Moriya interaction (DMI) stabilizes chiral spin textures Fert (1990). In most experiments, the size of the skyrmion falls into the 100 nm–1 m range, resulting in a skyrmionic bubble, whose reversed core is not point-like but extends over a diameter which is much larger than its chiral domain wall width Kiselev et al. (2011). Only in specific cases (with regard to the balance between the DMI Dupé et al. (2016); Yang et al. (2015), magnetic anisotropy Yu et al. (2016), dipolar interactions Boulle et al. (2016), etc. ) the DMI can stabilize sub-100 nm skyrmions as demonstrated recently Moreau-Luchaire et al. (2016); Soumyanarayanan et al. (), for example, in Pt/Co/Ir multilayers.

Skyrmions in such symmetry-breaking magnetic multilayers exhibit a fixed chirality defined by the sign of the chiral interaction Bode et al. (2007); Heinze et al. (2011), resulting in chirality-dependent motion Thiaville et al. (2012); Emori et al. (2013) that can be driven by spin-orbit torques Jiang et al. (2015); Woo et al. (2016)(SOTs). In addition to that, the spin texture of the skyrmions exhibits fascinating topology-related effects Nagaosa and Tokura (2013), such as the topological Hall effect Neubauer et al. (2009); Schulz et al. (2012); Kanazawa et al. (2015) and the skyrmion Hall effect Yin et al. (2015), both of which have been demonstrated experimentally Schulz et al. (2012); Kanazawa et al. (2015); Jiang et al. (2017); Litzius et al. (2017). From a technological perspective, the predicted low propagation current densities Iwasaki et al. (2013a) and highly controllable magnetization dynamics of skyrmions Ma et al. (2015) provide unique opportunities for developing high density, low power spintronic nanodevices for storage and logic operations.

In this prospect, adressing the skyrmion creation and controlling the motion of small skyrmions in technologically relevant materials is of much interest. However most of previous studies have been focusing on large, micrometer-sized skyrmionic bubbles Jiang et al. (2015, 2016, 2017); Yu et al. (2016). In this Letter, we investigate the manipulation of sub-100 nm skyrmions, with regard to both nucleation and spin-torque-induced motion. In particular, we show that such small compact skyrmions can be nucleated by applying a uniform current directly into nanotracks, and subsequently be moved through spin-orbit torques. For the present range of small skyrmion diameters, they reach a size comparable to the grains size in our sputtered grown magnetic multilayers, which results into a regime of motion significantly influenced by the magnetic inhomogeneities associated to these grains. With the support of micromagnetic simulations, we identify parameters influencing the perturbed or even impeded motion of skyrmions due to inhomogeneities, in order to provide directions towards better nanotrack designs for more efficient devices.

For the presented experiments, three multilayered stacks were deposited at room temperature by d. c. sputtering under an Ar pressure of 2.5 mbar. The structures are based on Pt/Co/Ir trilayers, which enable additive DMI at the Pt/Co and Co/Ir interfaces, as shown in a recent work Moreau-Luchaire et al. (2016). The complete stack is: thermally oxidised Si substrate/Ta 15/Co /[Pt 1/Ir 1/Co /Pt 3, where is 0.6, 0.8 and 1.0 nm, all other thicknesses are also given in nanometers. The Ta layer is used as a buffer layer for growth but also as a source of spin torques. The Pt capping layer prevents oxidation of the stack. The multilayers were subsequently patterned into 1 m-wide nanotracks by electron beam lithography and Ar ion milling, while unpatterned reference extended films were kept for standard magnetometry characterization.

Figure 1: a. Hysteresis loops of the perpendicular magnetization component as a function of the applied perpendicular field for the three samples of 0.6 nm Co (red triangles), 0.8 nm Co (blue diamonds) and 1.0 nm Co (black circles). The larger, filled points in each loop indicate the field set for the MFM images of Fig. 2, after the field was swept from the opposite magnetic saturation as indicated by the colored arrows. b. MFM image of the 0.6 nm sample just after growth (as-grown configuration). Dashed circles highlight some of the individual skyrmions. c. MFM image of the 0.6 nm Co sample at remanence after a field cycle. The blue/red contrast in b. and c. indicates areas with positive/negative magnetization along z.

The hysteresis loops of the out-of-plane (z) component of the magnetization (normalized to the saturation magnetization ) vs. perpendicular field are shown for the three multilayers in Fig. 1a. All the multilayers exhibit perpendicular magnetic anisotropy (PMA), and their magnetization breaks into domains during magnetization reversal. It can be seen on the hysteresis loops that the saturation field increases with increasing Co thicknesses, indicative of an increasing influence of the dipolar coupling between the thicker Co layers. Figures 1b,c are magnetic force microscopy (MFM) mappings of the phase, acquired in lift mode at room temperature on the unpatterned multilayer with 0.6 nm. They show the magnetic domains (magnetization up or down) respectively in the as-grown magnetic configuration and at remanence after having cycled the perpendicular field to saturation and back to zero. As can be seen from the circular shaped domains indicated by dashed circles in Fig. 1b, the deposited stack spontaneously exhibits skyrmionic bubbles and skyrmions, due to the strong DMI Moreau-Luchaire et al. (2016), together with worm domains. However, these skyrmions are not present anymore in the images acquired at remanence after cycling the field (Fig. 1c). This is due to the preferential formation of a labyrinthine demagnetized multidomain pattern after the nucleation of the first domains of reversed magnetization. The present example shows that it is not always possible to get spontaneous nucleation of skyrmions and their stabilization during field cycles, but only in some specific cases Moreau-Luchaire et al. (2016); Yu et al. (2016). Even if in the present case, cycling the field from the saturated state does not allow us to obtain skyrmions in the multilayers, the as-grown configuration shows their metastability once present. In this study, we thus investigate another method for skyrmions creation by applying a uniform current in narrow tracks. In the 1 m-wide nanotracks, the magnetic domain configuration is influenced by the confinement Sampaio et al. (2013); Rohart and Thiaville (2013), which actually allows the current-assisted nucleation of skyrmions, as shown below.

Before injecting the current, we first saturate out-of-plane the magnetization of the patterned tracks. Then, the applied perpendicular field direction is reversed and its value slightly shifted from zero in order to favor one magnetic orientation over the other. For the multilayers with 0.6, 0.8 and 1.0 nm, after saturation, the reversed field is set to 17, 20.5 and 44 mT respectively, as shown by the arrows in Fig. 1a. The MFM images of the tracks with 0.6, 0.8 and 1.0 nm are shown in Figs. 2a-c, respectively. The absolute values of the field are chosen to be increasing with the Co layer thickness in order to keep the same ratio. These values correspond to the large, filled symbols in Fig. 1a. The field for 1.0nm is purposively chosen negative to show the equivalence between the cases of positive and negative field bias (this is also the reason for the reversed contrast). For a Co thickness of 0.6nm, the track is almost fully saturated (Fig. 2a); for 0.8nm of Co thickness, it exhibits separated, isolated domains arising from the track edges (Fig. 2b); and for 1.0nm of Co thickness, it hosts extended, connected worm domains (Fig. 2c).

Figure 2: MFM images of the 1 m-wide tracks before any pulses in the a. 0.6 nm, b. 0.8 nm and c. 1.0 nm thick Co multilayers, and after the injection of 1000 current pulses of 200 ns in the d. 0.6 nm, e. 0.8 nm and f. 1.0 nm thick Co multilayers. The applied fields are respectively 17, 20.5 and -44 mT and the pulses amplitudes are respectively 1.43, 2.38 and 2.96  A m for the three different Co thicknesses. g. Dependance of , the square of the threshold current density for skyrmion nucleation, as a function of the Co thickness . The line is a linear fit to the data points. h. Dependance of the pulse amplitude as a function of the duty cycle of the excitation for changing the magnetic configuration with a bipolar pulsed excitation at 2.5 MHz. The line is an offset inverse function fit to the data points.

To further manipulate the magnetic configuration in the tracks, our approach is to leverage the action of electrical current pulses. The required current density for skyrmion nucleation is experimentally found by alternatively imaging the magnetic configuration in the track and applying pulses of increasing voltage, corresponding to increasing current densities , up to the observation of changes in the magnetic texture. In Figs. 2d-f we respectively show the same tracks as in Figs. 2a-c after the application of 1000 pulses of duration 200 ns each and current density 1.43, 2.38 and 2.96 A m, respectively. Unless otherwise stated, all pulses in our experiments are separated by more than 1 ms to avoid accumulated Joule heating between pulses. Whatever the initial magnetization (almost saturated, isolated domains or worms), a configuration that exhibits magnetic skyrmions can eventually be obtained at or above . In Fig. 2g we display the required current density squared for skyrmion formation as a function of . We note that the nucleation (and ) was found identical (not shown) for the two opposite polarities of the pulses. The linear fit (made allowing an offset due to room temperature) in Fig. 2g indicates a good correlation of with , that is, of the thermal energy with the magnetic volume, suggesting a thermal origin of the nucleation. Under the application of pulses, the Joule heating raises the temperature of the track. This allows to escape the magnetic states of Figs. 2a-c, which are determined by the magnetic history of the sample, to then reach the skyrmion state of Figs. 2d-f, similar to the as-grown state observed in Fig. 1b. In order to confirm the thermal origin of the skyrmion nucleation, we set a periodic, bipolar pulse excitation at the frequency of 2.5 MHz running continuously with different duty cycles, and find for each of them the required pulse voltage to affect the magnetic configuration. As shown in Fig. 2h, we obtain that the squared voltage scales inversely with the duty cycle (which is proportional to the accumulated heat), again consistent with Joule heating. It is important to pinpoint that the described nucleation process is reproducible and does not affect the magnetic properties of the tracks, as this method was used tens of times in the same tracks to initialize the skyrmion state after saturation, without noticeable changes in the skyrmions size, nor in the domain characteristics after a field cycle. The generation of skyrmions by homogeneous currents in nanotracks is another mechanism for skyrmion preparation, different from the methods using a constriction and inhomogeneous current injection Jiang et al. (2015, 2016); Heinonen et al. (2016); Hrabec et al. (), and may be more easily implementable when reducing the size of the devices towards the deep sub-micrometer range for applications.

Once the nanotracks are prepared in a magnetic state which exhibits isolated skyrmions, the next step is to study their current-induced motion. To limit the nucleation or deletion of skyrmions due to Joule heating, once the skyrmions are nucleated, we use shorter duration (100 ns) current pulses in order to induce motion. As a check method, we verify that the magnetic configuration remains unchanged, with no creation nor deletion of magnetic domains even after 1000 pulses are injected at the previously found nucleation current density . This again supports that the previous mechanism of skyrmion creation was driven by the temporary temperature rise in the track. As the 100 ns pulses have no effect on the magnetic configuration at , it is possible to slightly increase the current density to get stronger spin torques until motion can eventually be observed. In particular, for the multilayer with the 0.8 nm thick Co layers (on which we will focus from now on), increasing the current density from 2.38 A m to 2.85 A m for 100 ns pulses allows the motion of the skyrmions. An out-of-plane bias field 18.5 mT is applied to stabilize even more the skyrmions in these new conditions. In Figs. 3a-d we display the initial state in the track, and the magnetic configuration as observed by MFM after 50, 100 and 150 pulses respectively, for pulses of one polarity (charge current towards the left, as indicated on the figure). In Figs. 3e-h we show the corresponding images for pulses of the opposite polarity (charge current towards the right, as indicated on the figure). The whole sets of images are provided in animated images 1 and 2 as supplementary files in Supporting Information. In Figs. 3i,j the position of the center of each domain is tracked after each burst of 10 pulses for current respectively towards the left and towards the right, with the color of the points evolving from black, for the initial frame, to red, after 50 frames. The motion of the skyrmions can be seen from the circled skyrmions in Figs. 3a-h. As highlighted by the blue arrows following the traces of the skyrmions, their trajectories reveal motion of the skyrmions to the left (right) in Fig. 3i (Fig. 3j). The direction of the motion thus follows the direction of the charge current for both opposite polarities. As the motion is opposed to the electron flow, the SOTs must be the principal source of motion in this structure Thiaville et al. (2012). The SOTs in such multilayers are induced by the vertical spin currents due to the spin Hall effect in the thick Ta buffer, while the spin Hall effect in the repeated Pt layers is expected to be small due to their thickness which is small (only 1 nm) in comparison to their spin diffusion length (3–4 nm measured in our films) Rojas-Sánchez et al. (2014).

Figure 3: MFM images of the skyrmions in the track, in the initial state and after 50, 100 and 150 pulses, respectively, under 18.5 mT with a-d. current towards the left and e-h. towards the right. The colored circles connected by dashed lines highlight some moving skyrmions. Traces of the skyrmion positions for i. current towards the left, and j. current towards the right, with the iterations shown by the color, from black (initial) to red (after 50 bursts of 10 pulses), as shown on the color scale of t. MFM images of the skyrmions in the track, in the initial state and after 50, 100 and 150 pulses, respectively, k-n. under 11.5 mT and o-r.  21 mT, with current towards the left. Traces of the domain positions for s.  11.5 mT, and skyrmion positions for t.  21 mT, with the iterations shown by the color, from black (initial) to red (after 50 bursts of 10 pulses), as shown on the color scale of t. All images and graphs share the same x-y scale. For each group of images the direction of the current density is indicated by the large, black arrow. 2.85 A m and pulses are 100 ns long in each case. The blue arrows in i,j. and s,t. indicate the direction of motion for the most mobile skyrmions.

The images of Figs. 3a-h show a few worm domains among the skyrmions present inside the track. Without current applied, the isolated skyrmions are expected to keep a reduced size because the out-of-plane field is opposed to the magnetization inside the skyrmion core and confines it. However, it is to be reminded that the SOTs can act against this field that keeps the skyrmions compact, in turn allowing them to extend along the current direction during the pulses and thus creating worm domains. To demonstrate the controllability of this transition from skyrmions to worms, we can change the external out-of-plane field. In Figs. 3k-n (resp. Figs. 3o-r) we show the initial state in the track, and the magnetic configuration as imaged by MFM after 50, 100 and 150 pulses respectively, with charge current towards the left at 11.5 mT (resp. 21 mT). The whole sets of images are provided in animated images 3 and 4 as supplementary files in Supporting Information. Again, 2.85 A m. We see that after these sequences of pulses, the track is mainly populated by worm domains for the smaller external field 11.5 mT, whereas the skyrmions keep a circular, compact shape and move without extending into worms for the larger field 21 mT. Adjusting the external field thus allows to favor or suppress the worm domains under current.

In Fig. 3s (resp. Fig. 3t) we summarize the trajectories of the domains after each burst of 10 pulses under 11.5 mT (resp. 21 mT). It appears that at 11.5 mT, we do not observe motion but rather an elongation of the worm domains due to the current. On the contrary at 21 mT, a significant fraction (about half) of the skyrmions are moving, but it should be noted that some other skyrmions are completely pinned, as it has been shown in similar experiments by other groups Woo et al. (2016). From the analysis of the moving skyrmion trajectories, a mean velocity of about 20–40 mm s can be extracted. As will be discussed later, we also note that (i) the motion of the moving skyrmions is far from being uniform; and (ii) the trajectories tend to be aligned towards the direction of the current, without evidence for a preferred direction of transverse motion, that is, the absence of skyrmion Hall effect, usually expected from the gyrotropic motion of the skyrmions Jiang et al. (2017); Hrabec et al. ().

Upon further increasing the current density of the pulses, nucleation and/or deletion of some skyrmions occur and prevent the clear observation of the faster motion. For example in the images at 21 mT of Figs. 3o-r, one can notice the nucleation of more skyrmions in the track upon repeating the pulses, as we count 31 skyrmions in Fig. 3o and already 50 skyrmions in Fig. 3r. Note that for the multilayers with 0.6 nm and 1.0 nm, the current-induced motion for 100 ns pulses was found to begin at a current density slightly above the current density required for nucleation. For this reason, we have not been able to perform a quantitative analysis of the current induced skyrmion motion in these multilayers, nor at higher current densities in the 0.8 nm multilayer.

The mean velocities (20–40 mm s) and peak velocities (up to 0.5 m s) found in our present experiments are rather low as compared to those of other recent works where skyrmions in magnetic multilayers exhibit motion in the flow regime. Jiang et al. have found a skyrmion velocity Jiang et al. (2015) of 25 m s at 4.5 A m, but then 0.75 m s overcoming the pinning with the use of higher current densities Jiang et al. (2017) up to 6 A m. In two very recent works by Woo et al. Woo et al. (2016) and Hrabec et al. Hrabec et al. (), the velocity of the skyrmion reaches respectively up to 100 m s and 60 m s, at current densities around 5 A m. Up to 2–2.5 A m, however, these two studies find a velocity very close to zero due to pinning effects. Considering that we are using current densities close to those at which skyrmions are almost pinned in other studies Woo et al. (2016); Hrabec et al. (), and that no skyrmion Hall effect is observed here consistently suggests that the flow regime is not reached Jiang et al. (2015, 2017) even at the maximal current density in our tracks. To qualitatively understand the motion that we experimentally observed, we performed micromagnetic simulations using the solver MuMax3 Vansteenkiste et al. (2014). The magnetic parameters were taken from a previous work Moreau-Luchaire et al. (2016): saturation magnetization 956 kA m, uniaxial anisotropy 717 kJ m, exchange parameter 10 pJ m, and damping parameter 0.2. For simplicity, the torque considered here is only the antidamping torque originating from the spin Hall effect in the Ta layer (and possibly a small contribution from the Pt layers), with an effective spin Hall angle 0.033 Liu et al. (2012); Qiu et al. (2014), and neglecting the field-like torque.

Figure 4: a. Bottom. Velocities in the x direction (i. e. along the current, right-pointing triangles) and in the y direction (i. e. transverse to the current, up-pointing triangles) as a function of the free-standing skyrmion diameter . The different points correspond to different values of 0–130 mT. The open symbols denote the first simulated size for which the skyrmion is not stable and progressively expands non-circularly. Top. Skyrmion wall width as a function of the skyrmion diameter. The dashed line is and the horizontal line is 6.47 nm as expected from the micromagnetic parameters. b. Velocity in the x direction (i. e. along the current) as a function of the diameter of a skyrmion confined along the edge of a track. In each graph the corresponding values of in mJ m are denoted by the different colors and given by the numbers in parenthesis. The coordinates system with the current direction and an example skyrmion trajectory are depicted in schematics as insets of a. and b. where the dotted lines are edges with periodic boundary conditions. The curent density is set to 3 A m.

The current-induced motion of skyrmions was first simulated in a multilayer made of 11 repetitions, similar to our experiment, taking into account the dipolar coupling between each layer (of identical magnetization profile) by adding them at periodic positions in the z direction. We model the case in which each Co layer is subject to a damping-like torque from an adjacent heavy metal Liu et al. (2012) where the charge current flows in the x direction. We treat two different cases: (a) a single skyrmion in an extended film (x and y periodic boundary conditions) and (b) a skyrmion confined in a track (periodic boundary conditions along x only). The simulation area comprising one skyrmion, its trajectory in the coordinates systems and the definition of the current direction are shown for each case in the insets of Figs. 4a and 4b. For simplicity, temperature was not considered in our simulations, even if it is expected to influence the stochasticity of the motion Iwasaki et al. (2014) and the skyrmion stability Rohart et al. (2016). The velocity of an individual skyrmion has been extracted from the simulations under different magnitudes of the out-of-plane external field ( 0–130 mT) and values of the DMI for which a Néel skyrmion was stabilized ( 0.9–1.65 mJ m, to be compared to the threshold DMI strength without field 1.52 mJ m here). The current density is set at 3 A m.

In Figs. 4a (bottom graph) and 4b we display the expected velocities in the free and confined cases, respectively, as a function of the skyrmion diameter . We find that as long as the Néel skyrmion configuration is maintained, the velocity of the skyrmion is determined by its size. Two regimes can be observed on the graphs of Fig. 4a, and they can be distinguished by the ratio between the skyrmion diameter and the skyrmion wall width , which we display as a function of in the top graph of Fig. 4a. We define the skyrmion wall width as a usual domain wall width would be defined from the 1D magnetization profile in ), where is the distance from the center of the wall/skyrmion. asymptotically reaches the natural domain wall width value ( is the effective anisotropy also including the out-of plane field) for the large , as expected in a large skyrmionic bubble, which is shown by the solid line in Fig. 4a. On the contrary, for the small in the order of , the skyrmion is compact with its wall width compressed to , following the dashed line in Fig. 4a, hence the change of behavior for the smaller skyrmions.

For the sizes 50 nm, the skyrmion remains compact, that is, the extension of the part of reversed magnetization with pointing down is small as compared to the area covered by a tilted magnetization. For such compact skyrmions, both longitudinal and transverse speed increase monotonically with the diameter . For the diameters 50 nm, an area of uniformly reversed magnetization progressively appears and turns the skyrmion into a skyrmionic bubble for 100 nm (on the cut of the skyrmion along its diameter, the magnetization angle no longer rotates continuously but exhibits a plateau in the center, see Fig. 4b in Ref. 10, for example). The transition to skyrmionic bubbles significantly alters the dynamics. For increasing above 50 nm, the transverse velocity saturates and then decreases after 100 nm, while the longitudinal velocity in the longitudinal direction still increases but at a slower rate. A consequence of these results is that, in the sub-100 nm regime, as the force exerted by the SOTs on the skyrmion is proportional to the domain wall width Sampaio et al. (2013), smaller skyrmions thus move slower Ding et al. (2015). In the confined case, a similar transition between two regimes is observed in Fig. 4b, where the skyrmion velocity first increases almost linearly with the diameter as long as the skyrmion is compact, and then saturates. Note the significantly higher speed as compared to the free case, due to the effect of the confining potential Iwasaki et al. (2014), which allows a 5 to 15 times increase in the longitudinal velocity of the compact skyrmions ( 50 nm).

As explained above, the SOTs in the repeated Pt layers are expected to be much smaller than in the Ta layer, if not negligible. As a consequence, the torque originating in the bottom Ta layer exerted at on one single layer is diluted and become similar to a torque exerted on each of the repeated 11 Co layers at . As the velocities are proportional to the current density Iwasaki et al. (2013b); Sampaio et al. (2013), we thus need to rescale the velocities obtained above at 3 A m by a factor of up to 11 to compare with the presented experimental conditions. At the current density of 2.85 A m, for a diameter estimated Moreau-Luchaire et al. (2016) as 80 nm, the predicted velocities are thus around 1.4 m s and 2 m s if we neglect completely the torques from the repeated Pt layers. These size-dependance considerations, together with the less efficient spin torque in the structures, thus provide a first explanation why our experimental velocities for sub-100 nm skyrmions are found to be smaller that the ones observed by other groups for larger skyrmions Woo et al. (2016); Hrabec et al. () ( 200–350 nm). The size dependence could be thought as a disadvantage for applications, however, in the perspective of decreasing the device sizes, the distance a skyrmion has to travel in a device scales linearly with the size of the device/skyrmion. As a consequence, the almost linear dependence of the velocity as a function of the size of a compact skyrmion, as can be observed in Fig. 4a, does not affect the potential operation times in devices.

The Co-based multilayers with PMA grown by sputtering are known to exhibit local variations of the magnetic parameters due to the formation of material grains, which could also lead to pinning and apparently disordered motion Lemerle et al. (1998). The present samples have been characterized by transmission electron microscopy, which reveals the presence of grains of a characteristic size of 20 nm (typical for the present conditions of growth in our sputtering system). Moreover, the magnetic inhomogeneities in asymmetric magnetic multilayers similar to the present stacks have been characterized by Baćani et al. using calibrated MFM tips Baćani et al. (). They have found a mean value of the DMI very close to the one extracted from the field-dependance of the skyrmion size in our multilayers Moreau-Luchaire et al. (2016), but with strong local variations. In the images of Figs. 3a-d, there is a clear dispersion in the size of the skyrmions, with a standard deviation on the diameters 163 nm for a mean diameter 176 nm. Note that these sizes are apparent sizes after treatment of the MFM images for the above skyrmion tracking. The evaluation of the actual size of the skyrmion analysing the skyrmion cuts and taking into account the resolution of the tip rather suggests 80–100 nm, as in our previous work Moreau-Luchaire et al. (2016). In addition to that, some skyrmions on the same images have a non-ideal circular shape. These observations support that the magnetic properties in the system are not perfectly uniform within the track.

In order to achieve a more realistic modeling of the motion of skyrmions in our nanotracks, we have included such magnetic inhomogeneities in our micromagnetic simulation. We model a granular system in which the distribution of both positions and sizes of the grains, as well as the DMI value in these grains, are random. We have considered here only local variations of the DMI, but we note that variations of the other magnetic parameters, e. g. the magnetic anisotropy Kim and Yoo () or the exchange stiffness, would also result in alterations of the dynamics of the skyrmions. In Fig. 5a, we present an example of granular distribution used in these simulations, with a mean grain size 30 nm and a normal distribution of around 1.8 mJ m (value taken from our previous work Moreau-Luchaire et al. (2016) for a similar multilayer) with a standard deviation = 10%. For ultrathin layers of a thickness of a couple atoms, local variations of the thickness by one atom can dramatically change thickness-dependent anisotropy, DMI, etc. by up to 75% Baćani et al. (), so that our modeling here is realistic assumption (see the distribution ranging 55% of in Fig. 5a). For each grain size, the magnetic configuration is relaxed from a random magnetization distribution to obtain a population of skyrmions such as the one displayed in Fig. 5b (thus mimicking the mechanism of relaxation after Joule heating reinitialization of the magnetic order, similar to what has been shown before with MFM images) under a bias field 200 mT. The skyrmions remain compact and the inhomogenities of lead to a dispersion of their diameters, as can be observed in Fig. 5b. Then current-induced motion is simulated as above, for different current densities ranging 0.75–9 A m.

Figure 5: a. Map of the DMI in the simulation grid of 1024512 nm, where the values are indicated by the color scale on the right, for 30 nm. b. Initial state of the simulation after relaxing the skyrmion state, for 30 nm, showing on a scale ranging -1 (black) to 1 (white). Trajectories of around 80 different skyrmions simulated for c. no grains and 7.5 A m, d. no grains and 3.75 A m, e. grains with 30 nm and 7.5 A m, f. grains with 30 nm and 3.75 A m.

In Figs. 5c-f, we display the trajectories of all the skyrmions in the simulated nanotrack (1024512 nm, with x and y periodic boundary conditions) at two different current densities ( 7.5 A m and 3.75 A m), in the uniform case (no grains) and in the presence of inhomogeneities ( 30 nm). Comparing the length and angle of the trajectories, it appears that the local variations of cause an alteration of both the longitunal and transverse motion. As can be seen in Fig. 5e, the non-uniform leads to pinned skyrmions and reduces the total distance travelled by the skyrmions. Moreover, as can be seen in Figs. 5e,f, the non-uniform also causes the skyrmion Hall angle to drastically decrease, with the stronger effects at low current densities, where part of the skyrmions are even driven in a transverse direction opposite to the one of the uniform case. Such pinning behavior and perturbed trajectories are influenced by both grain size and current density. In Fig. 6a we summarize the simulated velocity (with combined x and y components) of the skyrmions for different mean grain sizes ( 15–60 nm and no grains) and for the whole range of simulated current densities. The velocity curves exhibit a clear pinning-like behavior, with the mean velocities at low currents being shifted towards zero due to the pinned skyrmions, consequently departing from the expected linear current-velocity relation of the uniform case. At high current densities however, above 6 A m, the skyrmion mobility (defined as ) is recovered and is comparable for all simulated grain sizes. Nevertheless, even at very high current densities (as shown after the break in the figure), the velocity is not totally recovered, different from the case of impurity potentials Iwasaki et al. (2013a). On each data point, the bar describes the standard deviation of the velocities for the 80 skyrmions simulated in the system. At low currents, in particular, these bars are highly asymmetric, meaning that many skyrmions are pinned and do not move at all, while a few of them only are able to move and to acquire a non-zero velocity. This behavior reproduces strikingly well the observations of the experiments presented above: at the corresponding current density indicated by a star in Fig. 6a, the velocity is expected to dramatically reduce by up to a factor of 5 due to the inhomogeneities of DMI, and by even more at lower currents, where motion is almost totally impeded due to pinning. In Fig. 6b we summarize the angle of the motion for the different mean grain sizes and current densities. At high current densities, the trajectories have nearly the same angle as in the uniform case (and actually converge to the same angle at very high currents, as shown after the break in the figure). However, at low current densities, the magnetic disorder affects the gyrotropic motion of the skyrmions and causes their trajectories to align with the current direction. At current densities comparable with our experiments (marked by a star in the graph), no significant angle between current and motion is expected. This effective cancellation of the skyrmion Hall effect is in perfect agreement with what we observe experimentally. These results showing the reduction of the skyrmion Hall effect are to be compared to those recently published by Reichhardt et al., which have simulated the presence of punctual defects in a skyrmionic system and obtained a comparable decrease of the apparent skyrmion Hall angle at low currents Reichhardt and Olson Reichhardt (2016).

Figure 6: a. Mean velocity and b. mean trajectory angle of the skymions vs. current density for different grain sizes 60 nm (blue hexagons), 30 nm (gold pentagons) and 15 nm (red diamonds). The uniform case (black squares) is fitted to a line showing in a. the proportionality of to and in b. the constant angle in this case. After the break in a., the scale is identical but the values are offset as indicated on the right part of the graph. The stars indicate the experimental current density. c. Mean velocity and d. mean trajectory angle of the skymions vs. current density for different grain sizes 15 nm (red diamonds), 10 nm (blue triangles) and 5 nm (green stars). The bars on top of each data point show the standard deviation for the around 80 skyrmions considered in each simulation. In a. and c. the colored lines are fits of the velocity for each grain size to a line of mobility identical to the mobility without grains. This allows us to extract the critical depinning current being the x-intercept, as indicated by the colored arrows. e. Critical depinning current as a function of , the skyrmion diameter relative to the grain size. The mean skyrmion diameter is here 24 nm and the curve is a log-gaussian fit showing the relevance of the relative size criterion.

We have further reduced the mean grain size in the simulations, down to 5 nm, in order to study the case where the grains are much smaller than the skyrmion diameter. As can be seen in Figs. 6c,d, for such a small scale of the inhomogeneities, the ideal behavior of steady motion is progressively recovered with decreasing grain size, because the variations in the magnetic energies affecting the skyrmions are averaged over their size. The influence of the grain size can be better seen in Fig. 6e which summarizes the extracted effective depinning current as a function of the diameter of the skyrmion , relative to the mean grain size . To find the critical depinning current density , we use linear fits to the velocity data points above 6 A m and extrapolate the current at null velocity (x-intercept, indicated by colored arrows in Figs. 6a,c). In this regime the skyrmion mobility is fully recovered, so that the mobility obtained from the uniform case can be used for all the fits. The error bars in Fig. 6e indicate the errors on the fit. The pinning is found to be strongest when the grain size matches the diameter of the skyrmions. As shown by the gaussian-log fit in Fig. 6e, the depinning current density can be roughly evaluated as . is the maximum depinning current density, which is here 2.7 A m, very close to our experimental current densities, while A is a width parameter. Because the Co-based multilayers studied in this work have a granular structure with grains at the 20 nm scale, the magnetic inhomogeneities are thus expected to largely influence the motion of sub-100 nm skyrmions. However, for larger skyrmions with diameters ranging 100 nm–1 m, the magnetic inhomogeneities are averaged over the size of the skyrmion, which results in a motion closer to the flow regime in an ideal, uniform material. However, it should be noted that significant pinning is still found in the experiments of other groups Woo et al. (2016); Hrabec et al. (). Magnetic inhomogeneities thus constitute another reason why the experimental velocities that we find for sub-100 nm skyrmions are much smaller than the ones observed for larger skyrmions, especially when other materials with possibly less pinning are used. In order to overcome the pinning to manipulate sub-100 nm skyrmions, better design of the material structure and control of material growth is necessary. Reducing the size of the skyrmions far below the grains size is another possibility to recover a non-disturbed motion.

In conclusion, we have experimentally studied the generation of sub-100 nm skymions in magnetic nanotracks observed by MFM, as well as their motion under spin-torques. Notably, the confinement to small widths of the tracks allows us to investigate a novel approach for the generation of skyrmions, under homogeneous current injection. Once the compact sub-100 nm skyrmions are obtained and stabilized directly in the track, the use of shorter and stronger pulses allows us to observe the current induced motion of the same skyrmions due to the spin-orbit torques arising from a bottom Ta layer. This way, the same track of simple design can be used for both nucleation and motion of the skyrmions. However we have found that the motion of small skyrmions in the sub-100 nm regime is largely influenced by the inhomogeneities of the magnetic Co layers. The disturbed and apparently disordered motion of the skyrmions is well reproduced by simulating a track made of grains of varying magnetic parameters, here a varying Dzyaloshinskii-Moriya interaction strength being taken as an example. Such micromagnetic simulations show the pinning behaviour of the skyrmion motion, leading to a depinning current which reaches a maximum for sizes of grains and skyrmions matching together, as well as the cancellation of the skyrmion Hall effect at low currents. These results highlight that creating skyrmionic structures in nanotracks is getting easier, but also highlight the importance of controlling the inhomogeneities of the magnetic properties in order to design skyrmion devices achieving a high degree of control even at low currents.

The authors acknowledge Cyrile Deranlot for his assistance and suggestions in the sample preparation as well as all the partners involved in MAGicSky consortium for fruitful discussions. Financial support from European Union grant MAGicSky No. FET-Open-665095. 103 is acknowledged.


  1. W. Jiang, P. Upadhyaya, W. Zhang, G. Yu, M. B. Jungfleisch, F. Y. Fradin, J. E. Pearson, Y. Tserkovnyak, K. L. Wang, O. Heinonen, S. G. E. te Velthuis,  and A. Hoffmann, Science 349, 283 (2015).
  2. G. Chen, A. Mascaraque, A. T. N’Diaye,  and A. K. Schmid, Appl. Phys. Lett. 106, 242404 (2015).
  3. C. Moreau-Luchaire, C. Moutafis, N. Reyren, J. Sampaio, C. A. F. Vaz, N. Van Horne, K. Bouzehouane, K. Garcia, C. Deranlot, P. Warnicke, P. Wohlhüter, J.-M. George, M. Weigand, J. Raabe, V. Cros,  and A. Fert, Nat. Nanotechnol. 11, 444 (2016).
  4. O. Boulle, J. Vogel, H. Yang, S. Pizzini, D. de Souza Chaves, A. Locatelli, T. O. Mentes, A. Sala, L. D. Buda-Prejbeanu, O. Klein, M. Belmeguenai, Y. Roussigné, A. Stashkevich, S.-M. Chérif, L. Aballe, M. Foerster, M. Chshiev, S. Auffret, I. M. Miron,  and G. Gaudin, Nat. Nanotechnol. 11, 449 (2016).
  5. G. Yu, P. Upadhyaya, X. Li, W. Li, S. K. Kim, Y. Fan, K. L. Wong, Y. Tserkovnyak, P. K. Amiri,  and K. L. Wang, Nano Lett. 16, 1981 (2016).
  6. S. Woo, K. Litzius, B. Krüger, M.-Y. Im, L. Caretta, K. Richter, M. Mann, A. Krone, R. M. Reeve, M. Weigand, P. Agrawal, I. Lemesh, M.-A. Mawass, P. Fischer, M. Kläui,  and G. S. D. Beach, Nat. Mater. 15, 501 (2016).
  7. A. N. Bogdanov and U. K. Rößler, Phys. Rev. Lett. 87, 037203 (2001).
  8. N. Nagaosa and Y. Tokura, Nat. Nanotechnol. 8, 899 (2013).
  9. A. Fert, V. Cros,  and J. Sampaio, Nat. Nanotechnol. 8, 152 (2013).
  10. N. S. Kiselev, A. N. Bogdanov, R. Schäfer,  and U. K. Rößler, J. Phys. D: Appl. Phys. 44, 392001 (2011).
  11. J. Sampaio, V. Cros, S. Rohart, A. Thiaville,  and A. Fert, Nat. Nanotechnol. 8, 839 (2013).
  12. N. Romming, C. Hanneken, M. Menzel, J. E. Bickel, B. Wolter, K. von Bergmann, A. Kubetzka,  and R. Wiesendanger, Science 341, 636 (2013).
  13. N. Romming, A. Kubetzka, C. Hanneken, K. von Bergmann,  and R. Wiesendanger, Phys. Rev. Lett. 114, 177203 (2015).
  14. W. Jiang, W. Zhang, G. Yu, M. B. Jungfleisch, P. Upadhyaya, H. Somaily, J. E. Pearson, Y. Tserkovnyak, K. L. Wang, O. Heinonen, S. G. E. te Velthuis,  and A. Hoffmann, AIP Adv. 6, 055602 (2016).
  15. W. Jiang, X. Zhang, G. Yu, W. Zhang, X. Wang, M. Benjamin Jungfleisch, J. E. Pearson, X. Cheng, O. Heinonen, K. L. Wang, Y. Zhou, A. Hoffmann,  and S. G. E. te Velthuis, Nat. Phys. 13, 162 (2017).
  16. K. Litzius, I. Lemesh, B. Krüger, P. Bassirian, L. Caretta, K. Richter, F. Büttner, K. Sato, O. A. Tretiakov, J. Förster, R. M. Reeve, M. Weigand, I. Bykova, H. Stoll, G. Schütz, G. S. D. Beach,  and M. Kläui, Nat. Phys. 13, 170 (2017).
  17. A. Hrabec, J. Sampaio, M. Belmeguenai, I. Gross, R. Weil, S. M. Chérif, A. Stachkevitch, V. Jacques, A. Thiaville,  and S. Rohart, arXiv:1611.00647 .
  18. A. Fert, Metallic MultilayersMater. Sci. Forum 59–60, 439 (1990).
  19. B. Dupé, G. Bihlmayer, M. Böttcher, S. Blügel,  and S. Heinze, Nat. Commun. 7, 11779 (2016).
  20. H. Yang, A. Thiaville, S. Rohart, A. Fert,  and M. Chshiev, Phys. Rev. Lett. 115, 267210 (2015).
  21. A. Soumyanarayanan, M. Raju, A. Oyarce, A. K. Tan, M.-Y. Im, A. Petrovic, P. Ho, K. Khoo, M. Tran, C. Gan, F. Ernult,  and C. Panagopoulos, arXiv:1606.06034v2 .
  22. M. Bode, M. Heide, K. von Bergmann, P. Ferriani, S. Heinze, G. Bihlmayer, A. Kubetzka, O. Pietzsch, S. Blügel,  and R. Wiesendanger, Nature 447, 190 (2007).
  23. S. Heinze, K. von Bergmann, M. Menzel, J. Brede, A. Kubetzka, R. Wiesendanger, G. Bihlmayer,  and B. Stefan, Nat. Phys. 7, 713 (2011).
  24. A. Thiaville, S. Rohart, E. Jué, V. Cros,  and A. Fert, Europhys. Lett. 100, 57002 (2012).
  25. S. Emori, U. Bauer, S.-M. Ahn, E. Martinez,  and G. S. D. Beach, Nat. Mater. 12, 611 (2013).
  26. A. Neubauer, C. Pfleiderer, B. Binz, A. Rosch, R. Ritz, P. G. Niklowitz,  and P. Böni, Phys. Rev. Lett. 102, 186602 (2009).
  27. T. Schulz, R. Ritz, A. Bauer, M. Halder, M. Wagner, C. Franz, C. Pfleiderer, K. Everschor, M. Garst,  and A. Rosch, Nat. Phys. 8, 301 (2012).
  28. N. Kanazawa, M. Kubota, A. Tsukazaki, Y. Kozuka, K. S. Takahashi, M. Kawasaki, M. Ichikawa, F. Kagawa,  and Y. Tokura, Phys. Rev. B 91, 041122 (2015).
  29. G. Yin, Y. Liu, Y. Barlas, J. Zang,  and R. K. Lake, Phys. Rev. B 92, 024411 (2015).
  30. J. Iwasaki, M. Mochizuki,  and N. Nagaosa, Nat. Nanotechnol. 8, 742 (2013a).
  31. F. Ma, Y. Zhou, H. B. Braun,  and W. S. Lew, Nano Lett. 15, 4029 (2015).
  32. S. Rohart and A. Thiaville, Phys. Rev. B 88, 184422 (2013).
  33. O. Heinonen, W. Jiang, H. Somaily, S. G. E. te Velthuis,  and A. Hoffmann, Phys. Rev. B 93, 094407 (2016).
  34. J.-C. Rojas-Sánchez, N. Reyren, P. Laczkowski, W. Savero, J.-P. Attané, C. Deranlot, M. Jamet, J.-M. George, L. Vila,  and H. Jaffrès, Phys. Rev. Lett. 112, 106602 (2014).
  35. A. Vansteenkiste, J. Leliaert, M. Dvornik, M. Helsen, F. Garcia-Sanchez,  and B. Van Waeyenberge, AIP Adv. 4, 107133 (2014).
  36. L. Liu, C.-F. Pai, Y. Li, H. W. Tseng, D. C. Ralph,  and R. A. Buhrman, Science 336, 555 (2012).
  37. X. Qiu, P. Deorani, K. Narayanapillai, K.-S. Lee, K.-J. Lee, H.-W. Lee,  and H. Yang, Sci. Rep. 4, 4491 (2014).
  38. J. Iwasaki, W. Koshibae,  and N. Nagaosa, Nano Lett. 14, 4432 (2014).
  39. S. Rohart, J. Miltat,  and A. Thiaville, Phys. Rev. B 93, 214412 (2016).
  40. J. Ding, X. Yang,  and T. Zhu, J. Phys. D: Appl. Phys. 48, 115004 (2015).
  41. J. Iwasaki, M. Mochizuki,  and N. Nagaosa, Nat. Commun. 4, 1463 (2013b).
  42. S. Lemerle, J. Ferré, C. Chappert, V. Mathet, T. Giamarchi,  and P. Le Doussal, Phys. Rev. Lett. 80, 849 (1998).
  43. M. Baćani, M. A. Marioni, J. Schwenk,  and H. J. Hug, arXiv:1609.01615 .
  44. J.-V. Kim and M.-W. Yoo, private communication .
  45. C. Reichhardt and C. J. Olson Reichhardt, New J. Phys. 18, 095005 (2016).
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Comments 0
Request comment
The feedback must be of minumum 40 characters
Add comment
Loading ...

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description