Advances in Condensed Matter Physics

Advances in Condensed Matter Physics / 2012 / Article
Special Issue

Low-Dimensional Magnetic Systems

View this Special Issue

Research Article | Open Access

Volume 2012 |Article ID 951976 |

M. Carpentieri, G. Finocchio, A. Giordano, B. Azzerboni, F. Lattarulo, "Micromagnetic Study of Synchronization of Nonlinear Spin-Torque Oscillators to Microwave Current and Field", Advances in Condensed Matter Physics, vol. 2012, Article ID 951976, 5 pages, 2012.

Micromagnetic Study of Synchronization of Nonlinear Spin-Torque Oscillators to Microwave Current and Field

Academic Editor: Giancarlo Consolo
Received30 Mar 2012
Revised17 Jul 2012
Accepted23 Jul 2012
Published30 Aug 2012


The nonautonomous dynamics of spin-torque oscillators in presence of both microwave current and field has been numerically studied in nanostructured devices. When both microwave current and field are applied at the same frequency, integer phase locking at different locking ratio is found. In the locking region, a study of the intrinsic phase shift between the locking force (current or field) and the giant magnetoresistive signal as a function of the bias current is also exploited.

1. Introduction

In the last years, the effects due to a direct transfer of the spin angular momentum [1โ€“4] in nanomagnets (magnetization reversal [5, 6] or persistent oscillation of the magnetization [7, 8]) have opened new perspectives in the field of nanotechnology. In particular, one of the most promising applications is the possibility to obtain a competitive generation of nanoscale microwave oscillators, namely, spin-transfer torque oscillators (STOs) [9, 10]. Nowadays, STO is promising from a technological point of view being one of the smallest auto-oscillators observed in nature. It exhibits properties such as frequency tunability on bias current and field and narrow linewidth. To find practical application, STOs have to improve their output power. Some years ago, some researchers demonstrated that synchronized oscillators provided increased output power [11, 12]. This effect, which is receiving a great deal of attention due to its potential applications in telecommunications, has been studied both experimentally and theoretically for the purpose of understanding and realizing mutual synchronization between two or more STOs for microwave source applications. Later, other studies showed the synchronization of serially connected STOs governed by phase locking to a microwave current offering a valid approach to fabricate STOs for output power levels closer to 1โ€‰ฮผW [13, 14].

The synchronization phenomenon is based on the well-known property that when the external frequency ๐‘“RF of a โ€œweakโ€ microwave current or field is close to the free running oscillation frequency ๐‘“0 of the STO, the self-oscillation mode moves and locks to the external frequency. Locking phenomena are also present for the ratio ๐‘Ÿ=๐‘“RF/๐‘“0 close to all integers (๐‘Ÿ=1, 2, 3, etc.) and several rational values [15].

We recently studied the injection locking phenomenon based on the application of a microwave field on perpendicular materials [16]. These systems mainly operate under โ€œweakโ€ microwave signal regime (the power of the microwave signal is negligible compared to the self-oscillation one) [9, 16โ€“18].

In this work, by means of a micromagnetic study, the nonautonomous dynamical behavior of STO in presence of microwave signal composed by the simultaneous application of microwave current density ๐ฝAC and field โ„ŽAC (both at the same frequency) is studied. Also, the influence of a static field on the frequency behavior is investigated.

2. Device and Numerical Details

