Newtonian photorealistic ray tracing of grating cloaks and correlation-function-based cloaking-quality assessment

doi:10.1364/OE.19.006078

« Show journal navigation

Newtonian photorealistic ray tracing of grating cloaks and correlation-function-based cloaking-quality assessment

Jad C. Halimeh, Roman Schmied, and Martin Wegener »View Author Affiliations

Optics Express, Vol. 19, Issue 7, pp. 6078-6092 (2011)
http://dx.doi.org/10.1364/OE.19.006078

View Full Text Article

Acrobat PDF (2105 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (31)
Cited By
Figures (7)
Metrics

Abstract

Grating cloaks are a variation of dielectric carpet (or ground-plane) cloaks. The latter were introduced by Li and Pendry. In contrast to the numerical work involved in the quasi-conformal carpet cloak, the refractive-index profile of a conformal grating cloak follows a closed and exact analytical form. We have previously mentioned that finite-size conformal grating cloaks may exhibit better cloaking than usual finite-size carpet cloaks. In this paper, we directly visualize their performance using photorealistic ray-tracing simulations. We employ a Newtonian approach that is advantageous compared to conventional ray tracing based on Snell’s law. Furthermore, we quantify the achieved cloaking quality by computing the cross-correlations of rendered images. The cross-correlations for the grating cloak are much closer to 100% (i.e., ideal) than those for the Gaussian carpet cloak.

1. Introduction

The concept of transformation optics [1

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 (2006). [CrossRef] [PubMed]

–5

M. Wegener and S. Linden, “Shaping optical space with metamaterials,” Phys. Today 63(10), 32–36 (2010). [CrossRef]

] has firmly established that arbitrary objects can be made to appear invisible by surrounding them with an appropriately shaped inhomogeneous and generally anisotropic optical magneto-dielectric structure. To bring invisibility cloaking “down to earth”, researchers now aim at obtaining good cloaking action for finite-size structures which are as small as possible compared to the object to be hidden, while obeying the following constraints: (i) avoid singularities of the optical properties, (ii) avoid metal structures necessary for magnetic resonances if possible (because of losses at optical frequencies), (iii) avoid very large anisotropies if possible (because this again requires resonances that lead to losses), and (iv) maintain operation over a wide frequency band.

The so-called carpet cloak [6

J. Li and J. B. Pendry, “Hiding under the carpet: a new strategy for cloaking,” Phys. Rev. Lett. 101(20), 203901 (2008). [CrossRef] [PubMed]

] introduced by Li and Pendry has already come amazingly close to this aim. Corresponding experiments have been presented in Refs [7

R. Liu, C. Ji, J. J. Mock, J. Y. Chin, T. J. Cui, and D. R. Smith, “Broadband ground-plane cloak,” Science 323(5912), 366–369 (2009). [CrossRef] [PubMed]

–12

H. F. Ma and T. J. Cui, “Three-dimensional broadband ground-plane cloak made of metamaterials,” Nat. Commun. 1(3), 1–6 (2010). [CrossRef]

]. However, the carpet cloak does suffer from certain limitations, such as a pronounced lateral beam displacement [13

B. Zhang, T. Chan, and B.-I. Wu, “Lateral shift makes a ground-plane cloak detectable,” Phys. Rev. Lett. 104(23), 233903 (2010). [CrossRef] [PubMed]

]. This artifact – which was already visible in previous publications [6

J. Li and J. B. Pendry, “Hiding under the carpet: a new strategy for cloaking,” Phys. Rev. Lett. 101(20), 203901 (2008). [CrossRef] [PubMed]

,14

J. C. Halimeh, T. Ergin, J. Mueller, N. Stenger, and M. Wegener, “Photorealistic images of carpet cloaks,” Opt. Express 17(22), 19328–19336 (2009). [CrossRef] [PubMed]

] – makes the cloaking structure detectable. We have recently pointed out [15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

] that the decay of the refractive-index profile away from the hidden object is related to the spatial-frequency distribution of the bump in the carpet. In particular, bump profiles with a cut-off at low spatial frequencies lead to a much faster decay of the refractive-index profile. Subsequently, the lateral beam displacement disappears almost completely for a single non-zero spatial frequency component. We have called this configuration the grating cloak [15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

However, it is difficult to assess the effect of the remaining imperfections of the finite-size cloaking structure from the mere inspection of selected rays in two dimensions (as done in our Ref [15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

].). In contrast, photorealistic ray tracing [16

A. S. Glassner, An Introduction to Ray Tracing (Morgan Kaufmann, 1989).

–18

A. J. Danner, “Visualizing invisibility: metamaterials-based optical devices in natural environments,” Opt. Express 18(4), 3332–3337 (2010). [CrossRef] [PubMed]

] is a well-known rendering technique based on ray optics, allowing to simulate images as they would appear to the human eye, and thus providing a direct and intuitive impression of their performance. These images can also be the basis for a quantitative assessment of the cloaking quality using cross-correlation functions (see section 5).

In this paper, we present corresponding numerical calculations. To obtain fully converged results under these conditions, we employ Newtonian ray tracing in a dedicated home-built software program.

2. The grating cloak

The refractive-index distribution n of the grating cloak has been given [15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

] in closed form. Replacing the two-dimensional coordinates (u,v) in Eq. (10) in Ref [15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

]. by (x = u, y = v, z) in three dimensions, we have

n (x, y, z) = \frac{1}{| 1 + W_{0} (i c k e^{i k (x + i y)}) |} .

(1)

Here, W ₀ is the Lambert function (or product logarithm). The constant wave number k is given by k = 2π /Λ with the real-space grating period Λ; the complex-valued constant c is given by c = iA with real-valued A, where 2A is the peak-to-peak vertical grating amplitude (see Eq. (2)). Thus, numerical errors in determining the index profile (as, e.g., possible for the quasi-conformal carpet cloak requiring numerical minimization of the modified-Liao functional) do not play any role at all in our present work. We do not consider any frequency dependence of n.

The profile y(x) of the metal grating off of which light rays are reflected and which is to be cloaked by the refractive-index profile given by Eq. (1) has the shape of a trochoid given by the parametric form [15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

]

\begin{matrix} x (ξ) & = & ξ - A \sin (k ξ) \\ y (ξ) & = & A \cos (k ξ), \end{matrix}

(2)

where the parameter ξ (corresponding to x in Ref [15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

].) runs from + ∞ to -∞. The grating profile given by Eq. (2) is translationally invariant along the z-direction and can be brought into the explicit form

x (y) = \pm (\frac{Λ}{2 π} \arccos (\frac{y}{A}) - \sqrt{A^{2} - y^{2}}) + N Λ,

(3)

with integer N = -∞, …, + ∞ for the different grating maxima. The “-“ sign in Eq. (3) applies for the left-hand-side slope of each maximum, the “+” sign for the right-hand-side slope. In the shallow-grating limit defined by 2π A/Λ <<1, hence x ≈ξ in Eq. (2), one gets the approximate explicit form

y (x) \approx A \cos (k x) .

(4)

However, for the numerical evaluation of the ray trajectories we will rather use the exact forms of Eqs. (2) and (3) throughout this work. Figure 1(a) shows an example for the normalized grating-cloak model parameters Λ = 1 and A = 3/(16π)≈0.06 that lead to experimentally accessible refractive indices. In Figs. 3 -7 , one normalized unit corresponds to an actual 19 cm.

Fig. 1 (a) False-color representation of the normalized refractive-index profile n/n _min of the conformal grating cloak in the xy-plane according to Eq. (1). The normalized grating-cloak parameters are Λ = 1 and A = 3/(16π)≈0.06 (leading to n _min = 0.727). The normalized height of the cloak is equal to 1.1. In Figs. 3-7, one normalized unit corresponds to 19 cm. (b) Same, but for the Gaussian conformal map with w = 0.306 and h = 0.13 (leading to n _min = 0.873). The normalized height of the cloak is 1.1. Note that these parameters have been chosen such that n/n _min in (a) and (b) varies in the same interval [1.0, 2.2]. This aspect allows for comparing the two setups in a meaningful manner.

Fig. 3 Color images rendered by the Newtonian ray-tracing approach. The corresponding scenery is illustrated in Fig. 2. (a) Flat mirror. (b) Grating profile with Λ = 1 (or 18.9 cm) and A = 3/(16π) ≈0.06 (or 1.1 cm), without cloak. (c) As (b), but with cloak. The cloak height is 1.1 normalized units (or 20.8 cm). Its refractive-index profile is depicted in Fig. 1(a). Note that this profile contains index values below unity with minimum n _min = 0.728.

Fig. 7 As Fig. 6, but for a yet smaller height of the grating-cloak structure of 0.41 normalized units (or 7.8 cm).

Since the grating cloak cannot be directly compared to the familiar carpet cloak [6

J. Li and J. B. Pendry, “Hiding under the carpet: a new strategy for cloaking,” Phys. Rev. Lett. 101(20), 203901 (2008). [CrossRef] [PubMed]

], we also generate images for a single-bump Gaussian conformal cloak [15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

], which is more amenable to such a comparison. In the shallow-bump limit its bump shape is given by

y (x) \approx h e^{- {(x / w)}^{2}},

(5)

with the height h and the width w. However, for the actual calculations, we use the exact implicit form given in Ref [15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

]. Figure 1(b) shows an example for the normalized model parameters w = 0.306 and h = 0.13.

3. Newtonian photorealistic ray tracing

We have previously [14

J. C. Halimeh, T. Ergin, J. Mueller, N. Stenger, and M. Wegener, “Photorealistic images of carpet cloaks,” Opt. Express 17(22), 19328–19336 (2009). [CrossRef] [PubMed]

,19

T. Ergin, J. C. Halimeh, N. Stenger, and M. Wegener, “Optical microscopy of 3D carpet cloaks:ray-tracing calculations,” Opt. Express 18(19), 20535–20545 (2010). [CrossRef] [PubMed]

] employed ray tracing based on Snell’s law for sceneries involving dielectric cloaks. These results have been fully converged. However, later we noticed that certain viewing directions, such as a perpendicular view onto a cloaked bump, tend to converge only very slowly or not at all within accessible computation times. Sometimes, stripes appeared as artifacts in the rendered images. Thus, we started investigating alternative ray-tracing approaches.

It is well established [18

A. J. Danner, “Visualizing invisibility: metamaterials-based optical devices in natural environments,” Opt. Express 18(4), 3332–3337 (2010). [CrossRef] [PubMed]

,20

J. L. Synge, Geometrical Mechanics and De Broglie Waves (Cambridge U. Press, 1954).

–29

K. Niu, C. Song, and M.-L. Ge, “The geodesic form of light-ray trace in the inhomogeneous media,” Opt. Express 17(14), 11753–11767 (2009). [CrossRef] [PubMed]

] that ray optics in inhomogeneous optical media can be formulated in analogy to Hamiltonian mechanics. Here we use a simple and intuitive formulation that is equivalent to Newton’s second law.

To emphasize the close analogy to mechanics and to aid the derivation, we briefly review three important Eqs. (6)-(8) from classical mechanics. We start from Hamilton’s variation principle for the action integral for the (generalized) coordinates and the (generalized) velocities for physical trajectories

δ \int_{t_{1}}^{t_{2}} L (\vec{r}, \vec{v}, t) d t = 0.

(6)

As described in many mechanics textbooks, this allows for deriving the Euler-Lagrange equations for the components of the coordinate vector r_i and the components of the velocity vector v_i

\frac{\partial L}{\partial r_{i}} - \frac{d}{d t} \frac{\partial L}{\partial v_{i}} = 0.

(7)

Identifying the first term as the force component F_i , the second term as the temporal derivative of the momentum component, and additionally assuming a constant particle mass m, one immediately arrives at Newton’s second law for the acceleration

\frac{d \vec{v}}{d t} = \frac{d^{2} \vec{r}}{d t^{2}} = \frac{\vec{F}}{m} .

(8)

We can now simply repeat this procedure for light rays, where Hamilton’s principle (Eq. (6)) has to be replaced by Fermat’s principle, i.e., by

δ \int_{t_{1}}^{t_{2}} n (\vec{r}) | \vec{v} | d t = 0.

(9)

In analogy to Eq. (7), this leads to

| \vec{v} | \vec{\nabla} n - \frac{d}{d t} (n \frac{\vec{v}}{| \vec{v} |}) = 0,

(10)

and upon inserting the modulus of the local phase velocity of light v = c ₀/n, which is known a priori at each point in space, to the Newtonian equation of motion

\frac{d \vec{v}}{d t} = \frac{d^{2} \vec{r}}{d t^{2}} = \frac{| \vec{v} |^{2} \vec{\nabla} n - 2 (\vec{\nabla} n \cdot \vec{v}) \vec{v}}{n} .

(11)

The effective light-ray acceleration on the right-hand side of Eq. (11) vanishes in homogeneous media and thus outside the cloak. Inside the cloak, we solve this ordinary second-order differential equation (Eq. (11)) numerically for each ray. At the interface between the cloak and the surrounding air (or embedding medium), Snell’s law and the Fresnel equations correctly describe the refraction and reflection of the ray due to the discontinuity in the refractive index. We neglect rays which are Fresnel reflected more than once, because their intensity will be negligible in practice.

The refractive-index profile of the grating cloak is given in the closed mathematical form of Eq. (1). Evaluating the gradient yields

\vec{\nabla} n = (\begin{matrix} (\frac{2 Re (1 + W_{0}) Im (W_{0})}{| 1 + W_{0} |^{2}} - Im (W_{0})) \times \frac{k}{| 1 + W_{0} |^{3}} \\ (\frac{{Im}^{2} (W_{0}) - {Re}^{2} (1 + W_{0})}{| 1 + W_{0} |^{2}} + Re (1 + W_{0})) \times \frac{k}{| 1 + W_{0} |^{3}} \\ 0 \end{matrix}),

(12)

with the abbreviation for the Lambert function W ₀ (see Eq. (1))

W_{0} = W_{0} (i c k e^{i k (x + i y)}) .

(13)

Equation (12) can be inserted on the right-hand side of Eq. (11). We mention that, even though the refractive-index profiles used in this work are translationally-invariant in the z-direction, and thus the z-component of the gradient of n is zero throughout space, the ray acceleration in the z-direction in Eq. (11) does in general not vanish.

4. Rendered images

The scenery used for the following renderings is illustrated in Figs. 2(a) and 2(b). In essence, a model, a physicist, looks at herself in a curved mirror in front of her and observes the images depicted in Figs. 3-7. Technically, this means that a virtual point camera with a “human” focal field of view (i.e., full opening angles of 50° horizontally and 42° vertically) is located between her eyes, which are centered in the photograph in Fig. 2(a). The distance between model and mirror is d = 0.5 m (or 2.64 normalized units) – a typical distance for standing in front of your bathroom mirror. The model is 1.7 m (or 9.0 normalized units) tall. In the coordinate system shown in Figs. 1(a) and 1(b), the model is located at the top and looks downward onto the (a) grating or (b) Gaussian bump. For case (a), she sees three grating maxima within her vertical field of view (FOV): one at the center and two at the edges. For case (b), the single maximum is at the center of her FOV.

Fig. 2 Illustration of the scenery used for all photorealistic images rendered in this work. (a) This color image actually serves as the input for the following Newtonian ray-tracing calculations. (b) In these calculations, the model looks at the mirror in front of her (compare Fig. 1(a)) and observes the images shown in Figs. 3-7. The mirror (gray) is located at a distance of d = 0.5 m from the model, the cloak is shown in green.

Corresponding photorealistic images rendered by the Newtonian ray-tracing approach described in the preceding section are shown in Fig. 3 for the grating profile and in Fig. 4 for the conformal Gaussian map. In both figures, (a) corresponds to a flat mirror and (b) to the corrugated mirror – a grating for Fig. 3 and a single bump for Fig. 4.

Fig. 4 As Fig. 3, but for the Gaussian conformal map with w = 0.306 and h = 0.13. The corresponding refractive-index profile is depicted in Fig. 1(b). Note that this profile contains index values below unity with minimum n _min = 0.873.

We see that panels (b) are distorted versions of panels (a). In both panels (b), a maximum of the profile is located in the vertical middle of the image. The maxima act as convex mirrors. This makes the model appear right-side-up in this part of the image, but strongly compressed along the vertical direction. For the grating mirror in Fig. 3, the minima of the profile act as concave mirrors. The next grating maxima appear at the very top and at the very bottom of the image in Fig. 3(b), respectively (i.e., tan(42°/2) = Λ/d). Between maxima and minima, inflection points of the mirror profile occur. At each inflection point, the image turns from right-side-up to upside-down or vice versa. The combination of these aspects leads to highly complex and pronounced distortions with respect to the reference, Fig. 3(a).

In panels (c) of Figs. 3 and 4, a graded-index profile according to Fig. 1 is added onto the corrugated mirror profile. For an ideal cloak, panels (c) should be identical to panels (a), i.e., all the distortions in the respective panels (b) should be compensated for. This is indeed very nearly the case on the vertical line through the center of the images, but increasing deviations occur towards the sides. The grating cloak in Fig. 3 appears to perform better than the Gaussian conformal map in Fig. 4. We will quantify these aspects by evaluating cross-correlation coefficients of panels (a) and (c) in the following section.

Figures 3 and 4 remain unrealistic in that we have used the unreferenced refractive-index profile, i.e., the index of air is taken as unity and the indices of the cloaks do exhibit values below unity, with the minimum values given in Fig. 1. This aspect makes experimental realization difficult. If one divides all refractive indices by the minimum value n _min, the results remain unchanged. However, it is more realistic to assume that we keep an air refractive index of n = 1 (rather than n = 1/n _min>1). It has been pointed out repeatedly [6

J. Li and J. B. Pendry, “Hiding under the carpet: a new strategy for cloaking,” Phys. Rev. Lett. 101(20), 203901 (2008). [CrossRef] [PubMed]

–14

J. C. Halimeh, T. Ergin, J. Mueller, N. Stenger, and M. Wegener, “Photorealistic images of carpet cloaks,” Opt. Express 17(22), 19328–19336 (2009). [CrossRef] [PubMed]

] that this step introduces additional imperfections. These imperfections are visualized for the grating cloak in Fig. 5 (same cloak height as in Fig. 3, i.e., 1.1 normalized units or 20.8 cm), in Fig. 6 (smaller cloak height of 0.55 normalized units or 10.4 cm), and in Fig. 7 (yet smaller cloak height of 0.41 normalized units or 7.8 cm). The “ideal” images as depicted in panels (a) of Figs. 5-7 are equivalent to adding a homogeneous dielectric plate with refractive index n = 1/n _min>1 and with the same height as the respective cloak onto the flat mirror. These images are all different and not identical to Fig. 3(a) because the homogeneous dielectric plate with a refractive index larger than that of air effectively reduces the field of view. This effect can be seen by comparing the books at the upper right-hand side corner of Fig. 3(a) and Fig. 5(a). It decreases with decreasing dielectric-plate thickness from Fig. 5(a) via Fig. 6(a) to Fig. 7(a). Furthermore, one can see very faint Fresnel reflections in the high-resolution versions of panels (a) of Figs. 5-7 (see, e.g., necklace on light blue shirt in Fig. 5(a)). The Fresnel transmissions also make the primary images very slightly dimmer than the reference image in Fig. 3(a).

Fig. 5 As Fig. 3, but the refractive-index profile of the grating cloak is divided by n _min (see Fig. 1(a)) such that the minimum refractive index in the cloaking profile becomes n = 1. In (a) and (b), we have added a homogeneous dielectric plate with index n = 1/n _min = 1.37 and with the same height as that of the cloak (see Fig. 1(a)), i.e., 1.1 normalized units (or 20.8 cm). The cloak height is successively reduced in Figs. 6 and 7.

Fig. 6 As Fig. 5, but for a smaller height of the grating-cloak structure of 0.55 normalized units (or 10.4 cm).

The lateral beam displacement mentioned above actually depends strongly on the asymptotic decay of the refractive-index profile to 1/n _min far away from the corrugation as discussed in Ref [15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

]. This effect is immediately seen in the compressed head of the model as in Fig. 4(c) while it is strongly mitigated in Fig. 3(c).

5. Correlation function as quantitative measure of cloaking quality

Early work on electromagnetic cloaking [1

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 (2006). [CrossRef] [PubMed]

–15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

] has been satisfied with some reasonable level of cloaking performance, judged by qualitative assessment. However, as the field of invisibility cloaking is getting more advanced, it is interesting to ask: Is there a quantitative measure for how good the cloaking performance of a particular cloaking device really is? With such a measure, one could quantitatively compare different cloaking approaches. Clearly there is no unique measure of this sort, but (cross-)correlation functions appear attractive to us as they have proven to be a useful approach in other fields of science and technology [30

M. James, Pattern Recognition (John Wiley & Sons, 1988). [PubMed]

,31

J. L. Rodgers and W. A. Nicewander, “Thirteen ways to look at the correlation coefficient,” Am. Stat. 42(1), 59–66 (1988). [CrossRef]

]. Let us briefly recall the underlying mathematics.

Suppose the functions f(x,y)≥0 and g(x,y)≥0 represent the intensity or gray levels of two monochrome two-dimensional finite-size photographs (like, e.g., the images in Figs. 3-7) on an arbitrary scale, where x and y are the orthogonal coordinates in the photographic plane. (More generally, you can consider three such functions for the three principal colors in color photographs.) For example, let f(x,y) be the unperturbed or reference image (i.e., no object and no cloak) and g(x,y) be the image taken with object and cloak. Thus, we will assume that the spatial average value of g is not larger than the spatial average value of f. First, we subtract their respective spatial average values

\begin{array}{l} f (x, y) \to f (x, y) - 〈 f 〉 \\ g (x, y) \to g (x, y) - 〈 g 〉 . \end{array}

(14)

Next, we define a two-dimensional image cross-correlation function, C, according to

C (Δ x, Δ y) : = \frac{\iint f (x, y) g (x + Δ x, y + Δ y) d x d y}{\iint f^{2} (x, y) d x d y} .

(15)

For the case of ideal cloaking, i.e., for f(x,y) = g(x,y) ∀ (x,y), we obtain the autocorrelation C(0,0) = + 100%. The value of C will decay away from the maximum at (Δx = 0,Δy = 0) towards 0% depending on the special properties of the photographic image f(x,y). For example, periodic patterns in f(x,y) will give rise to multiple maxima of C(Δx,Δy). Non-periodic images f(x,y) generally only lead to a single pronounced global maximum of C(Δx,Δy) at (Δx = 0,Δy = 0). If the image g(x,y) = a × f(x,y) with factor 0≤a≤1 is just dimmer than the reference image f(x,y), e.g., due to some absorption or scattering of the cloaking device, we get C(0,0) = a≤100%. If g(x,y) is the “inverse” or complement (i.e., black replaced by white and vice versa) of f(x,y), we get C(0,0) = –100%. If the image g(x,y) is identical to the original f(x,y), but shifted within the photographic xy-plane, a correlation maximum with a value approaching 100% (or strictly equal to + 100% for infinitely large photographs) appears at some shifted position (Δx≠0,Δy≠0). If g(x,y) has nothing to do with f(x,y), the correlation function will exhibit values near 0% for all (Δx, Δy).

Thus, we could either declare C(0,0) or the peak value of C(Δx≠0,Δy≠0) as the requested single-number measure for cloaking quality. The former is stringent, the latter is more generous with respect to possible shifts in the images (like, e.g., for the beam displacement of the carpet cloak). This shift could also be quoted. Note that either of the two measures based on Eq. (15) takes advantage of the entire images and is not based on individual points of the images only.

Ideally, one could finally average over the resulting cross-correlation functions for different reference images f_i (x,y) and corresponding cloaked images g_i (x,y) with index i to eliminate any dependence on the particular image scenery.

Let us now apply the cross-correlation function given in Eq. (15) to the images shown in the previous section. The results are summarized in Table 1 . There, “2D” refers to a near-zero field of view (FOV) in the horizontal direction, i.e., to a vertical line through the center of the images in Figs. 3-7. Note that the grating cloak (like the carpet cloak) is exact in this two-dimensional (2D) world only. Indeed, under these conditions, the correlation is C(0,0) = −20% without cloak and as large as 99% with cloak – essentially perfect cloaking is achieved. Going to a three-dimensional setting reduces the value for the cloaked case to 91% for 10° FOV, to 73% for 25° FOV and to a mere 51% for 50° FOV; whereas, the values for the uncloaked case are all very small or even negative, indicating strongly distorted images.

Table 1 Cross-Correlation Coefficients C: = C(0,0) According to Eq. (15) for Different Combinations of Images, f and g, As Indicated in the Two Left-Hand Side Columns for Different Horizontal Field of Views (FOV)

f	g	FOV = 0°	FOV = 10°	FOV = 25°	FOV = 50°
Figure 3(a)	Figure 3(b)	C = −20%	C = −1%	C = + 7%	C = + 6%
Figure 3(a)	Figure 3(c)	C = + 99%	C = + 91%	C = + 73%	C = + 51%
Figure 4(a)	Figure 4(b)	C = + 17%	C = + 46%	C = + 42%	C = + 43%
Figure 4(a)	Figure 4(c)	C = + 90%	C = + 87%	C = + 73%	C = 61%
Figure 5(a)	Figure 5(b)	C = + 17%	C = + 2%	C = + 8%	C = + 8%
Figure 5(a)	Figure 5(c)	C = + 97%	C = + 91%	C = + 76%	C = + 55%
Figure 3(a)	Figure 5(c)	C = + 74%	C = + 60%	C = + 40%	C = + 28%
Figure 6(a)	Figure 6(b)	C = + 25%	C = + 4%	C = + 10%	C = + 9%
Figure 6(a)	Figure 6(c)	C = + 71%	C = + 81%	C = + 73%	C = + 60%
Figure 3(a)	Figure 6(c)	C = + 75%	C = + 69%	C = + 52%	C = + 42%
Figure 7(a)	Figure 7(b)	C = + 25%	C = + 3%	C = + 10%	C = + 8%
Figure 7(a)	Figure 7(c)	C = + 40%	C = + 64%	C = + 54%	C = + 51%
Figure 3(a)	Figure 7(c)	C = + 40%	C = + 61%	C = + 49%	C = + 41%

FOV = 0 corresponds to a vertical line in the middle of the images or to an effectively 2D situation, whereas increasing FOV corresponds to a more and more three-dimensional viewpoint. The full horizontal (vertical) width in Fig. 3 corresponds to a FOV = 50° (FOV = 42°) full opening angle. The rows in white are reference results for which the correlation should ideally be C = 0%; the rows in light gray are cloaking results for which the correlation should ideally be C = 100%. Cloaking action better than C = 90% is highlighted in yellow.

For the referenced cloaks in Figs. 5-7, one can either compare the cloaked images in Figs. 5(c)-7(c) with the flat mirror in air in Fig. 3(a) or with the flat mirror plus a dielectric plate in Figs. 5(a)-7(a). The former comparison is stricter, the latter is more forgiving. The choice is ambiguous, however, because an observer will likely not be able to make a direct comparison between the two images. In any case, both options are listed in Table 1. Referencing to the flat mirror plus dielectric plate for the thick cloak in Fig. 5 leads to excellent correlations of 97% and 91% for 0° and 10° FOV, respectively. The correlation values are substantially lower when referencing to the flat mirror in air. The correlation values are yet smaller for the thinner cloaks in Figs. 6 and 7. Here, the lateral beam displacement [13

B. Zhang, T. Chan, and B.-I. Wu, “Lateral shift makes a ground-plane cloak detectable,” Phys. Rev. Lett. 104(23), 233903 (2010). [CrossRef] [PubMed]

,15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

], which also depends on the angle of incidence of the rays onto the mirror and hence on the model’s viewing direction, leads to large image distortions even for small vertical FOV and for referencing to the case with dielectric plate.

The Gaussian conformal map shown in Fig. 4 suffers from the same lateral-beam-displacement problem [15

R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]

] – even for large cloak heights – leading to much lower correlation values in Table 1 than for the grating cloak of the same height (Fig. 3) despite the fact that only the middle part along the vertical of the image in Fig. 4(b) is actually distorted due to the single bump (compare Fig. 1). The latter fact can also be seen from the fairly large correlation values for the Gaussian uncloaked case (Fig. 4(b)). If one narrows the cross-correlation analysis down to the middle parts of Fig. 4(a) and 4(c) only, the corresponding C further decreases.

6. Conclusion

We have visualized and quantified the performance of the grating cloak, which is an improved variation of the dielectric carpet cloak. While cloaking does exhibit cross-correlation coefficients near 100% (equivalent to perfect cloaking) for the finite-size grating cloak in 2D (or for small fields of view resembling 2D), the cross-correlation coefficients for the more realistic referenced 3D case with large fields of view are substantially smaller, although cloaking does appear reasonable at first sight. Our results should challenge others to obtain better cloaking as quantified by the cross-correlation coefficient under realistic conditions. We emphasize that the parameter choices in our work correspond to locally isotropic refractive indices between a minimum of 1.0 and a maximum of 2.2. These values are accessible throughout the entire visible spectral region with little dispersion, e.g., by nanostructuring titania along the lines of Ref [11

T. Ergin, N. Stenger, P. Brenner, J. B. Pendry, and M. Wegener, “Three-dimensional invisibility cloak at optical wavelengths,” Science 328(5976), 337–339 (2010). [CrossRef] [PubMed]

Acknowledgements

We thank our model, Tanja Rosentreter (LMU München), and the photographer, Vincent Sprenger (TU München), for help with the photograph used as input for the photorealistic ray-tracing calculations. We thank Tanja Rosentreter (LMU München) and Manuel Decker (KIT Karlsruhe) for help in the editing of the figures. We acknowledge the support of the Arnold Sommerfeld Center (LMU München), which allowed us to use their computer facilities for our rather CPU-time-consuming numerical ray-tracing calculations. J. C. H. gratefully acknowledges financial support by the Excellence Cluster “Nanosystems Inititiative Munich (NIM)”. M. W. acknowledges support by the Deutsche Forschungsgemeinschaft (DFG), the State of Baden-Württemberg, and the Karlsruhe Institute of Technology (KIT) through the DFG Center for Functional Nanostructures (CFN) within subproject A1.5. The project PHOME acknowledges the financial support of the Future and Emerging Technologies (FET) programme within the Seventh Framework Programme for Research of the European Commission, under FET-Open grant number 213390. The project METAMAT is supported by the Bundesministerium für Bildung und Forschung (BMBF).

References and links

1.	J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 (2006). [CrossRef] [PubMed]
2.	U. Leonhardt, “Optical conformal mapping,” Science 312(5781), 1777–1780 (2006). [CrossRef] [PubMed]
3.	V. M. Shalaev, “Physics. Transforming light,” Science 322(5900), 384–386 (2008). [CrossRef] [PubMed]
4.	H. Chen, C. T. Chan, and P. Sheng, “Transformation optics and metamaterials,” Nat. Mater. 9(5), 387–396 (2010). [CrossRef] [PubMed]
5.	M. Wegener and S. Linden, “Shaping optical space with metamaterials,” Phys. Today 63(10), 32–36 (2010). [CrossRef]
6.	J. Li and J. B. Pendry, “Hiding under the carpet: a new strategy for cloaking,” Phys. Rev. Lett. 101(20), 203901 (2008). [CrossRef] [PubMed]
7.	R. Liu, C. Ji, J. J. Mock, J. Y. Chin, T. J. Cui, and D. R. Smith, “Broadband ground-plane cloak,” Science 323(5912), 366–369 (2009). [CrossRef] [PubMed]
8.	J. Valentine, J. Li, T. Zentgraf, G. Bartal, and X. Zhang, “An optical cloak made of dielectrics,” Nat. Mater. 8(7), 568–571 (2009). [CrossRef] [PubMed]
9.	L. H. Gabrielli, J. Cardenas, C. B. Poitras, and M. Lipson, “Silicon nanostructure cloak operating at optical frequencies,” Nat. Photonics 3(8), 461–463 (2009). [CrossRef]
10.	J. H. Lee, J. Blair, V. A. Tamma, Q. Wu, S. J. Rhee, C. J. Summers, and W. Park, “Direct visualization of optical frequency invisibility cloak based on silicon nanorod array,” Opt. Express 17(15), 12922–12928 (2009). [CrossRef] [PubMed]
11.	T. Ergin, N. Stenger, P. Brenner, J. B. Pendry, and M. Wegener, “Three-dimensional invisibility cloak at optical wavelengths,” Science 328(5976), 337–339 (2010). [CrossRef] [PubMed]
12.	H. F. Ma and T. J. Cui, “Three-dimensional broadband ground-plane cloak made of metamaterials,” Nat. Commun. 1(3), 1–6 (2010). [CrossRef]
13.	B. Zhang, T. Chan, and B.-I. Wu, “Lateral shift makes a ground-plane cloak detectable,” Phys. Rev. Lett. 104(23), 233903 (2010). [CrossRef] [PubMed]
14.	J. C. Halimeh, T. Ergin, J. Mueller, N. Stenger, and M. Wegener, “Photorealistic images of carpet cloaks,” Opt. Express 17(22), 19328–19336 (2009). [CrossRef] [PubMed]
15.	R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]
16.	A. S. Glassner, An Introduction to Ray Tracing (Morgan Kaufmann, 1989).
17.	G. Dolling, M. Wegener, S. Linden, and C. Hormann, “Photorealistic images of objects in effective negative-index materials,” Opt. Express 14(5), 1842–1849 (2006). [CrossRef] [PubMed]
18.	A. J. Danner, “Visualizing invisibility: metamaterials-based optical devices in natural environments,” Opt. Express 18(4), 3332–3337 (2010). [CrossRef] [PubMed]
19.	T. Ergin, J. C. Halimeh, N. Stenger, and M. Wegener, “Optical microscopy of 3D carpet cloaks:ray-tracing calculations,” Opt. Express 18(19), 20535–20545 (2010). [CrossRef] [PubMed]
20.	J. L. Synge, Geometrical Mechanics and De Broglie Waves (Cambridge U. Press, 1954).
21.	M. Born, and E. Wolf, Principles of Optics (Pergamon, 1970).
22.	J. S. Desjardins, “Time-dependent geometrical optics,” J. Opt. Soc. Am. 66(10), 1042–1047 (1976). [CrossRef]
23.	J. Molcho and D. Censor, “A simple derivation and an example of Hamiltonian ray propagation,” Am. J. Phys. 54(4), 351–353 (1986). [CrossRef]
24.	P. S. J. Russell and T. A. Birks, “Hamiltonian optics of nonuniform photonic crystals,” J. Lightwave Technol. 17(11), 1982–1988 (1999). [CrossRef]
25.	C. Bellver-Cebreros and M. Rodriguez-Danta, “Eikonal equation from continuum mechanics and analogy between equilibrium of a string and geometrical light rays,” Am. J. Phys. 69(3), 360–367 (2001). [CrossRef]
26.	Y. Jiao, S. Fan, and D. A. B. Miller, “Designing for beam propagation in periodic and nonperiodic photonic nanostructures: extended Hamiltonian method,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 70(3), 036612 (2004). [CrossRef] [PubMed]
27.	D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14(21), 9794–9804 (2006). [CrossRef] [PubMed]
28.	D. A. Genov, S. Zhang, and X. Zhang, “Mimicking celestial mechanics in metamaterials,” Nat. Phys. 5(9), 687–692 (2009). [CrossRef]
29.	K. Niu, C. Song, and M.-L. Ge, “The geodesic form of light-ray trace in the inhomogeneous media,” Opt. Express 17(14), 11753–11767 (2009). [CrossRef] [PubMed]
30.	M. James, Pattern Recognition (John Wiley & Sons, 1988). [PubMed]
31.	J. L. Rodgers and W. A. Nicewander, “Thirteen ways to look at the correlation coefficient,” Am. Stat. 42(1), 59–66 (1988). [CrossRef]

OCIS Codes
(080.0080) Geometric optics : Geometric optics
(080.2710) Geometric optics : Inhomogeneous optical media
(160.3918) Materials : Metamaterials
(230.3205) Optical devices : Invisibility cloaks

ToC Category:
Geometric Optics

History
Original Manuscript: December 23, 2010
Revised Manuscript: February 9, 2011
Manuscript Accepted: February 10, 2011
Published: March 17, 2011

Citation
Jad C. Halimeh, Roman Schmied, and Martin Wegener, "Newtonian photorealistic ray tracing of grating cloaks and correlation-function-based cloaking-quality assessment," Opt. Express 19, 6078-6092 (2011)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-19-7-6078

Sort: Author | Year | Journal | Reset

References

J. B. Pendry, D. Schurig, and D. R. Smith, “Controlling electromagnetic fields,” Science 312(5781), 1780–1782 (2006). [CrossRef] [PubMed]
U. Leonhardt, “Optical conformal mapping,” Science 312(5781), 1777–1780 (2006). [CrossRef] [PubMed]
V. M. Shalaev, “Physics. Transforming light,” Science 322(5900), 384–386 (2008). [CrossRef] [PubMed]
H. Chen, C. T. Chan, and P. Sheng, “Transformation optics and metamaterials,” Nat. Mater. 9(5), 387–396 (2010). [CrossRef] [PubMed]
M. Wegener and S. Linden, “Shaping optical space with metamaterials,” Phys. Today 63(10), 32–36 (2010). [CrossRef]
J. Li and J. B. Pendry, “Hiding under the carpet: a new strategy for cloaking,” Phys. Rev. Lett. 101(20), 203901 (2008). [CrossRef] [PubMed]
R. Liu, C. Ji, J. J. Mock, J. Y. Chin, T. J. Cui, and D. R. Smith, “Broadband ground-plane cloak,” Science 323(5912), 366–369 (2009). [CrossRef] [PubMed]
J. Valentine, J. Li, T. Zentgraf, G. Bartal, and X. Zhang, “An optical cloak made of dielectrics,” Nat. Mater. 8(7), 568–571 (2009). [CrossRef] [PubMed]
L. H. Gabrielli, J. Cardenas, C. B. Poitras, and M. Lipson, “Silicon nanostructure cloak operating at optical frequencies,” Nat. Photonics 3(8), 461–463 (2009). [CrossRef]
J. H. Lee, J. Blair, V. A. Tamma, Q. Wu, S. J. Rhee, C. J. Summers, and W. Park, “Direct visualization of optical frequency invisibility cloak based on silicon nanorod array,” Opt. Express 17(15), 12922–12928 (2009). [CrossRef] [PubMed]
T. Ergin, N. Stenger, P. Brenner, J. B. Pendry, and M. Wegener, “Three-dimensional invisibility cloak at optical wavelengths,” Science 328(5976), 337–339 (2010). [CrossRef] [PubMed]
H. F. Ma and T. J. Cui, “Three-dimensional broadband ground-plane cloak made of metamaterials,” Nat. Commun. 1(3), 1–6 (2010). [CrossRef]
B. Zhang, T. Chan, and B.-I. Wu, “Lateral shift makes a ground-plane cloak detectable,” Phys. Rev. Lett. 104(23), 233903 (2010). [CrossRef] [PubMed]
J. C. Halimeh, T. Ergin, J. Mueller, N. Stenger, and M. Wegener, “Photorealistic images of carpet cloaks,” Opt. Express 17(22), 19328–19336 (2009). [CrossRef] [PubMed]
R. Schmied, J. C. Halimeh, and M. Wegener, “Conformal carpet and grating cloaks,” Opt. Express 18(23), 24361–24367 (2010). [CrossRef] [PubMed]
A. S. Glassner, An Introduction to Ray Tracing (Morgan Kaufmann, 1989).
G. Dolling, M. Wegener, S. Linden, and C. Hormann, “Photorealistic images of objects in effective negative-index materials,” Opt. Express 14(5), 1842–1849 (2006). [CrossRef] [PubMed]
A. J. Danner, “Visualizing invisibility: metamaterials-based optical devices in natural environments,” Opt. Express 18(4), 3332–3337 (2010). [CrossRef] [PubMed]
T. Ergin, J. C. Halimeh, N. Stenger, and M. Wegener, “Optical microscopy of 3D carpet cloaks:ray-tracing calculations,” Opt. Express 18(19), 20535–20545 (2010). [CrossRef] [PubMed]
J. L. Synge, Geometrical Mechanics and De Broglie Waves (Cambridge U. Press, 1954).
M. Born, and E. Wolf, Principles of Optics (Pergamon, 1970).
J. S. Desjardins, “Time-dependent geometrical optics,” J. Opt. Soc. Am. 66(10), 1042–1047 (1976). [CrossRef]
J. Molcho and D. Censor, “A simple derivation and an example of Hamiltonian ray propagation,” Am. J. Phys. 54(4), 351–353 (1986). [CrossRef]
P. S. J. Russell and T. A. Birks, “Hamiltonian optics of nonuniform photonic crystals,” J. Lightwave Technol. 17(11), 1982–1988 (1999). [CrossRef]
C. Bellver-Cebreros and M. Rodriguez-Danta, “Eikonal equation from continuum mechanics and analogy between equilibrium of a string and geometrical light rays,” Am. J. Phys. 69(3), 360–367 (2001). [CrossRef]
Y. Jiao, S. Fan, and D. A. B. Miller, “Designing for beam propagation in periodic and nonperiodic photonic nanostructures: extended Hamiltonian method,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 70(3), 036612 (2004). [CrossRef] [PubMed]
D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14(21), 9794–9804 (2006). [CrossRef] [PubMed]
D. A. Genov, S. Zhang, and X. Zhang, “Mimicking celestial mechanics in metamaterials,” Nat. Phys. 5(9), 687–692 (2009). [CrossRef]
K. Niu, C. Song, and M.-L. Ge, “The geodesic form of light-ray trace in the inhomogeneous media,” Opt. Express 17(14), 11753–11767 (2009). [CrossRef] [PubMed]
M. James, Pattern Recognition (John Wiley & Sons, 1988). [PubMed]
J. L. Rodgers and W. A. Nicewander, “Thirteen ways to look at the correlation coefficient,” Am. Stat. 42(1), 59–66 (1988). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper