# Injection locking at 2f of spin torque oscillators under influence of thermal noise

M. Tortarolo, B. Lacoste, J. Hem, C. Dieudonné, M.-C. Cyrille, A. Katine, D.

Mauri, A. Zeltser, Liliana D. Buda-Prejbeanu, U. Ebels

► To cite this version:

M. Tortarolo, B. Lacoste, J. Hem, C. Dieudonné, M.-C. Cyrille, et al. Injection locking at 2f of spin torque oscillators under influence of thermal noise. Scientific Reports, 2018, 8, pp.1728. 10.1038/s41598-017-18969-5. cea-01696935

HAL Id: cea-01696935

https://cea.hal.science/cea-01696935

Submitted on 24 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Injection locking at 2f of spin torque oscillators under influence of thermal noise

Received: 26 September 2017

Accepted: 18 December 2017

Published online: 29 January 2018

M. Tortarolo<fnref n="1" />, B. Lacoste²,⁵, J. Hem², C. Dieudonné², M.-C. Cyrille³, J. A. Katine⁴, D. Mauri⁴, A. Zeltser⁴, L. D. Buda-Prejbeanu² & U. Ebels²

Integration of Spin Torque Nano-Oscillators STNO's in conventional microwave circuits means that the devices have to meet certain specifications. One of the most important criteria is the phase noise, being the key parameter to evaluate the performance and define possible applications. Phase locking several oscillators together has been suggested as a possible means to decrease phase noise and consequently, the linewidth. In this work we present experiments, numerical simulations and an analytic model to describe the effects of thermal noise in the injection locking of a tunnel junction based STNO. The analytics show the relation of the intrinsic parameters of the STNO with the phase noise level, opening the path to tailor the spectral characteristics by the magnetic configuration. Experiments and simulations demonstrate that in the in-plane magnetized structure, while the frequency is locked, much higher reference currents are needed to reduce the noise by phase locking. Moreover, our analysis shows that it is possible to control the phase noise by the reference microwave current ($I_{RF}$) and that it can be further reduced by increasing the bias current ($I_{DC}$) of the oscillator, keeping the reference current in feasible limits for applications.

The increasing demand on miniaturization and the necessity of implementing more and more frequency standards on a single device (GPS, mobile phone, wifi) require new concepts to cover large frequency ranges at the nanoscale size, keeping costs low. In this sense the Spin STNOs where a spin polarized current passing through the magnetic multi-layered nanosystem can drive the magnetization into large amplitude periodic oscillations<fnref n="1" />⁻³ when the spin polarized current is large enough to compensate the natural damping, are an alternative path to current controlled microwave devices. Despite all their appealing features, as nano size, fast frequency tuning, high modulation speed and CMOS compatibility, one of the main issues that remains to be addressed is their low phase stability leading to a large phase noise⁴⁻⁹, which is the dominant contribution to linewidth in oscillators. Several options were proposed to stabilize the phase: electric coupling¹⁰,¹¹, dipolar coupling¹²⁻¹⁵ or spin wave coupled nanocontacts¹⁶,¹⁷. In order to understand the role of phase noise in electric synchronization of several oscillators by their own emitted microwave current, we studied the injection locking of an STNO to a reference microwave current $I_{RF}$ with known spectral specifications. Here we focus on standard uniform in plane magnetized oscillators (in-plane polarizer and in-plane free layer, IP), for which an in-plane precession (IPP) mode is stabilized. The injection locking of such an STNO to a reference current at two times the generated frequency ($f_{ext} = 2f_o$) was demonstrated both numerically and by experiments¹⁸. However, the linewidth in the locked regime was reduced only by a factor of seven, while a reduction to the linewidth of the microwave source (several Hz) was expected. These large linewidths are associated to the thermal noise that induces fluctuations which can drive the phase from an equilibrium state to a neighbouring one, with an associated phase slip of $\pm 2\pi$ which can be envisaged as non-synchronization and re-synchronization events. Zhou et al.¹⁹ demonstrated that the particular way the phase approaches its synchronized value has consequences in the transients that may limit the modulation of an STNO. Recent works investigated the mechanisms of the so called pure phase locking state in double vortex based STNO: Robust synchronization was experimentally shown, with a $10^5$ linewidth reduction²⁰ and the role of

