Magnonic Band Structure, Complete Bandgap and Collective Spin Wave Excitation in Nanoscale Two–Dimensional Magnonic Crystals

Magnonic Band Structure, Complete Bandgap and Collective Spin Wave Excitation in Nanoscale Two–Dimensional Magnonic Crystals

D. Kumar Thematic Unit of Excellence on Nanodevice Technology, Department of Condensed Matter Physics and Material Sciences, S. N. Bose National Centre for Basic Sciences, Block JD, Sector III, Salt Lake, Kolkata 700 098, India.    J. W. Kłos Faculty of Physics, Adam Mickiewicz University in Poznan, Umultowska 85, Poznań, 61-614, Poland.    M. Krawczyk Faculty of Physics, Adam Mickiewicz University in Poznan, Umultowska 85, Poznań, 61-614, Poland.    A. Barman Thematic Unit of Excellence on Nanodevice Technology, Department of Condensed Matter Physics and Material Sciences, S. N. Bose National Centre for Basic Sciences, Block JD, Sector III, Salt Lake, Kolkata 700 098, India.
August 28, 2019

We present the observation of a complete bandgap and collective spin wave excitation in two–-dimensional magnonic crystals comprised of arrays of nanoscale antidots and nanodots, respectively. Considering that the frequencies dealt with here fall in the microwave band, these findings can be used for the development of suitable magnonic metamaterials and spin wave based signal processing. We also present the application of a numerical procedure, to compute the dispersion relations of spin waves for any high symmetry direction in the first Brillouin zone. The results obtained from this procedure has been reproduced and verified by the well established plane wave method for an antidot lattice, when magnetization dynamics at antidot boundaries is pinned. The micromagnetic simulation based method can also be used to obtain iso–frequency countours of spin waves. Iso–frequency contours are analougous of the Fermi surfaces and hence, they have the potential to radicalize our understanding of spin wave dynamics. The physical origin of bands, partial and full magnonic bandgaps has been explained by plotting the spatial distribution of spin wave energy spectral density. Although, unfettered by rigid assumptions and approximations, which afflict most analytical methods used in the study of spin wave dynamics, micromagnetic simulations tend to be computationally demanding. Thus, the observation of collective spin wave excitation in the case of nanodot arrays, which can obviate the need to perform simulations may also prove to be valuable.

preprint: thanks: These authors have contributed equally to this work.thanks: These authors have contributed equally to this work.

I Introduction

Analogous to photonic crystals, magnonic crystals (MCs) Nikitov, Tailhades, and Tsai (2001); *Demokritov2013; *Kruglyak2010 are magnetic meta-materials designed for the propagation of spin waves (SWs).Hillebrands and Ounadjela (2002); Stancil and Prabhakar (2009); Neusser and Grundler (2009); *Lenk2011 Based on their design, MCs exhibit a characteristic SW dispersion relation complete with bands and, sometimes, band gaps which can be tuned by controlling material and structural parameters as well as the strength and orientation of the bias magnetic field.Mamica, Krawczyk, and Kłos (2012); Yang, Yun, and Cao (2011) This phenomenon makes MCs useful as potential candidates for the design of SW based signal processing and logic devices.Ding, Kostylev, and Adeyeye (2012)

