Efficient Synchronization of Dipolarly Coupled Vortex-Based Spin Transfer Nano-Oscillators.

Due to their nonlinear properties, spin transfer nano-oscillators can easily adapt their frequency to external stimuli. This makes them interesting model systems to study the effects of synchronization and brings some opportunities to improve their microwave characteristics in view of their applications in information and communication technologies and/or to design innovative computing architectures. So far, mutual synchronization of spin transfer nano-oscillators through propagating spinwaves and exchange coupling in a common magnetic layer has been demonstrated. Here we show that the dipolar interaction is also an efficient mechanism to synchronize neighbouring oscillators. We experimentally study a pair of vortex-based spin transfer nano-oscillators, in which mutual synchronization can be achieved despite a significant frequency mismatch between oscillators. Importantly, the coupling efficiency is controlled by the magnetic configuration of the vortices, as confirmed by an analytical model and micromagnetic simulations highlighting the physics at play in the synchronization process.

Beyond the traditional applications to data storage and field sensors, the recent progresses in spin transfer physics allows a widening of the application spectra for spintronics devices, notably toward multifunctional devices12 relying on their nonvolatile nature, scalability, and compatibility with existing CMOS processes. Among these novel opportunities, it is anticipated that spin-transfer nano-oscillators (STNOs) can represent a breakthrough in the next generation of information and communication technologies. These devices execute self-sustained magnetic oscillations in the GHz range that are induced by the spin-transfer torque and which can be efficiently converted into resistance and voltage oscillations through magneto-resistive effects3. To achieve these expectations, STNOs have been intensively studied in the last decade in order to improve the understanding of the mechanisms of the spin transfer torque (STT) as well as the physics of high frequency nonlinear magnetization dynamics4.

One of the most important properties of STNOs is their nonlinear nature and the associated (natural) tendency to adapt their own frequency to any change in their environment. On one hand, this effect contributes in certain limits to a broadening of the spectral linewidth, deteriorating the quality factor of the STNOs5. On the other hand, STNOs appear as model nanoscale systems for studying the effects of synchronization, notably in the regime of large nonlinearities6 and even their tendency to reach chaotic regimes under the influence of spin transfer forces7. Moreover, achieving synchronization of STNOs aims at improving the self-sustained oscillations stability that is crucial for radio-frequency applications as nanoscale tunable radiofrequency source or radiofrequency detector8910, and enables the development of innovative computing architectures, particularly oscillator-based associative memories based on information coding in the individual phases of large scale arrays of interacting STNOs1112.

The forced or non-autonomous dynamics of STNOs has been mainly investigated through the study of their ability to achieve phase locking to an external excitation being an external rf field4131415 or rf current1617181920212223. A further advance has been the observation of mutual synchronization between several STNOs. Different coupling mechanisms have been investigated, such as propagating spin-waves in an extended magnetic media242526 or magnetic coupling through exchange interactions27. However, these approaches to couple STNOs might be arduous for reaching synchronized states in large arrays of STNOs. Two other strategies have been recently investigated theoretically in order to achieve mutual synchronization of a large number of STNOs, namely the electric coupling by self-modulation of the current flowing through each device28 and the magnetic coupling between the oscillators arising from the dipolar influence of closely packed STNOs29303132. Our main objective is to propose a complete study combining experimental results and numerical as well as theoretical investigation of the mutual synchronization through dipolar coupling of two adjacent vortex-based STNOs. We demonstrate the ability of the two STNOs to achieve mutual synchronization despite a significant frequency mismatch. Moreover, we develop an analytical model describing the dynamics of the synchronization process and highlight the important role of the magnetic parameters (that are the core polarity and vortex chirality of the vortices) of the two STNOs for the optimization of the coupling efficiency.

Synchronization of two neighbouring vortex-based STNOs

Experimental results

Our system is composed of two 2R = 200 nm diameter spin valve nanopillars separated by L = 100 nm. These STNOs are patterned by e-beam lithography and ion etching from a magnetic trilayer deposited by sputtering: NiFe(15 nm)/Cu(8 nm)/NiFe(4 nm) (see Fig. 1). From our previous studies on a single STNO made from the same trilayer stack, we observed that the magnetic configuration in the thicker NiFe layer is a magnetic vortex at remanence. Depending on the applied dc current, the configuration in the thin NiFe layer can be either a quasi-uniform magnetization or a second magnetic vortex333435. In both cases, the resulting spin transfer forces on the thick layer vortex can lead to the sustained excitation of the vortex core gyration if the current sign is appropriately chosen. In the case of two vortices, an important outcome was to demonstrate that highly coherent gyrotropic oscillations of the coupled vortices can be achieved when their core polarities are in opposite directions33343536. Further description of the sample can be found in the methods section.

