Control of Néel-type Magnetic Kinks Confined in a Square Nanostructure by Spin-Polarized Currents

Magnetic skyrmion in chiral magnet exhibits a variety of unique topological properties associated with its innate topological structure. This inspires a number of ongoing searching for new topological magnetic textures. In this work, we used micromagnetic simulations and Monte Carlo simulations to investigate an exotic Néel-type magnetic kinks in square-shaped nanostructures of chiral magnets, which performs rather stably in the absence of magnetic field. The individual magnetic kink can reside in one of the four possible corners, and carry possibly upward or downward core polarity, constituting eight degenerate states. In addition, these kinks also exhibit unique behaviors of generation, stability and dynamics, as revealed by micromagnetic simulations. It was found that such kinks can be created, annihilated, displaced, and polarity-reversed on demand by applying a spin-polarized current pulse, and are easily switchable among the eight degenerate states. In particularly, the kinks can be switched toward the ferromagnetic-like states and backward reversibly by applying two successive current pulses, indicating the capability of writing and deleting the kink structures. These findings predict the existence of Néel-type magnetic kinks in the square-shaped nanostructures, as well as provide us a promising approach to tailor the kinks by utilizing the corners of the nanostructures, and control these states by spin-polarized currents. The present work also suggests a theoretical guide to explore other chiral magnetic textures in nanostructures of polygon geometries.


