{'fraction_non_alphanumeric': 0.024324670336704646, 'fraction_numerical': 0.013730636282166177, 'mean_word_length': 4.910896708286039, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
High-Directivity Emissions with Flexible Beam Numbers and Beam Directions Using Gradient-Refractive-Index Fractal Metamaterial

A three-dimensional (3D) highly-directive emission system is proposed to enable beam shaping and beam steering capabilities in wideband frequencies. It is composed of an omnidirectional source antenna and several 3D gradient-refractive-index (GRIN) lenses. To engineer a broadband impedance match, the design method for these 3D lenses is established under the scenario of free-space excitation by using a planar printed monopole. For realizations and demonstrations, a kind of GRIN metamaterial is proposed, which is constructed by non-uniform fractal geometries. Due to the non-resonant and deep-subwavelength features of the fractal elements, the resulting 3D GRIN metamaterial lenses have extra wide bandwidth (3 to 7.5 GHz), and are capable of manipulating electromagnetic wavefronts accurately, advancing the state of the art of available GRIN lenses. The proposal for the versatile highly-directive emissions has been confirmed by simulations and measurements, showing that not only the number of beams can be arbitrarily tailored but also the beam directions can be steerable. The proposal opens a new way to control broadband highly-directive emissions with pre-designed directions, promising great potentials in modern wireless communication systems.

M etamaterials, which contain a mass of artificial subwavelengh unit cells, have sparked exponentially increased interest among scientists and engieneers over the past decade due to the fantastic electromagnetic (EM) properties that are impossible to be achieved using natural materials. Since the initial motivation for scientific pursuit of negative refraction index in 2001 1 , physicists and engineers have pioneered new ways to guide and manipilate EM waves at will through appropriate choices and arrangements of meta-atoms (unit cells), giving rising to the directive 2 or enhanced 3 emissions, perfect imagings 4-6 , invisibe cloakings 7-9 , EM tunneling 10 , perfect absorbers [11][12][13][14] , and omnidirectional bending 15 , etc. In addition to the exhibited exciting phenomena in fundamental science, these conceptually new devices have broken the basic limitations for conventional devices, and thus provide an initial step towards unique functionalities. Nevertheless, the inherent drawbacks of large loss and limited bandwidth due to the resonant profiles of electric or magnetic inclusions as well as on occasion complicated fabrication process preclude applications of metamaterials in real world.
Despite these fruitful achievements, most previous reports were confined to the laboratory demonstration of the noted physics, for which the transverse-electric (TE) or transverse-magnetic (TM) polarized EM fields are examined extensively in a two-dimensional (2D) parallel-plate waveguide system [19][20][21][22][23][24][25][26][27][28][29][30][31][32][33] . Among them, the coaxial probe or point source was commonly employed as an ideal line source, which certainly conflicts with the freespace excitation in practice. Moreover, the homogenization of the GRIN materials can be further improved in terms of an effective medium. Most importantly, among the rarely reported three-dimensional (3D) lenses 34 the waveguide was exclusively adopted as the source and thus enabled only one desirable directive emission beam. This scheme to some extent deviates from the increasing demand on highly preferable directive emissions with flexible beam numbers and beam directions in a low-cost, low-weight and controllable manner. However, this 3D multi-beam highly-directive metamaterial lens antenna has not been reported yet so far. In this article, the effective excitation and broadband impedance match of a 3D GRIN metamaterial lens were comprehensively explored under a wave launcher of an actual printed monopole. The lens was designed in free space instead of 2D mimicking in the parallel-plate waveguide where the electric or magnetic fields perpendicular to the 2D plane are constant. A kind of GRIN metamaterial element constructed based on closed Sierpinski fractal ring was proposed, studied, and then employed as a fundamental block to build the 3D collimating lens with satisfied index gradients. Due to the omnidirectional feature of the monopole, a multi-beam antenna system can be readily conceived by employing more lenses. The proposed scheme has been verified by a far-field measurement of four GRIN lenses embedded with a printed monopole. To the best of our knowledge, this is the first realization of a 3D antenna system with flexible beam numbers and beam directions for broadband highly-directive emissions of more than an octave.

Results
Characterization of GRIN fractal metamaterials. The fractal scheme employed in this work is from the viewpoint of a much electrically smaller GRIN metamaterial element in virtue of the space-filling property, which is a new avenue advanced in material design 5 . Similar scheme has been employed to miniaturize the elements of composite right/left handed transmission line for a given resonant frequency 38,39 . The advantages of using rather a small element can be highlighted in twofold. First, the loading of one or two such elements may miniaturize considerably the planar circuits or antennas 5 , which is an intriguing characteristic that meets the increasing demand on high integration of modern communication systems. Second, the volumetric structures made of these periodically arranged elements are more feasible to be described by the constitutive parameters under the framework of an effective medium 40 and the resulting collective EM response is less influenced by the parasitic diffraction effects, enabling the elements in metamaterials to be analogous to the atoms or molecules in natural materials. Therefore, a desirable EM functionality can be engineered accurately by tailoring the element structures, which would benefit considerably the design accuracy and success rate. Fig. 1(a) shows the topology of the proposed GRIN metamaterial particle and its orientation with respect to the normally incident EM waves. The periodicity along x, y and z axis is p x , p y , and p z , respectively. The GRIN element is made of a fractal ring (FR) etched on one side of the F4B board (leaving the other side unoccupied) with dielectric constant e d 5 2.65, thickness h 5 0.3 mm, and loss tangent tans 5 0.001. Due to the four-fold rotational symmetry of the FR in xoy plane, the GRIN element presents a planar isotropy and thus supports dual-polarization operation. In this particular design, the GRIN element is illuminated by a normally incident plane EM wave with polarization in the x direction, and the Sierpinski ring of the second-iteration order is utilized by taking a tradeoff between miniaturization and fabrications. Indeed, any zigzag fractal structure with strong space-filling property may benefit the compact design.
In the simulation setup shown in Fig. 1(b), the top-and-bottom boundaries of the GRIN element are assigned as perfectly magnetic conductors (PMC), whereas the left-and-right boundaries are set as perfectly electric conductors (PEC). The impinging electric field drives the line inductor in the x direction, and the gap capacitor (inter-ring coupling) formed between neighbouring FR structures as well as the capacitor formed between close neighbouring sections of the FR itself, facilitating an electric resonance modeled by a LC resonator which restores and discharges alternately. This can be clearly shifts upwards as p x and p y increase simultaneously and appears to be saturated as p x and p y reach 9 mm. Although the decrease of p x and p y ensures increased gap capacitances, it either affords limited enhancement of capacitance or inherently brings about fabrication challenging and tolerances. Therefore, it is unadvisable to select this scheme to engineer a compact element.
Since the oscillating current forms and flows through the FR structure, the abnormal fractal zigzag boundary will succeed in extending the current path within a limited volume and thus efficiently lowers down the resonant frequency, enabling the element to be less than one-tenth of the working wavelength at f 0 . To afford an efficient control over f 0 in terms of variable inductors and capacitors in a simple manner, the Sierpinski ring is progressively scaled up and down by manipulating the scaling coefficient within a fixed lattice constant. In this regard, the effective parameters below f 0 change accordingly and thus enable us to achieve a controllable arithmetical progression of refractive index.
To illustrate the effects of fractal perturbation, the reflection and transmission coefficients of the proposed GRIN fractal particle and its conventional counterpart using the square ring (SR) 37 are characterized by the commercial software CST Microwave Studio. For fair comparisons, both rings are with the same scaling coefficient 0.77 and both elements have the same lattice constant p x 3 p y 3 p z 5 6 3 6 3 6 mm 3 . These S-parameters (not shown for brevity of contents) are then utilized to retrieve the effective constitutive parameters using the standard retrieval procedure 41 , as shown in Fig. 1(c). As is observed in both cases, the permittivity curve undergoes an electric resonance around f 0 , whereas the permeability curve suffers an antiresonance. Moreover, an obvious frequency reduction of f 0 near 2.9 GHz from 10.25 GHz to 7.34 GHz is inspected when the fractal geometry is introduced. Therefore, the GRIN fractal element obtains at least 28.3% size reduction for a given operation frequency. Moreover, the effective permittivity and permeability in both cases exhibit smooth slopes below f 0 and the permittivity of the fractal element suffers a sharper and stronger resonance at f 0 . As a consequence, the fractal element features larger index gradient with respect to its conventional counterpart.
Broadband Highly-Directive Emission System. Figs. 2(a)-(d) depict the topology and elaborated geometrical parameters of the proposed highly-directive emission system or in other words integrated lens antenna. As is shown, the antenna is composed of a broadband 3D GRIN collimating lens and an omnidirectional monopole to launch the lens. The distance from the monopole to the lens is D. Distinguished from the curved Fresnel lens which is often implemented by polystyrenes, the GRIN lens proposed here preserves flat surfaces and contains both core lens and coat lens (impedance matching layer) 37 to guarantee a wide operation band. Moreover, the GRIN lens is actualized by the proposed fractal metamaterials discussed previously, making current implementation distinguished from any earlier designs. Here, a planar microstrip trapezoid printed monopole deposited on a F4B substrate with e d 5 3.5, h 5 1.5 mm and tans 5 0.001 is explored as the source to achieve controllable uniform beams along several predicted directions with high integration, low weight, low cost, and broad band. To guarantee the bandwidth of the final antenna system, both the lens and monopole should operate in a wide band with fully overlapped frequencies. The wide bandwidth of the monopole is guaranteed by optimizing the geometrical parameters, whereas the design principle for a broadband lens is derived in the Methods section. Fig. 2(e) depicts the distributions of the theoretically calculated refractive index with discrete (solid line) and continuous (dashed line) profiles using the derived approach. As can been seen, the discrete index of the core lens varies gradually from 2.07 to 1.31, while that of the coat lens changes from 1.44 to 1.14. Notice that the discretized constitutive parameters have no penalty on lens function and performances 21 but enable the lens implementation more feasible and practical.
Although the inhomogeneous index gradient can be easily implemented by drilling isotropic holes in a dielectric slab 36 , the high fabrication cost, incompatibility with available printed circuit board (PCB) technology, and finite types of available aiguilles with discrete aperture dimensions limit the applications. Here, the proposed GRIN fractal particle is employed to realize the index gradient presented in Fig. 2(e) due to the unique EM response discussed in previous Section. Since the fractal element meets well the requirement for an effective medium and the material parameters exhibit a smooth dispersion in the off-resonant region, the undesirable distor- tion of optical quality and significant impedance change between two adjacent layers suffered in traditional approach will be mitigated or completely eliminated. Moreover, the EM wave outgoing from the monopole illuminates the lens with normal incidence along z axis and polarization along y axis. This EM condition coincides well with that presented in Fig. 1(a), and thus will render the desirable electric response discussed previously, and in turn the efficient excitation of the resulting lens.
Here, several guidelines are highlighted for an easy design. First, the larger the lattice constant along the propagation is, the smaller index gradient and lower index values become. Second, the thicker the dielectric board is, the larger permittivity while smaller permeability are, leading to smaller index and deteriorative impedance matching. Finally, the larger e d is, the larger n r is. The unit cells of the core and coat lens are built on two kinds of F4B boards with (e d , h, tans) 5 (2.65, 0.5 mm, 0.001) and (2.2, 0.5 mm, 0.001), corresponding to l 0 /12 3 l 0 /12 3 l 0 /20 and l 0 /12 3 l 0 /12 3 l 0 /10, respectively. Here, l 0 is the working wavelength in free space at 5 GHz.
For design and characterization, the GRIN element, the lens, and the integrated antenna system are numerically evaluated in CST. The retrieved effective material parameters for the core and coat lens at 5.5 GHz are portrayed in Figs. S2 and S3, respectively. In both cases, the dimensions of FR sweep from scale 5 0.4 to scale 5 0.6 in steps of 0.001. Therein, we observe that the permittivity and refractive index vary continuously, while the permeability almost maintains at unity. Therefore, the dimensions of FR reduce gradually from the center to the edge of the lens. Moreover, the refractive index of the core lens increases progressively from 1.37 to 2.09 while that of the coat lens increases from 1.17 to 1.57. These index gradients completely fulfill the theoretical requirements for the lens. Most importantly, the imaginary parts of the permittivity, permeability and refractive index in both cases are negligible across the entire sweeping range, indicating a low loss.
Given the parametric sweeping results, the lens can be exactly constructed based on the geometrical mapping process using a root-finding algorithm. Fig. 3 depicts the photograph of the fabricated 3D GRIN lens which is composed of ten core layers sandwiched by two coat layers on each side. In the radial direction, each concentric region shown in Fig. 2(d) is realized by two GRIN elements with identical geometrical parameters. Hence a total of 32 3 32 metallic patterns are positioned on either side of each F4B dielectric board. An additional 2.5-mm-thick foam plate with near-unity permeability and permittivity is placed between every two neighboring core layers to guarantee p z 5 3 mm, whereas an additional 5.5-mm-thick foam plate is inserted between every two coat layers to guarantee p z 5 6 mm. By tightly stacking these dielectric boards and foam plates alternatively and reinforcing them through adhesives, a 3D GRIN slab lens occupying an area of 170 3 170 3 54 mm 3 is implemented. The lens is packaged in hand-made foam to provide supporting frame and can be made more compact when smaller index gradients or less concentric layers are used. However, this would either degrade the accuracy of lens implementation or deteriorate the lens performances.
The single-beam highly-directive emission system can be directly generalized to a multi-beam one. Fig. 4 illustrates the schematic of various multi-beam emission systems and the fabricated prototype of a four-beam lens antenna for demonstration. As can be seen, two or four GRIN lenses are arranged around the source monopole with identical D. The highly-directive emissions can be controlled in following four manners. First, the number of beams can be engineered arbitrarily by appropriately distributing more GRIN lenses. Second, the beam directions can be arbitrarily manipulated by tailoring the included angles between neighboring lenses. Third, the beams can be engineered rotated by placing the lenses on a mechanically rotated platform isochronous with a time sequence and fixing the feeding antenna. Finally, the beam can be tailored nonuniform by varying D and index gradient of the lenses. Here, the direction of beam is defined corresponding to the orientation angle h with respect to z axis.  Numerical and experiemtal results. To validate the design method for the lens antenna, the gain is plotted as a function of D, see Fig. S4. The maximum gain of 14 dB occurs around D 5 90 mm which is used to design the index gradient. Consequently, the rationality of approximating the monopole as a point source and effective design of the lens are confirmed. Higher gain of more than 19.3 dB and larger aperture efficiency are achieved at 5 GHz when a horn antenna with length of D 5 90 mm and aperture of the lens is employed as the source at the cost of multiple beams. This level of gain approximates to the maximum gain of 20 dB according to G max~4 pL 2 lens =l 2 0 for which D should be sufficiently long. Therefore, the relative low gain of the lens antenna is attributable to the utilized omnidirectional monopole. In comparison with multi-beam emission system using several aperture antennas or arrays, our approach features moderate antenna gain, low loss, broad band, and requires only a single feed without using complicated feeding networks such as power dividers. Moreover, the distance D in our scheme can be significantly reduced relative to the bare horn for the same level of uniform phase distribution, enabling a compact multi-beam antenna.
To demonstrate the phase-collimating capability of the GRIN lens, the aperture fields in xoy plane for antennas with and without the GRIN lens are calculated and the electric field E y component is plotted, see Fig.  S5. Clearly, desirable phase-correction functionality can be identified from the uniform phase fronts across the entire aperture except for the edge of the lens. On the contrary, the phase front suffers significant fluctuations in the bare monopole case. Notice that the variation of phase front on the lens edge poses little influence on the radiation performances since most energy radiates through the central part.   5 illustrates the simulated 3D far-field radiation patterns and reflection coefficients S 11 of the monopole with and without GRIN lenses, aiming to provide an intuitive view on the controllable highly-directive emissions. Different scenarios of single beam, two beams, and four beams are extensively examined with different directions at 5.5 GHz. As is expected from Figs. 5(a)-(g), the omnidirectional patterns of the monopole in xoz plane (H plane) are transformed to several narrow-beam patterns with high directivity. With respect to the bare printed monopole, an average of 10 dB directivity enhancement is achieved along lens directions in all cases, which is a very attractive and challenging feature. Moreover, the monopolar pattern with null radiation at broadside still holds, indicating that these lenses only rearrange the energy distribution in H plane. Most importantly, it is demonstrated for the first time that the number of beams in free space can be tailored at will and the beams can be arbitrarily directed by simply steering the orientations of the lenses from an embedded actual printed monopole. The multi-beam directive emissions deviate from the multi-beam generation in 21,42,43 which are still constrained to surround a ideal probe with zero-refractive-index or GRIN metamaterials in a 2D waveguide. Our multi-beam control scheme provides more versatility, flexibility and stronger signal intensity, while also maintains easy implementation in practice. From Fig. 5(h), it is learned that the lenses introduced in all cases do not deteriorate the broadband impedance-matching performance of the printed monopole, for which the S 11 is better than 210 dB from 3.5 to 7.5 GHz. Fig. 6 gives the simulated far-field patterns in H plane and the realized gain at a number of frequencies from 3 to 8 GHz. Here, the antenna gain in all cases is calculated by taking losses into account. From Figs. 6(a)-(f), we learn that the sidelobes increase appreciably at the lower and upper edges of the passband, which is attributable to the dispersive GRIN material parameters. In this scenario, the constitutive parameters of the lens deviate slightly from the theoretically calculated ones when the frequencies deviate from the centre frequency 5.5 GHz at which the lens is exactly designed. Nevertheless, the slightly deteriorative sidelobe in the off-center band does not preclude the system from real-world applications since the directivity of main beams does not decay until 6.5 GHz. From Fig. 6 (g), an obvious gain enhancement is observed across the whole band in all cases and the minimum value is found to be 6.5 dB at 3 GHz. The gain decay at the lower frequency end is due to the size of the aperture in terms of the wavelength decreases significantly at that end, whereas the dispersive GRIN material parameters account for the raised return loss and material losses and thus give rise to the gain decay at the upper frequency end. Despite this, the radiation efficiency of the lens antenna g calculated as the proportion of the gain and directivity is larger than 90%, further demonstrating that the losses induced by the lens are extremely small. Moreover, the gain is very appealing and comparable with respect to aforementioned bare horn antenna, see the dashed blue curve. This is especially true for the gain at high frequencies.
To provide physical insight regarding the highly-directive emissions, the near-field mapping is performed in 3D free-space mea-  www.nature.com/scientificreports SCIENTIFIC REPORTS | 4 : 5744 | DOI: 10.1038/srep05744 surement apparatus. In the experimental setup, see Fig. S6, the broadband monopole placed on an aluminum plate is 90 mm away from the GRIN lens and is embedded in foam to provide mechanical support. It functions as a transmitter and moves together with an electronic step motor automatically along two dimensions by scanning an area of 200 3 200 mm 2 . At the front side of the lens, a 1.19 mm-diameter semi-rigid coaxial probe is fixed stationary to a ferric bracket to detect the electric fields. Both transmitting and receiving antennas are connected to the two ports of a phase-sensitive network analyzer (Agilent PNA-LN5230C) via two coaxial cables and are placed at the waist of the GRIN lens to avoid truncation effects of finite lens dimension. Fig. 7 plots the reflection coefficients and measured free-space electric-field distributions. The successful conversion of quasi-spherical waves to plane waves is clearly observed from 4.5 to 8 GHz in steps of 0.5 GHz by comparing the wavefronts outgoing from the lens antenna and bare monopole, respectively. Here, the results at the lower edge of passband are not given due to the large wavelength at these frequencies. In this regard, the long wavefront exceeds the scanning capability of available near-field measurement apparatus, making it impossible to clearly identify a complete wavefront. Nevertheless, the highly-directive emissions at these frequencies can be evidenced from the measured far-field patterns. Since the GRIN elements operate far from their resonant frequencies, the loss induced by the sample is extremely low and is hardly indicated from the measurements. Therefore, almost identical intensity can be inspected from the color maps of outgoing waves in all cases. The nonideal monopole with finite ground dimensions gives rise to the uneven and unsmooth wavefronts. Moreover, the measured reflection coefficients are in reasonable agreement with simulations, revealing that the return loss of the bare monopole and integrated lens antennas under different scenarios is better than 210 dB from 3.2 to 7.75 GHz, which is more than an octave. Therefore, the accuracy of employing GRIN fractal elements in constructing bulk GRIN lenses is verified.
To experimentally validate the appealing far-field characteristics of the 3D antenna system, Fig. 8 compares the simulated and measured H-plane radiation patterns of both single-beam and four-beam antennas from 3.5 to 7.5 GHz in steps of 2 GHz. The radiation patterns in other frequencies can be observed from Fig. S7. All patterns are normalized to the co-polarization of the designed monopole. As is inspected, the results between simulations and measurements are in reasonable consistency. In both single-beam and four-beam cases, the gain enhancement is observed more than 10 dB along certain orientations. The measured cross-polarization radiations are at least 15 dB lower than the co-polarization ones at most frequencies and become worse at both edges of the passband. The nonuniform intensity (amplitude fluctuations) of the four main beams is attributed to the slightly broken uniformity of the monopole pattern which is induced by the insufficiently large ground plane. The appreciably higher crosspolarization radiations in experiments, and slight discrepancies of copolarization patterns reflection coefficients (especially the shift of reflection zeros) between simulations and measurements are partially attributable to the tolerances that are inherent in fabrications, and partially to effects of supporting foam and adhesives that are not considered in simulations. Moreover, the misalignment of transmitting and receiving antennas and the nonideal background environment also account for the observed deviations. Nevertheless, these tolerances pose little influence on lens function.

Discussion
To conclude, we have succeeded in demonstrating a highly-directive emission system which possesses the merits of broad band, low loss, high gain, and flexible beam numbers and beam directions. Numerical and experimental results are in good consistency, showing that the emission system exhibits a near 10 dB gain enhancement, and moderate cross-polarized radiations which are 15 dB lower than the co-polarized radiations over an octave in super- extended C band. Therefore, the effective design for a broadband 3D slab lenses under the free-space excitation of a planar printed monopole and the advantages of using fractal strategy in 3D GRIN lens design without increasing the manufacturing complexity have been confirmed. The proposed scheme introduces an alternative avenue to improve the gain of microstrip antennas from a material perspective. Other omnidirectional source with uniform patterns can be employed for further improvements.

Methods
Suppose that the aperture of the monopole is negligible with respect to that of the 3D lens, the monopole can be approximated as a point source in design and the diffraction in terms of truncation effects can be eliminated. Despite this, the design process is more complicated than 2D design by using ideal point source, which will be discussed later. Denote the thickness and side length of the lens as t 5 t core 1 2t coat and L lens , respectively and assume that the wave inside the lens travels parallel to the optical axis without a lean. Here, t core and t coat are thicknesses of the core and coat lens, respectively. To guarantee the phase-collimating performance of the lens, the optical path from the source to any point in the upper surface of the GRIN lens should be the same, yielding the same phase delay.
where n 0 is the refractive index at the center of the lens, and n r is the radially distributed index as a function of radius r. Then, n r can be immediately obtained from Eq. (1).
Although the required index of the lens is obtained, it is impossible to achieve the rigorous criterion e r 5 m r 5 n r in broadband and there is a large index fluctuation at the interface of lens and air, resulting in a severe impedance mismatch. To relax the fluctuation and thus enable a smooth transition of index and desirable impedance matching, the aforementioned one-layer lens is substituted by a core lens sandwiched by two identical coat lenses 37 each functioning as a one-stage l/4 impedance transformer, see Fig. 2(c). To guarantee simultaneously the beam-collimating function, the optical path for the three-layer lens should meet the following criterion.
n core t core z2n coat t coat~nr t ð3Þ Since the effective permeability in this particular design maintains as unity across the entire band, we have n core~ffi ffiffiffiffiffiffiffi e core p and n coat~ffi ffiffiffiffiffiffiffi e coat p . Then e coat~ffi ffiffiffiffiffiffiffiffiffiffiffi e core e 0 p and n coat~ffi ffiffiffiffiffiffiffiffiffiffiffiffi ffi n core n 0 p can be immediately obtained by taking the one-stage l/4 transformer condition Z coat~ffi ffiffiffiffiffiffiffiffiffiffiffiffiffi Z core Z 0 p . By inserting these equations into Eq. (3), the radially distributed index of the coat and core lens can be achieved after some simplifications.
Here, the thickness of the coat lens can be roughly determined corresponding to l 0 /4, whereas that of the core lens should be designed as l 0 /2 to guarantee desirable phasecollimating performances. In actual realizations, the 3D collimating lens is discretized into eight concentric circular regions each having a specific refractive index. n 0 is cautiously chosen as n 0 5 1.8 to ease the fabrications in achieving large index gradients. Since the monopole is not an ideal point source, it is hard to strictly align the radiation center of the monopole with that of the lens, making the simple ray tracing do not always hold. To make the approximation more accurately, the aperture of the lens and D should be sufficiently large and should fulfill the condition ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi D 2 zL 2 {Dvs. The parameter s defines the accuracy criterion and is chosen as s 5 l 0 /6 in this particular design. Moreover, the relative position (xoy plane) of the monopole relative to the lens is optimized to have a maximum peak gain for a given D. Here, D plays a crucial role in determining the antenna performances such as gain when the position of the monopole is determined. To afford efficient excitation and maximum aperture efficiency, it is designed and optimized as D 5 90 mm by comparing the performances of antennas with different D and index gradients. In this connection, the radiation center of the outgoing wave can be tailored at the geometrical center of the lens, which guarantees a desirable gain and directivity.