Multibeam Characteristics of a Negative Refractive Index Shaped Lens

Narrow beam width, higher gain and multibeam characteristics are demanded in 5G technology. Array antennas that are utilized in the existing mobile base stations have many drawbacks when operating at upper 5G frequency bands. For example, due to the high frequency operation, the antenna elements become smaller and thus, in order to provide higher gain, more antenna elements and arrays are required, which will cause the feeding network design to be more complex. The lens antenna is one of the potential candidates to replace the current structure in mobile base station. Therefore, a negative refractive index shaped lens is proposed to provide high gain and narrow beamwidth using energy conservation and Abbe’s sine principle. The aim of this study is to investigate the multibeam characteristics of a negative refractive index shaped lens in mobile base station applications. In this paper, the feed positions for the multibeam are selected on the circle from the center of the lens and the accuracy of the feed position is validated through Electromagnetic (EM) simulation. Based on the analysis performed in this study, a negative refractive index shaped lens with a smaller radius and slender lens than the conventional lens is designed, with the additional capability of performing wide-angle beam scanning.


Introduction
In recent years, there has been a rapid development of antennas to meet the advanced mobile technology requirements. In the 5G mobile technology system, millimeter wave is used, which produces a smaller cell size and requires a multibeam radiation pattern in order to achieve massive Multiple Input Multiple Output (MIMO) operation [1,2] with the purpose to serve massive number of users, consistent interconnectivity and larger capacity. At millimeter wave, the base station size becomes less than 30cm. So, aperture antennas such as reflector antennas and dielectric lens antennas become promising candidates as compared to the present array antennas. For multibeam operation, the dielectric lens antenna can achieve very good performance by designing the lens using Abbe's sine condition [3], which is validated further in this paper.
There were several studies conducted to investigate the multibeam application using the lens antenna. The Luneburg lens antenna is one of the commonly used designs. However, major issues of the conventional Luneburg lens are large and heavy. In [4], a Luneburg lens with an operating frequency of 1-8 GHz was proposed. In the paper, multilayer Luneburg lenses were fabricated with the diameters of 0.6 and 1.0 m and weighing of 5 and 21 kg, respectively. A number of feed antennas was placed around the lens to produce multibeam with the same beam shapes. Beam scanning was achieved by switching the feed around the surface of the lens. The traditional Luneburg lenses weight has been reduced up to 8-10 times; however, in 5G base station application, this weight and size are not practical. A flat Luneburg lens antenna with a geometry of 80 × 5 mm 2 was designed to operate at Ka-band [5]. A linear antenna array with 11 E-shaped patch antenna elements is employed to feed the lens antenna. The lens is formed by rectangular shaped concentric six layered dielectric materials with a permittivity ranging from 1.74 to 10 to achieve higher aperture efficiency. This lens antenna produces a realized gain of 14.2 dBi. In total, 11 beams are produced for ±60 • beam scanning. However, a lens antenna employing natural dielectric material is known to produce a thicker lens size, and thus is not suitable for base station application. In order to reduce the lens thickness, a metamaterial lens concept was proposed [6]. Lens structure became concave and reduced the lens thickness. As for metamaterial lens fabrication, a lens of refractive index of n = −1 was shown in [7]. The structure is obtained by using 100 unit cells of dielectric resonator structure, 0.5 of f/D with lens aperture of 5.7 λ (156 mm) at an 11 GHz operating frequency. The achieved multibeam radiation pattern is shown in [8]. To achieve a good multibeam characteristic, Abbe's sine condition design is proposed in this paper, and comparison with the energy conservation technique is analyzed.
In this paper, the application of Abbe's sine design and energy conservation design are employed to develop a metamaterial lens. Based on the lens shaping equation, the possibility of the negative refractive index design is ensured. In order to solve the design equations, a MATLAB program is developed. The design accuracies of the MATLAB program are estimated using ray tracing results and calculation of the refraction angles on the lens surfaces. For validation, the multibeam radiation pattern is calculated using a commercial electromagnetic simulator, High Frequency Structure Simulator (HFSS). Figure 1 shows the proposed base station structure.
Sensors 2020, 20, x FOR PEER REVIEW 2 of 18 achieved by switching the feed around the surface of the lens. The traditional Luneburg lenses weight has been reduced up to 8-10 times; however, in 5G base station application, this weight and size are not practical. A flat Luneburg lens antenna with a geometry of 80 × 5 mm 2 was designed to operate at Ka-band [5]. A linear antenna array with 11 E-shaped patch antenna elements is employed to feed the lens antenna. The lens is formed by rectangular shaped concentric six layered dielectric materials with a permittivity ranging from 1.74 to 10 to achieve higher aperture efficiency. This lens antenna produces a realized gain of 14.2 dBi. In total, 11 beams are produced for ±60° beam scanning. However, a lens antenna employing natural dielectric material is known to produce a thicker lens size, and thus is not suitable for base station application. In order to reduce the lens thickness, a metamaterial lens concept was proposed [6]. Lens structure became concave and reduced the lens thickness. As for metamaterial lens fabrication, a lens of refractive index of n = −1 was shown in [7]. The structure is obtained by using 100 unit cells of dielectric resonator structure, 0.5 of f/D with lens aperture of 5.7 λ (156 mm) at an 11 GHz operating frequency. The achieved multibeam radiation pattern is shown in [8]. To achieve a good multibeam characteristic, Abbe's sine condition design is proposed in this paper, and comparison with the energy conservation technique is analyzed. In this paper, the application of Abbe's sine design and energy conservation design are employed to develop a metamaterial lens. Based on the lens shaping equation, the possibility of the negative refractive index design is ensured. In order to solve the design equations, a MATLAB program is developed. The design accuracies of the MATLAB program are estimated using ray tracing results and calculation of the refraction angles on the lens surfaces. For validation, the multibeam radiation pattern is calculated using a commercial electromagnetic simulator, High Frequency Structure Simulator (HFSS). Figure 1 shows the proposed base station structure.