INTRODUCTION
A magnetic skyrmion is a topologically stable configuration often observed in chiral magnets with broken inversion symmetry. The nanoscale skyrmion shows particle-like behavior, as it can be moved, created, and annihilated. These characters make it promising candidate as information carrier for future memory devices and logic devices [1][2][3][4][5][6]. It has thus aroused intense research efforts in recent years, leading to a series of breakthrough achievements in manipulation of skyrmion states, via injected spin-polarized currents, external electric-field, and so on [7][8][9][10][11][12][13][14]. For the application of skyrmions in spintronic devices, many feasible designs have been proposed for confining the skyrmions in geometric nanostructures, such as nanostripes and nanodisks, which may allow the precise control of individual skyrmions [8][9][10]. These achievements pave the way towards all-electrical manipulation schemes [2,3], and underpin the skyrmion-based information storage concepts, such as the skyrmion-based racetrack memory, high density magnetic random access memory, and logic gates, etc. [2-6, 8, 9].
These fascinating physical properties also inspired an ongoing search for new types of magnetic topological textures, for instance, fractional skyrmion emerging in various chiral and frustrated magnets [15][16][17][18][19][20][21], which exhibit some unusual physical phenomena distinctly different from skyrmion-host chiral magnetic materials such as MnSi and Fe 1-x Co x Si systems [22,23]. Moreover, some recent observations indicated that fractional skyrmions also emerge in geometric nanostructures due to the effects of geometric confinement and shape anisotropy [24][25][26][27][28]. For example, it was observed in Lorentz transmission electron microscopy (TEM) images of FeGe nanostripes that, the skyrmions survive when the width of the nanostripe much larger than the single skyrmion size, whereas some merons form at the edges that are not large enough for accommodating a complete skyrmion [24,25].
Generally, the presence of edges and corners in nanostructures can be utilized to tailor the magnetic textures and modify their dynamics behaviors. For the skyrmion confined in ultrathin film nanostructures with Dzyaloshinskii-Moriya (DM) interaction, the boundary constrictions naturally make the magnetization orientation undergo the 180°rotation at the edges, forming the so-called kink or π domain wall configurations [2,15,[29][30][31][32]. Depending on the crystal symmetry of chiral magnets, two distinct types of chiral kinks, namely, Néel-type and Blochtype kinks, are favorable in interfacial and bulk DM interaction systems respectively [29,30]. It had been theoretically proposed in chiral magnets that the skyrmions, antiskyrmions, and other magnetic configurations can be naturally interpreted in terms of chiral kinks. These kinks carry a topological charge and allow to construct new topological particle-like states [15,31,32]. In addition, the previous studies demonstrated that some special magnetic textures often appear in the corners of the polygon geometries like triangles, squares, rectangles, in which the corners may also act as pinning sites for the domain wall motion [33][34][35][36]. This indicates that the nanostructure with polygon geometries may provide a unique platform for studying the kink structures, favoring the control of the magnetic states, by utilizing their edges and corners. Certainly, there are more open questions to be explored along this line. For instance, is it possible to realize some new kink structures beyond the skyrmions and merons in nanostructures? Is it possible the kink structures stable in nanostructures without the assistance of external magnetic field? Is it possible to manipulate these kinks by pure electric currents? These issues are very critical for our understanding the physical properties of the new type chiral magnetic structures in confined nanostructures.
Inspired by this motivation, the purpose of our work is to explore new kink structures in nanostructures by using micromagnetic simulations and Monte Carlo simulations. In this work, we aimed at studying the square-shaped nanostructures with width smaller than the single skyrmion size. We demonstrated an exotic Néel-type kink structure, which can stably exist in the corners of the square-shaped nanostructures in the absence of external magnetic field. Such kinks possess eight switchable degenerate states and can be created, annihilated, displaced, and polarity-reversed by applying spin-polarized currents.

MODEL AND SIMULATION METHODS
In this work, the magnetic state in an ultrathin square nanostructure of chiral magnets is described by a classical Heisenberg model. The Hamiltonian consists of ferromagnetic exchange interaction (J FM ), DM interaction (D), anisotropy (A), and Zeeman term (B z ) [8,37,38]: where the magnetic moments are imposed on a two-dimensional L × L square lattice with free boundary conditions. m i denotes the magnetic moment at site i in the xy-plane, with fixed length |m i | 1.0. The magnetic field B z B z e z is applied normal to the xy-plane and easy-axis anisotropy (A > 0) is along the ± z-axis. In the simulations, we have neglected the dipole-dipole interactions. This treatment is suitable for chiral magnets with strong DM interaction and weak dipoledipole interactions such as Pd/Fe bilayers on Ir(111) substrate [39,40]. For Pd/Fe/Ir(111) the magnetic dipolar energy is very weak in energy scale in comparison with the ferromagnetic exchange energy and DM interaction, and thus it can be effectively included into the anisotropy energy. We choose typical parameters for Pd/Fe/ Ir(111) as J FM ≡ 2.95 meV, D/J FM 0.27, A/J FM 0.05, and a typical lattice constant a 5 Å [39,40].
To investigate the dynamics of the magnetic structures driven by the spin-polarized current, we numerically solved the Landau-Lifshitz-Gilbert (LLG) equation by using fourth-order Runge-Kutta method [8,10,41,42]: with the effective field B eff i −(1/Zc)(zH/zm i ), the gyromagnetic ratio c g s μ B /Z (here g s is the electron spin g-factor, and μ B is the Bohr magneton), and Gilbert damping α. The first and second terms describe respectively the gyromagnetic precession and the Gilbert damping, and the third term T denotes the spin transfer torque (STT) due to the spin-polarized current.
The skyrmion motion can be driven by spin-polarized current flowing in the nanostructure along either the in-plane or out-ofplane direction. For simulations of the skyrmion dynamics induced by the current-in-plane (CIP) injection, the corresponding torque T CIP is given by the following form [41]: where the first and second terms represent the coupling between magnetic moments and spin-polarized current j CIP via the spin transfer torque and via the non-adiabatic effects, with β the nonadiabaticity factor. M s is the saturation magnetization, P e is the polarization rate of the electric current, and e is the elementary charge.
For the current-perpendicular-to-plane (CPP) injection, the current-induced spin transfer torque T CPP includes an in-plane Slonczewski torque and an out-of-plane field-like torque [8,42,43]: where u Z e j CPP P e 2dMs is the Slonczewski torque coefficient, d the film thickness of FM layer, j CPP the current density, m p is the electron polarization direction, and ξ is the amplitude of the out-of-plane torque relative to the in-plane one. In the simulations, we fixed the coefficient ξ 0.2, the polarization rate P e 0.4, gyromagnetic ratio c 1.0, Gilbert damping α 0.3, and non-adiabaticity factor β 0.1, as the typical parameters for studying the currentinduced skyrmion dynamics [8,10,41]. The simulated time t is measured in the units of τ ≡ Z/J FM ∼ 0.2 ps, and current density j are scaled by κ CIP ≡ 2eJ FM /P e a 2 Z ∼ 1.4 × 10 13 Am −2 for the CIP case, and κ CPP ≡ 2edM s J FM /Z 2 P e c ∼ 8.5 × 10 13 Am −2 for the CPP case, respectively [41]. Here we chose the typical parameters for Pd/ Fe/Ir(111) as d 0.4 nm and M s 1.1 MA/m, and used these values to estimate the units of simulated time and current density.
To get the zero-temperature equilibrium state, we adopted a specific simulation scheme: the lattice was initialized as a paramagnetic phase at sufficiently high temperature, and annealed for obtaining the state at a low temperature, by using the Metropolis Monte Carlo simulation combined with overrelaxation algorithm [17,18,44]. Then the configuration was further relaxed for the equilibrium state, by solving the LLG equation under T 0 with long enough equilibration time.
In addition, the kink dynamics is controlled by injecting a CIP-type or CPP-type spin-polarized current pulse on demand. We introduced the schemes in our simulations for setting the amplitude and duration of current pulses as follows: 1) First, we tested the effect of amplitude of current density on the kink dynamics and chose the suitable amplitude of current pulses for control (generation, creation, annihilation, displacement or polarity reversal) of the kink states. 2) Then we carefully tracked its dynamics and turned off the current once the germinal kink appears or kink shifts (in this procedure, the current duration is determined), and further relaxed the system by solving LLG equation under zero-current condition to get an equilibrium kink state.

