1. Introduction

An ideal invisibility cloak makes any macroscopic object located inside of the cloak appear to an outside observer just like empty space – for any viewing direction in three dimensions, for any polarization of the light wave, for inspection in reflection as well as in transmission, and for a broad range of carrier frequencies of light. This ideal represents a demanding benchmark of the ideas of transformation optics [1–4] that has not yet been accomplished experimentally.

However, possible more down-to-earth applications of cloaks may impose less stringent requirements. For example, consider a solar cell. One may want to reduce the shadowing effect of some opaque electric cable or mechanical support structure located above the solar cell to increase the overall energy-conversion efficiency. On a clear day at noon, sunlight approximately impinges from one direction. Thus, the desired effect can be achieved partially by surrounding the cable with a free-space invisibility cloak that works for one direction only. In this case, unlike for the full cloaking problem, image distortions would also be of minor concern – as long as the light is not blocked by the cable and eventually hits the solar cell. As a related but different example, consider an advertising column located in front of a frosted-glass bathroom window, leading to little daylight in the room. A unidirectional free-space cloak around the advertising column would increase the level of daylight in the bathroom. A pedestrian walking by the house facade could still read the advertisement on the column though. Finally, a cylindrical unidirectional free-space cloak slowly rotating around the cylinder axis and containing some fixed good to be advertised might also serve as an appealing eye catcher for shop windows: The good can be seen clearly, slowly turns into strangely distorted, disappears and one can see straight through the arrangement, the good gradually appears again, etc. In this example, minimizing image distortions in the “visible state” would obviously be quite relevant.

It should be mentioned that one can construct rather trivial polarization-independent unidirectional “invisibility cloaks” (see, e.g., [5]). These devices are composed of a set of planar mirrors that guide light rays on a detour around an object to be hidden such that the direction of the emerging light rays is unchanged compared to empty space. These devices work only for one viewing direction because the auxiliary thin mirrors are only “invisible” if the incident rays lie in a plane parallel to these mirrors. However, due to the detour corresponding to the additional geometrical path length in air/vacuum, the optical path length through the apparatus becomes larger than in empty space. Such “ray cloaks” are thus distinct from true “wave cloaks”, which can also properly reconstruct the optical path length, hence reconstruct the time-of-flight (TOF) [6] of a light ray or a light pulse or, equivalently, reconstruct the phase of the wave [7]. To emphasize this aspect, we will show TOF difference images below for the polarization independent cloaks. For the polarization dependent (birefringent) cloaks, the concept of a single TOF is not meaningful as a single incident pulse would generally lead to several emerging main pulses, corresponding to the different paths of ordinary and extraordinary rays.

In two dimensions, unidirectional free-space wave cloaks have already been suggested in the pioneering work of Leonhardt [3,8] in 2006 and in [9,10]. Leonhardt showed by conformal coordinate transformations that a dedicated locally isotropic purely dielectric inhomogeneous material distribution is sufficient for accomplishing this goal. This implies that neither magnetism at optical frequencies nor material anisotropies are mandatory. In retrospect, intuitively speaking, this arrangement can be thought of as composed of two carpet cloaks [11,12] arranged back to back. The original proposal of the carpet cloak by Li and Pendry [11] used a quasi-conformal transformation. Later [12], a strictly conformal version was suggested. Unfortunately, Leonhardt’s particular 2006 suggestion [3,8] still contained singularities. It is critically discussed in [13]. More recently [14], demonstrated experimentally a “perfect” [14] unidirectional free-space cloak composed of piecewise homogeneous anisotropic resonant magnetic materials at microwave frequencies in a two-dimensional waveguide geometry and for one linear polarization of the electromagnetic wave.

