Curvature induced magnonic crystal in nanowires

A new type of magnonic crystals, curvature induced ones, is realized in ferromagnetic nanowires with periodically deformed shape. A magnon band structure of such crystal is fully determined by its curvature: the developed theory is well confirmed by simulations. An application to nanoscale spintronic devises with the geometrically tunable parameters is proposed, namely, to filter elements.


Introduction
Recent advance in fabrication of sub-micrometer sized magnetic nanowires of controlled shape [1][2][3] makes them promising construct elements for variety of spintronic devices. Curvilinear geometry opens new possibilities for the device miniaturization [4] and such curvilinear elements were suggested for using in computer memory and logic [5,6], in magnon waveguides [7]. That is why a clear understanding of physical picture of influence of the wire curvature on the magnetization dynamics is an actual problem of the modern micromagnetism.
Probably the most prominent manifestation of the curvature in magnetic systems is a magnetochiral effect [8]. In thin curved wires, the chiral symmetry breaking results in emergent geometry-induced Dzyaloshinskii-Moriya interaction (DMI), for a review see [1]. A number of the magnetochiral effects in dynamics of the nanowire magnetization were recently predicted. Spin waves can be bound by a local bending of the wire [9]. A wire twisting results in non-reciprocal spin-wave propagation [10]. Numerous curvature effects in the domain wall dynamics were found, namely the localized curvature defect (e.g. a bend of the wire) creates a pinning potential for the domain wall [11]; the curvature gradient acts as a driving force for the domain wall resulting in a fast translational motion free of the Walker limit [12]; torsion of the helix shaped wire can result in a negative mobility of the spin-torque driven domain wall [13], etc. These and others curvature induced phenomena may be explained within a general framework [14], which introduces a number of emergent geometry-induced interactions. The latter are effectively generated by the common energy terms comprising spatial derivatives, e.g. in presence of the curvature the exchange interaction generates the effective DMI and anisotropy [14,15], the intrinsic DMI generates effective anisotropy [16,17].
Here we consider the magnon spectrum of a wire with periodically deformed shape. We demonstrate that such geometry gives rise to the periodical potentials in equations for the magnon propagation, as a consequence, the band structure appears. The analysis is based on the previously developed approach [14] for curvilinear ferromagnetic wires. The full scale numerical analysis of the band structure is supplemented by the analytically obtained formulae valid for the limit cases of small and large curvatures. Additionally, we use the direct simulations of the Landau-Lifshitz equation in order (i) to verify the obtained results and (ii) to demonstrate that the periodically deformed wire can be utilized as a magnon filter. The latter is an important issue for the spintronic devices relied on the magnon crystals [7,[18][19][20][21][22][23][24].
A large area of modern spintronics deals with the data transfer and processing based on the magnon waves propagation [7]. In this case, magnonic crystals [18] are the key elements of the magnonic spintronic devices including resonators, generators, filters, wave-guides, see reviews Ref. [7,18,25] for more information. Magnonic crystal is a magnetic system with artificially introduced periodicity. The latter can be created by means of periodical arrangement of several different magnetic materials [18,23,26] or by using one material with periodically modulated geometrical parameters [18,24,27] including arrays of interacting magnetic strips [18,28,29] or dots [18,25]. The periodically arranged magnetization patterns can be also considered as a magnon crystal, e.g. the DMI induced helical structure [30,31] or skyrmion lattice [30,32].
Magnonic crystal, which we propose here, can be thought as a combination of the last two methods. On the one hand, we consider a curvilinear nanowire made of a single material. On the other hand, the periodically modulated curvature induces the periodical magnetization pattern, see Figs. 1, 2, resulting in periodical potentials. The magnon scattering on the periodical potentials induces the band structure. This process is completely analogous to formation of the energy zones in solid state crystals [33], but instead of the common Schrödinger equation we have its generalized form see Sec. 3 for details.
All analytical predictions are in good agreement with the spin lattice simulations, see App. C. In order to demonstrate that the meander-like periodical structure, see Fig. 1(a), can be used as a magnon crystal we consider a 14-periods part of such structure as a magnon filter. Using the spin-lattice simulations based on the Landau-Lifshitz equation we obtain the amplitude-frequency characteristic. The input signal was formed by a homogeneous rotation of the first spin with the input frequency. The output signal demonstrates the presence of the gaps in the system. Therefore the corresponding frequencies are effectively filtered. Position and width of the band gaps are determined by the curvature of the wire, see App. C.2, Fig. 1(d) and Supplemental movie 1 for details.

Model and the equilibrium magnetization pattern
A thin wire can be defined as a space domain r(s, ξ 1 , ξ 2 ) = γ(s) + ξ 1 e n (s) + ξ 2 e b (s). Here γ(s) is the central line of the wire parametrized by its arc length s, coordinates ξ 1 and ξ 2 parametrize the perpendicular cross-section of the wire, vectors e n and e b are the normal and binormal Frenet-Serret orts [34] of the line γ, respectively, see Fig. 1(a ). The term 'thin wire' means the simultaneous fulfilment of two conditions: geometrical one and magnetic one. The geometrical condition reads |∂ 2 s γ| 1/h max , where h max is the maximal lateral size of the perpendicular wire cross-section. The magnetic condition will be discussed latter. Together with the tangential unit vector e t = ∂ s γ the normal and binormal vectors compose the TNB basis {e t , e n , e b }. For the particular case of a planar wire lying within the xyplane the TNB basis vectors can be presented in the form e t (s) =x cos α(s) +ŷ sin α(s), e n (s) = −x sin α(s) +ŷ cos α(s), and e b =ẑ with alternating-sign curvature κ = ∂ s α.
Dynamics of the wire magnetization is governed by Landau-Lifshitz equation where m = M /M s is the unit magnetization vector normalized by the saturation magnetization, γ 0 > 0 is the gyromagnetic ratio, E is a total magnetic energy, and the damping is neglected. We model the magnetization subsystem by the Hamiltonian which takes into account only two magnetic interactions, namely: the isotropic exchange with the stiffness A (the first integrand) and the easy-tangential anisotropy with the coefficient K > Figure 1: Meander-like shape nanowire: Wire in form of periodically repeated semicircles (a),(a ) has the square-wave curvature κ(ξ) (b) with dimensionless amplitude κ 0 = /R and period 2ξ 0 = 2π/κ 0 . We demonstrate that the considered meander-like shape structure can be used as a magnon filter element using spin-lattice simulation based on Landau-Lifshitz equation with damping, for details see Appendix C.2. Colour scheme shows the binormal magnetization distribution in the wire (a) with curvature κ 0 = 1.2 and input signal of dimensionless frequency Ω = 1. Amplitude-frequency characteristics of input (c) and output (d) signals of a magnon filter element is composed of 14 periods of the meander shaped structure. The frequency ranges, where the magnons propagation is strongly depressed are shown by the grey shading. The band gaps position and width can be controlled by the parameter κ 0 . Decreasing the output amplitude with frequency takes place due to the natural damping (α = 0.01) included into the simulations. The terminator is placed after a magnon filter element to prevent signal reflection from the free boundary of the system. The terminator is composed of 5 periods of the meander shaped structure and has a spatially increasing damping coefficient, see Appendix C.2 for details. The magnetization dynamics are illustrated at the Supplemental movie 1 . 0 (the second integrand). Here S denotes the cross-section area, which is assumed constant along the wire. For the wires with the circular or square cross-sections the anisotropy constant reads K = πM 2 s +K a , where the first term comes from the magnetostatic contribution [35] and the second one represents the magnetocrystalline anisotropy. For the magnetically soft wires one has K a = 0. The model (1) determines the system and time-and length-scales presented by the frequency of the uniform ferromagnetic resonance for a straight wire ω 0 = 2Kγ 0 /M s , and magnetic length = A/K, respectively. The latter determines typical size of the magnetization nonuniformities, e.g. width of a domain wall. The model (1b) presumes that the magnetization is uniform within the wire cross-section: m = m(s, t). This assumption is valid for the case h max . The latter is the magnetic condition for the 'thin wire' approximation. In the following we use the dimensionless time τ = tω 0 and length ξ = s/ . This makes an analysis independent on the material parameters. We utilize the constraint |m| = 1 by means of the angular parametrization m = sin θ (cos φe t + sin φe n ) + cos θe b . For a planar wire the minimization of energy (1b) leads to the static equilibrium state θ = Θ = π/2 and φ = Φ(ξ), where the inclination angle Φ is determined by the driven pendulum equation [11,14] Φ − sin Φ cos Φ = −κ . ( Here and below prime denotes the derivative with respect to ξ, and κ = κ is the dimensionless curvature. The emergent geometry-induced DMI breaks the spatial symmetry of the solution: its contribution is determined by the gradient of the curvature in Eq. (2), the only geometrical parameter which influence the magnetization texture of any planar wire. In the following we consider a meander shaped wire in form of the periodically repeated semicircles of radii R, see Fig. 1. The spatial distribution of the curvature of such a wire is the square-wave function with period 2πR/ as shown in Figs. 1(c)-(d), see details in Appendix A.1. In this case the driving in Eq. (2) is determined by where κ 0 = /R is the curvature amplitude, ξ 0 = πR/ is half period of the curvature and δ(•) is the Dirac delta [36]. The corresponding solution of Eq. (2) reads where am(•, •) is Jacobi amplitude and • defines the integer part of • [36].
is elliptic integral of the first kind and the modulus k = k(Φ 0 ) is defined in (3b). Note that the curvature amplitude κ 0 is the only parameter which controls the system. The equilibrium solution (3b) is illustrated in Fig. 2. Magnetization of the equilibrium state always lies within the plane of the wire, but it is not tangential to it. The maximal deviation Φ 0 from the tangential direction takes place in points of junction of two semicircles, this is because of the curvature jump. In the limit case of small curvature κ 0 1 one has R and due to the easy-tangential anisotropy the wire is magnetized practically tangentially except the junction points, where the magnetization demonstrates the small deviations of amplitude Φ 0 ≈ κ 0 , see Fig. 2(b, c). In the intermediate case κ 0 ≈ 1 the curvature radius is close to the magnetic length (R ≈ ). This means the equilibrium state is determined by strong competition between the anisotropy and exchange interactions, see Fig. 2(d). The opposite limit case κ 0 1 corresponds to the regime R , when the dominating exchange interaction results in practically uniform magnetization aligned with x-axis. This corresponds to the maximal deviation amplitude Φ 0 π/2, see Fig. 2

Band structure
Here we consider linear dynamics of magnons excited over the equilibrium state described in the previous section. The technique which we use here was previously applied for a number of one-dimensional curvilinear magnetic problems: finding the magnon spectrum for the ring- [14] and helix-shaped [10] wires, localization and scattering of magnon waves on a localized curvilinear defect [9].
It is convenient to describe the linear magnetization dynamics by means of the complex valued function ψ = ϑ+iϕ constructed from deviations from the equilibrium state: θ = Θ+ϑ and φ = Φ + ϕ/ sin Θ. In this case, Landau-Lifshitz equation linearised in vicinity of the equilibrium state (Θ, Φ) has the form of generalized Schrödinger equation [9,37] −iψ = Hψ + W ψ * , for details see Appendix A. Here and below the overdot denotes the derivative with respect to the dimensionless time τ . In general case of an arbitrary planar wire the potentials V and  Table 1. The inset (c) shows detail of the band structure for small curvatures, the notations ate the same as in (b).
W have the form [9] where Φ is determined by Eq. (2). Application of Bogoliubov transformation ψ = u(ξ)e iΩτ + v * (ξ)e −iΩτ to Eq. (4) enables us to formulate an eigenvalue problem (EVP) for a Hamiltonian of Bogoliubov-de Gennes type Equation (6) is the general formulation of the spectral problem for a planar curvilinear wire described by the model (1). Generally, it can be solved numerically if the potentials V (ξ) and W (ξ) are known functions. Here we study a specific case, when the potentials (5) are periodic: Fig. 5. The periodic potential in a quantummechanical Schrödinger equation always produces the band structure. The band structure of the 1D system can be expressed quite simply in terms of the properties of scattering in the presence of a single-barrier potential [38]. Similar arguments can be used for the EVP (6), which results in the band structure with allowed frequency Ω(q) = K(q) 2 + v 2 − w 2 −1 and the quasi-wave vector q related by where v = max V (ξ) and w = min W (ξ). Using such an approach the band structure is determined by the scattering data for a single potential: the transmission coefficient T (K) and the scattering phase shift δ(K), see Appendix B for details. However even the scattering problem can be solved exactly only for few model potentials. Therefore below we also use another method.
Using the periodical properties of potentials, it is convenient to consider the EVP (6) in wave-vector space. For this case one can obtain the numerical solution for the band structure, see Appendix B.1, as well as analytical asymptotes for the limit cases of small and large curvatures, see Appendix B.2. The numerical solutions are shown in Fig. 3 by solid lines. In contrast to the straight wire, the dispersion relation for a given curvature κ 0 > 0 demonstrates the frequency reduction and appearance of the band structure, see Fig. 3(a). This is in a good agreement with the direct numerical simulation of the Landau-Lifshitz equation (1a). Evolution of the band structure with κ 0 is demonstrated in Figs. 3(b),(c). Generalizing the two-component approximation [33] for the case of the generalized Schrödinger equation (4), we approximate the band gap edges as Here Ω + ν and Ω − ν correspond to the top and bottom edges of the ν-th gap, respectively. Wavevector q ν = νκ 0 corresponds to the edge of the ν-th Brillouin zone. Coefficients V n , W n with n ∈ Z + denote Fourier components (B.6b) of the potentials. The approximation (8) is shown in Fig. 3(b) by thick dashed lines. The maximal relative deviation from the exact solution (solid lines) is 0.89%. Substituting the asymptotic approximations (B.14) and (B.17) for the Fourier coefficients into (8) we obtain the asymptotic behaviour of the gap edges for the limit cases of small and large curvatures, see Table 1. For the case κ 0 1 one has ∆Ω g ν ∝ κ 5 0 , where ∆Ω g ν = Ω + ν − Ω − ν is the width of the ν-th gap. In the opposite case all the gaps are closed by the law ∆Ω g ν≥2 ∝ κ −2 0 , except the first one, whose width goes to the constant value ∆Ω g 1 ≈ 3/4. Numerically found values of the first three band gaps are compared to the analytical asymptotes listed in Table 1 in Fig. 4(a). Interestingly that lim ν→∞ ∆Ω ν+1 /∆Ω ν = 1 for both cases of the vanishing and infinite curvatures, see Fig. 4(b).
The size of the main gap Ω 0 can be approximated by means of the one-component approximation The simple formula (9) demonstrates a good agreement with the numerical solution, see Fig. 3(b), therefore it is useful for practical estimations for the main gap width. However, it results in the wrong asymptotic for the limit case κ 0 1. Applying more exact threecomponent approximation (B.11) we obtain Ω 0 ≈ 1 − κ 2 0 /2 for κ 0 1 and Ω 0 ≈ 1/(2 √ 2κ 0 ) for κ 0 1. The asymptotic closing of the main gap for the large curvatures has the following physical sense. For large curvatures the curvature radius is much smaller than the magnetic length. In this case the wire magnetization is almost uniform and oriented along x-axis, see Fig. 2(c). On the other hand, the anisotropy axis frequently changes its direction within the xy-plane when moving along x-axis. Thus, in the long-wave approximation, such a wire is physically equivalent to the uniformly magnetized wire with easy-plane anisotropy, which is gapless [39].

Conclusion
We demonstrate that the periodical deformation of an uniaxial ferromagnetic wire results in appearance of the band structure in the magnon spectrum. As a case study we consider a meander-shaped planar wire composed of the semicircles of radius R. The normalized curvature κ 0 = /R is the only parameter which controls properties of the magnon spectrum. For small curvature κ 0 1 widths of all band gaps are ∆Ω g 1 ∝ κ 5 0 where ν ∈ N. In the opposite case κ 0 1 all the gaps closes by the law ∆Ω g ν≥2 ∝ κ −2 0 , except the first one, whose width goes to the constant value ∆Ω g 1 ≈ 3/4 (in units of the frequency of the uniform 2 Here we assume that Ω + 0 = Ω − 0 = Ω0. ferromagnetic resonance for a straight wire). The main gap shrinks as Ω 0 ≈ 1 − κ 2 0 /2 for small curvatures and it asymptotically closes Ω 0 ≈ 1/(2 √ 2κ 0 ) for large curvatures. In the latter case the curvilinear wire effectively behaves as an easy-plane magnet.
The curvature induced band structure is a direct physical evidence of the curvature effect on magnetic subsystem. Besides the fundamental meaning of this result the considered periodical system can be used in as a filter element in magnon spintronics applications. Here we consider a meander shaped wire with the periodically repeated semicircles of radii R, see Fig. 1. The curvature of this wire is the square-wave function with period 2ξ 0 :

A Equilibrium state and linearised equations of motion
For the considered geometry Eq. (2) is equivalent to the corresponding homogeneous equation supplemented by the derivative jumps in the junction points ξ = nξ 0 : where [Φ ] ξ = Φ (ξ + 0) − Φ (ξ − 0) and n ∈ Z. The opposite signs of the derivative jumps in the neighbouring junction points enable us to assume that the solution is a periodical function where the integration constant k is determined from the derivative jumps in (A.3). Note that Φ > 0 and Φ < 0 for the ranges where κ = −κ 0 and κ = κ 0 respectively. Thus, for all ξ one has Integration of Eq. (A.4) on the interval ξ ∈ [0, ξ 0 ] results in the solution The unknown constant Φ 0 can be found from the condition Φ(0) = −Φ 0 and Φ(ξ 0 ) = Φ 0 which is equivalent to the equation k(Φ 0 )F (Φ 0 , ik(Φ 0 )) = ξ 0 /2. Generalization of the half-period solution (A.6) for the all domain ξ ∈ R is presented by (3b). Asymptotic approximations for Φ 0 and Φ(ξ) for the cases of small and large curvatures are listed in Table. 1. Now we introduce the small deviations from the ground state as described in Sec. 3. The corresponding linearised Landau-Lifshitz equations have forṁ Introducing the complex valued function ψ = ϑ + iϕ we obtain the generalized Schrödinger equation (4) with potentials (5). Using (A.5) one presents the potentials in form where Φ is determined in (3b). Potentials (A.8) for different κ 0 are shown in Fig. 5. Note that V (ξ) = V (ξ + ξ 0 ) and W (ξ) = W (ξ + ξ 0 ) have twice as small period as Φ(ξ) has, this is because the dependence on ξ is reduced to the dependence on sin 2 Φ(ξ) in (A.8).
Looking for solution of (4) in form ψ = u(ξ)e iΩτ + v * (ξ)e −iΩτ one obtains the following set of homogeneous equations for functions u and v [9] which can be reformulated as EVP (6).

B Band structure for the Bogoliubov-de Gennes EVP in 1D
The general analysis of the band structure for the Bogoliubov-de Gennes EVP (6) can be performed, using an approach developed in Ref. [38] for the usual Scrödinger equation. We consider the periodic potentials V (ξ) and W (ξ) as a superposition of single-barrier potentials V and W of width ξ 0 : Let us start with the scattering problem for single-barrier potentials V(ξ) and W(ξ), which coincide with V (ξ) and W (ξ) within the interval ξ ∈ [−ξ 0 , ξ 0 ], respectively. Outside the interval we suppose, following [38], that potentials takes their limit values v = max V(ξ) and w = min W(ξ) when ξ / ∈ [−ξ 0 , ξ 0 ]. Then, general solution of eigenfunctions of the EVP (6) for single potentials in the range |ξ| > ξ 0 /2 has the following form, cf. [38] when ξ < − ξ 0 2 , Here we use the following notations: with transmission and reflection complex-valued coefficients being T = T e iδ and R = Re iδr , respectively, C j , j = 1, 4 being constants. Now using the Floquet theory, we conclude that u and v have to satisfy the following conditions: which are analogues of the Bloch's theorem for the EVP (6). By imposing conditions (B.4) and similar conditions for u and v at ξ = −ξ 0 /2 to (B.2), we find the band structure of EVP (6): cos (Kξ 0 + δ (K)) T (K) = cos qξ 0 .

B.1 Numerical band structure calculation in the wave-vector space
In the following it is convenient to utilize the periodicity of the potentials V and W and present them in form of the Fourier series where g 0 = 2π/ξ 0 = 2κ 0 is vector of the reciprocal lattice. Note that the potentials are even functions V (ξ) = V (−ξ), W (ξ) = W (−ξ). This has two consequences, namely V n = V −n , W n = W −n and V n = V * n , W n = W * n for all n. Substituting (B.6a) together with the representation u(ξ) = ∞ −∞û (q)e iqξ dq, v(ξ) = ∞ −∞v (q)e iqξ dq into (A.9) and using the orthogonality condition ∞ −∞ e iξ(q−q ) dξ = 2πδ(q−q ) we obtain the wave-vector representation of (A.9) For each q, equations (B.7) determine the tuple of numbers |Ψ q = | . . . ,û(q + ng 0 ), . . . ,v(q + ng 0 ), . . . which can be interpreted as eigenvector of the following eigenvalue problem HereĤ q is the block-matrix consisting of the blocksĤ q andŴ with elementŝ H nn q = δ n,n (q + ng 0 ) 2 + 1 + V |n−n | ,Ŵ nn = W |n−n | , (B.9) where n, n ∈ Z. Having the analytical expression (A.8) for the potentials V and W one can easily build matrixĤ q using (B.6b) for numerical calculation of the Fourier coefficientsV n and W n . The periodical band scheme for the first m zones is determined by m lowest eigenvalues of the matrixĤ q , see Fig. 6. In this case the size of the matrixĤ q must be much larger than m: −N 0 ≤ n, n ≤ N 0 , where N 0 m.

B.2 Analytical approximation
In vicinity of the ν-th Brillouin zone (ν ∈ N) the EVP (B.8) can be significantly simplified by applying the two-component approximation [33]. The solution is approximated by the standing wave with q = q ν , where the wave vector q ν = νg 0 /2 = νκ 0 corresponds to the edge Using the fact that matricesĤ qν andŴ commute, one obtains from (B.8) the following eigenvalues Ω qν = Ω ± ν .
Here Ω + ν and Ω − ν correspond to the upper and bottom edges of the ν-th gap, respectively, they are determined by Eq. (8). The value of the main gap Ω 0 , which corresponds to the case ν = 0 and q 0 = 0, can be roughly approximated by the one component approximation (8). Although, the latter works well in a wide range of κ 0 , it does not give the correct asymptotic for κ 0 1. For this case one should apply the three-component approximation: The straightforward calculation of eigenvalues of the matrixĤ q 0 and the subsequent substitution of the asymptotic approximations for the Fourier coefficients V 0,1 and W 0,1 listed in (B.14) and (B.17), enable us to estimate Ω 0 ≈ 1 − κ 2 0 /2 for κ 0 1 and Ω 0 ≈ 1/(2 √ 2κ 0 ) for κ 0 1.

C.1 Magnon dispersion
To verify our results, we numerically simulate the magnetization dynamics of the meander shape chain, Fig. 1(a) of the discrete magnetic moments m i , with i = 1, N . The magnetization dynamics is determined by the discrete Landau-Lifshitz equations with the relaxation term with ω 0 = 2Kγ 0 /M s , α being the damping parameter and E being the dimensionless energy normalized by K∆s 3 with ∆s being the lattice constant. The Hamiltonian of the magnetic wire reads with b i being a dimensionless external magnetic field, normalized by K/M s . The magnetization dynamics is described by a set of 3N ordinary differential equations (C.1) for unknown magnetization components m x i (t), m y i (t), m z i (t), i = 1, N . The set of equation (C.1) was integrated numerically.
We consider 20 periods of meander-like shape chain which corresponds to the length L = 2000∆s and the magnetic length l = 9.6∆s. During the simulation the discretization step is much smaller than the length of the single half-circle, ∆s π/κ 0 . The magnon dispersion relation is carried out in three steps. In the first step, we numerically simulate the magnetization dynamics for the system with Hamiltonian (C.2) and b i = 0 in overdamping regime, α = 0.1. The tangential distribution is chosen as initial one.
The obtained magnetization distribution corresponds to the ground state. Then the external magnetic field was applied along the shaped wire where q is the dimensionless wave vector, s i = (i − 1)∆s is the position of the magnetization vector m i and the amplitude of the external magnetic field b 0 = 0.005. In the last step, the magnetic field is switched off and the magnetization dynamics is simulated with damping coefficient α = 0.01. After that the space-time Fourier transform is performed for the binormal magnetization component. The frequency Ω corresponds to the maximum of the Fourier amplitude and the wave vector q is marked by a dot in Fig. 3.

C.2 Amplitude-frequency characteristic of the filter
We demonstrate that the magnon filter is a possible application for the proposed curvature induced magnonic crystal. We perform spin-lattice simulations of the magnetization dynamics, determined by Landau-Lifshitz equation (C.1). The meander shape chain of the discrete magnetic moments m i , with i = 1, N , lies in the plane x0y. The full energy of the magnetic system corresponds to the Hamiltonian (C.2) with b i = 0.
To form the input signal the first magnetic moment of the system is rotated with the input frequency Ω, m 1 =x cos Ωτ sin β +ŷ cos β +ẑ sin Ωτ sin β. We consider the 14 periods of the meander-like periodical structure with damping coefficient α = 0.01. To prevent signals reflection from the free boundary of the system the terminator is placed after the considered magnon crystal. The terminator is composed of 5 periods meander shape chain with N nodes and has the spatially inhomogeneous damping coefficient, α j = α exp (χs j ), j = N − N, N and the coefficient χ determines the rate of increase of the damping parameter.
The output signal is measured for the last node of the magnonic crystal. Due to non-linear processes, there is a transfer of energy density from the input frequency Ω to other spectral frequencies. To analyse the amplitude of the output signal on the input frequency the time Fourier transform is performed for one of the magnetization components.
Simulations are performed for curvature amplitudes κ 0 = 0.6 and κ 0 = 1.2, see Fig. 1(d). In the first case, the length of the periodic structure which corresponds to the magnonic crystal is 1400∆s and the first magnetic moment is rotating with β = π/2. The length of the terminator is 500∆s and χ = 0.01. In the second case, the length of the magnonic crystal is 700∆s and the first magnetic moment is rotating with β = π/12. The length of the terminator is 250∆s and χ = 0.025.