Size Effect of Nanostructures on the Magnetic Structures
We first investigated the size effect of the nanostructures on magnetic structures. The simulations were carried out on a L × L square lattice with free boundary conditions. The zero-temperature equilibrium states for the lattices with different sizes were obtained, and some of the typical configurations were presented in Figure 1. Here a magnetic field B z −0.008 J FM (corresponding to ∼ −0.1 Tesla, as estimated with the magnetic moment ∼3.0 μ B obtained from density functional theory for Pd/Fe/Ir(111) in Ref. [40]) is applied along the −z-axis, which is sufficient for creating the skyrmion states here. In the following, we used an external magnetic field to generate the exotic magnetic kink textures in the nanostructures, and then we demonstrated an alternative way for creating kink by applying spin-polarized current. As one may see in Figure 1A, a bubble domain arises in the square nanostructure at small lateral size L 20 (∼10 nm). As the size increases to L 24 (∼12 nm), a kink-like structure forms in the corner of the square nanostructure, surrounded by the ferromagnetic domain, as will be discussed below. For L 50 (∼25 nm), it was observed that some half-skyrmion-like (meron-like) structures with topological charge |Q| ≈ 1/2 and elongated stripes (fractional skyrmion) with topological charge 0 < |Q| < 1 emerges at the edge of the lattice, together with a kink-like structure that appears in the corner. Further increasing the size, the magnetic states evolve into multi-domains composed of spiral domains, edge-merons, and skyrmion as seen in Figure 1D for the relatively large size L 80 (i.e., ∼40 nm). It is interesting that the kink-like, meron-like, and skyrmion states tend to nucleate in the corner, the edge, and the inner region of the nanostructure, respectively. These are reasonable likely due to the different boundary constrictions on the formation of magnetic structures, which manifests the strong boundaries and geometric confinement effects in the nanostructures. It is noted that some similar edge states are also found experimentally and theoretically in diverse magnetic materials that host skyrmions with constricted geometry, as a consequence of the effects of geometric boundaries and confinements [37, 45].