The knowledge of dispersion relation of a wave propagating through a medium is necessary to understand its transmission characteristics. Although MCs have been a subject of intense study lately,Mamica et al. (2012); Au et al. (2012); Kłos et al. (2012); Mandal et al. (2012); *Ciubotaru2012; *Tacchi2012; *Gubbiotti2012; *Saha2013; Venkat et al. (2013); Kumar et al. (2013) reports on a time domain numerical calculation of dispersion relations of SWs propagating in two-dimensional (2D) MCs are very rare.Kim (2010); Schwarze et al. (2012); Dvornik, Au, and Kruglyak (2013) As other analytical methods are available, the use of time domain simulations and spatial Fourier transform to obtain the dispersion relation in a photonic or phononic crystal is rarely seenJafarpour et al. (2004) as well. We hope to fill that gap in research with this work. The underlying principles, over which the procedure described here is used, has been discussed more generally by us in Ref. Kumar et al., 2012. Here too, we essentially use a micromagnetic simulator called Object Oriented Micromagnetic FrameworkDonahue and Porter (2002) (OOMMF) to obtain magnetization , as a function of position , and time . Then we use a multi-domain discrete Fourier transform to obtain the desired dispersion relation: SW power as a function of wavevector , and frequency . However, while simulating the magnetization dynamics in large (ideally infiniteKruglyak et al. (2010)) 2D crystals, one can be expected to need far greater computational resources than during the simulations of their one-dimensional (1D) counterparts.Kłos et al. (2012) Using a finite sample size may produce some spurious modes in the obtained dispersion relation.Kruglyak et al. (2010) Thus, the use of 2D periodic boundary conditionWang et al. (2010a) (PBC) becomes mandatory in order to obtain good numerical resolution in wavevector and frequency domains while consuming finite computational resources. Also, 2D crystals have more high symmetry directions when compared to their 1D analogues. Different techniques will be required to obtain the results for different directions in the 2D reciprocal space covering the entire irreducible part of the Brillouin zone (BZ).Kittel (1976) Moreover, the signal which generates the waves will have to be carefully designed so that the resulting spectrum represents the physical dispersion relation of plane propagating SWs. Due to all these complications, a need to validate the results obtained here with a well established method, such as the plane wave method (PWM)Krawczyk and Puszkarski (2008) becomes very clear.

Both OOMMF and PWM solve the following Landau-Lifshitz-Gilbert (LLG) equation:Landau and Lifshitz (1935); *Gilbert2004


Here, is the saturation magnetization, is the gyromagnetic ratio and is the Gilbert damping constant. is the sum of the bias field , the excitation signal (used in micromagnetic simulations (MS) only), the exchange field, the demagnetizing field and the magneto-crystalline anisotropy. Other factors like magnetostriction, should also be considered here if applicable.

The details of MC considered here are presented in Sec. II.1. Simulation parameters and PWM are described further in Sec. II.2. OOMMF uses the finite difference method (FDM) to solve Eq. (1) as an ordinary differential equation in time and space (derivatives with respect to space are hidden away in ). PWM is based on the Bloch wave formalism. As these two methods are fundamentally different in approach, some quantitative differences in results are to be expected. The results from both the methods and their differences have been discussed in Sec. III for the antidot lattice (ADL). Due to small lattice constant, the considered system is an exchange dominated one and consequently, the differences in dispersion relations along the bias magnetic field and perpendicular to it are subtle. These differences have been explored by calculating the iso-frequency contours in the wavevector space using both MS and the PWM. The iso-frequency contours are the curves of the constant frequency plotted in the wavevector space, they are wave counterparts of the Fermi surfaces known from the theory of the solid state physics.Kittel (1976) The iso-frequency contours are very important tool for the analysis of the wave propagation phenomena, giving a deep insight into direction and velocity of propagating, reflected and refracted waves in artificial crystals. Such type of analysis, while widely explored in photonic and phononic crystals for designing their metamaterials properties,Luo et al. (2002); Enoch, Tayeb, and Gralak (2003); Li et al. (2011) is almost absent in magnonics. Thus, developing the ability to compute these iso-frequency contours using MS can be a breakthrough in exploring magnonic metamaterials based on MC; because the MSs can be performed without approximation limited applicability of the PWM (or other analytical methodsGrimsditch et al. (2004)), and thus yields experimentally realizable results even with complex magnetic configurations.

We also plot the energy spectral density and phase distributions associated with different modes in the SW spectrum in order to understand their physical origin and explain any observed partial or complete bandgaps. Finally, we use the method described here to obtain the SW dispersion relations in the case of 2D dot array where the SW propagation is mediated by inter-dot stray magnetic field as opposed to dipole-exchange interaction in ADL. This brings about an interesting change in the spectra, which is discussed in Sec. III along with their effective properties.

Ii Method

ii.1 Magnonic crystal lattice and material parameters

The structure considered here is an infinitely large square array of square antidots with their ferromagnet-air interface under pinned boundary conditions.Kłos et al. (2012) The geometrical structure of the sample is shown in Fig. 1 (a). The lattice constant nm and the antidots are square holes of edge length, nm. The material parameters of permalloy (Py: NiFe) are used during simulations and in PWM calculations: exchange constant, J/m, saturation magnetization, A/m, gyromagnetic ratio, m/As and no magnetocrystalline anisotropy. A saturating bias magnetic field of T points in direction.

