Maxwell Fisheye Lens Based Retrodirective Array

A Maxwell fisheye lens using parallel plate index grading is presented in this study to develop a passive retrodirective antenna array. As a proof-of-concept a design frequency of 10 GHz was selected for fabrication and experiment. The design principals of the lens are discussed, which enables 85% energy flow at the drain probe (also referred to as image point) of the lens. It is shown that the image in the Maxwell fisheye lens has a point symmetry with a reverse phase, which makes it possible to realize passive retrodirective action using the lens. This arrangement is significantly more practical than previous passive retrodirective topologies due to the un-constrained number of connections to radiating elements that it can support without the need for multi-layer technology. In the realization described here, a cross-polarized microstrip patch antenna array is connected to the source and drain probes of the lens structure in order to form the retrodirective array. The strategy for selecting the optimal transmission line lengths required to connect the antennas to the lens for maximum re-radiation power is described and implemented. Experimental results for a prototype high efficiency passive retrodirective array based on the theoretical design considerations presented in this paper are reported.


Maxwell Fisheye Lens Design Principles
Let x, y, z represent a point in a Cartesian coordinate system within which a ray trajectory can be written as r(t) = x(t), y(t), z(t) where t is the parametric representation throughout the trajectory. A set of rays in a spherically symmetric optical system can be projected into a plane without losing generality. In our case, the lens is cylindrical and lies in a disk in the xy-plane with the center of the lens located at z = 0. Here ray trajectories r(t) can be calculated using conformal mapping of the points on a sphere onto the aforementioned xy-plane 2 . For this case, the path length differential on the sphere is mapped as an optical length differential in the xy-plane. In this way, the geodesic points on the great circle of the sphere can be mapped directly to a fisheye refractive index profile. For a spherical coordinate system with θ = 0° defining the xy-plane the required mapping takes the form, sin cos cos2 and 1 sin cos sin2 (1) 2 2 This mapping is illustrated in Fig. 1(a). This is known as stereographic projection 2 in which southern hemisphere of a sphere is mapped to the interior, while the northern hemisphere is mapped to the exterior of a circle in xy-plane having radius equal to the radius of the sphere. The line element segments dx and dy can be expressed in terms of spherical coordinates, and for a cylindrical lens θ = 0° as open www.nature.com/scientificreports www.nature.com/scientificreports/ such that at the equator, n is equal to the index n 0 of the reference sphere. Typical lens profiles which emulate the fisheye lens 3 in a 2D plane can be described by: where, 2n 0 is the refractive index at the center of the lens, and R is the characteristic radius of the reference sphere. After translating 3D fisheye characteristics into 2D, the next step is its practical realization. Generally, one can implement a variety of technologies in order to create a refractive index profile that performs fisheye focusing e.g. thin plates 4 , holy parallel plates 5 . In this paper we take a different approach in which fisheye principle is realized using the mode theory of a parallel-plate waveguide 6 . Here between ideally conducting parallel plates with spacing d we introduce, a medium characterized by dielectric ε r (ω) given in Fig. 1(b). The general form for the propagation constant of the fields propagating parallel to the plates is 7 , when ∈  m {1, 2, } and ω ε = k c / r . The fields will propagate as long as k m is real, however when π ε > d k / r , k m becomes imaginary the propagating wave decays exponentially along the direction of propagation. The parallel plate waveguide then acts as a high pass filter with cut-off frequency Here m is considered the mode number. For the propagating wave with ω > ω c , the phase velocity will be v = ω/k m and the energy will transport at a group velocity v g = dω/dk m . We can control this group velocity by using a graded index profile following the fisheye Eq. (3). In this paper this is achieved by partially filling the parallel plate waveguide with dielectric material whose permittivity ɛ r is fixed but whose thickness varies across the lens we can create an effective dielectric ε ′ r such that.
where r is the radius of the parallel plates and t is the dielectric material thickness. Note that the notations ɛ r in this paper represents general relative permittivity of a material, whereas ɛ r ' represents the achieved effective permittivity inside partial dielectric filled metallic parallel plates. The method of controlling the permittivity in such a case is given in 8 . Considering (3), (5) and (6) we formulated the dielectric profile between parallel plates (blue disk in Fig. 1(c)). The next step is to excite EM energy in the parallel plate region, in our case by using coax probes. An important physical limitation needs to be considered when realizing the assembly in Fig. 1(c). The approach following (3) and (6) will eventually reduce the lens profile to the near-zero-index area, which to realize a perfect fisheye operation requires the lens to be infinitely large. We side-step this by applying a circular mirror (a cylindrical reflecting metal sheet) to confine the wave propagation within a finite zone as theoretically implied in 9 . To understand how the mirror functions, consider a conformal fisheye index profile as in Fig. 2, and its multiple circular ray trajectories. For a source placed within a confined space surrounded by a mirror (white line), the incident rays (solid lines) are reflected into ray trajectories (dotted line) that focus on the image point within the now confined space.