<footnote n="1">¹Centro Atómico Constituyentes CNEA, 1650 San Martín and Consejo Nacional de Investigaciones Científicas y Técnicas, C1033AAJ, Buenos Aires, Argentina. ²Univ. Grenoble Alpes, CEA, CNRS, Grenoble INP, INAC, SPINTEC, F38000, Grenoble, France. ³Univ. Grenoble Alpes, CEA-LETI MINATEC-CAMPUS, 38000, Grenoble, France. ⁴HGST, 3403 Yerba Buena Road, San Jose, California, 95135, USA. ⁵Present address: International Iberian Nanotechnology Laboratory, Braga, Portugal. Correspondence and requests for materials should be addressed to M.T. (email: tortarol@tandar.cnea.gov.ar)</footnote>
the phase slips in the synchronized state was investigated$^{21}$. In this work we study the injection locking at 2f to an external reference current of a uniform IP magnetized STNO under the influence of thermal noise. Contrary to most studies in injection locking, here the phase noise is investigated in detail to understand the non ideal locking reported on uniform IP magnetized STNO's$^{18}$, providing an analytical expression that relates the phase noise level with the intrinsic parameters of the STNO.

## Analytic model

The effect of thermal fluctuations on the transient behaviour of the injection locking state of an STNO is analyzed in the frame of a generic model of a non-linear auto oscillator$^{21}$ that we extended for the IPP mode synchronized by a reference microwave current at 2f (details in Supplementary material). Since STNO's are non-linear (non-isochronous) oscillators, the power and the phase of the oscillator are not independent, leading to a system of coupled equations.

$$ \frac{d\psi}{dt} = -\Delta\omega + 2N\delta p + \frac{1}{\sqrt{p_o}}\mathrm{Im}[\xi] \quad (1) $$

$$ \frac{d\delta p}{dt} \approx -2\Gamma_p \delta p + 2p_o F \cos(\psi) + 2\sqrt{p_o} \mathrm{Re}[\xi] \quad (2) $$

