Spectrum gaps of spin waves generated by interference in a uniform nanostripe waveguide.

We studied spin waves excited by two or more excitation sources in a uniform nanostripe waveguide without periodic structures. Several distinct spectrum gaps formed by spin waves interference rather than by Bragg reflection were observed. We found the center frequency and the number of spectrum gaps of spin waves can be controlled by modulating the distance, number and width of the excitation sources. The results obtained by micromagnetic simulations agree well with that of analytical calculations. Our work therefore paves a new way to control the spectrum gaps of spin waves, which is promising for future spin wave-based devices.

Magnonic crystals, counterparts of photonic crystals, are consisting of periodic magnetic materials1. The periodic structure affects the spin waves dispersion curve by causing a forbidden band at Brillouin zone boundaries due to Bragg reflection. The first evidence of Bragg reflection of spin wave was presented in a width-modulated nanostripe waveguide2. In periodic thin ferromagnetic films, spin waves with gigahertz frequencies can be confined123456, focused and guided789101112. Such spin waves have wavelengths that are a few orders of magnitude shorter than that of electromagnetic waves at the same frequency. These properties make magnonic waveguides (MWs) promising candidates for nanoscale signal processing and information transfer devices such as filters234, delay lines1314, and storage elements15.

In recent years, many MWs designs have been presented and tested to realize signal processing functions. Usually, these MWs are fabricated with different magnetic materials1617, with different structures236 or with a local Oersted fields13. However, fabrication of periodic nanostructures involves high-precision electron-beam lithography which is very complex and expensive. On the other hand, the spin waves in these periodic structures always reflect at the Brillouin zone boundaries, resulting in an increase of the transmission loss and a decrease of the propagation distance. Therefore, how to manipulate spin waves in uniform waveguides is still an intriguing project so far.

In this paper, we present a new way to manipulate spin waves by interference in a uniform waveguide rather than by traditional Bragg reflection in periodic structure waveguide. The simple uniform waveguide with two or more excitation sources (ESs) is studied. The spin waves, excited by different ESs, interfere constructively and destructively in the uniform waveguide and form conduction bands and forbidden bands. The physical origin of spectrum gaps in such multi-sources waveguide and the relations between the center frequency of spectrum gaps and the distance, number and width of ESs are studied in detail by numerical micromagnetic simulations and analytical calculations.

Results

The sketch of the present study is shown in Fig. 1. The initial magnetization along the longitudinal axis, the backward volume magnetostatic spin waves (BVMSW) are studied in the following and the analogous results in the Damon–Eshbach (DE) mode can be found as Supplementary Fig. S4 online. As a reference, we simulate a same dimension waveguide with a single ES. The frequency spectrum, transmission characteristic and dispersion curve of spin waves are displayed in Fig. 2 (a),(b), (c), respectively. The spin wave transmission isn't allowed below 12 GHz, which attributes to the nanostripe width confinement1819. For the waveguide with a single ES, there is no spectrum gap above 12 GHz. On the contrary, there are two spectrum gaps for the two ESs waveguide as shown in Fig. 2 (d) – (f), where the distance and width of the two ESs are D = 28 nm and W = 2 nm (the length of cell), respectively.

Figure 1

Schematic of waveguide with two excitation sources.

D and W are the distance and width of the ESs, respectively.

Figure 2

Characteristics of two uniform waveguides with different number ESs.

The top and bottom panels correspond to waveguides with one and two ESs, respectively. (a), (d) The frequency spectra obtained from FFT of Mz. (b), (e) Spin wave transmission characteristics obtained by meaning the spin waves intensity of a range of the waveguides (from 1000 nm to 1200 nm). The inset frequencies indicate the positions of the spectrum gaps. (c), (f) The dispersion curves of the first and third modes obtained by micromagnetic simulation and analytical calculation (the red and black line in Fig. 2(c)). The black dotted lines denote the positions of the spectrum gaps and the corresponding wave vectors kx. The inset wavelengths correspond to the wave vectors.

To understand such spectrum gaps features, we plot the dispersion curves of waveguide with single and two ESs in Fig. 2(c) and (f), respectively. The dispersion curves obtained by micromagnetic simulations agree with that of analytical calculations202122 (see Fig. 2(c)). We can observe two distinct spectrum gaps in Fig. 2(d) – (f). The wave vectors and center frequencies of the two spectrum gaps are k1 = 0.1135 nm−1, k2 = 0.3371 nm−1 and f1 = 19.53 GHz, f2 = 97.17 GHz, respectively. In terms of the dispersion characteristics, the spectrum gaps are not caused by Bragg reflection. The slopes of the dispersion curve at the edges of the spectrum gaps keep stable shown in Fig. 2(f), whereas the Bragg reflection caused the slopes decrease obviously21523.