In this paper, we aim at visualizing the performance of unidirectional free-space cloaks in three dimensions, i.e., for different viewing directions and for different linear polarizations of light. We consider four different types of designs in our visualizations based on rendering via ray tracing. The first two designs, strictly following [14] in terms of design and parameters, are based on piecewise homogeneous but anisotropic birefringent (magnetic) materials yielding in three dimensions two variations attributed to either using positive- or negative-uniaxial structures for the constituent homogeneous segments (both converge to the same material distribution in the center vertical plane). Design (i) is the positive- and (ii) the negative-uniaxial version. In both, birefringence inherently leads to a dependence on the direction of linear polarization. (iii) The third design is like the first two, but birefringence is eliminated by using anisotropic impedance-matched magneto-dielectric materials. In fact, the transformation that is behind this cloak allows for perfect cloaking in 3D for views along the transformation axis. Designs (i) and (ii) are merely uniaxial approximations of this cloak. (iv) The fourth design, following [3] and [12] in regard to using conformal mapping but not with respect to the particular transformation, uses locally isotropic (hence also polarization independent) but spatially inhomogeneous purely dielectric refractive-index distributions without singularities. Polymer structures according to design (iv) can be fabricated by three-dimensional direct-laser-writing optical lithography [15], possibly even at visible operation frequencies [16]. All designs (i)-(iv) include a metallic coating around the to-be-hidden region. Finally, virtues and limitations of all four approaches shall be discussed.

2. Ray-tracing approach

Ray tracing is an established technique for rendering photorealistic images of sceneries on the level of geometrical optics [2,17,18]. In our previous work [6,19,20], we have used a ray equation of motion that has been derived from Fermat’s principle. We believe that this approach is particularly well suited for treating transformation-optics problems. In [19] we have discussed distributions of locally isotropic refractive indices, in [6] distributions of anisotropic yet impedance-matched magneto-dielectrics, and in [20] spatial distributions of uniaxial birefringent dielectrics. Uniaxial purely magnetic materials are analogous. Interfaces to air are discussed in [20]. This discussion can be generalized to interfaces between two birefringent materials. We refrain from repeating and extending the underlying mathematics here and rather refer the reader to this previous work, noting only that, unlike in [20] where the cloak exhibited dielectric uniaxial anisotropy, here anisotropy is magnetic uniaxial and the ordinary and extraordinary refractive indices are given by $n_o={\sqrt{\in \mu }}_o$ and $n_e={\sqrt{\in \mu }}_e$ respectively, where ${\mu }_o$ and ${\mu }_e$ are the ordinary and extraordinary components of the diagonalized uniaxial permeability tensor and $\in$ is the (constant scalar) permittivity of each constituent segment. Moreover, in this case for extraordinary light, the magnetic-induction vector lies in the extraordinary main section [20], while for ordinary light it has to be perpendicular to the ordinary main section [20].

To allow for direct comparison with our previously rendered images for various other types of cloaks [6,19,20], we use the same scenery as previously. For the case of the above cloak designs (i)-(iii), the setting is illustrated in Fig. 1(a). Likewise, Fig. 1(b) depicts the setting for the cloak design (iv). A virtual point camera looks straight onto a model standing in front of a bookshelf. Curved shelf boards and/or book edges are an intuitive indicator for image distortions. Images shall be depicted for unpolarized detection as well as for detection of horizontally and vertically linear polarized light only.

 

Fig. 1 Illustration showing the virtual camera (blue eye) observing either (a) the piecewise homogeneous anisotropic unidirectional free-space cloak for designs (i)-(iii) or (b) the inhomogeneous isotropic unidirectional free-space cloak for design (iv). The model stands behind the cloak. The cloak in (a) has a square cross section of side length 14.1 cm enveloping a rhombic cloaked region of side length 11.51 cm. In the inset, the local u-axis (v-axis) makes an angle of 33.54 degrees counterclockwise with respect to the y-axis (x-axis). The cloak in (b) has a rectangular cross section of width 50 cm and height 57.14 cm, enveloping a roughly Gaussian cloaked region [12]. Each cloaked region has an area of 114.83 cm² in the center vertical plane. The axis of either cloak is 50 cm away from the model.

Download Full Size | PDF

3. Unidirectional piecewise homogeneous anisotropic square cloaks

In [14], a two-dimensional square cloak composed of a material with anisotropic magnetic permeability tensor and scalar constant electric permittivity is considered. Here, we take exactly the same material parameters and the same square and rhombic shapes of the cloak and corresponding cloaked region, respectively, in the xy-plane (see Fig. 1(a)). The absolute dimensions are quoted in the caption of Fig. 1. However, in three dimensions and for uniaxial constituent materials, one can choose the magnetic uniaxial optic axis to be either along the local u-axis or v-axis leading either to positive-uniaxial (version (i): ${\mu }_o=0.31$, ${\mu }_e=3.27$) or negative-uniaxial (version (ii): ${\mu }_o=3.27$, ${\mu }_e=0.31$) segments.