Figure 1: (a) The 2D antidot lattice under consideration. The lattice constant of the square lattice is nm. The thickness of the film is nm. The antidots are square (white) air holes of edge nm in ferromagnetic Py (black) medium. Dynamics is pinned at the edge of holes. The pinned region is marked with a different texture. Element geometry used in micromagnetic simulations extends to over hundred repetitions in length (horizontal dark arrows in (b), (c) and (d)) for good wavenumber resolution. 2D PBC is applied over these elements to mimic the infinite geometry. White arrows in (b), (c) and (d) show the direction of bias field used for simulations of SW dispersion for backward volume and Damon-Eshbach configuration. (d) shows the first BZ in the inverse lattice with typical symmetry point labels.

ii.2 Micromagnetic simulations and the plane wave method

The micromagnetic simulations involve solving the LLG equation (1) using a finite difference method based ordinary differential equation solver; and then, Fourier transforming the obtained space and time dependent magnetization data to get SW spectral density in wavevector and frequency domains.Kumar et al. (2012) Cell size nm along () axis was used during the FDM based simulations. The pinning in micromagnetic simulations was introduced by fixing magnetization vector in all cells of the discretization mesh, which border the antidots, i.e., in regions marked with different texture in Fig. 1. Figures 1 (b), (c) and (d) show parts of the elements over which 2D PBC are used to simulate the dispersion relation for different directions of the wave vector. These elements extend over (up to ) repetitions of unit cells in the horizontal direction to yield good resolution in the wavenumber domain. The 2D PBC are also implemented in order to improve the results with finite computational resources.Wang et al. (2010a) Figure 1 (e) shows the first BZ, the path in its irreducible part and typical symmetry points: , , and .Kittel (1976) Note that when the bias field is in the plane, an asymmetry is expected between the two orthogonal directions of SW propagation: (backward volume) and (Damon-Eshbach).Venkat et al. (2013) Thus, the triangle is no longer the irreducible BZ. However, in the forward volume arrangement when is perpendicular to the plane of the 2D MC, the symmetry is restored and dispersion is the same in the two orthogonal directions. The technique described here can be used independent of the direction of .

In order to get the results in the and directions, we use different excitation signals of the form , on elements shown in Figs. 1 (b) and (c), respectively. is horizontal along the axis (dashed white arrows). Similarly, dispersion along the and directions can be obtained when is across the width of the elements (vertical arrows along axis). Here, with mT and , and as given by Eqs. (2), (3) and (4), respectively:


See Eq. () in Ref. Kumar et al., 2013 for the detailed description of the terms involved in these equations. Here, the origin of coordinates is at the center of the considered geometry. It is due to and that the signal contains power between and in frequency and wavevector domains respectively.Kumar et al. (2012) should be asymmetric to ensure that both symmetric and antisymmetric modes are present in the resulting spectrum.Hillebrands and Ounadjela (2002) In Eq. (4), goes from to . While computing dispersion along and directions (Fig. 1 (b)), . However, for and directions (Fig. 1 (c)), . Both the elements in Figs. 1 (b) and (c) will span the same infinite 2D geometry under a 2D PBC; except, in the later case we can control whether the dynamics in the neighboring rows will be in phase or out of phase. Thus we can fix the wavevector component or to 0 or in the simulations. This is necessary to differentiate between the parallel directions and or and . Also, given by the expression


can be used instead of to selectively alter the amplitude of -th symmetric or -th antisymmetric mode. The freedom of choice of amplitudes and allows us to artificially control the statistical temperature of the magnons in the crystal and also helps in isolating a single mode in the case of a degeneracy. We can also sum over and to alter multiple modes in a single dynamic simulation. We also attempt to obtain the dispersion in direction by using the element shown in Fig. 1 (d). However, as there are two scattering centers (antidots) per cell in this arrangement, we can obtain the dispersion relations correctly only up to half of the BZ in that direction.Schwarze et al. (2012)

Until now, we could use a signal similar to the one we did in the case of an 1D lattice.Kumar et al. (2013) But, this limitation forced us to come up with a new signal


which has to be used in a larger 2D lattice of antidot array (with the cell size increased to nm to decrease time of computations). Here, is given by:


with analogous formula for . This signal is a point like source with the amplitude decay with distance as described by sinc function (in Eq. (3) along axis and in similar form for for dependence), having sharp cut-off in Fourier domain and able to excite symmetric and antisymmetric modes with respect to or axis. This signal was arrived upon largely by intuition, nevertheless, its agreement with the results obtained from PWM validates the usefulness of this signal. Spectral density, periodicity and asymmetry of the excitation signal (or source) should also be considered while developing similar techniques for other kinds of crystals (e.g. photonic or phononic crystals).

Three fold (one in time and two in space) Fourier transforms was needed to obtain the SW dispersion here. Magnetization was assumed to be uniform across the thickness of the film. We can now easily generalize that in the case of three-dimensional MCs, a signal of the form will be required followed by a four fold discrete Fourier transform.

We have also calculated the spatial distribution of energy spectral density (ESD), and phase, from the following equations:


Here, is the time domain Fourier transform, of a dynamical magnetization data. Unlike the new method used in Ref. Kumar et al., 2013, this gives us power from the entire wavevector domain for a selected frequency . However, if power is present for just one particular wavevector then both methods yield qualitatively identical results.

The PWM is a spectral method in which the eigenproblem is numerically solved in the frequency and wavevector domains by the standard numerical routines. We solve here LLG equation (1) in linear approximation without damping. The PWM calculations are performed with the assumption of the full magnetic saturation of the ADL along the bias magnetic field. As pinning during simulation will occur at the cell’s center, a hole size of was assumed during PWM calculations. Due to small thickness of the ADL, uniform SW profile across the thickness is assumed. The PWM in this formulation was already used in the calculations of the SW dynamics in 2D ADL and proved to give correct results.Neusser et al. (2011a); Tacchi et al. (2012b); Kłos et al. (2012, 2013) The detailed description of the method can be found in Refs. Tacchi et al., 2012b and Neusser et al., 2011b.

Iii Results and Discussions

Figure 2: (a) SW dispersion calculated using MSs (solid line) and PWM (dashed lines). ESD , distribution for the horizontal line ( GHz) shown in (a) in parts of the sample when the propagation direction is along for (b) and (c) . Corresponding phase , distribution is shown in (color online) (d) and (e), respectively.

The dispersion along the path in the first BZ shown in Fig. 1 (e) calculated with MSs by using the elements shown in Fig. 1 (b)-(d) is assembled as Fig. 2 (a) using solid lines. An overlay of dashed lines representing the SW dispersion relation obtained from the PWM is provided for comparison. Both these results appear to agree with each other except for the direction where the numerical method was able to yield results for only half of the total BZ extent. This is because we set to here (the spatial periodicity is ). Compared to the element shown in Fig. 1 (b), which can be used to produce results for or directions, the one in Fig. 1 (d) features two scatter centers per unit cell. And, if we artificially increase to , both scattering centers will be activated to produce additional spurious modes.Schwarze et al. (2012) To demonstrate the same we plot (normalized between and ) and (given by Eqs. (8) and (9), respectively), for frequency GHz in Figs. 2 (b) to (e). Note that the horizontal separation between regions of high ESD is about in Fig. 2 (b) for . This reduces to in Fig. 2 (c) for when both scattering centers in the unit cell (of the element shown in Fig. 1 (d)) are activated at once. The phase distributions also confirm that neighbouring locations of high ESD are about and radians out of phase with each other in former (: Fig. 2 (d)) and later (: Fig. 2 (e)) cases, respectively. Apart from incomplete result for the direction, we can also see that the modes here (shown by solid lines) do not match with those for direction at the point. This is because (cell size) nm was used while simulating for the direction as opposed to nm, which was used in the case of direction. Also, there are additional modes of lower amplitudes visible in the case of direction. This is due to the fact that becomes a stepped approximation of the right hand side of (3) by the use of the FDM; thus compromising the effectiveness of the cut off at , and exciting the second scattering centre to some extent (but not as well as ).

Figure 3: (a) SW dispersion calculated using MS (solid line) and the PWM (dashed lines). The full and partial magnonic bandgaps are marked and numbered by Roman numerals. The circled Arabic numerals indicate the points on the dispersion for which the mode profiles are calculated in Fig. 4. Iso-frequency lines from (b) GHz to GHz (c) GHz to GHz (it is around the top and bottom of the first and second magnonic band, respectively) using the PWM is shown with dashed lines. Iso-frequency lines for (b) GHz and (c) GHz calculated by the numerical method is shown using solid lines.

