Topological Magnon Modes in Patterned Ferrimagnetic Insulator Thin

Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE

Topological Magnon Modes in Patterned Ferrimagnetic Insulator Thin Films Yun-Mei Li, Jiang Xiao, and Kai Chang Nano Lett., Just Accepted Manuscript • DOI: 10.1021/acs.nanolett.8b00492 • Publication Date (Web): 20 Apr 2018 Downloaded from http://pubs.acs.org on April 20, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Topological Magnon Modes in Patterned Ferrimagnetic Insulator Thin Films Yun-Mei Li,†,‡ Jiang Xiao,¶,§ and Kai Chang∗,†,‡,k SKLSM, Institute of Semiconductors, Chinese Academy of Sciences, P.O. Box 912, Beijing 100083, China, College of Materials Science and Opto-Electronic Technology, University of Chinese Academy of Sciences, Beijing 100049, China, Department of Physics and State Key Laboratory of Surface Physics, Fudan University, Shanghai 200433, China, Collaborative Innovation Center of Advanced Microstructures, Nanjing 210093, China, and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China E-mail: [email protected],Phone:+8682304438

∗ To

whom correspondence should be addressed Institute of Semiconductors, Chinese Academy of Sciences, P.O. Box 912, Beijing 100083, China ‡ College of Materials Science and Opto-Electronic Technology, University of Chinese Academy of Sciences, Beijing 100049, China ¶ Department of Physics and State Key Laboratory of Surface Physics, Fudan University, Shanghai 200433, China § Collaborative Innovation Center of Advanced Microstructures, Nanjing 210093, China k Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China † SKLSM,

1

ACS Paragon Plus Environment

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract Manipulation of magnons opens an attractive direction in the future energy-efficient information processing devices. Such quasi-particles can transfer and process information free from the troublesome Ohmic loss in conventional electronic devices. Here, we propose to realize topologically protected magnon modes using the interface between the patterned ferrimagnetic insulator thin films of different configurations without the Dzyaloshinskii-Moriya interaction. The interface thus behaves like a perfect waveguide to conduct the magnon modes lying in the band gap. These modes are immune to back-scattering even in sharply bent tracks, robust against the disorders and maintain a high degree of coherence during propagation. We design a magnonic Mach-Zehnder interferometer, which realizes a continuous change of magnon signal with varying external magnetic field or driving frequency. Our results pave a new way for realizing topologically protected magnon waveguide and finally achieving a scalable lowdissipation spintronic devices and even the magnonic integrated circuit.

KEYWORDS: magnon, topological phase, interfacial state, Mach-Zehnder interferometer

2


Page 2 of 18

Page 3 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Currently there are significant interests to explore topological phases in various physical systems. These topological phases are characterized by integer topological invariants, 1–4 which can not be changed continuously. Therefore the physical properties associated to the topological invariant are robust against to perturbations. The integer topological invariants such as the Chern number indicate the nontrivial band structures of bulk crystals, and the existence of topologically protected edge states, which has been proposed for applications in, for example, spintronics and quantum computation. Very recently, most efforts in this field have been devoted to extend topological concepts to photon, 5 polariton 6 and even acoustic artificial systems. 7,8 Comparing with electron system, the photonic and acoustic systems can behave more clean because of lacking of interactions between photons and acoustic waves. For device purposes, it could be more easily to construct photonic and acoustic artificial crystals with flexible and tunable edge states since the characteristic size of artificial lattices in these systems is much larger than that in electron system. As collective excitations of magnetic order in magnetic materials, magnons, i.e., spin waves, are considered as a promising candidate as information carrier in the field of spintronics, aiming at reducing energy dissipation without compromising processing speed. 9–12 Comparing with traditional magnetic metals, magnetic insulators, e.g., yttrium-iron garnet (YIG), offer us a new platform to transfer and manipulate information by magnons without using conduction electrons, avoiding the Joule heating completely. 11,13,14 The frequencies of the magnons or spin waves range from gigahertz (GHz) to terahertz (THz), and the propagation distance is up to millimeters. 15–18 The dispersion and group velocity of magnons can also be engineered, which can form magnon waveguide 19–22 and give rise to the chiral magnon modes when considering the DzyaloshinskiiMoriya interaction 22–28 or dipole interactions. 29,30 The magnon dispersions can also be controlled straightforwardly by using the magnonic crystals, in which the magnetic properties are modulated periodically in real space. 11,12,31,32 In this Letter, we propose a topologically protected magnonic waveguide using the interface between two artificial magnonic crystals with different topological charge. The artificial magnonic crystal is constructed on a ferrimagnetic thin film with triangular-shaped holes periodically located

