## Abstract

A method to design isotropic inhomogeneous refractive index distribution is presented, in which the scalar wave field solutions propagate exactly on an eikonal function (*i.e.*, remaining constant on the Geometrical Optics wavefronts). This method is applied to the design of “dipole lenses”, which perfectly focus a scalar wave field emitted from a point source onto a point absorber, in both two and three dimensions. Also, the Maxwell fish-eye lens in two and three dimensions is analysed.

© 2010 OSA

## 1. Introduction

It is common to deduce Geometrical Optics as a limit case of Scalar Wave Optics, for instance, when the wavenumber *k* = *ω*/*c* is very large. In that limit, the propagation of the scalar field can be calculated with good approximation using rays. There are, however, some trivial cases in which the ray trajectories guide the scalar field in an exact manner (*i.e.* with no restriction to large *k*) so that the field is constant on the Geometrical Optics wavefronts. For instance, a (monopole) point source emitting from the center of a spherical symmetric refractive index distribution *n*(**r**) will generate a field which will depend on the radial coordinate only. In this paper we are going to discuss about isotropic inhomogeneous refractive index distributions that propagate nontrivial scalar fields exactly on eikonals (*i.e.*, remaining constant on the Geometrical Optics wavefronts), and we will find media that produce perfect focusing of rays and waves, perfect in the sense explained next.

In Geometrical Optics, an optical system is said to produce a perfect focus (or sharp image) of an object point **P** onto an image point **Q** when any ray trajectory emitted from **P** through the optical system will pass through **Q** in an exact way. Such points **P** and **Q** are said to be perfect conjugates [1]. A device is called an Absolute Instrument in Geometrical Optics if it produces perfect focusing of rays not just from a single object point, but of all points in a three-dimensional domain (*i.e.*, one with non-null volume) [1].

The ellipsoidal mirror is a well known example of perfect focusing of rays, but only for points **P** and **Q** coincident with the foci of the ellipsoid (*i.e.*, it is not an Absolute Instrument). Non-trivial examples of Absolute Instruments are based on non-homogeneous media [2], the most famous being the Maxwell fish-eye lens [1]. Unlike conventional imaging optics systems, in the Maxwell fish-eye lens the refractive index in the volume containing the object and image points is inhomogeneous, (*i.e.*, spatially varying) . There are, however, examples of Absolute Instruments in which the refractive index distribution in that volume is homogeneous [3].

An analogous concept to the performance of an Absolute Instrument in Scalar Wave and Electromagnetic Optics has been introduced by Pendry [4] in the field of metamaterials, under the rubric of Perfect Imaging, meaning image formation with unlimited resolution. It disclosed a specific device made of a negative refractive index slab with *n* = *ε* = *μ* = –1, the capacity of which for perfect imaging has been explained by the amplification of evanescent waves in the negative index material.

Recently, Leonhardt [5] has proven that Perfect Imaging for 2D Helmholtz scalar fields is also achieved by the cylindrical Maxwell fish-eye lens, and the mirrored version referred to as the Maxwell fish-eye mirror. In three dimensions, Leonhardt and Philbin have proven very recently [6] that the spherical Maxwell fish eye lens is perfect for focusing electromagnetic waves (not Helmholtz scalar wave fields) if the medium is impedance-matched (*n* = *ε* = *μ*). These results are especially relevant because it uses isotropic positive refractive index.

The fact that the cylindrical Maxwell fish-eye lens is perfect for focusing both rays and scalar field waves in two dimensions is not a general result. For instance, even though the elliptic mirror is ideal for rays when **P** and **Q** coincide with the foci, it is well known that it does not focus waves perfectly (neither in two nor in three dimensions) due to the coma of the mirror.

In this paper we are going to discuss systems that perfectly focus Helmholtz scalar field waves emitted from a point source at **P** onto a point drain at **Q**. This means that the local behavior of the field around **Q** will coincide asymptotically with a spherical converging wave in three dimensions. We will also consider the two dimensional case, valid for cylindrical waves in cylindrical symmetric optics.

Section 2 deals with the conditions that must be fulfilled for the Helmholtz scalar wave field to propagate exactly on an eikonal function, leading to a novel constructive design method. In section 3, this method is applied to cylindrical lenses, in particular to reproduce with a different approach the performance of the cylindrical Maxwell fish-eye in reference [5]. Section 4 will consider the three dimensional case, concluding with the design of a “3D dipole lens”, which perfectly focuses a point source **P** upon a point drain **Q** in three dimensions. Finally, Section 5 discusses resolution limits.

As a clarification, throughout this paper the term ‘Maxwell fish-eye’ will refer to the spatially unbounded refractive index distribution given by Eqs. (12) and (28), in the cylindrical and spherical symmetric cases, respectively, and not to the mirrored version discussed in references [3] and [5].

## 2. Helmholtz fields that propagate exactly on eikonals

#### 2.1 Statement of the problem

Consider a scalar field *U*(**r**) ∈ C, **r** ∈ *D* ⊂ R^{3} in a medium with refractive index distribution *n*(**r**):

*k*=

*ω*/

