Skyrmions in magnetic systems are particle-like objects predicted to occur in materials with chiral interactions [1]. The existence of a hexagonal skyrmion lattice in chiral magnets was subsequently confirmed in neutron scattering experiments [2] and in direct imaging experiments [3]. Since then, skyrmion states have been found in an increasing number of compounds [4–8], including materials in which skyrmions are stable at room temperature [9–14]. Skyrmions can be set into motion by applying an external current [15, 16], and effective skyrmion velocity versus driving force curves can be calculated from changes in the Hall resistance [17, 18] or by direct imaging of the skyrmion motion [9, 14]. Additionally, transport curves can be studied numerically with continuum and particle based models [19–23]. Both experiments and simulations show that there is a finite depinning threshold for skyrmion motion similar to that found for the depinning of current-driven vortex lattices in type-II superconductors [24–26]. Since skyrmions have particle like properties and can be moved with very low driving currents, they are promising candidates for spintronic applications [27, 28], so an understanding of skyrmion motion and depinning is of paramount importance. Additionally, skyrmions represent an interesting dynamical system to study due to the strong non-dissipative effect of the Magnus force they experience, which is generally very weak or absent altogether in other systems where depinning and sliding phenomena occur.
For particle-based representations of the motion of objects such as superconducting vortices, a damping term of strength aligns the particle velocity in the direction of the net force acting on the particle, while a Magnus term of strength
rotates the velocity component in the direction perpendicular to the net force. In most systems studied to date, the Magnus term is very weak compared to the damping term, but in skyrmion systems the ratio of the Magnus and damping terms can be as large as
[17, 19, 21, 29]. One consequence of the dominance of the Magnus term is that under an external driving force, skyrmions develop velocity components both parallel (
) and perpendicular (V⊥) to the external drive, producing a skyrmion Hall angle of
, where
. In a completely pin-free system, the intrinsic skyrmion Hall angle has a constant value
however, in the presence of pinning a moving skyrmion exhibits a side jump phenomenon in the direction of the drive so that the measured Hall angle is smaller than the clean value [22, 23, 30]. In studies of these side jumps using both continuum and particle based models for a skyrmion interacting with a single pinning site [22] and a periodic array of pinning sites [30], R increases with increasing external drive until the skyrmions are moving fast enough that the pinning becomes ineffective and the side jump effect is minimized.
Particle-based studies of skyrmions with an intrinsic Hall angle of moving through random pinning arrays show that
at small drives and that
increases with increasing drive until saturating at
for high drives [23]. In recent imaging experiments performed in the single skyrmion limit [31] it was shown that R=0 and
at depinning and that both quantities increase linearly with increasing drive; however, the range of accessible driving forces was too low to permit observation of a saturation effect. These experiments were performed in a regime of relatively strong pinning, where upper limits of
and
are expected. A natural question is how universal the linear behavior of R and
is as a function of drive, and whether the results remain robust for larger intrinsic values of
. It is also interesting to ask what happens in the weak pinning limit where the skyrmions form a hexagonal lattice and depin elastically. In studies of overdamped systems such as superconducting vortices, it is known that the strong and weak pinning limits are separated by a transition from elastic to plastic depinning and have very different transport curve characteristics [24, 26], so a similar phenomenon could arise in the skyrmion Hall effect. Noise fluctuations have also been used as another method to study the dynamics of magnetic systems [32]. In superconducting vortex systems, the plastic flow regime is associated with large voltage noise fluctuations of
form [33–36], while when the system dynamically orders at higher drives, narrow band noise features appear and the noise power is strongly reduced [26, 37–39]. Here we show that changes in the skyrmion Hall angle are correlated with changes in the skyrmion velocity fluctuations and the shape of the velocity noise spectrum. In the plastic flow region where R increases linearly with drive, there is a
velocity noise signal with
, while when R reaches its saturation value, there is a crossover to white noise or weak narrow band noise, indicating that noise measures could provide another way to probe skyrmion dynamics. In general, we find that the narrow band noise features are much weaker in the skyrmion case than in the superconducting vortex case due to the Magnus effect.
Simulation and system— We consider a 2D simulation with periodic boundary conditions in the x and y-directions using a particle-based model of a modified Thiele equation recently developed for skyrmions interacting with random [21, 23] and periodic [30, 40] pinning substrates. The simulated region contains N skyrmions, and the time evolution of a single skyrmion i at position is governed by the following equation:
Here, the skyrmion velocity is ,
is the damping term, and
is the Magnus term. We impose the condition
to maintain a constant magnitude of the skyrmion velocity for varied
. The repulsive skyrmion–skyrmion interaction force is given by
where
,
, and K1 is a modified Bessel function that falls off exponentially for large rij. The pinning force
arises from non-overlapping randomly placed pinning sites modeled as harmonic traps with an amplitude of Fp and a radius of
as used in previous studies [23]. The driving force
is from an applied current interacting with the emergent magnetic flux carried by the skyrmion [17, 29]. We increase
slowly to avoid transient effects. In order to match the experiments, we take the driving force to be in the positive x-direction so that the Hall effect is in the negative y-direction. We measure the average skyrmion velocity
(
) in the direction parallel (perpendicular) to the applied drive, and we characterize the Hall effect by measuring
for varied FD. The skyrmion Hall angle is
. We consider a system of size L=36 with a fixed skyrmion density of
and pinning densities ranging from np=0.006 25 to np=0.2.
In figures 1(a) and (b) we plot ,
, and R versus FD for a system with
,
, and
. In this regime, plastic depinning occurs, meaning that at the depinning threshold some skyrmions can be temporarily trapped at pinning sites while other skyrmions move around them. The velocity–force curves are nonlinear, and
increases more rapidly with increasing FD than
. The inset of figure 1(a) shows that
for
, indicating that just above the depinning transition the skyrmions are moving predominantly in the direction of the driving force. In figure 1(b), R increases linearly with increasing FD for
, as indicated by the linear fit, while for
R saturates to the intrinsic value of R=5.708 marked with a dashed line. The inset of figure 1(b) shows the corresponding
versus FD. From an initial value of 0°,
increases with increasing FD before saturating at the clean limit value of
. Although the linear increase in R with FD is similar to the behavior observed in the experiments of [31],
does not show the same linear behavior as in the experiments; however, we show later that when the intrinsic skyrmion Hall angle is small,
varies linearly with drive. We note that the experiments in [31] were performed in the single skyrmion limit rather than in the many skyrmion plastic flow limit we study. This could impact the behavior of the Hall angle, making it difficult to directly compare our results with these experiments.
In figure 2 we illustrate the skyrmion positions and trajectories obtained during a fixed period of time at different drives for the system in figure 1. At in figure 2(a), R=0.15 and the average drift is predominantly along the x-direction parallel to the drive, taking the form of riverlike channels along which individual skyrmions intermittently switch between pinned and moving states. In figure 2(b), for
we find R=0.6, and observe wider channels that begin to tilt along the negative y-direction. At
in figure 2(c), R=1.64 and
. The skyrmion trajectories are more strongly tilted along the
direction, and there are still regions of temporarily pinned skyrmions coexisting with moving skyrmions. As the drive increases, individual skyrmions spend less time in the pinned state. Figure 2(d) shows a snapshot of the trajectories over a shorter time scale at
where R=5.59. Here the plastic motion is lost and the skyrmions form a moving crystal translating at an angle of
with respect to the external driving direction, which is close to the clean value limit of
. In general, the deviations from linear behavior that appear as R reaches its saturation value in figure 1(b) coincide with the loss of coexisting pinned and moving skyrmions, and are thus correlated with the end of plastic flow.
In figure 3(a) we show R versus FD for the system from figure 1 at varied . In all cases, between the depinning transition and the free flowing phase there is a plastic flow phase in which R increases linearly with FD with a slope that increases with increasing
. In contrast to the nonlinear dependence of
on FD at
illustrated in the inset of figure 1(b), figure 3(b) shows that for
,
increases linearly with FD and
. To understand the linear behavior, consider the expansion of
For small
, as in the experiments,
, and since R increases linearly with FD,
also increases linearly with FD. In general, for
we find an extended region over which
grows linearly with FD, while for
, the dependence of
on FD has nonlinear features similar to those shown in the inset of figure 1(b). In figure 3(c) we plot R versus FD for a system with
for varied Fp. In all cases R increases linearly with FD before saturating; however, for increasing Fp, the slope of R decreases while the saturation of R shifts to higher values of FD. In general, the linear behavior in R is present whenever Fp is strong enough to produce plastic flow. In figure 3(d) we show R versus FD at
for varied pinning densities np. In each case, there is a region in which R increases linearly with FD, with a slope that increases with increasing np. As np becomes small, the nonlinear region just above depinning where R increases very rapidly with drive becomes more prominent.
For weak pinning, the skyrmions form a triangular lattice and exhibit elastic depinning, in which each skyrmion maintains the same neighbors over time. In figure 4(a) we plot the critical depinning force Fc and the fraction P6 of sixfold-coordinated skyrmions versus Fp for a system with and
. For
, the skyrmions depin elastically. In this regime,
and Fc increases as
as expected for the collective depinning of elastic lattices [25]. For
, P6 drops due to the appearance of topological defects in the lattice, and the system depins plastically, with
as expected for single particle depinning or plastic flow.
In figure 4(b) we plot R versus FD in samples with and np=0.1 in the elastic depinning regime for varied
. We highlight the nonlinear behavior for the
case by a fit of the form
with
and
. The dotted line indicates the corresponding clean limit value of R=5.708. We find that R is always nonlinear within the elastic flow regime, but that there is no universal value of β, which ranges from
to
with varying
. The change in the Hall angle with drive is most pronounced just above the depinning threshold, as indicated by the rapid change in R at small FD. This results from the elastic stiffness of the skyrmion lattice which prevents individual skyrmions from occupying the most favorable substrate locations. In contrast, R changes more slowly at small FD in the plastic flow regime, where the softer skyrmion lattice can adapt to the disordered pinning sites. In figure 4(c) we plot R versus FD at
and np=0.1 for varied Fp, showing a reduction in R with increasing Fp. A fit of the
curve in the plastic depinning regime shows a linear increase of R with FD, while for
in the elastic regime, the dependence of R on FD is nonlinear. Just above depinning in the elastic regime, the skyrmion flow direction rotates with increasing drive.
Correlations between noise fluctuations and the skyrmion Hall effect.
The power spectrum of the velocity noise fluctuations at different applied drives represents another method that can be used to probe the dynamics of driven condensed matter systems. In the superconducting vortex case, the total noise power over a particular frequency range or the overall shape of the noise power spectrum can be determined by measuring the voltage time series at a particular current. Both experiments and simulations have shown that in the plastic flow regime, where the vortex flow is disordered and consists of a combination of pinned and flowing particles, the low frequency noise power is large and the voltage noise spectrum has a character with
. In contrast, the low frequency noise power is considerably reduced in the elastic or ordered flow regime, where the noise is either white with
or exhibits a characteristic washboard frequency associated with narrow band noise. Based on these changes in the noise characteristics, it is possible to map out a dynamical phase diagram for the vortex system.
In the skyrmion system, the Hall resistance can be used to detect the motion of skyrmions, and therefore, in analogy with the voltage response in a superconductor, fluctuations in the Hall resistance at a specific applied current should reflect fluctuations in the skyrmion velocity. Since the recent experiments of [31] used imaging techniques to measure R, it is desirable to understand whether changes in R are correlated with changes in the fluctuations of other quantities. We measure the time series of and
in samples with
,
, and
, the same parameters used in figures 1 and 2. For these values, R increases linearly with increasing FD over the range
before saturating close to the clean value. For each value of FD, we then construct the power spectrum
where FD is held constant during an interval of simulation time steps.
In figure 5(a) we plot for
and
at
in the plastic flow regime. Here, the spectral shape is very similar in each case, while the noise power at low frequencies is slightly higher for
than for V⊥. At FD=0.4 in figure 5(b), deep in the plastic flow phase, the noise power is much higher for V⊥ than for
and can be fit reasonably well to a
form, while
also has an
shape over a less extended region. The two spectra have equal power only for high ω. Figure 5(c) shows
at
, which corresponds to the saturation region of R. The V⊥ signal still has the highest spectral power, but both spectra now exhibit a white or
shape. There is a small bump at low frequency which is more prominent in
that may correspond to a narrow band noise feature. We note that in the overdamped limit of
at this same drive, where the particles have formed a moving lattice, there is a strong narrow band noise feature, suggesting that the Magnus term is responsible for the lack of a strong narrow band noise peak in figure 5(c). In general, we find that the power spectrum for the skyrmions shows
noise in the plastic flow regime and white noise in the saturation regime.
Using the power spectrum, we can calculate the noise power S0 at a specific value of ω. In figure 6 we plot for
and V⊥ versus FD, along with the corresponding R curve from figure 1. At low FD, the value of S0 is nearly the same for both
and V⊥. The noise power for V⊥ increases more rapidly with increasing FD and both S0 curves reach a maximum near
before decreasing as R reaches its saturation value. In general S0 is large whenever the spectrum has a
shape. This result shows that noise power fluctuations could be used to probe changes in the skyrmion Hall effect and even dynamical transitions from plastic to elastic skyrmion flow.
We note that in real skyrmion systems, the skyrmions can also have internal modes of motion that could affect the noise power. Such internal modes are not captured by the particle model, and would likely occur at much higher frequencies than those of the skyrmion center of of mass motion that we analyze here. It would be interesting to see if such modes arise in experiment and to determine whether they can also modify the skyrmion Hall angle.
We have investigated the skyrmion Hall effect by measuring the ratio R of the skyrmion velocity perpendicular and parallel to an applied driving force. In the disorder-free limit, R and the skyrmion Hall angle take constant values independent of the applied drive; however, in the presence of pinning these quantities become drive-dependent, and in the strong pinning regime R increases linearly from zero with increasing drive, in agreement with recent experiments. For large intrinsic Hall angles, the current-dependent Hall angle increases nonlinearly with increasing drive; however, for small intrinsic Hall angles such as in recent experiments, both the current-dependent Hall angle and R increase linearly with drive as found experimentally. The linear dependence of R on drive is robust for a wide range of intrinsic Hall angle values, pinning strengths, and pinning densities, and appears whenever the system exhibits plastic flow. For weaker pinning where the skyrmions depin elastically, R has a nonlinear drive dependence and increases very rapidly just above depinning. We observe a crossover from nonlinear to linear drive dependence of R as a function of the pinning strength, which coincides with the transition from elastic to plastic depinning. We also show how R correlates with changes in the power spectra of the velocity noise fluctuations both parallel and perpendicular to the drive. In the plastic flow regime where R increases linearly with increasing FD, we find noise that crosses over to white noise at higher drives. The noise power drops dramatically as R saturates at high drives.
We gratefully acknowledge the support of the US Department of Energy through the LANL/LDRD program for this work. This work was carried out under the auspices of the NNSA of the US DoE at LANL under Contract No. DE-AC52-06NA25396.