3


Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 18

on a hexagonal lattice. The sign of the half-integer topological charge at each valley is determined by the rotation direction of the triangular holes (THs) with respect to the underlying lattice orientation. Such rotation breaks the mirror symmetry of the crystal, leading to a gap opening at the corners of the Brillouin zone. Thus an interface between two magnetic crystals with opposite rotation directions of THs is expected to hold an nontrivial topological invariant. According to the "bulk-edge correspondence", gapless chiral modes emerge in the band gap and propagate along the interface. These modes, due to the topological protection, are immune to the back-scattering, ideal for low-dissipation information transport. Utilizing the topological magnonic edge modes, we design a magnonic Mach-Zehnder interferometer, realizing a continuous change of magnon signal by changing external magnetic field or driving frequency. By proper design, the interferometer can be a magnetometer to detect magnetic field as small as the geomagnetic field. Our proposal paves the way for designing topological protected magnonic devices to achieve energy efficient information processing. As depicted in Figure 1 (a), the physical system under consideration is a hexagonal magnonic crystal consisting of triangular-shaped holes on the ferrimagnetic YIG thin film. The orientation of the THs has an angle α with the lattice vectors as defined in Fig. 1(a). An out-of-plane magnetic field Hext = Hextez (Hext > Ms ) is applied to guarantee the equilibrium magnetization along with the z direction (The magneto crystalline anisotropy energy in YIG thin film is equivalent to a perpendicular anisotropy field, 33,34 which can be absorbed into Hext ). The magnetization M(r,t) can be divided into the out-of-plane (static) and the in-plane (dynamical) components, M(r,t) ≃ Ms ez +m(r,t), where |m(r,t)| ≪ Ms . By defining a complex field m(r,t) = mx (r,t) −imy(r,t), the equation of motion for m(r,t) can be obtained from the standard Landau-Lifshitz-Gilbert equation as (see Supporting Information for the details), i(1 + iβ Ms)

∂ m(r,t) = ∂t

2 2 z ∇ m γ (Hext + Hms )m − γ Msλex

−γ Ms h(r,t),

4


(1)

Page 5 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

