Automatic Calibration Using a Modified Genetic Algorithm for Millimeter-Wave Antenna Modules in MIMO Systems

*is study proposes a method for designing and calibrating a millimeter-wave (mm-wave) multiple-input multiple-output (MIMO) antenna module. Herein, we adopt a design example involving a 64-element MIMO antenna array arranged in a triangular lattice (instead of the commonly used rectangular lattice) to achieve a 3°dB enhancement in effective isotropic radiated power. Analyzing a grating lobe diagram indicates a scan volume of ±60°/±45° in the azimuth/elevation direction. To calibrate the massive mm-wave MIMO antenna module, we propose a modified genetic algorithm to align the amplitude/phase of the transmitting/receiving signal of the module to reduce the time required for the calibration process. Finally, we conducted a simple experiment to validate the proposed method.


Introduction
Fifth-generation mobile communication technology is expected to introduce an advanced air interface in new radiofrequency bands. Millimeter-wave (mm-wave) frequency bands (i.e., ≥26 GHz) can provide a channel bandwidth of approximately 1 GHz (or higher) to enhance mobile broadband function with a throughput of up to 10 Gb, can support ultrareliable and low-latency communication capability with submillisecond latency under required packet loss, and can enable massive machine-type communication to connect billions of devices (i.e., Internet of ings (IoT)) [1]. Large-scale antenna arrays (also called multiple-input multiple-output (MIMO) arrays) with beam-forming technology for single-user MIMO or multiple-user MIMO communication between a base station and mobile devices can provide high antenna gain to overcome radiation propagation path losses at the aforementioned frequencies.
Moreover, sustaining an effective radio link budget can engender high-data-rate communications and high-capacity operation in networks. e effective isotropic radiated power (EIRP) is a crucial measure of transmitter performance.
Herein, we adopt a design example involving a 64-element MIMO array arranged in a triangular lattice to increase the antenna gain by 3 dB.
To equip a massive large MIMO system with beamsteering capability, the system must be adequately calibrated to generate the desired phase front to steer a dynamic threedimensional pencil beam in a specified direction. Consider, for example, MIMO antennas. After the radiating signals of the MIMO antennas are adequately calibrated and aligned uniformly across the antenna aperture, the antennas can generate a pencil beam pattern with a very high EIRP in the boresight direction; thus, the signal quality can be improved for high-data-rate applications.
Massive MIMO arrays can be calibrated through either over-the-air (OTA) field tests, including those that involve far-field [2][3][4][5] and near-field [6][7][8] test ranges, or built-in self-calibration tests [9]. Near-field techniques are advantageous for calibrating large-scale phased array systems [6] because the techniques provide useful information to antenna designers; specifically, the techniques provide information that is useful for diagnosing antenna illumination through the back-transformation of the computed far-field value to the aperture of the AUTusing an inverse fast Fourier transform. However, the signal level detected in the nearfield range is susceptible to noise interference, especially in mm-wave systems. In built-in self-calibration tests, phase toggling is applied to implement a closed-loop transmitting/ receiving coupler near the radiating element for signal detection. Although detected signals can be used for monitoring, calibration, and fault isolation in phased antenna array systems [9], signal deviations caused by parasitic circuit effects of the radiated structure, feeding network, and radome (i.e., housing of the system) cannot be considered in mm-wave systems. Moreover, the array systems must have a sufficient built-in circuit design for closed-loop testing and calibration to ensure adequate operation.
is can be achieved using radio-frequency integrated circuit (RFIC) technology for a medium-sized array system with low-power operation. However, such an arrangement is complex for an array with large antenna elements and high-power operation.
To solve the aforementioned MIMO calibration problems, this study proposes a calibration method that entails first applying a modified genetic algorithm (GA) in the farfield test range to calibrate an mm-wave massive MIMO system. e advantages of the proposed calibration method are as follows: (1) No additional closed-loop built-in circuits are required, thus reducing system design complexity (2) e overall system calibration is considered to ensure accuracy (3) e modified GA with machine learning capability based on big data can be easily applied to a mass production line for automatically calibrating 5G massive MIMO systems e remainder of this paper is organized as follows. Section 2 presents the theoretical approach of the proposed method and the numerical investigation of the simulated results. In Section 3, a simple experimental study is reported. Section 4 provides the study conclusion. Figure 1 illustrates the design concept of an mm-wave massive MIMO antenna module that can automatically undergo self-calibration through the modified GA. e MIMO antenna module comprises the following components: (1) 64-element microstrip patch antenna array; (2) 16 front-end modules, each of which has four-channel transmit/receive (TR) devices; (3) a power distribution network; and (4) a test antenna in the far-field range. is study adopted a four-element TRX half-duplex silicon front-end IC (F5280, IDT Inc.) [10], operating at 25-31 GHz. As presented in Figure 1 e far-field pattern of the simultaneously excited MIMO antenna array is expressed as follows:

System Configuration.
where where λ is the wavelength at frequency f; θ and ϕ are the azimuthal angle and angle of elevation of the far-field source, respectively; and T is the matrix transpose. ξ(θ, φ) is a function of the polarized pattern of an individual element, and w i (w i � c i e jψ i ) is the complex weighting coefficient at the ith element of weight vector W.
e Hermitian transpose, which combines transposition and conjugation, is denoted by H, and the radiation pattern of the array is denoted by f(θ, φ). To satisfy the scan volume requirement, the spacing arrangement of the antenna elements in the system should be specified, that is, either a rectangular or triangular lattice. For a rectangular lattice, the spacings in the x-axis, y-axis, and dx/dy can be expressed as follows: where θaz/θel is the scan volume specification in the azimuth/elevation direction. rough (4) and (5), the dx and dy values for elements arranged in a rectangular lattice can be derived as 0.536λ and 0.586λ, respectively, satisfying the scan coverage requirement of az � 60°and el � 45°. Moreover, for elements arranged in a triangular lattice, the dx and dy values can be derived as 1.0733λ and 0.6066λ, respectively, by analyzing a grating lobe diagram [11], as shown in Figure 1(c), where the solid and dashed lines represent the boundary of the scan volume and grating lobe, respectively. Since the aperture size is increased by arranging array elements in a triangular lattice, the MIMO provides relatively high antenna gain under the same scan volume requirement; nevertheless, such an increase is not suitable for miniaturization. On the basis of the derived spacing values, the maximum element gain (Ge � 4πdx dy/λ 2 ) corresponding to the rectangular and triangular lattice arrangements can theoretically be computed to be 5.96 and 9.13 dBi, respectively. is thus indicates that the EIRP is theoretically improved by approximately 3 dB. is improvement is crucial to MIMO design because power amplification results in poor efficiency in mm-wave systems. e use of additional power amplifiers for a 3 dB output power gain will increase power consumption levels and engender heat sink problems in such systems. According to the preceding theoretical calculations, the design of massive MIMO systems requires high EIRP in small-cell applications, which reduces volume constraints in space.

2
International Journal of Antennas and Propagation When the excitation coefficient vector W is calibrated uniformly with equal amplitude and phase distribution, the maximum antenna gain (Gmax) can be derived to be 27.19 dBi in the boresight direction (θ � 0°and ϕ � 0°). Calibration is a complex process that involves testing all states of both the phase shifter and the attenuator of each radiating element at various operational frequencies. Consider, for example, a front-end RFIC with 6-bit phase shifters and 4-bit attenuators; calibrating a 64-element MIMO antenna array with 10 frequencies would require approximately 655,360 (64 × 64 × 16 × 10 � 655,360) iterations of manual testing per each antenna element.
Accordingly, because of the lack of a built-in closed-loop calibration circuit in front-end RFICs (IDT, F5280) for precise phase and amplitude control of the excitation coefficient vectors, an artificial intelligence calibration method is essential for the automatic calibration of mass-produced massive MIMO systems.

