Research Article  Open Access
Propagation Effects in the SpinWave Spectrum of the Ferromagnetic Thin Film
Abstract
Exchange spin waves propagating in magnetic thin films exhibit some dynamic effects, that is, the effects caused solely by the wave propagation. In this paper we put our attention in four phenomena of such kind: the surface (and subsurface) localization, the collapse of the bulk band into a single energy level, the reversal of the mode order in the spinwave spectrum, and the dynamic separation of the thin film into two subsystems. We link these effects to properties of Hamiltonian matrix elements inherited from the spatial distribution of neighboring spins with respect to the film surface. We also provide necessary conditions for the occurrence of these effects.
1. Introduction
Magnetic thin films are the subject of intensive studies, both experimental [1–5] and theoretical [6–10], for more than half a century. This interest is stimulated by two main reasons: the huge progress in the fabrication technology of such structures [11] and the rich physics they exhibit. The latest is also inherited in more complicated systems such as magnetic multilayers [12–14] or thinfilm magnonic crystals [15–18]. A particular class of phenomena is related to the parallel to the surface propagation of spin waves (inplane propagation).
The propagation of spin waves in thin magnetic films leads to several effects including surface localization. For example, if magnetostatic interactions are taken into account famous DamonEshbach modes occur [19]. These modes are localized only at one surface, bottom or top, depending on the propagation direction. As a consequence nonreciprocity effect appears in the spinwave spectrum for thin films with one surface differing from another [20]. Both effects are dynamic and disappear for standing spin waves, that is, for inplane wave vector equal to zero. The dynamic localization appears also for purely exchange spin waves. Even in very simple model such as ferromagnetic thin film with natural surface and only nearest neighbor (NN) exchange interactions taken into account the surface localization occurs for propagating spin waves [21, 22]. In this case also another dynamic effect was observed, namely, the collapse of the bulk band, where all bulk modes are degenerated having the same energy [23].
Dynamic surface localization in thin magnetic films was studied theoretically in last decades [24–27] but it strongly gained in importance since spin polarized electron energy loss spectroscopy (SPEELS) was used to investigate propagating spin waves [28–30]. In contrast to the ferromagnetic resonance (FMR), which investigates standing spin waves, SPEELS gives possibility to scan the entire surface Brillouin zone (SBZ); thus it proves to be the powerful tool for the experimental study of propagation effects. The method allows experimental studies of the exchange interactions [31], magnetic anisotropy [32], or thermal properties [33]. For example, latest SPEELS researches show the interaction type in iron thin films to be strongly dependent on the substrate: from antiferromagnetic exchange [34] to DzyaloshinskiiMoriya one [35]. It allows also for the direct observation of spin waves localized dynamically at the surface [36].
In the planar model used in our work the thin film is treated as a set of layers representing lattice planes parallel to the surface. This allows us to reduce the problem to a onedimensional finite chain [7] with effective coupling between its elements (lattice planes). Propagation effects are strongly influenced by this effective interlayer coupling which reflects the total influence of exchange interactions from neighboring spins. The interlayer coupling was studied for standing spin waves in both dipolar and exchange coupled multilayers [37–42]. For exchange coupled layers the effective coupling depends on the wave vector of the propagating spin wave which may lead to its vanishing and the changing of its sign (between ferro and antiferromagnetic) [43].
The spinwave spectrum of different magnetic systems is strongly influenced by the range of considered interactions. For example, concurrence of the short range exchange and long range dipolar interactions generates menagerie of (meta)stable magnetic configurations [44–50] and causes rich spectrum of magnetic excitations [51–55] in magnetic dots, leads to the complete bandgap opening in magnonic crystals [56], or induces the spinwave pinning in strongly curveted wires [57]. In the case of thin films it is responsible for the occurrence of negative group velocity [58] and the splitting of the spinwave spectrum into subbands [14]. Also purely exchange spin waves exhibit new features if next nearest neighbors (NNN) interactions are considered, for example, subsurface localization [59] or separation of the film into subsystems [43].
These dynamic effects and their dependence on NNN interactions as well as on the crystallographic structure of thin film are the main objective of the present paper. We show that some of these effects could be useful for the determination of the NNN exchange interactions which is often a rather difficult task. One of more spectacular examples is the case of europium monochalcogenides (EuX). For EuS the literature data differ by the factor 2, and for EuO even the sign of remains in question [60–62]. The bulk band collapse, followed by the dynamic mode order reversal with changing wave vector, seems especially useful for this purpose.
In our work we use the planar model of the thin film with spin waves described by the Heisenberg Hamiltonian which contains NN and NNN interactions. The model has been known since the late sixties and its derivation is performed in several papers (see, e.g., [7]); thus we shall only recall its main assumptions in Section 2. In the next sections we use this model to explain dynamic (propagation) effects in the spinwave spectrum of thin films. In Section 3 we consider a surface localization of spin waves caused solely by their propagation and provide necessary conditions for this type of localization to occur. In Section 4 we study the dynamic collapse of the bulk band caused by the vanishing of the effective coupling between adjacent layers. For certain inplane wave vectors all bulk modes in the spinwave spectrum have equal energy, which makes the whole bulk band collapse into one energy level. The effect is usually followed by the reversal of the order of modes in two regions of the SBZ. In Section 5 we show that the thin film can separate into subsystems in which spin waves propagate independently for wave vectors from the border of the SBZ. Major results are summarized in Section 6. The paper is completed with two appendixes. In Appendix A we provide the neighbors distribution in the thin films considered in the paper. Appendix B contains the explicit dependence on the inplane wave vector of structural sums introduced in Section 2.
2. The Model
The schematic plot of a model is shown in Figure 1. The system under consideration is a thin magnetic film that consists of spins arranged in sites of a crystal lattice. The film can be treated as a set of crystallographic planes (layers) parallel to its surface ( plane). An external magnetic field strong enough to stabilize the homogenous ferromagnetic ground state is applied perpendicularly to the surface (along the direction). Under this assumption the Zeeman energy can be neglected as it only results in a uniform shift of the spinwave spectrum. The wave vector can be decomposed into a parallel to the film surface inplane component and perpendicular component . The first component is related to the inplane propagation of spin waves while the second one describes their quantization along the film thickness.
Considering uniform NN and NNN exchange interactions the Hamiltonian of the film takes the following form: where and are the NN and NNN exchange integrals, respectively. In the planar model of thin film the position of each spin is given by the layer number, , and the position vector within the layer, . The symbols in (1) refer to the neighbors of the spin , nearest in the first summation and next nearest in the second one. The spin from the plane can have neighbors in the same lattice plane (), in neighboring planes (), or in further planes.
The detailed description of the diagonalization procedure of Hamiltonian (1) is given in several papers (see, e.g., [7]). The first step is the HolsteinPrimakoff transformation, that is, the mapping to boson operators. Then the Fourier transformation in the plane of the film is performed and, subsequently, the BogolyubovTyablikov transformation in the direction perpendicular to it (mapping to the wavevector space). The assumptions the method is based on make it valid for low concentration of the spin waves (low temperatures in the case of thermal magnons) and for thin films large in the directions parallel to the surface.
The final outcome is an band matrix, where is the number of layers forming the thin film, in the following form:The number of nonzero elements depends on the spatial distribution of the NN and NNN (provided in Appendix A) and their explicit forms depend on the exchange integrals and , as well as on the inplane wave vector.
Matrix (2) is typical for onedimensional chain where the interaction between its links is described by the offdiagonal elements of the matrix. Other entries of the Hamiltonian matrix have the following meanings: represents the onsite energy; are surface (subsurface) parameters (see, e.g., [7] for more details). In other words, the method used transforms thin film into a onedimensional chain, where the role of chain elements is played by lattice planes parallel to the surface. This means the offdiagonal entries of the matrix describe the effective coupling between lattice planes: refers to the interaction between nearest layers, to the interaction between next nearest ones, and so forth. In our microscopic model this effective interlayer coupling reflects the total influence of exchange interactions from neighboring spins taking into account NN and NNN. The coupling depends not only on exchange integrals and the surface cut (NN and NNN distribution) but also on the inplane wave vector which is the reason for the occurrence of propagation (dynamic) effects in the spinwave spectrum. Thus the knowledge of the dependence of offdiagonal elements on the inplane wave vector allows forecasting the propagation effects. It is also worth noticing that the Zeeman energy, as an onsite energy, only affects the diagonal terms, shifting the whole spectrum of eigenvalues in a uniform way.
During the diagonalization of the Hamiltonian it is useful to introduce the socalled structural sums, which describe the spatial distribution of the NN and NNN [7], defined as follows:where summation over runs over the projections of (next) nearest neighbors of the spin in the plane belonging to the plane . The explicit forms of structural sums for surface cuts considered in this paper are given in Appendix B. Elements of the Hamiltonian matrix expressed in terms of structural sums are the following:where . Please notice that dependence of Hamiltonian matrix elements on the inplane wave vector is introduced only by structural sums. Therefore the evolution of these sums with determines propagation effects.
3. Dynamic Surface Localization
Surface localization of spin waves in thin magnetic films is widely studied in literature both experimentally and theoretically. For exchange spin waves it was found that the key factor for the surface localization is geometric disposition of interacting neighbors with respect to the surface rather than the interaction range: the emergence of surface states is due to interactions with neighbors the position vectors of which are oblique to the film surface, irrespectively of whether these “oblique” interactions are NN or NNN type [21, 22, 63]. Here we systemize the results and link them to the properties of Hamiltonian matrix elements. We also show that the range of interactions is important as well, however not the direct spinspin interactions but the effective interlayer ones.
The geometric disposition of neighboring spins is particularly important in the case of purely dynamical localization, that is, for the thin film with natural surfaces (no pinning). In such system there is no surface localization for the SBZ center but it could occur for nonzero inplane wave vectors (propagating spin waves). To obtain surface localization for (which can be detected by FMR experiments) one needs to introduce surface anisotropy to generate additional pinning at the surface. The anisotropy has large influence on the surface localization and can cover the dynamical localization effect. Therefore it is not taken into account in this paper. (The influence of the anisotropy on the spinwave localization is studied, e.g., in [64].)
The simplest case is sc(001) surface cut for which the Hamiltonian matrix is threediagonal only. This is a consequence of the NN and NNN disposition. Every spin has 4 NN and 4 NNN at the same lattice plane; at the adjacent plane it has 1 NN and 4 NNN (see Appendix A). This means the effective coupling is restricted only to the neighboring lattice planes even if we include NNN interactions. In such case, that is, if the Hamiltonian matrix is threediagonal, conditions for surface localization are [7]for symmetric and antisymmetric mode, respectively. Here is offdiagonal element and is the surface parameter. Please notice that the surface parameter according to (6) depends only on the number of NN and NNN at the adjacent plane so the dynamic surface localization can be induced solely by through its dependence on structural sums (and then on inplane wave vector).
According to the definition (3) if the position vector of the neighboring spin from adjacent layer is perpendicular to the surface (we will call it perpendicular neighbor) then the corresponding structural sum is equal to 1 (the number of adequate neighbors) independently of . In this case and neither condition (7) nor (8) is fulfilled and thus surface states cannot occur. In a consequence the necessary condition for the dynamic localization is an oblique disposition of the neighbors from adjacent lattice plane, which makes offdiagonal elements of the Hamiltonian matrix functions of the inplane wave vector. For example, for sc(001) surface cut there is only one “perpendicular” nearest neighbor at the layer and thus corresponding structural sum . The structural sum for NNN is given by which means these 4 “oblique” neighbors form layer are responsible for eventual surface localization. Thus, the offdiagonal element is given bywhere and are components of the inplane wave vector .
In Figure 2 we show exemplary spinwave spectra for sc(001) thin films consisting of lattice planes. The spectra are plotted along the highsymmetry path marked as a green solid line in the surface Brillouin zone (SBZ) shown in the inset in Figure 2(a). In the first case NNN interactions are assumed to be ferromagnetic (FM) with the ratio of NNN to NN exchange integral set to (Figure 2(a)). The second spectrum (Figure 2(b)) is obtained for antiferromagnetic (AFM) NNN interactions with . In both graphs the shaded area represents the bulk band. Below each spectrum we provide a graph with plotted as a green dotted line accompanied with horizontal blue lines for (dashed lines) and (solid lines). In Figures 2(c)–2(h) we show spinwave profiles for some particular wave vectors marked with purple lines in the corresponding spectra.
For case surface parameter and condition (7) means that ( should be between dashed lines in Figure 2(a)) which is fulfilled for any ; that is, there is a symmetric surface state within the whole SBZ except point . In Figure 2(c) we show spinwave profiles for the wave vector between points and in the SBZ. The lowest mode is localized at both surfaces while all other modes are bulk modes. The second surface mode, antisymmetric, occurs for wave vector large enough to fulfill condition (8) which in the case in question means ( between solid lines in Figure 2(a)). The points at which this mode becomes localized are marked by red arrows in Figure 2(a). In Figure 2(d) we show spinwave profiles for point in the SBZ where two lowest modes are localized. (The third profile given in panel (e) shows the mode order reversal whose effect is explained in Section 4.)
For case and the surface localization of the symmetric mode appears in some region in the SBZ around its center (again, except point ), that is, if is between dashed lines in Figure 2(b). For larger there is no surface state in the spectrum. Condition (6) for antisymmetric mode localization boils down to ( between solid lines in Figure 2(b)) which is fulfilled for wave vectors from some ring around point (between points marked by red arrows, cf. profiles in Figures 2(f) and 2(g)).
Other surface cuts for which the Hamiltonian matrix is threediagonal are bcc(110) and fcc(111) but now neighbors of both types, NN and NNN, which are located on the adjacent layer are disposed obliquely with respect to the surface plane. Therefore, both interactions are responsible for the dynamic surface localization.
Concerning the influence of the NN and NNN on the surface localization very instructive are surface cuts bcc(001) and fcc(001). In both cases there are 4 NN located “obliquely” at the adjacent plane () and one “perpendicular” NNN at nextneighboring plane (). This means only NN are responsible for spinwave localization. In other words, for the hypothetical situation when NN interactions are switched off () and only NNN interactions are taken into account () there will be no surface states in the spinwave spectrum. Exemplary spectra are shown in Figure 3 for bcc(001) thin film consisting of 11 lattice planes () with only one type of interactions taken into account: NN in panel (a) and NNN in panel (b). For the first case and the localization condition for symmetric mode is fulfilled for any inplane wave vector except 0 ( above dashed blue line in Figure 3(a)). For antisymmetric mode the condition is ; thus the second surface state occurs if is above solid blue line (between points marked by red arrows). In Figure 3(b), that is, for the hypothetical case of only NNN interactions switched on, there is no surface localization due to the perpendicular disposition of the NNN from the plane .
(a) ,
(b) ,
In the very last case the Hamiltonian matrix is fivediagonal; thus we should take into account both surface parameters and , and simple conditions (7) and (8) for the surface localization are not valid anymore. This also means that the effective coupling between layers now extends up to the secondneighboring layer which is described by offdiagonal element . For bcc(001) surface cut does not depend on because at the layer there is only one “perpendicular” neighbor. Different situation holds for fcc(110) surface cut for which there are two “oblique” NNN at the layer (and four “oblique” NN at ). Hence, we expect surface localization caused by NN and additional localization at one layer below surface (subsurface localization) caused by NNN. In Figure 4 we show two spinwave spectra for fcc(110) thin film consisting of 11 layers without (panel (a)) and with (panel (b)) NNN interactions. In the first case (Figure 4(a)) there are two surface states emerging from the bulk band in some region of the SBZ (please notice that even if NNN interactions are switched off the matrix is fivediagonal and conditions (7) and (8) are not valid). If NNN interactions are taken into account (Figure 4(b)) two other modes go out of the bulk band. The profiles of the lowest modes are shown in Figure 4(c): two lowest modes ( = 1, 2) are localized at the surface while the next two modes ( = 3, 4) are localized at the first plane below the surface layer. This subsurface localization is induced by the spinwave propagation due to the oblique disposition of NNN from the layer .
(a)
(b)
(c)
It is also worth noticing that each surface mode is localized at both surfaces. Moreover, in the case of symmetric thin film studied here the localization of particular mode is the same at both surfaces (the spinwave profile is symmetric or antisymmetric). To differ one surface from the other one should introduce additional pinning. It could be done by taking into account the surface anisotropy [64]. However, anisotropy terms in the elements of the Hamiltonian matrix do not depend on the wave vector (see (19) and (20) in [64]). Therefore we do not observe such effects like nonreciprocity found in the case of DamonEshbach modes [20].
4. Band Collapse and Mode Order Reversal
Let us go back to the case for which the Hamiltonian matrix is threediagonal only. Since offdiagonal element describes the effective coupling between neighboring layers only nearest layers are coupled. For the coupling is of ferromagnetic nature, that is, it favors the inphase precession of spins belonging to neighboring lattice planes. By analogy with phonons [65], this type of excitation is referred to as acoustic. On the other hand, if the coupling has antiferromagnetic nature energetically favoring antiphase precession (if the external field saturates the sample; if not the situation is opposite). This type of excitation is called optic. In the case of phonons the acoustic mode has a lower energy than the optic mode and the same pattern is considered the normal order of modes in the spinwave spectrum. If the energy of the optic mode is lower than that of the acoustic one, the spinwave spectrum has a reversed mode order, which occurs as a consequence of AFM interaction between adjacent layers. (This problem is widely discussed for standing spin waves for both types of interlayer coupling, exchange and dipolar; see, e.g., [37, 39–42].) However, for propagating spin waves depends on inplane wave vector which implies evolution of effective coupling within the SBZ. This means the character of effective coupling may change from FM to AFM for different wave vectors and, consequently, the order of modes in the spinwave spectrum may change as well.
Such dynamic (due to the wave propagation) mode order reversal we observe in spinwave spectra shown in Figure 2, that is, for sc(001) thin film composed of layers. In Figure 2(a) we show the spectrum for FM NNN interactions () accompanied by the dependence of . Going along highsymmetry path in the SBZ (green line in the inset in Figure 2(a)) at point the effective interlayer coupling is FM () until point is reached. In Figures 2(c) and 2(d) we show spinwave profiles for two points from this section of the highsymmetry path marked in panel (a) with purple lines. For both of them we observe normal order of modes with lowest mode acoustic and highest one optic. After crossing point offdiagonal element becomes greater than zero and effective coupling is AFM. This results in reversed mode order shown in Figure 2(e) where spinwave profiles for point in the SBZ are plotted. Going further along the highsymmetry path restores its negative value after some point between and . In other words, the reversed mode order occurs in some region of SBZ around point .
A similar feature we observe for AFM NNN interactions for which case the spinwave spectrum is presented in Figure 2(b). The main difference is that reversed mode order appears around point now (see the profiles provided in Figures 2(f) and 2(g)) because in this region . After crossing point effective coupling is FM () and the order of modes is normal (see spinwave profiles for point in Figure 2(h)).
These two regions, of normal and reversed order of modes, are separated by line at which (see dependence on Figures 2(a) and 2(b)). In such situation Hamiltonian matrix has nonzero elements only at its diagonal. In a consequence there are eigenvalues equal to (bulk modes) and two eigenvalues equal to (surface modes). This is what we call band collapse: whole bulk band is collapsed into one energy level. The collapse was discovered by Lévy et al. [23] in a semiinfinite ferromagnet of sc structure and (001) surface. Originally this effect was described as caused by NNN interactions. In [66] we studied this effect in thin films with different surface cuts and found it to occur even without NNN interactions. In [43] we show that the collapse may be followed by the reversal of the order of modes in the spinwave spectrum.
Since depends on exchange integrals for NN and NNN also two effects described above should. For sc(001) surface cut element is given by (9) and the equation has real solutions only if . Therefore, if NNN interactions are too weak there is no bulk band collapse nor dynamical change of the mode order in the spinwave spectrum of sc(001) films. For weak FM NNN interactions is negative and the order of modes in the spectrum is normal within the whole SBZ. On the other hand, for weak AFM NNN interactions the situation is opposite: and the order of modes is reversed for any . For stronger NNN interactions the solutions for different ratios are plotted in Figure 5 in a quarter of SBZ. For the critical value the equation is satisfied only by one inplane wave vector, corresponding to the highsymmetry point or for FM (Figure 5(a)) or AFM (Figure 5(b)) NNN interactions, respectively. This means that only at this one point does collapse occur and there is no dynamical change of the mode order. For NNN interactions stronger than the critical value solutions of the equation form a line in the SBZ which shifts from the corresponding highsymmetry point toward the center of the SBZ quarter. Thus, as we already mentioned, the region of the reversed mode order appears around the point for FM NNN interactions and around point for AFM ones under the condition that NNN interactions are strong enough.
(a)
(b)
A similar situation we have for the bcc(110) surface cut; the Hamiltonian matrix is also threediagonal. In this case the SBZ has the shape of an irregular hexagon which results in rather complicated formulas for the wave vector within the SBZ and, consequently, for the elements of the Hamiltonian matrix. To simplify the description of the problem the socalled uniform zone can be introduced [43], based on new coordinates, and , assumed to range from −0.5 to 0.5 within the SBZ. With such definition the uniform zone has the shape of a square regardless of the original SBZ geometry. For bcc(110) case the onetoone correspondence that maps each point in the SBZ to a point in the uniform zone is given by [43] The offdiagonal element of the Hamiltonian matrix expressed in coordinates and is given by Solutions of the equation for FM NNN interactions are shown in Figure 6. In contrast to sc(001) case there is no threshold of ratio and the collapse occurs even without NNN interactions. For it happens at the border of the SBZ (segment ) and is not followed by the mode order reversal. While FM NNN interactions increase the line is moving towards point and the region of reversed mode order appears. For AFM NNN interactions a necessary condition for the mode order reversal is , which means the NNN interactions must be stronger than NN ones.
The effect of band collapse followed by mode order reversal should be observable in experiments which gives possibility to look at spin waves over the whole SBZ, for example, by spin polarized electron energy loss spectroscopy (SPEELS) [30]. Observation of this effect gives possibility to determine precisely the ratio. For example, in sc(001) thin film the eventual collapse appears for some point in between and . From relation (9) for the wave vector from the segment  the condition for the collapse is fulfilled if where is the length of the wave vector along the  segment for which the collapse occurs and is the lattice constant.
For bcc(110) surface cut and FM NNN interactions if the wave vector changes along the line (the green dashed line in Figure 6) the condition for the collapse is . Thus the observation of the collapse for some particular directions in the SBZ gives automatically the ratio of NNN to NN interactions.
For other surface cuts the Hamiltonian matrix is threediagonal if NNN interactions are absent and then the collapse occurs only at the border of the SBZ. For example, for sc(110) thin films with only NN interactions taken into account the offdiagonal element of the Hamiltonian matrix is given by which means if , which is satisfied by wave vectors from the two borders of the SBZ (bold lines in Figure 7(a)). For fcc(001) surface cut and we have and thus for the whole SBZ border (, bold lines in Figure 7(b)). The same situation, that is, the collapse of the bulk band at the SBZ border, exists also for bcc(001) surface cut (see segment  in Figure 3(a)) for which
(a) sc(110)
(b) fcc(001)
5. Separation into Subsystems
If the Hamiltonian matrix is fivediagonal and the offdiagonal element depends on it could happen that and for some inplane wave vectors. In such case, according to the physical meaning of the offdiagonal elements, there is no effective coupling between neighboring layers while nextneighboring layers are coupled. In other words, from the point of view of propagating spin waves the system is separated into two subsystems consisting of the layers with odd and even numbers, respectively.
In Figure 8 we show the spinwave spectrum of bcc(001) thin film with FM (panel (a)) and AFM (panel (b)) NNN interactions. For bcc(001) surface cut there are 4 “oblique” NN at the adjacent layer () and one “perpendicular” NNN at layer . This means that depends on NNN interactions only and does not depend on the spinwave propagation. It has FM character for FM NNN interactions (, Figure 8(a)) and AFM character for AFM NNN interactions (, Figure 8(b)). Element depends only on NN and ; thus in both cases shown in Figure 8 its evolution with is the same: at the SBZ border, that is, between points and . For both cases we also show spinwave profiles (Figures 8(c)–8(f)) for two points in the SBZ: and (marked with red lines in the corresponding spectrum). For point both spectra are typical for FM thin film but for point we can clearly see the separation of the thin film into two subsystems. For two lowest modes, localized at surfaces, the nonzero dynamical magnetization (the amplitude of the spin precession) occurs only at layers with odd numbers. In the bulk band the spinwave excitations are restricted to only one type of layers: even or oddnumbered. For example, the third mode in Figure 8(d) is the uniform mode excited exclusively in the subsystem consisting of only evennumbered layers.
If the mode order in the spinwave spectrum is determined by the sign of (cf. Section 4). If the modes appear in normal order in the spectrum (for each subsystem, Figure 8(d)) while for the mode order is reversed (Figure 8(f)). Since the nextneighboring layer coupling is rather weak this effect is removed if .
If both and depend on they can vanish simultaneously for some inplane wave vectors. Such situation occurs for fcc(110) surface cut but only for one value of the ratio, namely, for 0.5 (Figure 9). In this case the collapse occurs at the SBZ border where both offdiagonal elements and are equal to zero due to the spinwave propagation.
6. Conclusion
In magnetic thin films some interesting phenomena occur just due to the spinwave propagation. In this paper we collect such dynamic effects in the case of exchange spin waves and link them to properties of Hamiltonian matrix elements. First of all, we study the dynamically induced (due to the wave propagation) surface localization which appears only if neighboring spins are located obliquely with respect to the surface of the film regardless of the circumstance of the interactions. We show that the same rule holds for the subsurface localization as well.
Another effect is the collapse of the bulk band into one energy level for some certain wave vectors caused by the vanishing of the effective coupling between adjacent layers of spins. The collapse, usually followed by the reversal of the mode order in the spinwave spectrum, can be utilized for precise determination of the next nearest neighbors exchange interactions. If the effective interlayer coupling ranges up to the nextneighboring layers the spin waves from the border of the surface Brillouin zone may exhibit another phenomenon: they can be excited in only every second lattice plane; that is, the thin film acts like that divided into two subsystems: one consists of oddnumbered layers and the second one consists of evennumbered ones.
Appendices
A. Nearest and Next Nearest Neighbors Distribution
Here we present an analysis of the nearest (NN) and next nearest (NNN) neighbors distribution in thin films discussed in this paper. In Figures 10–12 the 3D distribution is provided as well as 2D mapping onto the lattice plane parallel to the surface of the film. In each figure an arbitrary chosen spin is marked as the central black circle and the shaded plane represents the plane parallel to the film surface which contains the chosen spin. Below each 3D distribution the 2D mapping of respective neighbors is shown, that is, their projection onto the plane , where shaded circles (rings) represent NN (NNN) from the same plane while empty dashed ones represent the projections of respective neighbors from adjacent planes ().
(a) sc(001), NN
(b) sc(110), NN
(c) sc(111), NN
(d) sc(001), NNN
(e) sc(110), NNN
(f) sc(111), NNN
(a) bcc(001), NN
(b) bcc(110), NN
(c) bcc(111), NN
(d) bcc(001), NNN
(e) bcc(110), NNN
(f) bcc(111), NNN
(a) fcc(001), NN
(b) fcc(110), NN
(c) fcc(111), NN
(d) fcc(001), NNN
(e) fcc(110), NNN
(f) fcc(111), NNN
The numbers of neighbors within the individual lattice planes parallel to the surface are summarized in Table 1.

