A solution to the complement of the generalized Luneburg lens problem

Lenses are of interest for the design of directive antennas and multi-optics instruments in the microwave, terahertz and optical domains. Here, we introduce an optical problem defined as the complement of the well-known generalized Luneburg lens problem. The spherically symmetric inhomogeneous lenses obtained as solutions of this problem transform a given sphere in the homogeneous region outside of the lens into a virtual conjugate sphere, forming a virtual image from a real source. An analytical solution is proposed for the equivalent geodesic lens using the analogy between classical mechanics and geometrical optics. The refractive index profile of the corresponding inhomogeneous lens is then obtained using transformation optics. The focusing properties of this family of lenses are validated using ray-tracing models, further corroborated with full-wave simulations. The numerical results agree well with the predictions over the analyzed frequency bandwidth (10–30 GHz). This virtual focusing property may further benefit from recent developments in the fields of metamaterials and transformation optics. Spherically-symmetric lenses can create sharp virtual images, but a general expression for their refractive index profiles had not yet been developed. Here, this expression is provided via analogy between classical mechanics and geometrical optics, yielding solutions complementary to existing lenses obtained from the generalized Luneburg lens problem.

T he so-called generalized Luneburg lens problem was first formulated in its simplest form in the 1853 Cambridge and Dublin Mathematical Journal 1 , describing a "transparent medium [...] such that a ray of light within it is a given circle, the index of refraction being a function of the distance from a given point in the plane of the circle". The mathematical solution to that specific problem was published in the 1854 edition of the same journal 2 . An interesting feature of the described "transparent medium" is that "all the rays proceeding from any point in the medium will meet accurately in another point", a characteristic of absolute optical instruments, with no aberrations [3][4][5][6] . Although published anonymously, the problem and its solution are attributed to Maxwell 7 and the corresponding medium occupying the entire space was called Maxwell's fish-eye 3 . In the engineering community, the bounded version of this medium that is truncated at the normalized radius found more interest and is now generally called Maxwell's fish-eye lens. This latter design is a particular case of the more general problem formulated about a century later by Luneburg 8 where each point of a given sphere has a perfect image on another concentric sphere, the object and its image being in an homogeneous region outside or on the surface of the spherically symmetric inhomogeneous lens. For the particular case of Maxwell's fish-eye lens, both of the conjugate spheres are on the surface of the lens. In the case of the so-called Luneburg lens 8 , one sphere is on the surface of the lens and the conjugate sphere is at infinity, corresponding to the perfect collimation of a spherical wave emerging from a point source.
The definition of this problem was later extended by Eaton to include the possibility of having one of the spheres inside of the lens and to also have a point source and its image not aligned with the center of the lens 9 . A particular solution to this generalized problem is called the Eaton lens and works as an omnidirectional retro-reflector. Another particular case of interest is the Gutman lens 10 , which behaves as a Luneburg lens with its focal surface (or focal arc in the case of a planar implementation) inside the lens. Morgan extended that particular case by proposing a more general formulation of the Luneburg lens problem, which includes also the possibility of having one sphere inside the lens and the other outside the lens 11 . The formulation of the problem being reciprocal, the sphere inside the lens may correspond to either the source or the image. Besides the few specific cases mentioned above, the generalized Luneburg lens problem has no simple closed-form expression. Various methods had been proposed to solve the problem numerically, with residual aberrations depending on the achieved numerical accuracy [11][12][13] .
An equivalent to the planar Luneburg lens with a constant refractive index was proposed by Rinehart 14,15 and is known in its general form as the geodesic lens. The Rinehart-Luneburg lens has a simple closed-form expression when its rotationally symmetric profile is defined using the arc length measured from the symmetry axis as a function of the lens radius. This equivalence may be seen as an early implementation of conformal mapping, or transformation optics (TO) 16,17 , with Kunz 18 extending the particular transformation described by Rinehart to a more general equivalence between two rotationally symmetric surfaces. This property was recently exploited by Šarbort and Tyc 19 to provide a solution to the generalized Luneburg lens problem by deriving the analytical expression of its equivalent geodesic lens problem.
All the lenses discussed above are designed to produce a real image as a direct consequence of the formulation of the generalized Luneburg lens problem (see Fig. 1a). Miñano 4 noticed that an Eaton lens forms a virtual image diametrically opposite to the point source, corresponding to an absolute instrument with negative unit magnification. A similar property was obtained with the magnifying absolute instruments based on Eaton and Luneburg lenses described in 20 . Another class of spherically symmetric lenses providing virtual images is the so-called invisible lens, which produces an image collocated with the source, corresponding to a positive unit magnification 4,21,22 . Some less conventional geodesic and nonuniform refractive lenses were discussed in 23 to produce optical components equivalent to wellknown microwave components for use in microwave and optical circuitry. However, none of the spherically symmetric lenses reported above provide the more general transformation of a given sphere into a virtual conjugate sphere (see Fig. 1b). This problem is introduced here as the complement of the generalized Luneburg lens problem and is solved using the approach described in 19 . An analytical solution is also provided for its equivalent geodesic lens. This analytical formula can be used to evaluate the refractive index profile of any spherically and rotationally symmetric lens with the described virtual image property.
The Luneburg lens and its derivatives received a lot of attention since their inception because they provide unique scanning properties and require materials with relatively low refractive index, easy to engineer from existing materials. Early applications were mostly radar systems steering a highly directive beam over a wide angular range (up to full azimuthal coverage in some cases) 10,15,24,25 . With recent developments on metamaterials and TO [26][27][28][29][30] , there have been several advances reported in the microwave, terahertz, and optical domains [31][32][33][34][35][36] , with applications in terrestrial communications, satellite communications, and automotive radars to name a few [37][38][39][40][41] . Spherically and rotationally symmetric lenses producing sharp virtual images may be of interest in multi-optics systems, extending the range of existing solutions. They could be used in combination with reflector systems, improving for example the unifocal property of the hyperbolic sub-reflector in a classical dual-reflector Cassegrain geometry 42 . A unifocal lens with a virtual focal point was implemented in a dual-transmitarray configuration to reduce the overall height of the antenna system 43 . A similar unifocal duallens system based on homogeneous lenses is described in 44 . Fig. 1 Spherically symmetric inhomogeneous lens problems. A schematic representation of (a) the generalized Luneburg lens problem and (b) the proposed complementary problem is provided to define all design parameters. The circles delimit the boundaries between the inhomogeneous medium, n(r), and the homogeneous medium, n = 1. The inhomogeneous profile is such that any ray emerging from the source, represented by the point of coordinates (r s , ϕ s ) and hitting the lens, passes through a point image of coordinates (r i , ϕ i ) on the focal axis represented by dashed lines. In the case of the proposed complementary problem, it is as if the ray would emerge from the virtual image in a homogeneous medium, hence the partial representation of the ray in dotted line.
However, the extension of these specific solutions to the spherically and rotationally symmetric lenses described here is not straightforward. In its planar form, the general lens solution introduced in this paper could be combined with the doubly curved reflector geometry proposed in 45 .
This paper describes first the general problem, introduced as the complement of the generalized Luneburg lens problem. An analytical solution to the equivalent geodesic lens problem is derived, enabling the evaluation of different variations of the proposed lens supported by numerical results. An extension of the proposed general lens solution combined with simple mirrors is also discussed. Numerical results are provided to validate the proposed solution. Finally, some conclusions are drawn and perspectives are discussed.