Figure 1

Scheme of our coupled STNOs system.

Two spin valve nanopillars with 2R = 200 nm diameter separated by a distance L = 100 nm. Each pillar contains a NiFe(4 nm)/Cu(10 nm)/NiFe(15 nm) spin valve (NiFe = Ni81Fe19) and a vortex is present in each 15 nm NiFe layers. The current is injected in parallel into the two pillars.

Spintronic oscillators based on the dynamics of a magnetic vortex core have recently proven interest by achieving record quality factors and signal amplitudes373839. The main purpose here will be to further consider gyrotropic spin transfer vortex dynamics as a model system to investigate the synchronization of dipolarly coupled oscillators.

The dynamic behaviour of a magnetic vortex core is known to be well described through a single collective variable equation, the so-called Thiele equation describing the vortex core motion through its vector coordinates (Ref. 40, 41, 42, 43, 44, 45):

The gyrotropic force (left hand side of Eq. 1) is equal to the sum of the effective dissipative and conservative forces acting on the vortex core, which can be derived from the integration of the local magnetization dynamics. The three forces acting on the vortex core expressed in the right part of Eq. 1 are: the viscous damping force (proportional and opposed to the core velocity), the spring-like confinement force (opposed to the core displacement), and an additional force accounting for the effect of spin transfer torques. Details of the introduced parameters are given in supplementary informations. To achieve the dynamical regime of self-sustained gyration of the vortex core around the dot centre, the energy provided by the current and the associated spin transfer effect must fully compensate the energy dissipation. Two important conclusions arising from Thiele equation are that (i) the sense of gyration of the vortex core is directly related to the direction of the core polarity43 and (ii) the gyrotropic mode frequency can be tuned by the application of a perpendicular magnetic field Hperp with a quasi-linear evolution: where ω0 is the zero-field frequency and HS is the saturation field of the ferromagnet with a slope sign that is directly related to the sign of the vortex polarity P = ±1 (Ref. 46). In the following, we will use these two properties in order to reliably prepare the system of two neighbouring vortex-based STNOs in two different configurations being either parallel core polarities (Pc) or opposite core polarities (APc) between the two vortices in each STNO’s thick layer.

As shown in Fig. 1, the two STNOs are electrically connected in parallel. When the gyrotropic frequencies of the two STNOs are far from each other, we measure two independent peaks on the spectrum analyser, representative of the vortex dynamics of the two STNOs, as illustrated in Fig. 2(c). The oscillating in-plane stray field associated with the gyrotropic motion of each vortex will be the source of the dipolar dynamical field used to couple the two STNOs. Contributions of the out-of-plane components of the magnetization to the dynamical coupling are negligible because of the small volume of the vortex cores.

Figure 2

(a) Power spectrum map versus perpendicular field measured on Sample #1, in a case when vortices have opposite core polarities. The injected current through the two pillars is I = 50 mA. Insert shows evolution of frequencies versus field for the non-synchronized oscillators and their predicted linear evolution. Power spectra obtained at (b) 0.5 kOe and (c) 0.7 kOe, before (b) and after (c) synchronization is achieved.

In order to investigate the regime of potential synchronization through dynamical dipolar coupling, we have prepared the system of STNOs in order to have the oscillating vortices with opposite core polarities (APc) which is reflected by the two opposite slopes in their dependence on perpendicular field shown in Fig. 2(a), measured under a bias current I = +50 mA. As shown in Fig. 2, when starting from Hperp = 1.2 kOe, the frequency of the two STNOs are about 200 MHz apart from each other. This energy difference is very large compared to the coupling energy and the two STNOs almost do not feel each other. Then a mean to tune the frequency mismatch between the two STNOs is to sweep down the perpendicular field. As shown in Fig. 2(c) for Hperp = 0.7 kOe, the two STNOs are still behaving almost independently as two separate peaks can be measured. A drastic change is observed for a field slightly smaller field than Hperp = 0.6 kOe for which a single peak is recorded (see Fig. 2(b)). The observation of such a transition from two peaks associated to each auto-oscillator to a single peak is the evidence of the synchronization of the two STNOs. To discard the possibility that this transition might be related to a magnetic switching (e.g., reversal of a vortex polarity), we have verified that this phenomenon is reversible by sweeping the field back and forth around the critical field value. Indeed, the reversal of the core polarity of one of the two vortices would lead to the observation of an hysteretic transition.3547. Frequency pulling can also be observed, which is characteristic of a transition towards synchronized oscillations48. The measured critical frequency mismatch before synchronization is .

