Flat lens without optical axis: Theory of imaging

We derive a general theory for imaging by a flat lens without optical axis. We show that the condition for imaging requires a material having elliptic dispersion relations with negative group refraction, equivalent to an effective anisotropic refractive index n(theta). Imaging can be achieved with both negative (n<0) and positive (n>0) refractive indices. The Veselago-Pendry lens is a special case with isotropic negative refractive index of n(theta)=-1. Realizations of the imaging conditions using anisotropic media and inhomogeneous media, particularly photonic crystals, are discussed. Numerical examples of imaging and requirements for sub-wavelength imaging are also presented.

The concept of negative refraction has led to new fundamental approaches as well as applications in optics [6]. In 1968 Veselago [7] pointed out that for a material with refractive index n = −1 a flat surface would focus light. A decade later, Silin [8] discussed a more general case and obtained the lens equation of a flat slab with negative elliptic dispersion. Pendry's recent analysis [9] demonstrating the possibility of sub-wavelength resolution with a flat slab of such materials, as well as the experimental realization of the so-called left-handed materials using composite media [10] and photonic crystals (PhCs) [11,12] led to renewed interest in the unique electromagnetic properties of these artificial materials. The unique property of these flat lenses is the lack of optical axis. This type of flat lenses can be realized in a PhC using negative refraction [13] and has been demonstrated in microwave experiment [14]. Super-resolution imaging through single negative index media are also studied [15,16,17,18]. There is however no complete theory of imaging by a flat lens which properly describes the various features observed in the experiments and in numerical simulations.
In this paper, we present a general theory of imaging by a flat lens without optical axis, resulting in a proposal for a new material which is translationally invariant along the surface and has anisotropic refractive index n(θ). Defining an optical phase condition for imaging, we show that the condition requires a material having Ellip- * Electronic address: w.lu@neu.edu † Electronic address: s.sridhar@neu.edu tic dispersion relations with Negative Group Refraction (ENGRM) at the operating frequency. The ENGRM is defined by two parameters: the anisotropy parameter σ which is a measure of the ellipse eccentricity, and a phase factor κ which determines the center of the ellipse. Refraction laws for the wave vector and the group velocity for anisotropic media are derived. The theory shows that "perfect" images can be obtained, and imaging is possible with negative as well positive refractive indices. The required conditions and consequences of flat perfect and imperfect lenses made of ENGRMs are derived from a generalized Fermat's principle. Realization of flat lenses using PhCs is discussed. In real PhCs, σ is itself angledependent in most cases, leading to some limitations for image formation. Many of the proposals for PhC lenses are contained in the present theory. The theory is also applied to recent experiments using PhCs and is shown to successfully describe some intriguing aspects of the data. Numerical simulations are presented that describe visual details of the image formation in PhCs. We first consider the case of a single interface which is along the x-axis and at z = 0 as shown in Fig. 1(a). The wave vector k of an incident plane wave from vacuum towards the interface makes an angle θ with the z-axis.
Here we consider the case that for any incident k, there is only one refractive wave vector k r . Continuity along the interface requires that k rx = k x . A point source is at x = 0 and z = −u with u > 0. We now consider another point in the ENGRM at x = ∆x and z = v with v > 0. The phase difference between these two points for each incident k is Φ = Φ z + Φ x with Φ z = k z u+ k rz v and Φ x = k x ∆x. According to Fermat's principle, the formation of an image would require that the total phase Φ be stationary. This restriction can be relaxed for a flat interface. Since the surface is flat, any point of a finite-size source will be imaged to a point with the same x-coordinate. Thus ∆x = 0 and Φ x = 0. One only needs to consider Φ z . A generalized Fermat's principle states that an image will be formed if the phase Φ z is stationary, dΦ z /dk z = 0. The lens equation for a single interface is with a material constant This constant σ should be positive and independent of the incident angle for a focus without aberration. Thus the following rule for the wave vector refraction at the interface must be obeyed Here κ is the integration constant which is also an intrinsic property of the ENGRM. Since in the vacuum Note that this elliptic dispersion exists only at the operating frequency ω 0 [19]. In the neighborhood of this frequency ω ∼ ω 0 , (κ − k rz ) 2 /n 2 z + k 2 rx /n 2 x = ω 2 /c 2 with σ ≡ n z /n x ∼ σ. The deduction of Eq. (3) from an elliptic EFS requires negative group refraction together with the causality condition ∂ω/∂k rz > 0 that the energy must flow away from the interface which is illustrated in Fig. 2. The ray vector [20] which represents the direction of group velocity in a medium with elliptic EFS is s = (−σ sin αx + cos αẑ)/n xσ while k r − κẑ =n x (sin αx −σ cos αẑ) with sin α = n −1 x sin θ and θ the incident angle in the vacuum.
The Snell's law for wave vector loses its meaning in this anisotropic medium. However there is a law for group refraction, which is tan β = −σ tan α = −σ(n 2 x − sin 2 θ) −1/2 sin θ. At the operating frequency ω 0 of the flat lens, n x = 1,σ(ω 0 ) = σ, the group refraction law is simply Thus −σ can be regarded as the effective refractive index. for any σ as defined by Eq. (2). The proof is the following,