Néel-Type Magnetic Kinks
Now we focus on the intriguing features of kink-like structures and analyzing their current-induced dynamics in this paper. For this, we adopted a small square shaped nanomagnet consisting of 24 × 24 square lattices (∼12 × 12 nm) for studying the kink-like structure, in the rest part of this work. We first tested the stability of the kink-like structure in Figure 1B once the magnetic field was removed. For this, the magnetic structure was relaxed for the equilibrium state by solving the LLG equation, in which we used the kink-like state shown in Figure 1B as the initial state, and set B z 0 in the calculation. It was found that the kink-like structure remains stable without the assistance of an external magnetic field, as presented in Figure 2A. Note that the magnetic structure enclosed by the dash lines in Figure 2A is called Néel-type magnetic kink in this work, in which the magnetic moments undergo 180°rotation from the upward direction at its center to the downward direction in the periphery. The schematic magnetic configuration for a Néel-type kink is displayed in Figure 2D. For simplicity, we call the kink-like state as the kink state hereafter.
To explore the possible ground state for the nanomagnet of 24 × 24 square lattices at B z 0, the system was initialized from a the paramagnetic state at sufficiently high temperature T, and cooled down gradually until it reaches a very low T under B z 0, using Monte Carlo simulations with simulated annealing technique [17,18,44]. Then the system was further relaxed for the zero-temperature equilibrium state by solving LLG equation. As a result, a spiral structure forms as shown in Figure 2B. The numerical calculations indicate that the total Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 680698 energy of the spiral structure is a little smaller than that of the Néel-type kink structure. These results suggest the spiral structure to be the possible zero-field ground state, while the Néel-type kink structure to be a metastable state at B z 0, other than the minimum-energy state.
To further estimate the stability of the metastable Néel-type kink state under B z 0, we calculated the reduced local effective acting on the i-th magnetic moment m i . Interestingly, it was found in Figure 2C that the b i distribution and the magnetic kink state in Figure 2A are similar in spatial profile, and almost all the magnetic moments in the kink structure are oriented in parallel to the local field b i , making the metastable kink structure self-sustaining. In addition, it was noted that the anisotropy energy H ani −A i (m i · z) 2 contributes to the effective field in terms of −(1/Zc)(zH ani /zm i ) (2/Zc)Am z i e z . This suggests that the anisotropy plays an important role in guaranteeing the robustness of the metastable kink state, since the effective field (2/Zc)Am z i e z makes the stabilization of the out-of-plane magnetization component m z i [46,47]. In this regard, we reckon that the kink structure is a robust state after the withdrawal of the magnetic field, due to the kink state is trapped in an energy valley in the configurational energy landscape [47].
In fact, one may take account of the kink structures from the prototypical Néel skyrmion, whose configuration can be viewed as a coplanar spiral with magnetic moments lying in a plane perpendicular to the xy plane [1,38]: ⎧ ⎨ ⎩ m(r) sin q · r e q + cos q · r e z core − up, p 1 m(r) sin q · r + π e q + cos q · r + π e z core − down, p −1 , where the spatial variables of magnetization are defined in polar coordinates r (r cos φ, r sin φ), and the magnitude of spiral wave vector is q π/R, with the direction denoted by unit vector e q cos φe x + sin φe y . For the isolated Néel skyrmion configuration, the core magnetization points upward or downward, and smoothly changes to the opposite direction in the peripheral circle with radius R. We considered skyrmion structure with core-up (core-down) magnetization carries core polarity p 1 (p −1), as the skyrmion with p 1 depicted in Figure 2E. Here, R is used to define the radius of skyrmion, with 0 ≤ r ≤ R and the azimuthal angle 0 ≤ φ ≤ 2π for a single skyrmion. To proceed, one may divide an isolated skyrmion configuration into four quadrants along the dash lines, as illustrated in Figure 2E. The magnetic structures in quadrant I, II, III, and IV can be described by Eq. 5, with the azimuthal angle 0 ≤ φ I ≤ π/2, π/2 ≤ φ II ≤ π, π ≤ φ III ≤ 3π/2, 3π/2 ≤ φ IV ≤ 2π, respectively. The kink may appear in one of the four possible corners, i.e., the upper left (UL), upper right (UR), lower left (LL) or lower right (LR) corner of the nanostructure in the simulations, as shown in Figure 2F. The LL, LR, UR, and UL kink structures correspond to the well-defined kinks residing in the quadrant I, II, III, and IV, respectively. In addition, we define the polarity of the kink with core-up (core-down) as p 1 (p −1). Therefore, the kinks in the nanostructure can carry the core polarity  Figure 1A. Here the pie chart is generated from a set of data points for φ 2, 4, . . ., 360°at intervals of 2°.
Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 680698 p ± 1, and can reside in one of the four possible corners, constituting 8 degenerate kink states. By comparison, the single skyrmion often appears in circular shaped nanodisks, and possesses only two degenerate states with core polarity p ±1 [8,[48][49][50]. In this sense, the kinks in the square shaped nanostructure have the special feature of multiple degenerate states.

Generation and Switching of Kinks by In-Plane Current Pulse
In this section, we studied the generation of kinks by injecting into the nanostructure an in-plane spin-polarized current pulse.
Here, the direction of the CIP injection is defined by the angle φ with respect to the x-axis (see in Figure 3B), and φ is tunable. Note that in the following study, no external magnetic field was applied.
The simulations start from the initial ferromagnetic phase, with all magnetic moments aligning along z-axis at t 0.0 ps (see Figure 3A). We first tested the effect of CIP on the kink dynamics, using a moderate current density j CIP 0.1κ CIP ∼ 1.4 × 10 12 Am −2 with current duration t dur 250 τ ∼ 50 ps. Figure 3A shows the formation process of the UR kink for a typical case of φ 10°. At the beginning (t 20 ps) of the current duration, the embryo of kink pattern emerges in the UR corner of the nanostructures. This pattern gradually enlarges, and becomes a rough kink at t 50 ps. Subsequently, the current is turned off and the system evolves into an equilibrated kink state at t 120 ps.
Further simulations demonstrated that the generation of kinks is sensitive to injection angle φ, with the simulated results summarized in Figure 3C. It was found that the UR, UL, LL, and LR quarters can also be created respectively for −42°F

IGURE 4 | (A)-(H)
The gradual shift of kink from UR to LR corner, driven by in-plane current pulse with current density J CIP 1.4 × 10 12 Am −2 , current duration t dur 50 ps, and injection angle φ 135°. (I) Schematic illustrations of switching sequences for the kinks driven by in-plane current pulse. Here we used the straight red arrows to mark the injection direction of current J CIP , and the curved pink arrows to denote the switching sequences of magnetic kinks.
Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 680698 ≤ φ ≤ 14°, 48°≤ φ ≤ 104°, 138°≤ φ ≤ 194°, and 228°≤ φ ≤ 284°, while some bubble domain states form beyond these angles φ (see the while region in the Figure 3C). Note that the kinks with p −1 presented here are generated from the initial ferromagnetic phase with magnetization along the z-axis. If the initial ferromagnetic phase is magnetized along the −z-axis, the lattices will evolve to the kinks with p 1. Therefore, all the eight degenerate kink states can be created by tuning the direction of CIP injection with two different initial ferromagnetic orientations.
Next, we investigated the switching between these kinks using in-plane current pulse, which is fundamental to understand their dynamics properties. The current pulse (j CIP 1.4 × 10 12 Am −2 , current duration t dur 50 ps) with varying angle φ is exerted on the initial UR kink state with p −1. Simulated results reveal that the kink can transfer from UR corner to the LR corner in a clockwise (CW) direction, driven by current with 116°≤ φ ≤ 178°. However, it was noted that the UR kink cannot move to LL corner in a counterclockwise (CCW) direction or to diagonal UL corner with the adopted simulation parameters. Figures 4A-H show the dynamic process for the typical case of φ 135°. We can see clearly that the upper part of the UR kink is pushed towards the right-edge of the nanostructures, and the lower part is simultaneously dragged to the bottom of the nanostructures (see Figures 4A-D). Although the entire kink pattern deforms largely in this process, it gradually turns in a CW direction and eventually moves to the LR corner at t 50 ps (see Figures 4E,F). After that the current is switched off to zero for reaching the equilibrated kink state in the relaxation procedure. Moreover, simulations for the kinks with p −1 indicate that their switching sequences can be summarized as follows: UR kink → LR kink → LL kink → UL kink with suitable angle φ, as shown schematically Here two successional current pulses with injection angle φ 1 135°and tunable injection angle φ 2 , current density j CIP 2.8 × 10 12 Am −2 and duration t dur 20 ps is considered. These two pulses are used to annihilate and create the kink state, respectively. (D) The phase diagram for the equilibrated states at t 520 ps after relaxation as a function of angle φ 2 was summarized. The various kinks, spiral states, and bubble states are marked by the colored pieces, striped pattern and grid pattern, respectively. The color map below is used to scale the magnetization components along z-axis m z i . Here the pie chart is generated from a set of data points for φ 2 1, 2, . . ., 360°at intervals of 1°.
Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 680698 in Figure 4I. In this procedure, the switching sequences is in a CW direction, and the kink polarity has not been reversed by the CIP injection.

Annihilation and Creation of Kinks by In-Plane Current Pulses
In this section, we presented an effective way to annihilate and create the kink by applying two successive CIP pulses, as the results shown in Figure 5. Our simulations start from the initial UR kink state with p -1 in Figure 5A. We first found that the kink can be annihilated by applying a CIP current pulse, with a larger current density j CIP 2.8 × 10 12 Am −2 and an appropriate injection angle φ 1 135°. Here two successive current pulses are used to annihilate and create the kink states respectively, with injection angle φ 1 135°for the first pulse, and tunable injection angle φ 2 for the second one. The current duration t dur 200 ps is considered for both two pulses, as shown in Figure 5C. Figure 5A displays the annihilation of the kink driven by the first current pulse with angle φ 1 , whose pattern rotates along a CW direction and disappears at the LR corner of the nanostructure, forming the ferromagnetic-like state at t 40 ps. The ferromagnetic-like state remains at 40 ps ≤ t ≤ 200 ps, and we may see that most magnetic moments orientate along the z-axis in the ferromagnetic-like state, though some remaining magnetic moments align in the xyplane at the LR and LL corners of the nanostructure. More interestingly, the simulated results reveal that some other kink states can be created from the ferromagnetic-like state, by applying the second current pulse with varying injection angle φ 2 . The typical case for φ 2 230°is presented in Figure 5B, in which the domains first appear and gradually grow in the LL, and LR corners during 200 ps ≤ t ≤ 240 ps. The new embryonic kink forms at t 280 ps and remains during 280 ps ≤ t ≤ 400 ps, and it finally evolves into a core-up kink in the UL corners after the relaxation with j CIP 0. In this process, the kink has been annihilated and created by two successive current pulses, in which the switching between the kink and ferromagnetic-like states is reversible. This indicates the capability of reversible writing and deleting the kink states.
Further simulations generate the phase diagram for the equilibrated states at t 520 ps as a function of φ 2 , as summarized in Figure 5D. It was noted that current pulses with 37°≤ φ 2 ≤ 55°, 143°≤ φ 2 ≤ 145°, and 320°≤ φ 2 ≤ 325°m ay create the kinks with p −1 in the UL, LL, and UR corners of the nanostructure, while current pulses with 226°≤ φ 2 ≤ 238°, 252°≤ φ 2 ≤ 253°∪ 314°≤ φ 2 ≤ 315°, and φ 2 268°generate the kinks with p 1 in the UL, LL, and UR corners, respectively. In addition, some bubble or spiral states form beyond these angles φ (see the striped-pattern and grid-pattern region in Figure 5D), with some typical configurations shown at the bottom of Figure 5D. Regarding early investigations on magnetic vortex, it was demonstrated that the vortex polarity reversal may be triggered by a CIP pulse through the formation of a vortexantivortex pair [51], or by an alternating CIP through the resonant dynamics [52]. However, for the skyrmion in chiral magnets, CIP is usually used for displacing the skyrmion, while cannot change the skyrmion polarity in the dynamics [8,41,42]. In this sense, the core polarity reversal of kink can be achieved by CIP pulse, which is analogous to that of vortex system [51,52]. From the simulation results and analyses in Generation and Switching of Kinks by In-Plane Current Pulse and Annihilation and Creation of Kinks by In-Plane Current Pulses, one may create, annihilate, displace, and reverse the kink on demand by adjusting the CIP pulse injections. Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 680698

Manipulation of Kinks by Out-of-Plane Current Pulse
In this section, we investigated the core polarity reversal of the kinks with out-of-plane current, which is also a fundamental issue for understanding their dynamics properties. One may build a magnetic tunnel junction (MTJ) or spin valve to locally address the kinks [5,53]. As schematically shown in Figure 6A, a typical MTJ structure consists of the top and bottom two layers of ferromagnet, which is separated by an ultrathin spacer of insulator. The top ferromagnetic layer is a free layer which presents a kink state, while the bottom one is a fixed layer with the magnetization M fix fixed along the z direction. When a spin current is injected into the bottom fixed layer along the z direction, the spin is polarized along z direction. The polarized spin current then flows through the insulating layer to the free layer, acting on its magnetizations with CPP-type STT, described by Eq. 4. One may reverse the polarization direction of the spin current, by injecting the spin current into the top free layer along the −z direction. Here we defined the current density j CPP > 0 and j CPP < 0 in Eq. 4 corresponding to the injection of spin current flowing along z and −z direction, respectively. Our simulations start from an initial UR kink with polarity p −1 as a general representative. Note that the CPP-type STT induced by current j cpp < 0 tends to align the magnetic moments in the free layer along the −z direction, thus the kink may be reversed to its image structure with polarity p 1. Contrarily, the kink with polarity p 1 can be reversed, once the spin current is injected along the z direction. We first studied the dynamics of the kink driven by various current densities j CPP and duration time t dur , and the phase diagram for the equilibrated states obtained after the current pulses was summarized in Figure 6B. The results show that the kink is immovable under low current densities j CPP < 0.04κ CPP (∼3.4 × 10 12 Am −2 ), due to the weak STT strength. It was seen that the kink may be pulled at current densities larger than the critical value j CPP 0.05κ CPP (∼4.3 × 10 12 Am −2 ). At an intermediate current density j CPP ∼ −4.3 × 10 12 Am −2 and duration time t dur 160 ps − 180 ps, the kink moves in a CW direction to near LR corner of the nanostructure, accompanying with the kink polarity reverse. For the current pulse with a relatively large current density 0.06κ CPP ≤ j CPP ≤ 0.11κ CPP (i.e., 5.1 × 10 12 Am −2 ≤ j CPP ≤ 9.4 × 10 12 Am −2 ) and duration ∼20 ps ≤ t dur ≤ ∼ 60 ps, it was found that the kink polarity is also reserved, while the UR kink shifts diagonally to the LL corner of the nanostructure. It was also seen in the phase diagram that the large j CPP and large t dur usually lead to the bubble domains, and some spiral states appear at 0.06κ CPP ≤ j CPP ≤ 0.07κ CPP with a short pulse duration.
To further investigate the repeatability of polarity reversals of kink driven by CPP pulses, we took an examination on injecting some discontinuous current pulses to the system, as presented in Figure 7A. The insert shows the kink states at some typical time before and after the pulse injections in the evolutions. It was observed that all the four CPP pulses reverse the kink polarity, and the pulse with spin current j CPP < 0 (j CPP > 0) is used to reverse the kink with polarity p −1 (p 1). The first and second CPP pulses with j CPP ± 4.3 × 10 12 Am −2 and t dur 160 ps rotate the kink in CW and CCW directions respectively (see state transitions ① → ② → ③), which is a reversible switching. Note that the CW and CCW rotational directions here depend on the core polarity of the pre-evolutionary kink. The third and fourth CPP pulses with j CPP ± 0.07κ CPP ∼ ± 6.0 × 10 12 Am −2 and t dur 40 ps move the kink to the diagonal corner of the nanostructure (see state transitions ③ → ④ → ⑤), which is also repeatable. For tracking the variations of the magnetic structures, the out-of-plane magnetization M z i m z i /N is used to characterize the reversal of kink, where N is the total number of magnetic moments. Although some fluctuant variations appear in the M z − t curve in Figure 7B, it is clear that M z flips to the opposite direction after every current pulse, in response to the core polarity reversal of the kinks. Therefore, we may use the CPP current to reverse the kink polarity and displace their position.

DISCUSSION AND CONCLUSION
Before concluding this work, we briefly discuss the experimental observations and thermal stability of the kink states. The X-ray holography allow the imaging of ultrafast magnetization dynamics in magnetic nanostructure with sub-10 nm spatial resolution time-resolved [6]. Using time-resolved X-ray microscopy, a number of experimental studies reported the direct observation of nanoscale skyrmions, current-driven Frontiers in Physics | www.frontiersin.org July 2021 | Volume 9 | Article 680698 skyrmions dynamics, and the detailed evolution of magnetic configuration during the writing and deleting a skyrmion process [54,55]. On the other hand, all simulations in present study were calculated under zero temperature in which the thermal effect is neglected. However, the thermal stability of magnetic kink structure is also a crucial issue for a detailed understanding of their underlying physical properties. For this, we tested the effect of thermal fluctuation on the magnetic structure by using Monte Carlo simulation with ladder cooling protocol [44]. It was found that the magnetic structure may forms at a very low-temperature T ≤ 1.2J FM /k B (k B is the Boltzmann constant), under a magnetic field B z −0.008J FM . Note that temperature T in Monte Carlo simulations is scaled in the unit of J FM /k B . Thus the kink structures are estimated to be stable below a critical temperature of ∼ 40 K, with a typical exchange constant J FM 2.95 meV adopted for Pd/Fe/Ir(111). This suggests that the kink structures in nanostructures can only exist at cryogenic temperatures, in agreement with the latest experimental investigations on Néel skyrmion in Pd/Fe/Ir(111) [56]. In this regard, the previous studies on skyrmions [57][58][59][60] may provide us an enlightened approach to enhance the thermal stability of kinks, which would be an interesting challenge for our further studies.
In summary, we used micromagnetic simulations and Monte Carlo simulations to investigate an exotic Néel-type kink texture in square-shaped nanostructures, which may stably exist in the absence of magnetic field. It was interesting to find that individual kinks hold eight degenerate states as they can reside in one of the four possible corners of a nanostructure, and carry upward or downward polarity. In addition, we proposed some effective schemes to control their dynamics by means of injecting spinpolarization current pulses. It was found that kinks can be created, annihilated, displaced and reversed polarity on demand by applying an in-plane or an out-of-plane spinpolarized current pulse, and are facile switchable among the degenerate kink states. In particularly, the kinks can be switched toward the ferromagnetic-like states and backward reversibly by applying two successive current pulses, indicating the capability of writing and deleting the kink states. These findings predict the existence of Néel-type kinks in the squareshaped nanostructures, as well as provide us a promising approach to tailor the kinks by utilizing the corners of the nanostructures, and control these states by spin-polarized currents. This study also suggests a theoretical guide to explore other chiral magnetic textures in nanostructures of polygon geometries.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.