Here $\psi(t) = 2\Phi - \omega_{ext,t}$ is the phase difference between the STNO phase $\Phi$ and the phase of external source $\omega_{ext,t}$, $N$ is the coefficient of non-linear frequency shift, $F$ is a real parameter proportional to the reference current, $\Gamma_p$ is the damping rate of the power fluctuations and $\xi$ has the statistical properties of the Gaussian thermal noise ($\langle\xi(t)\rangle = 0$, $\langle\xi(t)\xi(t)\rangle = 0$, $\langle\xi(t)\xi(t)\rangle^* = \delta(t-t')$). Linearizing the equations (1) and (2) around a stable solutions $p_o$ and $\psi_o$ (without considering thermal noise) allows us to study the transient behaviour of the synchronized state, finding two solutions for the decay rate $\lambda$:

$$ \lambda_{1,2} = \Gamma_p \left[ 1 \pm \sqrt{1 - \frac{\varepsilon}{\varepsilon_c}} \right] \quad (3) $$

Here, $\varepsilon = I_{RF}/I_{DC}$ and $\varepsilon_c = \Gamma_p^2/(2Np_o^2\sin\psi_s P_x \Gamma_l |B|/A)$ (see Supplementary material for the definition of the parameters). When $\varepsilon > \varepsilon_c$, $\lambda$ is complex with a real part given by $\Gamma_p$, that is the decay rate to the phase locked state and an imaginary part that describes an oscillatory approach to the phase-locked state with a frequency given by:

$$ \Omega_t = \Gamma_p \sqrt{\varepsilon/\varepsilon_c - 1} \quad (4) $$

This is in agreement with Zhou et al.$^{19}$, where they found for out of plane (OP) magnetized STNO's that the phase approaches its locked state exponentially and oscillating above a certain critical reference current. These oscillations lead to sidebands of frequency $\Omega_t$ at both sides of the emission peak of the STNO as shown in Fig. 1d.

Before discussing in more detail the oscillatory transient, we first will provide an expression for the phase noise in the synchronized state. By taking into account the thermal noise, and supposing that phase and power are small deviation from the equilibrium we can calculate from eqs 1, 2 the power spectral density (or single side-band) $S_{\delta\phi}$ of the phase fluctuations of the synchronized state:

$$ S_{\delta\phi} = 2\pi\Delta f_o \frac{\left(\frac{\Gamma_p}{\pi}\right)^2 (1 + \nu^2) + f^2}{\left[-\left(\frac{\Gamma_p}{2\pi}\right)^2 \frac{\varepsilon}{\varepsilon_c} + f^2\right]^2 + 4\Gamma_p^2 f^2} \quad (5) $$

Here $\nu = Np_o/\Gamma_p$ is the dimensionless nonlinear frequency shift and $\Delta f_o$ is the free running linear linewidth. Eq. 5 is plotted in Fig. 1a with the parameters calculated from the analytical model (see Supplementary Material) for the free running state with a bias current $J_{DC} = 50 \times 10^{10} A/m^2$, which leads to an IPP stable precession mode around 4.7 GHz and a $\Delta f_o = 50$ MHz. In this configuration the system has a coefficient of nonlinear frequency shift $N = -3.16 \times 10^{10}$ rad/sec, a damping rate of the power fluctuations $\Gamma_p = 666$ Mrad/sec, a normalized dimensionless nonlinear frequency shift parameter $|\nu| = 16$, and $\varepsilon_c = 0.025$. Since this value of $\varepsilon_c$ for these uniform IP STNO's is small compared to $\varepsilon$ (~0.1 or higher) the phase locking always takes place via an oscillatory transient.

The characterization of the phase noise properties by the PSD in the Fourier space has the advantage that its inverse power law dependence on frequency $\text{PSD} \sim 1/f^\alpha$ provides information about underlying noise processes. The model predicts a crossover from $1/f^2$ to $1/f^\alpha$ with increasing reference current, with the two limit cases:

$$ S_{\delta\phi} = \begin{cases} f \gg f_{\text{roll-off}} : 2\pi\Delta f_o \left( \frac{1}{f^2 + 4\Gamma_p^2} \right) \\ f \ll f_{\text{roll-off}} : 2\pi\Delta f_o \frac{4(1 + \nu^2)}{(\Gamma_p/\pi)^2} \left( \frac{\varepsilon_c}{\varepsilon} \right)^2 \end{cases} \quad (6) $$

As already shown experimentally$^6$ for IPP STNO's, the free running oscillator ($\varepsilon = 0$) shows a $1/f^2$ dependence associated to a random walk of the phase (blue line, Fig. 1a). This behaviour is modified when applying a reference
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Phase noise calculated from the analytic model (a) and numerical simulations at 50 K for the phase (d) and amplitude (c) noise. Notice the peaks appearing ~200–300 MHz. (d) PSD of the signal as a function of ε for J<sub>DC</sub> = −50 × 10<sup>10</sup> A/m<sup>2</sup> and schematics of the oscillator (inset). (e) Linewidth vs. ε for low current J<sub>DC</sub> = −40 × 10<sup>10</sup> A/m<sup>2</sup> (red dots) and medium current J<sub>DC</sub> = −50 × 10<sup>10</sup> A/m<sup>2</sup> (blue dots). For ε > 0.3 (medium J<sub>DC</sub> regime), the linewidth falls below the resolution of the technique. (f) Sideband frequency for both low current regime (red set) and medium current regime (blue set) extracted from the PSD (stars), and from the peaks on the phase noise (Fig. 1b) (open circle). The solid line corresponds to the analytical model developed for the IPP geometry (eq. 4).</figcaption></figure>

current I<sub>RF</sub> at 2f: Even for a low external force (ε = 0.1, yellow line, Fig. 1a) below the roll of frequency f<sub>roll off</sub> ~1/Γ<sub>p</sub> down to the lowest (calculated) offset frequency the phase noise is constant. This corresponds to fluctuations of the phase around its locked value. The Eq. (6) shows that the phase noise level in this region can be decreased upon increasing the reference current ε. Nevertheless the phase noise falls as 1/ε² hence for higher reference currents the noise is reduced, but at a much smaller rate for further increasing. The lowest achievable phase noise is limited firstly by the STNO's voltage breakdown and secondly because a large reference current goes beyond the injection locking assumption, which is considered as a weak perturbation. Also the intrinsic parameters of the STNO as f<sub>o</sub>, Γ<sub>p</sub> and ν play a role in the lowest noise level suggesting that it can be improved by stack engineering. Last, above ε = 0.1 there is a peak around f<sub>roll off</sub> that is related with the oscillatory relaxation mechanism<sup>19,22-24</sup>, as will be discussed in the next section.

**Macrospin Analysis.** The S<sub>φφ</sub> at 50 K extracted from the numerical time integration of the LLG equation is shown in Fig. 1b. The corresponding evolution of the power spectral density of the *my* component of the magnetization with ε and its linewidth are displayed in Fig. 1d and e respectively. Both phase and amplitude noise (Fig. 1b and c) decrease with the reference current and a clear crossover from a 1/f<sup>2</sup> (random walk) to a 1/f<sup>0</sup> (white noise) is seen on the phase noise upon increasing ε.

Before addressing the phase noise level in comparison to the analytic results and the linewidth, we now discuss the PSD of the *m*<sub>*y*</sub> component of the magnetization, Fig. 1d. The peak of the free running state becomes very narrow and two symmetric sidebands appear upon increasing the reference current. These sidebands are also visible on S<sub>φφ</sub> extracted from the simulations (Fig. 1b), that shows a peak around f<sub>roll-off</sub> whose frequency depends on ε. Moreover, S<sub>φφ</sub> calculated with the parameters listed above (eq. 5, Fig. 1a) also shows the peaks associated to the sideband frequencies. Figure 1f shows the frequency of these sidebands extracted from the PSD of the *m*<sub>*y*</sub> component of the magnetization and from the numeric S<sub>φφ</sub> for two different bias currents J<sub>DC</sub> = −40 × 10<sup>10</sup> A/m<sup>2</sup> (red symbols) and J<sub>DC</sub> = −50 × 10<sup>10</sup> A/m<sup>2</sup> (blue symbols). The full line represents the model (eq. 4) for both bias currents. This comparison confirms that the peaks of the phase noise and the sidebands have the same physical origin arising from the oscillatory approach of the transient state. Furthermore the comparison supports the analytic model.

In the following we discuss the phase noise level of the numerical results that show two contributions to the locked state. The first one, as it was discussed for the analytical description, are phase fluctuations around the stable phase which is given by the external source plus a constant phase shift. The second contribution comes from the phase slips<sup>25</sup>, not considered in the analytical model but that are present in the numerical calculation. To understand their contribution to the phase noise and linewidth we extracted the phase from simulated time traces for different reference current values (Fig. 2a). The phase trace shows a drift in time, together with the appearance of the phase slips, which become well defined upon increasing the reference current. As can be seen the number of phase slips per sampling period (~40 µs) decreases with increasing ε. These phase slips are responsible for the
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Simulated (a) and experimental (c) phase temporal traces. Inset: detail of 5µs segments of the temporal trace. The phase slips decrease in number with increasing ε disappearing at ε = 0.3. (b) Phase noise analysis on the simulated temporal trace segments (inset) corresponding to no phase slips (black), 1 phase slip (red) and 2 phase slips (green). (d) Phase noise analysis of the experimental time trace from 3 s segments with 0, 1, 2, 10 phase slips.</figcaption></figure>

1/f² contribution of the phase noise at low offset frequencies. To demonstrate this, we compare the phase noise extracted from different 10 µs segments of the total 40 µs phase trace that contain respectively none, one or two phase slips. In Fig. 2b it is clearly seen that in presence of phase slips the phase noise has a 1/f² dependence at low frequency f < f<sub>roll off</sub>, while in absence of phase slips there is a constant phase noise level of around -100 dBc/Hz.

From the numerical analysis we can see that the drastic decrease in the linewidth with ε (Fig. 1e) can be related to the decreasing number of phase slips. Particularly, when the phase noise flattens for ε > 0.3 (ε > 2.75) for J<sub>DC</sub> = -50 × 10<sup>10</sup> A/m<sup>2</sup> (J<sub>DC</sub> = -40 × 10<sup>10</sup> A/m<sup>2</sup>), the phase slips are absent in the phase trace and the linewidth falls under the resolution limit of the numerical calculation (20 kHz). This fact shows that the so called “pure” synchronization<sup>20,26</sup> is due the absence of phase slips, where the STNO would reduce its linewidth ideally to the one of the reference source. In the case discussed here, this means that for values larger than ε = 0.3 (ε = 2.75) torques from the reference current on the magnetization are strong enough to stabilize the phase around a single value and the remaining noise is given by the one discussed within the analytical model, describing damped oscillations around the stable phase, for noise frequencies f < f<sub>roll off</sub>. We point out here that the absence of the phase slips depends on the observation time, i.e. the length of the temporal trace: longer observation times increase the probability of phase slips. These results evidence that even if the system is in the frequency locking regime, higher values of reference current are needed to achieve full linewidth reduction by phase locking. This observation is in agreement with Lebrun *et al.*<sup>26</sup> that highlighted the difference between the reported “frequency locked state<sup>18,27</sup>” and “pure phase locked state” in absence of phase slips for vortex oscillators with a free running frequency of ~200 MHz and ~100 kHz linewidth (free running).

**Experiment.** The analytical and numerical results explain qualitatively the experimental observations on the injection locking by a reference current. We present here results for a device with an autonomous, i.e. free running regime characterized by a free running frequency of f₀ = 7.5 GHz for a bias current I<sub>DC</sub> = -1.6 mA and an applied in plane field of 350 Oe, with a linewidth of 55 MHz. The PSD map of the output power for the STNO frequency f as a function of the source frequency f<sub>ext</sub> is shown in the Fig. 3a for ε = J<sub>RF</sub> / J<sub>DC</sub> = 0.7. In Fig. 3b it is clearly seen that for increasing reference current ε the linewidth reaches a 10 × reduction (8 MHz with a 1 MHz resolution bandwidth). The amplitude noise shows a 1/f⁰ behaviour both for locking (Fig. 3c, grey dashed line) and the free running state (Fig. 3c, black full line), characteristic of white noise fluctuations of the amplitude around its stable value. The experimental plots of the phase noise in the locking (Fig. 3c, red dashed line) and the free running state
<figure><img src="image_5.png" /><figcaption>**Figure 3.** PSD map of the output power at $I_{RF} = 1.12$ mA (a). Linewidth vs $\varepsilon = I_{RF}/I_{DC}$ (b) and amplitude and phase noise from the experiment (c) for the non-locked state, (continuous line, $\varepsilon = 0$) and locked state (dashed line, $\varepsilon = 0.7$). Notice that the injection locking mechanism is efficient to reduce the noise level (green arrows).</figcaption></figure>

(red full line) show that the injection locking mechanism is efficient to reduce the phase noise by 20 dBc/Hz with respect to the non-locked state. Both phase noise plots display a 1/f² behaviour, however the origin is different. In the free running state it results from a random walk of the phase, while in the locking state it is due to the phase slips as explained in the macrospin analysis. Our experiments show that for the maximum applied reference current, even though the oscillator is synchronized with the external source, the emission linewidth remains broad. The phase noise decreases but it does not reach the constant level for which the linewidth is expected to reduce to the source noise. This was not observed in our experiments because the voltage breakdown of the samples did not allow to continue increasing the reference current preventing the STNO from achieving a pure phase locked state. For the same reason, we were not able to observe sidebands in the experimental voltage output. Nevertheless, we have extracted the phase noise from shorter 3 µs segments of the temporal trace²⁸ of the phase (Fig. 2c), that include either 0, 1, 2 or 10 phase slips. As can be seen in absence of phase slips the phase noise is flat in a certain range of offset frequencies ~5–100 MHz in Fig. 2d. This demonstrates the lowest phase noise level that can be reached for the in-plane STNO, when phase slips would be completely suppressed. Note that this level of –100 dBc/Hz is the same for the analytic model (Fig. 1a) and simulations (Fig. 1b).

The large free running linewidth of ~50 MHz of our STNO's could be a drawback to linewidth reduction as discussed by Hamadeth et al.²⁰ where by decreasing 7 times the free running linewidth, the linewidth reduction in the locked regime goes from only 10 to 10⁵. Indeed, this is also witnessed in our numerical simulations, where the smaller free running linewidth at higher bias current ($J_{DC} = -50 \times 10^{10}$ A/m²) leads to a linewidth reduction for significantly lower reference currents.

**Perspectives.** We have studied the synchronization mechanism of a uniform IP magnetized STNO under thermal noise. The synchronization of these devices was demonstrated in several experiments, however no more than a 10 times reduction of linewidth was achieved. This is explained by numerical simulations including thermal noise. While the STNO can be synchronized by moderate reference currents, higher reference currents are needed for full linewidth reduction. With increasing reference current the number of phase slips is reduced resulting in a crossover from 1/f² to 1/f⁰ behaviour in the phase noise when the phase slips are suppressed. The simulations also show that it is possible to achieve linewidth reduction for lower reference currents by increasing the bias current of the oscillator. This study shows the role of the intrinsic parameters of an STNO on the phase noise level, which will be relevant for the design of STNO configurations of appropriate performances for microwave applications in the gigahertz range.

**Methods**

We performed macrospin simulations for the in-plane precession (IPP) mode of an in-plane magnetized polarizer and free layer magnetic tunnel junction MTJ, using a solver for the Landau-Lifshitz-Gilbert equation and taking
into account the damping like torque term (the field like term was neglected in this work, see Supplementary material). The simulation parameters are the following: free layer of size 90 × 80 × 3.9 nm³; spontaneous magnetization M_s = 1000 kA/m, damping parameter α = 0.02 and zero magneto-crystalline anisotropy. The polarizer is aligned in the plane at 165° from the free layer magnetization equilibrium position and a spin-polarization η = 0.37 is supposed. A static magnetic field of 40 mT was applied along the in plane easy axis (Ox). The continuous current density was set to J_Dc = -50.10¹⁰ A/m², leading to an IPP stable precession mode with f ~ 4.7 GHz. A white Gaussian thermal noise field was added, corresponding to 10 K, 20 K, 50 K, and 100 K¹⁸. The frequency of the RF current was set to two times the free running STNO frequency 2f₀ = 9,5396 GHz, which corresponds to the centre of the locking range. The phase and amplitude noise in the synchronized state as a function of the reference current were extracted from the simulated temporal traces of the mγ component of the magnetization (in-plane magnetization along the short axis of the pillar) using the Hilbert transform method²⁹,³⁰ which allows the reconstruction of an analytic signal from the voltage output:

$$V = V_0(1 + \delta a)\cos(2\pi ft + \delta \Phi)$$

The experimental studies were realized on the same type of devices presented in refs⁶,¹⁸,³¹, which are in plane magnetized MTJ, having a stack composition of IrMn/CoFeB/Ru/CoFeB/MgO/CoFe/CoFeB and nominal resistance area product RA = 1 Ωµm². The synchronization experiment was done varying the reference current frequency around two times the free running frequency (2f₀) of the oscillator, from 14 GHz to 16 GHz, and the source power was varied from -15 dBm to -5 dBm (corresponding to a reference current of ~0.3 to 1.3 mA), just before the sample starts to show signs of degradation. A detailed description of the experiment is available in ref.¹⁸. The temporal traces were measured using a single shot oscilloscope⁴,¹⁸,³², and amplitude and phase noise were extracted using the same protocol as for the simulated data.

## References

1. Slonczewsky, J. C. Current-driven excitation of magnetic multilayers. J. Magn. Magn. Mater 195, L1 (1996).

2. Berger, L. Emission of spin waves by a magnetic multilayer traversed by a current. Phys. Rev. B 54, 9353 (1996).

3. Tsoi., M. et al. Excitation of a magnetic multilayer by an electric current. Phys. Rev. Lett. 80, 4241-4284 (1998).

4. Keller, M. W., Kos, A. B., Silva, T. J., Rippard, W. H. & Pufall, M. R. Time domain measurement of phase noise in a spin torque oscillator. Appl. Phys. Lett. 94, 193105 (2009).

5. Keller, M. W., Pufall, M. R., Rippard, W. H. & Silva, T. J. Nonwhite frequency noise in spin torque oscillators and its effect on spectral linewidth. Phys. Rev. B 82, 054416 (2010).

6. Quinsat, M. et al. Amplitude and phase noise of magnetic tunnel junction oscillators. Appl. Phys. Lett. 97, 182507 (2010).

7. Grimaldi, E. et al. Response to noise of a vortex based spin transfer nano-oscillator. Phys. Rev. B 89, 104404 (2014).

8. Silva, T. J. & Keller, M. W. Theory of thermally induced phase noise in spin torque oscillators for a high-symmetry case. IEEE Trans. Magn. 46, 3555 (2010).

9. Sierra, J. F. et al. Influence of thermal fluctuations on the emission linewidth in MgO-based spin transfer oscillators. Appl. Phys. Lett. 101, 062407 (2012).

10. Georges, G. et al. Coupling efficiency for phase locking of a spin transfer nano-oscillator to a microwave current. Phys. Rev. Lett. 101, 017201 (2008).

11. Grollier, J., Cros, V. and Fert, A. Synchronization of spin-transfer oscillators driven by stimulated microwave currents. Phys. Rev. B 73, 060409(R) (2006).

12. Urazhdin, S., Tiberkevich, V. S. & A. Slavin, A. N. Parametric excitation of a magnetic nanocontact by a microwave field. Phys. Rev. Lett. 105, 237204 (2010).

13. Belanovsky, A. D. et al. Numerical and analytical investigation of the synchronization of dipolarly coupled vortex spin-torque nano-oscillators. Appl. Phys. Lett. 103, 122405 (2013).

14. Chen, H.-H. et al. Synchronization of spin torque nano-oscillators through dipolar interactions. J. App. Phys 115, 134309 (2014).

15. Chen, H.-H. et al. Phase Locking of spin-torque nano-oscillator pairs with magnetic dipolar coupling. Phys. Rev. B 93, 224410 (2016).

16. Sani, S. et al. Mutually synchronized bottom-up multi-nanocontact spin-torque oscillators. Nat Commun. 4, 2731 (2013).

17. Rippard, W. H. et al. Injection Locking and Phase Control of Spin Transfer Nano-oscillators. Phys. Rev. Lett 95, 067203 (2005).

18. Quinsat, M. et al. Injection locking of tunnel junction oscillators to a microwave current. Appl. Phys. Lett. 98, 182503 (2011).

19. Zhou, Y. et al. Oscillatory transient regime in the forced dynamics of a nonlinear auto oscillator. Phys. Rev. B 82, 012408 (2010).

20. Hamadeh, A. et al. Perfect and robust phase-locking of a spin transfer vortex nano-oscillator to an external microwave source. Appl. Phys. Lett. 104, 022408 (2014).

21. Slavin, A. & Tiberkevich, V. Nonlinear auto-oscillator theory of microwave generation by spin polarized current. IEEE Trans. Magn. 45, 4 (2009).

22. Bianchini, L. et al. Direct experimental measurement of phase-amplitude coupling in spin torque oscillators. Appl. Phys. Lett. 97, 032502 (2010).

23. d'Aquino, M., Serpico, C., Bertotti, G., Bonin, R. & Mayergoyz, I. D. Micromagnetic analysis of injection locking in spin-transfer nano-oscillators. Phys. Rev. B 84, 064415 (2010).

24. Finocchio, G., Carpentieri, M., Giordano, A. & Azzerboni, B. Non-Adlerian phase slip and nonstationary synchronization of spin-torque oscillators to a microwave source. Phys. Rev. B 86, 014438 (2012).

25. Pikovsky, A., Rosenblum, M. and Kurths, J. Synchronization: A universal concept in Nonlinear Sciences (Cambridge University Press, Cambridge 2001).

26. Lebrun, R. et al. Understanding of phase noise squeezing under fractional synchronization of a nonlinear spin transfer vortex oscillator. Phys. Rev. Lett. 115, 017201 (2015).

27. Demidov, V. E. et al. Synchronization of spin Hall nano-oscillators to external microwave signals. Nat. Commun. 5, 3179 (2014).

28. Houssameddine, D. et al. Not all the jumps on the experimental phase temporal trace (Fig. 3a) are phase slips. Only the ±2π jumps are compatible with phase slips, larger ones are related to extinctions in the oscillator, as seen in. Phys. Rev. Lett 102, 257202 (2009).

29. Picinbono, B. On instantaneous amplitude and phase of signals. IEEE Trans. Signal Process 45, 552 (1997).

30. Russek, S. E., Kaka, S., Rippard, W. H., Puffall, M. R. & Silva, T. Finite-temperature modeling of nanoscale spin-transfer oscillators. Phys. Rev. B 71, 104425 (2005).

31. Houssameddine, D. et al. Spin transfer induced coherent microwave emission with large power from nanoscale MgO tunnel junctions. Appl. Phys. Lett. 93, 022505 (2008).

32. Houssameddine, D., Ebels, U. & Deny, B. Temporal Coherence of MgO Based Magnetic Tunnel Junction Spin Torque Oscillators. Phys. Rev. Lett 102, 257202 (2009).
**Acknowledgements**

This work was supported in part by the French National Research Agency (ANR) under contract N° 2011 Nano 01607 (SPINNOVA) and in part by the EC under the FP7 program N° 317950 MOSAIC.

**Author Contributions**

M.T. and U.E. designed the experiment, analyzed the macrospin and experimental data and wrote the manuscript. L.D.B.P. performed the macrospin simulation, the analytical model and contributed to the writing of the manuscript. M.T. contributed to the analytic model and did the experiment. B.L. contributed to the analytic model. J.H. and C.D. contributed to the experiment and data analysis. M.-C.C., J.A.K., D.M. and A.Z. did the sample and microfabrication.

**Additional Information**

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-017-18969-5.

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018