Design Method
In designing the lens surface, the ray tracing method was employed. All rays passing through the aperture plane are designed to be parallel to the z-axis in order to achieve a flat wave front at the aperture plane. Figure 2 shows the proposed antenna configuration and the associated radiation parameters. The feed radiator radiates signal towards the negative refractive index shaped lens, and the result is observed at the aperture. The detailed description of the parameters is shown in Table 1.

Design Method
In designing the lens surface, the ray tracing method was employed. All rays passing through the aperture plane are designed to be parallel to the z-axis in order to achieve a flat wave front at the aperture plane.
Sensors 2020, 20, 5703 3 of 17 Figure 2 shows the proposed antenna configuration and the associated radiation parameters. The feed radiator radiates signal towards the negative refractive index shaped lens, and the result is observed at the aperture. The detailed description of the parameters is shown in Table 1.   Figure 2, at surface 1, S1 of the lens, Snell's law is given by Equation (1) [9].

Lens Configuration
At the lens surface 2, the expression for the slope can be derived from the condition that all exit rays after refraction are parallel to the z-axis, as shown in Equation (2). The expression can be separated into and , as shown in Equation (2) for variable change from dx to .
The equal condition of ray path length can be expressed in Equation (3) where indicates the total path length from the feed to the aperture plane.   Figure 2, at surface 1, S 1 of the lens, Snell's law is given by Equation (1) [9].
At the lens surface 2, the expression for the slope dz dx can be derived from the condition that all exit rays after refraction are parallel to the z-axis, as shown in Equation (2). The dz dx expression can be separated into dz dθ and dx dθ , as shown in Equation (2) for variable change from dx to dθ.
Sensors 2020, 20, 5703 4 of 17 The equal condition of ray path length can be expressed in Equation (3) where Lt indicates the total path length from the feed to the aperture plane.
By using this equation, the variable ϕ in Equations (1) and (2) can be expressed by the variable θ. Then, through simplification, the variable of Equations (1) and (2) becomes only θ. The electric power conservation at the ray is composed of dx and dθ The three differential Equations (1), (2) and (4) determine the lens shape in the MATLAB program. For the electric field intensity at feed radiator, Ep 2 (θ), Equation (5) is implemented. Meanwhile, Equation (6) is used for aperture illumination distribution, Ed 2 (x).

