# Spin-polarized Shapiro steps and spin-precession-assisted multiple Andreev reflection

###### Abstract

We investigate the charge and spin transport of a voltage-biased superconducting point contact coupled to a nanomagnet. The magnetization of the nanomagnet is assumed to precess with the Larmor frequency, , due to ferromagnetic resonance. The interplay between the ac Josephson current and the magnetization dynamics leads to spin-polarized Shapiro steps at voltages for and the subharmonic steps with are a consequence of multiple Andreev reflection (MAR). Moreover, the spin-precession-assisted MAR generates quasiparticle scattering amplitudes that, due to interference, lead to current-voltage characteristics of the dc charge and spin currents with subharmonic gap structures displaying an even-odd effect.

###### pacs:

75.76.+j, 74.50.+r, 75.50.Xx, 75.78.-nIntroduction. Control of the electron’s spin degree of freedom has led to a number of spintronics applications. tserkovnyak2005 () These applications often rely on using the spontaneously broken spin-rotation symmetry of a ferromagnet to act as a reference for the spin of itinerant electrons. In contrast to ferromagnets, the electron-electron interaction in a BCS superconductor leads to a spinless order parameter. While these two states of matter usually are incompatible in bulk materials, they can be combined in nanoscale junctions to create new phenomena. eschrig2011 (); bergeret2005 (); buzdin2005 () The interplay between superconductivity and ferromagnetism in such hybrid superconducting junctions allows for probing the spin degree of freedom as well as for coherently controlling spin-polarized currents. fogelstrom2000 (); zhao2004 (); andersson2002 (); zhao2007 (); cuevas2001 (); bobkova2006 (); keizer2006 (); senapati2011 (); khaire2010 (); huertas-hernando2002 (); shomali2011 () In addition, superconducting junctions coupled with magnetization dynamics have been studied volkov2009 (); zhu2003 (); michelsen2008 (). The lowest-energy magnetization excitation can be accessed by the application of an external magnetic field, which starts a precession of the magnetization around the direction of the field bell2008 (); morten2008 (); skadsem2011 (); richard2012 (). This ferromagnetic resonance (FMR) mode can also be achieved by a coupling between the magnetization and the Josephson effect: the ac Josephson current generates an oscillating magnetic field that resonantly excites the magnetization precession thirion2003 (); petkovic2009 (); barnes2011 ().

In Ref. petkovic2009 (), the coupling between a spatially dependent order parameter and the magnetization dynamics was studied in a tunnel junction whose width, , was larger than the superconducting coherence length, , and this coupling was detected as a rectification of the ac Josephson charge current. However, in the regime , which is the case for a superconducting quantum point contact (SQPC) muller1992 (); scheer1997 (), the spin-rotation symmetry of the superconducting order parameter may be broken locally and one may anticipate that spin-dependent modifications of the superconducting correlations will influence the transport properties teber2010 (); holmqvist2011 (). In this Letter, we study the coupling of the ac Josephson effect and magnetization dynamics of an SQPC containing a nanomagnet, e.g. a magnetic impurity thirion2003 () or a single-molecule magnet bogani2008 (); heersche2006 (); jo2006 (); kasumov2005 (); kahle2012 (). We show that the interplay between the Josephson and magnetization oscillations create Shapiro-like resonances that lead to a rectified spin current. In addition, we find a rich subgap structure of the current-voltage characteristics of the dc charge and spin currents displaying features related to the emission or absorption of energy quanta corresponding to the precession frequency. However, the features associated with an odd number of Andreev reflections are suppressed due to interference.

