# Zero-Field Fiske Resonance Coupled with Spin-waves in Ferromagnetic Josephson Junctions

AC Josephson current density in a Josephson junction with DC bias is spatially modulated by an external magnetic field, and induces an electromagnetic (EM) field inside the junction. The current-voltage (-) curve exhibits peaks due to the resonance between the EM field and the spatially modulated AC Josephson current density. This is called Fiske resonance. Such a spatially modulated Josephson current density can be also induced by a non-uniform insulating barrier and the Fiske resonance appears without external magnetic field. This is called zero-field Fiske resonance (ZFFR). In this paper, we theoretically study the ZFFR coupled with spin-waves in a superconductor/ferromagnetic insulator/superconductor junction (ferromagnetic Josephson junction) with a non-uniform ferromagnetic insulating barrier. The resonant mode coupled with spin-waves can be induced without external magnetic field. We find that the - curve shows resonant peaks associated with composite excitations of spin-waves and the EM field in the junction. The voltage at the resonance is obtained as a function of the normal modes of EM field. The ZFFRs coupled with spin-waves are found as peak structures in the DC Josephson current density as a function of bias voltage.

## 1 Introduction

The DC Josephson effect is characterized by the DC current flowing without a voltage-drop between two superconductors separated by a thin insulating barrier [1]. When a DC voltage is applied to the junction, the AC Josephson current with frequency flows in the junction driven by the difference of phases in two superconducting order parameters, i.e., Josephson-phase . If both the DC voltage and a magnetic field are applied to the junction, whose width is smaller than the Josephson penetration depth , the AC Josephson current density is spatially modulated and generates the electromagnetic (EM) field inside the junction. In this case, the current-voltage (-) curve exhibits peaks due to the resonance between the AC Josephson current density and the EM field. This is called Fiske resonance [2, 3, 4, 5, 6].

Josephson junctions composed of ferromagnetic metal (FM) and superconductors (Ss) are extensively studied for the last decade. The S/FM/S junctions exhibit fascinating phenomena which are not observed in the conventional Josephson junctions [7, 8, 9, 10, 11]. The interaction between Cooper pairs and spin waves in the FM is of importance in the transport properties in the S/FM/S and S/I/FM/S junctions. The dynamics of coupled with spin-waves in the FM has been investigated theoretically [12, 13, 14, 15, 16, 17, 18, 19, 20] and experimentally [21]. However, the Fiske resonance coupled with spin-waves is not yet observed experimentally.

Another type of Josephson junction with ferromagnetic insulator (FI) instead of the FM is also examined. It is reported that the dissipation effect in the S/FI/S junction is smaller than that the S/FM/S junction [22, 23, 24]. The damping of spin-waves is also very small in the FI compared to the case of the FM [25, 26, 27]. Therefore, the coupling between and spin-waves can be observed more clearly in the S/FI/S junction. In fact, in the S/FI/S junction, it is expected that the Fiske resonance has clear multiple structures associated with spin-wave excitation[28].

Here, we note that the Fiske resonance in the conventional Josephson junction is also induced by the non-uniform insulating barrier in the junction, since AC Josephson current density driven by a DC voltage is spatially modulated and then the EM field is generated inside the junction. In this case, the Fiske resonance occurs without external magnetic field. It is called zero-field Fiske resonance (ZFFR), which originates from the resonance between the EM field and the spatially modulated AC Josephson current density due to the non-uniform insulating barrier. This phenomenon has been widely studied experimentally and theoretically in the Josephson junction [6, 29, 30, 31, 32].

In this paper, we theoretically study the ZFFR coupled with spin-waves in an S/FI/S junction with a non-uniform FI. The merit of such a non-uniform geometry of junction is that the spatially modulated AC Josephson current density can be induced with no external magnetic field and thus the Fiske resonance occurs without external magnetic field. By solving the equation of motion of coupled with spin-waves, it will be found that the - curve shows resonant peaks. The voltage at the resonances is obtained as a function of the normal modes of EM field, which indicates composite excitations of the EM field and spin-waves in the S/FI/S junction. Dependence of those resonances on distributions of the Josephson critical current density is presented.

The rest of this paper is organized as follows. In Sec. II, we formulate the Josephson current in a Josephson junction with a non-uniform ferromagnetic insulator. In Sec. III, the DC component of Josephson current density at the ZFFR with spin-waves is calculated analytically and numerically. Summary is given in Sec. IV.

## 2 FORMULATION of JOSEPHSON CURRENT in FERROMAGNETIC JOSEPHSON JUNCTION with MAGNETIC INSULATOR