(i): Results for a configuration like the one shown in Fig. 1(a) are depicted in Fig. 2. By construction in two dimensions [14], for a view along the so-called transformation axis and for horizontal polarization of light, cloaking is perfect in the center vertical plane – where horizontal polarization of light corresponds to the magnetic-induction vector lying in the vertical plane – as shown in the rendered amplitude image in Fig. 2(d). Again by construction, for the same viewing direction but for vertical polarization of light, cloaking does not work at all. Consequently, for unpolarized detection, the image is an admixture of a cloaked and an uncloaked contribution.

 

Fig. 2 The piecewise homogeneous positive-uniaxial birefringent unidirectional free-space cloak under view along the transformation axis. (a) The bare field of view (FOV) as viewed by the virtual camera (see Fig. 1). (b) A rhombic mirror structure is placed in between the model and the virtual camera. (c) The cloak is placed around this structure and viewed for unpolarized light. (d) A horizontal linear polarizer is placed in front of the virtual camera, showing good cloaking behavior as extraordinary light dominates this setting, despite nontrivial ordinary contributions. (e) A vertical linear polarizer is placed in front of the virtual camera, showing bad cloaking behavior as ordinary light is dominant for this setting.

Download Full Size | PDF

General viewing directions towards the left- or right-hand side of the scenery are less obvious and our renderings give a first impression on the cloaking quality. For horizontal polarization of light, the overall cloaking performance is quite good, albeit not perfect. This is due to the fact that this uniaxial approximation of cloak (iii) utilized in versions (i) and (ii) cloaks perfectly only light whose wave vector and corresponding magnetic-induction vector both lie in the vertical plane. That happens only in the center vertical plane, but upon three-dimensional views, this is no longer the case and cloaking is no longer perfect. As to be expected, vertical polarization shows no cloaking at all.

Let us now turn to deviations in direction from the unidirectional design. This means that the cloak shown in Fig. 1(a) is rotated by a certain angle φ with respect to its own axis that is parallel to the z-axis. By design of the carpet cloak [11,12], we expect to see the projection of a flat mirror. This leads to an apparent mirror of width 19.94 cm rotated by φ. We only study the case of horizontal polarization as vertical polarization does not work anyway (see Fig. 2). Corresponding results are depicted in Fig. 3. Again, apart from these to-be-expected imperfections, cloaking is rather good in 3D except for significant Fresnel reflections and nontrivial ordinary-light contributions.

 

Fig. 3 The piecewise homogeneous positive-uniaxial birefringent unidirectional free-space cloak under views away from the transformation axis. The bare mirror structure is shown and the corresponding cloaking result for horizontally linear polarized light, respectively for rotations of (a), (b) 5 degrees, (c), (d) 10 degrees, and (e), (f) 20 degrees around the axis of the cloak. The behavior of the unidirectional cloak for views not along the transformation axis is akin to a carpet cloak. Note the strong Fresnel reflections (b), (d), and (f) and the significant ordinary-light contributions particularly present in (b).

Download Full Size | PDF

(ii): Fig. 4 is like Fig. 2 but for the negative-uniaxial version. Likewise, Fig. 5 is like Fig. 3 but for the negative-uniaxial version. This birefringent version shows significantly better cloaking than design (i) with much weaker Fresnel reflections and ordinary-light contributions.

 

Fig. 4 Same as Fig. 2, but for the negative-uniaxial variant of the piecewise homogeneous birefringent unidirectional free-space cloak. For horizontally polarized light, the cloaking behavior of this cloak is formidable as shown in (d). The step-like behavior of the vertically polarized light is due to ordinary light getting trapped in the corresponding segment of the cloak and escaping only after several total internal reflections. This is due to the fact that for this design ${\mu }_o=3.27\gg 1$ in the corresponding cloak segment. Fresnel-reflection coefficients for this cloak are of the order of less than 3%.

Download Full Size | PDF

 

