Discrimination of skyrmion chirality via spin–orbit and –transfer torques for logic operation

Recently many works on magnetic memories and logic circuits, which use a magnetic skyrmion have been reported. Previously we micromagnetically simulated a method to switch a chirality of a magnetic skyrmion formed in a magnetic thin film by introducing a pulsed heat spot. In this paper, we propose a method to discriminate the chirality of a skyrmion in a branched nanowire by using spin–orbit torque (SOT) and spin-transfer torque (STT), and confirm the validity of the method by using simulation. The simulated results show that the motion changes depending on the chirality when additional SOT is applied on a skyrmion moving in a branch by STT. This method can be used as a fundamental building block for electrical detection in memory and logic devices using the chirality of skyrmions as a data bit in addition to the presence (and polarity) of the skyrmions as conventionally used, which can be lead to multiple-valued operation.

Here, m, H eff , γ, and α are the unit magnetization vector that represents the direction of a local magnetic moment, the effective magnetic field acting on the magnetization, the gyromagnetic ratio, and the Gilbert damping constant, respectively. The velocity u is a vector along the direction of electron motion, with an amplitude of ( www.nature.com/scientificreports/ where j, P, g, μ B , e, and M s are the current density, spin polarization, the Lande g-factor, the Bohr magneton, the electron charge, and the saturation magnetization, respectively. β is the dimensionless parameter for STT. θ H , h, and n s are the spin Hall angle, the thickness of the film and the unit vector of the spin Hall torque, respectively. The effective magnetic field is calculated from the magnetic energy density by H eff = − 1 M s δε δm . For the magnetic energy, the exchange, anisotropy, demagnetizing, and Dzyaloshinskii-Moriya exchange interaction (DMI) energies are taken into account in the simulation 35,36 .
Here, A, K u , H dem , and D are the exchange stiffness constant, the uniaxial anisotropy constant, the demagnetizing field, which is calculated numerically, and the DM exchange constant, respectively. Here we assumed the interfacial DMI between the ferromagnetic thin film and the non-magnetic layer underneath. A free boundary condition was used at the edge of the film 36 .
Here, θ is a polar angle of a magnetic moment, n is the normal direction at the film edge, and ξ is a characteristic length determined by D 36 .
In this paper, two types of simulation were performed. In the first one, we investigate the effect of the D value on the direction of the skyrmion motion by SOT. The second one is the simulation of the skyrmion motion in the nanowire with a branch by SOT and STT. In the first simulation, we use the magnetic strip with 512 nm in length, 128 nm in width and 1.4 nm in thickness. In the second simulation, the branch structure is made by attaching the same strips used in the first simulation. The branch angle is varied from 15 to 90 degree from the nanowire direction. These strips are divided into rectangular prisms with their dimensions of 0.5 × 0.5 × 1.4 nm 3 in simulation.

Results and discussion
First, we show the change of the direction of the motion of the Néel and Bloch type skyrmions by SOT depending on the chirality by analytical model 31 . Figure 1 shows representative results. In the case of the Bloch type skyrmion, the direction of motion by SOT changes 180º depending on the chirality [see Fig. 1a]. This means that we can differentiate the chirality of the Bloch type skyrmions by detecting the opposite directions of motion. On the other hand, the direction of motion of the Néel type skyrmion is 90º, which requires further consideration as detailed below.
In the magnetic thin film, not only the diameter of the skyrmion, but also the direction of the magnetization on the domain wall surrounding the skyrmion changes by the D value. As shown in Fig. 1b, a Bloch type skyrmion in which the magnetization is parallel to the domain wall appears for the D = 0 erg/cm 2 case. The magnetization changes as the D value increases, and Néel type skyrmion appears for the D = 1.2 erg/cm 2 case. Figure 1c shows the simulated results of the change of the direction of the skyrmion motion by the D value and chirality. For the D = 0 erg/cm 2 cases, the same results with the analytical model are obtained. However, the directional angle dependent upon the chirality decreases with increasing the D value. In our previous paper 21 , in which we reported the control of the chirality by a heat pulse, we use the Bloch type skyrmion with D = 0.6 erg/cm 2 . In that case, the clockwise (CW) skyrmion moves to the + x direction by SOT, however the counter clockwise (CCW) skyrmion moves to the -y direction predominantly (Fig. 1c, D = 0.6 erg/cm 2 ). It shows that the branch structure to discriminate the chirality of the skyrmion should be designed considering the D value. Figure 2 shows the branch structure to discriminate the chirality of the Bloch type skyrmion. In this figure, the current flows to the -x direction (electron moves to the + x direction), and the skyrmion moves to the + x direction by STT. When SOT with the + y direction is also applied, the CW skyrmion with D = 0.6 erg/cm 2 moves to the + x direction as mentioned before. It moves to the + x direction at the branch, moving into the upper branch. On the other hand, the CCW skyrmion moves to the -y direction by SOT, moving to the below branch by the combination of STT and SOT.
We then investigate the change of the branch directions by chiralities using simulation. Figure 3 shows the results. Figure 3a,b show the motion of the CCW and CW skyrmions for 40 ns, respectively. Here, the branch angle is 45º, and the current density is j = 0.3 × 10 12 A/m 2 with P = 0.7 and θ H = 0.1 rad (red lines) or 0 rad (green dashed lines). These figures show that the direction of motion and the resulting branches for detection changes by the chirality and SOT (red line). The figures also show the simulated results without SOT (green dashed lines). The CW and CCW skyrmions move almost straight with minor departure by the branch. These results confirm that the branch direction does not change only by the chirality and SOT is required to control the branch. In all cases, the skyrmion moves to the center of the branches by the demagnetizing field and then moves along the branch. Figure 3c shows the demagnetizing field in the branch. The field becomes the maximum at the branch center because the width of the strip is maximum. That is the stable point for the skyrmion, and hence the skyrmion moves to the center of the branch first and then it moves along the branch.
(2) u = jPgµ B /(2eM s ),  Fig. 3e]. For the small current cases (j = 0.1 and 0.15 × 10 12 A/m 2 ), it takes long time to pass through the branch. As an example, the motions for 160 ns is shown in the figure. In Fig. 3d, the CW skyrmion moves to the upper branch for all cases. The CCW skyrmion moves to the bottom branch when the current is larger than 0.20 × 10 12 A/m 2 , however it crashes the branch center and disappears when j = 0.15 × 10 12 A/m 2 . In Fig. 3e, the CW skyrmion moves to the upper branch when the current is smaller than 0.2 × 10 12 A/m 2 , however it moves to the lower branch and crashes into the strip edge and disappears when the current is larger than 0.25 × 10 12 A/m 2 . The CCW skyrmion moves to the lower branch for all cases, however it crashes into the strip edge and disappear when the current is larger than 0.25 × 10 12 A/m 2 . Note that SOT is required to change the motion direction. However, when the spin Hall angle is large (~ 0.3 rad), the current range for the main branch decreases by  www.nature.com/scientificreports/ misdirectional motion or anihilation. These results show that the spin Hall angle of 0.1 rad is sufficient for the chirality discrimination, which is experimentally achievable 30,39 and larger angle is not necessary. Figure 3f-h show the phase diagram of the branch. Here, the branch angle and the current are varied between 15 and 90º, and between 0.1 and 0.4 × 10 12 A/m 2 , respectively. The circles show the branched cases and the crosses show the unbranched cases. The green (blue) symbols show the CW (CCW) skyrmion cases. The shaded regions show the angle and current ranges for the CW and CCW skyrmions to branch. Note that 0.1 × 10 12 A/m 2 is the minimum value we used in the simulation. Therefore, it is not the true minimum for normal operation. These figures show that the effect of the branch angle on the current range is small. The current range decreases because the maximum current decreases as the spin Hall angle increases. The misdirectional motion of the skyrmion does not depend on the branch angle, which is the same as that in Fig. 3d,e.
From these results, the skyrmion chiralities are proven to be discriminated electrically by SOT and STT. The branch direction of the Bloch type skyrmion moving by SOT and STT changes depending on the chirality. The spin Hall angle of 0.1 rad is sufficient for the chirality discrimination, which is experimentally achievable and larger angle is not necessary. The effect of the branch angle on the current range is minor, which is advantageous for device miniaturization.
In general, the trajectories of the skyrmion are not parallel to the electron flow direction by the Magnus force. The angle between the trajectories and the electron flow is determined by the ratio between the Gilbert damping constant, α, and the parameter for STT, β 9 . When β = α, the corresponding angle become zero, and the skyrmion moves parallel to the electron flow direction. The angle also changes by the chirality with SOT. This study shows that the chirality of the skyrmion can be determined by using the change of the angle by chirality. Note that in this study, we selected β = α for simplicity. For β ≠ α case, we can also determine the chirality by adjusting the branch angle.
Here, the influence of the Joule heating induced by an electrical current flow can be ignored due to the following reason. In the branched system here, the current density becomes strongly inhomogeneous. It is expected that the temperature induced by the Joule heating by the distribution of the current density increases and influences the chirality. This can be similar to the case we previously reported where the chiralities of the skyrmion can be switched by a local pulse heat spot with realistic parameters, a Gaussian shape with a diameter of about 60 nm, the maximum temperature of 550 K and a pulse width of 10 ns 21 . However in the branch geometry with the width of 128 nm and the length of 512 nm for simulation, temperature may be dissipated within large area and the corresponding temperature gradient is much smaller than that with the heat spot we reported previously. It is hence concluded that the chirality does not change by the thermal effect in this study.

Conclusions
We proposed a method to electrically discriminate the chirality of a skyrmion using a combination of STT and SOT. This proposed method was investigated and confirmed by micromagnetic simulation. This method is found to be highly effective to realize skyrmionics memory and logic devices with using the chirality of the skyrmion, allowing to achieve multiple-valued operation in combination with the presence and polarity of the skyrmion.

Data availability
The datasets generated during the current study are available from the corresponding author on reasonable request.