In order to quantitatively elucidate the physical origin of the spectrum gaps, we introduced the expression of dynamic magnetization of mz component for a stripe waveguide242526 where L is the thickness of waveguide, the magnitude h is the spatial Fourier component of microwave filed (h) corresponding to the wave vector k. In consideration of the interference between the two ESs, the resultant magnitude of the spin waves which are independently excited by different ESs is obtained: Note that we only consider the area far away from the antenna (x≫D) and the excitation fields are identical (amplitude, phase, direction and so on). The spin waves destructively interfere and disappeared (mz(sum-two) = 0) when the wavelength meets D = (2n + 1)/2 × λ, (n = 0,1,2,3,4……), where λ = 2π/k is the wavelength of spin wave, The condition for destructive interference is similar to double-slit interference in light. We calculated the wavelengths corresponding to the center frequencies of spectrum gaps. For 19.53 GHz and 97.17 GHz, the wavelengths are λ19.53GHz = 55 nm and λ97.17GHz = 18.6 nm, respectively. Taking into account the distance (D = 28 nm) between the ESs, we found D ≈ 1/2 × λ19.53GHz and D ≈ 3/2 × λ97.17GHz. Therefore, we can conclude that the spectrum gaps are mainly generated by the interference between the spin waves excited by different ESs. The spin waves, which meet the condition of destructive interference, can't propagate along waveguide and form spectrum gaps.

In order to obtain the dependence of the center frequency of spectrum gaps on the distance of ESs, we recorded and analyzed the frequency spectra at different distance of 10–56 nm. The results of these simulations agree with that of analytical calculations, as shown in Fig. 3. At D = 20 nm, there is only one spectrum gap (the center frequency 28.8 GHz) in the frequency range from 0 to 100 GHz, the wavelength of the center frequency is about 39.1 nm (D ≈ 1/2 × λ28.8GHz). More spectrum gaps can be found in this frequency range as the distance increases. The spectrum gaps at higher frequencies correspond to higher order interferences.

Figure 3

The center frequencies of spectrum gaps determined from micromagnetic simulation (square shape) and compared to the analytical calculation (line) for the different distances (D) of ESs.

The insets depict the frequency spectra for D = 20 nm, D = 38 nm and D = 52 nm. The black, blue and red line indicate the center frequencies of spectrum gaps caused by the first (D = 1/2 × λ), second (D = 3/2 × λ) and third (D = 5/2 × λ) order interferences, respectively.

The above discussions about two ESs waveguide prove that the generation of spectrum gaps of spin waves is caused by interference. We will further explore the effect of the number of ESs in this section. The transmission characteristic and dispersion curve show a distinct spectrum gap for the three ESs in Fig. 4 (a) and (b). The wave vector of the spectrum gap is 0.2094 nm−1, the corresponding wavelength can be calculated λ = 2π/k = 29.99 nm. The wavelength doesn't meet the principle obtained in two ESs interference. Taking multi-sources interference into account, the resultant magnitude of the spin waves can be expressed as: where N is the number of ESs. When the sum of spin waves intensity equal to zero (m = 0) the interference is destructive. Thus, we can drive a simple rule ND = nλ (n = 1,2,3,4,5……). If this rule is met, the interferences of the spin waves will be destructive and the spectrum gaps will be observed in the dispersion curves. For the three ESs (N = 3), the wavelength corresponding to the spectrum gap is λ = 29.99 nm. The result is closely correlated with the analytical calculations (λ = ND = 3 × 10 = 30 nm). The calculation (line) and micromagnetic simulations (dot) results are presented in Fig. 4 (c). All the data is obtained by changing the number of the ESs. The distance (D = 10 nm) and width (W = 2 nm) of the ESs are fixed during the process. Figure 4(c) shows that the developed theory above of multi-sources interference gives good quantitative description of the micromagnetic simulations.

Figure 4

(a) The transmission characteristic of spin wave obtained by meaning the spin wave intensity of a range of the waveguide (from 1000 nm to 1200 nm) for three ESs. The inset frequency indicates the position of the spectrum gap. (b) The dispersion curve of spin wave. The dashed black dot lines denote the frequency and wave vector of spectrum gap. (c) The numerical simulations (dot) and analytical calculations (line) results for the different number of ESs. The inset indicates the sketch of multi-sources waveguide.