Fig. 5 Same as Fig. 3, but for the negative-uniaxial variant of the piecewise homogeneous birefringent unidirectional free-space cloak. Note how here the Fresnel reflections are much diminished in contrast to Fig. 3. This cloak acts as a very good three-dimensional polarization dependent carpet cloak for views away from the transformation axis, leading to the illusion of a flat mirror rotated at the same angle the cloak is rotated at.

Download Full Size | PDF

(iii): The undesired polarization dependence of designs (i) and (ii) can be eliminated by considering impedance-matched magneto-dielectric materials for which the electric permittivity tensor equals the magnetic permeability tensor. This cloak is a perfect 3D unidirectional cloak, as shown in Fig. 6. When cloaked, the corresponding mirror structure is completely invisible, and the resulting view in Fig. 6(c) is identical to the bare view (no mirror, no cloak) of Fig. 6(a). Figure 6(d) shows the resulting TOF difference between Fig. 6(c) and Fig. 6(a). The cloaking is incredibly good, with the TOF values being limited only by machine precision. Thus, one can safely conclude that cloak (iii) is a perfect three-dimensional unidirectional free-space ray and wave cloak.

 

Fig. 6 Cloaking behavior of the piecewise homogeneous singly refracting unidirectional free-space cloak for a view along the transformation axis. This is a perfect cloak for this viewing direction. In (b), a rhombic mirror structure is introduced that hides the model and leads to a significant distortion of the bare FOV shown in (a). Upon introducing the singly refracting unidirectional cloak, the original bare FOV is retrieved perfectly, as shown in (c). Wave-cloaking behavior depicted in (d) by relative TOF difference is also perfect and limited only by machine precision, where the values are in the range of yoctoseconds (1 ys = 10⁻²⁴ s). Note the superiority of this design to its uniaxial approximations in cloaks (i) and (ii).

Download Full Size | PDF

Turning to deviations from views along the transformation axis, one sees that cloak (iii) is also a perfect three-dimensional carpet cloak. For any given angle φ of rotation, cloak (iii) produces the illusion of viewing a perfectly flat mirror of width 19.94 cm rotated at the same angle. In fact, simulations (not shown) of such a mirror rotated at an angle φ yields the exact same result shown in Figs. 7(b), 7(d), or 7(f) when φ = 5, 10, or 20 degrees, respectively.

 

Fig. 7 Cloaking behavior of the piecewise homogeneous singly refracting unidirectional free-space cloak for views away from the transformation axis. The bare mirror structure is visualized along with its corresponding cloaked view for rotations of (a), (b) 5 degrees, (c), (d) 10 degrees, and (e), (f) 20 degrees. It is observed that this cloak acts as a perfect carpet cloak for views away from the transformation axis.

Download Full Size | PDF

The anisotropic cloaks (i), (ii), and (iii) indeed exhibit very good to perfect unidirectional free-space cloaking. The disadvantage of design (iii), which is a perfect three-dimensional unidirectional free-space invisibility cloak, is that it requires an anisotropic tensor in both its permeability and permittivity, and hence, magneto-dielectric materials, and this has proven to be quite challenging experimentally. Cloaks (i) and (ii) [14], on the other hand, require only a permeability tensor, and as Figs. 2-5 show, the unidirectional cloaking of each is quite good. In fact, the negative-uniaxial version (ii) shows remarkable similarity in performance to the ideal cloak (iii), but is polarization dependent. Cloak (i) exhibits more impedance mismatch to the sides than cloak (ii). In fact, cloak (ii) shows very weak Fresnel reflections and ordinary-light contributions. Cloak (iii) is perfectly impedance matched to air/vacuum.

