arxiv:cond-mat/ v1 [cond-mat.other] 13 Apr 2006

Transcription

1 arxiv:cond-mat/060439v1 cond-mat.other] 13 Apr 006 Spin-wave instability for parallel pumping in ferromagnetic thin films under oblique field Kazue Kudo, Katsuhiro Nakamura Department of Applied Physics, Osaka City University, Osaka , Japan Abstract Spin-wave instability for parallel pumping is studied theoretically. The spin-wave instability threshold is calculated in ferromagnetic thin films under oblique field which has an oblique angle to the film plane. The butterfly curve of the threshold usually has a cusp at a certain value of the static external field. While the static field value of the cusp point varies as the oblique angle changes, the general properties of the butterfly curve show little noteworthy change. For very thin films, however, multiple cusps of the butterfly curve can appear due to different standing spin-wave modes, which indicate a novel feature for thin films under oblique field. Key words: spin-wave instability, parallel pumping, thin film, butterfly curve PACS: g, Ds 1 Introduction The first experiment for the parametric excitation of spin wave by parallel pumping was given by Schlöman et al. 1]. In a parallel pumping experiment, a microwave magnetic field is applied parallel to the external static field. When the microwave field amplitude h exceeds a certain spin-wave instability amplitude h crit, the parametric spin-wave excitation occurs. The excited spin waves have half the pumping frequency. The h crit curve plotted against the static field is called butterfly curve. The spin-wave instability for parallel pumping is often investigated for the purpose to study relaxation phenomena in ferromagnetic materials, since the instability threshold depends on the balance between the driving power and the damping of spin waves. Corresponding author. Tel.: ; fax: address: kudo@a-phys.eng.osaka-cu.ac.jp (Kazue Kudo). Preprint submitted to Elsevier Science 5 December 013

2 Theoretical studies have not been very successful to explain the instability threshold for parallel pumping: a combination of more than two equations (theories) is needed to obtain the butterfly curve. In other words, even a qualitative explanation for the butterfly curve cannot be successful so long as one tries to use only the Landau-Lifshitz (LL) equation governing the magnetization dynamics. By contrast the Suhl instability for perpendicular pumping has been well explained theoretically ]. The Suhl instability can be described by using only the LL equation. For parallel pumping, however, the LL equation is not enough to explain the instability, and the spin-wave line width H k is necessary to describe the relaxation of spin-waves as well as the instability. A typical butterfly curve for parallel pumping has a cusp at a certain external static field, which can be explained if the spin-wave line width H k is given. For example, Patton et al. proposed some trial H k functions and successfully explained butterfly curves 3]. Such a typical feature of butterfly curves usually does not depend on the shape of samples. For very thin films, however, more interesting features can appear. In fact, butterfly curves with multiple cusps were observed for 0.5-µm yttrium iron garnet (YIG) films under in-plane external field 4]. Those cusps were due to quantized standing-wave modes across the film. It is highly desirable to develop the theory of multi-cusp butterfly curves in the case of a general oblique field. In this paper, we examine the instability threshold for parallel pumping in ferromagnetic thin films theoretically. The butterfly curve of the threshold will be calculated in cases where the external field has an oblique angle to the film plane. For the calculation of the threshold, we use a trial H k function proposed by Patton et. al. 3], which was proposed originally for spherical samples but can be applied to thin films 5,6]. For usual thin films, the butterfly curve of the threshold has a single cusp at a certain value of the static external field. We will show how the cusp shifts when the oblique angle changes. Moreover, we will reveal that multi-cusp butterfly curves can be seen for very thin films. Equations of motion for the spin-wave amplitudes The dynamics of magnetization field M(r) is governed by the Landau-Lifshitz (LL) equation, 1 γ tm = M H eff λ M 0 M (M H eff ). (1) Here, γ is the gyromagnetic ratio (γ < 0 for spins); M 0 is the value of the magnetization in thermal equilibrium; λ is a damping constant parameter