Note that the diameters of two STNOs are not perfectly identical (about 3% difference) which results in two different frequencies at zero applied field as shown in the inset of Fig. 2 by extrapolating the unsynchronized frequency dependence. This difference in diameters explains why the synchronization is not observed around zero field, since the mode crossing occurs at positive fields . It is also noticed that the power amplitudes of the signals associated to each STNO when they are not synchronized are very different. Yet, this difference does not necessarily reflect a strong difference in the vortex gyration amplitudes as the associated voltage oscillations depend strongly on the actual magnetic distribution of the polarizing layers49.

Moreover, no notable modification of the linewidth was observed, being about 500 kHz before and after the synchronization. This differs from some of the previous observations, where STNOs synchronization was achieved through spin-wave coupling and a clear reduction of the peak linewidth was highlighted2425.

Our method to tune the frequency difference using a perpendicular field is obviously less efficient when the two core polarities are identical (Pc). In order to be able to study the synchronization in case of parallel core polarity, we have measured another pair of STNOs which have more separate natural gyrotropic frequencies. As shown in Fig. 3(b), a reversible transition from independent to synchronized auto-oscillations is observed. Note that the maximum frequency mismatch before synchronization is Δf = 20 ± 2 MHz for these measurements with identical core polarities. In order to make a quantitative comparison, we have also measured on the same system of coupled STNOs the synchronization in case of opposite core polarities (APc). As shown in Fig. 3(a), for the same bias current I = +35 mA, we find that the maximum frequency mismatch is Δf = 48 ± 2 MHz. On this measurement, the appearance of a secondary peak around 800 MHz before synchronization is again characteristic of the interaction of the two oscillators, corresponding to the frequency beating phenomenon48.

Figure 3

Power spectrum maps versus perpendicular field measured on Sample #2, in a case when vortices have opposite polarities (a) and when vortices have identical polarities (b). The injected current is I = 35 mA.

A striking result is thus that Δf > Δf, suggesting that the coupling energy between STNOs is more efficient to induce synchronization when the two interacting vortices have opposite core polarities. It is important to emphasize that the two cases mostly differ through the relative direction of gyration of the two vortices: the two vortex cores have the same sense of gyration when P1P2 > 0, whereas they gyrate in opposite directions when P1P2 < 0.

It is interesting to note that after synchronization, the resulting frequency and associated dependence to Hperp of the synchronized oscillations are close to one of the two auto-oscillators. Indeed, in each presented measurement, we can identify a leader oscillator, and a follower. We interpret that the latter has a stronger ability to adapt its frequency to external stimuli than the leader, and then experiences a bigger frequency shift to achieve synchronization.

Dipolar interaction between the oscillators

To establish a model, we consider a system composed of the two thick layers’ vortices, interacting through the mutual dipolar fields created by their rotating in-plane net magnetizations29303250515253. Contributions of the out-of-plane components of the magnetization to the coupling are neglected because of the small volume of the vortex cores. Given the large positive current flowing through each pillar, the two vortices have identical chiralities following the direction of the current-induced Oersted field. Moreover, given the smaller volume of the thin 4 nm layers, the possible influence of shallow dynamics occurring in these layers will not be considered and these polarizing layers’ magnetizations will be considered fixed.

We consider that the net in-plane magnetization of the shifted vortex increases linearly with the core displacement from the dot centre: , where is determined analytically54,Ci  = ± 1 is the vortex chirality, and is oriented along the pillar axis. While the vortex is gyrating, its net in-plane magnetization describes a 360° rotation. To model the dipolar coupling between the two vortices, we consider a simple model of two rotating in-plane macro-dipoles and at the disks’ centre positions293032. The interaction energy can then be expressed as:

where and are the polar coordinates of the 1 and 2 vortices cores respectively, and is a unit vector in the direction joining the two dots centres (Fig. 1). The coupling coefficients can then be rewritten as and so that μ(+) = 3μ(-), with D12 = 2R + L the interpillar centre-to-centre distance. In the following, the influence of the perpendicular field Hperp on these two coupling parameters will be neglected, considering that the applied field is a small portion of the saturating field , so that the vortices are not significantly tilted out-of-plane.

The vortex core direction of gyration is directly related to its polarity, which is translated in polar coordinates by: where ω is the absolute gyrotropic frequency of the vortex. As a consequence, the term relating to μ(+) should oscillate at the frequency (P1ω1 + P2ω2), while the second term relating to μ(−) should oscillate at the frequency (P1ω1 − P2ω2). The variation rates of the two terms depend then strongly on the relative polarities of the gyrating vortices.