It is interesting to note that each of designs (i) and (ii) show significant total internal reflection (TIR) behavior especially for ordinary light as shown in Figs. 2(e) and 4(e). Taking into consideration the first air-cloak interface, TIR occurs in air (Fresnel ordinary reflections) at that interface for design (i) since ${\mu }_{air}=1\gg {\mu }_o=0.31$, which in order to satisfy the continuity relations at the boundary leads to stronger Fresnel extraordinary reflections as seen in Figs. 2(d), 3(b), 3(d), and 3(f), while it occurs within the cloak segment at this interface for design (ii) as ${\mu }_o=3.27\gg {\mu }_{air}=1$ for this design. Moreover, in Figs. 2(d), 3(b), 3(d), and 3(f), one observes a decline in the cloaking performance as compared to Figs. 4(d), 5(b), 5(d), and 5(f). This is due to the fact that upon satisfying the continuity relations at the boundary of the positive-uniaxial cloak, horizontally polarized light has a significant ordinary contributions emerging from the cloak as horizontally polarized light unlike in the case of the negative-uniaxial cloak where such contributions are very weak.

4. Unidirectional inhomogeneous locally isotropic double-Gaussian cloak

In Section 3, piecewise homogeneous anisotropic birefringent (versions (i) and (ii)) and singly refracting (version (iii)) unidirectional free-space cloaks were investigated and their performance scrutinized. Even though version (iii) is an incredible unidirectional cloak that exhibits perfect cloaking in three dimensions, it is experimentally very difficult to realize. On the other hand, even though cloaks (i) and (ii) show remarkable cloaking in three dimensions and are experimentally feasible [14], one still faces the undesirable polarization dependence of the cloaking behavior. As such, we now turn our attention to an experimentally feasible polarization insensitive unidirectional cloak, namely the double-Gaussian unidirectional cloak. This design is based on the Gaussian carpet cloak [12,19] that is inhomogeneous and locally isotropic, and hence polarization independent. In the conformal map according to Eq. (4) of [12], we choose w = 10.02 cm and h = 4.25 cm. The truncating parameters are given in Fig. 1(b). This leads to a minimum (maximum) refractive index of 0.873 (1.920). Polarization independence is a very desirable feature in real-life applications, and particularly in the applications outlined in Section 1 of this manuscript. For example, in the case of solar-cell energy-conversion efficiency, cloaks (i) and (ii) would work for at most half the light impinging from the sun onto the solar panels. On the other hand, the double-Gaussian unidirectional cloak would, as will be shown and discussed shortly, transmit all of the impinging light for views along the transformation axis.

In fact, upon viewing this structure along the transformation axis, one sees very good cloaking in three dimensions. The bare view in Fig. 8(a) is severely distorted upon introducing the metal corrugation in Fig. 8(b). However, upon adding the double-Gaussian unidirectional cloak, the corrugation completely vanishes. Moreover, looking at the left- and right-hand sides of the image, one sees that there is a kind of cutting of the parts of the bare-view image around the middle horizontal line and the stretching of upper and lower parts of the image toward that line. This effect, present on the left- and right-hand sides but not the middle vertical line of the image, diminishes further away from the middle horizontal line of the image. Physically, this is due to the fact that the further from the center vertical plane the viewing direction is, the farther away light emerges from the corrugation on the other side of the cloak, rather than actually striking the corrugation itself. For an application such as solar-cell energy-conversion efficiency, for example, this is ideal, since this small cloaking imperfection in 3D is, for all purposes, perfect behavior for such an application. Also calculated and shown in Fig. 8(d) is the TOF difference between Figs. 8(c) and 8(a). The wave cloaking of this device is surprisingly good indicating a maximum relative error of 6%.

 

Fig. 8 The inhomogeneous locally isotropic double-Gaussian unidirectional free-space cloak for a view along the transformation axis. The Gaussian-like corrugation in (b) is cloaked amazingly well upon introduction of the cloak in (c). Worthy to note is that the cloaking performance in 3D is also very good. Moreover, the wave cloaking of this device is also formidable as shown in (d), where the TOF-difference map indicates a maximum relative error of 6% with values smaller than 0.1 nanosecond (1 ns = 10⁻⁹ s). Fresnel-reflection coefficients for this setting are of the order of 0.01% or less.

Download Full Size | PDF

Figure 9 shows views of the double-Gaussian structure for non-zero rotations around its axis, and hence, these views are not along the transformation axis. To achieve a comparable cloaked volume as its piecewise homogeneous anisotropic counterpart, the double-Gaussian unidirectional cloak exhibits a cloaked region more spread along the transformation axis. This leads the observer to see the illusion of a wider mirror upon rotations than in the case of the cloaks (i)-(iii). However, this can be avoided, in case unfavorable, by changing the values for w and h such that one has a thinner higher bump, and by cutting the edges in a proper manner. Such geometries are not considered here since the main interest lies in unidirectional-cloaking behavior.

 