Abbe's Sine Condition
In optic, a collimating lens can be designed to be a coma free for a limited scan by imposing the Abbe's sine condition. Coma refers to the aberration inherent to certain optical designs or due to imperfection in the lens or other components that results in off-axis point sources. This condition is automatically fulfilled if the inner surface of a conventional waveguide lens is spherical. The condition in the red circle in Figure 3 interprets the Abbe's sine condition. When the initial and the final ray are extended, these rays are intersecting inside the lens on a circle radius of f e [2]. Abbe's sine law in Equations (7) and (9) are applied in the MATLAB program to design the shaped lens.
Sensors 2020, 20, x FOR PEER REVIEW 4 of 18 By using this equation, the variable in Equations (1) and (2) can be expressed by the variable θ. Then, through simplification, the variable of Equations (1) and (2) becomes only θ. The electric power conservation at the ray is composed of and The three differential Equations (1), (2) and (4) determine the lens shape in the MATLAB program. For the electric field intensity at feed radiator, ( ), Equation (5) is implemented.

Abbe's Sine Condition
In optic, a collimating lens can be designed to be a coma free for a limited scan by imposing the Abbe's sine condition. Coma refers to the aberration inherent to certain optical designs or due to imperfection in the lens or other components that results in off-axis point sources. This condition is automatically fulfilled if the inner surface of a conventional waveguide lens is spherical. The condition in the red circle in Figure 3 interprets the Abbe's sine condition. When the initial and the final ray are extended, these rays are intersecting inside the lens on a circle radius of fe [2]. Abbe's sine law in Equations (7) and (9) are applied in the MATLAB program to design the shaped lens.

Program Flow Chart
A flow chart in Figure 4 represents the processes involved in negative refractive index lens design using MATLAB software. Program codes were developed based on the equations and formula that were previously described in Sections 2.2.1 and 2.2.2. The initial parameters (n, θ m , r o , d o ) determined at the lens edge are shown in Figure 5. The equations for feed radiation pattern, Ep 2 (θ), and aperture distributions, Ed 2 (x), are given. Next, differential Equations (1), (2) and (4) or (1), (2) and (9) are solved by the MATLAB routine of "ode45".
Sensors 2020, 20, x FOR PEER REVIEW 5 of 18 determined at the lens edge are shown in Figure 5. The equations for feed radiation pattern, ( ), and aperture distributions, ( ), are given. Next, differential Equations (1), (2) and (4) or (1), (2) and (9) are solved by the MATLAB routine of "ode45". The explanation for the initial input values is shown in Figure 5.
is an important parameter that determines the lens thickness where d is the initial thickness.  determined at the lens edge are shown in Figure 5. The equations for feed radiation pattern, ( ), and aperture distributions, ( ), are given. Next, differential Equations (1), (2) and (4) or (1), (2) and (9) are solved by the MATLAB routine of "ode45". The explanation for the initial input values is shown in Figure 5.

Simulation of Energy Conservation Law
is an important parameter that determines the lens thickness where d is the initial thickness.

Simulation of Energy Conservation Law
A negative refractive index lens is designed based on the energy conservation law by using Equations (3) and (4) in the MATLAB shape design program. The explanation for the initial input values is shown in Figure 5. θ 0 is an important parameter that determines the lens thickness where d 0 is the initial thickness.