We studied the dynamical behavior of exchange bias spin-valves composed by IrMn(8โ€‰nm)/Py(10โ€‰nm) (polarizer)/Cu(10โ€‰nm)/Py(4โ€‰nm) (free layer) with elliptical cross-sectional area (120โ€‰nmร—60nm) (see inset of Figure 1(a)). A Cartesian coordinate system has been introduced, where the x- and the y-axes are, respectively, related to the easy and hard in-plane axes of the ellipse. Our numerical experiment is based on the numerical solution of the Landau-Lifshitz-Gilbert-Slonczweski (LLGS) equation [1โ€“3]. In addition to the standard effective field (external, exchange, and self-magnetostatic), the Oersted field and the magnetostatic coupling with the polarizer are taken into account. The resulting equation (LLG+๐‘†) is expressed as ๐‘‘๐ฆ๐‘‘๐œ=โˆ’๐ฆร—๐กe๏ฌ€+๐›ผ๐‘€๐‘ ๎‚€๐ฆร—๐‘‘๐ฆ๎‚๐‘‘๐œโˆ’๐œ’๐ฆร—(๐ฆร—๐ฉ),(1) where ๐ฆ=๐Œ/๐‘€๐‘† and ๐ก๐ž๐Ÿ๐Ÿ=๐‡๐ž๐Ÿ๐Ÿ/๐‘€๐‘† are the dimensionless magnetization vector and the effective field, respectively, ๐œ=๐›พ0๐‘€๐‘ ๐‘ก is the dimensionless time, ๐›พ0 is the gyromagnetic ratio, and ๐‘€๐‘  is the saturation magnetization of the free layer. The first term on the right-hand side of (1) represents the precessional torque around the effective field, whereas the second one is the phenomenological dissipation term, ๐›ผ being the so-called Gilbert parameter. The third term in (1) is the dimensionless Slonczewski spin-transfer torque, where ๐ฉ=๐/๐‘€๐‘  represents the dimensionless magnetization in the pinned layer, and the prefactor ๐œ’ is given by ๐œ‡๐œ’=๐ต๐›พ0๐ฝapp๐‘€2๐‘ ๐‘’๐‘‘๐‘”(๐ฆ,๐ฉ),(2) where e and ๐œ‡๐ต are the electric charge and the Bohr magneton respectively, ๐ฝapp is the current per unit area (density current), and d is the thickness of the free layer, and g(m,p) is the polarization function depending on the relative orientation of the magnetizations [1โ€“3].

For a complete model description of the numerical techniques see also [19โ€“22]. Typical parameters for the Py have been used: saturation magnetization ๐‘€๐‘†=650ร—103โ€‰A/m,exchange constant ๐ด=1.3ร—10โˆ’11 J/m, damping parameter ๐›ผ=0.02, and polarization factor ๐œ‚=0.3 [1โ€“3]. The bias field is applied out-of-plane (z-direction) with a tilted angle of 10ยฐ along the x-axis. The polarizer is considered fixed along the x-direction. To study the locking, we consider a microwave current ๐ฝRF=๐ฝ๐‘€sin(2๐œ‹๐‘“AC๐‘ก+๐œ‹/2) (๐ฝ๐‘€โ‰ค2ร—107โ€‰A/cm2) and a microwave field linearly polarized at ๐œ‹/4 in the x-y plane โ„ŽRF=โ„Ž๐‘€sin(2๐œ‹๐‘“AC๐‘ก+๐œ‹/4)ฬ‚๐‘ฅ+โ„Ž๐‘€sin(2๐œ‹๐‘“AC๐‘ก+๐œ‹/4)ฬ‚๐‘ฆ (โ„Ž๐‘€โ‰ค3โ€‰mT). This microwave field can be generated by using the experimental technique developed in [25]. All the computations have been performed with no thermal effects.

3. Micromagnetic Results and Discussion

In order to characterize the device behavior, first of all we analyzed the STO in the free running regime. We observe dynamical regime in a wide range of current density for bias field larger than 180โ€‰mT. Here we discuss in detail data for a bias field of 250โ€‰mT, but qualitative similar results have been also observed for 200 and 300โ€‰mT.

