Arxiv:2106.09542V1 [Physics.App-Ph] 15 Jun 2021 Irwv Htnculn [ from Magnon- Beneﬁted Ultra-Strong Coupling and System

arXiv:2106.09542v1 [physics.app-ph] 15 Jun 2021 irwv htnculn [ from magnon- Benefited ultra-strong coupling and system. photon strong the magnon, density, microwave in called spin role high sphere Especially, vital the YIG a in plays excitation processing. which applica- spin potential information vast shown in the has prospects features, tion unique its from eirct [ reciprocity cino annt n(rica)ao [ atom (artificial) an to magnon of action † ∗ nonreciprocity. the is notable most the of Kerr one on and [ Based out, studies of series sphere. a YIG nonlinearity, magne- the in the is anisotropy from system originated tocrystalline the magnon of of nonlinearity advantage Kerr another merits, above [ to resonator microwave coplanar ducting [ phonons plat- interaction information complex [ dissipation low forms magnon, and the extensibility tunability, of high thanks Furthermore, the to experiment. [ in magnetometry demonstrated precision been and point ceptional ao,doe tmltdb h rgtpopc,con- prospect, bright the by Stimulated circu- isolator, diode. devices: lator, nonreciprocal manufac- chip-scale and com- communication turing For in signal one-way processing. applications examples, information valuable and and Nonrecipro- technology vast munication directions. has opposite physics two in cal signal of gation lcrncaddress: Electronic lcrncaddress: Electronic 3 eetyyas aiymgoissse,benefited system, magnonics cavity years, Recently Nonreciprocity eateto ahmtc n hsc,WhnPltcncU Polytechnic Wuhan Physics, and mathematics of Department 1 colo lcrclEgneig ua oyehi Unive Polytechnic Wuhan Engineering, Electrical of School 9 orcpoa ihodrsdbnsidcdb annKerr magnon by induced sidebands high-order Nonreciprocal a esrcue i salsigteinter- the establishing via structured be can ] 11 23 – e Wang, Mei ASnmes 22.s 03.65.Ta 42.25.Bs, numbers: PACS frequ nonrecipro in on-chip applications high-performance structuring detuning potential str optimal besides have the the in sidebands order, ratio sideband high-order isolation higher the the higher and by marked st YIG the inducing th in find thereby manipulating mode we strengthened, interestingly, magnon By be the can between nonlinearity sphere. Kerr detuning (YIG) frequency garnet the iron and sys yttrium magnonics cavity one a in and nonlinearity Kerr magnon on relies nldn h rt,scn-adtidodrsdbns at sidebands, third-order and second- first-, the including 15 , 4 epooea ffcieapoc o raigrbs nonreci robust creating for approach effective an propose We e aoaoyo polcrncCeia aeil n De and Materials Chemical Optoelectronic of Laboratory Key ,otclpoos[ photons optical ], 24 .INTRODUCTION I. soitdwt h raigo Lorentz of breaking the with associated , [email protected] wangmei ,i aietda smercpropa- asymmetric as manifested is ], 1, 2 ∗ [email protected] olg fPyisadEetoi cec,HbiNra Uni Normal Hubei Science, Electronic and Physics of College inhnUiest,Whn 304 epesRpbi fC of Republic People’s 430074, Wuhan, University, Jianghan u Kong, Cui 1 19 – 6 – unsi 302 h epesRpbi fChina of Republic People’s the 435002, Huangshi, ,osraino ex- of observation ], 16 22 , aebe carried been have ] 2 17 hoY Sun, Zhao-Yu 18 ,o supercon- or ], .I addition In ]. 10 ,acoustic ], 7 , Dtd ue1,2021) 18, June (Dated: 8 have ] 3 u Zhang, Duo otpoiigoe r ae nnniert [ nonlinearity on [ based engineering are reservoir ones The promising proposed. most been have mecha- frequencies mi- optical alternative at and inte-crowave propagation decades, nonreciprocal chip-scale few realizing and towards past nisms control the Over flexibly to poses gration. and hurdle that bulky crystal, serious the magneto-optic a the is of thing structure tricky heavy the However, reciprocity. non- the [ is nonre- watched fields of most electromagnetic the study of of the One reciprocity in made physics. been ciprocal has efforts siderable in[ tion both are [ crystal nonreciprocity magneto-optic the re- optical from achieved optical and typically and Microwave microwave the gions. of in propagation fields nonreciprocal electromagnetic the with concerned mainly eyn nFrdyrtto [ Faraday-rotation on relying we h rvn edadtemgo oeas affects also be- mode detuning magnon the frequency and hand, field driving other the tween the opposite two On de- in varying system directions. the induces of therefor responses It strength asymmetric grees effective nonlinearity. the Kerr exci- determines the affects of that On YIG magnon on field. be- of applied driving power tation detuning the driving frequency the and hand, YIG) in- the one (in power and mode sphere driving magnon tween YIG the sphere, on manipulating YIG jected ob- flexibly in be by can nonlinearity sidebands tained high-order Kerr sphere. of magnon nonreciprocity material) and strong the magnonic cavities suitable to microwave another Due coupled (or YIG two a system, of is magnonics consists approach cavity which three-mode This a in frequencies. implemented microwave third-order ro- and at achieve second-, sidebands first-, to the scheme of flexible nonreciprocity bust a propose we mechanism, st,Whn 300 h epesRpbi fChina of Republic People’s the 430040, Wuhan, rsity, iest,Whn 300 h epesRpbi fChina of Republic People’s the 430040, Wuhan, niversity, nprdb rvossuis ae nnonlinearity on based studies, previous by Inspired a devices. cal 47 ogtasiso o-eirct.More non-reciprocity. transmission rong nycm-iepeiinmeasurement, precision comb-like ency e opsdo w irwv cavities microwave two of composed tem 3 irwv rqece.Ti approach This frequencies. microwave , eie uhasre fequally-spaced of series a Such regime. ne h rnmsinnonreciprocity transmission the onger 48 uYn Wu, Yu-Ying ]. ie fMnsr fEducation, of Ministry of vices rvn oe ple nYIG on applied power driving e rct fhg-re sidebands, high-order of procity h rvn ed h effective the field, driving the 3 45 n iL Zheng Li-Li and versity, , hina 46 ,prmti iemodula- time parametric ], 32 – nonlinearity 34 obekteLorentz the break to ] 25 4, , † 26 ,wihis which ], 35 27 – – 44 31 ], ] 2 magnon excitation in YIG sphere. Thus the effective II. THE MODEL AND FORMULA DERIVATION Kerr nonlinearity will be influenced, as well as degree of the asymmetric responses in two directions. More interestingly, in the optimal parameters regime, the greater We consider a three-mode cavity magnonics sys- the order of the sideband, the stronger transmission non- tem [19] as shown in Fig. 1(a). It is composed of two reciprocity marked by the enhanced isolation ratio. Such three-dimensional rectangular microwave cavities (a and a series of equally-spaced high-order sidebands can be ap- b, frequency ωc) and one YIG sphere, where collective plied in frequency comb-like precision measurement, Be- spin excitation is quantified as magnon (m, frequency sides structuring nonreciprocal devices. In addition, fre- ωm). The cavities are linearly coupled to each other with quencies of the sidebands avoid the main cavity frequency coupling rate J; and each of them is driven and detected and this will strongly inhibit dissipation. The system pa- by a bichromatic microwave source, i.e., a strong control rameters are referenced from recent experiments [19, 49], field with amplitude εa (εb) at frequency ωd and a weak 2 specially, Kerr nonlinear coefficient (K reaches 10 nHz probe field sa (sb) at frequency ωp. The YIG sphere, refereing to [50]) in the system is highly dependent on the located in one cavity, is exposed to a uniformly DC mag- specific experimental implementation. This study broad- netic bias field (Hz) in the z direction. Meanwhile, it ens applications of the cavity magnonics system in infor- interacts with one cavity with strength g, and is pumped mation processing. by a microwave driving field with amplitude ε and fre- This paper is organized as follows: In Sec. II, we quency ωd. Note that amplitudes of the input microwave introduce theoretical scheme of the three-mode cavity fields can be converted into corresponding powers, i.e., ~ 2 ~ 2 ~ 2 magnonics system and its effective Hamiltonian with ex- Pj = ωdεj , P = ωdε and Psj = ωpsj (j = a,b) with ternal pump fields. Then we derive transmission co- control power Pj , driving power P and probe power Psj . efficients and isolation ratios of the first-, second- and third-order sidebands in two cases. In Sec. III, we study In the rotating frame, including all external microwave transmission nonreciprocity of the high-order sidebands sources, the system Hamiltonian (refering to [19, 49]) can through manipulating the driving field on the YIG sphere be written as and the frequency detuning between driving field and the magnon mode in YIG sphere. Conclusions are finally drawn in Sec. IV. † † † † † Hˆtot/~ = ∆c(â aˆ + ˆb ˆb) + ∆mmˆ mˆ + Kmˆ mˆ mˆ mˆ †ˆ ˆ† ˆ† ˆ † (a) Hz +J(â b +âb )+ g(b mˆ + bmˆ ) † ˆ† ˆ b +i√ηaκaεa(â aˆ)+ i√ηbκbεb(b b) YIG m − − † −i∆pt i∆pt J s +i√ηaκasa(â e aeˆ ) b − † −i∆pt i∆pt a +i√ηbκbsb(ˆb e ˆbe ) − s b † a cavity +i√ηmκmε(ˆm mˆ ), (1) x a − y cavity

here ∆m = ωm ωd, ∆p = ωp ωc and ∆c = ωc ωd (b) − − − Cavities is constant. In addition, the frequency detunings satisfy high-order sidebands relationship ∆p = ∆c + δ, where δ = ωp ωc refereing stokes YIG − (1) (2) to Fig. 1(b). The third term in Hamiltonian is magnon field ǻc į (3) (n) Kerr nonlinear term originating from the magnetocrys- ǻp ǻp ǻp Qǻ p Ȧ talline anisotropy in the YIG sphere. The Kerr nonlinear Ȧm /Ȧd Ȧc Ȧp 2 2 coefficient K = µ0Kanγ /(M Vm), where µ0 is the magnetic permeability of free space, Kan is the first-order FIG. 1: (Color online) (a) Schematic diagram of a three-mode anisotropy constant, γ is the gyromagnetic ratio, M is cavity magnonics system consisting of two coupled microwave the saturation magnetization, and V is the volume of cavities and a YIG sphere. (b) Frequency spectrogram of m the cavity magnonics system. Here the driving field (with the YIG sphere. κj (j = a,b,m) represents the total dis- amplitude ε) applied on YIG sphere and the control fields sipation of any microwave cavity or magnon mode with (with amplitudes εa, εb) injected from microwave cavities have cavity (magnon) coupling parameter ηj . the same frequency ωd. Frequencies of the control fields and the probe fields (with amplitudes sa, sb) are respectively de- In this work, we employ analytic solution method to tuned from the cavities’ frequency by ∆c and δ. And the solve Heisenberg-Langevin equations. It not only de- frequency detuning between the probe field and the control scribes the dynamical evolution of the cavity magnon- field is ∆p. There are higher-order sidebands at frequency ics system, but also responses the expectation values of ωd + n∆p, where n is an integer representing the sideband all operators, namely, oˆ = o (ô is a normal annihilation order. operator). Taking accounth i damping progresses in the system, Heisenberg-Langevin equations of the system can be 3 expressed as two parts: a steady-state mean value and a perturbation term, i.e., a = α + δa, b = β + δb, m = ̟ + δm. It means κa a˙ = (i∆c + )a iJb + √ηaκaεa that the steady-state mean values of the operators are − 2 − determined by the strong control fields, and the pertur- −i∆pt +√ηaκasae , (2a) bation terms are caused by the weak probe fields. When κb the probe fields are absent, with dα/dt = 0, dβ/dt = 0 b˙ = (i∆c + )b iJa igm + √ηbκbεb − 2 − − and d̟/dt = 0, the mean values of the operators from −i∆pt +√ηbκbsbe , (2b) the steady-state equations can be derived as κm 2 m˙ = (i∆m + )m igb i(2K m + K)m √ηaκaεa iJβ − 2 − − | | α = −κa , (3) i∆c + +√ηmκmε. (2c) 2 κa κa g(∆c i )̟ iJ√ηaκaεa + (i∆c + )√ηbκbεb β = − 2 − 2 , κa κb 2 Equations (2a) and (2b) describe the dynamics of cav- (i∆c + )(i∆c + )+ J ities a and b. Equation (2c) describes the dynamics 2 2 (4) of the magnon mode. The cubic term 2iK m 2m in κa − | | gJ√ηaκaεa + g(∆c i )√ηbκbεb √ηmκmε Equation (2c) presents the nonlinearity of the system. ̟ = − − 2 + , κa κb 2 It is derived from the Kerr nonlinearity term in Hamil- (i∆c + )(i∆c + )+ J Θ Θ 2 2 tonian (1) with factorizing averages approximation, i.e., (5) bc = b c . The input quantum and thermal noise h i h ih i terms are ignored here because their mean values are where we have defined zero. 2 κa g (i∆c + ) κ For existence of the nonlinear term, equations (2a)-(2c) Θ = 2 + m κa κb 2 are inherently nonlinear and cannot be solved directly (i∆c + 2 )(i∆c + 2 )+ J 2 and precisely. However, due to the strength of the con- 2 +i(∆m +2K ̟ + K). (6) trol fields (5 mW) are far stronger than the probe fields | | (1 nW) in the system, we can take perturbation method In Fig. 2, evolutions of the steady-state microwave- to solve equations (2). That is dividing each operator into photon numbers and the magnon number versus the resonant driving power on magnon mode are shown in two different cases. Figure 2(a) corresponds to case 1, only cavity a is driven by the bichromatic microwave fields and cavity b is free of driven. In case 2, the same bichromatic microwave fields are only applied on cavity b in Fig. 2(b). As we can see the steady-state microwave photon numbers and the magnon number in case 1 are larger than those in case 2 even with the same driving power P . In addition, in both cases the microwave-photon numbers in cavity a are smaller than that in cavity b, though the dissipation of cavity a is lower than that of cavity b and the magnon mode, i.e., κa < κb,κm (on the basis g 1 of J > 3 and K 1nHz). The result illustrates that the counterintuitive≥ asymmetric responds of the steady- state mean values are affected by the nonlinearity term in equations (2). Meanwhile, it also plays an important role in evolutions of the perturbation terms. Expanding equations (2a)-(2c) with the ansatz of operators and removing the steady-state terms, the perturbation terms can be expressed by the linearized Heisenberg- Langevin equations

κa δa˙ = (i∆ + )δa iJδb + √η κ s e−i∆pt,(7a) FIG. 2: (Color online) The steady-state microwave photon c a a a 2 2 2 − 2 − numbers |α| , |β| and magnon number |̟| versus with the κb δb˙ = (i∆c + )δb iJδa igδm resonant driving power P on the magnon mode in two cases. − 2 − − (a) case 1: Pa = 5 mW, Psa = 1 nW, Pb = Ps = 0; (b) −i∆pt b +√ηbκbsbe , (7b) case 2: Pb = 5 mW, Psb = 1 nW, Pa = Psa = 0. Other pa- 2 κm rameters are given as ωc/2π = 10.1 GHz, ωm/2π = ωd/2π = δm˙ = [i(∆m +4K ̟ + K)+ ]δm igδb 10.06 GHz, ∆c/2π = 40 MHz, ∆m = 0, κa/2π = 3.8 MHz, − | | 2 − 2 ∗ ∗ 2 ∗ κb/2π = 5 MHz, κm/2π = 20 MHz, δ/2π = −12 MHz, i2K(̟ δm + ̟ δm +2̟δm δm) −7 − K/2π = 3.9 × 10 Hz, J/2π = 12 MHz, g/2π = 8 MHz. i2Kδm∗δm2, (7c) − 4

the quadratic terms i2K̟∗δm2, i4K̟δm∗δm in with equation (7c) denote the− system nonlinearity− and can 2 κm = i(∆m + K +4K ̟ ) , (10) not be sneezed at. They inevitably cause asymmetric re- C − | | − 2 sponses of the perturbation terms, just as the responses κa κb 2 = (∆p ∆c + i )(∆p ∆c + i ) J ,(11) of the steady-state mean values in Fig. 2. However, the S(±) ± 2 ± 2 − ∗ 2 cubic term i2Kδm δm in equations (7) is so small that 2K ̟2∗ − = S(+) . (12) it will be ignored in the following calculation. 2 κa L g (∆p + ∆c + i ) (∆p i ) (+) To study asymmetric responses of the perturbations 2 − − C S from a new perspective, we explore the nonreciprocal For the second-order upper sideband, the coefficients are transmission of the output fields, mainly including the given as first-, second- and third-order sidebands, in case 1 and κa − g(2∆p ∆c + i )M case 2. To get analytical solutions of these high-order B− = − 2 2 , (13a) sidebands, we define the perturbation terms have the fol- 2 G − lowing forms − JB2 A2 = κa , (13b) 2∆p ∆c + i − 2 − −i∆pt + i∆pt − −2i∆pt + 2i∆pt 2 ∗ ∗ −2 − +∗ δa = A1 e + A1 e + A2 e + A2 e 2KG(̟ E + ̟ M +2̟M M ) M − = 1 1 1 , − −3i∆pt + 3i∆pt 2 − 2 ∗ 2 κa +A e + A e , (8a) (i 2∆p +2K̟ D )G + g (2∆p ∆c + i ) 3 3 C− − 2 − −i∆pt + i∆pt − −2i∆pt + 2i∆pt δb = B1 e + B1 e + B2 e + B2 e (13c) − −3i∆pt + 3i∆pt + −∗ +B3 e + B3 e , (8b) M2 = DM2 + E, (13d) − −i∆pt + i∆pt − −2i∆pt + 2i∆pt δm = M1 e + M1 e + M2 e + M2 e where +M −e−3i∆pt + M +e3i∆pt, (8c) 3 3 2̟2KF D = , (14) 2 κa g (2∆p + ∆c i 2 ) F (2∆p + i ) this ansatz is rooted from the internal mechanism of − − C 2̟∗KFM +2 +4̟KFM +M −∗ the system, where a series of high-order sidebands [51] E = 1 1 1 , (15) 2 κa g (2∆p + ∆c i ) F (2∆p + i ) with frequency ωd + n∆p (n =1, 2, 3...) are generated. n − 2 − C represents the sideband order as shown in Fig. 1(b). The κb κa 2 F = (2∆p + ∆c i )(2∆p + ∆c i ) J ,(16) negative and positive exponents of the powers are in pairs − 2 − 2 − and respectively denote the upper and lower sidebands. κb κa 2 G = (2∆p ∆c + i )(2∆p ∆c + i ) J .(17) − −i∆pt For examples, the term A1 e denotes the first-order − 2 − 2 − − + i∆pt upper sideband with coefficient A1 , and A1 e cor- To the third-order upper sideband, the coefficients have responds to the first-order lower sideband (the Stokes + the following forms field [52] in Fig. 1(b)) with coefficient A1 . For the sake κa − of simplicity, in what follows we only take the high-order g(3∆p ∆c + i )M B− = − 2 3 , (18a) upper sidebands into consideration, because of similar 3 κa κb 2 (3∆p ∆c + i )(3∆p ∆c + i ) J − 2 − 2 − characteristics in the lower sidebands. Of course the cor- − responding lower sidebands can also be studied by utiliz- − JB3 A3 = κa , (18b) ing the same method. 3∆p ∆c + i − 2 2 ∗ ∗ ∗ With the ansatz we can analytically solve equa- − 2K̟ H+ + H−(3∆c + gT i ) M3 = ∗ −∗ C 2 4 , tions (7), the coefficients of the high-order sidebands can (3i∆p + igQ)( 3∆p gT + i ) 4iK ̟ C− − − C − | | be derived, respectively. Firstly, the coefficients of the (18c) first-order upper sideband are derived as 2 −∗ + 2̟ KM3 + iH+ M3 = , (18d) κa 2 3∆p gB3 i i√ηbκbsb(∆p ∆c + i )(i ∆p +2̟ K ) − − − C B− = − 2 C− L 1 2 κa 2 g (∆p ∆c + i ) + (i ∆p +2̟ K ) in which − 2 C− L S(−) 2 κa iJ√ηaκasa(i ∆p +2̟ K ) g(3∆p ∆c + i ) + C− L , Q = − 2 , (19) 2 κa 2 κa κb 2 g (∆p ∆c + i ) + (i ∆p +2̟ K ) (−) (3∆p ∆c + i )(3∆p ∆c + i ) J − 2 C− L S − 2 − 2 − κa (9a) g(3∆p + ∆c i ) T = − − 2 , (20) − κa κb 2 JB + i√ηaκasa (3∆p + ∆c i )(3∆p + ∆c i ) J A− = 1 , (9b) − 2 − 2 − 1 κa ∗ + + + −∗ + −∗ ∆p ∆c + i H = 4i̟ KM M 4i̟K(M M + M M ), − 2 + 1 2 1 2 2 1 − − − − gB1 (21) M1 = − 2 , (9c) ∗ − − − +∗ − +∗ i ∆p +2K̟ H− = 4i̟ KM M 4i̟K(M M + M M ). C− L − 1 2 − 2 1 1 2 M + = ∗M −∗, (9d) (22) 1 L 1 5

− − − − It can be seen that the coefficients A1 , B1 , A2 , B2 , here o stands for either t, τ or ℓ corresponding to j = − − A3 and B3 of the high-order sidebands are closely de- 1, 2, 3. It, Iτ and Iℓ respectively denote isolation ratio pendent on the magnon mode including the steady-state of the first-, second- and third-order sidebands. Then amplitude and the perturbation terms. we use isolation ratio Io to quantitatively describe de- With above analytical solutions, we calculate the out- gree of the nonreciprocal transmissions in two cases. It put fields of each sideband in two cases. In case 1, the is decided by the value of Io: if Io = 0, equivalently, − 2 − 2 control and probe fields are injected from cavity a and √ηbκbB / √ηaκa = 1, it means transmissions of | j | | Aj | the output fields from cavity b are explored. With input- the high-order sidebands are reciprocal (symmetry) along output relation, the output fields of high-order sidebands two opposite directions. For other situations, the higher from cavity b can be expressed as the value of Io, the stronger is the nonreciprocity of the transmissions in two cases. For example, I 1, which − −i(ωd+∆p)t − −i(ωd+2∆p)t o saout = √ηbκbB1 e + √ηbκbB2 e illustrates transmission coefficient in one direction≥ is at − −i(ωd+3∆p)t +√ηbκbB3 e . (23) least 10 times that of the other direction. It also shows that transmittances along two opposite directions have a In the direction from cavity a to b, the transmission co- strong asymmetry. This situation can be approximated efficients of the first-, second- and third-order sidebands as one-way transmission of high-order sidebands. are respectively given as

− √ηbκbB1 III. NONRECIPROCAL TRANSMISSIONS OF tba = , (24a) sa THE HIGH-ORDER SIDEBANDS IN A CAVITY MAGNONICS SYSTEM − √ηbκbB2 τba = , (24b) sa Since the transmission coeﬃcients in equations (24)

− and (26) and the isolation ratio in equation (27) are η κ B √ b b 3 mainly dominated by coeﬃcients of each sideband, which ℓba = . (24c) sa are closely associated with the steady-state amplitude

and the perturbation terms of magnon mode. Under this In case 2, the same control and probe fields are ap- situation, we analyze influence of the driving power, applied on cavity b and the output fields from cavity a are plied on the magnon mode, and the frequency detuning − − researched. In order to avoid confusion, we use j , j between the driving field and the magnon mode on the − A B nonreciprocal transmission of each sideband. and j (j =1, 2, 3) to distinguish from the coefficients in caseM 1. Combining this condition and the input-output relation, output fields of the high-order upper sidebands 15 are given as follows

− −i(ωd+∆p)t − −i(ωd+2∆p)t sbout = √ηaκa e + √ηaκa e A1 A2 10 − −i(ωd+3∆p)t +√ηaκa e . (25) A3 Similarly, along the direction from cavity b to a, the transmission coefficients of the first-, second- and third- 5 order sidebands are defined as

− √ηaκa 1 tab = A , (26a) sb 0

− 0 1 2 3 4 √ηaκa 2 τab = A , (26b) sb

− FIG. 3: (Color online) The logarithmic amplitudes log10S(ω) √ηaκa 3 ℓab = A . (26c) of the transmitted high-order sidebands vs frequency ω/∆p sb in case 1 and case 2. Here P = 15 mW, ∆p/∆c = 0.7, other parameters are same as in Fig. 2. In order to more clearly describe nonreciprocal transmissions of the high-order sidebands in two cases, we For the same purpose, but utilizing a different method, introduce definition of the nonreciprocal isolation ratio, we first exhibit an overview comparison of the high-order which is expressed as sideband spectrums between case 1 and case 2. Here − 2 the magnon mode is resonantly driven by a 15 mW driv- √ηbκbB I = log | j | , (27) ing field. As shown in Fig. 3, spectrums of the adjacent o 10 − 2 √ηaκa j sidebands are equally spaced apart from each other at | A |

6 frequency ∆p in the rotating frame. Where ω/∆p = 0 sideband in two cases are reciprocal when the driving field denotes the control field, ω/∆p = n (n = 1, 2, 3) cor- on magnon mode is absent. Increasing driving power to responds to the nth-order sideband. It can be seen 20 mW in Fig. 4(b), intuitively, only a peak with value 2 that there exists amplitude gap between the sidebands 150 of tba is seen at δ/∆c = 0.3; and a peak value | 2| − with the same order in case 1 and case 2. Scilicet, 5.6 of tab is shown at δ/∆c = 0.25. However, if one | | − 2 the transmissions of the high-order sidebands are non- zooms in locally, a local peak of tba and a tiny one 2 | | reciprocal in two cases. What’s more, amplitudes of of tab can be found at the position of the blue-dotted the high-order sidebands decrease rapidly as the side- line.| This| illustrates that the first-order sideband (probe −3 band order increases. Specifically, S(∆p)/S(0) < 10 , field) is magnified hundreds of times compared with the −7 −12 S(2∆p)/S(0) < 10 and S(3∆p)/S(0) < 10 in both original input probe field and mainly transmitted at cases. This proves that it is valid to regard high-order δ/∆c = 0.3 in case 1. Whereas, in case 2, magnifica- sidebands as perturbation compared with the strong con- tion and− transmission of the first-order sideband mainly trol field. Then, we analyze in detail the transmission occurs at δ/∆c = 0.25. From another point of view, nonreciprocity of each sideband in the following content. the transmissions of− the first-order sideband in two cases are robustly nonreciprocal near δ/∆c = 0.3. Continue to enhance the driving power to 30 mW− in Fig. 4(c), we 2 find tba with two peaks is always smaller than 2; while A. Dependence of nonreciprocal transmissions for 2| | tab reaches its maximum value 15.5 at δ/∆c = 0.31, high-order sidebands on the resonant driving power |with| the other invisible peak at the blue-dotted line.− The P results illustrate that the first-order sideband is magnified less than twice in case 1, but more than ten times 2 Figure 4 shows transmission coefficients tba and in case 2; and the strong transmission nonreciprocity ex- 2 | | tab vary with frequency detuning δ/∆c under differ- ists near δ/∆ = 0.3. From above numerical analysis, | | c ent resonant driving powers P . In Fig. 4(a), we choose we find the transmission− nonreciprocity of the first-order 2 driving power P = 0. Under this situation, tba and sideband is effectively modulated by the driving power 2 | | tab are smaller than 0.2 and evolutions of them are syn- acted on the magnon mode. This stems from the fact | | chronous. It means the transmissions of the first-order that the increased driving power excites more magnon polarons and simultaneously enhances the effective Kerr

FIG. 5: (Color online) Transmission coefficients τba, τab and ℓba, ℓab of the second- and third-order sidebands vs frequency 2 FIG. 4: (Color online) Transmission coefficients |tba| and detuning δ/∆c. In order to clearly observe the transmission 2 |tab| of the first-order sideband vs frequency detuning δ/∆c. coefficients, the values of τab in (b) and ℓab in (e) have been 2 To observe transmission coefficients clearly, |tab| in (b) and respectively magnified by three and five orders of magnitude. 2 |tba| in (c) have been increased by one order of magnitude. Here P = 0 in panels (a) and (d), P = 20 mW in panels The driving powers P = 0, 20 mW, 30 mW in (a), (b) and (b) and (e), and P = 30 mW in panels (c) and (f). Other (c), respectively. Other parameters are same as that in Fig. 2. parameters are same as that in Fig. 2. 7 nonlinearity of the magnon, which eventually strengthens of equally spaced high-order sidebands also have poten- the nonlinear asymmetric responses of the system in two tial applications in frequency comb-like precision mea- cases. In addition, it needs to be emphasized that peaks surement through simply operating the driving power. on each curve should have been resonant with the effec- Then, we quantitatively describe the degree of trans- tive supermodes made up of two microwave cavity modes. mission nonreciprocity for each sideband by isolation ra- 2 κb−κa 2 tio. Figure 6 intuitively displays isolation ratio It varies Frequencies of the supermodes shift qJ ( 2 ) ± − with the resonant driving power P and frequency detun- (equal to 0.3∆ ) from the original cavity frequency ω . c c ing δ/∆ . As the driving power increases, one bright re- This is originated± from the strong cavity-cavity coupling c gion, composed of an upper and a lower parts, appears at rate, i.e., J >κ ,κ in the system. Practically, the ac- a b the center of δ/∆ = 0.3. In this region the optimal I tural resonant frequency exhibit shifts due to the joint c t is 5.6 dB obtained with− P = 21.8 mW. This can also be interference of the strong cavity-magnon coherent inter- explained by the result in Fig. 4, i.e., the prominent large action (g>J) and the magnon Kerr nonlinearity [20]. 2 2 differences of transmission coefficients tba and tab be- Figure 5 shows transmission coefficients τba, τab and | | | | tween the left resonant peaks near δ/∆c = 0.3. ℓba, ℓab of the second- and the third-order sidebands − For the second- and third-order sidebands, the cor- versus frequency detuning δ/∆ with different resonant c responding isolation ratios I and I versus the driv- driving powers. When the driving field is absent in τ ℓ ing power P and frequency detuning δ/∆ are given in Figs. 5(a) and 5(d), although the transmission coeffi- c Fig. 7(a) and 7(b). In Fig. 7(a), the value of I in most cients are very tiny, the transmission nonreciprocity is τ areas is lower than 3.5 dB, except two bright areas cen- very obvious in the regime δ/∆ [ 0.5 0.2]. Once c tered at δ/∆ = 0.3, where I reaches its maximum the driving field is considered, both∈ the− transmission− coef- c τ 9.9 dB. While the− distribution of I is quite different. ficients and nonreciprocity will be further enhanced. For ℓ As presented in Fig. 7(b), except a small dark area, the instance, P = 20 mW in Figs. 5(b) and 5(e), τ and ba values of I in other areas are higher than 10 dB. Es- ℓ reach respective unique maximum 0.7 and 1.5 10−3 ℓ ba pecially in the yellow area centered at δ/∆ = 0.4 at δ/∆ = 0.3. While τ and ℓ keep much smaller× c c ab ab with P [0 10] mW and the other one centered− at maximums.− Continue to increase the driving power to δ/∆ = ∈ 0.3 with P [22 28] mW, I can be higher 30 mW in Fig. 5(c) and 5(f). Both τ and ℓ have c ℓ ba ba than 30 dB.− By an overall∈ observing, what I and I have very tiny maximums. While τ and ℓ have soared ℓ τ ab ab in common is that the optimal value is distributed near to their peaks 4.8 10−3 and 4.6 10−7 at the point × × δ/∆c = 0.3, this is consistent with the result in Fig. 5. δ/∆c = 0.31. This illustrates that the resonant driving − field can− not only effectively manipulate the transmis- Taken together Fig. 6 and Fig. 7, we find that the isola- sion nonreciprocity of the first-order sideband, but also tion ratios It, Iτ and Iℓ show an overall increasing trend of the second- and the third-order sidebands. What’s with the increase of the sideband order. Specifically, op- more, the transmission nonreciprocity of the high-order timal It, Iτ and Iℓ are enhanced in turn from 5.6 dB, to 9.9 dB and lastly 30 dB near detuning δ/∆c = 0.3. sidebands have a cooperative and consistent characteris- − tic when facing different driving powers. Besides being The reason is rooted from the strengthened effective Kerr applied in high performance nonreciprocal directional- nonlinearity in each sideband. Besides the steady-state switching isolator [33, 53] and diode [54–56], such a series amplitude of the magnon mode, the first-order sideband is also associated with the first-order perturbation terms (in equations (9)) of the magnon mode; the second-order 50 sideband is simultaneously related to the first- and the second-order perturbation terms of the magnon mode (in 5 40 4 50 50 30 30 40 8 40 3 25 30 6 30 20 20 2 15 20 4 20 10 10 1 10 2 10 5

0 0 0 0 0 -0.5 -0.4 -0.3 -0.2 -0.1 -0.5 -0.4 -0.3 -0.2 -0.1 -0.4 -0.2 0 0.2 0.4

FIG. 7: (Color online) isolation ratios (a) Iτ of the second- FIG. 6: (Color online) The isolation ratio It of the first-order order sideband and (b) Iℓ of the third-order sideband vs the sideband vs the resonant driving power P and frequency de- resonant driving power P and frequency detuning δ/∆c. Pa- tuning δ/∆c. Parameters are the same as that in Fig. 2. rameters are same as that in Fig. 2. 8 equations. (13)); and the third-order sideband is totally non-resonant situations and results in much weaker effec- related to the first-, second- and third-order perturbation tive Kerr nonlinearity. This lastly induces much weaker terms (in equations (18)) of the magnon mode. Except transmission nonreciprocity and shifts the optimal de- for the steady-state amplitude of the magnon, the strong tuning regime. transmission nonreciprocity of each sideband originates For actual operation to obtain strong nonreciprocity from the aggregate nonlinear effects in each perturbation of the high-order sidebands, the bandwidth of the op- term. Therefor, the higher the order of the sideband, timal detuning regime is often more concerned. In the the stronger nonlinear asymmetry of the transmission, it resonant situation, we centralize It, Iη and Iℓ in Fig. 9. lastly leads to much higher isolation ratio. It can be obviously seen It, Iη and Iℓ increase in turn

FIG. 9: (Color online) The isolation ratios It, Iτ and (c) Iℓ vs δ/∆c when the magnon mode is resonant with the driving ﬁeld, i.e., ∆m = 0. Here P = 25 mW and the other parameters are the same as in Fig. 2.

and reach their optimal values in the gray shaded area. This means an overall controlling of the strong nonreciprocity for high-order sidebands can be realized in the FIG. 8: (Color online) The isolation ratios (a) It, (b) Iτ optimal detuning regime ranging from δ/∆c = 0.301 to − and (c) Iℓ vs δ/∆c with diﬀerent detunings ∆m between δ/∆c = 0.273. The bandwidth of the optimal detun- − the magnon mode and the correspong driving ﬁeld. Here ing regime is seven megahertz with ∆c =2π 40 MHz in P = 25 mW, and the other parameters are the same as in the system. It is achievable for experimental× operation to Fig. 2. capture strong transmission nonreciprocity of the high- order sidebands simultaneously in such a bandwidth.

B. Dependence of the nonreciprocal transmissions IV. CONCLUSION for high-order sidebands on frequency detuning ∆m In summary, we have theoretically provided a method In this section, we research influence of frequency de- to realize an overall controlling of the strong transmission tuning ∆m on the nonreciprocal transmission of the high- nonreciprocity of the first-, second- and third-order side- order sidebands. In Fig. 8, Isolation ratios It, Iτ and Iℓ vs bands in a cavity magnonics system. It is consist of two δ/∆c with different detunings ∆m are intuitively given. coupled microwave cavities and one YIG (or another suit- When magnon mode is resonant with the driving field i.e., able magnonic material) sphere. Our approach utilizes ∆m = 0, isolation ratios It, Iτ and Iℓ increase in turn. the self-Kerr nonlinearity of magnon in the YIG sphere. Especially in the optimal detuning regime (centered at We have shown that strong transmission nonreciprocity the resonant frequency of one effective supermode of the of the high-order sidebands can be achieved by respec- microwave cavities [20]), It, Iτ and Iℓ reach their max- tively regulating the driving power on the magnon mode imum values. The same evolution trend is also suitable and the frequency detuning between the magnon mode for non-resonant situations, i.e., ∆m/2π = 19.2 MHz, and the driving field. We also showed that the higher the 20 MHz, except for the reduced optimal isolation− ratios order of the sideband, the stronger is the nonreciprocity and the changed optimal detuning regimes. It is because marked by the enhanced isolation ratio in the optimal de- excitation of the magnetic polarons is suppressed in the tuning regime. We have illustrated that the bandwidth 9 of the optimal detuning regime can reach several mega- Acknowledgments hertz when magnon mode is resonant drived. This im- plies that it is experimentally feasible to capture robust transmission nonreciprocity of the high-order sidebands simultaneously. This study provides a promising route to We acknowledge professor Xin-You Lüfor his valuable realize strong nonreciprocity of the high-order sidebands suggestions on our work. And this work was supported at microwave frequencies and has potential applications by the National Natural Science Foundation of China in frequency comb-like precision measurement. (NSFC) under Grants 12005078, 11675124 and 11704295.

[1] M. Goryachev, W. G. Farr, D. L. Creedon, Y. Fan, M. berg, A. J. Ramsay, A. Nunnenkamp and A. J. Fer- Kostylev and M. E. Tobar, High-Cooperativity Cavity guson, Magneto-optical coupling in whispering-gallery- QED with Magnons at Microwave Frequencies, Phys. mode resonators, Phys. Rev. A 92, 063845 (2015). Rev. Applied. 2, 054002 (2014). [18] H. Huebl, C. W. Zollitsch, J. Lotze, F. Hocke, M. Greifen- [2] S. Sharma, B. Z. Rameshti, Y. M. Blanter and G. E. W. stein, A. Marx, R. Gross, and S. T. B. Goennenwein, Bauer, Optimal mode matching in cavity optomagnonics, High cooperativity in coupled microwave resonator ferri- Phys. Rev. B 99, 214423(2019). magnetic insulator hybrids, Phys. Rev. Lett. 111, 127003 [3] C. Braggio, G. Carugno, M. Guarise, A.Ortolan and G. (2013). Ruoso, Optical Manipulation of a Magnon-Photon Hy- [19] Y.-P. Wang, G.-Q. Zhang, D. Zhang, T.-F. Li, C.-M. Hu, brid System, Phys. Rev. Lett. 118, 107205 (2017). J. Q. You, Bistability of Cavity Magnon Polaritons, Phys. [4] X. Zhang, C.-L. Zou, L. Jiang and H. X. Tang, Strongly Rev. Lett. 120, 057202 (2018). Coupled Magnons and Cavity Microwave Photons, Phys. [20] L. V. Abdurakhimov, Yu. M. Bunkov and D. Konstanti- Rev. Lett. 113, 156401 (2014). nov, Normal-Mode Splitting in the Coupled System of [5] J. T. Hou, L. Liu, Strong coupling between microwave Hybridized Nuclear Magnons and Microwave Photons, photons and nanomagnet magnons, Phys. Rev. Lett. Phys. Rev. Lett. 114, 226402 (2015). 123, 107702 (2019). [21] X. Zhang, A. Galda, X. Han, D. Jin and V. M. Vinokur, [6] O.¨ O. Soykal and M. E. Flatté, Strong Field Interactions Broadband Nonreciprocity Enabled by Strong Coupling between a Nanomagnet and a Photonic Cavity, Phys. of Magnons and Microwave Photons, Phys. Rev. Applied Rev. Lett. 104, 077202 (2010). 13, 044039 (2020). [7] N. Crescini, C. Braggio, G. Carugno, A. Ortolan and G. [22] C. Kong, H. Xiong and Y. Wu, Magnon-Induced Nonre- Ruoso, Cavity magnon polariton based precision magne- ciprocity Based on the Magnon Kerr Effect, Phys. Rev. tometry, arXiv:2008.03062. Applied. 12, 034001 (2019). [8] Y. Cao and P. Yan, Exceptional magnetic sensitivity of [23] R. J. Potton, Reciprocity in optics, Rep. Prog. Phys. 67, PT -symmetric cavity magnon polaritons, Phys. Rev. B 717-754 (2004). 99, 214415 (2019). [24] C. Caloz, A. Alù, S. Tretyakov, D. Sounas, K. Achouri [9] D. Lachance-Quirion, Y. Tabuchi, A. Gloppe, K. Us- and Z.-L. Deck-Léger, Electromagnetic Nonreciprocity, ami, Y. Nakamura, Hybrid quantum systems based on Phys. Rev. Applied 10, 047001 (2018). magnonics, Appl. Phys. Express 12, 070101 (2019). [25] C. Caloz, A. Alù, S. Tretyakov, D. Sounas, K. Achouri, [10] Y. Tabuchi, S. Ishino, A. Noguchi, T. Ishikawa, R. Ya- and Z. Deck-Léger, Electromagnetic Nonreciprocity, mazaki, K. Usami, Y. Nakamura, Coherent coupling be- Phys. Rev. Applied 10, 047001 (2018). tween a ferromagnetic magnon and a superconducting [26] V. S. Asadchy, M. S. Mirmoosa, A. Daz-Rubio,´ S. Fan qubit, Science 349, 6246 (2015). and S. A. Tretyakov, Tutorial on electromagnetic nonre- [11] J. Li, S.-Y. Zhu and G. S. Agarwal, Squeezed states of ciprocity and its origin, arXiv:2001.04848v2. magnons and phonons in cavity magnomechanics, Phys. [27] J. Fujita, M. Levy, R. M. Osgood, L. Wilkens and Rev. A 99, 021801 (2019). H. Dötsch, Waveguide optical isolator based on Mach- [12] R. F. Sabiryanov and S. S. Jaswa, Magnons and Magnon- Zehnder interferometer, Appl. Phys. Lett. 76, 2158 Phonon Interactions in Iron, Phys. Rev. Lett. 83, 10 (2000). (1999). [28] T. R. Zaman, X. Guo and R. J. Ram, Faraday rotation in [13] X. Zhang, C.-L. Zou, L. Jiang, H. X. Tang, Cavity mag- an InP waveguide, Appl. Phys. Lett. 90, 023514 (2007). nomechanics, Sci. Adv. 2, e1501286 (2016). [29] F. D. M. Haldane and S. Raghu, Possible Realization [14] Y.-P. Gao, C. Cao, T.-J. Wang, Y. Zhang and C. Wang, of Directional Optical Waveguides in Photonic Crystals Cavity-mediated coupling of phonons and magnons, with Broken Time-Reversal Symmetry, Phys. Rev. Lett. Phys. Rev. A 96, 023826 (2017) 100, 013904 (2008). [15] L. Wang, Z.-X. Yang, Y.-M. Liu, C.-H. Bai, D.-Y. [30] Y. Shoji, T. Mizumoto, H. Yokoi, I. Hsieh and R. M. Wang, S. Zhang and H.-F. Wang, Magnon blockade Osgood, Magneto-optical isolator with silicon waveg- in a PT-symmetric-like cavity magnomechanical system, uides fabricated by direct bonding, Appl. Phys. Lett. 92, arXiv:2007.14645v1. 071117 (2008). [16] S. V. Kusminskiy, H. X. Tang and F. Marquardt, Cou- [31] Z. Wang, Y. Chong, J. D. Joannopoulos and M. pled spin-light dynamics in cavity optomagnonics, Phys. Soljaˇcić, Observation of unidirectional backscattering- Rev. A 94, 033821 (2016). immune topological electromagnetic states, Nature (Lon- [17] J. A. Haigh, S. Langenfeld, N. J. Lambert, J. J. Baum- don) 461, 772 (2009). 10

[32] M. Soljacic, C. Luo, J. D. Joannopoulos and S. Fan, Non- reciprocal Microwave Signal Processing with a Field- linear photonic crystal microdevices for optical integra- Programmable Josephson Amplifier, Phys. Rev. Applied tion, Opt. Lett. 28, 637 (2003). 7, 024028 (2017). [33] A. Alberucci and G. Assanto, All-optical isolation by di- [43] K. G. Fedorov, S. V. Shitov, H. Rotzinger and A. V. rectional coupling, Opt. Lett. 33, 1641 (2008). Ustinov, Nonreciprocal microwave transmission through [34] B. Anand, R. Podila, K. Lingam, S. R. Krishnan, S. S. S. a long Josephson junction, Phys. Rev. B 85, 184512 Sai, R. Philip and A. M. Rao, Optical diode action from (2012). axially asymmetric nonlinearity in an all-carbon solid- [44] N.R. Bernier, L.D. Tóth, A. Koottandavida, M. A. Ioan- state device, Nano Lett. 13, 5771 (2013). nou, D. Malz, A. Nunnenkamp, A.K. Feofanov, T. J. [35] M. Hafezi and P. Rabl, Optomechanically induced non- Kippenberg, Nonreciprocal reconfigurable microwave op- reciprocity in microring resonators, Opt. Express 20, tomechanical circuit, Nat. Commun. 8, 604 (2018). 007672 (2012); F. Ruesink, M.-A. Miri, A. Alù, E. Ver- [45] K. Fang, J. Luo, A. Metelmann, M. H. Matheny, F. Mar- hagen, Nonreciprocity and magnetic-free isolation based quardt, A. A. Clerk and O. Painter, Generalized non- on optomechanical interactions, Nat. Commun. 7, 13662 reciprocity in an optomechanical circuit via synthetic (2016); S. Manipatruni, J. T. Robinson and M. Lip- magnetism and reservoir engineering, Nat. Phys. 13, 465 son, Optical Nonreciprocity in Optomechanical Struc- (2017). tures, Phys. Rev. Lett. 102, 213903 (2009); Z. Shen, [46] A. Metelmann and A. A. Clerk, Nonreciprocal Photon et al. Experimental realization of optomechanically in- Transmission and Amplification via Reservoir Engineer- duced non-reciprocity, 10, 657-661 (2016); A. Seif and ing, Phys. Rev. X 5, 021025 (2015). M. Hafezi, Broadband optomechanical non-reciprocity, [47] D. L. Sounas and A. Alù, Non-reciprocal photonics based Nat. Photonics 12, 60-61 (2018); B. Peng, et al., Parity- on time modulation, Nat. Photonics 11, 774 (2017). time-symmetric whispering-gallery microcavities, Nature [48] L. Ranzani and J. Aumentado, A geometric description Phys. 10, 2927 (2014). of nonreciprocity in coupled two-mode systems, New J. [36] L. Du, Y.-T. Chen, J.-H. Wu and Y. Li, Nonrecipro- Phys. 16, 103027 (2014). cal quantum interference and coherent photon routing in [49] Y.-P. Wang et al. Magnon Kerr effect in a strongly cou- a three-port optomechanical system, Opt. Express 28, pled cavity-magnon system, Phys. Rev. B 94, 224410 3647-3659 (2020); R. Huang, A. Miranowicz, J.-Q. Liao, (2016). F. Nori and H. Jing, Nonreciprocal Photon Blockade, [50] G.-Q. Zhang, Y.-P. Wang, J. Q. You, Theory of the Phys. Rev. Lett. 121, 153601 (2018). magnon Kerr effect in cavity magnonics, Sci. China-Phys. [37] B. Megyeri, G. Harvie, A. Lampis, and J. Goldwin, Di- Mech. Astron. 62, 987511 (2019). rectional Bistability and Nonreciprocal Lasing with Cold [51] H. Xiong, L.-G. Si, X.-Y. Lü, X. Yang and Ying Wu, Atoms in a Ring Cavity. Phys. Rev. Lett. 121, 163603 Carrier-envelope phase-dependent effect of high-order (2018); L. D. Bino, J. M. Silver, M. T. M. Woodley, S. L. sideband generation in ultrafast driven optomechanical Stebbings, X. Zhao and P. Del’Haye, Microresonator iso- system, Opt. Lett. 38, 353-355 (2013). lators and circulators based on the intrinsic nonreciproc- [52] Q. Bin, X.-Y. Lü, F. P. Laussy, F. Nori and Y. Wu, N- ity of the Kerr effect, Optica 5, 000279 (2018). Phonon Bundle Emission via the Stokes Process, Phys. [38] A. E. Miroshnichenko, E. Brasselet and Y. S. Kivshar, Rev. Lett. 124, 053601 (2020). Reversible nonreciprocity in photonic structures infil- [53] Y. Kawaguchi, S. Guddala, K. Chen, A. Alè, V. Menon, trated with liquid crystals, Appl. Phys. Lett. 96, 063302 A. B. Khanikaev, All-optical nonreciprocity due to (2010); valley polarization in transition metal dichalcogenides, [39] J. Mund, D. R. Yakovlev, A. N. Poddubny, R. M. arXiv:2007.14934. Dubrovin, M. Bayer and R. V. Pisarev, Magneto- [54] C. Lu, X. Hu, Y. Zhang, Z. Li, X. Xu, H. Yang, Q. toroidal nonreciprocity of second harmonic generation, Gong, Ultralow power all-optical diode in photonic crys- arXiv:2005.05393v1. tal heterostructures with broken spatial inversion sym- [40] K. Wang, Q. Wu, Y. Yu, Z. Zhang, Nonreciprocal pho- metry, Appl. Phys. Lett. 99, 051107 (2011). ton blockade in a two-mode cavity with a second-order [55] D.-W. Wang, H.-T. Zhou, M.-J. Guo, J.-X. Zhang, J. Ev- nonlinearity, Phys. Rev. A 100, 053832 (2019); ers, S.-Y. Zhu, Optical diode made from a moving pho- [41] S. Toyoda, M. Fiebig, T.-h. Arima, Y. Tokura, N. Ogawa, tonic crystal, Phys. Rev. Lett. 110, 093901 (2013). Nonreciprocal second harmonic generation in a magneto- [56] L. Fan, J. Wang, L. T. Varghese, H. Shen, B. Niu, Y. electric material, arXiv:2006.01728. Xuan, A. M. Weiner, M. Qi, An All-Silicon Passive Op- [42] F. Lecocq, L. Ranzani, G. A. Peterson, K. Cicak, R. tical Diode, Science 335, 6067 (2012). W. Simmonds, J. D. Teufel and J. Aumentado, Non-