The system considered is a Josephson junction with a FI sandwiched by two superconductors with -wave symmetry as shown in Fig.1. The geometry of the junction is assumed to be a non-uniform junction to impose non-uniform Josephson current density without external magnetic field [33]. The magnetization in the FI is parallel to the -direction. Here, we assume a simple model of non-uniform Josephson current density given by,

(1) | |||||

(2) | |||||

(3) |

where and are the Josephson critical current density and Josephson frequency with bias voltage , respectively [34]. is the Josephson critical current density for the uniform geometry of junction, i.e., . The electromagnetic dynamics induces depending on space and time. The distribution of is determined by two parameters and , where is imposed. In the S/FI/S junction, spin-waves can be excited by the EM field inside the FI due to the AC Josephson current. In this situation, the equation of motion for coupled with spin waves is described by [28],

(4) | |||||

(5) |

The effective velocity of light in the FI is given by , Josephson penetration depth , dielectric constant and permeability . The London penetration depth is denoted by and means the damping factor caused by quasi-particle resistivity in the FI. The magnetic susceptibility of the FI in the linearized Landau-Lifshitz-Gilbert equation is given by [25],

(6) |

where , , and are the -component of the magnetization, Gilbert damping factor, and the gyromagnetic ratio, respectively. Magnetic susceptibility and spin-wave energy in a magnetic material are generally modified by geometry and thickness. On the other hand, Eq. (6) is obtained by assuming a uniform FI. This is justified, because the magnetic susceptibility and are insensitive to the thickness of FI, provided that the conformation of the ferromagnetic materials changes on a scale of nanometers [35]. Therefore, we adopt Eq. (6) and obtained in the uniform FI [36] for the non-uniform FI as an approximation, since we consider the thickness change of FI to be in a range of a few nanometers. In the FI, the dispersion relation of spin-waves with the frequency is given by

(7) |

where . The anisotropic field and the stiffness of spin-waves in the FI are denoted by and , respectively. The spin-wave having a finite wave number is neglected in the Fiske resonance because of the following reason: In Eq.(7), the first term is caused by the anisotropic and demagnetizing fields, and the wave number is given by . In a conventional FI, is about tens of eV [25]. On the other hand, is of the order peV due to the small stiffness of spin-waves [37] when is a few mm. Below, we only consider mode for spin-waves with the constant frequency .

## 3 DC Josephson current density with ZFFR and numerical results

In order to obtain the solution of Eq (4), we expand in terms of the normal modes of the EM field generated by the AC Josephson current as follows,

(8) |

where is a complex number and . This equation of satisfies , which is Kulik’s boundary condition [5, 6]. We consider to be a small perturbation and solve Eq.(8) by taking to be . Substituting Eq.(8) into Eq.(4), is determined as,

(9) | |||||

where .

Next, we calculate the DC Josephson current density coupled with spin waves as a function of . The function, , is expanded in terms of and is given by

(10) |

Introducing Eqs. (8) and (9) into Eq. (10), the analytic formula of without external magnetic field is obtained as,

(11) | |||||

(12) |

where , . Equation (11) clearly demonstrates that zero-field resonant modes depend on parameters and which determine the distribution of the Josephson critical current density flowing through the FI. Hence, one can easily find that three cases are possible for the zero-field resonance. When (), the zero-field resonance only appears at even (odd) numbers of . On the other hand, when , the zero-field resonance appears at all integers .

Next, we derive a condition for the ZFFR in the present system by analyzing Eq. (12). When the denominator of is minimum with respect to , takes a maximum, so that the DC Josephson current exhibits the resonant behavior. The DC voltage, at which the resonance occurs, is determined by neglecting the damping term of Eq. (12) as . Setting the denominator of to be zero, the voltage is given by

(13) |

We have two DC voltages, and , at which the ZFFR occurs for each . The integer is determined by the mode of the EM field in the junction. Eq. (13) clearly shows that there are two dispersion relations, which result from the coupling between the EM field and spin-waves in the FI. Note that the amplitude of ZFFR strongly depends on and as we will see in the next section.

At last, we numerically evaluate Eq. (11). Parameters are set to be T, [27], , [38], m/As[39], and GHz. Instead of plotting an - curve, at the resonances will be shown as a function of the voltage below. The amplitude of is associated with a height of resonant peak or a jump in the - curve (for instance, see Ref[6]). In the following numerical calculations, we exclude the contribution of in Eq.(11), since we discuss about the resonance between the spatially modulated AC Josephson current and standing wave of EM field.

