Writing and storing information in an array of magnetic vortex nanodisks using their azimuthal modes

Writing and storing information in an array of magnetic vortex nanodisks using their azimuthal modes


The switching of a vortex core of a single disk in an array of a multilayer system is investigated by micromagnetic simulation. We found that the perpendicular uniaxial anisotropy (PUA) decreases the frequencies of the azimuthal mode in disks with magnetic vortex configuration. We obtained a phase diagram of magnetic field intensity vs. frequency of the azimuthal mode, as a function of the value of perpendicular uniaxial anisotropy. We demonstrated that rotating magnetic fields (CW and CCW) with frequency equal to azimuthal modes can be used to switch the vortex core of single disks in a disk array. This allows obtaining different memory states with a single array of nanodisks, and therefore writing information through the application of rotating fields.

Magnetic vortex, azimuthal mode, vortex core switching, perpendicular anisotropy, memory states


1 Introduction

The magnetic vortex configuration is characterized by an in plane curling magnetization, and a core, where the magnetization points out of the plane. The curling direction defines the circulation C = +1 (counterclockwise (CCW)) and C = -1 (clockwise (CW)). The core has polarity p = +1 when it points along the +z direction and p = -1 in the -z direction (1).
A magnetic vortex presents a translation mode of low frequency, in the sub-gigahertz range, known as gyrotropic motion (1); (2) and other two modes of higher frequency ( 1GHz): azimuthal and radial modes, which have their origin in magnetostatic interactions and thus are also dependent on the dimensions of the disk (3); (4); (5).
Depending on the ratio of the thickness to the radius of the disk ( = L/R), the azimuthal modes have a splitting in the frequencies (6). These frequencies have CCW and CW senses of rotation (3).
Magnetic vortices have many potential applications in magnetic data storage devices (7); (2); (8); (9). For example, a vortex with p = +1 can store bit 1, and a vortex with p = -1 can store bit 0, or vice-versa. In these applications, the issue of switching vortex cores is a topic of great interest, that has been studied for a long time (10); (11); (12); (13); (14); (15).
Using rotating magnetic fields with a frequency equal to the gyrotropic frequency, it is possible to switch the vortex core polarity (15); (16), but with the downside that this only happens when the sense of rotation of the gyrotropic motion (which is determined by p) coincides with the sense of rotation of the magnetic field (10); (15); (17). Another proposal found in the literature is to use rotating magnetic fields with frequencies equal to the characteristic frequency of the azimuthal modes (11); (14). This method has the great advantage that magnetic rotating fields can be used with both directions of rotation (CCW and CW), and allow shorter switching times (14); (11).
Switching of vortex cores using azimuthal modes allows high switching critical velocities (of the order of 800m/s) (11), compared to those of the gyrotropic mode ( 330m/s) (16).
Nanodisks with a magnetic vortex configuration are generally produced by nanolithography in the form of arrays (matrices of nanodisks) on substrates that can influence their dynamic properties. For example, in a multilayer system, a perpendicular uniaxial anisotropy (PUA) can be induced due to the interface contribution, as has already been demonstrated by Garcia et al. (18). This PUA influences the processes of switching of the vortex core (19); (15).
An array of disks can be used to build an information storage device and/or build logic gate circuits (8); (9). In these arrays, the polarity has an important role, since it determines the type of logical gate to be obtained (8). Consequently, it is necessary to search for mechanisms to control the polarity of a single disk in an array, without altering the polarity of the neighboring disks.
The goal of this work is to propose a novel method for controlling the selectively switching of one single vortex core in a matrix of nanodisk multilayer system, in order to obtain the desired combinations of bits in this matrix. For this purpose, we have used micromagnetic simulation. All simulations were made using the open source software Mumax3 (20), with discretization cell size of 2 2 L nm, where L is the thickness of the disk. The material used is Permalloy (NiFe), with typical parameters (21); (1); (19): saturation magnetization M = 8.6 10 A/m, exchange stiffness = 1.3 10 J/m, and damping constant = 0.01. The perpendicular uniaxial anisotropy constant1 (K) varied from 0 to 200kJ/m. For larger values of K, skyrmion type magnetic structures emerge (15); (19).