Results and discussion
Problem formulation. The problem, as formulated by Luneburg himself, is to find a spherically symmetric partly inhomogeneous medium, corresponding to the desired lens, "such that two given spheres in the homogeneous part are perfect conjugate spheres" 8 . Including the extension formulated by Eaton, who suggested to remove "the restriction imposed by Luneburg that the emerging rays be parallel to the axis of symmetry of the system" 9 for the particular case of a conjugate sphere at infinity, one can schematically represent the problem to be solved with Fig. 1a. The spherical symmetry of this problem enables us to reduce its study to any plane passing through the center of the lens. Thus, it is convenient to use the polar coordinates (r, ϕ). The problem can also be normalized to the radius of the lens without loss of generality. A refractive index of 1 is assumed in the homogeneous region outside of the lens, here also without loss of generality as the problem may be scaled by any desired value. Note that the problem is described assuming a direction of propagation. However, the rays propagating in the forward direction could equally propagate in the backward direction, as a direct implication of the Stoke-Helmholtz reciprocity principle. Thus, no specific direction of propagation is highlighted in the ray-tracing representations provided here to emphasize that the source and the image may be swapped.
The problem is then to define the refractive index of the lens as a function of the radial coordinate, n(r), such that any ray emerging from a point source, at location (r s , ϕ s ), and hitting the lens, passes through a point image at location (r i , ϕ i ). As a consequence of the rotational symmetry, a distributed source on the sphere with radius r s produces a distributed image on the conjugate sphere with radius r i having a magnification of r i /r s . This property derives from the triangle similarity theorems. Note that with the generalization proposed by Eaton 9 , the image of a point source is in fact a circle with axis of symmetry the x-axis and radius r i sin ϕ i . This lens produces a conical wavefront in the particular case of a conjugate sphere at infinity. The image circle reduces to a point image when ϕ i À ϕ s ¼ 0 π ½ . The time to travel from the object to its image, or in other words the optical path length of the ray trajectory between the object and its image, is independent of the angular direction at the origin, α s . A convenient way to formulate and solve the problem is to consider the inverse problem, which consists in finding the ray trajectory assuming the refractive index, n(r), is known, and to make use of the well-known analogy between classical mechanics and geometrical optics. One may refer to 19 for a detailed derivation of the equations in the case of the generalized Luneburg lens problem and for the analytical expression of the corresponding equivalent geodesic lens.
Here a different problem is introduced. The problem remains to find a spherically symmetric medium such that two given spheres in the homogeneous part are perfect conjugate spheres. However, the problem is formulated such that the image is always virtual, as represented in Fig. 1b. The rays emerging from the lens are seen as if emerging from a point image (r i , ϕ i ) rather than converging to it. Equivalent parameter values are used in Fig. 1 to represent the two problems in order to emphasize their complementary nature. While the generalized Luneburg lens problem defines lenses focusing on one half of the focal axis (represented with dashed lines in Fig. 1), the reformulated problem defines lenses focusing on the complementary half of that focal axis. The complementary nature of the two problems is further highlighted through the specific examples discussed here. Because the proposed complementary problem maintains rotational symmetry, the magnification of the proposed lenses is also r i /r s .
Using the analogy between classical mechanics and geometrical optics, a ray trajectory is characterized by a constant value of the quantity L, labeled as an angular momentum in reference to its classical mechanics counterpart, and given by: where α is defined here as the angle measured towards the origin between the radial vector and the local tangent to the ray at any given point. Note that in some works, the supplementary angle is used instead, which does not change the value of L. This angle is represented in Fig. 1b for the particular cases of the point source and the point image with values α s and α i , respectively. This enables to evaluate L as a function of these angles and the radial position of the point source or point image as follows: There is another particular point of interest, also illustrated in Fig. 1, for which the ray trajectory has a minimum radius, r m . At this particular point, α = ± π/2. To simplify the analytical expressions, Luneburg 8 introduced the radial quantity ρ = rn(r), referred to as the turning parameter in recent literature 21 . This quantity corresponds to the radial coordinate transformation between a planar rotationally symmetric inhomogeneous lens and its geodesic equivalent 14,18 . Using differential geometry, the elementary angular displacement of a ray in a spherical medium with refractive index n(r) is given by: In ref. 19 , the parameter M is introduced to quantify the total change of polar angle as the ray propagates from the source to the image such that ϕ i − ϕ s = − Mπ, where M is a non-negative real number. Note that in the case of the nominal Luneburg lens problem 8 , where any point source and associated point image are aligned with the center of the lens, M = 1. In the proposed reformulated problem, the corresponding nominal case, where the point source, the point image, and the lens center are aligned, is obtained for M = 0. Using the same notation and counting the virtual portion of the ray trajectory negatively, the following equation is derived for the proposed complementary problem: The first and second terms on the right-hand side of the equation correspond to the real ray trajectory, these terms corresponding to the propagation outside the lens and inside the lens, respectively. The third term corresponds to the virtual propagation from the point image, and is the term counted negatively compared to the same equation in 19 . This leads to the following implicit equation, which defines r as a function of ρ: The quantity between brackets on the right-hand side is the scattering angle, χ, defined as the angular change to the ray trajectory in the homogeneous medium before and after the lens and counted positively when the ray is bent towards the x-axis. With this quantity, defined as a function of L, one can find the refractive index using the formula described in 8 : This formulation is based on a variation of the inverse Radon Transform introduced by Luneburg 8 , where the integration range is reduced to part of the lens region. Equation (6) provides the general solution to the proposed problem in the form (n(ρ), r(ρ)), defining the radial distribution of the refractive index of the lens, n(r), as a parametric function of ρ which may be solved numerically in the general case. In fact, Eq. (6) can be found directly from geometrical considerations giving the scattering angle, χ(L). However, the derivation leading to Eq. (5) is of interest as it helps defining an analytical formula for the proposed family of lenses using the equivalence with geodesic lenses, following an approach similar to the one described in 19 . A more general equivalence between two rotationally symmetric lenses on curved surfaces has been described 18,46 . For the purposes of this work, it is sufficient to consider the equivalence between a planar inhomogeneous lens and a homogeneous geodesic lens.
Before discussing the equivalent geodesic lens problem, one can already make some general remarks on the solutions to the complement of the generalized Luneburg lens problem formulated here. For the particular case χ = 0, Eq. (6) leads to n = 1, which is the trivial invisible lens solution. Another important remark, for the nominal case where M = 0, is that the sign of the scattering angle χ is independent of L and is positive when r i ≥ r s . This sign determines the sign of the integral in Eq. (6), meaning that the resulting refractive index is such that n(r) ≥ 1, ∀ r ∈ [0, 1], when r i ≥ r s . This property is true for any relative position of the source and its image in the case of the generalized Luneburg lens problem, as the scattering angle is always positive (see Fig. 1a). With the proposed lens, the scattering angle is negative for r i < r s and M = 0. This specific configuration results in a refractive index n(r) ≤ 1, ∀ r ∈ [0, 1], thus requiring to either adapt the refractive index of the homogeneous region outside the lens or to use adequate metamaterials. This opens opportunities for interesting developments.
Equivalent geodesic lens. The equivalence between a rotationally symmetry inhomogeneous lens and a geodesic lens is used to find an analytical solution to the proposed complementary problem 19 . The transformation, or conformal mapping, between the two lenses is defined by 14,18 : where s(ρ) is the arc length measured from the axis of symmetry of a rotationally symmetric geodesic lens defined in the cylindrical coordinate system (ρ, θ, z), as represented in Fig. 2. With this coordinate transformation, Eq. (3) becomes where L ¼ ρ sin α, with α being defined here as the angle between the ray trajectory and the meridian in the lens (see Fig. 2), and s 0 ðρÞ ¼ ds=dρ. This leads to the following implicit expression of the derivative of s(ρ), equivalent to Eq. (5): This equation differs from the one of the generalized Luneburg lens problem 19 in that the term with r i is negative. Thus, all terms that come from arcsinðL=r i Þ in 19 have to be taken with a negative sign. The term that comes from the generalization of the problem, Mπ/2, and which replaces Mπ=2 À arcsinðLÞ in 19 , needs to be reevaluated. The derivation of the associated term is detailed in the Methods and illustrates the use of Abel's integral equation to solve this problem. With these considerations, one can derive the following closed-form expression for the general equivalent geodesic lens solution to the proposed problem: Assuming r s = 1 and M = 0, this expression reduces to: The corresponding geodesic lens profile in cylindrical coordinates (ρ, θ, z) is obtained from the differential equation dz 2 = ds 2 − dρ 2 , while the inhomogeneous lens refractive index profile, n(r), is obtained from the coordinate transformation defined by Eq. (7). The following subsections provide numerical results, comparing the generalized Luneburg lens problem and the proposed complementary problem.
Discussion on complementarity. Numerical results are provided first to emphasize the complementarity of the two problems discussed above. The geodesic lens profiles are obtained using Eq. (10), while the inhomogeneous lens profiles are evaluated numerically from the equivalence defined by Eq. (7). Numerical results for the generalized Luneburg lens problem are obtained using corresponding formulas in 19 . The comparison of the geodesic lens profiles and inhomogeneous lens profiles are provided in Fig. 3 for the case r s = 1, corresponding to a source on the surface of the lens, and different values of the image radius, r i ≥ 1. A distinction is made using superscripts, with r r i corresponding to a real image lens design while r v i refers to a virtual image lens solution as proposed in this paper. The conventional Luneburg lens, corresponding to an image at infinity, is a common solution to both problems. In the case of the generalized Luneburg lens problem, the real image is located at + ∞, while in the case of the proposed complementary problem, the virtual image is at − ∞. This common solution corresponds to the yellow lines in Fig. 3. The generalized Luneburg lens problem covers the space between those curves and the ones of Maxwell's fish-eye lens, with an image on the surface of the lens and corresponding to the blue lines in Fig. 3. An example of intermediate lens design is provided for r r i ¼ 2:5 (orange lines). The profiles of the solutions with a real image inside of the lens 11 cover the space above Maxwell's fish-eye lens in Fig. 3.
The solutions to the proposed complementary problem occupy the space between the standard Luneburg lens and the trivial invisible lens solution, corresponding to a flat geodesic lens or an inhomogeneous lens with uniform refractive index equal to 1, both represented with green lines in Fig. 3. An example of intermediate lens design is provided for r v i ¼ 2:5 (purple lines). These results provide a good visualization of the complementary nature of the two problems. The proposed lens family leads to geodesic lens profiles lower than those of the generalized Luneburg lens problem, or similarly, lower refractive index profiles, which is expected to facilitate practical implementations and integration. As already observed in the subsection "Problem formulation", a virtual image inside of the lens would lead to solutions below the invisible lens case, r i = r s . This may be achieved in the case of inhomogeneous lenses with a refractive index below 1, using engineered materials. However, this is not achievable with geodesic lenses, as profiles with negative height values are in fact just symmetric designs of those with positive height values. This follows from the condition ds ≥ dρ which expresses the obvious fact that an elementary arc length, ds, cannot be shorter than the corresponding elementary radial displacement, dρ, in a geodesic lens. When combined with Eq. 7, this leads to the condition dn/dr ≤ 0 in the equivalent inhomogeneous lens. Thus, n(r) is a monotonically decreasing function that is always greater or equal to its value at the periphery of the lens and n(r) ≥ 1, ∀ r ∈ [0, 1], with the normalized notations used in this paper. The symmetric property of z(ρ) with respect to any plane perpendicular to the z-axis was exploited by Kunz 18 to fold the conventional geodesic lens, and was revisited recently to design low-profile geodesic lenses [39][40][41] . In particular, the numerical results reported in 40 confirm that the refractive index in the inhomogeneous lens equivalent to a modulated geodesic lens remains a monotonically decreasing function. Although the geodesic equivalent does not exist for inhomogeneous lenses with n(r) ≤ 1, ∀ r ∈ [0, 1], corresponding to the configuration with a virtual image closer to the lens center than the source, r i < r s , Eq. 10 still provides a numerical result that can be used to determine the refractive index, n(r), for that particular case, as demonstrated in the following subsection.
Ray-tracing models. Using the ray-tracing tool developed in 40 to analyze modulated geodesic lenses based on the non-Euclidian transformation of Eq. (7), some further numerical results are provided in this section exploring various configurations of the proposed family of lenses. The virtual image lens reported in Fig.  4a corresponds to the intermediate case also illustrated in Fig. 3 with r s = 1 and r i = 2.5. The rays have all the same length. A portion of an arc centered on the virtual image is illustrated with a red line in the case of the inhomogeneous lens, confirming the desired property. A ray-tracing representation of the equivalent geodesic lens is also provided in Fig. 4c. A second lens configuration is illustrated in Fig. 4b for the case of a source away from the lens, r s = 1.5, while keeping the virtual image at the same position, r i = 2.5. Moving the real source away from the lens and towards the virtual image reduces the maximum refractive index and corresponding height of the geodesic lens, as expected (see   . 4d). A ray-tracing representation is also provided for inhomogeneous lenses with r s > r i ≥ 1, thus leading to refractive index values below 1. Two configurations are illustrated in Fig. 5 for r i = 1 and r i = 1.5, with r s = 2.5 in both cases. The refractive index profiles are essentially symmetric to those obtained with r s < r i , having here the lower values at the center of the lens. Similarly to the case with refractive index above 1, moving the virtual image away from the lens and towards the real source reduces the refractive index variation across the lens, converging toward an invisible lens in the particular case where r i = r s .
All the results discussed above are obtained for M = 0. The general solution is also validated with M > 0, corresponding to a generalization of the problem equivalent to the one proposed by Eaton 9 . The ray-tracing representations reported in Fig. 6 are obtained with r i = 2 and M = 0.8. In Fig. 6a, the source is on the surface of the lens, while in Fig. 6b, the source is at infinity. A raytracing representation in the equivalent geodesic lens is provided in Fig. 6c and 6d for the configurations in Fig. 6a and 6b, respectively. The proposed complementary problem leads to the retro-reflective Eaton lens for r s = − ∞, r i = + ∞ and M = 1.
The configurations illustrated in Fig. 6 have a very high refractive index at the center of the lens. This may be significantly reduced using a half-lens design combined with a planar mirror. The concept of half-Luneburg lenses was introduced in the early 1950's to reduce the size and weight of Luneburg lenses 47 . With the addition of a planar mirror passing through the center of the  lens, the focusing properties are preserved and the source is transformed into a virtual source, still focusing the rays at infinity but in a symmetric direction with respect to the plane of the mirror. This property is also maintained in the case of planar 48 and geodesic 49 half-Luneburg lenses. This solution reduces the size of the lens by a factor of two, albeit with a reduced scanning range, as the omni-directional characteristic of the original Luneburg lens is lost. This solution may be of interest in some specific applications with limited field of view. This concept is also applicable to the virtual image lens introduced here. A raytracing representation of an inhomogeneous half-lens design is provided in Fig. 7. The lens is also designed with r i = 2, for comparison with the lens reported in Fig. 6a. The singularity is removed and the maximum for the refractive index is around 1.2.
Full-wave results. To validate further the theory, a full-wave model has been implemented in the frequency domain using the EM solver Ansys HFSS. The model is a geodesic lens in a parallel plate waveguide configuration with a homogeneous medium having a refractive index of 1. This modeling approach has already proved to be accurate for the characterization of geodesic lenses in the microwave domain 39,40,49 . A representation of the electric field distribution is provided in Fig. 8 for a virtual lens with parameters r s = 1, r i = 2.5 and M = 0, corresponding to the case illustrated in Fig. 4a, c. Numerical results are reported at 10, 20 and 30 GHz with a lens of 100 mm in diameter. The lens is excited with an isotropic source. These results confirm the frequency-independent response of the structure. The cylindrical  wavefront centered on the feed is clearly visible on the left side of the plots. The virtual image and expected cylindrical wavefront are highlighted in red. Besides some residual diffraction effects, the transformed wavefront is clearly visible on the right side of the plots and fully in line with predictions. The diffraction patterns are a consequence of the numerical implementation, requiring to discretize the lens profile, and of the transition between the surrounding planar parallel plate waveguide and the geodesic lens, leading to some residual aberrations. More generally, some limitations may result from the practical implementation of the lens. Nevertheless, these numerical results confirm the focusing properties of the virtual image lens solution detailed in this paper, even in the case of an electrically small lens as the analyzed model has a diameter of about three wavelengths only at the lowest reported frequency.
As mentioned in the introduction, the proposed family of lenses can also find applications in multi-optics systems. The analogy with dual-reflector configurations indicate that virtual image lenses may be used in a Cassegrain-like reflector geometry, which is known to provide a more compact layout 42 . A planar model of a lens-fed offset parabolic reflector geometry was implemented in Ansys HFSS. The numerical results for various positions of the source are reported in Fig. 9a-c. To further reduce the size of the proposed geometry, a half-lens geometry is implemented, corresponding to the case illustrated in Fig. 7. For these analyses, a more directive source is used to limit interference patterns coming from back radiation and to facilitate the visualization of the plane wave produced by the dual-optics system. Specifically, a rectangular waveguide with a broadwall dimension of 8.64 mm is used, corresponding to the standard size WR34 suitable for K a -band systems. Numerical results are reported at 20 GHz. The focal length of the offset parabolic reflector is set to 250 mm, while the lens diameter remains 100 mm. Due to the offset geometry, the distance that serves as reference for the magnification effect of this dual-optics system is the distance from the focal point to the reflector surface passing through the center of the lens. The reflector geometry was defined such that this distance is approximately 300 mm, thus magnifying the properties of the lens by a factor of 3 45 . The focal point of the parabolic reflector is set to coincide with the virtual image located at −45 ∘ in the coordinate frame with origin the center of the lens, corresponding to the configuration in Fig. 9b, where the virtual image is highlighted with a red dot. For sources placed at ± 15 ∘ around that nominal position, a tilt of the wavefront is observed while maintaining the plane wave characteristic. Due to the magnifying geometry, the expected tilt angle is equal to the source angular displacement divided by the magnification factor, corresponding to about ± 5 ∘ . The expected wavefront is  represented in Fig. 9 with a red solid line. The numerical results are well in line with predictions. These results provide a first demonstration of the use of the proposed virtual image lens in dual-optics systems.

Conclusions
We defined a family of spherically and rotationally symmetric lenses with a virtual image as the solution to the complement of the generalized Luneburg lens problem. An analytical formulation based on the equivalence between inhomogeneous lenses and geodesic lenses was proposed and numerical results were reported, emphasizing the complementary nature of these lenses with those obtained as a solution to the generalized Luneburg lens problem. A ray-tracing model was used to illustrate the focusing properties of various lens configurations. The solution was further validated using full-wave EM modeling, both as stand-alone and as feeding structure in a dual-optics system. The numerical results confirmed the expected virtual image focusing properties.
This complementary family of lenses is expected to extend the range of existing quasi-optical and optical systems, providing advanced focusing properties through the use of multi-optics configurations. Combined with the recent developments on metamaterials and transformation optics, this is expected to foster a number of new developments in the microwave, terahertz and optical domains.