GA for Far-Field Calibration.
On the basis of assumptions or simplifications, engineering problems can be mathematically modeled by imitating real-world principles. Subsequently, mathematical problems can be solved numerically. Finally, the derived problem solution can used to address a broader class of real problems. Several algorithms or systems have been demonstrated to be extremely helpful in addressing problems in various fields, including artificial neural networks [12,13], fuzzy control [14], simulated annealing [15,16], and GAs. Studies have successfully used GAs for antenna pattern synthesis and autocalibration [17,18] to obtain the optimum excitation coefficient vector (W opt ) for achieving high antenna performance (i.e., maximum antenna gain with a specified sidelobe level). Figure 2 presents the general design flow of a GA for the far-field calibration of a massive MIMO in an anechoic chamber; this procedure is used to mathematically imitate the evolution of the excitation coefficient vector (W) iteratively until the  Figure 1: (a) Schematic of mm-wave active massive MIMO system; (b) 8 × 8 element microstrip array (dx � 11.5 mm, dy � 6.5 mm) associated with a test antenna in the far-field range for calibration; (c) grating lobe diagram analysis (solid line: scan volume, dashed line: grating lobes) of the array with a triangular array lattice of dx � 11.5 mm and dy � 6.5 mm; (d) two-dimensional far-field patterns (E-H plane).
International Journal of Antennas and Propagation 3 optimal W opt is obtained. e stages in the iterative loop of the GA (Figure 2) are described as follows.

Evaluation.
is stage is also called initialization. It entails generating an initial population and its offspring and then ranking genes.

Selection.
is stage represents the mechanism for selecting individuals (populations) for reproduction according to their fitness (objective function) values.

Crossover.
is stage involves merging the genetic population of two individuals. Two carefully selected parents can produce excellent children.

Mutation.
is stage entails making random alterations to gene values. In real evolution, genetic information can be varied randomly by erroneous reproduction or other deformations of genes. In a GA, mutation can be realized as a random deformation of strings with a specified probability.
If the excitation coefficient vector of a MIMO system is adequately calibrated, the maximum measured EIRP would be attained in the boresight direction; this can be approximated as follows: EIRP � G A + 10 log 10 (N) + P elem + L + L scan + SP, (6) G A � 20 log 10 (|f(θ, ϕ)|), where G A is the antenna gain, P elem (13 dBm) is the transmitter power per element, N is the total number of antenna elements, SP (∼−64.91 dB) is the space loss at a 1.5 m range and 28-GHz frequency, L scan is the scan loss (0 dB at boresight), and L (∼−2 dB) is the total power loss in the antenna. When the weighting vector (W) is adequately calibrated with a uniform distribution, the maximum antenna gain Gmax is close to the antenna's area gain (10 log 10 ((4πA/λ 2 )) ∼27.19 dBi; A is the area of the antenna) in theory. us, according to (6), EIRP opt corresponding to a range R of 0 and 1.5 m can be estimated to be 56.26 and −8.65 dBm, respectively. During calibration, the measured EIRP value should be less than the EIRP opt value because of random amplitude/ phase errors of W caused by factors such as device variation, unequal transmission line trace lengths, and radome effects.
us, the cost function (f (W i )) at the ith iteration of the GA optimization loop can be defined as follows: where EIRP (i) represents the EIRP measured in the ith iteration of the GA optimization loop for a test antenna in the boresight direction in the far-field range (R � 1.5 m in this case). When (2) is substituted into (6) and (7), cost function (8) in a root-mean-square sense can be summarized as follows for numerical analysis: Accordingly, the GA design flow (Figure 2) can be applied to calibrate the massive MIMO system in order to achieve maximum EIRP and tune the phase and amplitude of W until the optimal excitation coefficient (W opt ) is attained. e error of the excitation weighting vector (δW � W opt − W) can be computed and written in the internal memory of the front-end ICs to correct the amplitude/phase error of the excitation weighting vectors. Phase errors engendered by the deviation of signal traces connecting from chips to antennas can be considered the determinant errors, but signal errors caused by chip-to-chip variations can be considered random errors, which is the major consideration in the following numerical analysis. e phase/amplitude of W is assumed to be randomly distributed uniformly around the interval (0, 2π)/(0.05, 0.2). e excitation coefficient of each radiating element (W i , i � 1, 2, . . ., 64) should be calibrated to the normalized value of 0.125 (1/ �� N √ ) in order to achieve maximum EIRP or antenna gain. Before the GA iterative loop (Figure 2) is executed for calibration, the excitation coefficient (W i , i � 1, 2, . . ., 64) of each element to the binary code should be encoded to represent each population of the generation. In this study, a 6-bit phase shifter associated with three LSBs of the attenuator was used for coding W.
In each population, the phase/amplitude error of the weighting vector (W i ) is also generated uniformly in the interval (0, 2π)/(0.05, 0.2) to compute the EIRP (or antenna gain) of each generation using (6) and (8) to rank genes. Subsequently, the inferior genes (poor EIRP) are discarded and superior genes (better EIRP) are crossed over to produce the next generation. Random mutation ensures the evolution of the generation to successfully attain the specified goal. After the GA iterative loop is executed, the massive MIMO system can be calibrated automatically to attain the maximum EIRP and W opt simultaneously. Figure 3 illustrates numerical simulation results regarding the autocalibration of the 64-element MIMO system through a GA with a population (M) of 20000 and 600 iterative loops (I c ). Figure 3(a) presents plots from a convergence analysis with respect to the iterative number, indicating the amplitude/phase standard deviation error of W to be less than 0.5% and antenna gain to be higher than 27.18 dBi after 300 iterative loops. In Figure 3(a), the solid line, dashed line, and dashed line with circles represent the simulated antenna gain, amplitude error (%), and phase error (%), respectively. Figure 3(b) shows a plot of the computed EIRP against the observation angle (θ) under W before and after calibration, revealing that the computed EIRP can be improved from −30 to −8.6 dBm in the boresight direction after the GA-based calibration. e corresponding phase and normalized magnitude of W used in the EIRP computation (Figure 3(b)) are plotted in Figure 3(c). e marginal deviation in the phase/amplitude of the excitation coefficient can be attributed to the quantization error.