2 Results and discussion

2.1 Isolated disk

We used disks with thickness L = 20nm and diameter D = 500nm. For these dimensions, the magnetic vortex configuration is stable (15). We assumed that the vortex core is initially at the equilibrium position at the center of the disk, and has polarity p = +1 and circulation c = +1.
In order to excite the azimuthal spin wave modes, we have applied an in-plane sinc pulse magnetic field B(t) = (B(x)/x,0,0), with x = 2f(t-t), centered on t = 1ns, where B = 1mT is the magnetic field amplitude and f = 50GHz is the frequency of the magnetic field pulse. The frequencies of the modes are obtained by fast Fourier Transform (FFT) from the time evolution of the x-component of the magnetization. These frequencies are shown in Fig. 1. We repeat the same procedure for each value of K.

Figure 1: (a) and (b) Values of the gyrotropic mode frequency and azimuthal spin wave frequencies for m = +1 (counterclockwise) and m = -1 (clockwise) obtained by a fast Fourier transform (FFT) from the time evolution of the x-component of the magnetization for each value of K, for p = +1. Orange triangles correspond to the counterclockwise frequencies, green squares to the clockwise frequencies, and blue stars to the gyrotropic mode.

There are three frequencies for each value of K. The lowest value frequency corresponds to the gyrotropic mode (f 0.35GHz) and the other two frequencies correspond to the m = -1 (clockwise) and m = +1 (counterclockwise) azimuthal modes, respectively (11); (14).
The frequency of the gyrotropic mode remains almost constant with the increase2 of K. The azimuthal frequencies decrease with increasing anisotropy, as shown in Fig. 1(a), because the influence of PUA modifies the configuration of the magnetic vortex (18); (19). It is important to note that the effect of PUA on the vortex configuration is totally different from that produced by a perpendicular magnetic field (PMF). Whereas PUA does not alter the gyrotropic mode, PMF does so in a manner proportional to the intensity of this field (14).
In Fig. 1(b) are shown all the frequencies (f = /2), considering negative values for m = -1 and positive values for m = +1 (14).
In order to switch the vortex core, we have used an in-plane rotating magnetic field ( for CCW and for CW ) bursts with duration of 24 periods, as suggested by Kammerer et al. (11). After the magnetic field is turned off, we have monitored the micromagnetic simulation for an additional 1.5ns with zero magnetic field, to observe possible switching that may occur due to delayed processes (22).
We start by exploring the switching vortex core using f = 0.35GHz (gyrotropic mode); we encountered a minimal magnetic field intensity = 1.2mT, to obtain the switching of vortex cores for the entire range of perpendicular uniaxial anisotropy constant (K) used in this work (see section 1).

Figure 2: Switching times trough the gyrotropic mode as a function of magnetic field intensity (B) for different values of K.

In Fig. 2 are shown the switching times3 (t) as a function of B, for each value of K used in this work, using the gyrotropic mode. We found a decrease of values of t with the increase of B.
For K = 0kJ/m it is obtained a switching time (t) of approximately 15ns (Fig. 2). This time is reduced by approximately 88 with the increase of B from 15ns (B = 1.2mT) to 1.82ns (B = 7mT). For larger intensity magnetic field, undesirable multiple switching events appear. The same behavior is observed for K 0.
Althought t decreases with the increase of B, the critical velocity (approximately 329m/s) that the vortex core reaches before switching, is the same for all values of B. This is known as the universal criterion of switching, as has already been demonstrated by Lee et al. (16). However, this critical velocity decreases when K 0, from 329m/s (K = 0kJ/m) to 200m/s (K = 200kJ/m), but is still independent of B. These values are consistent with those obtained by Fior et al. (15).
In order to obtain the magnetic field intensity to switch the vortex core using azimuthal modes, we have varied B from 1mT to 6mT for m = -1 mode, and from 1mT to 8mT for the m = +1 mode, in 0.2mT steps for both m = +1 (CCW) and m = -1 (CW) modes. The switching phase diagram (B vs. frequency) is shown in Fig. 3. We have slightly varied the values of the frequencies shown in Fig. 1, as suggested in ref. (14), in order to obtain lower values of the threshold magnetic field intensity.