In order to characterize the oscillator regime of the device, we swapped the dc current exceeding the critical current value (to obtain dynamics regime, threshold current density was ๐ฝ=3ร—107โ€‰A/cm2) up to current values, where oscillation regime is degraded by noise. Frequency and power behavior of the nano-oscillator with respect to the current density is shown in Figures 1(a) and 1(b) (no microwave signal). The frequency curve ๐‘“0 as a function of J presents red shift from the critical current up to ๐ฝ1=โˆ’3.5ร—107โ€‰A/cm2, where the dynamics is characterized by an in-plane oscillation axis. For |๐ฝ|>|๐ฝ1| the magnetization precesses around an out-of-plane axis and the blue shift is achieved. The discontinuities observed in the oscillation frequency are related to jumps of the oscillation axis that correspond to transitions between strongly nonlinear oscillation modes. Similar discontinuities have been also observed in some experimental work [23]. In addition, we performed micromagnetic simulations increasing the out-of-plane static field up to 600โ€‰mT, and we observed that for out-of-plane bias field larger than 400โ€‰mT the red-shift zone disappears, jumps are softer, and the synchronization data achieved are similar to results already published in the literature (see e.g., [10] for a review). At larger field (>600โ€‰mT), the polarizer is moved from the x-direction towards the out-of-plane z-direction that might generate additional noise showing not coherent precessional states.

Figure 1(b) shows the power versus current behavior. Nonlinear power strongly increases at low current and for current values greater than 4ร—107โ€‰A/cm2 is about constant.

We systematically studied the locking to the first harmonic (the same of the self-oscillation) in the blue shift region as a function of the ๐ฝ๐‘€ and โ„Ž๐‘€. Figure 2 shows the precession frequency of the GMR signal as a function of the microwave source frequency in two cases: RF current only and RF current and field together (for ๐ปDC=250โ€‰mT and ๐ฝDC=โˆ’5ร—107โ€‰A/cm2). We found different locking regions at the locking ratio 1โ€‰:โ€‰1, 2โ€‰:โ€‰1, and 3โ€‰:โ€‰1 (in the last case only when the microwave field component is applied, not shown here). Typically, the locking region is much larger when the microwave force is a field (or a combination of current and field). In fact, in the case of current we found a locking region of about 150โ€‰MHz (1โ€‰:โ€‰1) and 50โ€‰MHz (2โ€‰:โ€‰1), whereas no locking on the third harmonic is found. Microwave field provides a locking region larger than 1โ€‰GHz and, since driving force breaks the oscillation symmetry, the 1โ€‰:โ€‰1 synchronization region has a specific asymmetric shape. Then, whereas for small forcing signal the synchronization region can be described by an analytical theory (symmetric tongue where the locking region increases linearly with force amplitude) [24], it cannot be described analytically when increasing the forcing signal and their precise determination requires necessarily specific numerical techniques.

Figure 3(a) summarizes the Arnold tongue (๐ฝ=โˆ’5ร—107โ€‰A/cm2) computed up to ๐ฝ๐‘€=1ร—107โ€‰A/cm2 (โ„Ž๐‘€=0โ€‰mT) and then increasing โ„Ž๐‘€ up to 3โ€‰mT with ๐ฝ๐‘€=1ร—107โ€‰A/cm2 held unchanged. The border lines have been computed considering the lower (in the left part) and higher (in the right part) microwave frequency where the phase locking is achieved. In the low regime of microwave source, the Arnold tongue is related to the only application of the microwave current, which can be considered as a โ€œweakโ€ microwave signal. In fact, such a signal gives rise to symmetric synchronization region (no hysteresis is observed) with a locking band linearly dependent on the force locking.

When both microwave current and field are applied simultaneously at the same frequency, the nonautonomous response becomes more complicated. The presence of an additional weak microwave field gives rise to increasing of the locking region from 150โ€‰MHz at โ„Ž๐‘€=0โ€‰mT to 1.3โ€‰GHz for โ„Ž๐‘€=1โ€‰mT. As can be observed the locking region is strongly asymmetric. This is caused by the strong nonlinearity of the dependence of the auto-oscillation frequency on the oscillation power [25].

