We investigate a special design of two-dimensional magnonic crystal, consisting of two superimposed lattices with different lattice constants, such that spin waves (SWs) can propagate either in one or the other sublattice, depending on which of the two frequency bands they belong to. The SW bands are separated by a very large bandgap (in our model system, 6 GHz), easily tunable by changing the direction of an applied magnetic field, and the overlap of their spatial distribution, for any frequency of their bands, is always negligible. These properties make the designed system an ideal test system for a magnonic dual band waveguide, where the simultaneous excitation and subsequent propagation of two independent SW signals are allowed, with no mutual interference.
1. Introduction
In recent years the field of magnonics has received considerable attention and has developed very fast: the idea of manipulating spin wave (SW) propagation, through its Bragg diffraction in artificial lattices in which magnetic properties are periodically varied (magnonic crystals, MCs), has been extensively investigated in numerous systems, both theoretically and experimentally, with numerous applications [1, 2]. Depending on geometry, material, magnetization, and applied field, collective SWs in MCs can be controlled at will for specific purposes and applications [3–8]. Versatile devices have been proposed, in which SWs can propagate at a tunable group speed, which can turn from waveguides to memories and vice-versa under the action of some external field [9–16] and logic gates where the output signal is determined by specific interactions between input propagating SWs [17–19]. The use of SWs as information carriers is particularly interesting in the context of limiting energy dissipation, with reference to ohmic losses, since no direct electron motion is involved [20, 21], and hence magnonic devices are considered key ingredients in the context of enhancing the device efficiency, even when integrated with the traditional ones.
In this paper we investigate a special two-dimensional (2D) MC, designed in such a way that SWs can propagate in distinct regions of it, depending on the frequency band they belong to. We characterize the two spin wave types, by determining the two localization areas and the extension of the corresponding bands.
2. Calculation Methods and System Design
We study the system within the micromagnetic framework, computing the equilibrium magnetization by the software OOMMF [22] and the spin wave dynamics by the dynamical matrix method (DMM) [23]. The equilibrium maps provided by OOMMF are the input of the DMM software, which calculates the normal mode frequencies and spatial profiles within a Hamiltonian framework, after computing a matrix, on which an eigenvalue/eigenvector problem is solved for any wavevector value and direction. Eigenvalues and eigenvectors provide spin wave frequency and spatial profile, respectively. We used the following magnetic parameters: saturation magnetization Ms=800 kA/m, exchange stiffness constant A=1.3×10-11 J/m, and gyromagnetic ratio γ = 185 rad GHz/T. An external magnetic field of constant magnitude B0=0.1 T is applied, at different directions, as shown in Figure 1.
Primitive cell of the 2D magnonic crystal used for calculations. The angle θ is the direction at which the external magnetic field B0 is applied. Distance between the small squares is a1=15 nm. This cell is unitary only for the secondary sublattice (a2=4a1). The red line square is the ideal unitary primitive cell of the primary sublattice (which was however not directly used in the calculations).
Due to the long-range nature of dipolar interactions, it is not trivial to find out a 2D periodic magnetic system, in which the magnetization dynamics involve two specific and distinct regions of it for all modes in two complete frequency bands, and, at the same time, the two bands are enough separated in frequency. To this purpose, the basic idea is using shape anisotropy to differentiate in energy the two sublattices within the primitive cell, and this is attained with the following MC design (Figure 1): a main (primary) square sublattice, with a short lattice constant (a1=15 nm), is made of cubes with a side of 5 nm, that is slightly lower than the exchange correlation length (λex=5.5 nm, in our case), and exactly corresponding to the micromagnetic elemental cell. A secondary square sublattice, with a long lattice constant (a2=60 nm), is made of larger prism particles, with identical thickness and width, but with a double length: 10×5×5 nm3 (hence consisting of two micromagnetic cells). In this way, the dipolar energy in smaller particles is larger than that in bigger particles, and this is a condition necessary to both increase the frequency of the modes localized in larger particles and determine a fast decay of the mode oscillation outside that region. Actually, the two lattices that we have chosen are not geometrically disjoint, but an element of the primary sublattice lays always over the left part of any element of the secondary sublattice. However, this is not altering the overall purposes and results, because the intersection area is always not involved by dynamic magnetization of the primary lattice in no cases: the two lattices can be seen as always dynamically disjoint.
We would like to comment our choice in the perspective of an experimental investigation. Apparently, due to the reduced dimensions of the nanoparticles, our MC system might look rather academic: as already remarked, the design of the present MC is intended to test the feasibility of dual propagation of independent signals within the same medium. However, the lowest limit of current technology (either by self-assembly techniques or by bioengineered crystallization) is a crystal dimension within 10÷5 nm, that is, rather close to our choice, 5 nm [24–26]. Moreover, in recent past other systems based on very small particles were investigated ([27] and references therein). At present, the investigated systems, usually composed of bigger elements, sometimes with different magnetic materials, allow spin waves with different localization features, which however belong to the same frequency interval: in these cases, when spatial separation is considered, frequency bands are not separated (i.e., dispersion curves of different modes cross and mix), and when a complete frequency gap is occurring between two bands, SW spatial separation is not occurring for all the frequencies of the corresponding bands, and hence there is no duality [4, 28]. For example, in a typical dot lattice [15], end mode bands are separated by a significant gap from extended wave bands; however, in the upper bands the dynamic magnetization involves even the region of the end modes: a frequency band separation does not strictly meet a SW spatial separation. The same can be said for other kinds of localized SWs in dot or antidot lattices and even bicomponent systems [4, 8, 29].
3. Results
The dynamic magnetization of each mode, resulting as output of DMM calculations, can be interpreted as a Bloch wave, using the following expression:(1)δmkr=δm~kreik·r,where r is the position vector in the direct space, k is the Bloch wavevector (k1 or k2, when relevant to the primary or secondary sublattice, correspondingly), and δm~k is the cell function, which has the periodicity of the array. In the following, we will consider only the case of uniform oscillation in the primitive cell (of either of the sublattices) and study the variation of frequency with k. Due to the fact that we used the primitive cell of Figure 1 for calculations, we could directly control k2 and only get an indirect evaluation of k1, at each frequency, after numerical inspection of the mode phase profiles, by extracting the effective wavelength concerning the primary sublattice elements and computing the corresponding wavevector value (k1).
The dual band property of the system allows the simultaneous propagation of two independent signals, of course, within the linear regime, which in our calculations was used to transform the Landau-Lifshts equation into a linear eigenvalue/eigenvector problem, while in real experiments it can be provided by a proper limitation of the excitation power [8]. The two signals could be excited by either ordinary or polarized current flowing on antennas or point contacts that are preliminary located on one or more elements of either of the lattices [30–33]. Another possibility is an electric excitation of a ferroelectric substrate under the MC, by means of the ferroelectric/ferromagnetic interaction [9, 34]. Recently, even the possibility of exciting a MC by means of optical/microwave frequencies has been explored, resulting in the need for defects in order to get a local excitation and opening the door to studies of defects in MCs [35–38]. In this way, we would like to contribute to the research demanding, from one side, ever smaller devices and, from the other side, ever higher density of information, either stored or delivered at frequencies within the GHz regime: in our model, these two apparently competing features meet in a simple, experimentally challenging answer.
In Figure 2, we plot the internal effective field Heff=B0/μ0+Hex+Hdip (i.e., the sum of the external field H0=B0/μ0 and the exchange Hex and dipolar Hdip=-|Hdip| internal fields) as a function of position along x-axis (with reference to Figure 1), for different directions of the applied field. From the plot, it is clear that the field direction, θ, critically affects Heff values only in the bigger elements (secondary sublattice), while variations of Heff in the smaller elements (primary sublattice) are rather small. Moreover, the mean values in the two sublattices are the farthest at θ=0° and the closest at θ=90°: actually, in the bigger elements, the dipolar field significantly increases as the external field is gradually rotated, because it forces the magnetization to progressively align to the minor axis of the particle (the secondary lattice element), suffering an increasing demagnetizing effect. Since spin wave frequency is critically dependent on the internal effective field, just from this observation it is possible to predict that both the spin waves will occur within either one sublattice or the other, and the corresponding mode frequency difference (magnonic gap) will be the largest at 0°.
Spatial profiles of the internal effective field, for different directions of the applied field of constant magnitude μ0H0=0.1 T (degrees, indicated in figure). Note that the internal field is zero in correspondence with the void spaces. At 0° the discrepancy between the region 0–10 nm and the other regions is the largest, while at 90° it is the minimum.
In Figure 3 we show the spatial maps of the fundamental modes [k=(0,0)] that we found: clearly, the dynamic magnetization is localized in either one sublattice or the other, and this property is rigorously hold even for any other wavevector value and direction. We note here that the profile of the fundamental mode relative to the primary sublattice (Figure 3(a)) has a nonuniform amplitude distribution, differently from the secondary one (Figure 3(b)): this is due to the fact that calculations are performed with the primitive cell of Figure 1 (superposition of both sublattices), which is in fact nonunitary for the primary sublattice and contains an element of the secondary sublattice (acting as a sort of “defect” for the primary sublattice, though periodic as well): this fact gives origin to a nonuniform distribution of the dipolar fields among the elements of the primary sublattice (squares) that are contained in the nonunitary primitive cell and, in turn, a nonuniform distribution of the dynamic magnetization of the corresponding quasi-uniform mode, shown in the figure. For this reason, when trying to extract from Figure 3(a) the “ideal” unitary cell function δm~kr relevant to the primary sublattice only, it is necessary to forget about the “real” intensity modulation and take the lowest periodic area (a1×a1) containing a single square element only (with reference to the red line square of Figure 1), while for the secondary sublattice the “ideal” cell function δm~kr can be extracted from Figure 3(b) considering the area (a2×a2) corresponding to the extension shown in Figure 1, with a single rectangle and no square elements. However, the realistic cell function δm~kr, which we calculate and associate to the dual band MC fundamental modes, is the superposition of the two “ideal” ones and can be found in Figures 3(a) and 3(b) by simply looking at the area (a2×a2) exactly corresponding to what is shown in Figure 1.
Out-of-plane real component of the dynamic magnetization (arbitrary units) for the two fundamental modes localized in complementary regions of the array (with B at θ = 0°), consisting of the superposition of two sublattices with different lattice constant. Here, we plot an illustrative 3×3 periodic cell extract of the ideally infinite 2D array. (a) Fundamental spin wave mode (k1=0) localized within the elements of the primary sublattice only, showing an irregular intensity distribution discussed in the text; (b) fundamental spin wave mode (k2=0) strictly localized in the secondary sublattice.
Moreover, in Figure 4 we show the magnonic gap for the two modes at k=(0,0) as a function of the applied field direction: and, as predicted from Figure 2, the gap is continuously decreasing as the applied field direction moves from 0° to 90°.
Frequency gap between the two spin wave bands as a function of the propagation direction. In this case SW propagation is always parallel to the applied field direction θ (defined in Figure 1). The gap decreases as a consequence of the frequency decrease of the upper band, due to the increasing of the demagnetizing field in the elongated elements of the lattice (secondary sublattice).
In Figure 5 we plot the dispersion curves of the SW modes propagating across either the primary (a) or secondary (b) lattice, when the external field is applied at θ=0°. Red and blue lines correspond to propagation perpendicular (with reference to thin film nomenclature, Damon-Eshbach-like) and parallel (backward-like) to the applied field, respectively [39–41]. The former increase as the wavevector magnitude is increased, and saturate, due to the full dipolar nature of the interaction among the lattice elements, in both cases. The latter, instead, decrease with increasing the wavevector magnitude, once more due to the dipolar interaction. From these curves, we compute the average group velocity vg for SWs of either of the sublattices (a is the proper lattice constant), being(2)vg=dωdk=ωπ/a-ω0π/awhich is an indication of the propagation speed of the wave along a given direction. With reference to Figure 5, we calculate different values of vg, depending on which lattice type (primary/secondary) and propagation direction k with respect to the applied field B are considered:(3)vg=2.07 m/sprimary sublattice, k⊥B11.6 m/sprimary sublattice, k∥B0.600 m/ssecondary sublattice, k⊥B1.68 m/ssecondary sublattice, k∥B.These values are rather low, considering that often SW propagation speed can get the km/s order of magnitude: however, the low speed is here a consequence of a low interaction among the MC elements, and, moreover, in MC devices, due to miniaturization, it is the SW speed difference that is crucial, in order to associate different SWs to different signals or digits. Of course, any other orientation of k with B can be, in principle, computed: the above values are shown just to give a test characterization of our choice of MC. Moreover, any other direction θ of the applied field B can be considered, which we already know to have little effects on the primary sublattice and major ones on the secondary sublattice.
Dispersion relations of SW localized within the primary sublattice ((a), a1=15 nm), and secondary sublattice ((b), a2=60 nm); in (a) and (b), the red line corresponds to the DE-like mode (k⊥B) and blue one to the BA-like mode (k∥B). Zone boundary is at k1=2.09×108 rad/m (a) and k2=5.23×107 rad/m (b). Note that in (b) the wavevector k2 was a parameter of the calculation, while in (a) k1 was extrapolated by phase map inspection.
4. Comments and Remarks
The key element in the choice of our system is twofold: from one side, the superposition of two sublattices each made up of nanoparticles with different aspect ratio, that is, with large anisotropy differences, and from the other side, the small (nanometric) size of the lattice constituents. The first aspect (anisotropy) determines the presence of two fundamental modes at k=0, localized in different sublattices; each one of these modes can be in principle independently excited, due to the large bandgap. The other aspect (nanometric size) wipes many modes out of the overall low-frequency dynamics, making the occurrence of crossing dispersion curves unlikely: with our choice of taking the micromagnetic cell as large as the nanoparticle (size 5 nm), we actually are under the so-called macrospin approximation [42], since when the dot size is lower than the exchange length (in our case, 5.5 nm), each nanoparticle acts as a single magnetic moment. With this choice, clearly we will not consider modes with a nonuniform cell function profile [43], that is, modes with nonuniform oscillations within the single nanoparticle, or localized modes (like the end modes in dots with rectangular cross section) [44, 45]: in our system these modes would occur at a very high frequency (close to or larger than 100 GHz), so with our choice we do not miss any realistic (low-frequency) solution.
Concerning the dual band feature, anisotropy differences at this small length scale are to some extent more important than the filling fraction of the system, that is, the ratio between the actual magnetic area in the periodic cell and the total area of the periodic cell itself. The filling fraction, actually, has proven to be important in determining the opening of a complete magnonic bandgap in magnonic crystals [46, 47]. Within a micromagnetic model, we can increase the filling fraction by keeping the lattice elements closer to each other, and this would increase the interdot interaction and decrease the internal field differences between the elements of each lattice: as a direct consequence, nanodots belonging to different sublattices would be less independent, so that excitations of one sublattice would interfere with the other one, breaking the dual property. Conversely, if we decrease the filling factor by keeping the lattice elements far apart, the system would behave more and more as a set of degenerate (noninteracting) dots, with very flat dispersion curves, less attractive for information delivery issues.
On the other side, for any given filling fraction, different systems characterized by increasing dot size (larger than the exchange length) would more and more tend to break the macrospin picture, due to the increasing rise of modes with complex (nonuniform) cell function, with increasing chances of crossing among different mode dispersions, so that the frequency bandgap between the fundamental modes of the two sublattices can be extremely reduced or not present at all; moreover, due to the consequent increasing interdot distance, the interdot interaction would be lower, with resulting flatter dispersion curves, with smaller group velocities.
A remarkable feature of truly dual band MCs is that it is not possible to interpret the long wavelength behavior (k=0) of the SW dynamics within the usual effective medium model [48–51]: this model approximates the system by a single effective medium with effective (intermediate) properties, with a resulting single effective frequency at k=0. But the existence of a very large fluctuation of the internal field values between the two different sublattices in the real system (at the origin of the existence of two fundamental-like SW modes with a very large frequency gap even at k=0) makes any intermediate (“effective”) value a nonphysical solution. Similar situations were discussed in other systems, with similar large fluctuations in the internal field [52–54]: in particular, in dots with triangular cross section [52] it was shown that an effective Kittel (quasi-uniform) mode was impossible to consider as an actual solution, because of the very strong variation of the internal field from the apex to the basis of the triangle.
5. Conclusion
In this paper, we proposed and characterized a model system of MC, in which SWs occur in two bands, separated by a very large frequency gap. SWs of either of the bands propagate in different regions of the MC, with no superposition, and can be considered completely independent. This allows considering the system as ideal for simultaneous propagation of two independent signals, at very different frequencies. This property, which was designed by exploiting different shape anisotropy of different MC elements, is of great importance in technology of miniaturized devices for high-density information storage/delivery, where, usually, the need for miniaturization contrasts with the other important need of signal independence, due to magnetic interactions. This work is intended to be introductory for further investigations, both theoretical and experimental, on combined systems, made of superposed lattices with special design, geometry and, possibly, different magnetic materials, all tailored in order to get a very sharp change of the internal field profile across the medium: this feature was found to be the key ingredient for the occurrence of dual band SW dynamics.
Competing Interests
The authors declare that they have no competing interests.
StampsR. L.BreitkreutzS.ÅkermanJ.The 2014 magnetism roadmap2014473333300110.1088/0022-3727/47/33/333001KrawczykM.GrundlerD.Review and prospects of magnonic crystals and devices with reprogrammable band structure2014261212320210.1088/0953-8984/26/12/123202GubbiottiG.MontoncelloF.TacchiS.Angle-resolved spin wave band diagrams of square antidot lattices studied by Brillouin light scattering20151062626240610.1063/1.4923271KrawczykM.MamicaS.MruczkiewiczM.KlosJ. W.TacchiS.MadamiM.GubbiottiG.DuerrG.GrundlerD.Magnonic band structures in two-dimensional bi-component magnonic crystals with in-plane magnetization2013464949500310.1088/0022-3727/46/49/4950032-s2.0-84890285310MamicaS.Tailoring of the partial magnonic gap in three-dimensional magnetoferritin-based magnonic crystals2013114404391210.1063/1.48166842-s2.0-84882318169ToppJ.HeitmannD.KostylevM. P.GrundlerD.Making a reconfigurable artificial crystal by ordering bistable magnetic nanowires20101042020720510.1103/PhysRevLett.104.2072052-s2.0-77953105148YuH.DuerrG.HuberR.BahrM.SchwarzeT.BrandlF.GrundlerD.Omnidirectional spin-wave nanograting coupler20134, article 270210.1038/ncomms37022-s2.0-84889246186Van de WieleB.MontoncelloF.A continuous excitation approach to determine time-dependent dispersion diagrams in 2D magnonic crystals2014473131500210.1088/0022-3727/47/31/315002Van de WieleB.HämäläinenS. J.BalážP.MontoncelloF.van DijkenS.Tunable short-wavelength spin wave excitation from pinned magnetic domain walls201662133010.1038/srep21330KlinglerS.PirroP.BrächerT.LevenB.HillebrandsB.ChumakA. V.Design of a spin-wave majority gate employing mode selection20141051515241010.1063/1.48980422-s2.0-84908052388KimS.-K.LeeK.-S.HanD.-S.A gigahertz-range spin-wave filter composed of width-modulated nanostrip magnonic-crystal waveguides200995808250710.1063/1.31867822-s2.0-69549086642ChumakA. V.VasyuchkaV. I.SergaA. A.KostylevM. P.TiberkevichV. S.HillebrandsB.Storage-recovery phenomenon in magnonic crystal20121082525720710.1103/PhysRevLett.108.2572072-s2.0-84862543365SergaA. A.ChumakA. V.AndréA.MelkovG. A.SlavinA. N.DemokritovS. O.HillebrandsB.Parametrically stimulated recovery of a microwave signal stored in standing spin-wave modes of a pagnetic film2007992222720210.1103/PhysRevLett.99.2272022-s2.0-36448958754MontoncelloF.GiovanniniL.Bandwidth broadening and asymmetric softening of collective spin waves in magnonic crystals20141042424240710.1063/1.4883878MontoncelloF.TacchiS.GiovanniniL.MadamiM.GubbiottiG.CarlottiG.SirotkinE.AhmadE.OgrinF. Y.KruglyakV. V.Asymmetry of spin wave dispersions in a hexagonal magnonic crystal20131022020241110.1063/1.48076572-s2.0-84878376557MontoncelloF.GiovanniniL.Vortex mode dynamics and bandwidth tunability in a two-dimensional array of interacting magnetic disks20121001818240610.1063/1.4710549KhvalkovskiyA. V.ApalkovD.WattsS.ChepulskiiR.BeachR. S.OngA.TangX.Driskill-SmithA.ButlerW. H.VisscherP. B.LottisD.ChenE.NikitinV.KrounbiM.Erratum: basic principles of STT-MRAM cell operation in memory arrays2013461313960110.1088/0022-3727/46/13/1396012-s2.0-84874878220ChumakA. V.SergaA. A.HillebrandsB.Magnon transistor for all-magnon data processing20145, article 470010.1038/ncomms57002-s2.0-84907359510KhitunA.BaoM.WangK. L.Magnonic logic circuits2010432626400510.1088/0022-3727/43/26/2640052-s2.0-77953859602ChumakA. V.VasyuchkaV. I.SergaA. A.HillebrandsB.Magnon spintronics201511645346110.1038/nphys33472-s2.0-84930340720HeydermanL. J.StampsR. L.Artificial ferroic systems: novel functionality from structure, interactions and dynamics2013253636320110.1088/0953-8984/25/36/3632012-s2.0-84882734272DonahueM. J.PorterD. G.1999Gaithersburg, Md, USANISTVersion 1.0GiovanniniL.MontoncelloF.NizzoliF.Effect of interdot coupling on spin-wave modes in nanoparticle arrays200775202441610.1103/physrevb.75.0244162-s2.0-33846460092KasyutichO.DesautelsR. D.SouthernB. W.Van LieropJ.Novel aspects of magnetic interactions in a macroscopic 3D nanoparticle-based crystal20101041212720510.1103/physrevlett.104.1272052-s2.0-77950207503KasyutichO.TatchevD.HoellA.OgrinF.DewhurstC.SchwarzacherW.Small angle X-ray and neutron scattering study of disordered and three dimensional-ordered magnetic protein arrays2009105707B52810.1063/1.30758652-s2.0-65249121818KasyutichO.SaruaA.SchwarzacherW.Bioengineered magnetic crystals2008411313402210.1088/0022-3727/41/13/1340222-s2.0-48249127239MontoncelloF.GiovanniniL.KrawczykM.Spin wave localization and softening in rod-shaped magnonic crystals with different terminations2012112303391110.1063/1.4743003TacchiS.GruszeckiP.MadamiM.CarlottiG.KłosJ. W.KrawczykM.AdeyeyeA.GubbiottiG.Universal dependence of the spin wave band structure on the geometrical characteristics of two-dimensional magnonic crystals201551036710.1038/srep103672-s2.0-84934882054ZivieriR.TacchiS.MontoncelloF.GiovanniniL.NizzoliF.MadamiM.GubbiottiG.CarlottiG.NeusserS.DuerrG.GrundlerD.Bragg diffraction of spin waves from a two-dimensional antidot lattice201285101240310.1103/physrevb.85.0124032-s2.0-84856454978ChumakA. V.SergaA. A.WolffS.HillebrandsB.KostylevM. P.Scattering of surface and volume spin waves in a magnonic crystal2009941717251110.1063/1.31272272-s2.0-65449184202ChumakA. V.SergaA. A.WolffS.HillebrandsB.KostylevM. P.Design and optimization of one-dimensional ferrite-film based magnonic crystals2009105808390610.1063/1.30982582-s2.0-65449136624MadamiM.BonettiS.ConsoloG.TacchiS.CarlottiG.GubbiottiG.MancoffF. B.YarM. A.ÅkermanJ.Direct observation of a propagating spin wave induced by spin-transfer torque201061063563810.1038/nnano.2011.1402-s2.0-80455174647DemidovV. E.UrazhdinS.UlrichsH.TiberkevichV.SlavinA.BaitherD.SchmitzG.DemokritovS. O.Magnetic nano-oscillator driven by pure spin current201211121028103110.1038/nmat34592-s2.0-84870064673Van De WieleB.LaursonL.FrankeK. J. A.Van DijkenS.Electric field driven magnetic domain wall motion in ferromagnetic- ferroelectric heterostructures2014104101240110.1063/1.48609632-s2.0-84892172521AuY.AhmadE.DmytriievO.DvornikM.DavisonT.KruglyakV. V.Resonant microwave-to-spin-wave transducer20121001818240410.1063/1.47110392-s2.0-84862529632KoeneB.SavoiniM.KimelA. V.KirilyukA.RasingT.Optical energy optimization at the nanoscale by near-field interference2012101101311510.1063/1.47332922-s2.0-84863685670DingJ.KostylevM.AdeyeyeA. O.Magnonic crystal as a medium with tunable disorder on a periodical lattice2011107404720510.1103/physrevlett.107.047205RossN.KostylevM.StampsR. L.Effect of disorder studied with ferromagnetic resonance for arrays of tangentially magnetized submicron Permalloy disks fabricated by nanosphere lithography2011109101390610.1063/1.3526307DamonR. W.EshbachJ. R.Magnetostatic modes of a ferromagnet slab1961193-430832010.1016/0022-3697(61)90041-52-s2.0-33751271066WhiteR. M.2007chapter 8Berlin, GermanySpringer10.1007/978-3-540-69025-2VenkatG.KumarD.FranchinM.DmytriievO.MruczkiewiczM.FangohrH.BarmanA.KrawczykM.PrabhakarA.Proposal for a standard micromagnetic problem: spin wave dispersion in a magnonic waveguide201349152452910.1109/tmag.2012.22068202-s2.0-84871823717BerkovD. V.MiltatJ.Spin-torque driven magnetization dynamics: micromagnetic modeling200832071238125910.1016/j.jmmm.2007.12.0232-s2.0-39949085772BayerC.KostylevM. P.HillebrandsB.Spin-wave eigenmodes of an infinite thin film with periodically modulated exchange bias field2006881111250410.1063/1.2184767PuszkarskiH.KrawczykM.LévyJ.-C. S.Localization properties of pure magnetostatic modes in a cubic nanograin200571101442110.1103/physrevb.71.0144212-s2.0-16844368858GrimsditchM.LeafG. K.KaperH. G.KarpeevD. A.CamleyR. E.Normal modes of spin excitations in magnetic nanoparticles2004691717442810.1103/physrevb.69.1744282-s2.0-37649027758KrawczykM.PuszkarskiH.Plane-wave theory of three-dimensional magnonic crystals20087751305443710.1103/PhysRevB.77.054437MamicaS.KrawczykM.SokolovskyyM. L.Romero-VivasJ.Large magnonic band gaps and spectra evolution in three-dimensional magnonic crystals based on magnetoferritin nanoparticles2012861414440210.1103/physrevb.86.1444022-s2.0-84867132787AlmeidaN. S.MillsD. L.Effective-medium theory of long-wavelength spin waves in magnetic superlattices198838106698671010.1103/physrevb.38.66982-s2.0-0000081405WangX.-Z.TilleyD. R.Magnetostatic modes on lateral magnetic superlattices1994187432533010.1016/0375-9601(94)90009-4WangX.-Z.TilleyD. R.Magnetostatic surface and guided modes of lateral-magnetic-superlattice films19945018134721347910.1103/PhysRevB.50.134722-s2.0-0001488146WangX.-Z.TilleyD. R.Magnetostatic modes of lateral-magnetic-superlattice films in a transverse field19979275777578610.1088/0953-8984/9/27/009MontoncelloF.NizzoliF.Spin modes of triangular magnetic nanodots in the vortex, Y, and buckle states2010107202390610.1063/1.3276228MontoncelloF.GiovanniniL.NizzoliF.VavassoriP.GrimsditchM.Dynamic origin of first and second order phase transitions in magnetization reversal of elliptical nanodots2008772121440210.1103/physrevb.77.2144022-s2.0-44649188761MontoncelloF.GiovanniniL.NizzoliF.TanigawaH.OnoT.GubbiottiG.MadamiM.TacchiS.CarlottiG.Magnetization reversal and soft modes in nanorings: transitions between onion and vortex states studied by Brillouin light scattering20087810710442110.1103/PhysRevB.78.104421