Modified GA for Far-Field Calibration.
Generally, information regarding the randomness of the error of the excitation weighting vector is not available. In the previous section, to ensure the accuracy of the calibration results, we assume that the phase/amplitude of the error of weighting vectors is uniformly distributed in the interval (0, 2π)/(0.05, 0.2); moreover, we use a large population and high number of iterations in the GA optimization loops. is is a timeconsuming process and is unsuitable for mass production. If the randomness of the error of the vector can be adequately characterized and set in advance for the GA iterative loop, the calibration time can be reduced considerably.
Accordingly, the modified GA proposed in this study solves the aforementioned problem of high time complexity. In the modified GA, a large historical database with a reward mechanism is established for machine learning of the characteristics of random errors. erefore, in mass production, the random behavior of each population can be estimated according to the parameters established in the historical database such that the calibration time can be reduced considerably. Figure 4 displays the flowchart of the modified GA. To assess the feasibility of the modified GA, we can use a preset Gaussian error of the vector (δW s ) with an amplitude (v/v) and phase (rad) mean/standard deviation of m � 0/σ � 0.05 and m � 0/σ � 1.57, respectively.
We intentionally ignore this background information and set the population of random genes (M) with a uniform distribution of the phase/amplitude of the error of the vector in (0, 2π)/(0, 0.2) to 20000; moreover, we set number of iterative loops (I c ) to 500. We execute the GA-based calibration process 60 times (NL � 60). When the specified goal is achieved such that W opt (i) (i � 1,. . ., NL), the error of the vector can be simply computed using formula δ W (i) � W opt (i) − W (i) to produce 3840 (NL × N � 64 × 60) random amplitude/phase error numbers. For 3840 random variables, we can fit the probability distribution function of this event after running 60 GA calibration procedures. e simulation of the historical database after the 60 calibration runs indicates that the mean (µ) and standard deviation (σ) of the amplitude/phase of the error vector (δW) with a Gaussian distribution are −0.0097/−0.00117 and 0.048/0.151, respectively. Hence, we adaptively reset the estimation parameters (μ and σ) for use in the iterative loops of our modified GA (Figure 4) to reduce the calibration time.
e convergence analysis results obtained from the simulation of our GA are plotted in Figure 5 and compared with those obtained using the conventional GA. e solid lines with symbols represent the simulated results for our modified GA and the dashed lines with solid symbols represent the results for the conventional GA. Figure 5 reveals that our modified GA reduces the required number of iterative loops by 75 (from 275 to 200 loops), representing a nearly 27% reduction, compared with the conventional GA. A 3-byte SPI write operation for the transmission of write commands would require 1200 ns on the basis of a 20-MHz clock rate; this would result in a latency of at least 76.8 ms for one beam-steering command as well as additional computation time.
erefore, the proposed GA would require at least 15.36 secs (76.8 ms × 200) to execute calibration processes.

