High-Performance Magnetoinductive Directional Filters

: Multiport magnetoinductive (MI) devices with directional ﬁlter properties are presented. Design equations are developed and solved using wave analysis and dispersion theory, and it is shown that high-performance directional ﬁlters can be realised for use both in MI systems with complex, frequency-dependent impedance and in conventional systems with real impedance. Wave analysis is used to reduce the complexity of circuit equations. High-performance MI structures combining directional and inﬁnite rejection ﬁltering are demonstrated, as well as multiple-passband high-rejection ﬁltering. A new method for improving ﬁltering performance through multipath loss compensation is described. Methods for constructing tuneable devices using toroidal ferrite-cored transformers are proposed and demonstrated, and experimental results for tuneable MI directional ﬁlters are shown to agree with theoretical models. Limitations are explored, and power handling sufﬁcient for HF RFID applications is demonstrated, despite the use of ferrite materials.

Here we show that a pair of coupled right-and left-handed MI waveguides can form a directional filter and investigate its potential for signal blocking in MI-based high-frequency radio frequency identification (HF RFID), where MI antennas may also find application [46]. The filters are designed to effectively operate at high RF power, since operation is based on complementary outputs rather than absorption, and to ensure effective blocking, frequency selectivity is tuneable. Furthermore, a new method for realising four-port directional Figure 1a is a schematic of a magnetoinductive waveguide, which consists of magnetically coupled resonators with inductance L, capacitance C, resistance R, and nearestneighbour mutual inductance M. Although non-nearest-neighbour coupling is a common problem, we assume that steps have been taken to avoid it. Assuming that the loop current in resonator n is I n , the circuit equations relating the currents in adjacent elements are: jωL + 1 jωC + R I n + jωM(I n−1 + I n+1 ) = 0. (1) onics 2022, 11, x FOR PEER REVIEW 2 of 16 effective blocking, frequency selectivity is tuneable. Furthermore, a new method for realising four-port directional filters with tuneable infinite-rejection responses is presented. Wave theory is developed to explain the principle of multipath loss compensation, and high attenuation is confirmed by experiment. Designs are extended to allow multiple stopbands with improved rejection. Three-port directional filters are also shown to be feasible and can be made using a section of MI waveguide with the addition of as few as three resonators. In each case, wave theory is used to reduce the complexity of the circuit equations. Simple multiplexers can be implemented by cascading directional filters. Together, these features can introduce additional functionality and complexity in magnetoinductive systems, since the filter components can be integrated directly with MI waveguides. Additionally, such devices can be successfully used in conventional systems. Essential background is first reviewed in Sections 2.1 and 2.
2. An analytic model for MI directional filters is presented in Section 2.3, the circuit equations are solved using discrete travelling wave analysis, and methods for tuning are proposed. Conditions for infinite rejection are identified in Section 2.4, and the theory is extended to multiple passbands. Methods for constructing tuneable devices with high rejection and multiple passbands are described. Theoretical predictions are shown to be in good agreement with experimental results in the HF band in Section 3. Performance parameters including power handling and frequency scaling are also considered, with conclusions in Section 4. Figure 1a is a schematic of a magnetoinductive waveguide, which consists of magnetically coupled resonators with inductance , capacitance , resistance , and nearest-neighbour mutual inductance . Although non-nearest-neighbour coupling is a common problem, we assume that steps have been taken to avoid it. Assuming that the loop current in resonator is , the circuit equations relating the currents in adjacent elements are:

Magnetoinductive Waveguides
(   For convenience, we define the self-impedance Z = jωL + 1 jωC + R, mutual impedance Z m = jωM, and magnetic coupling coefficient κ = 2M L . Assumption of a travelling current wave in the form I n = I 0 e −jkna , where k = k − jk is the complex propagation constant and a is the element spacing, yields the well-known dispersion relation [34]: Here is the angular resonant frequency and Q = ω 0 L R is the quality factor. For lossless waveguides, the Q-factor is infinite, and the left-hand side of Equation (2) is real with k = 0. In real systems, Q must be finite. If losses are low (k k ), dispersion and loss may be approximated using: Equation (3) is the dispersion relation for lossless guides. It implies that MI waveguides are band limited, with propagation in the range Guides with positive M are right-handed and have positive phase velocity, while negative M yields left-handed guides with negative phase velocity. At the resonant frequency, k a = π 2 when M > 0 and − π 2 when M < 0. Equation (4) implies that losses are inversely proportional to |κQ|, are lowest at the resonant frequency, and rise rapidly at the band edges. Figure 1b shows the frequency dependence of k a and k a for Q = 200 and two different values of |κ|, which clearly show the band-limited nature of propagation. Crosses show three experimentally measured values along the dispersion relation for κ = 0.3 and are clearly in agreement with Equation (3).

Impedance and Power Flow
It is simple to show that the characteristic impedance of an MI waveguide is [47]: If the waveguide is lossless, the characteristic impedance Z 0M at resonance is purely real, with Z 0M = ω 0 M. In general, the impedance is complex valued and frequency dependent, which implies that care must be taken when matching to real-valued loads. An important consequence is that conventional scattering parameters cannot be used; instead, the formulation of 'power waves', originally developed by Kurokawa [50] and recently adapted to MI waveguides [51], is required. For MI lines extending to the left, type 1 power waves A1 and B1 are used at element n, as shown in Figure 1a, defined as: Here Z * 0 represents the complex conjugate of Z 0 , and Re(Z 0 ) its real part. For MI lines extending to the right, type 2 power waves A2 and B2 must be used, defined as: The squared moduli of these coefficients provide the correct definition of power.

Magnetoinductive Directional Filters
A directional filter has four ports connected to lines of characteristic impedance Z 0 . A voltage source connected to port 1 then gives rise to incident and reflected waves A1 1 and B1 1 and transmitted waves B2 2 , B2 3 , B1 4 from ports 2, 3, and 4. Port 2 has the response of a band-reject filter, and Port 4 that of a bandpass filter. In a matched and lossless device, S 21 and S 41 are complementary, S 31 = 0 for all ω, and S 11 = 0 at resonance. Equivalent properties apply no matter which port is used as input. We now develop the analysis for MI directional filters and consider their frequency response, bandwidth, and tuneability. Figure 2a shows the equivalent circuit of a filter based on a pair of infinite horizontal MI waveguides consisting of identical resonators, but with opposite coupling polarity so that the lower guide is right-handed (blue) and the upper guide left-handed (red). The mutual inductances are equal in magnitude, with M > 0, so the coupling coefficients are κ f = 2M L and κ b = − 2M L . The waveguides are connected vertically via two resonators, which couple the guides via mutual inductances M C and coupling coefficients κ c = 2M c L , such that κ c < κ f . The network has four ports, 1-4, formed by sections of an MI waveguide. The lower line supports currents I n , and the top line supports J n , while the two coupling resonators have currents K and P. Note that guide terminations can also be represented by lumped impedances, as shown in Figure 2b. There are three main possibilities: an infinite MI waveguide, represented by Z 0 (1); a resistive termination, represented by Z 0M (2); and a broadband termination to Z 0M , achieved by halving the inductance and doubling the capacitance in the resonant element (3). Termination 2 provides an impedance match at ω 0 , while termination 3 provides matching at two frequencies, ω 0 and ω 0 √ 1−κ 2 , with low reflection in between [49]. develop the analysis for MI directional filters and consider their frequency response, bandwidth, and tuneability. Figure 2a shows the equivalent circuit of a filter based on a pair of infinite horizontal MI waveguides consisting of identical resonators, but with opposite coupling polarity so that the lower guide is right-handed (blue) and the upper guide left-handed (red). The mutual inductances are equal in magnitude, with > 0, so the coupling coefficients are = 2 and = − 2 . The waveguides are connected vertically via two resonators, which couple the guides via mutual inductances and coupling coefficients = 2 , such that < . The network has four ports, 1-4, formed by sections of an MI waveguide. The lower line supports currents , and the top line supports , while the two coupling resonators have currents and . Note that guide terminations can also be represented by lumped impedances, as shown in Figure 2b. There are three main possibilities: an infinite MI waveguide, represented by 0 (1); a resistive termination, represented by 0 (2); and a broadband termination to 0 , achieved by halving the inductance and doubling the capacitance in the resonant element (3). Termination 2 provides an impedance match at 0 , while termination 3 provides matching at two frequencies, 0 and 0 √1 − 2 , with low reflection in between [49]. For now, we assume type 1 terminations and that an MI wave is incident on port 1 from the left. Experience suggests that this wave will be partially reflected at = 0 due to the discontinuity caused by the coupling resonator, and that the remaining power, less any loss, will appear as MI waves from ports 2 to 4. We write expressions for the currents and in terms of incident, reflected, and transmitted MI waves as: Here and are the amplitudes of the incident and reflected waves at port 1; 2 , 3 , 4 are the wave amplitudes from ports 2, 3, and 4; = ′ − ′′ is the propagation constant of the forward waveguide; and = ′ − ′′ is the propagation constant of the backward one. From Figure 1b and Equation (2), we can assume: Combining Equation (1) with circuit equations for the currents and in the coupling resonators, a set of coupled equations may be obtained at the six central resonators: For now, we assume type 1 terminations and that an MI wave is incident on port 1 from the left. Experience suggests that this wave will be partially reflected at n = 0 due to the discontinuity caused by the coupling resonator, and that the remaining power, less any loss, will appear as MI waves from ports 2 to 4. We write expressions for the currents I n and J n in terms of incident, reflected, and transmitted MI waves as: Here I and R are the amplitudes of the incident and reflected waves at port 1; T 2 , T 3 , T 4 are the wave amplitudes from ports 2, 3, and 4; k f = k f − jk f is the propagation constant of the forward waveguide; and k b = k b − jk b is the propagation constant of the backward one. From Figure 1b and Equation (2), we can assume: Combining Equation (1) with circuit equations for the currents K and P in the coupling resonators, a set of coupled equations may be obtained at the six central resonators: Substituting the assumed solutions and using the dispersion equation, K and P may be eliminated, allowing the reflection and transmission coefficients to be calculated: Here µ = M c M is the mutual inductance ratio. Note that the coefficients above are complex valued and frequency dependent. They are expressed only in terms of µ and k f , which may be found from Equation (2). However, care must be taken to ensure that k f ≥ 0 and 0 ≤ k f ≤ π. The results imply that no signal leaves port 2 at resonance, provided the guides are lossless. Similarly, the transmission from port 4 is maximum at resonance, while the reflectance from port 1 is zero. Port 3 is isolated at all frequencies, both within and outside the MI passbands. This behaviour is inherent in the way signals are combined. The phase shift at resonance along path P1 in Figure 2a is π 2 , while the phase shift along path P2 is 3π 2 . If the device is lossless, the signal from port 2 is then zeroed by cancellation. If Q-factors are finite, losses are different along P1 and P2 and cancellation is partial. Two signal paths, P3 and P4, to port 3 are in antiphase for all frequencies.
Although the results above provide valid amplitude reflection and transmission coefficients, they offer no insight on power flow in devices that are either lossy or operating out of band. We therefore define three sets of discrete power wave pairs, {A1 1 , B1 1 }, {A2 2 , B2 2 }, and {A1 4 , B1 4 }, for ports 1, 2, and 4, as shown in Figure 2a, as follows: Here Z 0 f = jωMe −jk f a and Z 0b = −jωMe −jk b a are the characteristic impedances of the forward and backward guides. In fact, Z 0 f = Z 0b whenever Equation (9) applies. Power wave scattering parameters with respect to port 1 may then be found as: Using the transmission and reflection coefficients, we may obtain analytic expressions for the scattering parameters of a directional filter with MI line terminations (1) as: Moduli squared of Equations (11) and (14) are equal whenever k f = k b = 0.

Frequency Response
To illustrate the above, Figure 3a shows the frequency dependence of magnitude S-parameters for different terminations, with L = 1 µH, κ = 0.6, and κ c = 0.3 (so that µ = 1/2). Results for lossless and lossy elements with magnetoinductive terminations are Electronics 2022, 11, 845 6 of 16 obtained using Equation (14), and results for lossy elements with narrowband resistive terminations using numerical solution of circuit equations.

Frequency Response
To illustrate the above, Figure 3a shows the frequency dependence of magnitude Sparameters for different terminations, with = 1 , = 0.6, and = 0.3 (so that = 1 2 ⁄ ). Results for lossless and lossy elements with magnetoinductive terminations are obtained using Equation (14), and results for lossy elements with narrowband resistive terminations using numerical solution of circuit equations. In each case, frequencies near 0 are diverted into port 4, while the remaining power is sent to port 2. Port 4 is therefore bandpass, while port 2 is a bandstop. The MI passband is as estimated from previous arguments. The isolation of port 3 is unaffected by loss; however, power emerging from port 4 is reduced and port 2 rejection diminished. Since losses are inversely proportional to and , high coupling coefficients and Q-factors are required to minimise loss. The effect of a narrowband resistive termination is to reduce the bandwidth. However, this can be partly restored by using a broadband termination.

Bandwidth
An analytic expression for the 3-dB bandwidth of | 41 | 2 may be obtained in the lossless case. The half-power frequencies can be found as follows: Within the passband, there are two solutions for ′ . If these are expressed in terms of , , and 0 , the half-power bandwidth is given by: Thus, the width of the bandpass channel increases with and 0 . It is simple to calculate the Q-factor of the bandpass response, as = 0 . Figure 3b shows a contour map of the variation of with and thus obtained. High-filter Q-factors require high values of and low values of ; however, the achievable Q-factor is limited since the use of real resonators with finite Q-factors will result in high losses if is too low. In each case, frequencies near ω 0 are diverted into port 4, while the remaining power is sent to port 2. Port 4 is therefore bandpass, while port 2 is a bandstop. The MI passband is as estimated from previous arguments. The isolation of port 3 is unaffected by loss; however, power emerging from port 4 is reduced and port 2 rejection diminished. Since losses are inversely proportional to κQ and κ c Q, high coupling coefficients and Q-factors are required to minimise loss. The effect of a narrowband resistive termination is to reduce the bandwidth. However, this can be partly restored by using a broadband termination.

Bandwidth
An analytic expression for the 3-dB bandwidth of |S 41 | 2 may be obtained in the lossless case. The half-power frequencies can be found as follows: Within the passband, there are two solutions for k f . If these are expressed in terms of κ, µ, and ω 0 , the half-power bandwidth B is given by: Thus, the width of the bandpass channel increases with µ and ω 0 . It is simple to calculate the Q-factor of the bandpass response, as Q bp = ω 0 B . Figure 3b shows a contour map of the variation of Q bp with κ and κ c thus obtained. High-filter Q-factors require high values of κ and low values of κ c ; however, the achievable Q-factor is limited since the use of real resonators with finite Q-factors will result in high losses if κ c is too low.

Frequency Tuning
So far, the filter centre frequency has been determined by ω 0 , which also determines the operation of the MI waveguides. It would clearly be beneficial to introduce tuning within the MI passband. Figure 3c shows the result of tuning the resonance of the coupling resonators to different values, namely, (a) ω 0 = 0.9ω 0 , (b) ω 0 = 0.95ω 0 , and (c) ω 0 = 1.1ω 0 , assuming lossless resonators and lossless MI terminations. Although the best performance is obtained when ω 0 = ω 0 , tuning around ω 0 is feasible for the variations in ω 0 shown.
From the above, we conclude that the network shown in Figure 2a with lossless terminations possesses all the properties of a directional filter except that it has null reflectivity only within the MI passband. S 11 is zero at the resonant frequency and very small in its vicinity. Lossy systems and systems with narrowband terminations have similar properties, the main difference being a reduction in the magnetoinductive bandwidth.

Advanced Filters
In this section, we consider the properties of more advanced filters, which offer either infinite rejection or multiple bandstop frequencies.

Infinite Rejection
In the absence of loss, complete rejection occurs at ω 0 at port 2 when port 1 is used as an input. However, for finite Q-factors, the two signal paths from port 1 to port 2 in Figure 2a incur unequal losses, and some leakage into port 2 then occurs even at resonance. To equalise signal amplitudes along the two paths, the magnetic coupling strength between resonators 0 and 1 in the forward waveguide may be adjusted. We define a new positive coefficient, κ f = λκ f , between these elements and reconsider the port reflection and transmission coefficients. The circuit equations must be modified to: As before, we may substitute the assumed solutions and use the dispersion relation to obtain modified reflection and transmission coefficients. The transmission from port 2 may then be expressed as a rational polynomial in λ as follows: For infinite rejection at midband, we require the numerator of T 2 /I to equal zero: At resonance, the propagation constant can be rewritten in terms of κ and Q alone: Substituting into Equation (19), we obtain the condition for infinite rejection as: Figure 4a shows the variation of λ with µ for a range of Q-factors. Lower values of λ are required for higher losses, as well as for lower values of µ. In each case, λ → 1 as µ increases, with sharper increases at lower losses. When λ is set as above, port 2 is isolated at resonance. However, there is an impact on filter performance; namely, port 3 is no longer fully isolated, the passband of port 2 reduces, and reflectance at port 1 is improved.  Figure 4b shows the frequency dependence of scattering parameters for filters with lossy MI (1) and resistive narrowband (2) terminations for κ = 0.6, κ c = 0.3, Q = 200, and λ = 0.8746 from Equation (21). Both cases demonstrate excellent directional filtering with reflectance mostly <−30 dB within the passband, as well as infinite rejection at resonance. If isolation is important, λ may be adjusted to balance rejection with leakage into port 3. The mutual inductance discontinuity generated by tuning λ implies that the complementary waveguide pair is no longer symmetrical, meaning the ports cannot be used interchangeably. Only the two ports immediately adjacent to the λ-tuned resonator pair would yield an infiniterejection response. The remaining two ports would suffer from reduced performance if used as inputs. Figure 4a shows the variation of with for a range of Q-factors. Lower values of are required for higher losses, as well as for lower values of . In each case, → 1 as increases, with sharper increases at lower losses. When is set as above, port 2 is isolated at resonance. However, there is an impact on filter performance; namely, port 3 is no longer fully isolated, the passband of port 2 reduces, and reflectance at port 1 is improved. Figure 4b shows the frequency dependence of scattering parameters for filters with lossy MI (1) and resistive narrowband (2) terminations for = 0.6, = 0.3, = 200, and = 0.8746 from Equation (21). Both cases demonstrate excellent directional filtering with reflectance mostly <−30 dB within the passband, as well as infinite rejection at resonance. If isolation is important, may be adjusted to balance rejection with leakage into port 3. The mutual inductance discontinuity generated by tuning implies that the complementary waveguide pair is no longer symmetrical, meaning the ports cannot be used interchangeably. Only the two ports immediately adjacent to the -tuned resonator pair would yield an infinite-rejection response. The remaining two ports would suffer from reduced performance if used as inputs.

Multiple Bandstop Frequencies
A directional filter with complementary bandstop and bandpass responses can be realised by increasing the number of coupling resonators from 1 to , as shown in Figure  5a. Here MI waveguides with coupling coefficients ± are connected using vertical lines containing weakly coupled resonators with coupling . The resonator rows are labelled from 0 (bottom) to top ( − 1) and have the currents and flowing in the left and right columns, respectively. The lines of coupling resonators can be considered as two identical uniform MI waveguides with a passband of  The solutions for currents and remain unchanged from Equation (8). Once again, we may write down the equations at the junctions as: 0 + ( −1 + 1 ) + 0 = 0, 1 + ( 0 + 2 ) + 0 = 0 The solutions for currents I n and J n remain unchanged from Equation (8). Once again, we may write down the equations at the junctions as: Equation (23) can again be solved for the reflection coefficient at port 1 and transmission coefficients at ports 2-4 by substituting the dispersion relations and eliminating the terms P T , P R , and K T , K R . When N = 1, Equation (24) is identical to Equation (11). R I = (e j2k c a(N+1) − 1) e j2k c a 1 + e j2k f a − µe j(k f +k c )a 1 + e j2k c a + µ 2 e j2(k f +k c )a /D The coefficients can be converted to scattering parameters as before. Figure 5b shows the frequency dependence of S-parameters for N = 2 and N = 5. The additional coupling resonators generate multiple bandpass responses in S 41 and complementary notches in S 21 . The value of |κ c | determines the sharpness of the peaks in S 41 , as well as the notch frequencies, with a larger separation between peaks for higher values. The quality factors of the peaks are higher for increasing N, and the effect of loss is to reduce rejection in S 21 and increase insertion loss in S 41 .
The frequencies of peaks and notches can be found by computing the maxima of S 41 or the minima of S 21 within the MI passband. The central frequencies ω k simplify to: The notch spacing is nonuniform and can be adjusted by varying ω 0 or κ c . The values for ω k are independent of κ, provided that κ > κ c . When N is odd, there is always a peak at the resonant frequency. As previously shown, infinite rejection can be achieved at resonance for N = 1. However, since path losses cannot be equalised at several frequencies, infinite rejection cannot be obtained for all notches simultaneously when N > 1. However, simulations show that it is simple to improve performance by adjusting λ.

Three-Port Filters and Multiplexers
In the absence of infinite or improved rejection (λ = 1), the current in port 3 is identically zero at all frequencies with port 1 as an input. This implies that port 3 can be physically removed from either the single-or multiple-notch device, along with a resonator at n=1 in the backward-wave guide in Figure 2a or Figure 5a. Hence, a threeport directional filter can be generated with the addition of just three weakly coupled resonators, and a backward-wave resonator at n = 0 loaded with one of the terminations from Figure 2b. Loss compensation, however, is then impossible due to the removal of signal paths. Nevertheless, filter performance would be indistinguishable from Figure 3a or Figure 5b, as long as λ is unity. Multiplexers can also be implemented by cascading filters, with the notch frequencies of the constituent filters defining the frequencies of extracted signals.

Results and Discussion
In this section, we present experimental results for filters constructed for f 0 = 13.56 MHz to remove residual carrier from tag responses in HF RFID. A tuneable design was developed to allow the demonstration of filter operation and extension to complex filters.

Single-Notch Filter
The overall design was as shown in Figure 2a; however, resistive terminations were used, and ports were floated using paired capacitors. To allow frequency tuneability, all capacitors included varicaps. To confine magnetic fields and ensure nearest-neighbour coupling, all inductors used low-loss ferrimagnetic toroidal cores with inner and outer diameters of 13.2 and 21 mm (Fair-Rite 5967000601, based on "67" Ni-Zn ferrite material with cutoff at 50 MHz). The initial relative permeability was ≈40, implying that the magnetic flux was almost completely contained in the core. Two windings (with inductances L C and L T , such that L C + L T = L) were used for the waveguide elements, and three (with inductances L C , L C , and L S , such that 2L C + L S = L) for the coupling resonators. Adjacent windings were reversed to obtain negative coupling coefficients. A method was used to allow the tuning of mutual inductance using components constructed by 3D printing with a polycarbonate filament. Figure 6a shows this mechanism, allowing control over the relative position of two copper windings. The coupling coefficient was a nonlinear function of the winding separation with minimum separation yielding maximum coupling. signal paths. Nevertheless, filter performance would be indistinguishable from Figure 3a or 5b, as long as is unity. Multiplexers can also be implemented by cascading filters, with the notch frequencies of the constituent filters defining the frequencies of extracted signals.

Results and Discussion
In this section, we present experimental results for filters constructed for 0 = 13.56 MHz to remove residual carrier from tag responses in HF RFID. A tuneable design was developed to allow the demonstration of filter operation and extension to complex filters.

Single-Notch Filter
The overall design was as shown in Figure 2a; however, resistive terminations were used, and ports were floated using paired capacitors. To allow frequency tuneability, all capacitors included varicaps. To confine magnetic fields and ensure nearest-neighbour coupling, all inductors used low-loss ferrimagnetic toroidal cores with inner and outer diameters of 13.2 and 21 mm (Fair-Rite 5967000601, based on "67" Ni-Zn ferrite material with cutoff at 50 MHz). The initial relative permeability was ≈40, implying that the magnetic flux was almost completely contained in the core. Two windings (with inductances and , such that + = ) were used for the waveguide elements, and three (with inductances , , and , such that 2 + = ) for the coupling resonators. Adjacent windings were reversed to obtain negative coupling coefficients. A method was used to allow the tuning of mutual inductance using components constructed by 3D printing with a polycarbonate filament. Figure 6a shows this mechanism, allowing control over the relative position of two copper windings. The coupling coefficient was a nonlinear function of the winding separation with minimum separation yielding maximum coupling.   Figure 6b shows a photograph of the completed filter constructed as a copper-clad printed circuit board using low-loss mica capacitors. Mechanical tuning was provided for all mutual inductances. The value of M needed for impedance matching to a 50 Ω load at resonance was first identified. Inductors L s were then chosen to equalise the self-inductance of all resonators, and capacitances were tuned to achieve resonance at 13.56 MHz. Table 1 shows the values of the most important parameters after tuning. Four-port scattering parameters were measured using an electronic network analyser (Agilent E5061B). To begin with, λ was set to unity so that responses were independent of the port used as input. Figure 7a compares the frequency dependence of S-parameters with the predictions of the theoretical model for κ c = 0.28. Excellent agreement with experimental results can be seen. Larger values of µ were investigated and found to reduce losses and increase the passband of the bandpass channel and the notch depth.

Single-Notch Filter with Infinite Rejection
Stopband rejection can be greatly improved by fine-tuning according to Equation (21). For parameters in Table 1, the value of required for infinite rejection is 0.79, well within range for the coupling tuning mechanism. Figure 7b compares the simulated and measured scattering parameters obtained in this case. A performance improvement of over 50 dB is achieved compared with that in Figure 7a, with 70 dB rejection at port 2. This improvement is at the expense of reduced isolation of port 3 with a minimum attenuation of 18 dB at resonance between ports 1 and 3. However, losses in 21 are still extremely low.

Double-Notch Filters
A two-notch directional filter with parameter values shown in Table 1 was implemented in the same way by increasing the number of coupling resonators in each vertical branch from one to two. The magnitude of determines the frequency locations of the notches in 21 according to Equation (25). However, smaller values of also

Single-Notch Filter with Infinite Rejection
Stopband rejection can be greatly improved by fine-tuning λ according to Equation (21). For parameters in Table 1, the value of λ required for infinite rejection is 0.79, well within range for the coupling tuning mechanism. Figure 7b compares the simulated and measured scattering parameters obtained in this case. A performance improvement of over 50 dB is achieved compared with that in Figure 7a, with 70 dB rejection at port 2. This improvement is at the expense of reduced isolation of port 3 with a minimum attenuation of 18 dB at resonance between ports 1 and 3. However, losses in S 21 are still extremely low.

Double-Notch Filters
A two-notch directional filter with parameter values shown in Table 1 was implemented in the same way by increasing the number of coupling resonators in each vertical branch from one to two. The magnitude of κ c determines the frequency locations of the notches in S 21 according to Equation (25). However, smaller values of κ c also increase propagation losses through the coupling resonators, reducing notch depth and increasing losses in the bandpass channel. Figure 8a compares theoretical and experimental results for κ c = 0.3. Agreement between simulation and experimental results is again excellent. Attenuation in S 21 is the same as for the single-notch filter, implying that two signals can be extracted simultaneously without hindering performance.

Double-Notch Filter with Improved Rejection
Although infinite rejection appears impossible in multinotch filters, rejection can be improved by adjusting . Figure 8b compares theoretical and measured S-parameters for a two-notch filter with = 0.3 with chosen to equalise rejection at the two notches. A significant performance improvement in a rejection of over 15 dB is obtained, once again

Double-Notch Filter with Improved Rejection
Although infinite rejection appears impossible in multinotch filters, rejection can be improved by adjusting λ. Figure 8b compares theoretical and measured S-parameters for a two-notch filter with κ c = 0.3 with λ chosen to equalise rejection at the two notches. A significant performance improvement in a rejection of over 15 dB is obtained, once again at the expense of port 3 isolation. As regards single-notch filters, lower values of λ are needed for higher values of µ. Low insertion loss is observed around the resonance frequency but increases elsewhere as matching between the guides and load deteriorates.

Applications
Toroid transformers based on soft ferrites allow magnetoinductive designs to be realised with low non-nearest-neighbour coupling. However, core materials must have low loss and high linearity. Some applications also require high power handling. In HF RFID, a modulated AC magnetic field delivers power and data to tags, which couple to the field using inductive loops. For the ISO/IEC 14443 type A standard, the carrier frequency is f c = 13.56 MHz, while tags use an internally generated subcarrier at f m = f c 16 to return data, generating sideband spectra with central frequencies, f c ± f m . Since the carrier power can be orders of magnitude greater than that of the sidebands, receiver protection is required. The single-notch filter presented here can be integrated into existing MI systems via direct connection to the waveguides, has high rejection at resonance, and may be used to transfer residual carrier to a load. However, the filter must be effective at the representative carrier powers, despite potential nonlinearities introduced by ferrite cores. To emulate generic HF RFID signals, we modulated a high-power carrier at f c using a simple cosine message waveform at f m . The waveform of the input signal s(t) is: Here A c and A m are the carrier and message amplitudes. The signal was derived from a first signal generator (Agilent E4433B) modulated by a second signal modulator (Agilent N5181A) with amplification by a 27 dB power amplifier (Minicircuits ZFDWHA-1-20+). A c was adjusted to deliver a carrier power, P c , into a 50 Ω load, while the modulation index A m /A c was chosen to generate typically small sidebands. Detection was performed using a digital storage oscilloscope (Keysight InfiniiVision DSOX3024T, Santa Rosa, CA, USA) using a 30 dB attenuator for unfiltered signals. Figure 9a shows the power spectral density (PSD) of the signal input to port 1 of the single-notch filter for P c = 1 W. Amplifier nonlinearities have resulted in minor harmonics of the AM spectrum, which have been partially suppressed at high frequencies with a low-pass filter. Figure 9b shows the PSD at port 2, the stopband channel. The output shows a carrier attenuation of 73 dB and a sideband attenuation of 1.2 dB, implying an output signal with almost intact message and negligible carrier. No significant additional harmonics are introduced. Figure 9c shows the attenuation of the two sidebands and the carrier over a power range between 0.1 and 1 W. Attenuation is almost independent of P c , implying effective power handling despite the use of ferrites. partially suppressed at high frequencies with a low-pass filter. Figure 9b shows the PSD at port 2, the stopband channel. The output shows a carrier attenuation of 73 dB and a sideband attenuation of 1.2 dB, implying an output signal with almost intact message and negligible carrier. No significant additional harmonics are introduced. Figure 9c shows the attenuation of the two sidebands and the carrier over a power range between 0.1 and 1 W. Attenuation is almost independent of , implying effective power handling despite the use of ferrites.
(a) (b) (c) Figure 9. Power spectral densities of (a) input signal to port 1 of MI single-notch directional filter; (b) output signal at port 2 normalised by the input carrier power . Dashed lines: filter bandwidth; red: harmonics at multiples of carrier frequencies; green: harmonics at message sidebands. (c) Attenuation of sidebands and carrier.

Frequency Scaling
The potential for scaling to higher frequencies is largely determined by ferrite performance. Table 2 shows performance parameters for low-cost NiZn soft ceramic ferrite materials and toroids designed for RF applications, such as EMI suppressors (upper group) and broadband transformers and baluns (lower), including their initial relative permeability , Curie temperature TC, and loss factor. The main difference between groups is the maximum usable frequency Fmax. Following Snoek's limit [52], falls as Fmax is raised. High-frequency performance is therefore accompanied by a reduction in magnetic field confinement and an increase in loss. Small toroids may be used to minimise the volume of magnetic material. However, careful design will be required to optimise performance, power handling will be reduced, and automated coil winding may be required for construction [53]. Despite this, the data suggest that response may be extended to the low UHF band. Above this, air-cored inductors may be used, but nonnearest-neighbour coupling and radiation will inevitably complicate the design and reduce performance.

Frequency Scaling
The potential for scaling to higher frequencies is largely determined by ferrite performance. Table 2 shows performance parameters for low-cost NiZn soft ceramic ferrite materials and toroids designed for RF applications, such as EMI suppressors (upper group) and broadband transformers and baluns (lower), including their initial relative permeability µ i , Curie temperature T C , and loss factor. The main difference between groups is the maximum usable frequency F max . Following Snoek's limit [52], µ i falls as F max is raised. High-frequency performance is therefore accompanied by a reduction in magnetic field confinement and an increase in loss. Small toroids may be used to minimise the volume of magnetic material. However, careful design will be required to optimise performance, power handling will be reduced, and automated coil winding may be required for construction [53]. Despite this, the data suggest that response may be extended to the low UHF band. Above this, air-cored inductors may be used, but non-nearest-neighbour coupling and radiation will inevitably complicate the design and reduce performance.

Future Work
MI directional filters can also be designed for use at UHF. Potential research could concentrate on counteracting ferrite losses or, alternatively, minimising non-nearest-neighbour coupling and radiation in open-loop designs. Furthermore, the device wave operating principles can be extended to other metamaterial types, such as elastic or acoustic [54]. In particular, it has already been shown by us that MEMS could be used to generate highperformance two-port notch and comb filters [55]. Once again, wave analysis was used to simplify the coupled dynamical equations. Extension to multiport devices can be a fitting research target. Higher performance is expected due to inherently higher Q-factors of MEMS resonators.

Conclusions
New configurations of a magnetoinductive device with directional filter properties have been introduced based purely on magnetically coupled LC resonators. Design rules have been established, and methods for calculating scattering parameters when filters are used in MI systems (which have complex, frequency-dependent impedance) or conventional systems (with real impedance) have been clarified. The circuit equations have been simplified and solved using wave analysis. Analytic expressions have been developed for S-parameters and bandwidth, a simple method for introducing tuneability has been proposed, and extensions to allow infinite rejection at a chosen frequency or multiple stopbands with high rejection have been described.
Device construction and operation have been verified at HF. Ceramic ferrite cores have been used to ensure nearest-neighbour coupling in a compact layout, and simple mechanical methods of tuning resonant frequency and coupling (essential in blocker applications) have been proposed and demonstrated. Experimental results for tuneable filters, filters with infinite rejection, and filters with multiple bandstop frequencies have all been shown to agree with theoretical models. Power handing capability sufficient for blocker applications in HF RFID has been demonstrated, with little harmonic generation due to ferrite nonlinearity. Further development will involve direct connection to a magnetoinductive antenna, and this work is in progress.