B. Structural Sums
In this appendix we provide the explicit dependence on the inplane wave vector of structural sums for thin films considered in the paper for NN (Table 2) and NNN (Table 3). For the detailed calculations please see [67].
 
Symbol: . 
 
Symbol: . 
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The author thanks JeanClaude Serge Lévy for valuable discussions. This study received support from the National Science Centre (NCN) of Poland, Project DEC2012/07/E/ST3/00538. The author acknowledges the financial support of the Operational Programme Human Capital, cofinanced by the European Social Fund in Poland, Grant no. POKL.04.01.0100019/10.
References
 M. J. Klein and R. S. Smith, “Thin ferromagnetic films,” Physical Review, vol. 81, no. 3, pp. 378–380, 1951. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. E. Wigen, C. F. Kooi, M. R. Shanabarger, and T. D. Rossing, “Dynamic pinning in thinfilm spinwave resonance,” Physical Review Letters, vol. 9, no. 5, pp. 206–208, 1962. View at: Publisher Site  Google Scholar
 P. Grünberg, J. Barnas, F. Saurenbach, J. A. Fuss, A. Wolf, and M. Vohl, “Layered magnetic structures: antiferromagnetic type interlayer coupling and magnetoresistance due to antiparallel alignment,” Journal of Magnetism and Magnetic Materials, vol. 93, pp. 58–66, 1991. View at: Publisher Site  Google Scholar
 J. Dubowik, K. Załȩski, H. Głowiński, and I. Gociańska, “Angular dependence of ferromagnetic resonance linewidth in thin films,” Physical Review B, vol. 84, no. 18, Article ID 184438, 2011. View at: Publisher Site  Google Scholar
 P. Kuświk, P. L. Gastelois, M. M. Soares et al., “Effect of CoO/Ni orthogonal exchange coupling on perpendicular anisotropy of Ni films on Pd(001),” Physical Review B, vol. 91, no. 13, Article ID 134413, 2015. View at: Publisher Site  Google Scholar
 A. R. Ferchmin, “A quantum theory of spinwave resonance in cubic ferromagnets,” Physics Letters, vol. 1, no. 7, pp. 281–283, 1962. View at: Publisher Site  Google Scholar
 H. Puszkarski, “Theory of surface states in spin wave resonance,” Progress in Surface Science, vol. 9, no. 56, pp. 191–247, 1979. View at: Publisher Site  Google Scholar
 I. Zasada, B. Busiakiewicz, and L. Wojtczak, “Shortrange order in Valenta model of thin films,” Journal of Magnetism and Magnetic Materials, vol. 312, no. 1, pp. 58–64, 2007. View at: Publisher Site  Google Scholar
 D.T. Hoang, M. Kasperski, H. Puszkarski, and H. T. Diep, “Reorientation transition in molecular thin films: potts model with dipolar interaction,” Journal of Physics Condensed Matter, vol. 25, no. 5, Article ID 056006, 2013. View at: Publisher Site  Google Scholar
 H. Puszkarski and P. Tomczak, “Spinwave resonance in (Ga,Mn)As thin films: probing inplane surface magnetic anisotropy,” Physical Review B, vol. 91, no. 19, Article ID 195437, 2015. View at: Publisher Site  Google Scholar
 C. A. F. Vaz, J. A. C. Bland, and G. Lauhoff, “Magnetism in ultrathin film structures,” Reports on Progress in Physics, vol. 71, no. 5, Article ID 056501, 2008. View at: Publisher Site  Google Scholar
 M. Krawczyk, H. Puszkarski, J.C. S. Lévy, S. Mamica, and D. Mercier, “Theoretical study of spin wave resonance filling fraction effect in composite ferromagnetic [ABA] trilayer,” Journal of Magnetism and Magnetic Materials, vol. 246, no. 12, pp. 93–100, 2002. View at: Google Scholar
 W. P. Zhou and G. H. Yun, “The energy bands and the dispersion relations of the spin waves in a ferromagnetic bilayer system,” Surface Science, vol. 553, no. 1–3, pp. 75–81, 2004. View at: Publisher Site  Google Scholar
 S. Pal, B. Rana, S. Saha et al., “Timeresolved measurement of spinwave spectra in CoO capped [Co(t)/Pt(7Å)]_{n−1} Co(t) multilayer systems,” Journal of Applied Physics, vol. 111, no. 7, Article ID 07C507, 2012. View at: Publisher Site  Google Scholar
 V. V. Kruglyak, P. S. Keatley, A. Neudert, R. J. Hicken, J. R. Childress, and J. A. Katine, “Imaging collective magnonic modes in 2D arrays of magnetic nanoelements,” Physical Review Letters, vol. 104, no. 2, Article ID 027201, 2010. View at: Publisher Site  Google Scholar
 S. Tacchi, G. Duerr, J. W. Klos et al., “Forbidden band gaps in the spinwave spectrum of a twodimensional bicomponent magnonic crystal,” Physical Review Letters, vol. 109, no. 13, Article ID 137202, 2012. View at: Publisher Site  Google Scholar
 S. Mamica, M. Krawczyk, and J. W. Kos, “Spinwave band structure in 2D magnonic crystals with elliptically shaped scattering centres,” Advances in Condensed Matter Physics, vol. 2012, Article ID 161387, 6 pages, 2012. View at: Publisher Site  Google Scholar
 M. L. Sokolovskyy, J. W. Klos, S. Mamica, and M. Krawczyk, “Calculation of the spinwave spectra in planar magnonic crystals with metallic overlayers,” Journal of Applied Physics, vol. 111, no. 7, Article ID 07C515, 2012. View at: Publisher Site  Google Scholar
 R. W. Damon and J. R. Eshbach, “Magnetostatic modes of a ferromagnet slab,” Journal of Physics and Chemistry of Solids, vol. 19, no. 34, pp. 308–320, 1961. View at: Publisher Site  Google Scholar
 M. Mruczkiewicz, M. Krawczyk, G. Gubbiotti et al., “Nonreciprocity of spin waves in metallized magnonic crystal,” New Journal of Physics, vol. 15, Article ID 113023, 2013. View at: Publisher Site  Google Scholar
 H. Puszkarski, J. C. S. Lévy, and S. Mamica, “Does the generation of surface spinwaves hinge critically on the range of neighbour interaction?” Physics Letters A, vol. 246, no. 34, pp. 347–352, 1998. View at: Google Scholar
 S. Mamica, R. Józefowicz, and H. Puszkarski, “The role of obliquetosurface disposition of neighbours in the emergence of surface spin waves in magnetic films,” Acta Physica Polonica A, vol. 94, no. 1, pp. 79–91, 1998. View at: Google Scholar
 J. C. Lévy, E. Ilisca, and J. L. Motchane, “Influence of surface anisotropy and nextnearestneighbor coupling on surface spin waves,” Physical Review B, vol. 5, no. 3, pp. 1099–1105, 1972. View at: Publisher Site  Google Scholar
 H. Puszkarski, “Surface spin waves in hexagonal cobalt thin films,” Physica Status Solidi (B), vol. 22, no. 2, pp. 355–359, 1967. View at: Publisher Site  Google Scholar
 J. C. S. Lévy, “Surface and interface magnons: magnetic structures near the surface,” Surface Science Reports, vol. 1, no. 2, pp. 39–119, 1981. View at: Publisher Site  Google Scholar
 B. X. Xu, M. Mostoller, and A. K. Rajagopal, “Spin waves at surfaces and interfaces in cubic Heisenberg systems,” Physical Review B, vol. 31, no. 11, pp. 7413–7423, 1985. View at: Publisher Site  Google Scholar
 J. M. Wesselinowa, “Influence of magnetic surface anisotropy on the dynamic properties in ferromagnetic thin films,” Journal of Physics Condensed Matter, vol. 18, no. 35, pp. 8169–8178, 2006. View at: Publisher Site  Google Scholar
 M. Plihal, D. L. Mills, and J. Kirschner, “Spin wave signature in the spin polarized electron energy loss spectrum of ultrathin Fe films: theory and experiment,” Physical Review Letters, vol. 82, no. 12, pp. 2579–2582, 1999. View at: Publisher Site  Google Scholar
 H. Ibach, D. Bruchmann, R. Vollmer, M. Etzkorn, P. S. A. Kumar, and J. Kirschner, “A novel spectrometer for spinpolarized electron energyloss spectroscopy,” Review of Scientific Instruments, vol. 74, no. 9, pp. 4089–4095, 2003. View at: Publisher Site  Google Scholar
 R. Vollmer, M. Etzkorn, P. S. Anil Kumar, H. Ibach, and J. Kirschner, “Spinpolarized electron energy loss spectroscopy of high energy, large wave vector spin waves in ultrathin fcc Co films on Cu(001),” Physical Review Letters, vol. 91, no. 14, p. 147201/4, 2003. View at: Google Scholar
 K. Zakeri, T.H. Chuang, A. Ernst et al., “Direct probing of the exchange interaction at buried interfaces,” Nature Nanotechnology, vol. 8, no. 11, pp. 853–858, 2013. View at: Publisher Site  Google Scholar
 T.H. Chuang, K. Zakeri, A. Ernst et al., “Magnetic properties and magnon excitations in Fe(001) films grown on Ir(001),” Physical Review B, vol. 89, no. 17, Article ID 174404, 2014. View at: Publisher Site  Google Scholar
 K. Zakeri, J. Prokop, Y. Zhang, and J. Kirschner, “Magnetic excitations in ultrathin magnetic films: temperature effects,” Surface Science, vol. 630, pp. 311–316, 2014. View at: Publisher Site  Google Scholar
 Y. Meng, Kh. Zakeri, A. Ernst et al., “Direct evidence of antiferromagnetic exchange interaction in Fe(001) films: strong magnon softening at the highsymmetry M point,” Physical Review B, vol. 90, no. 17, Article ID 174437, 2014. View at: Google Scholar
 K. Zakeri, Y. Zhang, J. Prokop et al., “Asymmetric spinwave dispersion on Fe(110): direct evidence of the dzyaloshinskiimoriya interaction,” Physical Review Letters, vol. 104, no. 13, Article ID 137203, 2010. View at: Publisher Site  Google Scholar
 K. Zakeri, Y. Zhang, and J. Kirschner, “Surface magnons probed by spinpolarized electron energy loss spectroscopy,” Journal of Electron Spectroscopy and Related Phenomena, vol. 189, pp. 157–163, 2013. View at: Publisher Site  Google Scholar
 P. Grünberg, “Magnetostatic spinwave modes of a heterogeneous ferromagnetic double layer,” Journal of Applied Physics, vol. 52, no. 11, pp. 6824–6829, 1981. View at: Publisher Site  Google Scholar
 J. Barnaś, “Interlayer exchange coupling in ultrathin layered structures,” Journal of Magnetism and Magnetic Materials, vol. 128, no. 12, pp. 171–178, 1993. View at: Publisher Site  Google Scholar
 R. W. Wang and D. L. Mills, “Magnetic properties of finite antiferromagnetic superlattices: statics, dynamics, and the surface spinflop phase,” Physical Review B, vol. 50, no. 6, pp. 3931–3941, 1994. View at: Publisher Site  Google Scholar
 S. M. Rezende, C. Chesman, M. A. Lucena, A. Azevedo, F. M. De Aguiar, and S. S. P. Parkin, “Studies of coupled metallic magnetic thinfilm trilayers,” Journal of Applied Physics, vol. 84, no. 2, pp. 958–972, 1998. View at: Publisher Site  Google Scholar
 R. Zivieri, L. Giovannini, and F. Nizzoli, “Acoustical and optical spin modes of multilayers with ferromagnetic and antiferromagnetic coupling,” Physical Review B, vol. 62, no. 22, Article ID 14950, 2000. View at: Publisher Site  Google Scholar
 R. Zivieri, “Metamaterial properties of onedimensional and twodimensional magnonic crystals,” Solid State Physics, vol. 63, pp. 151–216, 2012. View at: Publisher Site  Google Scholar
 S. Mamica, “Dynamic effects on the spinwave spectrum of the bcc thin film,” The European Physical Journal B, vol. 87, no. 12, 2014. View at: Publisher Site  Google Scholar
 E. Y. Vedmedenko, A. Ghazali, and J.C. S. Lévy, “Magnetic structures of Ising and vector spins monolayers by MonteCarlo simulations,” Surface Science, vol. 402404, pp. 391–395, 1998. View at: Publisher Site  Google Scholar
 T. Shinjo, T. Okuno, R. Hassdorf, K. Shigeto, and T. Ono, “Magnetic vortex core observation in circular dots of permalloy,” Science, vol. 289, no. 5481, pp. 930–932, 2000. View at: Publisher Site  Google Scholar
 K. L. Metlov and Y. P. Lee, “Map of metastable states for thin circular magnetic nanocylinders,” Applied Physics Letters, vol. 92, no. 11, Article ID 112506, 2008. View at: Publisher Site  Google Scholar
 S. Mamica, J.C. S. Lévy, P. Depondt, and M. Krawczyk, “The effect of the singlespin defect on the stability of the inplane vortex state in 2D magnetic nanodots,” Journal of Nanoparticle Research, vol. 13, no. 11, pp. 6075–6083, 2011. View at: Publisher Site  Google Scholar
 S. Mamica, J.C. S. Lévy, M. Krawczyk, and P. Depondt, “Stability of the Landau state in square twodimensional magnetic nanorings,” Journal of Applied Physics, vol. 112, no. 4, Article ID 043901, 2012. View at: Publisher Site  Google Scholar
 L. Rondin, J.P. Tetienne, S. Rohart et al., “Strayfield imaging of magnetic vortices with a single diamond spin,” Nature Communications, vol. 4, article 2279, 2013. View at: Publisher Site  Google Scholar
 X. Li, W. V. Liu, and L. Balents, “Spirals and skyrmions in two dimensional oxide heterostructures,” Physical Review Letters, vol. 112, no. 6, Article ID 067202, 2014. View at: Publisher Site  Google Scholar
 O. G. Heinonen, D. K. Schreiber, and A. K. PetfordLong, “Micromagnetic modeling of spinwave dynamics in exchangebiased permalloy disks,” Physical Review B, vol. 76, no. 14, Article ID 144407, 2007. View at: Publisher Site  Google Scholar
 L. Giovannini, F. Montoncello, R. Zivieri, and F. Nizzoli, “Spin excitations in nanometric magnetic dots: calculations and comparison with light scattering measurements,” Journal of Physics Condensed Matter, vol. 19, no. 22, Article ID 225008, 2007. View at: Publisher Site  Google Scholar
 Z. Liu, R. D. Sydora, and M. R. Freeman, “Shape effects on magnetization state transitions in individual 160nm diameter Permalloy disks,” Physical Review B, vol. 77, no. 17, Article ID 174410, 2008. View at: Publisher Site  Google Scholar
 F. Montoncello, L. Giovannini, F. Nizzoli, R. Zivieri, G. Consolo, and G. Gubbiotti, “Spinwave activation by spinpolarized current pulse in magnetic nanopillars,” Journal of Magnetism and Magnetic Materials, vol. 322, no. 16, pp. 2330–2334, 2010. View at: Publisher Site  Google Scholar
 S. Mamica, J.C. S. Lévy, and M. Krawczyk, “Effects of the competition between the exchange and dipolar interactions in the spinwave spectrum of twodimensional circularly magnetized nanodots,” Journal of Physics D: Applied Physics, vol. 47, no. 1, Article ID 015003, 2014. View at: Publisher Site  Google Scholar
 S. Mamica, M. Krawczyk, M. L. Sokolovskyy, and J. RomeroVivas, “Large magnonic band gaps and spectra evolution in threedimensional magnonic crystals based on magnetoferritin nanoparticles,” Physical Review B, vol. 86, no. 14, Article ID 144402, 2012. View at: Publisher Site  Google Scholar
 V. S. Tkachenko, A. N. Kuchko, M. Dvornik, and V. V. Kruglyak, “Propagation and scattering of spin waves in curved magnonic waveguides,” Applied Physics Letters, vol. 101, no. 15, Article ID 152402, 2012. View at: Publisher Site  Google Scholar
 M. Krawczyk, S. Mamica, M. Mruczkiewicz et al., “Magnonic band structures in twodimensional bicomponent magnonic crystals with inplane magnetization,” Journal of Physics D: Applied Physics, vol. 46, no. 49, Article ID 495003, 2013. View at: Publisher Site  Google Scholar
 S. Mamica, H. Puszkarski, and J. C. S. Lévy, “The role of nextnearest neighbours for the existence conditions of subsurface spin waves in magnetic films,” Physica Status Solidi B, vol. 218, no. 2, pp. 561–569, 2000. View at: Publisher Site  Google Scholar
 S. Mühlbauer, P. Böni, R. Georgii, A. Schmehl, D. G. Schlom, and J. Mannhart, “Field and temperature dependence of the magnetization in ferromagnetic EuO thin films,” Journal of Physics Condensed Matter, vol. 20, no. 10, Article ID 104230, 2008. View at: Publisher Site  Google Scholar
 W. Söllinger, W. Heiss, R. T. Lechner et al., “Exchange interactions in europium monochalcogenide magnetic semiconductors and their dependence on hydrostatic strain,” Physical Review B—Condensed Matter and Materials Physics, vol. 81, no. 15, Article ID 155213, 2010. View at: Publisher Site  Google Scholar
 P. Liu, J. A. C. Santana, Q. Dai, X. Wang, P. A. Dowben, and J. Tang, “Sign of the superexchange coupling between nextnearest neighbors in EuO,” Physical Review B—Condensed Matter and Materials Physics, vol. 86, no. 22, Article ID 224408, 2012. View at: Publisher Site  Google Scholar
 R. F. Wallis, A. A. Maradudin, I. P. Ipatova, and A. A. Klochikhin, “Surface spin waves,” Solid State Communications, vol. 5, no. 1, pp. 89–92, 1967. View at: Publisher Site  Google Scholar
 H. Puszkarski and S. Mamica, “Interfacelocalized mode in bilayer film ferromagnetic resonance spectrum,” Horizons in World Physics, vol. 244, pp. 1–43, 2004. View at: Google Scholar
 C. Kittel, Introduction to Solid State Physics, Wiley, 2004.
 S. Mamica, “Effects of spinwave energy degeneracy induced by inplane wave propagation in ferromagnetic thin films,” Journal of Magnetism and Magnetic Materials, vol. 322, no. 20, pp. 2994–2999, 2010. View at: Publisher Site  Google Scholar
 S. Mamica and H. Puszkarski, “Nearest and nextnearest neighbourhood surface mapping in thin cubic films. Properties of structural sums,” Acta Physicae Superficierum, vol. 9, no. 9, pp. 17–70, 2006. View at: Google Scholar
Copyright
Copyright © 2015 Sławomir Mamica. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.