The transport properties of nonmagnetic SQPCs can be understood in terms of Andreev reflection. An applied bias voltage leads to the occurrence of multiple MAR octavio1983 () and the ac Josephson effect, which is characterized by the Josephson frequency, . The ac current of SQPCs also includes higher harmonics of the Josephson frequency, i.e. (setting ) bratus1995 (); averin1995 (); cuevas1996 (). In addition, MAR processes produce a subharmonic gap structure of the dc current that can be seen in the current-voltage characteristics as steps at voltages , The ac Josephson effect can be probed by the application of microwave radiation as their interplay leads to Shapiro steps and a rich subgap structure due to photon-assisted MAR cuevas2002 (); uzawa2005 (); chauvin2006 (). In the case of SQPCs with magnetic interfaces, the transport properties are modified due to spin-dependent effects fogelstrom2000 (); andersson2002 (); cuevas2001 (). E.g., a magnetic interface may produce different transmission probabilities for the spin-up and spin-down bands, and , which may lead to the occurrence of spin currents bobkova2006 (); zhao2007 ().

Here, we consider a point contact consisting of two superconducting leads coupled over a nanomagnet. The left (L) and right (R) leads are assumed to consist of s-wave superconductors in the clean limit cuevas1996 (). We further assume that the magnetization of the nanomagnet can be treated classically within the macrospin model and hence be described by single spin vector, tserkovnyak2005 (). Under FMR conditions, external dc and rf fields are applied to create an effective magnetic field, , in which the spin precesses with the Larmor frequency , where is the gyromagnetic ratio. This effective field also includes any effects of the resonant coupling between the ac Josephson current and the magnetization dynamics. The nanomagnet’s spin, whose dynamics can be described by the Landau-Lifshitz-Gilbert equation LLG (), generates a time-dependent exchange field that is felt by the quasiparticles tunneling across the junction and can be incorporated into a phenomenological tunnel Hamiltonian that allows for spin flip scattering () which is accompanied by absorption (emission) of an energy quantum by the quasiparticle (see Ref. holmqvist2011 () for details).

The transport characteristics depend on the time-dependent superconducting phase difference across the contact, , where is the initial phase difference, in addition to the Larmor frequency. The spin traces out an in-plane angle described by in analogy with . The scattering resulting from the MARs combined with the precession-induced spin flips transfers quasiparticles into sidebands with energies , where and . The charge () and spin (-component) currents in lead can then be written in terms of all harmonics as

(1) |

Model. The currents are obtained using nonequilibrium Green’s function techniques and the calculation is based on the Hamiltonian approach cuevas1996 () in which the leads are described by Keldysh Green’s functions in the quasiclassical approximation eliashberg1972 () and the interface containing the nanomagnet is treated as a strong impurity serene1983 (). The effects of are accounted for in the standard way by a transformation of the unperturbed lead Green’s functions, , where ”” denotes a matrix in Keldysh space, that is expressed as in Nambu space (””). The voltage-biased lead Green’s functions, , are then expressed in terms of the equilibrium ones, , as where refers to one of the retarded (), advanced () or Keldysh () components. The equilibrium Green’s functions have a simple Fourier transformation given by . Since the leads consist of s-wave superconductors, the retarded component of is given by , where , , . The advanced component is while the Keldysh Green’s function is given by .

The magnetic interface is characterized by a tunnel Hamiltonian, , in which the hopping amplitudes, , depend on the instantaneous direction of the spin, . Within the quasiclassical approximation, the Fermi surface limits of the hopping amplitudes can be expressed as . We can define the spin-up (spin-down) transmission probability as , where and is the angle between and .

The hopping elements’ influence on the lead Green’s functions is captured by the t-matrix equation , where . The time dependence of the superconducting phase of the Green’s functions can be transferred to the hopping elements by the transformation according to . The -matrix equation then becomes an algebraic equation in energy space and

(2) |

By only considering currents on the left side of the junction (thereby dropping the index ), one can write the Fourier components as and the current components of Eq. (1) as

(3) |

where , , for . The -matrix components can be found from the equation where the retarded components are given by

(6) | |||||

(7) |

where we have defined , and if . The advanced and Keldysh matrices are similarly defined.