2 ∇)m is the exchange field with exchange coupling strength parametrized by λ 2 . where ∇ · (λex ex z h(r,t) = hx − ihy , hx , hy are the x and y component of the demagnetization field, respectively, Hms

is the z component of the magnetostatic field (see Supporting Information for the definitions). Here

γ is the electron gyromagnetic ratio, β is the Gilbert damping constant. We model the boundaries of the THs with the exchange boundary condition n · ∇m(r,t) = 0, 35 where n is the normal direction of the boundary. Similar as solving the Schrödinger equation in periodic structure, we employ the ansatz m(r,t) = mk (r)ei(k·r−ω t) , h(r,t) = hk (r)eik·r−iω t 31 and neglect the damping term. The Eq. (1) becomes a secular equation 2 ω mk (r) = γ [H0 − Msλex (2ik · ∇ + ∇2 )]mk (r)

− γ Ms hk (r),

(2)

z + M λ 2 k2 and the exchange boundary condition reads n · ∇m (r) + in · where H0 = Hext + Hms s ex k

kmk (r) = 0. Due to the non-trivial boundary configurations, we calculate the band structures numerically using the finite element method (see Supporting Information for the details). The band structures of the hexagonal magnonic crystal with rotation angle α = 0◦ and 10◦ are shown in Figure 1 (b). The dispersion shows a pair of well-defined extrema (i.e., the valleys) in each corner of the first Brillouin zone (FBZ) , i.e., in K and K′ valleys, which are related by the time-reversal symmetry. For α = 0, the magnon band structure exhibits a gapless feature at the corners of the FBZ. For α 6= 0, the broken mirror symmetry lifts the degeneracies of the valley states, inducing a finite band gap. The gapless (α = 0) and the gapped (α 6= 0) band structures are similar to the cases of graphene without and with staggered potentials on A and B sublattices. 36–38 The in-plane magnetization of the valley states, the Kv and Kc states in the upper inset of Figure 1(b), localized at the points with C3 symmetry, i.e. the A and B points, as shown in Figure 1(c). We denote the two orthogonal states as mA (r) and mB (r) with the corresponding eigenfrequencies

ωA and ωB . The reversal of the rotation direction leads to the sign change of the effective band

5


Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 18

gap ∆ = ωA − ωB = −sgn(α )∆K , where ∆K = |ωA − ωB | is the absolute band gap. The state in the q

q

vicinity of K (K ′ ) point can be expanded as: mq (r) = cA mA (r) + cB mB (r), where q = k − K(′) . By Substituting the state mq (r) into Eq. (2), one can obtain the effective low-energy Hamiltonian near K and K ′ points 39 (see Supporting Information for the derivation): Heff = ω0q + vK (qx σx + τz qy σy ) − sgn(α )

∆K σz , 2

(3)

where σ = (σx , σy , σz ) are the Pauli matrix, τz = ± denotes the valley index. ω0 = (ωKc + ωKv )/2+ q

R

2 and px = m∗ (r)(−i∂ )m (r)d 3 r. The integral η0 (q · K + 12 q2 ), vK = η0 |pxAB |, η0 = 2γ Msλex x B A AB

region is the unit cell. In order to describe the topological property of the magnonic crystal, we define the Berry curvature and Chern number in the second band as Ω (k) = ∇k × hmk |i∇k |mk i and C =

1 R 2 2π BZ d kΩ(k).

The effective Hamiltonian in Eq. (3) has been investigated in graphene system 36 and the topological charge in each valley is given by C = 21 τz sgn(α ). The Chern number of the band is zero because of the time reversal symmetry. In the artificial magnonic crystal made of THs, the sign of the topological charge in each valley can be controlled by the rotation direction of the THs. Thus we may easily create interfaces between two regions with different topologies simply by drilling THs with opposite rotation angles on the two sides of the interface, as shown in the insets in Figure 2 (a) and (b). The topological charge is C = C1 − C2 = τz , where C1 and C2 are the topological charges at the two sides of the interface, respectively. Then each valley carries an integer topological charge. As a result, gapless chiral modes are expected at their interface. To confirm the existence of the gapless states, we numerically calculate the band structures for two types of interface configurations, i.e., a horizontal and a vertical direction interface, as shown in Figure 2(a) and (b). Both configurations contain two sets of gapless states with opposite group velocities. Therefore, such interfaces form ideal waveguides for the gapless states. For the case of horizontal interface as in Figure 2 (a), the two sets of gapless states are bounded to different

6


Page 7 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

valleys, and they move towards opposite directions along the interface, similar to the graphene system with stagger potentials. 37,40 These chiral states are localized at the interface. Figure 2(c) and (d) show the propagation of these chiral states in a patterned 4-nm-thick film. 41 In our simulations, we use a radio frequency antenna to excite magnons. For α = ±10◦ , the magnons are found to be localized and propagate along the interface as shown in Figure 2(c) (d). The effective width of the waveguide is determined by the gap. A larger gap (i.e. a larger rotation angle) gives a narrower effective width. The effective width of the waveguide is about 200 nm for the angle α = ±10◦ while 150 nm for the angle α = ±20◦ and 100 nm for the angle α = ±30◦ . The values are larger than the value 5.7 nm reported in Ref. 19 and smaller than 340 nm reported in the recent experiment. 20 The chiral states lies in the band gap, i.e. near K and K′ points. The states thus carry a crystal momentum K (or K′ ). In the vertical interface (see Figure 2 (d)), only the chiral states near the red and blue points in Brillouin zone (see the down-left inset of Figure 2(d)) are allowed to propagate. In order to excite ◦

the magnons, the antenna should be placed titled respect to the horizontal axis with the angle 30 . In our simulations, the Gilbert damping effect is included in our micromagnetic simulation. The path length of the interface in Figures 2(c) and (d) is 1.2 µ m, much smaller than the decay length in YIG thin film, approximately 860 µ m , 17 therefore the damping effect is quite small in this case. The existence of the gapless states are independent of the interface configurations. Meanwhile, the states moving in opposite directions hold different topological charges, which suppresses the back-scattering during the propagation of magnons. Consequently, the magnon propagates nicely through a bent interface, as shown in Figure 3 (a). The rotation angle here is ±20◦ . The waveguide shows a narrower effective width because of a larger band gap. Because the two valleys have opposite topological charge, the inter-valley scattering is forbidden for gapless magnons. Based on this feature, we design a magnonic beam splitter using a cross-junction interface as shown in Figure 3 (b). The thin film is partitioned into four quadrants, where the THs in the first and third quadrants are rotated clockwise (negative α ), and those in the second and fourth quadrants are rotated counter-clockwise (positive α ). Such arrangement of THs forms a cross-road type

7


Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 18

interfaces. The value of topological charge associated with each path along the interfaces are marked with the solid (blue) and dashed (red) lines in the inset of Figure 3 (b) . It is evident that the east-bound magnon on path 1 cannot make its way to path 3, but can be diverted into path 2 and path 4, i.e. a beam splitter. Because of the locking of the valley degree of freedom to the propagation direction, the chiral modes are topologically protected, thus robust against the disorders. Here we consider two kinds of disorders: the uncertainty in rotation angles and sizes of the THs. In the first case, the rotation angles of THs randomly distribute between 10◦ and 20◦ but the rotation direction is kept unchanged. In the second case, the length of sides of the THs randomly distributes between 16 nm and 24 nm. It is well-known that in disordered systems, the waves are more sensitive to disorders at lower dimension. 42,43 The chiral modes are effectively one-dimensional states, therefore one would expect these states to experience Anderson localization. Nevertheless, this behavior appears to be completely evaded by the topological protected states, just like the edges states in the disordered quantum Hall systems. Figure 3 (c) and (d) show that the waveguide effect of the interface remains with either type of disorder, no localized states or standing waves formed along the interface. Above, we have proved the existence of the topological states and its robustness regardless of the shape and disorder. Now we investigate the possibility to construct the magnonic device. As previously proved, the crossroad of the interface is like a beam splitter. The beam is split into two beams at the cross. In coherent optics, the coherence between the two beams split from one beam is guaranteed when the optical path difference is smaller than the coherent length. In coherent optics, the coherence between the two beams split from one beam is guaranteed when the optical path difference is smaller than the coherent length. Different from optics, the coherence between the split two magnon beams is ensured by the back scattering immunity feature of the topological states. Or equivalently speaking, the phase of magnons keep undisturbed by the bend tracks and disorders. Therefore, we can also say that topological magnons in YIG has a very long coherent length as the damping is quite small. The interference of the split two beams in Figure 3 (b) can

8


Page 9 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

be observed by proper path design. We construct a Mach-Zehnder interferometer (MZI), as shown in Figure 4 (a). The beam is split at the left crossroad and then meets at the right crossroad. The path difference between two path is δ l ≃ 20a, where a is lattice constant. At the right crossroad, there are also two allowed propagation paths, i.e., path 3 and 4, which are marked by the blue dashed arrows. Interestingly, there is an intrinsic π phase difference between the two paths (See Supporting Information for the proof). By changing the excitation frequency or external bias magnetic field, we can change the wave vector and hence change the phase difference between the two paths. Let q be the wave vector away from K points. When q = 0 (magnetic field 344 mT) , the path 3 experiences constructive interference while the path 4 experiences destructive interference [see Fig. 4 (b)]. While q · δ l = π (magnetic field 362 mT), the interference results are reversed, as clearly shown in Fig. 4 (c). Experimentally, via a continuous adjustment of the external magnetic field Hext , one can achieve continuous change of q, resulting in a phenomenal continuous change of magnon signal at the detector. The topological modes are linear dispersed. q is proportional to the change of the magnetic field. By increasing the path difference δ l, a smaller variation of the magnetic field will remarkably change the magnon signal at the detector. For example, by extending the length of the path 2 in Figure 4 (a) (path 2 can be any shape, a battlement-shaped interface can largely increase the path length), the path difference will be increased significantly by orders of a few hundred times larger. When δ l = 4000a, the magnetic field sensitivity of the signal is as small as 1 Gs, comparable to the magnitude of the geomagnetic field. The above band structure analysis indicates that only the change of out-of-plane component of the magnetic field will affect the interference pattern. To detect the direction and magnitude of the unknown external magnetic field, we can rotate the MZI with the fixed magnetic field Hext = 344 mT kept always perpendicular to the YIG film. The interference pattern of MZI will be changed when the MZI is rotated respect to an unknown external magnetic field. From the variation of the interference pattern, one can determine the magnitude and direction of external magnetic fields. Since the interference patterns are sensitive to the magnetic fields. The magnonic MZI can thus be used as a magnetometer. The system we proposed should be quite practical to be realized experimentally using the

9


Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 18

nanofabrication techniques. 31,44–46 The topological modes are tunable by changing the lattice constant the size of THs, or by tuning the external magnetic field. In our model, the spin wave mode we discussed is in the exchange-dominated region. Since the topological mode is inherently related to the symmetry and symmetry breaking, therefore we believe the topological interfacial magnon modes still can exist for the dipolar-dominated magnons, i.e., the large lattice constant case (usually a few hundred of nanometers to a few micrometers). In summary, we proposed a topologically protected magnonic waveguide constructed from the interface between two artificial magnonic crystals with opposite rotation directions of the THs. The rotation of the THs breaks the mirror symmetry of the crystal, leading to the gap opening and a half-integer topological charge at the corners of the FBZ. The interface between the two crystals holds an integer topological charge, giving rise to the topologically protected modes that lie in the band gap and propagate along the interface. These modes are immune to back-scattering in bent tracks and disordered interfaces. We design a magnonic MZI, achieving a continuous change of magnon signal with varying external magnetic field or driving frequency. Although we consider YIG ferrimagnetic thin film here, our conclusions also hold for ferromagnetic insulators. Our proposal paves a completely new way to construct magnon-based spintronic device and even magnonic integrated circuit.

Competing interests: The authors declare no competing financial interest.

Acknowledgement Y.-M. Li thanks Weichao Yu for his helpful discussions. Y.-M. Li and K. Chang were supported by the MOST of China (Grants No. 2015CB921503, No. 2016YFE0110000 and No.2017YFA0303400 ), and the NSFC (Grants No. 11434010). J. X. was supported by the National Natural Science Foundation of China under Grant No. 11722430 and National Key Research Program of China under Grant No. 2016YFA0300702.

10


Page 11 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Supporting Information Available The calculation of the band structures, the derivation of effective Hamiltonian, the details of micromagnetic simulation and the phase calculation in Mach-Zehnder interferometer. This material is available free of charge via the Internet at http://pubs.acs.org/.

References (1) Kane, C. L.; Mele, E. J. Phys. Rev. Lett. 2005, 95, 146802. (2) Kane, C. L.; Mele, E. J. Phys. Rev. Lett. 2005, 95, 226801. (3) Bernevig, B. A.; Hughes, T. L.; Zhang, S.-C. Science 2006, 314, 1757–1761. (4) König, M.; Wiedmann, S.; Brüne, C.; Roth, A.; Buhmann, H.; Molenkamp, L. W.; Qi, X.-L.; Zhang, S.-C. Science 2007, 318, 766–770. (5) Lu, L.; Joannopoulos, J. D.; Soljaˇci´c, M. Nat. Photon. 2014, 8, 821–829. (6) Karzig, T.; Bardyn, C.-E.; Lindner, N. H.; Refael, G. Phys. Rev. X 2015, 5, 031001. (7) Yang, Z.; Gao, F.; Shi, X.; Lin, X.; Gao, Z.; Chong, Y.; Zhang, B. Phys. Rev. Lett. 2015, 114, 114301. (8) Lu, J.; Qiu, C.; Ye, L.; Fan, X.; Ke, M.; Zhang, F.; Liu, Z. Nat. Phys. 2017, 13, 369–374. (9) Chumak, A. V.; Vasyuchka, V. I.; Serga, A. A.; Hillebrands, B. Nat. Phys. 2015, 11, 453–461. (10) Kruglyak, V. V.; Demokritov, S. O.; Grundler, D. J. Phys. D: Appl. Phys. 2010, 43, 264001. (11) Serga, A. A.; Chumak, A. V.; Hillebrands, B. J. Phys. D: Appl. Phys. 2011, 43, 264002. (12) Lenk, B.; Ulrichs, H.; Garbs, F.; Münzenberg, M. Phys. Rep. 2011, 507, 107–136. (13) Yan, P.; Bauer, G. E. W. Phys. Rev. Lett. 2012, 109, 087202.

11


Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 18

(14) Barker, J.; Bauer, G. E. W. Phys. Rev. Lett. 2016, 117, 217201. (15) Schneider, T.; Serga, A. A.; Leven, B.; Hillebrands, B.; Stamps, R. L.; Kostylev, M. P. Appl. Phys. Lett. 2008, 92, 022505. (16) Kajiwara, Y.; Harii, K.; Takahashi, S.; Ohe, J.; Uchida, K.; Mizuguchi, M.; Umezawa, H.; Kawai, H.; Ando, K.; Takanashi, K. Nature 2010, 464, 262–266. (17) Yu, H.; Kelly, O. D.; Cros, V.; Bernard, R.; Bortolotti, P.; Anane, A.; Brandl, F.; Huber, R.; Stasinopoulos, I.; Grundler, D. Scientific Reports 2015, 4, 6848–6848. (18) Maendl, S.; Stasinopoulos, I.; Grundler, D. Applied Physics Letters 2017, 111, 012403. (19) Garcia-Sanchez, F.; Borys, P.; Soucaille, R.; Adam, J.-P.; Stamps, R. L.; Kim, J.-V. Phys. Rev. Lett. 2015, 114, 247206. (20) Wagner, K.; Kákay, A.; Schultheiss, K.; Henschke, A.; Sebastian, T.; Schultheiss, H. Nature Nanotechnology 2016, 11, 432. (21) Duerr, G.; Thurner, K.; Topp, J.; Huber, R.; Grundler, D. Phys. Rev. Lett. 2012, 108, 227202. (22) Lan, J.; Yu, W.; Wu, R.; Xiao, J. Phys. Rev. X 2015, 5, 041049. (23) Zhang, L.; Ren, J.; Wang, J.-S.; Li, B. Phys. Rev. B 2013, 87, 144101. (24) Mook, A.; Henk, J.; Mertig, I. Phys. Rev. B 2014, 90, 024412. (25) Chisnell, R.; Helton, J. S.; Freedman, D. E.; Singh, D. K.; Bewley, R. I.; Nocera, D. G.; Lee, Y. S. Phys. Rev. Lett. 2015, 115, 147201. (26) Chernyshev, A. L.; Maksimov, P. A. Phys. Rev. Lett. 2016, 117, 187203. (27) Kim, S. K.; Ochoa, H.; Zarzuela, R.; Tserkovnyak, Y. Phys. Rev. Lett. 2016, 117, 227201. (28) Laurell, P.; Fiete, G. A. Phys. Rev. Lett. 2017, 118, 177201.

12


Page 13 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

(29) Shindou, R.; Ohe, J.-i.; Matsumoto, R.; Murakami, S.; Saitoh, E. Phys. Rev. B 2013, 87, 174402. (30) Shindou, R.; Matsumoto, R.; Murakami, S.; Ohe, J.-i. Phys. Rev. B 2013, 87, 174427. (31) Krawczyk, M.; Puszkarski, H. Phys. Rev. B 2008, 77, 054437. (32) Krawczyk, M.; Grundler, D. J. Phys. Condens. Matter 2014, 26, 123202. (33) Wang, H.; Du, C.; Hammel, P. C.; Yang, F. Phys. Rev. B 2014, 89, 134404. (34) Jianbo, F.; Hua, M.; Wen, X.; Xue, M.; Ding, S.; Wang, M.; Yu, P.; Liu, S.; Han, J.; Wang, C.; Du, H.; Yang, Y.; Yang, J. Applied Physics Letters 2017, 110, 202403. (35) Rado, G. T.; Weertman, J. R. J. Phys. Chem. Solids 1959, 11, 315–333. (36) Xiao, D.; Yao, W.; Niu, Q. Phys. Rev. Lett. 2007, 99, 236809. (37) Yao, W.; Yang, S. A.; Niu, Q. Phys. Rev. Lett. 2009, 102, 096801. (38) Castro Neto, A. H.; Guinea, F.; Peres, N. M. R.; Novoselov, K. S.; Geim, A. K. Rev. Mod. Phys. 2009, 81, 109–162. (39) Koshino, M.; Morimoto, T.; Sato, M. Phys. Rev. B 2014, 90, 115207. (40) Semenoff, G. W.; Semenoff, V.; Zhou, F. Phys. Rev. Lett. 2008, 101, 087204. (41) d’Allivy Kelly, O. et al. Applied Physics Letters 2013, 103, 082408. (42) Anderson, P. W.; Thouless, D. J.; Abrahams, E.; Fisher, D. S. Phys. Rev. B 1980, 22, 3519– 3526. (43) Azbel, M. Y. Phys. Rev. B 1982, 26, 4735–4738(R). (44) Gates, B. D.; Xu, Q.; Stewart, M.; Ryan, D.; Willson, C. G.; Whitesides, G. M. Chem. Rev. 2005, 105, 1171–1196, PMID: 15826012. 13


Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 18

(45) Duerr, G.; Madami, M.; Neusser, S.; Tacchi, S.; Gubbiotti, G.; Carlotti, G.; Grundler, D. Appl. Phys. Lett. 2011, 99, 202502. (46) Tacchi, S.; Duerr, G.; Klos, J. W.; Madami, M.; Neusser, S.; Gubbiotti, G.; Carlotti, G.; Krawczyk, M.; Grundler, D. Phys. Rev. Lett. 2012, 109, 137202.

14


Page 15 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Figure 1: (a) Schematic of hexagonal magnonic crystal consisting of rotated triangular-shaped holes, with the rotation angle α . a1 and a2 are lattice vectors and the side length of the THs is L. A and B are the positions with C3 symmetry. (b) The magnon energy dispersions. The lattice constant |a1 | = 34.6 nm and L = 20 nm. µ0 Hext = 344mT. The red dashed line shows the band structure when α = 0, while the blue solid line indicates the one when α = 10◦ . The film thickness is 4 nm. The up-middle inset is the zoom in around the K point. The right inset shows the first Brillouin zone of the hexagonal magnonic crystal. (c) The in-plane magnetization distributions of Kv and Kc states. (d) Band gap dependence on the rotation angle.

15


Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 18

Figure 2: (color online). Magnon chiral states in the presence of the interfaces. (a) (b). Dispersions of chiral states (red curves) with horizontal direction interface (a) and vertical direction interface (b). The widths of the ribbon are 600nm(a) and 693nm(b), respectively. The rotation angle is α = ±10◦ . (c) (d). Propagation of chiral state at frequency f = 14.3 GHz in the horizontal direction interface (c) and vertical direction interface (d). The rotation angle is α = ±10◦ . The blue strip represents the radio frequency antenna generating an alternating magnetic field that excites magnons. The top right inset is the zoom in figures of the interfacial states. In (d), down right inset shows the propagation direction. The size for the lattice in (c) and (d) is 1.2 µ m×1.2 µ m.

16


Page 17 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Figure 3: (color online). Magnon propagation in various types of interfaces. The excitation frequency is f = 14.3 GHz. (a) The interface with two 90◦ corners. The symbol "-" and "+" in the down left inset denote the rotation directions of the THs. "-" means clockwise while "+" represents counter-clockwise. The size of the square is 450 nm×450 nm. The same for the subsequent figures. (b) Cross-junction interface. The rotation angle for (a) and (b) is α = ±20◦ . (c) Interfacial states with rotation α randomly distributed between 10◦ and 20◦ . The length of THs is 20 nm. (d) Interfacial states with the length of THs randomly distributed between 16 nm and 24 nm. The rotation angle is α = ±20◦ . The size of the lattice in (c) and (d) is 1.33 µ m×1.23 µ m.

Figure 4: (color online). Magnon Mach-Zehnder interferometer. (a). Schematic of Mach-Zehnder interferometer. (b) (c), Interference pattern with external magnetic field 344 mT (b) and 362 mT (c). The excitation frequency is 14.4 GHz. The rotation angle α = ±30◦ . The cell size is 1.5 µ m×1.26 µ m

17


Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

For TOC only

18


Page 18 of 18

Topological Magnon Modes in Patterned Ferrimagnetic Insulator Thin

Recommend Documents