Experimental Study
Because of budget limitations, we conducted a simple experimental study (Figure 6(a)) using a 2 × 2 subarray antenna module of the massive MIMO system (Figure 1(b)). In the experiment, four microstrip antennas were printed on the top-side and the front-end IC (IDT5280) associated with one I/O and one power connector mounted on the back side of the PCB.
Considering the lack of the dedicated testing software to command the front-end IC, we performed a bench test (Figure 6(b)) by using IDT's testing software (Timing Commander) to generate the database for calibration through the proposed method. e database included the amplitude/phase information of 64-phase shifter command states (6 bits) changing corresponding to each element at operational frequencies. On the basis of the bench test results, the signal vector error (W) engendered by changes in the phase command states was simulated to calibrate the phase distribution of the 2 × 2 MIMO system using the proposed method (Section 2.2). We performed the calibration at a frequency of 27 GHz only because the chipset vendor already provided the default setting of 28 GHz. To simplify the experiment, we ensured that all channels constantly radiated at the maximum power gain level. Table 1 lists the calibration data of phase shifter command states corresponding to the four-element system at 27 and 28 GHz. Running a bench test would require 2.5 h; nevertheless, our proposed calibration method could shorten the duration of this time-consuming process to a few seconds subject to the availability of dedicated testing software. Subsequently, as depicted in Figure 7, we applied the default setting of the front-end IC (f � 28 GHz) for the experimental International Journal of Antennas and Propagation study by conducting OTA measurements at a compact antenna test range at Oriental Institute of Technology. Figure 8 displays plots of the measured copolarization radiation patterns at 27/28 GHz, indicating a measured peak gain of 44.2/48.65 dBi. e measured patterns indicate the simultaneous achievement of both maximum peak gain and adequately calibrated phase distributions, as illustrated by the symmetrical radiation patterns in the boresight direction ( Figure 9). We subsequently applied the data calibrated by the proposed method (Table 1 at 27 GHz) to the F5280 IC register to perform OTA measurements again. e measurement results demonstrated that the gain could be improved by 3.3 dB at 27 GHz but degraded by 2.85 dB at 28 GHz, indicating that excitation coefficient calibration against operating frequencies is essential for MIMO system performance. Figure 9 illustrates plots of the measured EIRP and TRP against the input power level, with the power being swept from 0 to 10 dBm at VNA and being attenuated by approximately 34.5 dB cable loss. Finally, we changed the phase shifter states to be out of phase along the x-axis to form difference patterns. e measurement results ( Figure 10)       International Journal of Antennas and Propagation confirmed the beam-steering capability of the MIMO antenna module that was adequately calibrated using the proposed method.

Conclusion
e deployment of the mm-wave massive MIMO in the form of a small cell or customer-provided equipment is essential for 5G applications. Reliable and time-efficient calibration and test methods are crucial for developing mmwave massive MIMO systems. In this study, we developed an 8 × 8 massive MIMO system with a triangular lattice to achieve a theoretical maximum EIRP of 56.26 dBm. We first propose a modified GA for the far-field calibration of the mm-wave massive MIMO system. We then present a numerical analysis to validate the effectiveness of the proposed GA. Finally, we executed a simple experiment to validate that a well-calibrated phase distribution can enable the achievement of maximum array gain in the boresight direction. erefore, our proposed method can enable automatic calibration of massive MIMO systems [19].
Data Availability e figures' data used to support the findings of this study are currently under embargo while the research findings are commercialized. Requests for data, 12 months after publication of this article, will be considered by the corresponding author.

Conflicts of Interest
e authors declare that they have no conflicts of interest.