# Numerical Extraction of the Equivalent Circuit for a Basic Magnetoelectric Dipole Antenna

## Article information

## Abstract

Magnetoelectric dipoles have attracted global research attention due to its broadband, unidirectional, and high front-to-back ratio characteristics. This study implemented a co-simulation between a basic magnetoelectric dipole and its front feeding circuit through the step-by-step numerical extraction of its equivalent circuit model equipped with lumped and frequency-independent components. First, the series resonance subcircuit was derived from the series resonance point in the impedance of the magnetoelectric dipole. Second, the parallel resonance sub-circuit was achieved based on the parallel resonance point. By combining the series and parallel sub-circuits according to the sequence of their resonance frequency, the final form of the equivalent circuit for the basic magnetoelectric dipole was realized. Furthermore, to obtain the component values of the proposed circuit, a numerical fitting technique was adopted to accurately match the input impedance of the antenna and its equivalent circuit. A comparison of the circuit and antenna electromagnetic simulations showed that they agreed well with each other. Hence, the correctness and feasibility of the extraction process were verified. The overall results showed that the proposed circuit model can easily substitute for a basic magnetoelectric dipole in the implementation of antenna/circuit co-simulation in circuit simulators.

## I. Introduction

Among the many analysis methods for examining the impedance properties of antennas, the theory of equivalent circuits has been the most widely adopted [1–5]. An equivalent circuit model can substitute for a real antenna structure in conducting antenna/circuit co-simulation in system-level simulations, making the process both easy and time-saving. For the above reason, several equivalent circuit models have been proposed in the literature [6–12]. For instance, in the early days, Chu [6] examined the impedance property of an antenna from the purview of an equivalent circuit. In 1979, Prony’s method was adopted to conduct an equivalent circuit synthesis [7]. A few years later, in 1981, Pearson and Wilton [8] analyzed the physical realizability of broadband equivalent circuits. Following this, Hamid and Hamid [9] analyzed the equivalent circuit of a dipole antenna characterized by an arbitrary length. Furthermore, in 2012, an equivalent circuit for modeling the radiation impedance of tightly coupled spirals was proposed [10]. Notably, this circuit was able to describe dielectric superstrates and higher-order grating lobes. Additionally, a slot antenna fed by a microstrip was described by two short-circuited slot lines parallel with conductance [10]. In 2013, an equivalent model featuring a transmission line accurately reflected the impedance at frequencies above the first resonant mode [11]. Moreover, to accurately compute circuit element values with the least effort, Simpson analyzed a four-element circuit topology for Hertz’s dipole using the Argand diagram [12].

The magnetoelectric (ME) dipole is a unidirectional antenna primarily proposed by Chlavin [13] and developed by Luk and Wong [14]. In recent years, several novel designs for this antenna, characterized by dual polarization [15], low profile [16], wide beam [17], reconfigurable [18], and ultra-wideband [19], have been reported. Notably, the implementation of a co-simulation between the basic ME dipole and its front-feeding circuit necessitates an equivalent circuit model of the antenna. Thus far, many equivalent circuit models have been proposed for this kind of antenna [20–23]. However, these models failed to reveal the real impedance characteristics of the antenna. Therefore, to accurately reproduce the input impedance of a basic ME dipole antenna in system-level simulations, an equivalent circuit model equipped with lumped and frequency-independent components is urgently required.

The objective of this work was to extract the equivalent circuit for a basic ME dipole. For this purpose, the component values of the circuit had to be lumped and frequency independent. Although the proposed circuit could not reveal the physical principles of the antenna, it was able to reproduce the antenna’s impedance characteristics in a relative bandwidth of 85%. Therefore, it can serve as a substitute for the real antenna structure in implementing an antenna/circuit co-simulation. The equivalent circuit extraction process was organized as follows: first, the form of the equivalent circuit model was outlined, after which suitable methods for determining the initial value of the components were discussed. Following this, a numerical fitting technique was implemented. Finally, a comparison was conducted to verify the correctness and feasibility of the model.

## II. Equivalent Circuit with Frequency-Independent Elements

This section describes the process by which an equivalent circuit with frequency-independent lumped components was extracted for a basic ME dipole.

### 1. Antenna Geometry