*c*. As with [5] and [7], Eq. (1) is herein referred to as Helmholtz equation for inhomogeneous media (or just the Helmholtz equation for short) and its solution

*U*(

**r**) is referred to a Helmholtz scalar wave field. Note that the name Helmholtz equation is sometimes reserved to the case

*n*= constant (i.e., the case of homogeneous media), while Eq. (1) is formally equivalent to the time-independent Schrödinger equation. This equation is relevant in other areas of physics, such as acoustics or optics. In optics, this equation in 2D is exact for describing TE polarized light in cylindrical media (in which electric field vector

**E**points orthogonal to the cross section of the cylinder). It is not exact but approximate, however, for describing electromagnetic fields in 3D.

Consider a function *S*(**r**) ∈ R that is a particular solution of the eikonal equation in the domain *D*:

The integral curves of the vector field ∇*S*(**r**) are the geometrical optics rays associated with the eikonal function S(**r**), which measures the advance of the optical path length along the rays. We can create the tri-orthogonal curvilinear coordinates (*S*,*u*,*v*) where *u* and *v* are the coordinates defined on the surfaces *S*(**r**) = constant. This coordinate system is used to express the transport equations in the Geometrical Optics approximation, and has been specifically used in Electromagnetic Optics by Stavroudis [8]. Additionally, it is also commonly applied in the method of characteristics for solving the time-dependent wave equation. In these new coordinates *S = S*(**r**), *u = u*(**r**), *v* = *v*(**r**), Eq. (1) becomes:

Here *U* = *U*(*S*,*u*,*v*) and *n* = *n*(*S*,*u*,*v*). If there were solutions of *U* depending on *S* only, then (3) reduces to:

For this to be possible, it must be that

Given this, the field *U*(*S*) will be a solution of the following ordinary differential equation:

Note that since ${\left(\nabla S\right)}^{2}={n}^{2}(r)$we can state that there are solutions of *U* depending solely on *S* if (and only if) S fulfills:

Equivalently, Eq. (7) can be written as the following system of two third-order partial differential equations:

Note that any solution of Eq. (8) [or equivalently Eq. (7)] provides wavefronts for which *U* = *U*(*S*) exists and the refractive index is calculated as$n(r)=\left|\nabla S\right|$. The energy flux of the scalar wave field is given by [9]:

Here the asterisk denotes complex conjugation and the subindex denotes the derivative with respect to *S*. Therefore, in these fields **I** is tangent to the rays (as $\nabla S$ is).

As can be easily checked, Eq. (7) is also valid for two dimensions, *i.e.*, when the previous reasoning is followed with the Helmholtz Eq. (1) in two dimensions.

#### 2.2 Solutions of Eq. (7)

We will prove next that a function *S* is a solution of functional Eq. (7) in the domain *D* if and only if there is a real increasing function *V*(*S*) which is harmonic in *D*, *i.e.*, whose Laplacian is zero (∆*V* = 0):

The function *F*(*S*) of Eq. (7) is thus

Since *V _{S}* is the argument of a logarithm, it must be positive, which implies that

*V*(

*S*) is monotonically increasing (and therefore invertible).

This solution gives a constructive method to calculate the refractive index distributions that propagate fields *U* exactly on eikonals. The design method comprises the following steps:

*U*(

**r**) can finally be obtained by solving Eq. (6) for

*U*(

*S*), and then calculating

*U*(

*S*(

*V*(

**r**))).

#### 2.3. Equivalent statements

We have found a scalar wave field *U*(**r**) solution of the Helmholtz equation (in both two and three dimensions) that propagates with the help of two associated functions *S*(**r**) and *V*(**r**), which are solutions of the eikonal equation and Laplace equations, respectively. The three functions are functionally interrelated, that is, *U = U*(*S*), *S = S*(*V*), and thus *U = U*(*V*). It is remarkable that the problem could have been stated in two alternative ways, different from that in Section 2.1, and the results would have been the same. These two alternative statements and their solutions are:

- a) Find a scalar wave field
*U*(**r**) particular solution of the Helmholtz equation which is a function*U = U*(*V*) of a solution*V*(**r**) of Laplace equation. The necessary and sufficient condition is that the solution takes the form*V = V*(*S*), where*S*(**r**) is a solution of the eikonal equation. - b) Find a particular solution
*S*(**r**) of the eikonal equation which is a function*S = S*(*V*) of a solution*V*(**r**) of Laplace equation. When that solution is found, then a field*U*(**r**) is a solution of the Helmholtz equation of the form*U = U*(*S*) and can be found by solving Eq. (6).

The formulation has provided us a tool for the design of refractive-index distributions that propagate Helmholtz scalar wave fields. In Sections 3 and 4 we will apply this tool for perfect focusing of waves from a point source in two and three-dimensions, respectively.

## 3. Perfect focusing in two dimensions

#### 3.1 The cylindrical Maxwell fish-eye lens

Consider first the cylindrical Maxwell fish-eye lens case with refractive index distribution:

Here *L* and *a* are positive real constants and *ρ*, *z* are cylindrical coordinates (*ρ*
^{2} = *x*
^{2} + *y*
^{2}). As mentioned in the introduction, the Maxwell fish-eye lens is an Absolute Instrument in geometrical optics, and thus the rays contained in planes *z* = *z*
_{0} passing through any point **P** will focus again at its conjugate point **Q**. The position of **Q** in the x-y plane can be obtained as the transformation from **P** by an inversion with respect to the circle of radius *a* centered at **O,** followed by central symmetry with respect to **O**. Therefore, **P**·**Q** = −*a*
^{2}.

Leonhardt has recently proven in [5] that such ideality for focusing rays of the cylindrical Maxwell fish-eye lens also occurs for 2D Helmholtz fields. Leonhardt solved it by using the classical stereographic projection between the sphere and the plane, which Luneburg used as well to explain the ideality of the Maxwell fish-eye lens in ray optics [10]. We will confirm Leonhardt’s result here by using the approach of the previous section: showing that the function *S* of the rays emitted from **P** in the cylindrical Maxwell fish-eye lens fulfills Eq. (5).

Without loss of generality (due to the cylindrical symmetry), let us consider **P** = *a*(−cosh α − sinh α,0,0) and **Q** = *a*(cosh α − sinh α,0,0) for a given real constant *α*. Note that **P**·**Q** = −*a*
^{2}cosh^{2}α + *a*
^{2}sinh^{2}α = −*a*
^{2}. In the case α = 0, **P** = (−*a*,0,0) and **Q** = (*a*,0,0).

Apply the following change of variables to the following decentered cylindrical bipolar coordinates (σ, τ, *z*) [11]:

The cross section of surfaces of constant *σ* and constant *τ* are circles that intersect at right angles. As shown in Fig. 1
, curves of constant *σ* (in blue) correspond to circles that intersect at the two points **P** and **Q**, while the curves of constant *τ* (in red) are non-intersecting circles of different radii that surround the points **P** and **Q**.

In this coordinate system, the gradient and Laplacian operators take the forms [11]:

The refractive index distribution of the Maxwell Fish Eye [Eq. (12)] can be written in these coordinates (σ, τ, *z*) as:

We can confirm that the *τ* isosurfaces coincide with Geometrical wavefronts by the fact that the eikonal equation (∇*S*)^{2} = *n*
^{2} has particular solutions depending only on τ, *i.e.*, *S* = *S*(*τ*), which are calculated with Eq. (14) and Eq. (16):