Results. – The dc charge and spin currents of Eq. (1) consist primarily of contributions from the components and can be described as background currents in the spirit of Ref. cuevas2002 (). There is also the possibility of Shapiro-like resonances between the Josephson and the Larmor frequencies at voltages , . The Larmor precession produces a spin structure of the matrices in which the index corresponds to a spin-raising (spin-lowering) matrix. This spin structure causes the Shapiro resonances to give contributions to the spin-polarized current only and thus result in a Shapiro spin current given by . This behavior is in contrast to that of a weak link coupled to a nanomagnet whose magnetization performs a precessional motion. In that case, the resulting time-dependent magnetic field produces Shapiro-type resonances that give contributions to the charge current cai2010 ().

Figure 1 shows the current-voltage characteristics as well as the conductance for the charge background current, . As can be seen, the most prominent feature is the subharmonic gap structure appearing at voltages , where This structure can clearly be seen in the conductance where it is displayed as distinct satellite peaks for junctions in the tunnel limit. The subharmonic gap structure is similar to that of the current-voltage characteristics generated by photon-assisted MAR which shows features due to the absorption or emission of photons cuevas2002 (). However, in the case of MAR processes influenced by Larmor precession, a quasiparticle can only absorb or emit one quantum of since these processes are accompanied by a spin flip (see e.g. Fig. 2 for the absorption processes with and ). As a result, there is only one satellite peak on each side of the feature corresponding to a tunnel process of order . In addition, only the processes for which is an odd number display side peaks. One can show that this suppression is due to destructive interference of the Andreev reflection amplitudes corresponding to the side peaks of processes with Defining , the height of the current steps at can be approximated in the tunnel limit by , where is the Heaviside step function, and and are the complete elliptic integrals of the first and second kinds. In the limit , one obtains and the heights of the current steps at are thus .

The component of the background spin current, , is shown in Fig. 3 (the in-plane components are zero). For transmission probabilities and , our results reproduce those of Ref. zhao2007 () for zero spin mixing. In the case of nonzero Larmor frequency, the spin current is finite even for , as can be seen in panel (b). In the tunnel limit, the spin current can be divided into a spin-filter current and a spin-pump current, . The component of the spin-filter current is given by for zhao2007 (). The spin-pump current can correspondingly be approximated by . For and , the total spin-pump current is given by if and is hence and does not depend on the bias voltage.

Now, we turn to the Shapiro-like resonances in the spin current. Since the Larmor precession requires to take the values for the Shapiro resonance condition, is spin polarized in the plane (see Ref. holmqvist2011 () for a detailed analysis of the matrices’ spin structure). Defining where , one can write the th component as , where and the transformation rotates the vector around the direction of the external magnetic field, , through an angle . The existence of a dc current that is spin polarized in the plane implies that the rotational symmetry around the axis is broken and that the magnitude of the Shapiro currents depend on the initial angle of nanomagnet’s magnetization direction. This behavior is analogous to the dependence on the initial superconducting phase difference which also is present in the behavior of microwave-irradiated Josephson junctions cuevas2002 (). The Shapiro contributions strongly depend on the transmission probability as can be seen in Fig. 4.

The values in Fig. 1 are chosen for clarity but may be smaller in practice. Consider an SQCP consisting of Nb, whose superconducting gap is meV, containing a nanomagnet that under FMR conditions reaches precession frequencies of GHz which is well below the critical magnetic field, , of Nb bell2008 (). Then, typically and the temperature is restricted to mK. Increasing closer to the critical field value decreases and hence increases the ratio which allows for a better resolution of the subgap structures of and . Alternatively, using the ac Josephson current to resonantly excite , values corresponding to an effective magnetic field can be achieved petkovic2009 (); barnes2011 (); thirion2003 (). Detection of the spin-polarized Shapiro currents would then be a measurement of the FMR frequency of the nanomagnet.

Conclusions. – In conclusion, we have calculated the dc charge and spin currents through a voltage-biased superconducting point contact coupled to the spin of a nanomagnet that under ferromagnetic resonance conditions precesses with the Larmor frequency . We have shown that coherent multiple Andreev reflection in the presence of the Larmor dynamics leads to a subharmonic gap structure of the dc charge and spin currents with features at originating from the absorption or emission of a precession quantum. However, destructive interference suppresses these features for even values of . In addition, the Larmor dynamics combined with broken spin-rotation symmetry generates Shapiro-like resonances resulting in an additional spin current polarized in the plane.