Fisheye Lens Development
The fisheye lens schematic presented in Fig. 3(a) comprises of two perfectly electrically conducting, PEC, plates spaced with d = 5 mm. The substrate profile used to fill in the parallel plate waveguide is Rexolite 10 , ɛ r = 2.53, dispersion factor = 0.00066, and low coefficient of linear thermal expansion = 3.8 × 10 −5 1/°F in inches. A circular band of copper of diameter = 100 mm is used to define a closed boundary for the lens, see Fig. 2. A coaxial probe www.nature.com/scientificreports www.nature.com/scientificreports/ with h = 4.3 mm (selected for 50Ω match) is used to excite the lens while the same probe length is used to form the drain point. A Triumph Dulex milling machine is used to machine the lens profile from a cylindrical Rexolite stock, and then a Tech-Gen precision finisher is deployed in order to polish the surface, making sure to have a surface roughness of <λ/8 across the lens profile. Full-wave EM simulations using CST Microwave Studio 11 show the maximum electric fields within this confined space. Energy focusing at the image point is illustrated via the field strength plot in Fig. 3(b) confirming the operation of the device. Figure 4(a) shows the focusing capabilities of the device, here ~85% simulated E-field can be extracted from the drain point. Since the lens structure is confined, we believe the remaining field loss is associated with imperfect imaging in the active drain located at the image point of the lens, as argued in 12,13 . A wave leaving the source in the fisheye lens follows a circular path and focuses at the image point. The real component of the electric field in time lapse images are shown in Fig. 4(b), here the difference in the propagating EM waves along the edges and centre of the lens can be clearly seen. In simulation we noted that if we locate the mirror very close to the probe location, energy loss to impedance mismatch occurs, which is undesirable for the device operation. The waves propagating between the outer edge of the lens (mirror) and the start of the Rexolite dielectric are found to have the same v g which aids in image formation at the drain probe through phase alignment. Due to multiple transitions, from mirrored parallel plate, to dielectric filled parallel plate and coax-type probe (drain) it is not straight forward to identify which modes are propagating in these portions of the fisheye lens structure. The wave equation solution for two parallel plates are characterized by m (Eq. (5)) and TE m and TM m modes are defined on discrete frequencies with specific cut-off frequencies. It is important to note that a solution of wave equation with no magnetic-field in the direction of propagation (TM m ) has a fundamental mode (TM 0 ) that has no cut-off frequency. Considering an infinite parallel plate waveguide where wave is propagating along z-direction, at the TM 0 mode, the propagation constant is equal to k. This means that the TM 0 mode has neither electric nor magnetic field along the direction of propagation, making it a transverse electro-magnetic (quasi-TEM) mode, also supported by the lens structure.  www.nature.com/scientificreports www.nature.com/scientificreports/ Drain probe location in our work verifies the location identified by the time reversed scattering and transformation optics techniques [13][14][15][16] , however it is different from the optical fisheye lens in 17 . Note that the operation frequency has a direct link to the optimal drain location in a mirrored fisheye lens. Since in this work we focus only on a narrow band around 10 GHz frequency operation, further details of this relation are not emphasized here and can be found in 18 . Other ways of close-to-perfect drain from fisheye lens is with the waveguide feeds [19][20][21] .
To confirm the phase conjugation, we simulated E z and H x 2.5 mm inside the lens structure using a 2D mesh with pixel size = 1.15 × 1.15 mm 2 . The results are shown in Fig. 5(a-d). We used these results to compute the wave impedance phase as shown in Fig. 5(e). 180-degree transformation of the phase is observed when comparing the semi-circular side contain the source to the semi-circular side containing the image, confirming P* at the drain probe. The wave impedance map Fig. 5(e) allows us to choose the impedance location for the coaxial probe insertion position and its depth such that 50 Ω circuit impedance match to the external environment is ensured.