This utilizes the principal value of the arctangent function. Choose the solution with positive sign (so the optical path length *S* increases with *τ*) and *C* = 0. The minimum value of *S*(*τ*) is reached at the **P** where *τ→ −∞* and *S*(*−∞*) = 0, while the maximum value is achieved at the **Q** where *τ→∞* and *S*(*∞*) = *πL*. Therefore 0 ≤ *S* ≤ *π L*.

We can invert *S*(*τ*) to obtain:

By (15):

Then we can easily check that Eq. (7) is fulfilled, since:

This uses Eq. (18). Therefore, there is a solution of the field *U* which is function of S only, and it fulfils the ordinary differential Eq. (6):

Using the change of variables $p=\mathrm{tan}\left(\frac{S}{2L}\right)$(*i.e.*, $p=\mathrm{exp}(\tau -\alpha )$), we obtain:

This is Eq. (9) of reference [5], with the wavenumber scaled by the factor *L* (which was normalized to 1 in that reference). The resolution was described in detail there but is not pertinent and so will not be continued here.

*V*(*S*) is obtained by integrating Eq. (11) with *F*(*S*), given by Eq. (20) as:

Here the last equality is valid for *A* = 1/*L* and *B* = *α* (see Eq. (18)). This indicates that the function *τ*(*x*,*y*) is harmonic, which is in fact well known, since the mapping of the planes x-y to τ-σ produced by the bipolar coordinate transformation (13) is conformal [11], and therefore both *τ*(*x*,*y*) and *σ*(*x*,*y*) are harmonic functions.

#### 3.2 The general cylindrical case: The “2D dipole lens”

The method of section 1.2 can find a more general solution of a device that focuses perfectly the scalar waves emitted by a point source a **P** at a point drain at **Q**. For that task, consider again the following harmonic function in R^{3}, valid except at the lines x = -*a*; y = 0 and x = -*a*; y = 0

This harmonic function is (up to a multiplicative constant) equal to the electrostatic potential created by an electric dipole with cylindrical symmetry along the z axis, *i.e.*, two line charges of equal magnitude but opposite sign crossing the z = 0 plane on the points **P** = (-*a*,0,0) and **Q** = (*a*,0,0). The function (24) coincides with function *τ*(*x*,*y*) of the previous section for the case *α* = 0 (as can be checked by inverting the mapping (13)), which is well known in electrostatic theory [11].

According to the method described in section 1.2, we can select an arbitrary monotonic function *S*(*V*). We will call this family of solutions the “2D dipole lens”. Different functions *S*(*V*) lead to different refractive index distributions, but all of them have the same wavefronts and rays defined by the coordinate lines of the bipolar reference system (13). One distinguished solution of this family is the cylindrical Maxwell fish-eye lens discussed in the previous section. In general, the refractive index distribution will not be a function of *ρ* only, as was the case for the cylindrical Maxwell fish-eye lens. Also unlike the Maxwell fish-eye lens, a general 2D dipole lens will not produce perfect focusing of the cylindrical scalar waves emitted by a point source different than **P** or **Q** (which are focused one into the other). Therefore, no perfect imaging for a volumetric region is expected in the general case.

In the selection of *S*(*V*) care must be taken with its asymptotic values when |*V|*→∞. This care is needed for the resulting refractive index (given by |∇*S*|) to be bounded, considering that *V* is unbounded around **P** and **Q**. In that case we must choose that when |*V|*→∞, *i.e.* for **r** close to **P** (or **Q)**:

Here r denotes the distance from **r** to **P** (or **Q**).

## 4. Perfect focusing in three dimensions

Note that, in order to simplify the nomenclature, in this section the points **P** and **Q** are placed on the *z* axis, whereas in the previous section they were on the *x* axis.

#### 4.1 The spherical Maxwell fish-eye lens

In the spherical Maxwell fish-eye lens, the rays trajectories (and their orthogonal surfaces) from **P** = *a*(0,0,cosh α − sinh α) and **Q** = *a*(0,0,−cosh α − sinh α) define the bispherical coordinate system (σ,τ,φ) given by [11]:

In this coordinate system, the eikonal equation takes the form:

Here the refractive index of the spherical Maxwell fish-eye lens is written in the (σ,τ,φ) coordinates as:

The particular solution *S* of Eq. (27) of the function *S* = *S*(*τ*) coincides with the cylindrical Maxwell fish-eye lens [Eq. (17)], as was expected. Using the expression for the Laplacian in the (σ,τ,φ) coordinates [11] gives:

From (28) and (29) it can be seen that (except in the trivial asymptotic cases *α*→ ± ∞, in which **P** or **Q** are at the origin) *S* does not fulfill Eq. (7) for any value of *α*, which means that there is no solution of the field *U* that is a function of *S* only, as occurred in 2D. This means that the spherical Maxwell fish eye lens does not belong to the class of media discussed here, but this does not mean that it cannot perfectly focus Helmholtz scalar wave fields of a point source **P** onto its image point **Q** in three dimensions (that ability or inability is has still to be proved). Leonhardt and Philbin, however, have proven very recently [6] that the spherical Maxwell fish eye lens can perfectly focus electromagnetic waves (not Helmholtz scalar wave fields) in three dimensions if the medium is impedance-matched (*n* = *ε* = *μ*).

We will see next that the “dipole lens” concept introduced in section 3.1 can be analogously defined in 3D, and it will lead to refractive index distributions which are not spherical symmetric but that produce the perfect focusing of the Helmholtz scalar wave field from **P** to **Q** in three dimensions.

#### 4.2 The general rotational case: The “3D dipole lens”

Consider the following harmonic function in R^{3} except at the points (*x*,*y*,*z*) given by **P** = (0,0,-*a*) and **Q** = (0,0,*a*) (so that *D* = R^{3}\{*P*,*Q*})

Here *ρ*, *z* are again cylindrical coordinates (*ρ*
^{2} = *x*
^{2} + *y*
^{2}). This harmonic function is (up to a multiplicative constant) equal to the electrostatic potential created by an electric dipole (two point charges of equal magnitude but opposite sign) located on the points **P** and **Q**. That potential distribution has rotational symmetry with respect to the z axis, and the cross sections of the *V* = constant surfaces and electrostatic field lines are shown in Fig. 2
(see also an animation in [13]). In our case, the *V* = constant surfaces are the wavefronts and the electrostatic field lines correspond to the ray trajectories and energy flux lines. Note that both families of curves are not circumferences, as occurred in the 2D case. Particularly, the *V* = constant surfaces are octics (8-degree algebraic curves) that belong to the family of Generalised Cayley’s ovals [14].

As in the 2D case, for the refractive index *n* = |∇*S*| to be bounded, we must select *S*(*V*) with the appropriate asymptotic behavior when |*V|*→∞, which in this 3D case is:

This is achieved for instance by the function:

Here again the principal value of the function arctan is used. The minimum value of *S*(*ρ*,*z*) is only reached at **P**, *S*(0,–*a*) *=* 0, while the maximum value is achieved at **Q**, *S*(0,*a*) *= πL*. Therefore 0 ≤ *S* ≤ *π L*.

The refractive index distribution can be computed as:

Here $s=+\sqrt{{a}^{4}+2{a}^{2}\left({\rho}^{2}-{z}^{2}\right)+{\left({\rho}^{2}+{z}^{2}\right)}^{2}}$. Such a rotationally symmetric distribution is even in *z* (as expected by the symmetries of *V*(**r**) and of *S*(*V*)) and is shown in Fig. 3
. At **P** and **Q**, $n=L/(\sqrt{2}a)$; at (*ρ*,*z*) = (0,0), $n=2L/a$, which is the maximum value of n. For large values of |**r**|, *n*(**r**) is smaller than one (as also occurs in the Maxwell fish eye lens).

In order to calculate the field *U*(*S*), first calculate *F*(*S*) using Eq. (11) and (32):

The field *U*(*S*) is a solution of Eq. (6), which in this case is:

This is close to Eq. (21) except for the factor of 2. The general solution of Eq. (35) has the following simple closed form [12]:

Here *C*
_{1} and *C*
_{2} are arbitrary constants and:

Note that that *κ* →*k* when *k*→∞. Select the particular solution *C*
_{2} = 1/(4*πa*); *C*
_{1} = i/(4*πa*), which leads to

This field (38) represents a wave propagating towards increasing values of *S*, that is, from point **P** to point **Q**, with phase advancing in proportion to *κ*. The energy flux vector of this field can be computed from Eq. (9) as:

The field can be expressed in cylindrical coordinates as *U*(*S*(*V*(*ρ*,*z*)), so that from Eqs. (38), (32) and (30):

Figure 4
shows the real part of *U*(*ρ*, *z*)exp(ω*t*) for ω*t* = 1.033*π*, *a* = 1, *L* = 1,*κ* = 9.5. Time evolution of the field is shown in Media 1.

The field diverges at both points **P** and **Q,** since at *S*(**P**) *=* 0 and *S*(**Q**) *= πL* the sine function vanishes in (38) and asymptotically behaves as:

This is because function *S* given by Eq. (32) approximates to:

Calling *r* the distance from a point (*ρ*, *z*) to **P** and *r’* the distance to **Q**, *i.e.*

**P**and

**Q**can be written as:

These expressions are identified as the field of a point source at **P** emitting radiation and point drain at **Q** receiving said radiation (because the energy flux (39) follows the ray trajectories), so no flux is radiated towards infinity.

## 5. Resolution limit

Using small spheres of radius *r* and *r’* around **P** and **Q**, respectively, we can confirm from (44) that:

**S**points towards the spheres exterior.

This confirms that the field (40) is a solution of the Helmholtz equation defined for the whole R^{3} as:

The addition of the Dirac-delta at **Q** is necessary for mathematical consistency. In the absence of the drain at **Q**, the field amplitude there will not diverge, because the wave will pass through the focus **Q**, expand from it and converge back upon **P**. We can calculate the field *W*(**r**) in absence of the drain at **Q** as:

Here *S(ρ,z)* is given by Eq. (30) and (32), (note that 0≤*S*≤*πL*). Around **Q** (*i.e.*, *ρ* = 0, *z* = *a*) we can calculate from Eq. (48) the asymptote:

This is non-singular at **Q** and has amplitude:

Around **P** (*i.e.*, *ρ* = 0, *z* = −*a*) the asymptote is:

As can be checked by direct calculation, the field *W*(**r**) of Eq. (48) is then the solution of the following Helmholtz equation:

Here *κL = κ*
_{1}
*a*≠*l*, with *l* a positive integer. At the frequencies for which *κL = κ*
_{1}
*a* = *l*, the system is an oscillator, so that Eq. (52) becomes homogeneous and the field amplitude at **P** equals that at **Q** of Eq. (50).

Equation (49) is the well known [11] spherical wave converging to and diverging from **Q**. It resembles an Airy pattern in three dimensions, with peak amplitude given by Eq. (50), and the radius of the sphere of the first null is:

Here *λ* is the wavelength of the field in vacuum, and *n*
_{max} = 2*L*/*a* is the maximum value of *n* [see Eq. (33)] obtained for (*ρ*,*z*) = (0,0). This seems to indicate that both the amplitude at focus and the resolution capacity of the 3D dipole lens (and the Maxwell fish eye in two dimensions [5]) is limited by the wavelength. This is not the case, however, when an absorber much smaller than the wavelength is placed at focus. It is expected that if a good absorber with absorption cross section much smaller than the wavelength (for instance, a McDonalds body in the case of acoustic waves [15]) is located at **Q**, it will present a high field amplitude on its surface. Using these subwavelength absorbers as detectors, it should be possible to resolve two point sources separated by a distance much smaller than the wavelength.