Figure 3: Switching phase diagrams for (a) m = +1 (CCW) and (b) m = -1. Red triangles indicate no switching, green squares indicate single switching and blue squares indicate multiples switching.

We encountered three regions for both modes m = +1 and m = -1: 1) no switching, 2) single switching, and 3) multiple switching. We are interested in the region of single switching, with the purpose of having total control in the selectivity of the resulting polarity p.
In Fig. 3 we do not show the region for 1mT B 2mT for m = +1 mode, and the region for 1mT B 3.6mT for m = -1 mode, since there is no switching for all the values of K used in this work.
The threshold of the magnetic field intensity (B) and the range of single switching of vortex core are different for each value of K (see Fig. 3), and for modes m = +1 and m = -1. A wider range of B values is found for m = -1 mode that result in a single switching (green squares), in comparison with m = +1. For m = +1 mode, multiple switching events are dominant in the phase diagram, whereas for m = -1, single switching events are more frequent.
Multiple switching events appear because the applied magnetic field pumps enough energy to reverse repeatedly the vortex core between p = +1 and p = -1 (13).
For the m = +1 mode, we used a magnetic field (B) of up to 8mT to obtain single switching for K = 200kJ/m. This is different for the m = -1 mode, where a threshold of magnetic field B = 3.6mT is necessary to obtain a single switching of vortex core.
All these differences between modes m = +1 and m = -1 are due to the fact that the modes act differently in the creation of a dip, which is the first step in the switching process (23); (11). Whereas that m = -1 mode leads to the formation of a single dip, m = +1 mode leads to the formation of a double dip (11).
Fig. 4 shows the switching times (t) for both modes, and for the case K = 0kJ/m and K = 100kJ/m. These times decrease with increasing intensity of the magnetic field for m = +1 mode, however, for m = -1, it is observed that for some values of B, t does not have the same behavior4. This can also be attributed to nolinear dynamics.
Similar behaviors were obtained by Kammerer et al. (11) for Permalloy disks, and using K = 0kJ/m.

Figure 4: Switching times versus magnetic field intensity for (a-b) m = +1 (CCW) and (c-d) m = -1 (CW), for K = 0kJ/m and K = 100kJ/m.

It is important to note that although the switching time can be reduced increasing the magnetic field intensity, and using the gyrotropic mode (Fig. 2), we obtain shorter switching times with lower field intensities using the azimuthal modes (Fig. 4). As mentioned earlier, we get t = 1.82ns (B = 7mT) using the gyrotropic mode for K = 0kJ/m, but using the azimuthal mode we get shorter times, of approximately t = 1.29ns for the m = +1 mode, and t = 0.85ns for m = -1 mode. Both values of t were obtained with magnetic field intensity smaller that 7mT (Fig. 4).
The nolinear dynamics breaks the universal criterion for switching of vortex cores found for the gyrotropic mode. We have obtained different values for critical velocities. For example, for K = 0kJ/m, we found an average critical velocity (v) for the entire region where there is single switching (Fig. 3) of approximately 816m/s and 400m/s, for modes m = +1 and m = -1, respectively. These values are higher than those found in the gyrotropic mode, and similar to those found in ref. (11).
In Table 1 are shown the values of the average critical velocities for modes m = +1 and m = -1.

K (kJ/m) v (m/s) v (m/s)
m = +1 m = -1
0 816 400
50 913 380
100 686 381
150 806 491
200 601 375
Table 1: Average critical velocity for different values of K and modes m = +1 and m = -1.