Simulation of Energy Conservation Law
A negative refractive index lens is designed based on the energy conservation law by using Equations (3) and (4) in the MATLAB shape design program. Figure 6a shows the horn feed radiation, Ep 2 (θ). The lens area is about −55 • to 55 • at −10 dB. It means that the area that is covered by the horn radiation beam is from −55 • to 55 • . The aperture distribution, Ed 2 (x), is shown in Figure 6b, which represents the edge level of the aperture distribution.

Aperture Distribution Design
Sensors 2020, 20, x FOR PEER REVIEW 6 of 18 3.1.1. Aperture Distribution Design Figure 6a shows the horn feed radiation, ( ). The lens area is about −55° to 55° at −10 dB. It means that the area that is covered by the horn radiation beam is from −55° to 55°. The aperture distribution, ( ), is shown in Figure 6b, which represents the edge level of the aperture distribution. Considering all the design parameters stated previously, the shaped lens was obtained as illustrated in Figure 8. For the designed lens in Figure 8, the accuracy of the MATLAB program was validated by calculating the incident angle, and refracted angle, at each surface of the lens. The refractive index value, , was manually determined using the Snell's Law in (10). The calculated value is expected to be similar to the value ( =−√2), which is used in the MATLAB design program. The lens parameters are tabulated in Table 2.

=
(11) Aperture distribution can be expressed as in (10). When Equation (10) is compared with theoretical Equation (6) in Figure 7, it is observed that both equations have a good agreement with each other. Considering all the design parameters stated previously, the shaped lens was obtained as illustrated in Figure 8. For the designed lens in Figure 8, the accuracy of the MATLAB program was validated by calculating the incident angle, and refracted angle, at each surface of the lens. The refractive index value, , was manually determined using the Snell's Law in (10). The calculated value is expected to be similar to the value ( =−√2), which is used in the MATLAB design program. The lens parameters are tabulated in Table 2. Considering all the design parameters stated previously, the shaped lens was obtained as illustrated in Figure 8. For the designed lens in Figure 8, the accuracy of the MATLAB program was validated by calculating the incident angle, θ i and refracted angle, θ r at each surface of the lens. The refractive index value, n, was manually determined using the Snell's Law in (10). The calculated n value is expected to be similar to the n value (n = − √ 2), which is used in the MATLAB design program. The lens parameters are tabulated in Table 2.
Sensors 2020, 20, x; doi: FOR PEER REVIEW www.mdpi.com/journal/sensors validated by calculating the incident angle, and refracted angle, at each surface of the lens. The refractive index value, , was manually determined using the Snell's Law in (10). The calculated value is expected to be similar to the value ( =−√2), which is used in the MATLAB design program. The lens parameters are tabulated in Table 2.   Table 3 represents the calculated values considering the four first rays coming from the shaped lens in Figure 9a,b. The 1 and 2 values were calculated based on Equation (11) and compared to the exact value that is used in MATLAB which is = −√2. It is noticed that there is a small difference between the calculated values due to some errors that occurred during angle determination.

Parameters Value
Focal length (mm) 46.56 Table 3 represents the calculated n values considering the four first rays coming from the shaped lens in Figure 9a,b. The n 1 and n 2 values were calculated based on Equation (11) and compared to the exact n value that is used in MATLAB which is n = − √ 2. It is noticed that there is a small difference between the calculated values due to some errors that occurred during angle determination. Table 3. Calculated n-value.    Table 3 represents the calculated values considering the four first rays coming from the shaped lens in Figure 9a,b. The and values were calculated based on Equation (11) and compared to the exact value that is used in MATLAB which is = −√2. It is noticed that there is a small difference between the calculated values due to some errors that occurred during angle determination.

Abbe's Sine Condition
The same process was repeated for Abbe's sine Law by using Equations (7) and (9) in the MATLAB design program. The shaped lens is shown in Figure 10. The lens parameters are tabulated in Table 4.