## 6. Conclusions

We have found non-trivial isotropic non-homogenous refractive-index distributions that propagate Helmholtz fields *U*(**r**) exactly on eikonals in that medium, in the sense that that *U*(**r**) = *U*(*S*(**r**)), where *S*(**r**) is a particular solution of the eikonal equation in that medium.

As an example, we have found a (positive) refractive index distribution *n*(**r**), the “3D dipole lens”, which in three dimensions perfectly focuses a point source at **P,** emitting a Helmholtz scalar wave field, onto a point drain at **Q**. To our knowledge this is the first example of such a system with a positive refractive index. Apart from the difficulty of making the 3D dipole lens (because eventually *n*(**r**) goes below 1), they have a clear theoretical interest. Perfect focusing of Helmholtz fields in 3D is particularly relevant to acoustics (for example, to ultrasound imaging), turning the sound waves of an explosion into an implosion.

The cylindrical Maxwell fish eye lens in two dimensions has shown to be a particular example of the “2D dipole lens”, although it has the very remarkable property that produces the perfect focusing for the point source **P** located in any arbitrary position, which makes it a perfect imaging device in two dimensions [5]. We have also shown that the spherical Maxwell fish eye lens does not belong to the family of refractive indices discussed here, that is, there is no solution of the Helmholtz equation in three dimensions that propagates exactly on eikonals in the spherical Maxwell fish eye medium.