All values of critical velocities shown in Table 1, using the azimuthal modes, are higher than the universal critical velocity found using the gyrotropic mode, of approximately 330m/s. These higher velocities are responsible for the shorter switching times.
The average critical velocities do not show a linear behavior with the increase of K; they increase for K = 50kJ/m (m = +1), then decrease for K = 100kJ/m and then increase again. This behavior contrasts with the case where PMF increases or decreases v depending on whether the PMF is parallel or antiparallel to p (14).
This difference is expected, since PUA does not modify the gyrotropic mode, whereas PMF does. Moreover, PUA and PMF act differently on the vortex core, leading to totally different behaviors in the switching processes, as has already been demonstrated by Fior et al. (15).
Next, we used the influence of PUA and rotating magnetic fields in order to obtain different final states in a 22 matrix of disks.

2.2 2 2 matrix

We will now describe a matrix of four vortex disks, where we can write four bits of information using the rotating magnetic fields.
We used an array of four identical disks, as shown in Fig. 5, with thickness L = 20nm, diameter D = 500nm and separated by an edge to edge distance x = 500nm. Each disk has its own K, with n = 1,2,3,4 and K K K K. The initial configuration is that all disks have polarity and circulation p = C = +1. It is important to mention that the magnetostatic interation between the disks can change the values of the intensity of magnetic fields for switching (24) shown in Fig. 3. However, this does not alter the principle that the magnetic field having frequency equal to either modes (m = +1 or m = -1) only reverses the vortex core of the disk to which these modes correspond.

Figure 5: Array of disks with magnetic vortex configuration and different values of K.
Figure 6: Initial configuration of the matrix of disks and the resulting configurations obtained after applying the rotating fields. Blue number 1 and red number 0 correspond to polarity p = +1 and p = -1, respectively

In order to obtain different final states, as shown in Fig. 6, we used a global CCW rotating magnetic field acting on the entire matrix, with duration of 24 periods.
We used the following convention: blue digit 1 to indicate positive polarity p = +1, and red digit 0 to indicate negative polarity p = -1 (see Fig. 6).
In order to obtain any of the final states that have a single red digit, we choose a frequency f equal to the azimuthal m = +1 mode of the disk of interest to be switched. For example, to switch the vortex core of disk 1 (K = 0kJ/m), we used a frequency f = 8.4GHz (see Fig. 1), and B = 4.6mT. The magnetic field will efficiently excite the disk 1, to which the f frequency corresponds, switching the vortex core only in this disk (see videos in the Supplementary material). It is important to remark that this would be impossible using a frequency equal to the gyrotropic mode, because disks with higher K would switch before those of smaller values of K (15).
A final state with two red digits can be obtained in two steps using the azimuthal modes: first switching one of the disks, and next the second disk. This could be done in one step using the gyrotropic mode, depending on which disks the user wants the switching. For example, if one desires to switch only disks 3 and disk 4, the duration of applied magnetic field will have to be that necessary for the switching to occur on disk 3, since the switching on disk 4 would occur earlier (15), due to K K, thus t t. However, if one desired to switch disk 1 (K) and disk 4 (K), the application of magnetic field will switch the disks with intermediate values of K, such as K and K.
Final states with three red digits can be obtained in three steps, following similar procedure, as in the case of two red digits.
Full switching of all disks is trivial and can be obtained using a frequency equal to that of the gyrotropic mode.5
Switching from negative polarity p = -1 to positive polarity p = +1 is also possible using the azimuthal modes, but changing the sense of rotation of these modes. For p = +1, we have m = +1 (CCW) and m = -1 (CW), and for p = -1 we have m = +1 (CW) and m = -1 (CCW) (11).

3 Conclusion

In this work, we initially studied the influence of PUA on azimuthal modes in a matrix of disks with magnetic vortex configuration, using micromagnetic simulations. Our results show that azimuthal mode frequencies decrease with increase of PUA, and modified the intensity of the magnetic field necessary to switching the vortex core.
Based on this initial study, we then demonstrated that the azimuthal modes can be used to selective switching in arrays of disks, therefore obtaining several different final state configurations using a single array of disks. This shows the great advantage of using the azimuthal modes in comparison to the use of the gyrotropic mode.
This work also shows that when the intrinsic variable K is considered, the universality of the value of the critical velocity is broken, even for the gyrotropic mode.
Our proposal addresses a subject not very studied, the influence of the PUA in the magnetic vortex dynamics, allowing to write information in an array of disks through selective switching of vortex cores. This simple system can be expanded to larger arrays.