Abbe's Sine Condition
The same process was repeated for Abbe's sine Law by using Equations (7) and (9) in the MATLAB design program. The shaped lens is shown in Figure 10. The lens parameters are tabulated in Table 4.     Table 5 represents the calculated n value considering the four first rays from the shaped lens in Figure 11a,b. A small difference between the calculated values and the exact value used in MATLAB may be due to errors during angle determination.

Abbe's Sine Condition
The same process was repeated for Abbe's sine Law by using Equations (7) and (9) in the MATLAB design program. The shaped lens is shown in Figure 10. The lens parameters are tabulated in Table 4.   Table 5 represents the calculated value considering the four first rays from the shaped lens in Figure 11a,b. A small difference between the calculated values and the exact value used in MATLAB may be due to errors during angle determination. Figure 11. Rays in and rays out. Figure 11. Rays in and rays out.

Horn Feed
In order to calculate a negative refractive index material, the HFSS Simulator was used. In MATLAB, a point source was used to represent the feed, but in HFSS, a horn antenna was used to represent the feed radiator. Table 6 shows the simulation parameters using HFSS. Table 6. Simulation parameters using High Frequency Structure Simulator (HFSS).

Parameters Description/Value
Boundary Condition Radiation Boundary Refractive index, −1.4142 ( A conical horn antenna is designed to operate at 28 GHz as shown in Figure 12a with the radiation pattern illustrated in Figure 12b.

Horn Feed
In order to calculate a negative refractive index material, the HFSS Simulator was used. In MATLAB, a point source was used to represent the feed, but in HFSS, a horn antenna was used to represent the feed radiator. Table 6 shows the simulation parameters using HFSS. A conical horn antenna is designed to operate at 28 GHz as shown in Figure 12a with the radiation pattern illustrated in Figure 12b.

On-Focus
The designed shaped lens structure was simulated for performance evaluation. The lens structure has a permittivity of = −2 and permeability = −1. Figure 13a,b show the electric field distribution and the antenna gain, respectively. From Figure 13a, it is observed that the lens transformed the spherical wave into the plane wave. The obtained parameters are shown in Table 7. The theoretical value shows the uniform aperture case. This lens antenna structure produced a 27.55 dB gain with a beamwidth of 7.81°. The efficiency of the structure is about 66%.

On-Focus
The designed shaped lens structure was simulated for performance evaluation. The lens structure has a permittivity of ε r = −2 and permeability µ r = −1. Figure 13a,b show the electric field distribution and the antenna gain, respectively. From Figure 13a, it is observed that the lens transformed the spherical wave into the plane wave. The obtained parameters are shown in Table 7. The theoretical value shows the uniform aperture case. This lens antenna structure produced a 27.55 dB gain with a beamwidth of 7.81 • . The efficiency of the structure is about 66%.

Off-focus
In order to prove that the designed lens shape is suitable for wide-angle beam scanning, multibeam characteristics are investigated in this paper. The feed horn position shifts spherically by 10° using the formula of = , where is the distance between the feed position and the center of the lens, while is the focal length. The off-focus position is illustrated in Figure 14. The results of antenna gain for different feed positions is illustrated in Figure 15 and Table 8. It is observed from the graph that there is no significant gain reduction for all the scanning angles. The scanning losses are all less than 3 dB, which indicates that beamwidths are consistent for all scanning angles. It can be concluded that this structure is suitable for wide-angle beam scanning.

Off-Focus
In order to prove that the designed lens shape is suitable for wide-angle beam scanning, multibeam characteristics are investigated in this paper. The feed horn position shifts spherically by 10 • using the formula of F = f s, where F is the distance between the feed position and the center of the lens, while f s is the focal length. The off-focus position is illustrated in Figure 14.

Off-focus
In order to prove that the designed lens shape is suitable for wide-angle beam scanning, multibeam characteristics are investigated in this paper. The feed horn position shifts spherically by 10° using the formula of = , where is the distance between the feed position and the center of the lens, while is the focal length. The off-focus position is illustrated in Figure 14. The results of antenna gain for different feed positions is illustrated in Figure 15 and Table 8. It is observed from the graph that there is no significant gain reduction for all the scanning angles. The scanning losses are all less than 3 dB, which indicates that beamwidths are consistent for all scanning angles. It can be concluded that this structure is suitable for wide-angle beam scanning. The results of antenna gain for different feed positions is illustrated in Figure 15 and Table 8. It is observed from the graph that there is no significant gain reduction for all the scanning angles. The scanning losses are all less than 3 dB, which indicates that beamwidths are consistent for all scanning angles. It can be concluded that this structure is suitable for wide-angle beam scanning.    Figures 16 and 17 show the results of the electric distribution in magnitude and phase, respectively, performed through near-field calculation. For on-focus condition, both intensity and phase distributions show uniform and symmetrical characteristics as shown in Figures 16a and 17a, respectively. In off-focus condition, the uniform characteristics are expected to be disturbed by some distortions due to the non-linearity in phase; thus, the behaviors are non-symmetrical, as observed in Figures 16b and 17b. Both electrical intensity and phase are slightly distorted when the feed is displaced from its original position.  Figures 16a and 17a, respectively. In off-focus condition, the uniform characteristics are expected to be disturbed by some distortions due to the non-linearity in phase; thus, the behaviors are non-symmetrical, as observed in Figures 16b and  17b. Both electrical intensity and phase are slightly distorted when the feed is displaced from its original position.   Figures 16 and 17 show the results of the electric distribution in magnitude and phase, respectively, performed through near-field calculation. For on-focus condition, both intensity and phase distributions show uniform and symmetrical characteristics as shown in Figures 16a and 17a, respectively. In off-focus condition, the uniform characteristics are expected to be disturbed by some distortions due to the non-linearity in phase; thus, the behaviors are non-symmetrical, as observed in Figures 16b and 17b. Both electrical intensity and phase are slightly distorted when the feed is displaced from its original position.

On-Focus
For Abbe's sine condition, a shaped lens structure with the permittivity value of ε r = −2 and permeability value of µ r = −1 was simulated as shown in Figure 18a. The performance of the designed lens is illustrated in Figure 18b. It is observed that this Abbe's sine shaped lens produced a gain of 27.48 dB. The beamwidth is 7.64 • . This lens also has 65% efficiency, as shown in Table 9.

On-focus
For Abbe's sine condition, a shaped lens structure with the permittivity value of = −2 and permeability value of = −1 was simulated as shown in Figure 18a. The performance of the designed lens is illustrated in Figure 18b. It is observed that this Abbe's sine shaped lens produced a gain of 27.48 dB. The beamwidth is 7.64°. This lens also has 65% efficiency, as shown in Table 9.

Off-focus
For multibeam investigation, the feed position of the horn antenna is varied from 0° to 40° with 10° intervals as illustrated in Figure 19. The same condition in energy conservation ( = ) was applied for Abbe's sine shaped lens. The results for all scanning angles are represented in Figure 20 and Table 10. The gain is slightly lower as compared to the energy conservation shaped lens.

Off-Focus
For multibeam investigation, the feed position of the horn antenna is varied from 0 • to 40 • with 10 • intervals as illustrated in Figure 19. The same condition in energy conservation (F = f s) was applied for Abbe's sine shaped lens. The results for all scanning angles are represented in Figure 20 and Table 10. The gain is slightly lower as compared to the energy conservation shaped lens. However, there is no significant gain reduction for every scanning angle. The beamwidths are also consistent for all scanning angles.
Sensors 2020, 20, x FOR PEER REVIEW 13 of 18 However, there is no significant gain reduction for every scanning angle. The beamwidths are also consistent for all scanning angles.           Figures 22a and 23a, respectively. Figures 22b and 23b present off-focus condition, which shows that both electrical intensity and phase are slightly distorted. The feed displacement from its original position affects the electrical intensity and the phase distribution linearity. Figures 22 and 23 show the results of the electric distribution in magnitude and phase, respectively, for Abbe's sine lens. For on-focus condition, both intensity and phase distributions exhibited uniform characteristics as shown in Figure 22a and Figure23a, respectively. Figure 22b and Figure 23b present off-focus condition, which shows that both electrical intensity and phase are slightly distorted. The feed displacement from its original position affects the electrical intensity and the phase distribution linearity.   exhibited uniform characteristics as shown in Figure 22a and Figure23a, respectively. Figure 22b and Figure 23b present off-focus condition, which shows that both electrical intensity and phase are slightly distorted. The feed displacement from its original position affects the electrical intensity and the phase distribution linearity.   The performance for both types of lens are compared in Figure 24a,b. In terms of gain and beamwidth, the performances of both lenses are similar: there is not much difference, as shown in Figure 24a below. The performance for both types of lens are compared in Figure 24a,b. In terms of gain and beamwidth, the performances of both lenses are similar: there is not much difference, as shown in Figure 24a below.

Feed Position Analysis
In this section, the performance of the feed position for scanning angles of 20° until 40° is analyzed. The performance of the designed energy conservation lens and Abbe's sine lens is evaluated when the feed is located closer and farther from the lens. The original position is = as shown in Figure 25. θ is the scanning angle which is varied from 20° to 40°. The original position is the calculated position based on the shaped lens. Position 2 is where the feed is placed farther from the lens, while position 3 is where the feed is placed closer to the lens.

Feed Position Analysis
In this section, the performance of the feed position for scanning angles of 20 • until 40 • is analyzed. The performance of the designed energy conservation lens and Abbe's sine lens is evaluated when the feed is located closer and farther from the lens. The original position is F = f s as shown in Figure 25. θ is the scanning angle which is varied from 20 • to 40 • . The original position is the calculated position based on the shaped lens. Position 2 is where the feed is placed farther from the lens, while position 3 is where the feed is placed closer to the lens.

Feed Position Analysis
In this section, the performance of the feed position for scanning angles of 20° until 40° is analyzed. The performance of the designed energy conservation lens and Abbe's sine lens is evaluated when the feed is located closer and farther from the lens. The original position is = as shown in Figure 25. θ is the scanning angle which is varied from 20° to 40°. The original position is the calculated position based on the shaped lens. Position 2 is where the feed is placed farther from the lens, while position 3 is where the feed is placed closer to the lens.      Table 11 shows the overall performance for both types of lens at 20 • , 30 • and 40 • scanning angle. From all the results shown in Figure 26, Figure 27 and Table 11, it can be seen that all the positions along the lines satisfied the F = f s condition with optimum performance. It shows that the F = f s condition and both lens structures are suitable for multibeam applications. Hence, this structure can be applied for 5G mobile base station.   Table X shows the overall performance for both types of lens at 20°, 30° and 40° scanning angle. From all the results shown in Figure 26, Figure 27 and Table 11, it can be seen that all the positions along the lines satisfied the = condition with optimum performance. It shows that the = condition and both lens structures are suitable for multibeam applications. Hence, this structure can be applied for 5G mobile base station.

Conclusions
Two types of negative refractive indexes lens were designed in MATLAB software and simulated by using the electromagnetic simulator ANSYS HFSS. Both the energy conservation lens and Abbe's sine lens offer high gain and narrow beam width characteristics with 66% and 65% efficiency, respectively. A 27.55 dB maximum gain is achieved for the energy conservation lens and 27.48 dB for Abbe's sine lens. Moreover, both shaped lenses provide optimum results for beam scanning up until 40 • . High gain, narrow beam width and wide-angle beam scanning capability are the key elements for 5G application. Thus, it can be concluded that these two types of refractive index lens are good candidates for 5G mobile base station application.