Retrodirective Array Operation
The Van Atta array requires multiple crossovers 22,23 hence multilayer implementation which can be lossy and are difficult to rescale for larger arrays. In the fisheye lens, when multiple sources and drain probes pairs are placed along the outer boundary of the profiled Rexolite dielectric (Fig. 3), phase conjugation analogous to Van Atta phase aligned transmission lines (e.g. 24,25 ) is possible without the requirement of crossovers. Next, we show how to create a retrodirective reflector using the fisheye lens.
We first designed a λ/2 spaced 6-element microstrip patch antenna array operating at 10 GHz. The antenna unit cell used in the array is similar to a dual-feed microstrip patch antenna configuration 26 and is optimized for 10 GHz operation. Also, the antenna used in this work have an optimized location of coax feed such that the distance of co-polarization (co-pol) and cross-polarization (x-pol) feed is same from the closest edge of the antenna. This makes the antenna structure square. The patch antenna unit cell is designed and fabricated on a 0.813 mm thick Rogers 4003 C grounded substrate (relative permittivity ɛ r = 3.55 and electric loss tangent tanδ = 0.0027). Same unit cell is replicated to form 6-element array. Based on the incoming wave's direction of arrival (DoA), we deduced the re-radiated and scattered patterns by qualitative and physical considerations when standard phased aligned transmission lines are replaced by the fisheye lens, e.g. Figure 6 shows the patterns when a plane incident wave at DoA = 30° hits the patch antenna array. As required for retrodirective array functionality the resulting magnitude patterns evaluated in the plane normal to the E-field component of the incident wave shows that there is a re-radiated portion of field toward the DoA and also a scattered field which has maximum field strength specular to the direction of the incoming wave. We find that the overall re-radiated energy is affected by the length of the transmission lines connecting the fisheye lens to the array antenna. Consider all the lens to array interconnect transmission line lengths to be equal, and when ∈  p {0, 1, 2, } and < < q 0 1 . The length of a transmission line in (7) is deconstructed in two parts. First part is an integer multiple of wavelength, i.e. 1λ, 2λ, 3λ … and second is an additional section of a transmission that is subwavelength long represented by q. As observed in 27 length q has an impact on the re-radiation and this impact is most pronounced along the DoA. Another observation is that the phase of the re-radiated field changes with a change in q, while it has negligible impact on the phase of the field along the direction of reflected waves (in this case about 330°). We also observed that for the microstrip patch antenna array, maximum re-radiation occurs when q = 0.6λ which is close to the conclusions in 27 for a classical Van Atta configuration. Optimum transmission line length governed by q has an impact of approximately 1.5 dB in the re-radiated field strength along the DoA. With an optimal q = 0.6λ, retrodirective array performance for several representative DoA's is shown in Fig. 7. As the DoA is skewed away from the broadside direction of the array, the re-radiated field strength weakens such that the difference between re-radiated field at DoA = 60° is ~13 dB below compared to the re-radiated field at DoA = 0°. To isolate the retro and scattered patterns we used orthogonal polarization diversity such that the www.nature.com/scientificreports www.nature.com/scientificreports/ co-polarization component of the antenna is connected to the source probe while the cross-polarization component is attached to drain probe. It is important to mention that mutual coupling will have an impact on the results, and in the results below this is taken into consideration. The impact is reduced by designing the antennas with low return loss (<−10 dB), high port isolation between co-and cross-polarization ports of the same antenna (i.e. < −20 dB) and minimum achievable coupling between neighbouring array elements. Fabricated microstrip patch antenna unit cell, co-and cross-polarization feed, and the 6-element array is presented in Fig. 8(a). Note that we used phase aligned coax cables with SMA type RF connectors as transmission line of length l. Simulated   www.nature.com/scientificreports www.nature.com/scientificreports/ and measured results of the array are shown in Fig. 8(b). Measured |S 11 | <−10 dB bandwidth of the antenna unit cell is from 9.88 GHz to 10.11 GHz. Although the proposed fisheye lens structure is frequency dependent 18 but since our application demanded the bandwidth limitation to <300 MHz, we consider the same limit for the fisheye lens design. Results for the wave with a DoA = 15° are presented in Fig. 8(c,d) when E θ and E φ represent the cross-and co-polarization components of the fields. Undesirable scattered fields are shown to be nullified by the proposed approach. Results evaluated by experiments shown in Fig. 9 (E θ component in far-field) verifies the theoretical prediction.
The antenna unit cell number in passive retrodirective array in our work is 6, while this number is different in other retrodirective array examples in literature. Uniform linear array (ULA) configuration, for instance, are shown to have 2 28 , 4 and 29,30 16 31 antenna units. In addition to azimuth direction, elevation direction are also captured using uniform rectangular array (URA) configurations, shown to have 25 32 , 32 33,34 , even 40 35 antenna unit cells. This makes it very difficult to directly compare our approach with different retrodirective mechanisms. However, the measured E θ -field results in Fig. 9 show a normalized comparison with the simulated scenario, proving the proposed system to be a practically viable solution. conclusion A design approach for a cylindrically bounded Maxwell fisheye lens connected to a microstrip antenna array for retrodirective operation is presented. The operation of the lens connected to a linear array is shown to produce retrodirective action, wherein the re-radiated retrodirected beam has high isolation from un-wanted specular reflections by using polarisation diversity. Engineering details including the optimum fisheye lens design and transmission line lengths connecting the lens to the array are discussed. The approach suggests significant enhancement in the degrees-of-freedom for the implementation of high efficiency passive retrodirective array.

Data availability
The dataset based on measurements is available at go.qub.ac.uk/fisheye-lens.