The authors would like to thank the support of the Brazilian agencies CNPq and FAPERJ.



  1. These values of K can be obtained experimentally increasing the thickness of disk as shown in ref. (18).
  2. There is a variation of approximately 3 for the maximum value of K, but this change is negligible. This too was demonstrated by Fior et al. (15).
  3. See Supplementary material for details of how these values were obtained these values.
  4. Supplemental material shows t for all values of K used in this work.
  5. The same methodology can be used for the case of having the disks located in the form of nanopillars.


  1. K. Y. Guslienko, Magnetic vortex state stability, reversal and dynamics in restricted geometries, J. Nanoscience Nanotechnol. 8 (2008) 2745–2760. doi:10.1166/jnn.2008.003.
  2. A. P. Guimarães, Principles of Nanomagnetism, 2nd Edition, Springer, Cham, 2017.
  3. K. Y. Guslienko, A. N. Slavin, V. Tiberkevich, S.-K. Kim, Dynamic origin of azimuthal modes splitting in vortex-state magnetic dots, Phys. Rev. Lett. 101 (2008) 247203. doi:10.1103/PhysRevLett.101.247203.
  4. A. A. Awad, K. Y. Guslienko, J. F. Sierra, G. N. Kakazei, V. Metlushko, F. G. Aliev, Precise probing spin wave mode frequencies in the vortex state of circular magnetic dots, Applied Physics Letters 96 (1) (2010) 012503. doi:10.1063/1.3268453.
  5. M. Sproll, M. Noske, H. Bauer, M. Kammerer, A. Gangwar, G. Dieterle, M. Weigand, H. Stoll, G. Woltersdorf, C. H. Back, G. Schütz, Low-amplitude magnetic vortex core reversal by non-linear interaction between azimuthal spin waves and the vortex gyromode, Appl. Phys. Lett. 104 (1) (2014) 012409. doi:10.1063/1.4861779.
  6. J. P. Park, P. A. Crowell, Interactions of spin waves with a magnetic vortex, Phys. Rev. Lett. 95 (2005) 167201. doi:10.1103/PhysRevLett.95.167201.
  7. S. Bohlens, B. Krüger, A. Drews, M. Bolte, G. Meier, D. Pfannkuche, Current controlled random-access memory based on magnetic vortex handedness, Appl. Phys. Lett. 93 (14) (2008) 142508. doi:10.1063/1.2998584.
  8. H. Jung, Y.-S. Choi, K.-S. Lee, D.-S. Han, Y.-S. Yu, M.-Y. Im, P. Fischer, S.-K. Kim, Logic operations based on magnetic-vortex-state networks, ACS Nano 6 (5) (2012) 3712–3717, pMID: 22533663. doi:10.1021/nn3000143.
  9. H. Vigo-Cotrina, A. P. Guimarães, Single array of magnetic vortex disks uses in-plane anisotropy to create different logic gates, Journal of Magnetism and Magnetic Materials 441 (2017) 14 – 20. doi:10.1016/j.jmmm.2017.05.027.
  10. M.-W. Yoo, K.-S. Lee, D.-E. Jeong, S.-K. Kim, Origin, criterion, and mechanism of vortex-core reversals in soft magnetic nanodisks under perpendicular bias fields, Phys. Rev. B 82 (2010) 174437. doi:10.1103/PhysRevB.82.174437.
  11. M. Kammerer, M. Weigand, M. Curcic, M. Noske, M. Sproll, A. Vansteenkiste, B. Van Waeyenberge, H. Stoll, G. Woltersdorf, C. H. Back, G. Schuetz, Magnetic vortex core reversal by excitation of spin waves, Nat. Comm. 2 (279). doi:10.1038/ncomms1277.
  12. M.-W. Yoo, J. Lee, S.-K. Kim, Radial-spin-wave-mode-assisted vortex-core magnetization reversals, Appl. Phys. Lett. 100 (17) (2012) 172413. doi:10.1063/1.4705690.
  13. M. Kammerer, H. Stoll, M. Noske, M. Sproll, M. Weigand, C. Illg, G. Woltersdorf, M. Fähnle, C. Back, G. Schütz, Fast spin-wave-mediated magnetic vortex core reversal, Phys. Rev. B 86 (2012) 134426. doi:10.1103/PhysRevB.86.134426.
  14. M.-W. Yoo, S.-K. Kim, Azimuthal-spin-wave-mode-driven vortex-core reversals, J. Appl. Phys. 117 (2) (2015) 023904. doi:10.1063/1.4905689.
  15. G. B. M. Fior, E. R. P. Novais, J. P. Sinnecker, A. P. Guimarães, F. Garcia, Indirect switching of vortex polarity through magnetic dynamic coupling, J. Appl. Phys. 119 (2016) 093906. doi:10.1063/1.4942534.
  16. K.-S. Lee, S.-K. Kim, Y.-S. Yu, Y.-S. Choi, K. Y. Guslienko, H. Jung, P. Fischer, Universal criterion and phase diagram for switching a magnetic vortex core in soft magnetic nanodots, Phys. Rev. Lett. 101 (2008) 267206. doi:10.1103/PhysRevLett.101.267206.
  17. M. Noske, A. Gangwar, H. Stoll, M. Kammerer, M. Sproll, G. Dieterle, M. Weigand, M. Fähnle, G. Woltersdorf, C. H. Back, G. Schütz, Unidirectional sub-100-ps magnetic vortex core reversal, Phys. Rev. B 90 (2014) 104415. doi:10.1103/PhysRevB.90.104415.
  18. F. Garcia, H. Westfahl, J. Schoenmaker, E. J. Carvalho, A. D. Santos, M. Pojar, A. C. Seabra, R. Belkhou, A. Bendounan, E. R. P. Novais, A. P. Guimarães, Tailoring magnetic vortices in nanostructures, Appl. Phys. Lett. 97 (2) (2010) 022501. doi:10.1063/1.3462305.
  19. E. R. P. Novais, S. Allende, D. Altbir, P. Landeros, F. Garcia, A. P. Guimarães, Effect of perpendicular uniaxial anisotropy on the annihilation fields of magnetic vortices, J. Appl. Phys. 114 (15) (2013) –. doi:10.1063/1.4824803.
  20. A. Vansteenkiste, J. Leliaert, M. Dvornik, M. Helsen, F. Garcia-Sanchez, B. Van Waeyenberge, The design and verification of Mumax3, AIP Advances 4 (2014) 107133. doi:10.1063/1.4899186.
  21. K. Y. Guslienko, B. A. Ivanov, V. Novosad, Y. Otani, H. Shima, K. Fukamichi, Eigenfrequencies of vortex state excitations in magnetic submicron-size disks, J. Appl. Phys. 91 (2002) 8037–8039. doi:10.1063/1.1450816.
  22. M. Kammerer, M. Sproll, H. Stoll, M. Noske, M. Weigand, C. Illg, M. Fähnle, G. Schütz, Delayed magnetic vortex core reversal, Applied Physics Letters 102 (1) (2013) 012404. doi:10.1063/1.4773592.
  23. K. Y. Guslienko, K.-S. Lee, S.-K. Kim, Dynamic origin of vortex core switching in soft magnetic nanodots, Phys. Rev. Lett. 100 (2008) 027203. doi:10.1103/PhysRevLett.100.027203.
  24. Y. Lu, Z. Zhang, Y. Liu, Magnetic interaction effect on the critical switching current in vortex arrays, J. Appl. Phys. 109 (10) (2011) 103906. doi:10.1063/1.3590333.
This is a comment super asjknd jkasnjk adsnkj
The feedback cannot be empty
Comments 0
The feedback cannot be empty
Add comment

You’re adding your first comment!
How to quickly get a good reply:
  • Offer a constructive comment on the author work.
  • Add helpful links to code implementation or project page.