Analytical model for synchronization

The purpose is now to define the synchronization criterion. We propose here to model the coupled system through two coupled Thiele equations:

where the interaction force is then expressed as: . In the general case of a non-uniform polarizer, with magnetization distributed in plane or out-of-plane, the spin transfer force can be described by the general expression: where J is the current density flowing through the pillar and λ(X, J) expresses the spin transfer efficiency414244455556. By using complex number coordinates, these equations can then be written as:

where are the complex coordinates of the core positions and are their complex conjugates.

The coupling action is considered as a perturbation to the equilibrium auto-oscillations of the isolated pillars. In order to extract a differential equation for the relative phase dynamics, we expand these equations around their equilibrium values through the oscillations power: , and the subsequent deviations of the instantaneous frequencies: , and dissipative forces: where N’s are the coefficients of the non-linear frequency shifts, and are the damping rates for small power deviations6.

Details of the developments of these equations are proposed in supplementary informations. To capture the dominant mechanisms responsible for synchronization, the following approximations are used: , , , , , given the small applied field and the small differences between the two pillars. We further suppose that the two auto-oscillators only differ through their frequencies ω1≠ω2. We then extract the differential equation governing the dynamics of the signed phase difference , in the cases of identical polarities (P1 = P2 = +1) and opposite polarities (P1 = −P2):

where is the normalized dimensionless nonlinear frequency shift. Equations 6 and 7 are similar to the typical Adler equation48. Synchronization can be achieved only if a stable solution exists for this equation. Given that C1 = C2 in our experiment, the conditions on the frequency mismatch Δω for synchronization are:

Equation 8 states that there are two mechanisms contributing to synchronization between the two STNOs. The first mechanism, associated with the αη multiplying coefficient, corresponds to the direct action of the dipolar force on the oscillator phase. The second mechanism, associated with the ν coefficint, corresponds to the action of the dipolar force leading to the modification of the vortex orbits, and subsequently their gyration frequency, to achieve synchronization. For the STNOs under study, we evaluate that and v ≳1 (Ref. 39), highlighting that the second mechanism is clearly the dominant one.

From this analytical study, it appears that the coupling action efficiency depends on the relative polarities. Notably, from the expression of dipolar coupling presented in Eq. 3, this model predicts that when polarities are anti-parallel (APc) the synchronization can occur for a frequency mismatch that is three times larger than when polarities are parallel (Pc): .

These results are in good agreement with the experimental observations we presented. Indeed, the critical frequency mismatches before synchronization measured on the second sample verify . This confirms that the effective coupling energy is stronger when the vortices are gyrating in opposite directions, and that this configuration is preferable to achieve efficient synchronization.

Finally, another prediction from the model is that, depending on the relative polarities, the equilibrium phase shift after synchronization is different, being either close to for identical core polarities or to for opposite core polarities. Changing the relative chiralities signs would reverse this statement, but without affecting the coupling amplitude and the associated synchronization features.

Note that in Eq. 8, we neglected interaction terms associated to fast variations of the coupling energy which do not participate in the synchronization process. However, it should be noted that the high frequency oscillations of the coupling energy still exist and are certainly perturbing the synchronized oscillations. Notably, a major consequence will be to increase the phase fluctuations of the synchronized oscillators, thus increasing the synchronized signal linewidth. This phenomenon might prevent the synchronization phenomenon to increase the coherence of the oscillations.

Study of the dipolar energy versus relative polarities

To gain further insight on the relative efficiency of the dipolar action for synchronization depending on the relative polarities, we present in Fig. 4 a sketch of the temporal evolution of the coupling energy in the ideal condition of phase locking ω1 = ω2 for the two cases of identical polarities (φ2 = φ1) and opposite polarities (φ2 = π − φ1). In the case of opposite polarities, when vortices gyrate in opposite directions, the relative alignment of the net magnetizations is always favourable during the motion, so that the mean interaction energy is large and negative. In contrast, when vortices gyrate in identical direction (polarities are identical), the synchronized gyration brings the system in situations when net magnetizations are in a favourable alignment (Wint < 0) but also in a unfavourable alignment (Wint > 0).

Figure 4

(a) Illustration of the synchronized motion of the two auto-oscillating vortices cores, and of their net in-plane magnetizations (black arrows), for the cases of identical polarities (left) and opposite polarities (right) during one period. The sign of the associated dipolar coupling energy Wint is also given. (b) Predicted evolution of the dipolar coupling energy Wint versus part of oscillation period for identical (blue) and opposite polarities (red). The dotted lines show the average interaction energy in both cases.

