3. Results
The function between skyrmion velocity and the position in the circular ring nanotrack enables us to determine the dynamics of skyrmion. We calculate the skyrmion velocity and decompose it into two components, parallel velocity
and perpendicular velocity
, where
(
) is the velocity component along(perpendicular) the current direction. The positive directions of
and
are along
and
, respectively. It can be seen from
Figure 1(b) that
shows a periodic trend of increasing first and then decreasing in both CCW and CW directions, where
of CCW (CW) with
Q=-1 overlaps with CW(CCW) with
Q=1. This fluctuation also occurs in
, as shown in
Figure 1(c), which fluctuates over a non-zero value and has the same periodicity as
. This fluctuation of velocity is significantly different from the case in a straight track. In a straight track, the velocity remains constant in parallel and zero in perpendicular direction for
. However, in a circular ring nanotrack, the velocity fluctuates significantly. To compare the movement behavior of a skyrmion in both tracks, we define
Rsk as the radius of the real moving trajectory and
Rm as the radius of the middle line of the circular ring
. Here we let
Rd=
Rm-
Rsk, which is the distance between the middle line and the real moving trajectory.
Figure 1(d) shows the variation of
Rd with the angle position. The initial non-zero value of
Rd indicates that the skyrmion undergoes a transverse drift from the middle line. It means that the skyrmion initially stabilizes inside the middle line and then moves towards the inner or outer boundary of the track driven by the current flow.
Upon observing skyrmion movement in the circular ring nanotrack, we found that the skyrmion moves at a certain angle relative to the driving current, similar to the drift caused by the skyrmion Hall effect (SkHE). We so calculate the Hall angle
of the skyrmion in a circular ring track and show it in
Figure 2(a), where
. It is evident that the circular ring nanotrack exhibits a non-zero Hall angle, with a maximum value of about 2°. We discovered that this type of SkHE still exists when the track thickness is changed to
h=0.4 nm and
h=1 nm, and it is not affected by the thickness of the circular track. In traditional SkHE, the skyrmion topological number
Q and the movement direction are the key factors affecting the process, as shown in the gray section of
Figure 2(a). Similar to the traditional SkHE, the Hall angle also demonstrates a notable reliance on the direction of motion and the skyrmion topological number
Q in a circular ring nanotrack. For a skyrmion with
Q=-1, a CCW current drives it towards the inner side, while the CW current drives it towards the outer side of the circular ring nanotrack. The situation is exactly the opposite for a skyrmion with
Q=1.
Based on the obtained results, it can be concluded that the movement of skyrmion in circular ring nanotracks is similar to SkHE. However, SkHE does not occur in the straight nanotrack when
. This leads us to assume that the SkHE observed in circular rings may be due to the asymmetry inner and outer boundaries of the circular ring. In
Figure 2(b), we show the variation of
with the drift distance
Rd for different
Rm. It can be found that
are linearly related to
Rd in both directions, it is just a negative correlation in CW and a positive correlation in CCW. The situations of
Q=1 are the opposite. Furthermore, special attention should be paid to the fact that the stable positions of the skyrmion, represented by the intersection of the velocities in two directions, are significantly different for circular ring nanotracks with different sizes. The stable position of a skyrmion is directly influenced by the size of the circular ring nanotrack. As the size of the ring decreases, its curvature becomes more prominent, leading to reduced symmetry between its inner and outer boundaries. Conversely, larger circular rings appear more symmetrical with a smaller difference between the inner and outer boundaries. Due to the asymmetric boundary of the circular ring nanotrack, the stable position of the skyrmion drifts farther from the middle line as the circular ring nanotrack size decreases. When a skyrmion moves along a straight track, it remains stable on the centerline without deviation for a certain distance. However, when a skyrmion moves along a curved track, the collision of the skyrmion with the asymmetric boundaries of the track results in periodic changes in the skyrmion's Hall angle. This, in turn, leads to periodic increases and decreases in velocity. Therefore, the period of such collisions is uncertain. When we modify the size of the nanotrack or the driving current density, the SkHE still exhibits an oscillation mode, but the period changes accordingly.
To deeply understand the influence of asymmetric boundaries on the movement of the skyrmion in the circular ring nanotrack, the Thiele equation is used. Since the size of skyrmion changes is not obvious and the topological structure and stability are maintained, we still treat the skyrmion as a rigid object in asymmetric track. This rigid approximation allows the Thiele equation for the skyrmion is written as:
where the first term is the Magnus force,
G is the gyrovector related to the topological number
Q of skyrmion, and
G=(0, 0, 4π
Q) is a constant for Néel skyrmion. The second term of Eq. (2) stands for the dissipative force and the
D comes from a dimensionless matrix related to the dissipative force, being its components
Dxx=
Dyy=
D.
is the velocity of the electron, and
l is the lattice constant of the material. The
is the velocity of skyrmion that can be decomposed into two orthogonal parts (the tangential and radial directions of the circular-ring nanotrack). The third term
represents the force acting on the moving skyrmion and
, including the forces from the boundary and driving current.
For
, we can get the simplification form of skyrmion velocity from Eq. (2):
Considering
D=|
G| and
the function of velocity can be further simplified as:
It can be concluded from Eq. (4) that the
is determined by the sign of
and the product of
G (
Q) and
. The force
acting on the skyrmion mainly comes from the repulsive force from the boundary, it is negative when the skyrmion moves towards the inner boundary and positive when the skyrmion moves towards the outer boundary. When a skyrmion with
Q=-1 moves in a counterclockwise direction and approaches the inner boundary, its velocity increases. However, when it approaches the outer boundary, its velocity decelerates. For skyrmions with
Q=+1, the situation is the opposite since
G is positive. This is consistent with the variations in skyrmion velocity described in
Figure 1(b).
The SkHE caused by the asymmetry of the inner and outer boundaries is also reflected in the drift in the skyrmion stable. We so calculated the relationship between the initial drift distance
Rd and the radius
Rm of the circular ring nanotrack under a range of values
. The results, as shown in
Figure 3(a), revealed that the skyrmion's stable position tends to drift further away from the middle line as the size of the circular ring decreases. As shown in the inset of Fig.3(a), even if the values
are different,
Rd is linearly correlated with the curvature
of the nanotrack. Because the symmetry of the boundary is directly affected by the curvature of the circular ring. Additionally, both
Figure 3(a) and the illustration in it show that the dependence of skyrmion size
(the red line) on the
Rm is the same as that of
Rd. Essentially, the force from the asymmetric boundary changes the size of the skyrmion and thus changes its stable position on the track.
In addition to the size
Rm of the circular ring nanotrack, the width also affects the force exerted by the boundary on skyrmion [
46,
47].
Figure 3(b) illustrates the drift distance
Rd changes with the track’s width when the
Rm is fixed at 80 nm, 100 nm, and 120 nm. Generally, it is still true that the smaller the
Rm, the more obvious the skyrmion stable state drift from the initial position. However, for the same
Rm values, the stable state can be divided into four regions at different width values, as shown in
Figure 3(b). In region I (0<
W<40 nm), the tracks are too narrow for a stable skyrmion. In region II (40<
W<70 nm), although the asymmetric boundary reduces the size of the skyrmion, the increase in width significantly increases the size of the skyrmion, causing
Rd to gradually increase until the increase in width has little effect on the size of skyrmion. When
W>70 nm, the skyrmion size does not increase significantly, so the
Rd in region III and IV is consistent with the effect of the size of a circular ring. For larger
W in region IV (
W>140nm), the corresponding
Rd ≤0.15 nm and the asymmetric boundary drift can be almost ignored. This is because as the track becomes wider, although the force caused by the asymmetric boundary still exists because the boundary is too far away from the skyrmion, the stable position will hardly move much.
In order to confirm our analysis, we also calculate the skyrmion velocity on a wider nanotrack. Our findings reveal that as the track width is very wide, the periodic oscillation of
(depicted in
Figure 1(b)) diminishes or entirely disappears, and
(depicted in
Figure 1(c)) comes close to zero. To verify this, we tested the case of a width of
W=150 nm. As shown in
Figure 4(a) and (b),
no longer displays an obvious fluctuation, but remains close to the velocity in the long straight track, and
is also close to zero, as expected. The reason for this is for a fixed
Rm, the boundary’s symmetry for the skyrmion's stable position remains the same. However, the repulsive force from the boundary weakens as the track gets wider. This force is related to the size of the skyrmion and its distance from the boundary. As a result, the skyrmion remains unaffected by the boundary.
Since the asymmetric inner and outer boundaries of a circular ring can generate the SkHE, we have also studied the impact of asymmetric boundaries on the circular rings that already have the SkHE. As illustrated in
Figure 4(c) and (d), when
the Hall angle in a circular ring track is greater than that in a straight track, the SkHE is enhanced. On the other hand, when
the Hall angle in the circular ring nanotrack is smaller than that in a straight track, the SkHE is weakened. Therefore, the asymmetric boundaries may be used to strengthen or weaken the SkHE. That is utilizing a circular ring nanotrack offers a promising way to manipulate the SkHE in devices.
According to research, SkHE can also be induced by inhomogeneous DMI [
48]. In this section, we compared the effect of asymmetric boundary with that of inhomogeneous DMI. We created straight nanotracks with inhomogeneous DMI, which had the same width as the circular ring nanotracks. Firstly, we compared the effect of non-uniform DMI on the skyrmion’s stable position with the circular ring. The graph in
Figure 5(a) illustrates the relationship between the DMI change values and the drifted distance of the skyrmion's stable position for different initial DMI values. When there was no change in DMI, the skyrmion stayed in the center of the track. However, when there was a slight DMI change, the skyrmion drifted from its middle line. The
DMI gradient in can cause a similar drift distance of the skyrmion’s stable position with the circular ring with
Rm=80 nm.
We also conducted an analysis to determine the factors that influence the Hall angle in the inhomogeneous DMI distribution tracks. The trajectory of the skyrmion in the modeled nanotrack is shown in
Figure 5(b). It was observed that the presence of an inhomogeneous DMI distribution in the nanotrack causes the skyrmion to move toward the boundaries as it moves in a circular ring. Additionally, we found that the sign of the Hall angle depends on the sign of
and skyrmion topological number
Q.
Figure 5(c) illustrates the relationship between the Hall angle and the value of the DM changes value for various current densities. Compared with the circular ring,
DMI gradient can only yield a Hall angle of about
for
driving current. For high current density, a larger DMI gradient is required to cause a
Hall angle. Therefore, the SkHE caused by the circular ring is much more significant than that caused by the inhomogeneous DMI.