3 (λ < 0 in this case); H eff is the effective magnetic field: H eff = D M + H d + H a + H 0 + h cosωt. () The first term comes from the exchange interaction; the third term is the anisotropy field, which is omitted below for convenience; the fourth and fifth terms are the external static and pumping fields, respectively. For parallel pumping, H 0 //h and they are parallel to the z axis. The second term on right hand side of Eq. () is the demagnetizing filed given by the gradient of the magneto-static potential φ: H d = φ. (3) The magneto-static potential obeys the Poisson equation: 4π M, inside the sample, φ = 0, outside the sample. (4) The number of independent components of M is two, since the length of the magnetization vector is invariant ( M = M 0 ). The normalized magnetization S = M/M 0 can be represented by a point on a unit sphere. It is convenient to project the unit sphere stereographically onto a complex variable ψ(r, t) 7]: ψ = S x + is y 1 + S z, (5) where S x = ψ + ψ 1 + ψψ, S y = i(ψ ψ) 1 + ψψ, S z = 1 ψψ 1 + ψψ. (6) In terms of ψ, the LL equation (1) is rewritten as t ψ i(1 iλ) γ {DM 0 ψ ψ ( ψ) ] ψψ (1 ψ ) x φ i } (1 + ψ ) y φ + ( z φ H 0 h cosωt)ψ = 0. (7) Inside the sample, φ satisfies φ= 4πM 0 (1 + ψψ ) { (1 ψ ) x ψ + (1 ψ ) x ψ 3

4 d ~ z α z H0 // h ~ x x Fig. 1. Schematic picture of sample setting. The sample extends on the x-ỹ plane. The z axis corresponds to the direction of the external filed and has an angle α to the z axis. + i(1 + ψ ) y ψ (1 + ψ ) y ψ] (ψ z ψ + ψ z ψ) }. (8) Now let us consider the boundary conditions, which affect the demagnetizing field. We assume a film of thickness d infinitely extended in the x-ỹ plane as shown in Fig. 1. The external field is at an angle α to the z axis and corresponds to the z axis. We assume unpinned surface spins, which satisfy Neumann-like boundary conditions, z S z=±d/ = 0, (9) namely, z ψ = z z=±d/ ψ = 0. (10) z=±d/ Before proceeding to the linear instability of spin waves, we introduce the following dimensionless time and space units 8]: t t 4π γ M 0, r rd. (11) Then we obtain the linearized equations of motion corresponding to Eqs. (7) and (8): t ψ + i(1 iλ) l ψ 1 ] ( x + i y )Φ + ( z Φ ω H ω h cos ω p t)ψ = 0,(1) 4

5 Φ = ( x i y )ψ + ( x + i y )ψ when 1 < z < 1, (13) where l = D 4πd, Φ = φ 4πM 0 d, ω H = H 0 4πM 0, ω h = h 4πM 0, ω p = ω 4πM 0 γ.(14) Now we consider the undriven case (i.e., ω h = 0), and expand ψ(r, t) and ψ (r, t) so that they fulfills the boundary conditions (10). For even modes, k z = mπ (m: integer), ψ( r, t) = k ψ ( r, t) = k a k(t) e i( k x x+ k yỹ) cos k z z a k (t) ei( k x x+ k yỹ) cos k z z. (15) For odd modes, k z = (m + 1)π (m: integer), ψ( r, t) = k a k(t) e i( k x x+ k yỹ) sin k z z ψ ( r, t) = k a k (t) ei( k x x+ k yỹ) sin k z z. (16) Here, k x = k x cosα + k z sin α, k y = k y, and k z = k z cosα k x sin α. Using the expansions (15) and (16), we obtain solutions of Eq. (13): for even modes, Φ( r) = k e i( k x x+ k yỹ) k i( k a k + k + a k ) { cos k z z ( 1) m e k / cosh k z } for odd modes, + k z sin α (a k + a k ) { sin k z z ( 1) m e k / ( k z / k ) sinh k z }], (17) Φ( r) = k e i( k x x+ k yỹ) k i( k a k k + a k ) { sin k z z ( 1) m e k / sinh k z } k z sin α (a k a k ) { cos k z z + ( 1) m e k / ( k z / k ) cosh k z }],(18) where k + = k x cos α+i k y, k = k x cos α i k y and k = k x + k y. The detailed derivation of these solutions is shown in Appendix A. Here we approximate z Φ = H d z /4πM 0 with the value of uniform magnetization, k = 0: H d z = H d x sin α + H d z cosα = 4π(N x M x sin α + N z M z cosα), where N x and N z are 5