To summarize, in case of identical polarities, Wint oscillates with large amplitude with a mean energy , while in case of opposite polarities, Wint oscillates with small amplitude but with , three times larger. In other words, the dipolar interaction is more efficient for stabilizing synchronization when in AP configuration.

Finally, the validity of the macro-dipoles and Thiele approach was evaluated by conducting a micromagnetic study. Using full micromagnetic simulations, we studied the synchronization dynamics between the two identical oscillators, with no frequency mismatch, at zero temperature and zero field. This specific case allowed us to monitor the synchronization phenomenon whatever the amplitude of the coupling. Following the approach described in (Ref. 29,30,32) we extracted the effective coupling coefficient μeff, defined as , as a function of the centre-to-centre distance. The result is presented in Fig. 5, where the ratio μeff/G is plotted for the two polarities configurations. This study allows us to get a quantitative prediction of the critical frequency mismatch for synchronization as a function of the pillars interdistance.

Figure 5

Average interaction energy versus interpillar distance extracted from micromagnetic simulations for the case of two identical synchronized oscillators with radii R = 100 nm, when vortices have identical polarities (blue line and filled squares) and opposite polarities (red line and open squares).

It appears from this study that the prediction is a fair approximation in the present experimental case where D12 = 300 nm. Yet, this ratio (see insert in Fig. 5) is predicted to decrease as the pillars interdistance decreases. Comparison to experimental results suggests that the real interpillar distance is actually smaller than the nominal one.

Discussion

In this study, both experiments and analytical modelling demonstrated that the effective dipolar coupling energy, and the subsequent critical frequency mismatch for synchronization, is stronger for STNOs with interacting vortices gyrating in opposite direction i.e., with opposite polarities, than for interacting vortices gyrating in the same direction i.e., with identical polarities. Synchronization of STNOs with frequency mismatches corresponding, respectively, to approximately 10% and 6% of the mean frequency of the two oscillators have been observed in this optimal configuration. These numbers fall reasonably in the range of observed frequency spreading in networks of STNOs, and the results can be extrapolated for potential synchronization of a large number of oscillators. Applications of such networks have been proposed recently for the realization of neuro-inspired STNO-based associative memories1112. In such architectures, the control over the coupling between neighbouring oscillators is a major asset as it will be the support of the information. Considering the implementation of the coupling between oscillators, these results can be safely extended to other spin-transfer oscillators architectures, such as the ones based on magnetic tunnel junctions, which recently demonstrated higher power with low linewidth signal3857.

Methods

The samples are composed of two 2R = 200 nm circular spin-valve nanopillars separated by a L = 100 nm edge-to-edge distance, patterned in a Cu(60 nm)/NiFe(15 nm)/Cu(10 nm)/NiFe(4 nm)/Au(25 nm) stack (NiFe = Ni81Fe19) by standard e-beam and ion etching lithography process.

Samples with single pillars prepared from the same stack have already been intensively studied in previous works1533343947. Applying a positive current I > 10 mA ensures that the two vortices have identical positive chilarities (direction of the curling magnetization) following the direction of the current induced Oersted field. Control over the vortices polarities can be done by applying a perpendicular magnetic field47.

We note that the induced Oersted field from each pillar is also felt by the neighbouring pillar. This crosstalk field was evaluated between 100 Oe and 200 Oe for a current I = 20 mA flowing in each pillar. The consequence will be an offset of the vortex cores equilibrium positions by a few nanometers, being much smaller than the pillar radius. However, this field is not sufficient to affect the vortex stability.

The two self-sustained oscillators are supplied in parallel through single top and bottom electrodes for current injection (see Fig. 1), and the amplitudes of the current flowing through each pillar are hence not independently controlled. Despite possible small size and shape deviations between the two pillars, we assume that equal current flows through each pillar. In our convention, a positive current corresponds to electrons flowing from the thick to the thin layer. Given the parallel supply, it is not possible to independently measure the magnetoresistive signals associated to each pillar, but only to measure the voltage across the two pillars.

The ac-voltage is isolated from biasing circuit by a bias-T and then measured on a spectrum analyzer after a 30 dB amplification.

Additional Information

How to cite this article: Locatelli, N. et al. Efficient Synchronization of Dipolarly Coupled Vortex-Based Spin Transfer Nano-Oscillators. Sci. Rep. 5, 17039; doi: 10.1038/srep17039 (2015).