Figure 2 shows induced by ZFFR as a function of [40] for and . With these parameters, is linearly distributed in the junction (see the inset of Fig. 2). The black (solid) and red (dashed) lines are and the component with in Eq. (11), respectively. This result clearly demonstrates that the Fiske resonance occurs without external magnetic field, i.e., ZFFR. The additional resonance peak around arises from the presence of spin-wave excitation in the FI. This resonance comes from the inhomogeneity of Josephson critical current density induced by the non-uniform geometry of junction. Moreover, in the present case, Eq. (11) becomes

(14) |

In Eq. (14), it is found that the ZFFR only occurs at odd number of . Since resonant peaks of ZFFR with are much smaller than that with , main contribution to the ZFFR as depicted in Fig. 2 is the mode of .

Figure 3 is the case for and . The black (solid) line is . Red (solid), blue (dashed), and green (chain) lines are each component with in Eq. (11). It is found that ZFFR peaks of clearly appear at in Fig. 3 in contrast to Fig. 2. The reason is simply due to the non-linear Josephson critical current density to contain both symmetric and antisymmetric components with respect to . The present non-linear distribution of the Josephson critical current density will be more realistic. Therefore, we can expect that multiple resonant peaks such as Fig. 3 is practically observed without external magnetic field.

## 4 Summary and Discussion

We have theoretically studied the zero-field Fiske resonance (ZFFR) in the S/FI/S junction with several patterns of spatial variation in the Josephson critical current density, which is induced by a non-uniform ferromagnetic insulating barrier. Such a non-uniform AC Josephson current density can excite the EM field inside the FI without external magnetic field. It is found that the current-voltage (-) curve shows two resonant peaks without external magnetic field in the present system, i.e., the ZFFR coupled with spin-waves occurs. Voltage at the resonances is obtained as a function of the normal modes of EM field, which indicates composite excitations of the EM field and spin-waves in the S/FI/S junction.

The present study will provide a platform to study the dynamics of Josephson phase and the magnetic excitation. Furthermore, in the non-uniform S/FI/S junction, several applications such as spin-current emitter by utilizing spin-wave excitation in the FI [41] may be also possible in analogy with the emission of coherent THz radiation in the high-T cuprate [42, 43, 44], although Josephson junctions based on the high-T cuprate are usually laminated structures differently from the single Josephson junction discussed here. In fact, the inhomogeneity of the junction was one of essential factors to realize the emission without external magnetic field [42, 43, 44]. However, novel devices using the S/FI/S junction are beyond the scope of the present paper and will be studied elsewhere.

###### Acknowledgements.

This work is supported by Grant-in-Aid for Research Activity Start-up (No. 25887053) from the Japan Society for the Promotion of Science and Grant-in-Aid for Scientific Research from MEXT (Grant No.24540387, No.24360036, No.23340093, and No.25287094), Center for Computational Science and e-Systems of JAEA, and the inter-university cooperative research program of IMR, Tohoku University.## References