In pursuit of our quest to close the gap in the direction we eventually decided to simulate the SW dynamics in a large 2D MC with signal defined by Eq. (6) and perform a three-fold Fourier transform in contrast with the two-fold transforms done earlier. We transformed time to frequency domain and and dimensions to the 2D wavevector domain. The resulting dispersion relation as calculated from numerical method is shown in Fig. 3 (a) using solid lines. The agreement with the PWM results (shown by dashed lines) is poorer in comparison with Fig. 2 (a). This is due to the fact that cell size in the later attempt was increased from nm to nm. The complete and partial bandgaps width and center frequency, as seen from the dashed lines in Fig. 3 (a), are extracted in Tab. 1. Here, values for partial bandgaps depend upon the path, which has been used to plot the dispersion. Bandgap I is the only complete bandgap observed here with the maximum width of GHz.

Label Extent Center (GHz) Gap Width
I Complete Bandgap 76.39 15.36
II 100.18 9.24
III 114.45 7.7
IV 155.85 1.9
V 107.75 5.9
VI 149.6 8.4
VII 100.68 10.24
VIII 114.2 8.2
IX 85.80 3.01
X 108.1 4.6
XI 150.1 9.4
Table 1: Magnonic bandgap widths and center frequencies across different directions as calculated by the PWM and labeled in Fig. 3 (a).

Most bands observed in Fig. 3 (a) increase or decrease almost monotonously along any high symmetry direction. Consequently, the width of bandgap I too, appears to decrease monotonously as we move either along or . Both upper and lower limits of bandgap I are present at point which suggests an anti-crossing of bands at that point. This can also be regarded as the cause of the gap formation. Narrower bandgap widths have been observed by different techniques before.Wang et al. (2010b) The relatively high width of GHz of bandgap I here can be attributed to small lattice dimensions and edge pinning.Kłos et al. (2012) Bandgaps II to XI (Fig. 3 (a) and Tab. 1) are direction dependent partial bandgaps. This is mainly because bands approaching point from other high symmetry directions (with the exception of the band starting at 5⃝) tend to show greater slopes. As 5⃝5⃝ is a relatively flatter line, bandgap IV survives for three high symmetry directions. In a more isotropic forward volume arrangement,Schwarze et al. (2012) bandgap IV might also have qualified as a complete bandgap if the dispersion in the direction was also calculated. On the other hand, if wavevector dependent anisotropy is overlooked,Kim (2010) partial bandgaps (e.g. bandgap IV, or II and VII, or III and VIII) will appear as a complete bandgap. Partial bandgaps IV, V, VIII, X and XI are direct, while II, III, VI, VII and IX are indirect. Direct bandgaps are formed when the minimal and the maximal frequency of the magnonic bands surrounded the bandgap, from the top and bottom, respectively, are characterized by the same wavevector. While two different wavevectors are involved in the formation of indirect bandgap. In Fig. 3 (a) the minimal and maximal frequencies appear at high symmetry points. Occasionally, a bandgap may form between two high symmetry points due to anti-crossing of modes in a folded BZ,Kłos et al. (2013) but that is not observed here.

Figure 4: (Color online) ESD , and phase , for high symmetry points , , and at points 1⃝ through 5⃝ marked on Fig. 3 (a).