Acknowledgments. – We wish to thank J.C. Cuevas for useful discussions. C.H. and W.B. were supported by the DFG and SFB 767. M.F. acknowledges support from the Swedish Research Council (VR).

## References

- (1) Y. Tserkovnyak, A. Brataas, G.E.W. Bauer, B.I. Halperin, Rev. Mod. Phys. 77, 1375 (2005).
- (2) A.I. Buzdin, Rev. Mod. Phys. 77, 935 (2005).
- (3) F.S. Bergeret, A.F. Volkov, and K.B. Efetov, Rev. Mod. Phys. 77, 1321 (2005).
- (4) M. Eschrig, Physics Today 64, 1, 43 (2011).
- (5) M. Fogelström, Phys. Rev. B 62, 11812 (2000).
- (6) J.C. Cuevas and M. Fogelström, Phys. Rev. B 64, 104502 (2001).
- (7) D. Huertas-Hernando, Yu.V. Nazarov, and W. Belzig, Phys. Rev. Lett. 88, 047003 (2002).
- (8) M. Andersson, J.C. Cuevas, and M. Fogelström, Physica C 367, 117 (2002).
- (9) E. Zhao, T. Löfwander, and J.A. Sauls, Rev. B 70, 134510 (2004).
- (10) R.S. Keizer, S.T.B. Goennenwein, T.M. Klapwijk, G. Miao, G. Xiao, and A. Gupta, Nature (London) 439, 825 (2006).
- (11) I.V. Bobkova, Phys. Rev. B 73, 012506 (2006); I.V. Bobkova and A.M. Bobkov, ibid. 76, 094517 (2007).
- (12) E. Zhao and J.A. Sauls, Phys. Rev. Lett. 98, 206601 (2007).
- (13) T.S. Khaire, M.A. Khasawneh, W.P. Pratt, and N.O. Birge, Phys. Rev. Lett. 104, 137002 (2010).
- (14) K. Senapati, M.G. Blamire, and Z.H. Barber, Nature Materials 10, 849 (2011).
- (15) Z. Shomali, M. Zareyan, and W. Belzig, New J. Phys. 13, 083033 (2011).
- (16) J.X. Zhu and A.V. Balatsky, Phys. Rev. B 67, 174505 (2003); J.X. Zhu, Z. Nussinov, A. Shnirman, and A.V. Balatsky, Phys. Rev. Lett. 92, 107001 (2004).
- (17) J. Michelsen, V. S. Shumeiko, and G. Wendin, Phys. Rev. B 77, 184506 (2008).
- (18) A.F. Volkov and K.B. Efetov, Phys. Rev. Lett. 103, 037003 (2009).
- (19) C. Bell, S. Milikisyants, M. Huber, and J. Aarts, Phys. Rev. Lett. 100, 047002 (2008).
- (20) J.P. Morten, A. Brataas, G.E.W. Bauer, W. Belzig, and Y. Tserkovnyak, EPL (Europhysics Letters) 84, 57008 (2008).
- (21) H.J. Skadsem, A. Brataas, J. Martinek, and Y. Tserkovnyak, Phys. Rev. B 84, 104420 (2011).
- (22) C. Richard, M. Houzet, and J.S. Meyer, Phys. Rev. Lett. 109, 057002 (2012).
- (23) C. Thirion, W. Wernsdorfer, and D. Mailly, Nature Materials 2, 524 (2003).
- (24) I. Petković, M. Aprili, S. E. Barnes, F. Beuneu, and S. Maekawa, Phys. Rev. B 80, 220502 (2009).
- (25) S.E. Barnes, M. Aprili, I. Petković, S. Maekawa, Superconductor Science and Technology 24, 024020 (2011).
- (26) C.J. Muller, J.M. van Ruitenbeek, and L.J. de Jongh, Phys. Rev. Lett. 69, 140 (1992).
- (27) E. Scheer, N. Agrait, J.C. Cuevas, A. Levy-Yeyati, B. Ludolph, A. Martin-Rodero, G.R. Bollinger, J.M. van Ruitenbeek, and C. Urbina, Nature (London) 394, 154 (1998); E. Scheer, W Belzig, Y. Naveh, M.H. Devoret, D. Esteve, and C. Urbina, Phys. Rev. Lett. 86, 284 (2001).
- (28) S. Teber, C. Holmqvist, and M. Fogelström, Phys. Rev. B 81, 174503 (2010).
- (29) C. Holmqvist, S. Teber, and M. Fogelström, Phys. Rev. B 83, 104521 (2011); C. Holmqvist, W. Belzig, and M. Fogelström, ibid. 86, 054519 (2012).
- (30) A. Yu. Kasumov, K. Tsukagoshi, M. Kawamura, T. Kobayashi, Y. Aoyagi, K. Senba, T. Kodama, H. Nishikawa, I. Ikemoto, K. Kikuchi, V. T. Volkov, Yu. A. Kasumov, R. Deblock, S. Guéron, and H. Bouchiat, Phys. Rev. B 72, 033414 (2005).
- (31) H. B. Heersche, Z. de Groot, J. A. Folk, H. S. J. van der Zant, C. Romeike, M. R. Wegewijs, L. Zobbi, D. Barreca, E. Tondello, and A. Cornia, Phys. Rev. Lett. 96, 206801 (2006).
- (32) M.-H. Jo, J. E. Grose, K. Baheti, M. M. Deshmukh, J. J. Sokol, E. M. Rumberger, D. N. Hendrickson, J. R. Long, H. Park, and D. C. Ralph, Nano Lett. 6, 2014 (2006).
- (33) L. Bogani and W. Wernsdorfer, Nat. Mater. 7, 179 (2008).
- (34) S. Kahle, Z. Deng, N. Malinowski, C. Tonnoir, A. Forment-Aliaga, N. Thontasen, G. Rinke, D. Le, V. Turkowski, T. S. Rahman, S. Rauschenbach, M. Ternes, and K. Kern, Nano Lett. 12, 518 (2012).
- (35) M. Octavio, M. Tinkham, G.E. Blonder, and T.M. Klapwijk, Phys. Rev. B 27, 6739 (1983).
- (36) E.N. Bratus, V.S. Shumeiko, and G. Wendin, Phys. Rev. Lett. 74, 2110 (1995).
- (37) D. Averin and A. Bardas, Phys. Rev. Lett. 75, 1831 (1995).
- (38) J.C. Cuevas, A. Martín-Rodero, and A.Levy Yeyati, Phys. Rev. B 54, 7366 (1996).
- (39) J.C. Cuevas, J. Heurich, A. Martín-Rodero, A. Levy Yeyati, and G. Schön, Phys. Rev. Lett. 88, 157001 (2002).
- (40) Y. Uzawa and Z. Wang, Phys. Rev. Lett. 95, 017002 (2005).
- (41) M. Chauvin, P. vom Stein, H. Pothier, P. Joyez, M.E. Huber, D. Esteve, and C. Urbina, Phys. Rev. Lett. 97, 067006 (2006).
- (42) L. Landau and E. M. Lifshitz, Phys. Z. Sowjetunion 8, 153 (1935); T. L. Gilbert, Phys. Rev. 100, 1243 (1955).
- (43) G. Eilenberger, Z. Phys. 214,195 (1968); A.I. Larkin and Y.N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 55, 2262 (1968) [Sov. Phys. JETP 28, 1200 (1969)]; G.M. Eliashberg, Zh. Eksp. Teor. Fiz. 61, 1254 (1971) [Sov. Phys. JETP 34, 668 (1972)].
- (44) J.W. Serene and D. Rainer, Phys. Rep. 101, 221 (1983).
- (45) L. Cai and E.M. Chudnovsky, Phys. Rev. B 82, 104429 (2010).