Note that this refraction law is very general and is valid
The space inside the medium can be considered to be "optically stretched" in the z-direction by a factor of σ −1 . The waves inside the medium are represented by ψ(r ′ ) = exp(ik x x − ik z z ′ + iκ ′ z ′ ) with z ′ = σz and κ ′ = κ/σ. Thus in the case that κ = 0, the waves at this frequency will experience a periodicity in the z-direction with length scale 2π/|κ| since the waves are modulated by exp(iκz). This periodicity is the effective periodicity and may not be the actual physical periodicity such as the lattice spacing of a PhC.
We now consider the imaging of an ENGRM slab with thickness d as shown in Fig. 1(b). The first surface is at the origin and the second at z = d. A point source is placed at z = −u. We consider another point outside the ENGRM slab at z = d + v. If the value of v satisfies the following lens equation the phase difference Φ z = k z (u + v) + k rz d = κd is stationary and an image will be formed without aberration. Note that for an ordinary thin curved lens with focal length f , the lens equation is With given value of σ and slab thickness d, the distance between the object and the image outside the slab is fixed, Once the location of the object is fixed, the position of the image is also fixed, no matter where the slab is placed! We remark that a flat slab can focus with negative as well as with positive wave vector refraction. The necessary condition for flat lens imaging is negative group refraction. In general, there is not much meaning of Snell's law for wave vectors for anisotropic media. Only the re-fraction law for the direction of group velocity Eq. (6) is meaningful and does correspond to ray diagrams. The Veselago-Pendry lens is a special case of our flat lens with σ = 1, κ = 0 and is the only case that the ray diagram is applicable for the wave vectors.
The lensmaker's formula Eq. (1) or (7) only provides a necessary condition for imaging. An additional constraint is required for the formation of a "perfect" image, viz. the surface reflection coefficient r should vanish, r = 0, for far field and should diverge, r → ∞, for near field. These set the conditions for a perfect flat lens without optical axis. The requirement of no reflection r = 0 is not essential for a slab to focus far field. The principal conditions for the flat lens to focus light are that Eq. (3) and |r| ≪ 1 should be satisfied for all or at least a large range of incident angles. The presence of reflection will make the image dim and may give rise to multiple images. The indefinite indices medium [21] can have an elliptic dispersion relation, but in general will not satisfy the flat lens imaging conditions.
The field inside the ENGRM flat lens can also be described by a partial differential equation [σ −2 (iκ + ∂ z ) 2 + ∂ 2 x + (ω/c) 2 ]Ψ = 0. The medium can be regarded as a distorted space with the metric ds 2 = dx 2 + σ 2 dz 2 . Thus there is an analogy with gravitational lensing where refraction occurs due to the warping of space caused by general relativistic effects, such as in the vicinity of a massive object. It is worth noting that general relativistic effects are observable due to optical refraction.
In the rest of the paper, we will focus on S-polarized electromagnetic waves with the electric field in the ydirection. In general, the magnetic permeability µ of the ENGRM should be a tensor, and its relationship to the other material parameters σ and κ can be deduced. Here we assume that effective indices ǫ and µ can be used to describe the medium of interest. Due to its symmetry, the ellipsoid can always be transformed to its principal axes [2]. For the S-polarized waves, the H field is in the xzplane, only µ x and µ z will enter our discussion. Furthermore, µ x plays a more active role than µ z since the surface reflection coefficient is r = (µ x k z −k rz )/(µ x k z +k rz ). Thus to achieve perfect imaging for far field, the effective permeability of ENGRM should be Only in the case κ = 0 that µ x = −σ is isotropic while for κ = 0, µ x has to be anisotropic and diverges at k x = k 0 . For κ = 0, any finite µ x at k x = k 0 results in reduced transmission. Notice that in PhCs including metamaterials, the effective µ x should be anisotropic in general. To see how the propagating waves are focused, we consider a S-polarized point source −(i/4)H field is Note that k z = k 0 cos θ and k x = k 0 sin θ. The image outside the slab is located at z = (1 + σ)d − u. One can see that except for a global phase exp(iκd), the far field of the image is exactly that of the source. The far field inside the slab is for 0 ≤ z ≤ d. Another perfect image is formed inside the slab at z = u/σ. The first surface acts as a mirror with an extra phase exp(iκz). The space inside the ENGRM slab is optically stretched by a factor of σ −1 . A flat lens imaging is shown in Fig. 3. Evanescent waves and sub-wavelength imaging For evanescent waves k x > k 0 , the complex extrapolation of Eq. (4) gives In order to amplify evanescent waves, a singularity must exist in T to compensate to certain extent the decay of evanescent waves in the vacuum exp[−q(u + v)]. In the case of single-interface resonance [9,24], r → ∞, T exp[−q(u + v)] = exp(iκd), all the evanescent waves will be amplified only if The images both inside and outside have the same subwavelength features of the source. If µ x couldn't take the above form, sub-wavelength imaging is still possible. At the so-called overall resonance condition [24], 1 − r 2 exp(2ik rz d) = 0, thus r = ± exp(−ik rz d), one gets For κ = 0, µ 0 x and µ ± x are all complex with negative real parts. One notices that µ + x (q) is flatter than µ − x (q) as functions of q. In terms of effective indices, the existence of surface modes ω 0 (k x ) requires that ℜµ x be negative for S-polarization. To amplify evanescent waves, the ω 0 (k x ) curve must be very flat [24]. This is equivalent to say that the curve µ + x (q) should be flat so that the contribution of evanescent waves are constructive for certain window of q. The closer µ x is to µ + x (q), the more sub-wavelength features the images will have. Note that for large κd, the evanescent waves are located in the vicinity of the interfaces for constant µ x .
Realization in real materials We note that no known natural material is found to have negative dispersion. Negative group refraction can be achieved in periodic or quasi-periodic media with nonzero κ [22]. For most media, Eq. (8) and (12) are unlikely to be fulfilled, thus a perfect flat lens is unattainable. Realization using real materials is discussed next.
Ordinary material may not have the EFS described by Eq. (4) for all k x . For an uniaxial crystal, since the group refraction is positive, thus k rz = κ − σk z with κ = 0 and σ < 0. There will be no focus though negative refraction for certain incident angles can be easily achieved [23]. Instead a virtual image will be formed satisfying the same lens equation. However, a PhC could have both an elliptic EFS and negative group refraction for certain windows of k x at certain frequency. Here we use some general features of PhCs to explain the imaging mechanism. To this end, we consider the first band of a square lattice PhC of unit lattice spacing. For small k x along the ΓM direction, the dispersion can be approximated elliptically as k rz ≃ κ − σ 0 k z with appropriate constants κ and σ 0 (see the dashed line in Fig. 4). At ω 0 = ω u , the upper limit of all-angle negative refraction (AANR) in the first band [13], κ = √ 2π − k 0 (1 − σ 0 ). Since for AANR with ω 0 away from ω u , σ 0 will be smaller, we only concentrate here on the discussion for ω 0 = ω u . To enhance the transmission, the slab thickness can be selected to satisfy the Fabry-Perot resonance condition (κ − σ 0 ω/c)d = 0 (mod 2π). The slab is on resonance only for k x → 0. For large k x , k rz ≃ √ 2π − σ −1/2 0 k z due to the C 4v symmetry of the EFS. A desired large transmission would require µ x to diverge. Thus the image location is bounded σ 0 ≤ (u + v)/d ≤ σ −1/2 0 and the Fourier components of the object with k x → k 0 will be partially lost. This leads to a so-called self-collimation effect for small σ 0 . Far field images both inside and outside the slab with reasonably good quality will be formed. However, the image inside the PhC is stretched and modulated due to the partial transmission besides the Bloch wave modulation. Even with the removal of Bloch wave modulation, the image inside a PhC will be visible only for large thickness. For evanescent waves at ω = ω u with k x → k 0 , one has ℜµ 0 To have substantial sub-wavelength feature of the image, the real part of µ x must take large negative value for small σ 0 . An example is shown in Fig. 4, which is very similar to Fig. 5 in Ref. [13]. The image is confined in the vicinity of the second surface of the PhC slab due to a small value of σ 0 . In this paper, we have discussed the group refraction of an anisotropic medium characterized by two materials parameters σ and κ. When σ is a positive constant, the anisotropic medium is an ENGRM and a flat slab of such medium can be used as a focal lens without optical axis and leading to images free of aberration. Our theory of flat lens imaging is a generalization of the Vesalago-Pendry perfect lens and beyond Silin's formula. The theory is valid for both real and virtual images. This flat lens theory leads to a clear understanding of negative group refraction and flat lens imaging in electromagnetism [13,14], acoustics [25], and electron waves.
Numerical simulations of homogeneous ENGRMs and of PhCs were carried out supporting our theory.
The lack of optical axis for the flat lens has very broad applications and confers important advantages in optics. Clearly there are no aperture restrictions. There are two characteristics that might be viewed as limitations. The lensmaker's formula clearly dictates upper bounds to u, v < σd for real image formation. The other limitation is that magnification is always unity. For any medium the lens properties σ, κ, and also the working frequency ω 0 can be obtained by inspection of the EFS. To have high quality image, σ should be constant and close to unity. The present theory can be used to design tailor-made flat lenses. Extension of our theory to three-dimensions is straightforward.
Work supported by the National Science Foundation and the Air Force Research Laboratories, Hanscom, MA.