The effect of a neighboring metal layer on the high-frequency characteristics
of a thin magnetic stripe
Marina Vroubela兲and Behzad Rejaei
Delft Institute of Microelectronics and Submicron Technology (DIMES), Faculty of Electrical Engineering, Mathematics, and Computer Science, Delft University of Technology, Feldmannweg 17, 2628 CT
Delft, The Netherlands
共Received 17 August 2007; accepted 8 April 2008; published online 9 June 2008兲
The spin-wave spectrum of a ferromagnetic stripe placed above a metallic layer with finite conductivity is studied by using the magnetostatic Green’s function formalism. It is shown that the frequency and linewidth of the resonances are uniquely determined by complex, mode-dependent demagnetization factors. The formalism developed is used to analyze the resonance characteristics of the magnetic stripe as a function of its width and separation from the metallic layer. © 2008 American Institute of Physics.关DOI:10.1063/1.2937211兴
I. INTRODUCTION
The resonance spectrum of spin waves in thin, finite fer-romagnetic 共FM兲 patterns has been intensively studied over the last decade.1–7Such studies are particularly important in view of the potential application of FM elements in magnetic data storage and microwave devices.8–11In many cases, how-ever, magnetic elements do not stand alone, but are 共par-tially兲 surrounded by metallic conductors. One, for instance, can think of microwave transmission lines with FM cores9–11 or inductive techniques utilizing coplanar waveguides for probing magnetization dynamics in thin films.12,13 Neverthe-less, the effect of neighboring metallic layers on the spwave spectra of confined magnetic samples has not been in-vestigated, although the former’s influence on the dispersion relation of spin waves in infinite films is well known.14,15
This paper presents an analysis of the spin-wave spec-trum of a magnetic stripe placed above a metal layer. The analysis is carried out by using the thin-film approximation of the magnetostatic Green’s function formalism.4,16 It is shown that the presence of the metal ground increases the overall effective anisotropy of the magnetic stripe, shifting the spin-wave resonances to higher frequencies. Our results agree with experimental data obtained from thin Permalloy 共NiFe兲 stripes built on top of a dielectric/metal substrate. Furthermore, we show that the finite conductivity of the metal ground yields a mode-dependent extrinsic damping constant, in addition to the intrinsic damping constant of the film. Finally, by using the formalism developed, we study the dependence of the resonance frequency and damping of spin-wave modes on the stripe-ground plane separation and stripe width.
II. THEORY
Consider the structure shown in Fig. 1, consisting of a FM rectangular stripe placed above a metallic ground plane 共infinitely extended in the x-z plane兲 covered by an isolating 共e.g., dielectric兲 layer with the thickness d. The width 共w兲,
thickness 共t兲, and length 共L兲 of the magnetic stripe are de-fined along the x-, y-, and z-directions, respectively. For sim-plicity, we assume uniformity along the z-direction.
For a FM stripe magnetized along its length共z-direction兲 the ac-magnetization m = mxxˆ + myyˆ gives rise to an ac-demagnetization field hm= hm,xxˆ + hm,yyˆ inside the stripe, where
hm共r兲 = − ⵜ
冕
Sⵜ
⬘
G共r,r⬘
兲 · m共r⬘
兲dS. 共1兲 Here r =共x,y兲, ⵜ=共x,y兲, and S is the cross section of the magnetic strip in the x-y plane. In the above equation G is Green’s function for the magnetic potential in two dimen-sions in the presence of a conductive ground plane and is given by共see Appendix兲G共r,r
⬘
兲 = G0共r,r⬘
兲 + F共兲 + F共*兲, G0共r,r⬘
兲 = − 1 4ln关共x − x⬘
兲 2+共y − y⬘
兲2兴, F共兲 = − 1 4ln+ 2关H1共兲 − Y1共兲兴 − 1 22, =共x − x⬘
兲 + i关y + y⬘
+ 2共d + t/2兲兴, 共2兲 where =冑
i0, with the 共angular兲 frequency, 0 the vacuum permeability, andthe electrical conductivity of the ground layer. In the above equation, H1and Y1are the Struvea兲Electronic mail: m.vroubel@ewi.tudelft.nl. FIG. 1. Cross section of a FM stripe with a width w and thickness t, placedat a distance d above a metallic ground plane.
and Bessel functions of the first kind, respectively. Note that, in the limit of a perfectly conducting ground plane共→⬁兲, Eq.共2兲 is reduced to G→⬁共r,r
⬘
兲 = − 1 4ln关共x − x⬘
兲 2+共y − y⬘
兲2兴 − 1 4ln兵共x − x⬘
兲 2+关y + y⬘
+ 2共d + t/2兲兴2其, 共3兲 which is the free space field of a magnetic line charge at x⬘
, y⬘
, plus the field of its image 共with the same sign兲 at x⬘
, −y⬘
− 2共d+t/2兲. Therefore, a perfect ground actually in-creases the demagnetization field inside the stripe. As we shall see later, this conclusion remains valid for imperfect grounds whereis finite.The excitation spectrum of the magnetic stripe can be determined by solving Eq. 共1兲 together with the linear-response equation,17 m共r兲 =ញ · hm共r兲, ញ = M ⍀H2 −2
冋
⍀H i − i ⍀H册
, ⍀H=H+ i␣, M=␥M,H=␥Ha, 共4兲where ␥ is the gyromagnetic constant, M is the saturation magnetization, Hais the static共dc兲 anisotropy field along the z-direction, and␣is the Gilbert damping constant. For a thin FM strip共tⰆW兲, the variation of the magnetization along the film thickness共y-direction兲 can be neglected. Hence, one can formulate the problem in terms of the average magnetization vector,4
m˜共x兲 =1 t
冕
−t/2t/2
m共x,y兲dy. 共5兲
Upon combining Eqs.共1兲and共4兲, and performing averaging over the film thickness t, one arrives at the one-dimensional matrix integral equation,
Qញ · m˜共x兲 +
冕
−W/2 W/2 Gញ 共x,x⬘
兲 · m˜共x⬘
兲 = 0, 共6兲 where Qញ =ញ−1= 1 M冋
⍀H − i i ⍀H册
, 共7兲and Gញ is matrix Green’s function,
Gញ 共x,x
⬘
兲 =冋
␦共x − x⬘
兲 − g1共x,x⬘
兲 ig2共x,x⬘
兲ig2共x,x
⬘
兲 g1共x,x⬘
兲册
, 共8兲 where␦共x−x⬘
兲 is the Dirac delta function. The functions g1 and g2 are symmetric, i.e., g2共x,x⬘
兲=g2共x⬘
, x兲 and g1共x,x⬘
兲 = g1共x⬘
, x兲 so that Gញij共x,x⬘
兲=Gញji共x⬘
, x兲 for i, j=1,2. 共The ex-pressions for g1 and g2are quite complicated and, therefore, are given in the Appendix.兲Equation共6兲possesses nontrivial solutions only for cer-tain frequencies, i.e., the magnetostatic resonance frequen-cies of the stripe. However, before trying to solve Eq.共6兲let us first consider the eigenvalue problem,
冕
−W/2 W/2 Gញ 共x,x⬘
兲 · ⌿k共x⬘
兲 = k⌿k共x兲, ⌿k共x兲 =冋
k1共x兲 k 2共x兲册
. 共9兲Because of the special form of Gញ 关Eq.共8兲兴, it can be shown that if an eigenvaluekand eigenfunction⌿k共x兲 satisfy Eq.
共9兲, then so do the combination ¯k= 1 −k,⌿¯k共x兲 =
冋
k2共x兲 −k
1共x兲
册
. 共10兲Returning to Eq. 共6兲, we next try a solution of the type m˜共x兲=Ak⌿k共x兲+Bk⌿¯k共x兲, where Ak, Bk are constants. This yields the equation
冋
⍀H+kM − ii ⍀H+共1 − k兲M
册
·冋
AkBk
册
= 0, 共11兲
which allows a nonzero solution only for
2= k 2
=关⍀H+共1 − k兲M兴共⍀H+kM兲, 共12兲 i.e., the spin-wave resonance frequencies of the stripe. Note that kand 1 −kcan be interpreted as modependent de-magnetization factors in the x- and y-directions, respectively.16
The solutions of Eq. 共6兲 correspond to standing spin waves caused by the reflection of the waves by the lateral edges of the stripe. Figure 2 shows the distribution of the in-plane magnetization m˜x共x兲 for the first four modes 共k = 1 , 2 , 3 , 4兲 of a 50-m-wide and 0.1-m-thick Permalloy stripe, placed 1m above a conducting ground plane. Like in the case of a free magnetic stripe,1 the distribution of m˜x共x兲 for different modes roughly resembles standing sinu-soidal waves arising from imposing the quantization condi-FIG. 2. The共real part of the兲 in-plane magnetization m˜xfor the first four
eigenmodes 共k=1,2,3,4兲 of a 50-m-wide and 0.1-m-thick Permalloy stripe, placed 1m above a conducting ground plane. The conductivity of the metal ground is= 3.3⫻107S/m. The corresponding eigenvalues were
1= 0.003 96+ 0.000 23i, 2= 0.008 77+ 0.000 67i, 3= 0.012 78+ 0.001 18i,
and 4= 0.016 33+ 0.001 67i. The results were obtained by the numerical
tion kx= k/w on spin waves propagating with a wave num-ber kxalong the x-direction in an unbounded magnetic film. Although expression共12兲for the resonance frequency is identical to that of a free magnetic stripe, i.e., in the absence of a ground plane4,16 one should bear in mind that Gញ and, therefore, its eigenvalues k depend on the stripe-ground plane distance d, the conductivityof the ground, and even the frequency共through the parameter兲. The dependence of k on the frequency implies that Eq. 共12兲 should be viewed as a self-consistent equation for the resonance fre-quenciesk. Furthermore, because the components of Gញ are complex quantities, k are complex as well. Note that the imaginary part ofk=k
⬘
+ ik⬙
is, in fact, induced by the finite conductivity of the ground. A perfect ground conductor would yield a real matrix Gញ and, through symmetry proper-ties of Gញ, real values ofk.The imaginary part of k contributes to the imaginary part of the resonance frequencyk=k
⬘
+ ik⬙
. If bothk⬙
and the Gilbert damping constant ␣ 共representing the intrinsic magnetic relaxation loss of the stripe兲 are small, then it fol-lows from Eq.共12兲that共k
⬘
兲2⬇ 关H+共1 − k⬘
兲M兴共H+k⬘
M兲, 共13兲 k⬙
⬇冉
H+ M 2冊
␣+ 共1 − 2k⬘
兲M 2 2k⬘
k⬙
. 共14兲Note that in an actual experimentk
⬘
corresponds to the cen-tral frequency of the resonance observed while the resonance linewidth is given by⌬k⬇2k⬙
.16Thus, the Ohmic losses in the ground conductor cause additional absorption of electro-magnetic energy, increasing the magnetostatic resonance linewidth by an amount proportional tok⬙
. In analogy with the intrinsic magnetic loss, the extra loss induced by the ground plane can be expressed in terms of the mode-dependent, extrinsic damping constant␣k ex = 共1 − 2k
⬘
兲M 2 k⬘
共2H+M兲 k⬙
. 共15兲The overall damping constant for any given mode is then simply␣+␣kex.
To obtain an experimental verification of the model pre-sented, we measured the resonance frequency of the lowest mode 共k=1兲 Permalloy 共NiFe兲 stripes with a width of 100m and thicknesses of 100 and 200 nm. The stripes were built on top of a 2-m-thick aluminum layer with a conductivity of= 3.3⫻107 S/m, covered by a 1-m-thick SiO2 layer for isolation. A second SiO2 layer was next de-posited to cover the stripes, followed by the fabrication of microstrip probe lines on top of the structure. High-frequency impedance measurement of the microstrip lines was performed by using a vector network analyzer. Magnetic parameters共saturation magnetization, intrinsic magnetocrys-talline anisotropy兲 of the Permalloy film were found from M-H loop measurements performed by using a Princeton Mi-cromag 2900 AGM. Based on the values obtained 共M = 1.15 T , Ha= 25 Oe兲, we calculated the first spin-wave
reso-nance frequency of the stripes. The results, given in TableI, are in good agreement with the experimental values. III. NUMERICAL EXPERIMENTS
In this section, we use the formalism discussed in Sec. II to numerically study the dependence of the magnetostatic resonance spectra of FM stripes on their width and distance from a conducting ground plane. For our simulations we use a saturation magnetization of M = 1 T and an intrinsic 共mag-netocrystalline兲 anisotropy field of Ha= 5 Oe 共these are the most common values reported for Permalloy films in the literature兲. The electrical conductivity of the ground plane is taken to be= 3.3⫻107 S/m.
Figure3shows the resonance frequencyk
⬘
/2 and the extrinsic damping constant␣kex共k=1,2,3,4兲 of the first four spin-wave resonances calculated as function of the ground plane separation d, for a 50-m-wide and 0.1-m-thick FM stripe. The resonance frequency increases if the ground plane is brought closer to the stripe by reducing d. This is because, as mentioned in Sec. II, the ground plane strengthens the ac-demagnetization field induced by the magnetization dis-tribution on the stripe. This leads to larger values of the demagnetization coefficientkand, in turn, higher resonance frequencies. In fact, it can be shown than if the magnetic stripe is directly put on top of a perfect ground共d=0兲, then the resulting demagnetization factors are twice larger than that of a free magnetic stripe. Reducing the distance d also leads to an increase in the effective extrinsic damping con-stant ␣kexwhich is also shown in Fig. 3. This, of course, is due to the increase in the flow of eddy currents in the ground plane caused by the ac magnetic field around the stripe.
By increasing the distance d, the resonance frequencies saturate to their “free space” values, as one would have ex-pected. The saturation occurs earlier for modes with higher numbers. This result can be understood by noting that the magnetic field accompanying a spin-wave decays exponen-tially outside the magnetic film in the y-direction共Fig.1兲. A
ground plane at a distance larger than this decay length is not “felt” by the stripe. For spin waves propagating with a wave number kxin an unbounded film, the decay length is given by
= 1/kx. Since the eigenmodes of the stripe roughly corre-spond to standing waves formed by imposing the quantiza-tion condiquantiza-tion kx= k/w 共see previous section兲, the decay length corresponding to the kth resonance is k⬇w/k. TABLE I. Measured and calculated values of the frequency of the first spin-wave resonance of a 100-m-wide Permalloy stripe. The 共vertical兲 distance between the stripes and the aluminum ground plane was 1m. For the calculation we used a saturation magnetization of M = 1.15 T and an internal anisotropy of Ha= 25 Oe, as obtained from M-H loop measurements
on nonpatterned magnetic films. An experimentally obtained gyromagnetic ratio of␥⬇1.84⫻1011rad/s T was used in the simulations. For
compari-son, results of the calculation for free magnetic stripes are also included. Film thickness 共nm兲 Measured frequency 共GHz兲 Calculated frequency共GHz兲 Strip above Al ground Free stripe
100 2.15–2.38共16 devices兲 2.19 1.92
Thus, fields associated with higher modes decay faster in the vertical direction and start to behave as free-stripe modes at smaller values of d. To provide an independent verification of our method we also calculated the resonance frequency and linewidth of the stripe of Fig.2 for d = 1m using AN-SOFT HFSS, a commercial full-wave electromagnetic simula-tor. The results are shown as circles in Fig.3.
Figure4 presents the dependence of the first spin-wave resonance frequency 1
⬘
/2 and its associated extrinsic damping constant␣1ex, on the width of a FM strip. The most interesting feature of this figure concerns the behavior of␣1ex. For the particular configuration considered, the value of␣1ex reaches a maximum at w⬇8m and becomes zero as w →0 or w→⬁. This result can be qualitatively understood as follows. The extrinsic damping of spin waves is caused by the electric currents induced on the surface of the ground conductor. Those currents flow in the z-direction with a sheet density equal to the tangential magnetic field Hxon the con-ductor surface. The total loss induced in the concon-ductor is Pl=共1/2兲Rs兩Hx兩2, where Rs= Re共Zs兲 is the surface resistance of the conductor关see Eq.共A10兲兴. As spin waves are bound to the magnetic film, the magnetic field decays exponentiallyinside the dielectric layer as it reaches the ground surface. Thus, for the lowest mode, Hx⬀exp共−d/1兲=exp共−d/w兲. On the other hand, the surface resistance of the conductor is given by Rs⬃共
⬘
10/2兲1/2. The observed maximum in␣1ex as function of W is due to the competition between Rs and Hx. Increasing the stripe width leads to a lower resonance frequency and, therefore, a smaller Rs. On the other hand, it results in larger surface currents due to a larger magnetic field Hx reaching the surface of the ground conductor. In-creasing the distance d leads to less rapid rise in Hx as a function of w, shifting the maximum observed to higher val-ues of the width.
IV. CONCLUSION
The effect of the neighboring metal layer on the netic characteristics of a FM strip is analyzed by using mag-netostatic Green’s function formalism. It is shown that the metallic layer strengthens the demagnetizing field inside the FM strip and, consequently, increases the frequency of the spin-wave resonances. The finite conductivity of the metal layer leads to extra broadening of the resonance peaks, which is described in terms of a mode-dependent extrinsic damping constant. The extrinsic magnetic damping increases by re-ducing the distance between the metal ground and the FM stripe. However, it shows a more complicated behavior as function of the width of the stripe, reaching a maximum for a particular value of the stripe width.
ACKNOWLEDGMENTS
We would like to thank Y. Zhuang, P. Khalili, and J. N. Burghartz for many usefull discussions. This work was sup-ported by the Foundation for Technical Science共STW兲. FIG. 3. Resonance frequencyk⬘/2共a兲 and extrinsic damping constant␣k
ex
共b兲 vs distance d for the first four spin-wave modes 共k=1,2,3,4兲. The circles correspond to the results obtained fromANSOFT HFSS, a commercial three-dimensional full-wave electromagnetic simulator. The saturation mag-netization of the stripe is M = 1 T and the internal anisotropy is Ha= 5 Oe.
The width and the thickness of the stripe are w = 50m and t = 0.1m, respectively. The conductivity of the metal ground is= 3.3⫻107S/m.
FIG. 4. Resonance frequencyk⬘/2and extrinsic damping constant␣kexvs
the width w for the first resonance mode 共k=1兲 of a magnetic strip of thickness t = 0.1m, placed above a conductive ground plane at the distance d = 1m.
APPENDIX
In what follows, we derive Green’s function G for the magnetic potential in the presence of a ground conductor. In the regions of space filled with media with zero electric con-ductivity 共e.g., air, dielectrics, nonconductive magnetic ma-terials兲, one can employ the magnetostatic approximation and represent the magnetic field h as
h =ⵜ, 共A1兲
where is the magnetic potential. This approximation is justified if electromagnetic propagation effects are negli-gible, i.e., if the dimension of the structure to be studied is far less than that of the electromagnetic wavelength. From Maxwell’s equationⵜ·b=0, where b=0共h+m兲 is the mag-netic induction, one then arrives at the equation
ⵜ2= −ⵜ · m. 共A2兲
The magnetostatic approximation looses its validity in-side the conductive ground layer. Therefore, we solve Eq.
共A2兲outside the ground conductor only, taking into account the latter by imposing appropriate boundary conditions on at its surface. To obtain those conditions, we first note that the time-dependent variation of the magnetic field induces an electric field which, assuming uniformity along the stripe, only has a component ezin the z-direction. Using Faraday’s law and Eq.共A1兲just above the conducting plane共inside the dielectric layer兲, we have
ez
x = i0hy= i0
y. 共A3兲
Furthermore, the tangential components of the electric and magnetic fields are continuous across the dielectric-ground plane interface and are related by
ez= − ZShx, 共A4兲
where ZSis the surface impedance of the conductor.18 Com-bining Eqs.共A1兲,共A3兲, and共A5兲, we arrive at the boundary condition at the metal surface,
2 x2 = − i0 ZS y, y = − d − t/2. 共A5兲
For any distribution of magnetization m, the solution of Eq. 共2兲 can be expressed in terms of Green’s function G satisfying the equation
ⵜ2G共r,r
⬘
兲 = −␦共r − r⬘
兲. 共A6兲The functionG can be evaluated by applying a Fourier trans-form in the x-direction,
⌫共y,y
⬘
;kx兲 =冕
−⬁⬁
G共x,y;0,y
⬘
兲exp共ikxx兲dx, 共A7兲 which results in2⌫
y2 − kx
2⌫ = −␦共y − y
⬘
兲. 共A8兲The boundary condition for⌫ reads
kx2⌫ = i0
ZS
⌫
y, y = − d − t/2, 共A9兲
where the surface impedance in the Fourier domain is given by
ZS=
i0
冑
kx2+ i0. 共A10兲
Solving Eqs.共A8兲–共A10兲yields
⌫共y,y
⬘
;kx兲 = 共2兩kx兩兲−1兵exp关− 兩kx共y − y⬘
兲兩兴+ R exp关− 兩kx兩共y + y
⬘
+ 2d兲兴其, 共A11兲 where R =冑
kx 2 + i0−兩kx兩冑
kx 2 + i0+兩kx兩 . 共A12兲Green’s function 共4兲 in the space domain is finally deter-mined by applying the inverse Fourier transform,
G共r,r
⬘
兲 =冕
−⬁⬁
⌫共y,y
⬘
;kx兲exp关− ikx共x − x⬘
兲兴dkx2. 共A13兲
The one-dimensional matrix Green’s function Gញ in Eq.
共8兲 is computed from Gញ 共x,x
⬘
兲 =1 t冕
−t/2 t/2冕
−t/2 t/2 ⵜ ⵜ⬘
G共r,r⬘
兲dydy⬘
共A14兲 which results in the elementsg1共x,x
⬘
兲 = 1 t冕
−t/2 t/2冕
−t/2 t/2 y y⬘
G共r,r⬘
兲dydy⬘
=1 t再
1 2ln冋
共x − x⬘
兲2 共x − x⬘
兲2+ t2册
+ U共x − x⬘
兲 + V共x − x⬘
兲冎
g2共x,x⬘
兲 = 1 t冕
−t/2 t/2冕
−t/2 t/2 x y⬘
G共r,r⬘
兲dydy⬘
=1 t关U共x − x⬘
兲 − V共x − x⬘
兲兴, 共A15兲 where U共x − x⬘
兲 = F共1兲 + F共2兲 − 2F共3兲, V共x − x⬘
兲 = F共1*兲 + F共2*兲 − 2F共3*兲, 1= x − x⬘
+ i2共d + t兲, 2= x − x⬘
+ i2d, 3= x − x⬘
+ i2共d + t/2兲. 共A16兲The function F is defined by Eq.共2兲.
1P. H. Bryant, J. F. Smyth, S. Shultz, and D. R. Fredkin,Phys. Rev. B47,
2J. Jorzick, S. O. Demokritov, C. Mathieu, B. Hillebrands, B. Bartenlian, C.
Chappert, F. Rousseaux, and A. N. Slavin,Phys. Rev. B60, 15194共1999兲.
3Y. Roussigne, S. M. Cherif, C. Dugautier, and P. Moch,Phys. Rev. B63,
134429共2001兲.
4K. Yu. Guslienko, S. O. Demokritov, B. Hillebrands, and A. N. Slavin,
Phys. Rev. B66, 132402共2002兲.
5J. P. Park, P. Eames, D. M. Engebretson, J. Berezovsky, and P. A. Crowell,
Phys. Rev. Lett.89, 277201共2002兲.
6M. Bailleul, D. Olligs, and C. Fermon,Phys. Rev. Lett.91, 137204共2003兲. 7C. Bayer, J. P. Park, H. Wang, M. Yan, C. E. Campbell, and P. A. Crowell,
Phys. Rev. B69, 134401共2004兲.
8K. Kim and G. H. Koh, 24th International Conference on
Microelectron-ics, 2004共unpublished兲, IEEE Catalog No. 04TH8716, 377 vol. 1.
9N. Cramer, D. Lucic, R. E. Camley, and Z. Celinski,J. Appl. Phys.87,
6911共2000兲.
10Y. Zhuang, B. Rejaei, E. Boellaard, M. Vroubel, and J. N. Burghartz, IEEE
Microw. Wirel. Compon. Lett. 12, 473共2002兲.
11E. Saluhun, P. Queffelec, G. Tanne, A.-L. Adenot, and O. Acher,J. Appl.
Phys.91, 5449共2002兲.
12T. M. Crawford, M. Covington, and G. J. Parker,Phys. Rev. B67, 024411
共2003兲.
13M. Bailleul, R. Hollinger, and C. Fermon, Phys. Rev. B 73, 104424
共2006兲.
14T. Ykawa, J. Yamada, K. Abe, and J. Ikenoue,Jpn. J. Appl. Phys.16, 2187
共1977兲.
15D. F. Vaslow, Proceedings of the IEEE, Proceedings January 1973
共un-published兲, p. 142.
16B. Rejaei and M. Vroubel,J. Appl. Phys.101, 074904共2007兲.
17A. G. Gurevich and G. A. Melkov, Magnetisation Oscilations and Waves
共CRC, New York, 1996兲.