- [1] B. D. Josephson, Phys. Lett. (1962) 251.
- [2] M. D. Fiske, Rev. Mod. Phys. (1964) 221.
- [3] R. E. Eck, D. J. Scalapino, and B. N. Taylor, Phys. Rev. Lett. (1964) 15.
- [4] D. D. Coon and M.D Fiske, Phys. Rev. (1965) A744.
- [5] I. O. Kulik, JETP Lett. (1965) 84.
- [6] A. Barone and G. Patern, (Wiley, New York, 1982).
- [7] A. A. Golubov, M. Yu. Kupriyanov, and E. Ilfichev, Rev. Mod. Phys. (2004) 411.
- [8] A. I. Buzdin, Rev. Mod. Phys. (2005) 935.
- [9] V. V. Ryazanov, V. A. Oboznov, A. Yu. Rusanov, A. V. Veretennikov, A. A. Golubov, and J. Aarts, Phys. Rev. Lett. (2001) 2427.
- [10] T. Kontos, M. Aprili, J. Lesueur, F. Genet, B. Stephanidis, and R. Boursier, Phys. Rev. Lett. (2002) 137007.
- [11] J. W. A. Robinson, S. Piano, G. Burnell, C. Bell, and M. G. Blamire, Phys. Rev. Lett. (2006) 177003; J. W. A. Robinson, S. Piano, G. Burnell, C. Bell, and M. G. Blamire, Phys. Rev. B (2007) 094522.
- [12] Z. Nussinov, A. Shnirman, D. P. Arovas, A. V. Balatsky, and J. X. Zhu, Phys. Rev. B (2005) 214520.
- [13] S. Takahashi, S. Hikino, M. Mori, J. Martinek and S. Maekawa, Phys. Rev. Lett. (2007) 057003.
- [14] C. Bell, S. Milikisyants, M. Huber, and J. Aarts, Phys. Rev. Lett. (2008) 047002.
- [15] M. Houzet, Phys. Rev. Lett. (2008) 057009.
- [16] S. Hikino, M. Mori, S. Takahashi, and S. Maekawa, J. Phys. Soc. Jpn. (2008) 053707.
- [17] F. Konschelle and A. Buzdin, Phys. Rev. Lett. (2009) 017001.
- [18] T. Yokoyama and Y. Tserkovnyak, Phys. Rev. B (2009) 104416.
- [19] A. F. Volkov and K. B. Efetov, Phys. Rev. Lett. (2009) 037003.
- [20] S. Mai, E. Kandelaki, A. F. Volkov, and K. B. Efetov, Phys. Rev. B (2011) 144519.
- [21] I. Petkovi, M. Aprili, S. E. Barnes, F. Beuneu, and S. Maekawa, Phys. Rev. B (2009) 220502.
- [22] S. Kawabata, S. Kashiwaya, Y. Asano, Y. Tanaka, and A. Golubov: Phys. Rev. B (2006) 180502(R).
- [23] S. Kawabata, Y. Asano, Y. Tanaka, A. A. Golubov, and S. Kashiwaya Phys. Rev. Lett. (2010) 117002.
- [24] Both the S/FM/S and the S/FI/S junctions are called ferromagnetic Josephson junction, that means the latter one in this paper.
- [25] B. Hillebrands and K. Ounadjela, Spin dynamics in confined magnetic structures II (Springer-Verlag Berlin Heidelberg, New York, 2003)
- [26] D. D. Stancil and A. Prabhakar, Spin Waves Theory and Applications (Springer Science + Business Media, LLC 2009)
- [27] Y. Kajiwara, K. Harii, S. Takahashi, J. Ohe, K. Uchida, M. Mizuguchi, H. Umezawa, H. Kawai, K. Ando, K. Takanashi, S. Maekawa, and E. Saitoh Nature (2010) 262.
- [28] S. Hikino, M. Mori, S. Takahashi, and S. Maekawa, J. Phys. Soc. Jpn. (2011) 074707.
- [29] M. Russo and R. Vaglio, Phys. Rev. B (1978) 2171.
- [30] T. C. Wang, J. Appl. Phys. (1979) 2859.
- [31] C. Camerlingo, M. Russo, and R. Vaglio, J. Appl. Phys. (1982) 7609.
- [32] C. Nappi, E. Sarnelli, M. Adamo, and M. A. Navacerrada, Phys. Rev. B (2006) 144504.
- [33] The non-uniform geometry of the junction is the essential point, since such a geometry introduces the non-uniform current distribution. In this paper, the Josephson critical current density have a linear current distribution with in this geometry when and .
- [34] Note that we use the voltage bias model for simplicity, since we only focus on the DC component of Josephson current[6].
- [35] For instance, Y. Sun, Y. Song, H. Chang, M. Kabatek, M. Jantz, W. Schneider, M. Wu, H. Schultheiss, and A. Hoffmann, Appl. Phys. Lett. (2012) 152405.
- [36] The distribution of Josephson current density with (see Eq. (2)) can be obtained by changing the thickness of insulating barrier along direction, since the Josephson current density exponentially decreases with increasing the thickness of insulating barrier. For instance, see Ref [6].
- [37] M. Pajda, J. Kudrnovsk, I. Turek, V. Drchal, and P. Bruno, Phys. Rev. B (2001) 174402.
- [38] M. P. Lisitskiy and M. V. Fistul, Phys. Rev. B (2010) 184505. We adopted the value of damping in a SIS junction as reported in this paper.
- [39] S. Chikazumi, Physics of Magnetism (Oxford University Press, New York, 1997), p.559-560.
- [40] For mm and , the magnitude of DC voltage () is about tens of .
- [41] S. Maekawa, S. O. Valenzuela, E. Saitoh, and T. Kimura, Spin Current (Oxford University Press, Oxford, 2012).
- [42] L. Ozyuzer, A. E. Koshelev, C. Kurter, N. Gopalsami, Q. Li, M. Tachiki, K. Kadowaki,T. Yamamoto, H. Minami, H. Yamaguchi, T. Tachiki, K. E. Gray, W.-K. Kwok, U. Welp, Science 318 (2007) 1291.
- [43] K. Kadowaki, H. Yamaguchi, K. Kawamata, T. Yamamoto, H. Minami, I. Kakeya, U. Welp, L. Ozyuzer, A. Koshelev, C. Kurter, K.E. Gray, W.-K. Kwok, Physica C 468 (2008) 634.
- [44] A. E. Koshelev and L.N. Bulaevskii, Phys. Rev. B 77 (2008) 014530.