The constructive procedure developed here opens the possibility of designing other interesting novel refractive-index distributions that propagate fields exactly on eikonals. For instance, perfect scalar wave focusing from a point source into multiple point drains (or from multiple point source into a single point drain) is expected when the harmonic function *V*(**r**) is selected as that produced by several electric charges.

## Acknowledgments

The authors thank the Spanish Ministries MCEI (Consolider program CSD2008-00066, DEFFIO: TEC2008-03773), MITYC (OSV: TSI-02303-2008-52), and the Madrid Regional Government (LED-TV: 130/2008 TIC, ABL: PIE/466/2009, F3: PIE/469/2009 and CAM/UPM-145/Q060910-103) for the support given in the preparation of the present work. The authors also thank Jesús López for creating the video and Bill Parkyn for editing the manuscript.

## References and links

**1. **M. Born, and E. Wolf, *Principles of Optics*, (Pergamon, Oxford, 1975, 5th ed)

**2. **S. Cornbleet, *Microwave and Geometrical Optics*, (Academic, London, 1994)

**3. **J. C. Miñano, “Perfect imaging in a homogeneous threedimensional region,” Opt. Express **14**(21), 9627–9635 (2006), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-21-9627. [CrossRef]

**4. **J. B. Pendry, “Negative refraction makes a perfect lens,” Phys. Rev. Lett. **85**(18), 3966–3969 (2000). [CrossRef]