Fig. 9 The inhomogeneous locally isotropic double-Gaussian unidirectional free-space cloak for views away from the transformation axis. The mirror corrugation and its cloaking are shown for rotations of (a), (b) 5 degrees, (c), (d) 10 degrees, and (e), (f) 20 degrees. It is seen that the double-Gaussian unidirectional cloak behaves similarly to the Gaussian carpet cloak [12,19] for views away from the transformation axis. Fresnel-reflection coefficients for the above settings are 0.2% or less.

Download Full Size | PDF

Even though the conformal map that leads to the carpet cloak [12,19] utilized in the construction of the double-Gaussian unidirectional cloak is over the entire (infinite) half-plane, the actual carpet cloak itself has a finite size. This finite-size truncation leads to imperfect impedance matching. However, our calculations show that Fresnel reflections for all settings shown in Figs. 8 and 9 at the air-cloak interface are very weak, in all less than 0.2%. This is not as good as the ideal design (iii) in terms of impedance matching, but still significantly better (by orders of magnitude) than its birefringent approximations (i) and (ii). Indeed, as illustrated in Figs. 8 and 9, one sees no specular reflections on the surface of the cloak, and for all intents and purposes, for a human observer, this cloak appears to be impedance matched.

5. Conclusion

We have rendered photorealistic images of sceneries in three dimensions including unidirectional free-space invisibility cloaks. Compared to omnidirectional free-space cloaks, the material requirements can be relaxed for the unidirectional case. The birefringent designs (i) and (ii) work well for one particular viewing direction (along the transformation axis), but the performance deteriorates for angles departing from this view. Furthermore, due to the use of birefringence, the devices clearly work for one linear polarization of light only. This aspect is eliminated for the corresponding anisotropic but impedance-matched design (iii), which uses magneto-dielectric materials – at the expense of being much more demanding to realize.

Design (iv) only exploits locally isotropic refractive indices, some of which are below unity. However, as done previously, one can reference this structure to unity via dividing the index distribution by its minimum value [19,21], thus achieving a graded-index minimum refractive index of unity. The unreferenced conformal carpet cloak by design gives the illusion of a flat mirror in vacuum. The referenced conformal cloak will by design give the illusion of a flat mirror in a homogeneous isotropic medium of refractive index 1/n_min, where n_min is the smallest refractive index of the unreferenced conformal carpet cloak. Thus, the referenced unidirectional cloak (not shown) will still achieve similar cloaking quality but will additionally give the illusion that one is looking through a cuboid slab composed of a homogeneous isotropic medium of refractive index n_min.

The free-space cloaking performance of (iv) in terms of range of viewing angles is not as good as (i)-(iii), and this is mainly due to the fact that to achieve a comparable cloaked volume, (iv) requires a cloaked region that is more spread along the transformation axis than in the case of (i)-(iii), thus leading to wider apparent mirrors upon views away from the transformation axis. In terms of performance in three dimensions while viewing along the transformation axis (Figs. 2, 4, 6, and 8), we see that (iii) is a perfect ray and wave cloak. However, its birefringent versions (i) and (ii) show a deterioration in cloaking performance away from the center vertical plane, although design (ii) still exhibits great 3D cloaking behavior for horizontally polarized light. (iv) shows amazingly good performance in three dimensions as well, making it a formidable experimentally feasible three-dimensional unidirectional free-space invisibility cloak. In miniature form, polarization insensitive three-dimensional structures like (iv) could for instance be fabricated in polymer form by direct-laser-writing optical laser lithography [15], combined with stimulated-emission depletion even at visible operation frequencies [16].

Table 1 aims at helping the reader with comparing the photorealistic renderings shown in this paper with all of those of our [6,19–22,24] and other groups’ [23,25,26] previous publications on three-dimensional transformation-optics architectures, including sceneries involving negative-index materials, which can be seen as resulting from a piecewise linear one-dimensional coordinate transformation [8].

Table 1. List of Photorealistic Renderings of Various Devices

View Table