Now we calculate mode profiles ESD and phase , at the high symmetry points, using the PWM, for the first five modes as marked in Fig. 3 (a). The results are tabulated as Fig. 4 where is represented by color saturation and is represented by hue. A general trend of higher frequency mode profiles limiting themselves to smaller regions in real space is observed. This trend has been seen for 1D systems as well.Kłos et al. (2013) Here, mode profiles appear similar in size at points 5⃝, 5⃝ and 5⃝. Although, the distribution at 5⃝ is vastly different due to a (nearby) mode–crossing in the direction (see Fig. 3 (a)). Mode profile at may be obtained by rotating the mode profiles at by . Modes with negligible group velocity are trapped and forbidden to move in specific high symmetry directions. Also, the number of nodal lines, which controls the spatial quantization of modes, generally increases with mode number i⃝. No nodal lines are evident for 1⃝. Vertical and horizontal nodal lines are seen at 1⃝ and 1⃝, respectively; while 1⃝ features both vertical and horizontal nodal lines. From Fig. 3 (a), we can see that points 2⃝, 2⃝, 3⃝ and 3⃝ belong to the same mode and points 3⃝, 3⃝, 2⃝ and 2⃝ belong to a different mode. As the crossing between these modes occurs along the direction, the mode profiles at 2⃝ and 3⃝ are comparable. Similarly, mode profiles at 3⃝ and 2⃝ are also comparable except, 3⃝ has higher frequency and consequently, is more confined is space. In general, vertical and horizontal nodal lines dominate at points and , respectively; while a more isotropic distribution is observed at point and . Modes 1⃝, 2⃝, and 5⃝ are isotropic along and axes for the point. However, modes 3⃝ and 4⃝ are disposed along rows and columns, respectively. There local shape and size is comparable and accordingly, they are also degenerate as seen in Fig. 3 (a). Going from to either or , 4⃝ maintains its size and frequency; except the Damon-EshbachLeite and Filho (2001) geometry is evident in the later case. Similarly, the expanses of mode profiles at M2⃝ and 3⃝ are comparable (as there frequencies are within GHz of each other), and yet their orientations are mutually orthogonal.

Iso-frequency lines are shown in Figs. 3 (b) and (c), using both the PWM (dashed lines) and the MSs (solid lines). Iso-frequency contours calculated using the proposed method are thicker because small yields a low wavevector resolution. The agreement between the results obtained from the two methods as GHz line calculated using the MSs and the GHz line calculated using the PWM is clear, but the 2 GHz difference in frequencies is due to the shift of the dispersion curves calculated with both methods shown in Fig. 3 (a). In contrast to Fig. 3 (b) the two methods appear to give identical results for the GHz iso-frequency line, where the results of MS and PWM coincide. The shapes of iso-frequency lines control the direction of the propagating waves and consequently also alter the shapes of their wavefronts. Thus, although the dispersion along and directions may appear comparable, the wavefronts of the propagating SWs from the first band will quickly uncover the underlying anisotropy, because of slightly different group velocity and curvature of different iso-frequency contours in two orthogonal directions, which is easily noticeable by the inspection of the contours for 63 and 65 GHz in Fig. 3 (b). This anisotropy is a manifestation of dipolar interactions hardly visible in this size and frequency regime in the magnonic band structure shown in Fig. 3 (a). Backward volume modes are characterized by negative group velocity in the case of dipole dominated or dipolar-exchange SW propagating in a ferromagnetic thin film.Stancil and Prabhakar (2009) This is not seen in Fig. 3 (a) as due to weakness of the dipolar interactions the exchange field makes a significant contribution with increasing wavevector already near the BZ center.

Figure 5: First mode in a permalloy nano–dot array with varying angle , between the bias field ( = 1 T), and wavevector , showing the transition from magnetostatic backward volume mode to Damon-Eshbach configuration. The dashed lines are calculated using the analytic expressions for these two configuration with a reduced saturation magnetization. The structure considered here is given in the top left corner with nm, nm and thickness nm. Material parameters remain the same as before.

The developed method is not limited to the antidot lattices nor exchange dominated SWs. To prove this and better understand the properties of dipolar waves in MCs we take a look at the dispersion of SWs in the case of 2D MC composed of a square array (of lattice constant nm) of square dots (of edge nm and nm thick). This structure is shown in the top left panel of Fig. 5 along with the dispersion relations of the first mode with increasing angle (from to ), between  and in the subsequent panels. Here the wave propagation is mediated by the dipolar field only. We note how the mode’s group velocity gradually increases from negative (backward volume) to positive (Damon-Eshbach) as goes from to .Dvornik, Au, and Kruglyak (2013) The transition appears to occur at a critical angle . Note that here the direction of  is being changed as opposed to that of in the previous case. It is interesting to note that the dispersion relations obtained here for the array of nano–dots reminds us of the dispersion of magnetostatic waves in thin ferromagnetic film. To verify this hypothesis we calculate the dispersion relation of magnetostatic waves in the thin ferromagnetic magnetic film (3 nm thick) with reduced magnetization, i.e. with the effective value of the saturation magnetization . The dashed lines overlaid in Fig. 5 are computed using the analytical expression for backward volume and Damon-Eshbach configuration in the case of thin filmStancil and Prabhakar (2009) with a reduced saturation magnetization . A good agreement between the dispersion in the array and the effective thin film is found. A minor disagreement is introduced by the presence of the BZ boundaries but only near these boundaries. Further, the critical angle,  Hurben and Patton (1995) in the case of such thin film is also . This implies that one should also be able to use the analytical expression to calculate the SW manifold between backward volume and Damon-Eshbach geometries. This also shows, that a thin film MC composed of an array of saturated ferromagnetic nanodots can be used as a magnonic metamaterial, i.e., an artificial crystal with tailored effective properties of spin wave dynamics.Kruglyak et al. (2012); Zivieri and Giovannini (2012); Mruczkiewicz et al. (2012); Dmytriiev et al. (2012) Further studies are necessary to elucidate the limits of the effective saturation magnetization approach presented here, the influence of the dot–shape, their arrangement and inter–dot separation (mode-splitting has been experimentally demonstrated for nano–dot arraysRana et al. (2011)), but this is outside the scope of this paper.