6 demagnetizing factors. In this case, N x = Nỹ = 0 and N z = 1. Assuming M z M 0, namely, M z M 0 cosα, we obtain z Φ = H d z /4πM 0 = cos α. Then Eq. (1) for ω h = 0 is rewritten as t ψ + i(1 iλ) l ψ 1 ] (cosα x + i ỹ sin α z )Φ + (cos α ω H )ψ = 0.(19) Combining Eqs. (15)-(19), we obtain t a k + i(1 iλ)a k a k + i(1 iλ)b k a k = 0, (0) where A k = C k + Re C k, B k = ( 1) n ( D k e iϕ k + C k ). (1) Here, n is an integer, cos ϕ k = k x /k, sin ϕ k = k y /k, k = k x + k y, and C k = cos α l k ω H 1 f k sin θ k, C k = k z k sin α ( k e iϕ k + k z sin α ) f k, D k = 1 f k sin θ k, ( f k = 1 ( δ 0n ) 1 exp k k ]) k k z z, k sin θ k = k k, kz = πn. () The detailed derivation is given in Appendix B. Equation (0) represents two coupled harmonic oscillators a k, a k. We diagonalize Eq. (0) by means of the Holstein-Primakoff transformation: where a k =ν k b k µ k b k a k =ν kb k µ k b k, (3) ν k =cosh χ k, µ k =e iβ k sinh χ k, 6

7 C k + Re C k cosh χ k = (C k + Re C k ) (1 + λ ) D k e iϕ k ]. (4) 1/ + Ck Substituting Eq. (3) into Eq. (0), we obtain t b k i(ω k + iη k )b k = 0, (5) where ωk = (C k + ReC k ) (1 + λ ) D k e iϕ k + C k, (6) η k = λ(c k + Re C k ) (7) Equations (6) and (7) express the dispersion relation and a damping rate, respectively. 3 The instability threshold Now we consider the case where the microwave field h cosωt is applied. The equation of motion of a k corresponding to Eq. (0) is t a k + i(1 iλ)a k a k + i(1 iλ)b k a k i(1 iλ)ω h cosω p t a k = 0. (8) After the the Holstein-Primakoff transformation (3), we substitute the following equations into Eq. (8), b k (t) =b o k (t) expi(ω p/)t η k t], b k (t) =bo k (t) exp i(ω p/)t η k t], (9) since the resonance occurs at ω k = ω p /. Then the slowly-varying variable b o k satisfies ( t b o k + ω k ω ) ] p ρk b o k = 0, (30) where D k e iϕ k + Ck ρ k = ω h 1 + λ. (31) ω k 7

8 Therefore, the exponentially increasing solution for b k is b k exp ( ρ k η k )t + i(ω p /)t], (3) where ρ k > η k. (33) The instability threshold ωh crit ω crit h = ω p 1 + λ min k is now given as η k D k e iϕ k +. (34) Ck The threshold, Eq. (34), is obtained by using only the LL equation, but cannot explain the experimentally-observed instability threshold for parallel pumping. In fact, Eq. (34) proves to give an instability curve totally different from the real butterfly curve. Since the relaxation phenomenon is essentially nonlinear one, the linearizion analysis is not sufficient to discuss the instability threshold. The problem is beyond the purpose of this paper, and it will be discussed in a forthcoming paper 9]. This difficulty, however, can be overcome by replacing the damping rate η k by a suitable spin-wave line width. Using the spin-wave line width H k, we rewrite Eq. (34) as ω crit h = ω p min k H k D k e iϕ k. (35) Ck Here we adopt a simple trial H k function 3], H k = A 0 + A 1 sin (θ k ) + A k, (36) where A 0, A 1 and A are adjustable parameters. 4 Butterfly curves Typical butterfly curves of the threshold ωh crit have a cusp at a certain static field: as the static field increasing, ωh crit decreases below the cusp point and increases above that. For static fields below the cusp point, the minimum threshold modes corresponds to a spin wave propagating with θ k = π/. As the static field increases, the wave vector k of the threshold modes decreases, 8

9 hc Oe] hc Oe] α= H0 Oe] α=π/ H0 Oe] hc Oe] hc Oe] α=π/ H0 Oe] α=π/ H0 Oe] Fig.. Theoretical butterfly curves for d = 15µm for some values of the oblique angle α. and k 0 at the cusp. For static fields above the cusp point, the wave number remains at k 0 and θ k decreases from π/ to 0. Figure shows theoretical butterfly curves calculated by using Eqs. (35) and (36). The material parameters used in the calculation are typical values for yttrium iron garnet (YIG) materials: γ = rad/(s Oe), 4πM 0 = Oe, D = Oe cm /rad. The other parameters are as follows: ω/(π) = Hz, d = 15µm, A 0 = /(4πM 0 ), A 1 /A 0 = 0.4, A /A 0 = /d. We calculate the static field value for the cusp point with k = 0 and θ k in a bulk approximation: we assume f k = 1 in Eq. (). Setting ω k = ω p / in Eq. (6), we obtain the static field for the cusp point ω crit H = cos α ω p + 1. (37) While the cusp point shifts to lower static fields as the oblique angle α increases, the typical shape of curves is found to show little noteworthy change. In some cases, however, curious butterfly curves can emerge. Figure 3 is a theoretical butterfly curve for the case of d = 5µm and α = π/. The other parameters used in the calculation are the same as used in Fig.. There appear 9

