Abstract
A framework for the design and optimization of large dual-linear polarized, shaped-beam reflectarrays for advanced wireless communications is presented. The methodology is based on the generalized intersection approach (IA) algorithm for both phase-only synthesis (POS) and direct optimization of the reflectarray layout, as well as on the use of a method of moments in the spectral domain assuming local periodicity. A thorough description of the design and optimization procedures is provided. To demonstrate the capabilities of the proposed framework, two examples are considered. The first example is a shaped-beam reflectarray for future 5G base stations working in the millimeter waveband, radiating a sectored-beam pattern in azimuth and squared-cosecant pattern in elevation to provide constant power in the coverage area. The second example is a very large contoured-beam reflectarray for direct-to-home (DTH) broadcasting based on real mission requirements with Southern Asia coverage.
Keywords
- array pattern synthesis
- reflectarrays
- optimization
- wireless communications
- 5G
- base station
- space communications
- shaped-beam
- contoured-beam
- generalized intersection approach
1. Introduction
Wireless communication technologies have experienced a constant and rapid development over the past few decades. This has resulted in communication systems that need to fulfill increasing tighter requirements with the goal of improving their performance and quality. In particular, future developments and integration of 5G technologies for terrestrial and space communications [1] represent a great challenge. Specifically, the antenna is an important subsystem for wireless communications, since it is the device that converts the guided waves into propagating waves in free space and vice versa. Different parameters of the antenna may be optimized depending on the application, such as size, radiation pattern, matching, etc. In many cases, a shaped-beam pattern is necessary to adequately redirect power to the desired area. For instance, direct-to-home (DTH) applications need a contoured-beam footprint to match some specific geographic area on the surface of the Earth [2]. Also, an interesting feature for base stations for wireless communications is to provide constant power over a certain angular range. This may be achieved with a shaped-beam squared-cosecant pattern [3]. Traditionally, shaped parabolic reflectors or phased arrays have been employed for these applications [2, 3]. However, shaped parabolic reflectors are bulky and expensive to manufacture, while phased arrays require complex feeding networks which introduce high losses. Nonetheless, with the popularization of the microstrip technology, reflectarray antennas have become a potential substitute to parabolic reflector dishes and phased arrays.
The concept of reflectarray antenna was first introduced in 1963 [4] as a type of antenna that combines the simplicity of reflectors and the versatility of arrays, using waveguides as the reflecting element. This resulted in a bulky and expensive structure. However, reflectarrays were not widely studied until the development of low-profile printed antennas in the 1980s, when the printed planar reflectarray was developed [5]. It consists of an array of radiating elements that are spatially fed by a primary feed, which is usually a horn antenna. Its working principle is based on altering the properties of the electromagnetic field impinging from the feed. By adjusting the dimensions of the reflectarray elements, a phase shift is introduced in the impinging field [6], allowing to obtain the desired radiation pattern.
Although designing reflectarrays for high-gain pencil beam patterns at a certain direction may be achieved with analytical equations [5], the synthesis of noncanonical beams is a challenging task and requires the use of an optimization algorithm, especially in cases with tight requirements, such as space applications [2]. Since reflectarrays are usually comprised of hundreds or even thousands of elements, the employed algorithm must be computationally efficient. Until recently, the dominant approach was the phase-only synthesis (POS) [5], which employs a simplified analysis of the unit cell. This results in an extremely efficient synthesis [7, 8, 9] but has no control over the cross-polarization performance. The first approach to the crosspolar direct optimization of reflectarray antennas was presented in [10], using a method of moments based on local periodicity (MoM-LP) for the analysis of the unit cell. However, the algorithm was slow and only handled 1 polarization and small reflectarrays (225 elements). Other approaches for the minimization of the crosspolar component of the far field include a proper arrangement of the elements [11] and the minimization of the undesired tangential field adjusting the dimensions of the element [12] or through rotation [13]. These techniques are faster, but they work at the element level and thus provide suboptimal results.
In this chapter, we present a general framework for the efficient and accurate pattern optimization of reflectarray antennas for advanced wireless communications, including copolar and crosspolar specifications. It is based on the use of the generalized intersection approach (IA) algorithm [14] for the optimization and a MoM-LP [15] for the accurate characterization of the reflectarray unit cell. The design procedure is divided in several stages. First, a phase-only synthesis (POS) is carried out, to efficiently obtain the desired copolar pattern. Then, by using a zero-finding routine and the MoM-LP, the layout of the reflectarray is obtained adjusting the dimensions of each unit cell. Finally, an optional stage to improve the cross-polarization performance may be carried out. It employs the MoM-LP directly in the optimization loop to accurately characterize the crosspolar pattern. Both the POS and direct layout optimization are carried out with the generalized IA, demonstrating the versatility of the algorithm. Two relevant examples are provided to demonstrate the capabilities of the proposed framework. First, a shaped-beam reflectarray for future 5G base stations at millimeter waveband is proposed. It radiates a sectored-beam pattern in azimuth and a squared-cosecant pattern in elevation. The second example is a very large contoured-beam, spaceborne reflectarray for direct-to-home (DTH) broadcasting based on a real space mission.
The rest of the chapter is divided as follows. Section 2 introduces the optimization framework based on the generalized IA algorithm. Section 3 describes the design and optimization methodology using the generalized IA. Section 4 contains the results regarding the two reflectarray designs for advanced wireless communications. Finally, Section 5 contains the conclusions.
2. Optimization framework for electrically large reflectarrays
2.1 Pattern requirements in the optimization procedure
Before describing in detail the optimization algorithm, we will establish the different pattern requirements that can be imposed in the optimization procedure and how they are implemented in the generalized IA. For the case of radiation pattern optimization, the requirements may be imposed in the copolar and crosspolar components. When performing a POS, only copolar requirements are considered due to the simplifications in the analysis of the unit cell [9]. However, a direct optimization of the layout may consider both copolar and crosspolar requirements. In the generalized IA, the copolar requirements are given by means of two mask templates, which impose the minimum (
where

Figure 1.
Typical requirement templates for (a) squared-cosecant pattern and (b) sectored-beam pattern.
On the other hand, there are several methodologies to implement crosspolar requirements. A typical approach is to minimize the crosspolar far field component by means of templates [16], similarly to the procedure followed with the copolar pattern. However, the crosspolar pattern does not need a lower bound in the optimization. Thus, only the maximum mask
However, there are applications in which the figure of merit for cross-polarization performance is not the crosspolar pattern. In particular, some space missions [2] give the requirements for the crosspolar discrimination (XPD) and/or crosspolar isolation (XPI).
The XPD is defined for a certain coverage zone as the difference (in logarithmic scale) point by point of the copolar gain and the crosspolar gain. Mathematically it is expressed as
where
Similarly, the XPI is defined for a certain coverage zone as the difference (in logarithmic scale) of the minimum copolar gain and the maximum crosspolar gain:
where XPI is in dB and

Figure 2.
Graphical definition of the parameters for co- and cross-polarization performance: The crosspolar discrimination (XPD), which is defined point by point as the difference between the copolar gain and the crosspolar gain, and the crosspolar isolation (XPI), which is defined for the coverage zone as the difference between the minimum copolar gain and the maximum crosspolar gain. The copolar and crosspolar patterns are in dBi, while the XPD and the XPI are in dB.
The optimization procedure should maximize the XPDmin and/or XPI. Thus, if
2.2 Generalized intersection approach
The framework for the optimization of the radiation pattern of reflectarray antennas is based on the generalized intersection approach (IA) [4]. A flowchart of the algorithm is shown in Figure 3. It is an iterative algorithm which performs two operations at each iteration

Figure 3.
Flowchart of the generalized intersection approach algorithm as applied to the far field in gain.
where
where
is the matrix of reflection coefficients which define the electromagnetic behaviour of the unit cell. These coefficients are complex numbers and are computed by a full-wave analysis tool assuming local periodicity [5].
In (8),
This operation is also applied to the crosspolar pattern when performing a direct optimization of the reflectarray layout. If the cross-polarization performance is improved by means of the XPDmin or XPI optimization, a similar expression to Eq. (11) is used but only taking into account the minimum masks, as in Eqs. (6) and (7).
The second operation of the generalized IA, denoted by
The latter operation is performed by a general minimizing algorithm [14]. In addition, as a distance definition, we employ the Euclidean norm for square-integrable functions [8], which is implemented by the weighted Euclidean metric:
where it was taken into account that the result of the forward projection is the trimmed gain in Eq. (11);
This sum can be minimized by the Levenberg-Marquardt algorithm (LMA) [9].
Finally, the generalized IA can be applied to perform a phase-only synthesis (POS), where the optimizing variables are the phase shift introduced by each reflectarray element corresponding to the phases of the direct coefficients in Eq. (10), or a direct layout optimization, where the optimizing variables are the geometrical features of the unit cell.
3. Design and optimization methodology
This section briefly describes the design methodology employing the optimization framework presented in the previous section. It is applied to a reflectarray in single-offset configuration, as shown in Figure 4. The procedure is divided into three stages: first, a phase-only synthesis to obtain the desired radiation pattern; then, a design procedure to adjust the element dimensions yielding a reflectarray layout; and the last and optional stage is the optimization of the cross-polarization performance of the reflectarray antenna.

Figure 4.
Diagram of the planar reflectarray antenna optics and the considered unit cell based on two sets of parallel and coplanar dipoles. © 2018 IEEE. Reprinted, with permission, from [
3.1 Phase-only synthesis for the copolar pattern
The first step in the design of a shaped-beam reflectarray antenna is a phase-only synthesis (POS). The aim of the POS is to obtain a phase-shift distribution that generates the desired shaped radiation pattern, which in general cannot be obtained through analytical means since that approach presents some limitations [9]. Since we are interested in dual-linear polarized reflectarrays, two phase-shift distributions are necessary, one for each linear polarization. In addition, the generalized IA is a local search algorithm. Thus, a good starting point is of utmost importance. It has been demonstrated that a properly focused pattern is sufficient for the POS [7]. In that case, the initial phase distribution for the POS may be obtained analytically [5]:
where
Then, the generalized IA is employed to synthesize the desired pattern. For the POS, the elements are modelled as ideal phase shifters, in which there are no losses (
where
3.2 Obtaining a reflectarray layout from a phase-shift distribution
The procedure to obtain a reflectarray layout from the two phase-shift distributions obtained after the POS is summarized in the flowchart of Figure 5. It requires the use of a full-wave technique based on local periodicity (FW-LP) to analyse the unit cell. Here, we employ the MoM-LP described in [18] to analyse the unit cell shown in Figure 4. In this step, a common procedure in the literature is to use a design curve obtained at normal incidence to seek the size of the reflectarray element that matches the required phase shift. However, it is recommended to consider the real angle of incidence to increase accuracy, especially for very large reflectarray antennas, since the phase shift varies with the angle of incidence [5].

Figure 5.
Flowchart of the procedure to obtain a reflectarray layout from the synthesized phase-shift distribution for two linear polarizations.
This procedure is divided into three steps. Firstly, a phase-shift table is generated, increasing the size of the element (for instance, the patch size or dipole length) in little intervals. For the case at hand and using the unit cell based on two sets of parallel dipoles of Figure 4, two variables,
3.3 Improvement of the cross-polarization performance through direct optimization
The third and final stage is optional and consists in improving the cross-polarization performance of the synthesized reflectarray by directly optimizing its layout using a FW-LP tool. This is especially important for applications with tight cross-polarization requirements, such as space missions [2], since the layout obtained in the previous stage most likely will only comply with copolar specifications. As a starting point, the layout obtained in the previous stage is employed. Also, the copolar specification masks are maintained to keep the copolar pattern within specifications while the cross-polarization performance is improved.
There are a number of approaches that can be followed in this stage depending on the application. A common approach in the literature is to directly minimize the crosspolar component of the far field [19, 20]. This is done by applying Eq. (2) for the crosspolar pattern masks in the forward projector of the generalized IA. Another approach is to impose Eqs. (6) and (7) in order to maximize the XPD
Nevertheless, this stage requires the use of a FW-LP tool to obtain the full matrix of reflection coefficients in Eq. (10) in order to correctly characterize the crosspolar radiation pattern. Thus, the improvement in cross-polarization performance will be slower than the POS in the first stage.
4. Examples of application
Here, we present two examples of application of the optimization framework presented in the previous sections. First, a medium-sized reflectarray is designed to work in a base station for future 5G application in the millimeter band at 28 GHz. The second example is a very large contoured-beam reflectarray for direct-to-home broadcasting in the Ku-band at 12.5 GHz, based on a real mission with Southern Asia coverage.
4.1 Reflectarray for 5G base station
4.1.1 Antenna specifications
For the first example, the considered reflectarray is circular and comprised of 912 unit cells (34 elements in the main axes). The periodicity is 5.36 mm in both axes, which is half a wavelength at the working frequency, 28 GHz, in order to avoid grating lobes [5]. The feed is placed at (−79.3, 0.0, 200.2) mm with regard to the centre of the reflectarray (see Figure 4), and it is modelled as a
The unit cell shown in Figure 4 is used here. The separation between dipoles is set to
The same substrate is used in both layers of the unit cell, with

Figure 6.
Unit cell study for the reflectarray for 5G base station at 28 GHz showing the phase shift (left) and the magnitude (right) for several angles of incidence. Unit cell presents a good angular stability with low losses.
Regarding the far field specifications, the chosen pattern for the 5G base station has a 30° sectored beam in azimuth and a squared-cosecant beam in elevation to provide constant power flux in an elevation span of 50°.
4.1.2 Results of the antenna design
The starting point for the POS is a pencil beam pointing at

Figure 7.
For polarization X: (a) starting phase distribution (in degrees) obtained with
The obtained layout was simulated with a MoM-LP [18], and the resulting radiation pattern for polarization X is shown in Figure 8, where the copolar and crosspolar components of the far field are shown in the

Figure 8.
Radiation pattern in the whole visible region radiated by the reflectarray designed for a 5G base station for polarization X. (a) Copolar component of the far field. (b) Crosspolar component of the far field.
On the other hand, Figure 9 represents the main cuts in elevation and azimuth for both linear polarizations along with the mask requirements. Here, it can be better appreciated how the specifications are met, with side lobes lower than −2 dB, which represent a SLL better than 20 dB for this shaped pattern. In addition, Figure 10 shows the radiation pattern in 3D perspective, along with a sketch of the reflectarray panel.

Figure 9.
Main cuts for both linear polarizations in (a) elevation and (b) azimuth with the mask requirements for the reflectarray designed for a 5G base station.

Figure 10.
3D representation of the copolar component of the radiation pattern for a 5G base station.
Regarding the cross-polarization performance, the initial design presents maximum crosspolar values of −6.1 and −6.8 dBi for polarizations X and Y, respectively, while the maximum copolar gain is 19.6 dBi for both polarizations. This gives a maximum copolar gain/maximum crosspolar gain ratio (CP
After the crosspolar optimization, the radiation pattern shown in Figure 11 was obtained for polarization X. When compared with the far field of Figure 8, it can be seen how the crosspolar pattern maximum value has been considerably reduced while keeping the copolar pattern within specifications. In fact, the maximum copolar gain is now 19.7 dBi and 19.6 dBi for polarizations X and Y, respectively. At the same time, the maximum crosspolar values are −15.7 and −15.4 dBi, with CP

Figure 11.
Radiation pattern in the whole visible region radiated by the reflectarray designed for 5G base station for polarization X after the optimization to improve the cross-polarization performance: (a) copolar component of the far field and (b) crosspolar component of the far field.
Polarization X | Polarization Y | |||||
---|---|---|---|---|---|---|
CPmax | XPmax | CPmax−XPmax | CPmax | XPmax | CPmax−XPmax | |
Initial design | 19.6 | −6.12 | 25.7 | 19.6 | −6.84 | 26.4 |
Optimized design | 19.7 | −15.75 | 35.4 | 19.6 | −15.40 | 35.0 |
Table 1.
For the reflectarray for 5G base station, summary of the performance of the initial and optimized designs regarding the maximum copolar gain (CPmax), the maximum crosspolar gain (XPmax) and the difference between them (CPmax-XPmax) for both linear polarizations.
CPmax and XPmax are in dBi, while CPmax−XPmax is in dB.

Figure 12.
Layout of the reflectarray designed for 5G base station that generates a radiation pattern with a 30° sectored beam in azimuth and a squared-cosecant beam in elevation: (a) bottom layer and (b) upper layer.
4.2 Reflectarray for direct-to-home satellite application
4.2.1 Antenna specifications
For the second example, an elliptical reflectarray with axes
A similar unit cell as in the previous example is used with different dimensions and materials. The separation between dipoles is now set to

Figure 13.
Unit cell study for the reflectarray for DTH at 12.5 GHz showing the phase shift (left) and the magnitude (right) for several angles of incidence. Unit cell presents a good angular stability with low losses while providing more than 600° of linear phase shift.
Figure 14 shows the contour requirements for the Southern Asia mission, similar to that provided by the SES-12 satellite. Zone 1 includes India, Nepal, Bhutan, Bangladesh and Sri Lanka, while zone 2 includes Pakistan and Afghanistan. According to the official specifications [22], the satellite provides an EIRP of 52 dBW for zone 1 and 48 dBW for zone 2. The EIRP can be converted into gain using the following expression:

Figure 14.
Footprint of the southern Asia coverage for direct-to-home broadcasting application. Zone 1 includes India, Nepal, Bhutan, Bangladesh and Sri Lanka, while zone 2 includes Pakistan and Afghanistan. This coverage mimics the one provided by the SES-12 satellite, placed in geostationary orbit at 95°E. (u,v) coordinates are in the satellite coordinate system.
where
4.2.2 Results of the antenna design
For the first step, a POS is carried out to obtain the desired copolar pattern in dual-linear polarization. Figure 15a shows the initial phase shift for the POS obtained with Eq. (15). It generates a focused beam in the direction

Figure 15.
For polarization X: (a) starting phase distribution (in degrees) obtained with

Figure 16.
Radiation pattern of the initial layout with southern Asia coverage for polarization X: (a) copolar pattern and (b) crosspolar pattern.
Space missions usually impose very stringent cross-polarization requirements in the form of crosspolar discrimination (XPD) and crosspolar isolation (XPI) for the transmit and receive bands, respectively. Notice that according to the definitions of minimum XPD in Eq. (4) and the XPI in Eq. (5), the XPI is a more stringent parameter than the XPD
Zone 1 | Zone 2 | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
Polarization X | Polarization Y | Polarization X | Polarization Y | |||||||||
CP | XPD | XPI | CP | XPD | XPI | CP | XPD | XPI | CP | XPD | XPI | |
Initial design | 31.52 | 32.50 | 32.08 | 31.36 | 30.51 | 29.56 | 28.63 | 32.99 | 31.13 | 28.91 | 31.76 | 29.87 |
Optimized design | 30.07 | 38.99 | 38.11 | 30.04 | 39.77 | 37.51 | 29.15 | 39.01 | 37.53 | 29.17 | 39.79 | 38.23 |
Table 2.
For the reflectarray with southern Asia coverage, comparison of the performance of the initial design after the POS and the optimized layout for improved cross-polarization performance.
CP is the minimum copolar gain in dBi in a coverage area, XPD
After the direct layout optimization, the cross-polarization performance of the reflectarray antenna significantly improved. The worst parameter is the XPI for zone 1 and polarization X, which has a value of 37.5 dB. It improved to 8 dB over the value for the initial design. The minimum improvement was 6 dB for the XPI for zone 1 and polarization X and XPD

Figure 17.
Radiation pattern of the optimized layout with southern Asia coverage for polarization X with improved cross-polarization performance: (a) copolar pattern and (b) crosspolar pattern.
5. Conclusions
A framework for the design and optimization of large dual-linear polarized reflectarray antennas has been presented. It is based on the generalized intersection approach (IA) algorithm, which is used for both a phase-only synthesis to obtain the initial design and a direct optimization of the reflectarray layout for a subsequent optimization of the cross-polarization performance. The IA employs the Levenberg-Marquardt algorithm in the backward projector and a method of moments based on local periodicity (MoM-LP) to accurately characterize the electromagnetic response of the unit cell.
In order to demonstrate the capabilities of the proposed framework, two examples for advanced wireless communication applications are provided. For the two designs, a unit cell that consists in two sets of parallel dipoles is employed. Each set of dipoles controls the phase shift of a linear polarization. The first example is a shaped-beam reflectarray for future 5G base stations, radiating a sectored beam in azimuth and a squared-cosecant beam in elevation. A circular reflectarray with a diameter of 182 mm and comprised of 912 elements was proposed. The shaped-beam reflectarray achieves a maximum gain of 19.6 dBi and a SLL better than 20 dB. In addition, the gain smoothly decreases its value in elevation a total of 15 dB over a tilt of 50°. This feature makes the synthesis of a squared-cosecant pattern a challenging task. Finally, the maximum value of the crosspolar pattern was reduced more than 8.5 dB for both linear polarizations after a direct layout optimization using MoM-LP directly in the optimization loop, while maintaining the copolar pattern within requirements.
The second example consists of a 1.1-metre reflectarray for direct-to-home (DTH) broadcasting application. A Southern Asia coverage footprint that emulates the requirements of the SES-12 satellite has been selected. This coverage presents two zones with different gain requirements: one comprising India, Nepal, Bhutan and Bangladesh with a 30 dBi of minimum copolar gain requirement and another for Pakistan and Afghanistan with a requirement of 26 dBi. After a phase-only synthesis, the layout of the reflectarray was obtained using a zero-finding routine and simulated with a MoM-LP tool. The minimum copolar gain achieved in both linear polarizations is better than 31 dBi for zone 1, while it is better than 28 dBi for zone 2. Then, a direct optimization of the layout with MoM-LP was carried out to improve the cross-polarization performance. Both the minimum crosspolar discrimination and crosspolar isolation improved at least 6 dB for both zones and linear polarizations while keeping the minimum copolar gain within requirements.
The results shown here demonstrate the versatility of the proposed framework for the design and optimization of reflectarrays, as well as the feasibility of this type of antenna for advanced wireless communications.
Acknowledgments
This work was supported in part by the Ministerio de Ciencia, Innovación y Universidades under the project TEC2017-86619-R (ARTEINE); by the Ministerio de Economía, Industria y Competitividad under the project TEC2016-75103-C2-1-R (MYRADA); by the Gobierno del Principado de Asturias/FEDER under the project GRUPIN-IDI/2018/000191; by the Gobierno del Principado de Asturias through the Programa “Clarín” de Ayudas Postdoctorales/Marie Curie COFUND under the project ACA17-09; and by Ministerio de Educación, Cultura y Deporte/Programa de Movilidad “Salvador de Madariaga” (Ref. PRX18/00424).
Thanks
The authors would like to thank Dr. R. Florencio, Prof. R. R. Boix and Prof. J. A. Encinar for providing the MoM-LP software for the analysis of the reflectarray cell.