Figure 3(b) shows the phase difference between the magnetization oscillation and the microwave source (negative angle means magnetization in delay with respect to microwave source, ฮจ being the phase of the natural precession of the magnetization and ฮจ๐‘’ the force phase) when the microwave component is applied at the same frequency of the free precession one. As shown, typically in both cases (microwave field or current), the phase increases with current and decreases after a maximum value. The phase shift depends on the initial detuning, and it goes to zero only if there is no detuning, that is, in the center of the synchronization region. It should be remembered that the phase difference depends on the initial phase of the force and the way it affects the oscillator. We emphasize here that the phase locking implies that phase difference will be kept bounded inside a finite range of detuning, that is, within the synchronization region. In this case the microwave current and field give the same qualitative behavior to the oscillator, but the phase difference is translated by an angle of about 120ยฐ. As expected, phase difference between microwave source and the magnetization precession follows the frequency behavior. It initially decreases with dc current amplitude (precession frequency decreases) and then increases with dc current (precession frequency increases) up to current values of the order of 7ร—107โ€‰A/cm2. After that current value, phase shift jumps about 180ยฐ, typically this is due to the different oscillation mode from in-plane to out-of-plane mode. Lastly, the phase difference gradually increases following the frequency slow rising with dc current magnitude [26].

In the locking region an intrinsic phase shift ฮจ๐‘–, computed as the difference between the phase of the self-oscillation ฮจ and the phase of the microwave current ฮจ๐‘’, is found.

Figure 4 summarizes ฮจ๐‘– as a function of the microwave frequency, as can be observed a linear relationship between ฮจ๐‘– and ๐‘“AC is achieved with a range of ฮจ๐‘– which can cross 0 or ๐œ‹/2 depending on the bias current density. As reported by Slavin and Tiberkevich [10], analytical formulation of the phase difference in the locking region is given by ฮฆ0๎‚ต๐‘“=arcsinACโˆ’๐‘“0๎‚ถฮ”๐‘“โˆ’arctan(๐œ),(3) where ๐œ=(๐‘/(๐บ+โˆ’๐บโˆ’)) is the nonlinear frequency shift. We found N which is characterized from two different values of nonlinear frequency shift ๐‘=2๐œ‹(๐‘‘๐‘“/๐‘‘๐‘), ๐‘ is the oscillation power, and ๐บ+ and ๐บโˆ’ are the non-linear damping coefficients as in [10]. Phase difference is strongly dependent on the frequency in the locking region [27]. For all the cases, the phase difference decreases from left to right in the Arnold tongue. In particular, the slope of the curve is strong for RF current locking whereas it is more soft for RF field, a result also in substantial agreement with the analytical theory. Considering (3) for RF applied field, the locking bandwidth and nonlinear frequency shift are larger with respect to ac current. A good agreement between analytical and numerical data is found.

4. Conclusions

In summary, we have studied micromagnetically the nonlinear behavior of spin-torque nano-oscillators in locking regime driven by microwave current and field. We found a large locking region at different harmonics when RF field is applied. The effect of the static applied field is also studied. Finally, we showed and explained the intrinsic phase shift due to the difference between microwave source and magnetization precession inside the locking region, also comparing our numerical data with a recent analytical theory.