10 .5 hc Oe] H0 Oe] Fig. 3. Theoretical butterfly curve for d = 5µm. The external field is applied parallel to the film plane. multiple cusps in the low field region of the butterfly curve. These cusps are due to different standing spin-wave modes across the film thickness. Let us recall that the z component of wave vectors are quantized. Effects of quantization are essential when the boundary conditions become important. Comparing the results in Figs. and 3, one may say multi-cusp feature of a butterfly curve can be seen when the thickness of a film is very small. This finding is supported by some experimental studies 4,5]. In Ref. 5], a standard butterfly curve of the instability was shown for 15.9-µm-thick yttrium iron garnet (YIG) film under in-plane external field. However, in Ref. 4], multiple butterfly curves appeared for a YIG film under in-plane external field. The difference between the two is just the thickness. The thickness of the film in Ref. 4] was 0.5µm. 5 Discussion We have revealed how butterfly curves depend on the oblique angle between the external field and the film plane. We have calculated theoretical butterfly curves by using the LL equation together with the spin-wave line width H k. The calculation was performed under the assumption that parameters A 0, A 1 and A for H k do not depend on the oblique angle. These parameters were originally introduced to fit a theoretical butterfly curve to experimental data 3], and might depend on the oblique angle. We have also shown qualitative features of the novel aspect of butterfly curves with multiple cusps for parallel pumping. Those multi-cusp curves come from quantized standing-wave modes across the film. However, it is not clear how such cusps appear, since there remains an ambiguity about how to evaluate A 0, A 1 and A. Figure 3 is one of possible theoretical curves. To proceed to quantitative evaluations, the experimental test is highly desirable to confirm the features predicted here. 10

11 Acknowledgments The authors thank to Prof. Mino of Okayama university for useful discussion. One of the authors (K. K.) is supported by JSPS Research Fellowships for Young Scientists. A Demagnetization field First, we should use the ( x, ỹ, z) coordinate system rather than (x, y, z) to use Eqs. (15) and (16) in the Poisson equation (4), and employ the transformation:, x cosα x sin α z. We consider only the case of 1 < z < 1. For even modes, k z = mπ (m: integer), the Poisson equation is rewritten as Φ = k e i( k x x+ k yỹ) (i k cos k z z + k z sin α sin k z z)a k + (i k + cos k z z + k z sin α sin k z z)a k], (A.1) where k + = k x cosα + i k y and k = k x cosα i k y. For odd modes, k z = (m + 1)π (m: integer), Φ = k e i( k x x+ k yỹ) (i k sin k z z k z sin α cos k z z)a k (i k + sin k z z k z sin α cos k z z)a k]. (A.) Here we note a well known fact: the Poisson equation, φ(r) = 4πρ(r), (A.3) has the solution as ρ(r ) dv φ(r) = r r. (A.4) Applying Eq. (A.4), we obtain the solution of Eqs. (A.1) and (A.): for even modes, Φ( r) = 1 4π k i( k a k + k + a ) k dṽ r r ei( k x x + k yỹ ) cos k z z 11

12 + k z sin α (a k + a k ) dṽ ei( k x x + k yỹ ) r sin r k z z ;(A.5) for odd modes, Φ( r) = 1 4π k i( k a k k + a ) k k z sin α (a k a k ) After integration, Eqs. (A.5) and (A.6) become dṽ ei( k x x + k yỹ ) r sin r k z z dṽ ei( k x x + k yỹ ) r cos r k z z.(a.6) and Φ( r) = k e i( k x x+ k yỹ) { i( k k a k + k + a ) cos k k z z ( 1) m e k / cosh k ] z + k z sin α (a k + a ) sin k k z z ( 1) m e k ]} k / z sinh k k z, (A.7) Φ( r) = k e i( k x x+ k yỹ) { i( k k a k k + a ) sin k k z z ( 1) m e k / sinh k ] z k z sin α (a k a ) cos k k z z + ( 1) m e k ]} k / z cosh k k z,(a.8) respectively. B Equation of motion for a k First, let us calculate the derivative of Φ in Eq. (19). For even modes, k z = mπ (m: integer), (cos α x + i ỹ sin α z )Φ = k e i( k x x+ k yỹ) k { ( k + k a k + k + a k ) cos k z z ( 1) m e k / cosh k ] z + i k + kz sin α (a k + a k ) sin k z z ( 1) m e k / ( k z / k ) sinh k z ]} 1

