Abstract
The collective dynamics in populations of magnetic spin torque oscillators (STO) is an intensely studied topic in modern magnetism. Here, we show that arrays of STO coupled via dipolar fields can be modeled using a variant of the Kuramoto model, a wellknown mathematical model in nonlinear dynamics. By investigating the collective dynamics in arrays of STO we find that the synchronization in such systems is a finite size effect and show that the critical coupling—for a complete synchronized state—scales with the number of oscillators. Using realistic values of the dipolar coupling strength between STO we show that this imposes an upper limit for the maximum number of oscillators that can be synchronized. Further, we show that the lack of long range order is associated with the formation of topological defects in the phase field similar to the twodimensional XY model of ferromagnetism. Our results shed new light on the synchronization of STO, where controlling the mutual synchronization of several oscillators is considered crucial for applications.
Introduction
The emergence of coherent phases of interacting oscillators is at the foundation of the cooperative functioning of a wealth of different systems in nature^{1}. Examples of collective behavior can be chosen within a wide range of systems such as laser arrays^{2}, Josephson junctions^{3}, chemical reactions^{4}, synchronously flashing firefly populations^{1}, disease spreading^{5}, or cortical oscillations in the brain^{6,7}. Science has sought mathematical models for understanding collective phenomena in large populations of oscillators that were tractable both analytically and numerically.
The Kuramoto model is a well known mathematical model in nonlinear dynamics that describes large systems of coupled phase oscillators^{8}. The model, with a remarkable simplicity, has been used to describe the essential features of collective excitations in a vast set of biological and physical phenomena^{8,9,10,11,12,13,14,15,16}. Although the Kuramoto model originally described oscillators interacting alltoall with the same strength, variations of the model have been used to describe systems with phase offset and time delays in the couplings, other topologies like onedimensional structures with local couplings etc. (see e.g. ref. 9 for an overview of extensions of the Kuramoto model). In particular, twodimensional Kuramoto networks with diffusive local coupling accept solutions consisting in waves, spirals and many other patterns^{17}.
Understanding the collective behavior in oscillator networks is also an intensely studied topic in modern magnetism: the synchronization of spin torque oscillators (STO). STO are strongly nonlinear magnetic oscillators that can be implemented into nanoscale devices working at microwave frequencies and can be frequency and phase locked to external oscillatory signals or other STO^{18,19,20,21,22,23,24,25,26,27,28,29,30,31,32}. They are envisaged to be useful for a variety of advanced magnetic nanodevices, as microwave sources and for signal processing in telecommunication technologies (see e.g. refs 33, 34, 35 and references therein). STO have also been proposed as possible candidates for a full spintronic implementations of neural networks, based on nanodevices emulating both neurons and synapses^{35,36}. Building artificial neural networks for computation is an emerging field of research within bioinspired computing^{33,34,35,36,37,38,39,40}, where controlling the collective behavior in oscillator networks is crucial.
In both experimental and theoretical studies, most of the work has been performed for limited number of oscillators. Experimentally, the synchronization of STO has proven to be difficult and the synchronization of only a few oscillators has been demonstrated^{30,32}. Theoretically, the magnetization dynamics of STO is modeled with the LandauLifshitzGilbertSlonzewski (LLGS) equation^{41,42}, but large number of STO lead to challenging computations caused by the nonlocal dipolar fields. It is important to consider that in these nonlinear systems “more is different” and that the collective behavior can not be derived simply from the behavior of its individual elements. Thus, a theoretical framework capable to capture the essential dynamics would be ideal to explore those systems.
Here, we show that twodimensional arrays of STO coupled via dipolar fields can be modeled by a variant of the Kuramoto model. We begin with describing two coupled STO with the Thiele equation^{43} and show that for smallamplitude oscillations the system can be described as a simple phase oscillator model. Next, we model the interactions for the case of a twodimensional array of oscillators based on the dipolar coupling and obtain a modified Kuramoto model. Finally we compare the results from our model to the micromagnetic solution of the LLGS equation.
We find that the synchronization in twodimensional arrays of dipolar coupled STO is purely a finite size effect and the critical coupling strength for obtaining a globally synchronized state scales with the number of oscillators N as λ_{crit} ∝ log(N). Using realistic values of the dipolar coupling strength between STO we show that this imposes an upper limit for the maximum number of STO that can be synchronized. Further, we study the synchronization transition between the initial formation of locally synchronized clusters and the globally synchronized and phase coherent state and correlate it with a transition in the local order of the system. We also observe the emergence of topological defects and the formation of patterns in the phase field similar to the twodimensional XYmodel of magnetism—suggesting a connection between arrays of STO, systems described by a 2d Kuramoto model and the 2d XY model of statistical mechanics.
Results
From the Thiele equation to the Kuramoto model
We are considering STO whose free layer ground state configuration is a magnetic vortex. The vortex state is characterized by inplane curling magnetization and a small (~10 nm) region of the vortex core with outofplane magnetization^{44}. The gyrotropic motion of the vortex core is driven by the injection of a DC spin polarized current through the STO stack and can be described by a gyration radius r and phase θ, as illustrated in Fig. 1a. We first model the interaction of two vortices with the Thiele equation with an extra term that accounts for the vortex interaction. These equations describe the vortices motion given by their coordinates X_{1,2} in their self induced gyrotropic mode and include the spintransfertorque (STT) as well as a coupling term^{43,45}.
Here, G is the gyroconstant, k(X_{1,2}) the confining force, D_{1,2} the damping coefficient and F_{STT} the STT. The interaction between the neighboring STO illustrated in Fig. 1b is summarized by a dipolar coupling term given by F_{int} = −μ(d)X_{2,1}, where μ(d) describes the interaction strength as a function of the separation d between the two STO.
Assuming a small difference in the nominal frequencies of two coupled STO described by Eq. (1), one can linearize the set of equations following the approach by Belanovsky et al.^{27}, showing that the dynamics of the phase difference between the STO can be described by Adler’s equation^{46}. Following these approximations, the set of equations reduce to that of two coupled phase oscillators θ_{1} and θ_{2}: (see ‘Supplementary information’ for details).
where ω_{1,2} are eigenfrequencies of oscillators θ_{1,2} respectively and λ describes the interaction strength trying to synchronize them.
To check the validity of the approximations, one can compare the results obtained using a simplified phase oscillator model to a numerical solution of Eq. (1) as well as a micromagnetic solution of the full system using the LLGS equation. This was done by Belanovsky et al.^{27}, where they found that for a small difference in nominal frequencies, in their case given by a difference in STO disc diameter ΔD/D_{0} ≤ 5%, the synchronization can be qualitatively described using the simplified model. Assuming the error in stateoftheart fabrication processes is below this limit, the simplified equations are a valid description of the system.
The functional form of Eqs (2, 3) is the same as that of the well known Kuramoto model^{8,9}, which is a generalization for the case of an ensemble of weakly coupled phase oscillators. Considering the interaction between several STO, we obtain a Kuramoto model where the single oscillator state is described through the dynamic equation of its phase θ_{i} due to the interaction with its surrounding oscillators θ_{j}:
The coupling term is here generalized to include the interaction between several oscillators, determined by the interaction strength λ_{ij} between oscillators θ_{i} and θ_{j}. This determines the nature of the interaction, ranging from a global alltoall coupling where λ_{ij} = λ for all oscillators, to a local interaction where λ_{ij} = 0 for all but the nearest neighbors. Here, we are considering the intermediate case of a nonlocal coupling to mimic the dipolar interaction between neighboring STO. Starting from a macrodipole approximation for the dipolar energy between two magnetic dipoles μ_{1} and μ_{2}, the average interaction strength is found to decay as ^{47}, where d_{ij} is the distance between oscillators θ_{i} and θ_{j}. We thus set the coupling strength to
where we include interactions within a coupling radius R, indicated by the blue circle in Fig. 1c. The network model for the STO array is implemented with biperiodic boundary conditions and the time evolution of the oscillator phases given by Eq. (4) is solved numerically. A small random disorder in the oscillator eigenfrequencies is included by setting ω_{i} = ω_{0} ± δω_{i}. Here, ω_{0} = 1 GHz and δω_{i} represents a uniformly distributed random disorder where δω_{i}/ω_{0} ≤ 2.5%. The interaction strength is determined by the STO spacing and size, as well as the magnetic material properties^{31}. Here, the interaction strength has been varied in the range λ = 1–20 MHz and is in the same range as the interaction strength extracted from micromagnetic simulations for similar STO^{31}.
In order to evaluate the Kuramoto model as a valid description for arrays of STO, we compare it to a micromagnetic solution of the complete system, accounting for all dynamic dipolar terms (see ‘Methods’ section). To compare the Kuramoto model and the micromagnetic solution, we define a suitable order parameter to distinguish disordered and synchronized states. The phase of the individual oscillators θ_{i} is used to define the order parameter ρ, describing the phase coherence in a system of N oscillators:
The case ρ = 0 corresponds to the maximally disordered state, whereas ρ = 1 represents the state where all oscillators are perfectly synchronized and phase coherent. In addition to the global order parameter ρ, we define a local correlation function β:
where the brackets indicate a summation over neighboring oscillators and n is the number of neighbors. β_{i} is a measure of the phase correlation of oscillator θ_{i} and its neighbors, indicated by the blue and red oscillators in Fig. 1c respectively. If oscillator θ_{i} is located within a synchronized cluster, β_{i} → 1. Calculating β thus allows for investigating the formation of locally synchronized clusters and the emergence of patterns of synchronized states, which can not be obtained simply from the global order parameter ρ.
Kuramoto model versus micromagnetic simulations
We now compare the Kuramoto model given by Eq. (4) and the full micromagnetic solution of the LLGS equation (see ‘Methods’ section for details). Starting from a disordered initial state, we investigate the synchronization dynamics by calculating the time evolution of the phase distribution θ_{i} and local correlation β_{i}. As an example, we show in Fig. 2 snapshots of θ_{i} and β_{i} for the Kuramoto model and the micromagnetic solution for a system of 45 × 45 oscillators at times t_{1} and t_{2} > t_{1}, with random initial phases. At time t_{1}, one notices the initial formation of small locally synchronized clusters, as seen through both the phase distribution and the bright areas in the correlation maps, where β → 1. As time progress to t_{2}, these clusters grow in size and merge with neighboring clusters.
Comparing the two models, we find that they both show the same behavior. For weak interaction strengths, the system tends to be in a disordered state with no correlation between neighboring oscillators. By increasing the interaction strength above a certain threshold, synchronized clusters begin to form. The oscillators within each cluster are synchronized, but might not be phase coherent with other clusters. This can be seen in the phase maps in Fig. 2, where the individual clusters have different phases. As time progresses, the transition from a disordered to a synchronized state is governed by the growth and merging of neighboring clusters, reaching a globally synchronized and phase coherent state for sufficiently strong interactions.
Depending on the interaction strength, the system ends up in either a disordered, partially synchronized or globally synchronized state. Controlling the interaction strength is thus the key parameter to determine the system behavior. The interaction strength needed to obtain synchronization will depend on the differences in the nominal frequencies of the oscillators^{31}. Another important consideration, is whether the critical interaction strength also depends on the number of oscillators. Lee et al.^{48} have studied the synchronization in a 2d Kuramoto model with a nearest neighbor interaction. They showed that the transition to a synchronized state depends strongly on system size and that the critical coupling strength needed to synchronize scales with the number of oscillators N as λ_{crit} ∝ log(N). This raises the question if such a scaling law can also be observed in our models: observing the same scaling laws in both the Kuramoto model and the micromagnetic solution would strengthen the suggestion of the Kuramoto model as a valid description of arrays of STO.
We first consider our Kuramoto model, which has a nonlocal interaction to mimic the dipolar interaction in arrays of STO. Due to the increased complexity compared to the nearest neighbor model studied by Lee et al.^{48}, an analytical derivation of the scaling behavior with system size is to our knowledge still an open question. To investigate the scaling behavior we thus resort to a numerical solution. We performed simulations with the number of oscillators ranging from N = 9 to N = 2500, gradually increasing the interaction strength between each simulation until the system reaches a synchronized state at a critical coupling strength, λ_{crit}. 100 simulations were performed for each system size, with different initial oscillator phases and eigenfrequencies. In Fig. 3a we show a plot of λ_{crit} vs. number of oscillators, N. The results indicate that the critical coupling strength scales as λ_{crit} ∝ log(N), same as the nearest neighbor Kuramoto model investigated by Lee et al.^{48}. The main difference compared with the nearest neighbor model is that we include interactions within a coupling radius R, as indicated by the blue circle in Fig. 1c. The imposed cutoff radius has a physical justification when considering a realistic system, which would inevitably include thermal noise. As the STO we consider are weakly coupled and the dipolar interaction decay with distance, there will be a limiting spacing where the thermal noise level is comparable to the coupling strength. In our model we observe that the results do not depend qualitatively on the range of the cutoff radius and that the large scale behavior is dominated by the diffusive coupling. This suggests that insight from studying the analytically tractable nearest neighbor model of Lee et al.^{48} might provide valuable insight into the behavior of STO arrays.
The coupling in the Kuramoto model is defined simply as an interaction strength given by λ in Eq. (4). In the micromagnetic model the coupling comes from dipolar interactions, determined by the magnetic material properties and the STO spacing. In order to investigate the scaling behavior in the micromagnetic solution, one thus needs to relate the critical interaction strength to a critical STO spacing. Following the aforementioned macrodipole approximation that the effective dipolar interaction decay as 1/d^{3}, one can relate the coupling strength λ_{crit} to a critical spacing d_{crit} between neighboring STO as d_{crit} ∝ [log(N)]^{−1/3}. Micromagnetic simulations were then performed to obtain d_{crit} vs. number of oscillators in the range N = 3 × 3 to N = 15 × 15. The results are shown as the red datapoints in Fig. 3b. For comparison, we show as a solid line the expected scaling from the Kuramoto model: d_{crit} ∝ [log(N)]^{−1/3}. The good agreement between the micromagnetics result and the Kuramoto model indicate that they both follow the same scaling law, strengthening the suggestion of the Kuramoto model as a valid description for arrays of STO.
Scaling of output power in arrays of STO
To investigate the implications of the scaling with system size, we consider a model calculation of the output power as we increase the number of STO. For applications of STO as e.g. nanoscale microwave generators, the power output of a single STO is not competitive. Decisive improvement is expected from the synchronization and phase locking of several STO, as this would result in a quadratic increase of the output power, P ∝ N^{2} for N synchronized oscillators.
The output of a single STO can be described by its amplitude and phase, . In our model we have assumed a constant amplitude for all STO and the total amplitude A for an array of STO is given by: . The power output is proportional to A^{2} and for N oscillators we obtain , from the definition of the order parameter ρ in Eq. (6). A quadratic scaling in the power output, P ∝ N^{2}, implies a perfectly synchronized and phase coherent state, given by ρ = 1. However, as the coupling strength to obtain a synchronized state scales with the number of STO, this will affect the power output when scaling up to large arrays.
As an example, we consider a system of STO composed of 200 nm diameter spin valve nanopillars with 15 nm thick Permalloy as the ferromagnetic layer. The average interaction energy can be extracted from micromagnetic simulations and for an edgeedge spacing of 150 nm the interaction strength is found to be 7–8 MHz^{31}. To increase the output power, we are now interested in scaling up to a large number of STO. We keep the same STO spacing when scaling up, e.g. keeping λ fixed. From Fig. 3a, we see that for an interaction strength , there is a corresponding number of STO, , where . This can be illustrated by calculating the total power output as the number of STO is increased, as shown in Fig. 3c. For a small number of STO, we see that the power output follows close to the ideal N^{2} scaling. However, as the number of STO is increased, the scaling with system size becomes increasingly important. For a certain number of STO, indicated by in Fig. 3a, the interaction is no longer strong enough to obtain a synchronized state. This is also illustrated in the inset of Fig. 3c, where we plot the order parameter ρ for arrays of 3 × 3 and 13 × 13 STO respectively, showing the transition from a synchronized state to chaotic behavior as the number of STO is increased above .
This illustrates the importance of our findings that synchronization in such 2d arrays is purely a finite size effect. The interaction strength is limited by the material properties and STO spacing and using realistic values of the coupling strength we start to see significant deviations from the ideal P ∝ N^{2} scaling for array sizes larger than 10 × 10 STO (see Fig. 3c). This means that in a physical realizable system, the scaling with system size imposes an upper limit for the maximum number of STO that can be synchronized.
Topological defects
That the synchronized and phase coherent state is purely a finite size effect, is similar to that of the classical 2d XY model of magnetism. The Kuramoto model is indeed similar to the 2d XY model^{49}, where the direction of spin in the XY model corresponds to the oscillator phase in the Kuramoto model. In the 2d XY model a long range ordered phase is absent due to the presence of spin wave fluctuations and topological defects. The lack of long range order is a specific case of the MerminWagner theorem in spin systems^{50}, stating that continuous symmetries cannot be spontaneously broken in systems with sufficiently short range interactions in dimensions d ≤ 2. The fluctuations preventing long range order in the 2d XY model diverge logarithmically with system size^{49}, in agreement with the logarithmic scaling observed in our system of STO.
Similar topological defects at the boundaries between locally synchronized clusters have previously been observed in the nearest neighbor Kuramoto model^{48} as well as in other twodimensional oscillator network models^{17,51,52,53}. In oscillator networks this is associated with the appearance of topological defects in the oscillator phase field, θ_{i}. In the continuum limit this is expressed as:
where dl is an integration path enclosing the defect and n is the topological charge.
Such topological features are observed also in our Kuramoto model for arrays of STO. The presence of vortices in the phase field is more pronounced as the system size increases. As an example we here consider an array of 50 × 50 oscillators. The disorder in the system is kept constant (given by the difference in the nominal frequencies of the oscillators) and the interaction strength λ is varied, acting as the inverse temperature: as coupling increases, the system becomes more ordered. Starting from a disordered initial state and varying the interaction strength between each simulation, we observe 4 different regimes: For weak interaction strengths we observe the formation of locally synchronized clusters, where cluster sizes increase with interaction strength. Apart from the localized clusters there is no long range order in the system (indicated as regime 1 in Fig. 4). Increasing λ above a certain threshold, we enter regime 2. Here we observe the formation of vortices in the phase field and as an example we show in Fig. 4c a state with 4 vortices. The topological charge is conserved in the system and two vortices of charge ±1 respectively is present.
In both regime 1 and 2 long range order in the system is absent and ρ ≈ 0, as indicated in Fig. 4a. (that ρ > 0 here is a result of fluctuations due to the finite array size). Increasing λ further we enter regime 3, where the transition from regime 2 → 3 is governed by vortex annihilation processes (see ‘Supplementary information’). Here there are no topological defects in the phase field and the lack of global phase coherence is due to spin waves in the phase field where the oscillator phases change smoothly across the array (regime 3 in Fig. 4b). For sufficiently strong interactions we enter regime 4. Increasing the interaction strength is analogous to increasing the exchange coupling in a Ferromagnetic system, resulting in a more ordered state. The result here is a gradual suppression of the spin waves in the phase field observed in regime 3 as the interaction strength is increased. Regime 4 is thus characterized as the globally phase coherent state where all oscillators are synchronized and phase coherent (regime 4 in Fig. 4b).
The growth of the order parameter ρ with increasing coupling strength λ in Fig. 4a resembles that of a phase transition. Previous work have shown that the synchronization transition in the globally coupled Kuramoto model can be described as a phase transition, where the nature of the transition can be of first or secondorder depending on the frequency distribution and coupling topology^{8,54}. The Kuramoto model with finite range couplings is less studied, as these systems are difficult to analyze and solve analytically. A study of the locally coupled Kuramoto model on a ddimensional lattice have shown that the synchronization transition depends strongly on the lattice dimensionality and indicates d = 4 as the lower critical dimension for phase synchronization^{55}. This is in agreement with the observed scaling with system size in our model, which indicates that the synchronization transition is purely a finite size effect.
In order to investigate the synchronization transition in our model further, we calculate the spatial correlation function for the array of oscillators. The correlations decay with distance and asymptotically the correlation function is given by: 〈θ(r) ⋅ θ(R)〉 ∝ e^{−r−R/ξ}/r − R^{η}. This describes the correlation between oscillators at positions r and R respectively and the correlation length ξ is obtained by averaging over all positions r and R in the array (an example is shown in ‘Supplementary information’ Fig. S2). From the decay of the correlation function, we obtain the correlation length ξ as a function of the interaction strength λ. Conventional phase transitions are accompanied by a diverging correlation length close to the transition. Here, we do not observe a diverging ξ going from the disordered to the phase coherent state (1 → 4 in Fig. 4) and the correlation length remains finite. As inset in Fig. 4a we show a loglog plot of the correlation length normalized to the system size, ξ/L, (N = L × L). The results indicate a power law relating the correlation length and the interaction strength as ξ ∝ λ^{ν}, where the exponent for this case was found to be ν = 2.1 ± 0.1. This means that the correlation length simply scales with the coupling strength and the transition between regimes 1–4 in Fig. 4 correspond to structures of ever increasing length scales. The transition to the phase coherent state (ρ → 1) occurs when the correlation length approaches the system size L, underlining the finite size effects on the synchronization transition and that the system is not undergoing a conventional phase transition. Further investigations of finite size effects, the lack of long ranger order in the Kuramoto model and the connection to the 2d XY model will be the subject of future work.
Discussion
To summarize, we have shown that the Kuramoto model provides a good description of arrays of STO. It provides a simple theoretical model to study large populations of coupled STO, which were previously unaccessible due to the long computation time for a full micromagnetic solution. By investigating the collective dynamics in large arrays of STO, we observed a scaling with system size indicating that the synchronization in arrays of dipolarcoupled STO is purely a finite size effect. The critical coupling strength to obtain a globally synchronized state scales with the number of oscillators, as λ_{crit} ∝ log(N), preventing global synchronization for large system sizes. As a consequence of the scaling with system size, we showed that for realistic values of the dipolar coupling strength between STO this imposes an upper limit for the maximum number of oscillators that can be synchronized. Further, we showed that the lack of long range order and scaling with system size is associated with the emergence of topological defects and the formation of patterns in the phase field, similar to that of the 2d XYmodel of magnetism.
In the present study we considered dipolarcoupled STO, where the short time delay in the coupling between neighboring oscillators compared to the oscillator frequency means that phase delay in the couplings can be neglected. However, for other coupling mechanisms, time delay can become significant. Interaction mediated by spin waves provide a different mechanism to obtain synchronization of STO^{32}, where the finite propagation speed of the spin waves results in a phase offset in the couplings. Another recent proposal includes the use of nonlocal electrical couplings, where the coupling phase can be externally tuned through an electrical delay line^{56}.
From a dynamical systems point of view, the study of time delay induced modifications to the couplings is of fundamental interest, as well as of practical relevance for modeling of physical, biological and chemical systems. In such systems, time delay is associated with finite propagation velocity of the couplings via e.g latency times of neuronal excitations, reaction times in chemical systems etc. (see e.g. ref. 9 and references therein). The possibility of designing STO arrays with a defined phase offset in the couplings suggests a real world analog to the more general SakaguchiKuramoto model on a 2d lattice, which allows for a phase lag in the couplings^{57}.
Our study suggests, on the one hand, that the use of models from nonlinear dynamics can be useful for describing synchronization of magnetic oscillators and, on the other hand, arrays of STO as a a physical realizable model system for the Kuramoto model on a 2d lattice.
Methods
Micromagnetic model
The micromagnetic model is defined as arrays of discs, where the system is divided into a grid with a mesh size of 5 nm, limited by the exchange length of the ferromagnetic material (here Py). The volumetric quantities such as the magnetization M and effective field H_{eff} are treated at the center of each cell, whereas coupling quantities such as exchange strength are considered at the faces between cells. The numerical solution was obtained using the micromagnetic solver Mumax^{3} ^{58}, which uses a RKF 45 method to solve the LandauLifshitzGilbertSlonzewski (LLGS) equation^{41,42} given by:
Here, γ is the gyromagnetic ratio, α the damping parameter and M_{s} the saturation magnetization. The spintransfer torque term is given by , the charge current density J and the free layer thickness d. P(θ) is a polarization function assumed to increase with the relative angle θ between the magnetization of the free layer and the fixed layer and m_{f} is a unit vector in the direction of the magnetization of the free layer.
In the model, each disc is composed by a magnetic free layer and a fixed polarizer which generates a perpendicular spin polarization p_{z}. The free layer in the STO is 30 nm thick Py with a disc diameter of 150 nm and the damping parameter α was set to 0.01. A small disorder in the eigenfrequencies of the individual STO is included through a random distribution in the saturation magnetization in the range [865, 885] ⋅ 10^{3} A/m, resulting in a slight variation of STO eigenfrequencies. All discs were initialized with a vortex of same polarity and chirality and the centercenter spacing of the discs was varied to change the interaction strength. The polarizing layers are not included in the model as these layers, being uniformly magnetized in z direction have almost no influence on the vortices motion. The vortex gyration was driven by a DC spin current with a polarization p_{z} = 0.3 and current density J ≈ 4.3 ⋅ 10^{7} A/m^{2}. During the simulations, a static magnetic field of 150 mT was applied along the z direction to set the vortex core polarity.
Additional Information
How to cite this article: Flovik, V. et al. Describing synchronization and topological excitations in arrays of magnetic spin torque oscillators through the Kuramoto model. Sci. Rep. 6, 32528; doi: 10.1038/srep32528 (2016).
References
Strogatz, S. H. & Stewart, I. Coupled oscillators and biological synchronization. Sci. Am. 269(6), 102–109 (1993).
Kourtchatov, S. Yu., Likhanski, V. V., Napartovich, A. P., Arecchi, F. T. & Lapucci, A. Theory of phase locking of globally coupled laser arrays. Phys. Rev. A 52, 4089 (1995).
Wiesenfeld, K., Colet, P. & Strogatz, S. H. Synchronization Transitions in a Disordered Josephson Series Array. Phys. Rev. Lett. 76, 404 (1996).
Singh, R. & Sinha, S. Spatiotemporal order, disorder and propagating defects in homogeneous system of relaxation oscillators. Phys. Rev. E. 87, 012907 (2013).
Kuperman, M. & Abramson, G. Small World Effect in an Epidemiological Mode. Phys. Rev. Lett. 86, 2909 (2001).
Breakspear, M., Heitmann, S. & Daffertshofer, A. Generative models of cortical oscillations: neurobiological implications of the Kuramoto model. Front. Hum. Neurosci. (2010).
Izhikevich, E. M. Which Model to Use for Cortical Spiking Neurons? IEEE transactions on neural networks. 15.5, 1063–1070 (2004).
Kuramoto, Y. & Nishikawa, I. Statistical macrodynamics of large dynamical systems. Case of a phase transition in oscillator communities. J. Stat. Phys. 49, 569 (1987).
Acebrón, J. A. et al. The Kuramoto model: A simple paradigm for synchronization phenomena. Rev. Mod. Phys. 77, 137 (2005).
Filatrella, G., Nielsen, A. H. & Pedersen, N. F. Analysis of a power grid using a Kuramotolike model. The European Physical Journal B. 61(4), 485–491 (2008).
Daniels, B. C., Dissanayake, S. T. M. & Trees, B. R. Synchronization of coupled rotators: Josephson junction ladders and the locally coupled Kuramoto model. Phys. Rev. E 67, 026216 (2003).
Wiesenfeld, K., Colet, P. & Strogatz, S. H. Frequency locking in Josephson arrays: Connection with the Kuramoto model. Phys. Rev. E 57, 1563 (1998).
Néda, Z., Ravasz, E., Vicsek, T., Brechet, T. & Barabási, A. L. Physics of the rhythmic applause. Phys. Rev. E 61, 6987 (2000).
Heinrich, G., Ludwig, M., Qian, J., Kubala, B. & Marquardt, F. Collective Dynamics in Optomechanical Arrays. Phys. Rev. Lett. 107, 043603 (2011).
Cumin, D. & Unsworth, C. P. Generalising the Kuramoto model for the study of neuronal synchronisation in the brain. Physica D: Nonlinear Phenomena, 226(2), 181–196 (2007).
Timms, L. & English, L. Q. Synchronization in phasecoupled Kuramoto oscillator networks with axonal delay and synaptic plasticity. Phys. Rev. E. 89, 032906 (2014).
Lauter, R., Brendel, C., Habraken, S. J. M. & Marquardt, F. Pattern phase diagram for twodimensional arrays of coupled limitcycle oscillators. Phys. Rev. E. 012902 (2015).
Rippard, W. H. et al. Injection Locking and Phase Control of Spin Transfer Nanooscillators. Phys. Rev. Lett. 95 067203 (2005).
Kaka, S. et al. Mutual phaselocking of microwave spin torque nanooscillators. Nature, 437 389–392 (2005).
Mancoff, F. B., Rizzo, N. D., Engel, B. N. & Tehrani, S. Phaselocking in doublepointcontact spintransfer devices. Nature 437 393–395 (2005).
Pufall, M. R., Rippard, W. H., Russek, S. E., Kaka, S. & Katine, J. A. Electrical Measurement of SpinWave Interactions of Proximate Spin Transfer Nanooscillators. Phys. Rev. Lett. 97, 087206 (2006).
Bonin, R., Bertotti, G., Serpico, C., Mayergoyz, I. D. & d’Aquino, M. Analytical treatment of synchronization of spintorque oscillators by microwave magnetic fields. Eur. Phys. J. B. 68, 221 (2009).
Zhou, Y., Person, J. & Åkerman, J. Intrinsic phase shift between a spin torque oscillator and an alternating current. J. Appl. Phys. 101, 09A510 (2007).
Recende, S. M., de Aguiar, F. M. & Azevedo, A. SpinWave Theory for the Dynamics Induced by Direct Currents in Magnetic Multilayers. Phys. Rev. Lett. 94, 037202. (2005).
Rezende, S. M., de Aguiar, F. M. & Azevedo, A. SpinWave Theory for the Dynamics Induced by Direct Currents in Magnetic Multilayers. Phys. Rev. B. 73, 094402 (2006).
Slavin, A. N. & Tiberkevich, V. S. Theory of mutual phase locking of spintorque nanosized oscillators. Phys. Rev. B. 74, 104401 (2006).
Belanovsky, A. D. et al. Numerical and analytical investigation of the synchronization of dipolarly coupled vortex spintorque nanooscillators. Appl. Phys. Lett. 103, 122405 (2013).
Belanovsky, A. D. et al. Phase locking dynamics of dipolarly coupled vortexbased spin transfer oscillators. Phys. Rev. B 100409(R) (2012).
Georges, B., Grollier, J., Cros, V. & Fert, A. Impact of the electrical connection of spin transfer nanooscillators on their synchronization: an analytical study. Appl. Phys. Lett 92, 232504 (2008).
Ruotolo, A. et al. Phaselocking of magnetic vortices mediated by antivortices. Nature Nanotechnology 4, 528–532 (2009).
Locatelli, N. et al. Efficient Synchronization of Dipolarly Coupled VortexBased Spin Transfer NanoOscillators. Sci. Rep. 5, 17039 (2015).
Housang, A. et al. Spinwavebeam driven synchronization of nanocontact spintorque oscillators. Nature Nanotechnology 11, 280–286 (2016).
Bonetti, S. & Åkerman, J. Magnonics: From Fundamentals to Applications (2013).
Macià, F., Kent, A. D. & Hoppensteadt, F. C. Spinwave interference patterns created by spintorque nanooscillators for memory and computation. Nanotechnology 22(9), 095301 (2011).
Locatelli, N., Cros, V. & Grollier, J. Spintorque building blocks. Nature Materials, 13, 11 (2014).
Csaba, G. & Porod, W. Computational Study of SpinTorque Oscillator Interactions for NonBoolean Computing Applications. IEEE Trans. Magn. 49(7), 4447–4451 (2013).
Hoppensteadt, F. C. & Izhikevich, E. M. Oscillatory Neurocomputers with Dynamic Connectivity. Phys. Rev. Lett. 82, 2983 (1999).
Maffezzon, P., Bahr, B., Zheng, Z. & Daniel, L. Oscillator Array Models for Associative Memory and Pattern Recognition. IEEE Trans. Circuits Syst. I, Reg. Papers, 62(6), 1591–1598 (2015).
Locatelli, N. et al. Spintronic devices as key elements for energyefficient neuroinspired architectures. Proceedings of the 2015 Design, Automation & Test in Europe Conference & Exhibition.
Vincent, A. F. et al. Spintransfer torque magnetic memory as a stochastic memristive synapse for neuromorphic systems. IEEE Transactions on Biomedical Circuits and Systems. 9(2) (2015).
Landau, L. & Lifshitz, E. On the Theory of the Dispersion of Magnetic Permeability in Ferromagnetic Bodies. Phys. Z. Sowjetunion 8, 153 (1935).
Slonczewski, J. C. Currentdriven excitation of magnetic multilayers. J. Magn. Mater, 159, 1–2, L1–L7 (1996).
Thiele, A. A. SteadyState Motion of Magnetic Domains. Phys. Rev. Lett. 30, 230 (1973).
Shinjo, T., Okuno, T., Hassdorf, R., Shigeto, K. & Ono, T. Magnetic Vortex Core Observation in Circular Dots of Permalloy. Science, 289, 930–932 (2000).
Khvalkovskiy, A. V., Grollier, J., Dussaux, A., Zvezdin, K. A. & Cros, V. Vortex oscillations induced by spinpolarized current in a magnetic nanopillar: Analytical versus micromagnetic calculations. Phys. Rev. B 80, 140401(R) (2009).
Adler, R. A Study of Locking Phenomena in Oscillators. Proc. IRE 34, 351 (1946).
Araujo, F. A. et al. Optimizing magnetodipolar interactions for synchronizing vortex based spintorque nanooscillators. Phys. Rev. B 92, 045419 (2015).
Lee, T. E., Tam, H., Refael, G., Rogers, J. L. & Cross, M. C. Vortices and the entrainment transition in the twodimensional Kuramoto model. Phys. Rev. E 82, 036202 (2010).
Kosterlitz, J. M. & Touless, D. J. Ordering, metastability and phase transitions in twodimensional systems. J. Phys. C. 6. 1181 (1973).
Mermin, N. D. & Wagner, H. Absence of Ferromagnetism or Antiferromagnetism in One or TwoDimensional Isotropic Heisenberg Models. Phys. Rev. Lett. 17, 1133 (1966).
Davidsen, J. & Kapral, R. Phase synchronization and topological defects in inhomogeneous media. Phys. Rev. E. 66, 055202(R) (2002).
Singh, R., Xu, J., Garnier, N. G., Pumir, A. & Sinha, S. SelfOrganized Transition to Coherent Activity in Disordered Media. Phys. Rev. Lett. 108, 068102 (2012).
Groβmann, G., Peruani, F. & Bär, M. Superdiffusion, largescale synchronization and topological defects. Phys. Rev. E. 93, 040102(R) (2016).
Basnarkov, L. & Urumov, V. Phase transitions in the Kuramoto model. Phys. Rev. E. 76, 057201 (2007).
Hong, H., Park, H. & Choi, M. Y. Collective synchronization in spatially extended systems of coupled oscillators with random frequencies. Phys. Rev. E. 72, 036217 (2005).
Lebrun, R. et al. Mutual synchronization of spin torque nanooscillators through a nonlocal and tunable electrical coupling. arXiv:1601.01247.
Sakaguch, H. & Kuramoto, Y. A Soluble Active Rotater Model Showing Phase Transitions via Mutual Entertainment. Prog. Theor. Phys. 76, 576 (1986).
Vansteenkiste, A. et al. The design and verification of MuMax3. AIP Advances 4, 107133 (2014).
Acknowledgements
This work was supported by the Norwegian Research Council (NFR), project number 216700. V.F. acknowledge partial funding obtained from the Norwegian PhD Network on Nanotechnology for Microsystems, which is sponsored by the Research Council of Norway, Division for Science, under contract no. 221860/F40. F.M. acknowledges financial support from the Ramón y Cajal program through RYC201416515 and from the MINECO through the Severo Ochoa Program for Centers of Excellence in R&D (SEV 20150496).
Author information
Affiliations
Contributions
V.F. initiated the project, developed the model, performed the calculations/simulations and wrote the manuscript. F.M. and E.W. supervised the project and provided valuable input during the analysis/discussion of results and writing of the manuscript.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Electronic supplementary material
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Flovik, V., Macià, F. & Wahlström, E. Describing synchronization and topological excitations in arrays of magnetic spin torque oscillators through the Kuramoto model. Sci Rep 6, 32528 (2016). https://doi.org/10.1038/srep32528
Received:
Accepted:
Published:
Further reading

Emulating the local Kuramoto model with an injectionlocked photonic crystal laser array
Scientific Reports (2021)

Higherorder and longrange synchronization effects for classification and computing in oscillatorbased spiking neural networks
Neural Computing and Applications (2021)

Chimeras and complex cluster states in arrays of spintorque oscillators
Scientific Reports (2017)

A NanotechnologyReady Computing Scheme based on a Weakly Coupled Oscillator Network
Scientific Reports (2017)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.