This paper was supported by Spanish Project under Contract no. MAT2011-28532-C03-01. The authors would like to thank Sergio Greco for his support with this paper.


  1. J. C. Slonczewski, โ€œCurrent-driven excitation of magnetic multilayers,โ€ Journal of Magnetism and Magnetic Materials, vol. 159, no. 1-2, pp. L1โ€“L7, 1996. View at: Publisher Site | Google Scholar
  2. J. C. Slonczewski, โ€œExcitation of spin waves by an electric current,โ€ Journal of Magnetism and Magnetic Materials, vol. 195, no. 2, pp. L261โ€“L268, 1999. View at: Publisher Site | Google Scholar
  3. J. C. Slonczewski, โ€œCurrents and torques in metallic magnetic multilayers,โ€ Journal of Magnetism and Magnetic Materials, vol. 247, no. 3, pp. 324โ€“338, 2002. View at: Publisher Site | Google Scholar
  4. L. Berger, โ€œEmission of spin waves by a magnetic multilayer traversed by a current,โ€ Physical Review B, vol. 54, no. 13, pp. 9353โ€“9358, 1996. View at: Publisher Site | Google Scholar
  5. I. N. Krivorotov, N. C. Emley, J. C. Sankey, S. I. Kiselev, D. C. Ralph, and R. A. Buhrman, โ€œTime-domain measurements of nanomagnet dynamics driven by spin-transfer torques,โ€ Science, vol. 307, no. 5707, pp. 228โ€“231, 2005. View at: Publisher Site | Google Scholar
  6. D. C. Ralph and M. D. Stiles, โ€œSpin transfer torques,โ€ Journal of Magnetism and Magnetic Materials, vol. 320, no. 7, pp. 1190โ€“1216, 2008. View at: Publisher Site | Google Scholar
  7. S. I. Klselev, J. C. Sankey, I. N. Krivorotov et al., โ€œMicrowave oscillations of a nanomagnet driven by a spin-polarized current,โ€ Nature, vol. 425, no. 6956, pp. 380โ€“383, 2003. View at: Publisher Site | Google Scholar
  8. A. A. Tulapurkar, Y. Suzuki, A. Fukushima et al., โ€œSpin-torque diode effect in magnetic tunnel junctions,โ€ Nature, vol. 438, no. 7066, pp. 339โ€“342, 2005. View at: Publisher Site | Google Scholar
  9. W. H. Rippard, M. R. Pufall, S. Kaka, S. E. Russek, and T. J. Silva, โ€œDirect-current induced dynamics in Co90Fe10/Ni80Fe20 point contacts,โ€ Physical Review Letters, vol. 92, no. 2, Article ID 027201, 4 pages, 2004. View at: Google Scholar
  10. A. Slavin and V. Tiberkevich, โ€œNonlinear auto-oscillator theory of microwave generation by spin-polarized current,โ€ IEEE Transactions on Magnetics, vol. 45, no. 4, pp. 1875โ€“1918, 2009. View at: Publisher Site | Google Scholar
  11. S. Kaka, M. R. Pufall, W. H. Rippard, T. J. Silva, S. E. Russek, and J. A. Katine, โ€œMutual phase-locking of microwave spin torque nano-oscillators,โ€ Nature, vol. 437, no. 7057, pp. 389โ€“392, 2005. View at: Publisher Site | Google Scholar
  12. F. B. Mancoff, N. D. Rizzo, B. N. Engel, and S. Tehrani, โ€œPhase-locking in double-point-contact spin-transfer devices,โ€ Nature, vol. 437, no. 7057, pp. 393โ€“395, 2005. View at: Publisher Site | Google Scholar
  13. J. Grollier, V. Cros, and A. Fert, โ€œSynchronization of spin-transfer oscillators driven by stimulated microwave currents,โ€ Physical Review B, vol. 73, no. 6, Article ID 060409(R), 4 pages, 2006. View at: Publisher Site | Google Scholar
  14. W. H. Rippard, M. R. Pufall, S. Kaka, T. J. Silva, S. E. Russek, and J. A. Katine, โ€œInjection locking and phase control of spin transfer nano-oscillators,โ€ Physical Review Letters, vol. 95, no. 6, Article ID 067203, pp. 1โ€“4, 2005. View at: Publisher Site | Google Scholar
  15. S. Urazhdin, P. Tabor, V. Tiberkevich, and A. Slavin, โ€œFractional synchronization of spin-torque nano-oscillators,โ€ Physical Review Letters, vol. 105, no. 10, Article ID 104101, 4 pages, 2010. View at: Publisher Site | Google Scholar
  16. M. Carpentieri, G. Finocchio, B. Azzerboni, and L. Torres, โ€œSpin-transfer-torque resonant switching and injection locking in the presence of a weak external microwave field for spin valves with perpendicular materials,โ€ Physical Review B, vol. 82, no. 9, Article ID 094434, 8 pages, 2010. View at: Publisher Site | Google Scholar
  17. G. Consolo, V. Puliafito, G. Finocchio et al., โ€œCombined frequency-amplitude nonlinear modulation: theory and applications,โ€ IEEE Transactions on Magnetics, vol. 46, no. 9, pp. 3629โ€“3634, 2010. View at: Publisher Site | Google Scholar
  18. R. Bonin, G. Bertotti, C. Serpico, I. D. Mayergoyz, and M. D'Aquino, โ€œAnalytical treatment of synchronization of spin-torque oscillators by microwave magnetic fields,โ€ European Physical Journal B, vol. 68, no. 2, pp. 221โ€“231, 2009. View at: Publisher Site | Google Scholar
  19. A. Romeo, G. Finocchio, M. Carpentieri, L. Torres, G. Consolo, and B. Azzerboni, โ€œA numerical solution of the magnetizationreversal modeling in a permalloy thin film using fifth order Runge-Kutta method with adaptive step size control,โ€ Physica B, vol. 403, no. 2-3, pp. 464โ€“468, 2008. View at: Publisher Site | Google Scholar
  20. E. Martinez, L. Torres, L. Lopez-Diaz, M. Carpentieri, and G. Finocchio, โ€œSpin-polarized current-driven switching in permalloy nanostructures,โ€ Journal of Applied Physics, vol. 97, no. 10, Article ID 10E302, 3 pages, 2005. View at: Publisher Site | Google Scholar
  21. Martinez et al., โ€œstandard problem #4 report,โ€ The free layer has been discretized in computational cells of 5×5×4 nm3. The time step used was 32 fs, View at: Google Scholar
  22. M. Carpentieri, L. Torres, B. Azzerboni, G. Finocchio, G. Consolo, and L. Lopez-Diaz, โ€œMagnetization dynamics driven by spin-polarized current in nanomagnets,โ€ Journal of Magnetism and Magnetic Materials, vol. 316, no. 2, pp. 488โ€“491, 2007. View at: Publisher Site | Google Scholar
  23. I. N. Krivorotov, D. V. Berkov, N. L. Gorn et al., โ€œLarge-amplitude coherent spin waves excited by spin-polarized current in nanoscale spin valves,โ€ Physical Review B, vol. 76, no. 2, Article ID 024418, 14 pages, 2007. View at: Publisher Site | Google Scholar
  24. P. Maffezzoni, โ€œComputing the synchronization regions of injection-locked strongly nonlinear oscillators for frequency division applications,โ€ IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, vol. 29, no. 12, pp. 1849โ€“1857, 2010. View at: Publisher Site | Google Scholar
  25. P. Tabor, V. Tiberkevich, A. Slavin, and S. Urazhdin, โ€œHysteretic synchronization of nonlinear spin-torque oscillators,โ€ Physical Review B, vol. 82, no. 2, Article ID 020407, 4 pages, 2010. View at: Publisher Site | Google Scholar
  26. Y. Zhou, J. Persson, S. Bonetti, and J. Akerman, โ€œTunable intrinsic phase of a spin torque oscillator,โ€ Applied Physics Letters, vol. 92, no. 9, Article ID 092505, 3 pages, 2008. View at: Publisher Site | Google Scholar
  27. G. Finocchio, G. Siracusano, V. Tiberkevich, I. N. Krivorotov, L. Torres, and B. Azzerboni, โ€œTime-domain study of frequency-power correlation in spin-torque oscillators,โ€ Physical Review B, vol. 81, no. 18, Article ID 184411, 6 pages, 2010. View at: Publisher Site | Google Scholar

Copyright © 2012 M. Carpentieri et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.