Discussion

All the results and discussions above are based on the fixed width of ESs (W = 2 nm) which is less than the exchange length of the simulated magnetic material. However, it is very difficult to be fabricated using the existing lithography tools. Hence, the different widths of the ESs must be considered. Figure 5(a) shows the sketch of a single ES with width W. The single ES can be viewed as many narrower ESs with width λex, where λex is the exchange length. Each portion of the ES will act as a source of spin waves. The multi-sources waveguide have been discussed in above sections. In this condition, the center distance of the nearest neighbor ESs is D = λex and the number of the ESs is N = W/λex. The condition of destructive interference is We compare the calculated (line) and micromagnetic simulated (dot) data in Fig. 5(b). The results show that the developed approximate theory above of wider ES gives good quantitative description of the data of micromagnetic simulations.

Figure 5

(a) The sketch of a single ES with width W. (b) The numerical simulations (dot) and analytical calculations (line) results dependencies of the width of ESs. The insets indicate the frequency spectra of spin wave for W = 26 nm and W = 42 nm, respectively.

At last, we consider general situations in a multi-sources waveguide. Figure 6 shows the frequency spectra and transmission characteristics for different ESs. The frequency spectra with only one spectrum gap for the single ES with W = 20 nm is shown in Fig. 6(a) and with two spectrum gaps for the two ESs with W = 2 nm and D = 40 nm is shown in Fig. 6(b). As we expected, there are three spectrum gaps for two ESs with W = 20 nm and D = 40 nm, as shown in Fig. 6(c). The highest one is caused by the interference of the single ES with W = 20 nm and the lower two are caused by the interference of two ESs with W = 2 nm and D = 40 nm. The transmission characteristics of spin waves for different ESs are shown in Fig. 6(d). The blue line (W = 20 nm and D = 40 nm) is the superposition of the black (single ES W = 20 nm) and red (W = 2 nm and D = 40 nm) lines. We have previously discussed the effect of the distance and width of the ESs on spectrum gaps, respectively. Therefore, we could easily tune the spectrum gaps for a general situation in multi-sources waveguide.

Figure 6

The frequency spectra of spin wave for (a) single ES with W = 20 nm, (b) two ESs with W = 2 nm and D = 40 nm, (c) two ESs with W = 20 nm and D = 40 nm. The insets indicate the sketch of different ESs. (d) The transmission characteristics of spin waves for different ESs.

In conclusion, we performed micromagnetic simulations of the spin wave in a uniform waveguide with several ESs. The appearance of spectrum gaps was clearly observed in multi-sources waveguide. The spectrum gaps are originated from the interference between the different spin waves independently excited by different ESs. Moreover, the position and the number of spectrum gaps can be manipulated by the distance, number and width of ESs. Although our calculations were done in nanometer, similar interference conditions could be obtained in micrometer (see Supplementary Fig. S1–S3). This work is expected to provide an alternative way to manipulate spin waves besides traditional magnonic crystals.

Methods

The waveguide used in our numerical simulations is sketched in Fig. 1. The 10 nm thick permalloy waveguides with width ranging from 40 to 1000 nm and length ranging from 1300 to 6000 nm were studied in our works. Considering the slow calculation speed, we mainly studied the waveguide with width of 40 nm and length of 1300 nm. Others sizes were provided in the Supplementary Fig. S1–S3 online. Simulations were performed using OOMMF code [A version of the OOMMF code used is 1.2a5 (2012)] to solve the Landau-Lifshitz-Gilbert equation. The parameters are as follows: the exchange constant A = 13 × 10−12 J/m, the saturation magnetization Ms = 8.6 × 105 A/m, the damping constant α = 0.01 and the cell size 2 × 2 × 10 nm3. The damping parameter α in the boundary areas (~20 nm on each side) was increased more than fifty times to prevent the reflection at the ends of the structure. In order to excite spin waves with a wide frequency range, we applied a sinc function excitation with frequencies ranging from 0 to 100 GHz and an amplitude of 1.0 T. The excitation, after through the input antenna, was separated into two identical ESs (amplitude, direction and phase).

Author Contributions

Q.W., H.W.Z. and Z.Y.Z. conceived and supervised the research. Q.W. conducted the micromagnetic simulations and wrote the manuscript. G.K.M., Y.L.L. and X.L.T. made the theoretical analysis. All the authors discussed the results.