Iv Conclusions

We have described a numerical algorithm to calculate the dispersion of plane propagating SWs in a 2D MC using multi-domain Fourier transform of results obtained from micromagnetic simulations. At the core of this technique is a new excitation signal, which is capable of generating SWs whose energy spectral density corresponds to the characteristic dispersion relation of the 2D MC. The lack of such signal has been discussed before in the case of 1D MCs.Di et al. (2013); *Lee2013 The results obtained from this procedure were verified by the plane wave method when magnetization dynamics at antidot boundaries is pinned. We noted that both methods were in qualitative agreement with each other. The fact that better quantitative agreement was observed while using 2D PBC over 1D elements was due to lower cell size.Kłos et al. (2013)

Apart from a new numerical algorithm to compute the dispersion relation in any given direction of a two- or three-dimensional inverse lattice, this method will also allow for the numerical computation of iso-frequency contours from micromagnetic simulations. Thus the numerical tool for study metamaterials properties of MCs was provided. It gives the possibility for design the properties of SWs relevant to technological applications and potentially exceeding these known from the homogeneous ferromagnetic thin films. The negative refraction, unidirectional media or caustic propagations are only some of examples here.Veerakumar and Camley (2006); Kostylev, Serga, and Hillebrands (2011); Kim et al. (2008) Further, this method can be generalized to aid the numerical computation of dispersion or iso-frequency contours in the case of two- or three-dimensional phononicAlleyne and Cawley (1991) and photonicBurt et al. (2004); *Gersen2005 crystals as well.

The dispersion here appeared to be similar in and directions. However, a noticeable anisotropy between the backward volume and Damon-Eshbach geometries was very evident from the study of the mode profiles and the iso-frequency contours. As dipole field mediates the SW propagation in a 2D dot array we were able to obtain the negative group velocity associated with the first mode in the case of a backward volume magnetostatic configuration. We were also able to analyze the nature of bands and complete and partial bandgaps that were obtained from the dispersion calculations in the case of an MC. This can be useful in the design of attenuators,Sekiguchi et al. (2012) phase-shifters,Au et al. (2012) filtersKim, Lee, and Han (2009) and logic gates.Khitun, Bao, and Wang (2008)

Recent advances in lithography techniquesMandal et al. (2012); Rahman et al. (2012); Sidorkin et al. (2009); Siegfried et al. (2011) have made it possible to fabricate dot and antidot lattices with a resolution below 10 nm. Thus, one can fabricate samples with dimensions comparable to the systems considered here. Experimental techniques similar to Brillouin light scattering spectroscopyPerzlmaier et al. (2005) can be used to explore the SW dispersion relation.

We acknowledge the financial support from the Department of Science and Technology, Government of India (Grant nos. INT/EC/CMS (24/233552), Department of Information Technology, Government of India (Grant no. 1(7)/2010/M&C), the European Community’s FP7/2007-2013 (GA nos. 233552 (DYNAMAG) and 228673 (MAGNONICS)) and NCN of Poland (DEC-2-12/07/E/ST3/00538). D.K. would like to acknowledge financial support from CSIR - Senior Research Fellowship (File ID: 09/575/(0090)/2011 EMR-I) and fruitful discussions with A. Mookerjee at the S. N. Bose National Centre for Basic Sciences.


Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description