Fig. 1 presents the structure of a basic ME dipole antenna, whose optimized sizes are listed in Table 1, with a 10-mm long coaxial feeding line. To validate the feasibility of the model, a prototype was fabricated, as shown in Fig. 1(e). The magnitude of the reflective coefficient *S*_{11} of the antenna is plotted in Fig. 2. The simulation results of the antenna achieved a relative bandwidth of 86% for |*S*_{11}| < −10 dB, which ranged from 2.2 GHz to 5.5 GHz. Notably, the simulation and measurement results agreed well with each other, although a slight discrepancy was observed owing to fabrication and measurement errors.

### 2. Extraction of the Equivalent Circuit Model

The first step of the extraction procedure involved simulating the input impedance of the basic ME dipole using the Ansys High Frequency Structural Simulator (HFSS), as plotted in Fig. 3. The operating band (2.2–5.5 GHz) comprised three series resonance points (marked by the solid green lines) and three parallel resonance points (marked by the solid blue lines). The series resonance points were located at 2.3 GHz, 3.3 GHz, and 5.1 GHz, with impedances of *Z*_{s}_{1} = 33 Ω, *Z*_{s}_{2} = 34 Ω, and *Z*_{s}_{3} = 70 Ω. Meanwhile, the parallel resonance points were located at 2.7 GHz, 4.3 GHz, and 5.5 GHz, with impedances of *Z*_{p}_{1} = 65 Ω, *Z*_{p}_{2} = 80 Ω, and *Z*_{p}_{3} = 100 Ω. Based on these series/parallel resonance points, the equivalent circuit was extracted, as depicted in the flow diagram in Fig. 4. First, the series resonance was corresponded to a series *RLC* circuit, as shown in Fig. 5(a). Second, the parallel resonance was corresponded to a parallel *RLC* circuit, as shown in Fig. 5(b). Effectively, three series circuits and three parallel circuits were obtained. Third, based on the order of the resonance frequency, the previously mentioned series resonance circuits and parallel resonance circuits were connected in sequence, as shown in Fig. 5(c), generating three series branches (*R*_{s}_{1}/*L*_{s}_{1}/*C*_{s}_{1}, *R*_{s}_{2}/*L*_{s}_{2}/*C*_{s}_{2}, and *R*_{s}_{3}/*L*_{s}_{3}/*C*_{s}_{3}) and three parallel branches (*R*_{p}_{1}/*L*_{p}_{1}/*C*_{p}_{1}, *R*_{p}_{2}/*L*_{p}_{2}/*C*_{p}_{2}, and *R*_{p}_{3}/*L*_{p}_{3}/*C*_{p}_{3}). To obtain the values of the components in Fig. 5(c), the initial values for all *R*/*L*/*C* had to be obtained in advance according to the impedance values at the corresponding resonant points. The detailed procedure to determine these initial values is discussed in the next section. Finally, a numerical fitting procedure was implemented based on the initial element values. After setting reasonable iterations and tolerance, the final equivalent circuit was realized.

### 3. Initial Value Determination

To obtain the final values of the components in Fig. 5(c), a fitting iteration process was implemented. This section demonstrates the method and the problems encountered in determining the initial values during the fitting iteration. For simplicity, only the first series resonance branch (*R*_{s}_{1}, *L*_{s}_{1}, and *C*_{s}_{1} in Fig. 5(c)) and the first parallel resonance branch (*R*_{p}_{1}, *L*_{p}_{1}, and *C*_{p}_{1} in Fig. 5(c)) are discussed in detail. Notably, the determination process for the other branches was the same as that for the first series and parallel branches.

To attain the first series resonance frequency *f*_{s}_{1} (indicated by the first solid green line in Fig. 3), it was assumed that *R*_{s}_{1} is equal to the resistance (33 Ω) at the frequency point. The first extreme point in reactance after *f*_{s}_{1} was termed *f*_{f}_{s1}, marked by the dashed green line in Fig. 3. The figure further shows that *f*_{s}_{1} and *f*_{f}_{s1} correspond to 2.3 GHz and 2.5 GHz, respectively. To obtain the *L*_{s}_{1} and *C*_{s}_{1} values for the first resonance subcircuit, two equations pertaining to *L*_{s}_{1} and *C*_{s}_{1} were formulated. The first equation was based on the assumption that the series resonance at frequency *f*_{s}_{1} is caused by *L*_{s}_{1} and *C*_{s}_{1}, while the second was based on the assumption that the total reactance value at *f*_{f}_{s1} is the sum of the reactance of *L*_{s}_{1} and *C*_{s}_{1}. Hence, the initial values for *L*_{s}_{1} and *C*_{s}_{1} were calculated using the following equations:

where *ω*_{s}_{1} and *ω*_{fs}_{1} are the angular frequencies of *f*_{s}_{1} and *f*_{f}_{s1}, respectively. Furthermore, *X*_{f}_{s1} refers to the reactance at *f*_{f}_{s1}, which is equal to 23.3 Ω, as indicated in Fig. 3 by the intersection of the first dashed green line with the black line. By computing Eq. (1), the following values were acquired: *L*_{s}_{1} = 9.28 nH and *C*_{s}_{1} = 0.52 pF. Therefore, the initial element values for the first *RLC* series branch were calculated according to the first resonant frequency *f*_{s}_{1} and the first extreme point after *f*_{s}_{1}.

Next, the initial component values for the first parallel branch were calculated according to the first parallel resonance frequency *f*_{p}_{1} and the first reactance extreme point frequency *f*_{f}_{p1} after *f*_{p}_{1}. Notably, *f*_{f}_{p1} is denoted by the dashed blue line in Fig. 3. The results showed that *f*_{p}_{1} = 2.7 GHz and *f*_{f}_{p1} = 2.9 GHz, with the parallel branch composed of *R*_{p}_{1}, *L*_{p}_{1}, and *C*_{p}_{1}. Similar to the calculations for the series resonance circuit, the initial *R*_{p}_{1} in the current estimation was equal to the resistance at *f*_{p}_{1}, as a result of which *R*_{p}_{1} = 65 Ω in this branch. Moreover, two additional equations were formulated to obtain the initial values of *L*_{p}_{1} and *C*_{p}_{1}. In contrast to the series resonance circuit, the first equation in this case was a resonance equation that accounted for the resonance frequency *f*_{p}_{1} at which a parallel resonance occurred. Meanwhile, the second equation was a susceptance equation, where the susceptance at *f*_{f}_{p1} was assumed to be contributed by *L*_{p}_{1} and *C*_{p}_{1}. Effectively, Eq. (2) were constructed to calculate *L*_{p}_{1} and *C*_{p}_{1}, as follows:

where *ω*_{p}_{1} and *ω*_{f}_{p1} are the angular frequencies at *f*_{p}_{1} and *f*_{f}_{p1}, respectively. Meanwhile, *B*_{f}_{p1} refers to the susceptance at *f*_{f}_{p1}, with a value equal to 0.0045 S. By computing Eq. (2)*L*_{p}_{1} = 1.73 nH and *C*_{p}_{1} = 2 pF were obtained.

The other series (*R*_{s}_{2}/*L*_{s}_{2}/*C*_{s}_{2} and *R*_{s}_{3}/*L*_{s}_{3}/*C*_{s}_{3}) and parallel branches (*R*_{p}_{2}/*L*_{p}_{2}/*C*_{p}_{2} and *R*_{p}_{3}/*L*_{p}_{3}/*C*_{p}_{3}) were calculated in the same way as described above. For brevity, the detailed deduction process is not presented here. After conducting a deduction process for each branch, the initial values of all components in Fig. 5(c) were obtained, as listed in Table 2.

### 4. Curve Fitting and Validation

Before curve fitting was performed, the input impedance expression for the circuit presented in Fig. 5(c) was estimated in advance. The impedance of the series branch and the admittance of the parallel branch were calculated using Eqs. (3) and (4), respectively, as follows:

Based on the above formulations, the total input impedance for the equivalent circuit depicted in Fig. 5(c) can be expressed as follows:

where *Z*_{2} can be calculated using the following equation:

To obtain a good fitting result, the objective function in the fitting algorithm was considered the difference between the input impedance at the antenna port and the equivalent circuit port, which can be expressed as:

where *n* refers to the total number of sample frequency points within the concerned frequency range and *i* indicates the *i*-th sample point. Therefore, *delta_fun* represents the algebraic sum of squares of the difference between the input impedances of the antenna and the circuit at all frequency points.

To obtain the values of the circuit components, this study adopted a numerical fitting technique using the MATLAB program in the least square sense. To achieve an accurate result, *delta_fun* was set to 0.05, while the maximum number of iterations was set to 8,000. The initial element values are listed in Table 2. Notably, the lower boundaries for *R*, *L*, and *C* were found to be zero, while the upper boundaries were 2,000 Ω, 500 nH, and 1,000 pF, respectively. The final element values were calculated by implementing the fitting process, the results of which are listed in Table 3. Based on the values in Table 3, a circuit model was built in the Advanced Design System for validation, as shown in Fig. 6. In this circuit, the source impedance was set at 50 Ω to maintain consistency with the characteristic impedance of the antenna port. The input impedance for the circuit illustrated in Fig. 6 is plotted in Fig. 7. Furthermore, the input impedance at the antenna port is depicted in Fig. 7 for comparison. It is obvious that the impedance curves for the equivalent circuit and the basic ME dipole antenna are in good agreement with each other. This agreement validates the reasonability and feasibility of the equivalent model and the robustness of the deduction procedure. Moreover, the proposed method is also applicable for extracting the equivalent circuit model for other types of antennas.

## III. Conclusion

This paper numerically extracted an equivalent circuit model for a basic ME dipole. First, the structure of the basic ME dipole and its impedance were obtained. Second, the equivalent circuit form of the antenna was extracted based on series resonance and parallel resonance. Subsequently, the proposed circuit was composed of three series and three parallel resonance sub-circuits. Third, based on the resonance points and extreme points in the reactance, the initial values of the circuit components were calculated. The final component values were obtained by accurately implementing a numerical fitting. Notably, the extraction procedure involved both model establishment and numerical calculations. Finally, the equivalent circuit simulation and antenna EM simulation were observed to be in good agreement. As a result, it is established that the proposed equivalent circuit model offers good circuit representation for basic ME dipole antennas and exhibits great potential for use in antenna/circuit co-simulation.

## References

## Biography

Zhiyi Li, https://orcid.org/0000-0003-2119-1443 was born in Sichuan, China, in 1990. He received his B.S. degree in electronics information engineering from Qiqihar University, Qiqihar, China, in 2013 and his M.E. degree in electromagnetic field and microwave technology from Xidian University, Xi’an, China, in 2016. From 2016 to 2018, he worked as an antenna designer at CEC Jinjiang INFO Industrial Co. Ltd., and as an integration and verification engineer at Ericsson in Chengdu, China. He is currently pursuing a Ph.D. in radiophysics at Sichuan University in Chengdu, China. His research interests include electromagnetic metamaterials, antenna designs, optimization algorithms, numerical simulations, matching networks, and wireless power transfer.

Yuzhu Tang, https://orcid.org/0000-0002-0516-1642 was born in Sichuan, China, in 1990. He received his bachelor’s degree in communication engineering from the University of Electronic Science and Technology of China in 2012 and his master’s degree in communication and information systems from Chengdu University of Information Engineering, China, in 2017. He is currently pursuing a Ph.D. in radiophysics at Sichuan University. His research interests include rectifiers, antennas, and functional metamaterialss.

Zhifan Zhao, https://orcid.org/0000-0001-7012-9588 was born in Ordos, Inner Mongolia, China, in 1996. He received his B.S. degree from Southwest Minzu University in 2019. He is currently pursuing a master’s degree at Sichuan University. His research interests include the design of magnetoelectric dipole antennas and TEM horn antennas.

Linwan Deng, https://orcid.org/0000-0003-0795-4811 was born in Ordos, Inner Mongolia, China, in 1996. He received his B.S. degree from Southwest Minzu University in 2019. He is currently pursuing a master’s degree at Sichuan University. His research interests include the design of magnetoelectric dipole antennas and TEM horn antennas.

Hongzheng Zeng, https://orcid.org/0000-0002-6692-6936 received his M.S. degree in electronic science and technology from Chongqing University of Posts and Telecommunications, China, in 2014, and his Ph.D. degree in radio physics from Sichuan University, China, in 2020. He is currently an assistant researcher at the Key Laboratory of Flight Technology and Flight Safety of Civil Aviation at the Civil Aviation Flight University of China. His research interests include numerical methods applied in electromagnetics and electromagnetic compatibility and interference. He is an assistant electronic design engineer at the Chinese Institute of Electronics.

Xing Chen, https://orcid.org/0000-0002-3757-0619 received his M.S. degree in radiophysics and Ph.D. degree in biomedical engineering from Sichuan University, Sichuan, China, in 1999 and 2004, respectively. In 1991, he joined the College of Electronics and Information Engineering, Sichuan University, as a teaching staff member, where he is currently a professor. His research interests include antennas, microwave imaging, global optimization, numerical methods applied in electromagnetics, and parallel computation. Dr. Chen is also a senior member of the Chinese Institute of Electronics.