**5. **U. Leonhardt, “Perfect imaging without negative refraction,” N. J. Phys. **11**(9), 093040 (2009). [CrossRef]

**6. **U. Leonhardt and T. G. Philbin, “Perfect imaging with positive refraction in three dimensions,” Phys. Rev. A **81**(1), 011804 (2010). [CrossRef]

**7. **M. A. Alonso and G. Forbes, “Stable aggregates of flexible elements give a stronger link between rays and waves,” Opt. Express **10**(16), 728–739 (2002), http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-10-16-728. [PubMed]

**8. **O. N. Stavroudis, *The Mathematics of Geometrical and Physical Optics*, (Willey-VCH, Weinheim, 2006)

**9. **Yu. A. Kratsov, and Yu. I. Orlov, *Geometrical Optics of Inhomogeneous Media,* p.23, Springer Verlag, Berlin Heidelberg, 1990.

**10. **R. K. Luneburg, *Mathematical Theory of Optics*, (University of California Press, Los Angeles 1964)

**11. **P. M. Morse, and H. Feshbach, *Methods of Theoretical Physics* (New York, McGraw-Hill, 1953)

**12. **A. D. Polyanin, and V. F. Zaitsev, *Handbook of Exact Solutions for Ordinary Differential Equations*, 2nd Edition, Chapman & Hall/CRC, Boca Raton, 2003. It also coincides with the change of variables *S*’ = *S/L* + *π*/2 in the equation described in http://eqworld.ipmnet.ru/en/solutions/ode/ode0235.pdf.

**13. **http://www.youtube.com/watch?v=bG9XSY8i_q8

**14. **http://www.mathcurve.com/courbes2d/cayleyovale/cayleyovale.shtml

**15. **Yu. I. Bobrovnitskiĭ, “Impedance theory of sound absorption: The best absorber and the black body,” Acoust. Phys. **52**(6), 638–647 (2006). [CrossRef]