13 + sin α { i( k a k + k + a k ) k z sin k z z ( 1) m e k / k sinh k z ] + k z sin α (a k + a k ) cos k z z ( 1) m e k / cosh k z ]}], (B.1) where k + = k x cosα + i k y and k = k x cosα i k y. For odd modes, k z = (m + 1)π (m: integer), (cos α x + i ỹ sin α z )Φ = k e i( k x x+ k yỹ) k { ( k + k a k k + a k ) sin k z z ( 1) m e k / sinh k ] z i k + kz sin α (a k + a k ) cos k z z ( 1) m e k / ( k z / k ) cosh k z ]} + sin α { i( k a k k + a k ) kz cos k z z ( 1) m e k / k cosh k z ] + k z sin α (a k a k ) sin k z z ( 1) m e k / sinh k z ]}]. (B.) Let us project Eqs. (B.1) and (B.) onto cos k z z and sin k z z, respectively. For 1/ < z < 1/, the projection F( z) of a function f( z) onto cos k z z is given by, F( z) = 1/ d zf( z) cos k z z / 1/ d z cos kz z. (B.3) 1/ 1/ After the projection, Eqs. (B.1) and (B.) are reduced to (cosα x +i ỹ sin α z )Φ = k e i( k { x x+ k yỹ) k 1 ( δ m0 ) ( 1 e k ) k k } { ( k + k + k z sin α)a k + ( k + + k z sin α)a k} cos kz z, (B.4) where δ ij is the Kronecker delta, and (cosα x +i ỹ sin α z )Φ = k e i( k { x x+ k yỹ) k 1 ( 1 e k ) k k } { ( k + k + k z sin α)a k ( k + + k z sin α)a k} sin kz z, (B.5) respectively. Let us impose a restriction on the kinds of parameters about wave vectors to describe the equation of motion: we only use k, k = k x + k y, k z = πn (n: integer), θ k and ϕ k. Here, sin θ k = k /k; cosϕ k = k x /k ; sin ϕ k = k y /k. The restriction leads to the following: 13

14 k + k + k z sin α = k + (k cosϕ k + k z sin α) k z sin α; k + + k z sin α = k eiϕ k +(k e iϕ k + k z sin α) k z sin α; k = k k z. (B.6) From the above two equations and Eq. (19), we obtain Eq. (0): t a k + i(1 iλ)a k a k + i(1 iλ)b k a k = 0. (B.7) For even modes, A k = C k + Re C k, B k = ( D k e iϕ k + C k ), (B.8) and for odd modes, A k = C k + Re C k, B k = ( D k e iϕ k + C k ). (B.9) Here we note that C k, Ck and D k are described in Eq. (). The difference between Eqs. (B.8) and (B.9) is just the sign of B k. Therefore, for both even and odd modes, we may write B k = ( 1) n ( D k e iϕ k + C k ). (B.10) References 1] E. Schlömann, J. J. Green and U. Milano, J. Appl. Phys. 31, 386S (1960). ] H. Suhl, J. Phys. Chem. Solids 1, 09 (1957). 3] See, for example, M. Chen and C. E. Patton, in Nonlinear Phenomena and Chaos in Magnetic Materials, edited by P. E. Wigen (World Scientific, Singapore, 1994), pp ] B. A. Kalinikos, N. G. Kovshikov and N. V. Kozhus, Sov. Phys. Solid State 7, 1681 (1986). Fiz. Tverd. Tela (Leningrad) 7, 794 (1985).] 5] G. Wiese, L. Buxmzn, P. Kabos and C. E. Patton, J.Appl. Phys. 75, 1041 (1994). 6] P. Kabos, M Mendik, G. Wiese and C. E. Patton, Phys. Rev. B 55, (1997). 7] M. Lakshmanan and K. Nakamura, Phys. Rev. Lett. 53, 497 (1984). 8] F. J. Elmer, Phys. Rev. B 53, 1433 (1996). 9] K. Kudo and K. Nakamura, in preparation. 14