## 1.

## Introduction

Transformation optics (TO) is an emerging technique for the design of advanced electromagnetic (EM) media. It is based on the concept that Maxwell’s equations can be written in a form-invariant manner under coordinate transformations, such that only the permittivity and permeability tensors are modified.^{1}^{–}^{3} With the coordinate transformation applied to the constitutive parameters, EM waves in one coordinate system can be described as if propagating in a different one. The geometric interpretation of Maxwell’s equations utilized in the TO approach provides a powerful and intuitive design tool for the manipulation of EM fields on all length scales.

In the past, the form invariance of Maxwell’s equations has been exploited as a computational tool to simplify numerical electrodynamic simulations. In 1996, a transformation from Cartesian to cylindrical coordinates was applied to solve for the modes of an optical fiber with circular cross section.^{1} This transformation allowed for efficiently solving a cylindrical geometry using a finite-difference computer code implemented in Cartesian coordinates. The TO technique, however, extends well beyond the domain of computational approaches and has gained a great deal of relevance over the past decade in conjunction with the emerging field of metamaterials. Metamaterials are artificially structured media whose effective material parameters can be engineered to have, in principle, any combination of anisotropic electric and magnetic responses,^{4}^{–}^{8} making them an enabling path for transformation optical structures.

## 2.

## Basic Theory

The physical meaning of coordinate transformation can be given as follows. We start from a Cartesian system with a given set of electric and magnetic fields and their associated Poynting vectors. Next, imagine that the coordinates are continuously distorted into a new system. TO was born of the realization that as the system is distorted it carries with it all the associated fields. Hence, to guide the trajectory of light, only a distortion in the underlying coordinate system is needed, automatically taking with it the light ray. Knowledge of the transformation in turn provides the values of $\epsilon $ and $\mu $ required to steer the light in this way.

To give an intuitive view of the TO scheme, we consider a very simple distortion of space: a section of the $x$ axis is compressed as shown in Fig. 1. To find the values of $\epsilon $ and $\mu $ that would lead to the ray trajectory given by the red line, we now illuminate the compressed region with test EM waves of two different polarizations (see the right panel).

For a ray propagating along the $x$ axis, the condition ${k}^{\prime}\alpha d={k}_{0}d$ must hold so that the phase is preserved at the far side of the compressed region. (Here, ${k}_{0}$ is the free-space wave vector, ${k}^{\prime}$ is the wave vector in the compressed region, $\alpha $ is the compression factor, and $d$ is the thickness of the original uncompressed layer.) This condition requires ${k}^{\prime}={k}_{0}\sqrt{{\epsilon}_{y}^{\prime}{\mu}_{z}^{\prime}}={\alpha}^{-1}{k}_{0}$, and thus ${\epsilon}_{y}^{\prime}{\mu}_{z}^{\prime}={\alpha}^{-2}$. Since Maxwell’s equations are symmetrical for $\epsilon $ and $\mu $, any transformation must apply equally to these two parameters.^{2} Therefore, we have

## Eq. (1)

$${\epsilon}_{y}^{\prime}={\mu}_{y}^{\prime}={\epsilon}_{z}^{\prime}={\mu}_{z}^{\prime}={\alpha}^{-1}.$$Note that $\epsilon $ and $\mu $ appear on the same footing because of the symmetry between electric and magnetic fields in the transformed space. For a ray propagating along the $z$ axis, since there is no compression in this direction, the refractive index must remain unchanged $\sqrt{{\epsilon}_{y}^{\prime}{\mu}_{x}^{\prime}}=\sqrt{{\epsilon}_{y}{\mu}_{x}}=1$. Hence, we can deduce

To conclude, if we compress a coordinate system along a certain axis, both $\epsilon $ and $\mu $ are decreased by the compression factor in the direction of distortion. On the other hand, perpendicular to the direction of distortion, $\epsilon $ and $\mu $ are increased by the inverse of the compression factor. For a general compression, the formula is applied successively along each of the three axes.

It is possible to follow this intuitive approach of compressing and expanding space to the design of much more complex and functional devices. However, leveraging the formal structure of electromagnetism, we can follow a systematic, general approach that allows the consideration of arbitrary transformations. Under a general spatial operation, EM fields are distorted in a way that is exactly equivalent to a transformation of the electric permittivity and magnetic permeability tensors of the form

## Eq. (3)

$${\overline{\overline{\epsilon}}}^{\prime}=\frac{\overline{\overline{\mathrm{\Lambda}}}\xb7\overline{\overline{\epsilon}}\xb7{\overline{\overline{\mathrm{\Lambda}}}}^{\mathrm{T}}}{\mathrm{det}(\overline{\overline{\mathrm{\Lambda}}})},\phantom{\rule[-0.0ex]{2.0em}{0.0ex}}{\overline{\overline{\mu}}}^{\prime}=\frac{\overline{\overline{\mathrm{\Lambda}}}\xb7\overline{\overline{\mu}}\xb7{\overline{\overline{\mathrm{\Lambda}}}}^{\mathrm{T}}}{\mathrm{det}(\overline{\overline{\mathrm{\Lambda}}})},$$## Table 1

Summary of transformations of different physical quantities.

Physical quantities | Before transformation | After transformation |
---|---|---|

Scalar potential | $\varphi $ | ${\varphi}^{\prime}=\varphi $ |

Charge density | $\rho $ | ${\rho}^{\prime}=\rho /\mathrm{det}(\overline{\overline{\mathrm{\Lambda}}})$ |

Electric field | $\overline{E}$ | $\left[\begin{array}{l}{\overline{E}}^{\prime}\\ {\overline{H}}^{\prime}\\ {\overline{k}}^{\prime}\end{array}\right]={({\overline{\overline{\mathrm{\Lambda}}}}^{\mathrm{T}})}^{-1}\xb7\left[\begin{array}{l}\overline{E}\\ \overline{H}\\ \overline{k}\end{array}\right]$ |

Magnetic field | $\overline{H}$ | |

Wave vector | $\overline{k}$ | |

Magnetic flux density | $\overline{B}$ | $\left[\begin{array}{l}{\overline{B}}^{\prime}\\ {\overline{D}}^{\prime}\\ {\overline{j}}^{\prime}\\ {\overline{S}}^{\prime}\end{array}\right]=\overline{\overline{\mathrm{\Lambda}}}\xb7\left[\begin{array}{l}\overline{B}\\ \overline{D}\\ \overline{j}\\ \overline{S}\end{array}\right]/\mathrm{det}(\overline{\overline{\mathrm{\Lambda}}})$ |

Electric displacement | $\overline{D}$ | |

Current density | $\overline{j}$ | |

Poynting vector | $\overline{S}$ |

TO has provided a powerful tool for the design of structures capable of controlling the flow of light. The most well-known example is the EM invisibility cloak proposed in 2006.^{2}^{,}^{10} Figure 2(a) shows a ray of light traveling in free space. Suppose we wish to hide the object lying in a spherical region with radius ${R}_{1}$ by directing the rays around this region while confining the distorted rays within a larger sphere radius ${R}_{2}$ [Fig. 2(b)]. In this way an external observer would be aware neither of the presence of the cloak nor its contents. In other words, any object hidden in the cloak would be invisible to outer observers [see Fig. 2(d)].

The corresponding space distortion can be realized through the following coordinate transformation:

## Eq. (4)

$${r}^{\prime}=\left(\frac{{R}_{2}-{R}_{1}}{{R}_{2}}\right)r+{R}_{1},\phantom{\rule[-0.0ex]{1.0em}{0.0ex}}{\theta}^{\prime}=\theta ,\phantom{\rule[-0.0ex]{1.0em}{0.0ex}}{\phi}^{\prime}=\phi .$$This transformation compresses the space in the radial direction. Therefore, along the angular coordinates the values of $\epsilon $ and $\mu $ are increased by the inverse of the compression factor, whereas along the radial direction the values are reduced. The corresponding cloak parameters are

## Eq. (5)

$${\epsilon}_{r}^{\prime}={\mu}_{r}^{\prime}=\frac{{R}_{2}}{{R}_{2}-{R}_{1}}{\left(\frac{r-{R}_{1}}{r}\right)}^{2},\phantom{\rule[-0.0ex]{2.0em}{0.0ex}}{\epsilon}_{\theta}^{\prime}={\mu}_{\theta}^{\prime}={\epsilon}_{\phi}^{\prime}={\mu}_{\phi}^{\prime}=\frac{{R}_{2}}{{R}_{2}-{R}_{1}}.$$If we take a closer look at the cloak parameters, we find that a spherical cloak is singular on its inner surface, i.e., the values of $\epsilon $ and $\mu $ either become infinity or zero. To implement such a cloak, it requires that the metamaterials are in general extremely anisotropic and the difficulty of manufacturing such materials increases severely at higher frequencies. Therefore, any practical implementation will always involve a degree of approximation. A two-dimensional (2-D) cylindrical cloak was later proposed as a simplified scheme, and reduced parameters are also suggested to further simplify the design and realization process.^{11}^{–}^{14} The cylindrical cloak was first demonstrated at microwave frequencies in 2006 using metamaterials,^{15} followed by a huge surge of attentions from different research groups implementing their own versions of the technology.^{16}^{–}^{25}

Another example of TO device is the electromagnetic field concentrator that can focus the incident electromagnetic waves to the central area and enhance the electromagnetic energy density.^{26}^{,}^{27} It can be designed by applying the transformation function

## Eq. (6)

$$f(r)=\{\begin{array}{ll}\alpha r,& \text{when}\text{\hspace{0.17em}\hspace{0.17em}}r<{R}_{1}\\ \frac{{R}_{2}-\alpha {R}_{1}}{{R}_{2}-{R}_{1}}(r-{R}_{1})+\alpha {R}_{1},& \text{when}\text{\hspace{0.17em}\hspace{0.17em}}{R}_{1}<r<{R}_{2}\end{array},$$^{28}

^{–}

^{30}Space transformation with folded region can also lead to perfect imaging with negative refractive index superlens,

^{31}

^{–}

^{34}and optical illusion devices that optically cancel the scattering from original objects and restore the optical paths of light to generate illusions.

^{35}The realization of such devices requires layered complex metamaterial structures to achieve the spatially gradient and anisotropic material parameters. Thus, a simplification process is usually taken to reduce the difficulty in designs and fabrications.

^{36}

^{–}

^{38}Note that devices with extreme-value parameters, including spherical/cylindrical cloaks, superscatterers, superlens, etc., only work within limited frequency band with inevitable absorptions due to the resonant nature of metamaterials.

## 3.

## Quasi-Conformal Mapping

To mitigate the material parameter constraints, the flexibility of the coordinate transformations has been explored. Attentions have been turned to the so-called quasi-conformal mapping that allows the design of devices with isotropic dielectric materials or materials with very small anisotropy. As a natural extension from the conformal mapping, the quasi-conformal mapping relaxes the severe restrictions on the conformal mapping while remaining orthogonal such that the permittivity and permeability tensors can be easily realized. These mappings allow for transformations between domains with different conformal modules in two steps. To compensate for the mismatched conformal modules, the virtual domain is first mapped to an intermediate domain with the same conformal module as the physical domain. This can be simply achieved using a uniform compression/expansion ${y}^{\prime}=\alpha \gamma $, where $\alpha $ is the compression/expansion factor. The intermediate domain is then conformally mapped to the physical domain, giving

## Eq. (7)

$$\overline{\overline{\epsilon}}=\overline{\overline{\mu}}=\mathrm{Diag}[{\alpha}^{-1},\alpha ,{(\alpha |\mathrm{\Lambda}|)}^{-1}].$$In contrast to conformal mapping, the material tensors produced by quasi-conformal mapping are not equal to each other due to perturbations to the conformal module, as can be seen from Eq. (7). However, these perturbations are generally small and thus the resulting anisotropy can typically be ignored. In general, this technique does not have a closed form analytical solution. However, the quasi-conformal map can be approximated by solving Laplace equation on the coordinates.^{39}

One example of quasi-conformal mapping is the design of “carpet cloak,” which provides an alternative form of invisibility cloak.^{40} Instead of completely hiding objects in a free space, a carpet cloak itself appears as a flat reflecting surface so that the objects and the bump made by the objects underneath the carpet are undetectable by light. The carpet cloak can be designed by transforming a rectangular region in the virtual space to an arbitrary region with a defined curved bottom boundary [see Fig. 3(a)]. To minimize the anisotropy factor, the relative size of the cloaked region has to be very small compared to the size of the whole cloak device. For the case in Fig. 3(a), the anisotropy factor becomes a constant of 1.042, and ${n}^{2}$ ranges from 0.68 to 1.96.^{40} Thus, the anisotropy of the carpet cloak can be neglected. Figure 3(b) shows the ${E}_{z}$ distribution when a Gaussian beam is launched at 45 deg toward the ground plane covered by a rectangular carpet cloak with a PEC bottom boundary. The cloak cancels out the original scattering by the bump itself so that it appears as a flat mirror. However, there is an unavoidable lateral shift comparable to the height of the cloaked object in the reflected wave due to the approximations made to the anisotropic parameters, as shown in Fig. 3(c), making the object eventually detectable.^{41}

The material parameters of the carpet cloak are nonsingular and can in principle be realized without the need of resonant features. Furthermore, the isotropic material profile simplifies the fabrication process, making its realization in optical spectrum possible. A number of experimental realizations have been reported, including 2-D and 3-D carpet cloaks working at microwave,^{42}^{–}^{44} near-infrared,^{45}^{–}^{48} and visible spectra.^{49}^{,}^{50} However, in all these realizations, the ratio of the size of the hidden region to that of the cloak is fairly small (normally $<0.01$), due to the isotropy approximation taken in the quasi-conformal mapping. Furthermore, the spatially gradient material parameters require long design and implementation cycles. Despite these limitations, quasi-conformal mappings still hold great utility in bringing TO designs closer to realization.

## 4.

## Linear Transformation

Linear transformation is another type of coordinate transformation used to simplify the device design. In 2-D cases, an arbitrary linear transformation can be described as ${x}^{\prime}={\alpha}_{1}x+{\beta}_{1}y+\gamma $, ${y}^{\prime}={\alpha}_{2}x+{\beta}_{2}y+{\gamma}_{2}$, ${z}^{\prime}=z$, which transfer a line in the initial coordinate to another one in the physical space. Thus, we can find a unique linear transformation to map an arbitrary triangular region to another one (see Fig. 4). The corresponding permittivity/permeability tensor in the transformed geometry is calculated as

## Eq. (8)

$$\overline{\overline{\epsilon}}=\overline{\overline{\mu}}=\left[\begin{array}{cc}\frac{{\alpha}_{1}^{2}+{\beta}_{1}^{2}}{{({\alpha}_{1}{\beta}_{2}-{\alpha}_{2}{\beta}_{1})}^{2}}& \frac{{\alpha}_{1}{\alpha}_{2}+{\beta}_{1}{\beta}_{2}}{{({\alpha}_{1}{\beta}_{2}-{\alpha}_{2}\beta )}^{2}}\\ \frac{{\alpha}_{1}{\alpha}_{2}+{\beta}_{1}{\beta}_{2}}{{({\alpha}_{1}{\beta}_{2}-{\alpha}_{2}\beta )}^{2}}& \frac{{\alpha}_{1}^{2}+{\beta}_{1}^{2}}{{({\alpha}_{1}{\beta}_{2}-{\alpha}_{2}{\beta}_{1})}^{2}}\end{array}\right],$$^{51}

^{,}

^{52}

Therefore, linear transformation opens the way for the design of a homogeneous carpet cloak by linearly compressing the triangular region along the $y$ axis into a polygonal region [see Fig. 5(a)],

## Eq. (9)

$${x}^{\prime}=x,\phantom{\rule[-0.0ex]{1.0em}{0.0ex}}{y}^{\prime}=\frac{{H}_{2}-{H}_{1}}{{H}_{2}}y+\frac{D-x\text{\hspace{0.17em}}\mathrm{sgn}(x)}{D}{H}_{1},\phantom{\rule[-0.0ex]{1.0em}{0.0ex}}{z}^{\prime}=z.$$In the transformed geometry, the permittivity tensor of the nonmagnetic carpet cloak takes the form of

## Eq. (10)

$${\overline{\overline{\epsilon}}}_{x-y}^{\prime}=\epsilon \left[\begin{array}{cc}{\left(\frac{{H}_{2}}{{H}_{2}-{H}_{1}}\right)}^{2}& -\frac{{H}_{1}{H}_{2}^{2}}{{({H}_{2}-{H}_{1})}^{2}d}\text{\hspace{0.17em}}\mathrm{sgn}(x)\\ -\frac{{H}_{1}{H}_{2}^{2}}{{({H}_{2}-{H}_{1})}^{2}d}\text{\hspace{0.17em}}\mathrm{sgn}(x)& 1+{\left(\frac{{H}_{2}}{{H}_{2}-{H}_{1}}\right)}^{2}{\left(\frac{{H}_{1}}{d}\right)}^{2}\end{array}\right].$$Consider a transverse-magnetic (TM, magnetic field perpendicular to the cloak device) polarized Gaussian beam incident obliquely upon such a carpet cloak on top of a flat surface. The anisotropic cloak layer guides the beam around the bump, making the output beam propagate in exactly the same way as that reflected from a flat surface, as shown in Fig. 5(b).

The permittivity tensor in Eq. (10) can be diagonalized by rotating the optical axis, and may be realized with natural birefringent crystals^{53}^{–}^{55} or metamaterials with effective uniaxial material profiles, such as dielectric grating structures.^{56}^{–}^{59} The experimental demonstration of this kind of carpet cloak has been done in visible spectrum, with the calcite crystal^{53}^{,}^{54} and in near infrared spectrum, with the silicon grating structure.^{56} Note that the use of natural crystals enables the realization of optical cloak in macroscopic scale, whereas the metamaterial solution offers more freedom as we can engineer the geometrical parameters, such as the filling factor of the gratings, to achieve anisotropy not found in natural materials. The linear transformation has been further applied to design a unidirectional cloak,^{60} omnidirectional cloaks,^{61} and other photonic devices, such as waveguide adapters^{62}^{–}^{65} and magnifying lens.^{34}

## 5.

## Conformal Transformation

Although quasi-conformal mapping and linear transformation make the realization of devices more feasible, the inhomogeneous or anisotropic distributions of permeability and permittivity are still challenging to achieve, especially at optical spectrum. Conformal transformation is a scheme that can completely eliminate the requisite anisotropy of the material parameters, thus it is especially useful in the device design at optical frequencies.

A conformal mapping is an analytic transformation that preserves local angles. If we consider an analytic function $\varphi (z)$ in the complex plane with $z=x+iy$, it must satisfy the Laplace’s equation:

## Eq. (11)

$$\frac{{\partial}^{2}\varphi}{\partial {x}^{2}}+\frac{{\partial}^{2}\varphi}{\partial {y}^{2}}=0.$$If we make a coordinate transformation ${z}^{\prime}={x}^{\prime}+i{y}^{\prime}=f(z)$, ${\varphi}^{\prime}({z}^{\prime})$ also satisfies the Laplace’s equation in the new coordinate system:

## Eq. (12)

$$\frac{{\partial}^{2}{\varphi}^{\prime}}{\partial {x}^{\prime 2}}+\frac{{\partial}^{2}{\varphi}^{\prime}}{\partial {y}^{\prime 2}}=0,$$Moreover, the preservation of local angles ensures that the boundary conditions in the transformed space remain unchanged. Thus, the dielectric constant of each material is also conserved:

Since both the electrostatic potential and the material permittivity are preserved under the 2-D conformal mapping, the delicate design of a metamaterial with a spatial variation in its constitutive parameters is no longer necessary. Thus, conformal transformation not only simplifies the fabrication process but also provides an easy route to engineering the plasmonic properties of the transformed nanostructures. In recent years, optical conformal transformation has been exploited extensively to treat subwavelength fields occurring in plasmonic nanosystems, which are difficult to study analytically with traditional theoretical methods, giving a precise design tool.^{66}^{–}^{69} By applying transformations to simple plasmonic structures that are well understood but do not have the properties we desire, a variety of complex plasmonic structures, such as sharp edges,^{70} nearly touching spheres,^{71} and nonconcentric core-shell structures,^{72} can be generated and studied analytically through the optical properties of the original simple structures.

## 5.1.

### Singular Plasmonic Structures

From traditional concepts, it is usually believed that a metallic structure should have a large physical size (as compared to the wavelength) to allow for a broadband light-harvesting process and a nanoparticle of finite size usually sustains localized surface plasmon resonances at discrete, rather than continuous, frequencies. However, there are exceptions to these rules. Some finite nanostructures containing sharp edges (or corners) can behave like infinite plasmonic systems and show a continuous interaction with light over a broad frequency range.^{73}^{,}^{74} The theory of TO enables the understanding of the physical mechanism behind this phenomenon. The general strategy starts with a well understood canonical plasmonic system whose analytical description is possible, and then a conformal transformation is applied to deduce the solution for a much more complex geometry. This strategy can explain through the two examples shown in Fig. 6.

A thin slab of metal can support surface plasmon excitations with a lower bound cutoff at the zero frequency and an upper bound cutoff at the surface plasmon frequency. However, the energy is dispersed to infinity and cannot be collected for an efficient light harvesting process. By applying a 2-D inverse transformation, which converts the infinite metal slab in the original space into a crescent-shaped cylinder depicted in Fig. 6(a), a dipole source in addition to the slab is transformed to an incident plane wave, and the induced surface plasmons are now incoming waves focused onto the cusp of the crescent. Another example is shown in Fig. 6(b), where the inverse transformation bridges a plasmonic system consisting of two semi-infinite metal slabs and a pair of touching metallic nanowires. The sharp geometrical features in these nanostructures act as singularities for surface plasmons, causing them to propagate toward the sharp points, slowing down as they progress, but never reaching these sharp points. Consequently, light energy builds up around the singularities. In other words, the transformed geometries can harvest light from infinity and ideally concentrate energy to a nanoscale region around the cusp or touching point.

Detailed calculations show that absorption cross sections for the two geometries depicted in Fig. 6 take the following forms, respectively:

## Eq. (15)

$$\sigma =4{\pi}^{2}{k}_{0}{\left(\frac{\rho}{1-\rho}\right)}^{2}{D}_{o}^{2}\text{\hspace{0.17em}}\mathrm{Re}\left[\frac{\epsilon}{1-{\epsilon}^{2}}{\left(\frac{\epsilon -1}{\epsilon +1}\right)}^{-\frac{2{D}_{i}}{{D}_{o}-{D}_{i}}}\text{\hspace{0.17em}}\mathrm{ln}\right(\frac{\epsilon -1}{\epsilon +1}\left)\right],$$^{75}

## Eq. (16)

$$\sigma ={\pi}^{2}{k}_{0}\frac{{D}_{1}^{2}{D}_{2}^{2}}{{({D}_{1}+{D}_{2})}^{2}}\text{\hspace{0.17em}}\mathrm{Re}\left\{\right[{\left(\frac{\epsilon -1}{\epsilon +1}\right)}^{\frac{{D}_{1}-{D}_{2}}{{D}_{1}+{D}_{2}}}+{\left(\frac{\epsilon -1}{\epsilon +1}\right)}^{\frac{{D}_{2}-{D}_{1}}{{D}_{1}+{D}_{2}}}+2\left]\mathrm{ln}\right(\frac{\epsilon -1}{\epsilon +1}\left)\right\},$$^{76}Here, $\epsilon $ is the permittivity of the metal, ${D}_{i}$ and ${D}_{o}$ denote the inner and outer diameters of the crescent; ${D}_{1}$ and ${D}_{2}$ correspond to the diameter of the two nanowires.

As shown in Figs. 7(a) and 7(b), which display the absorption cross sections as a function of the physical cross section of the crescent and two touching cylinders, both structures exhibit broadband spectra where the absorption cross section is of the order of the overall physical cross section even for such small particle sizes. This broadband light harvesting behavior results from the fact that the induced surface plasmons propagating along the nanostructure surface are compressed toward the singularity, where the group velocity vanishes and energy accumulates. As a consequence, the wavelength is reduced while the energy is increased by the compression factor, leading to a huge increase in the electric field close to the singularity. As illustrated in Figs. 7(c) and 7(d), the maximum field enhancements induced in the crescent and touching-cylinder configurations reach about 1500 and 15,000, respectively.

Based on this presented approach, an analytical relationship between a canonical metallodielectric system and a variety of singular plasmonic structures, including wedges, crescents, rough surfaces, touching cylinders, etc., has been established.^{70}^{,}^{75}^{,}^{77}^{–}^{82} Compared with traditional computational methods, this conformal transformation approach does not require the implementation of adaptive meshes around the sharp geometrical boundaries (such as edges or surface protrusions), thereby allowing for a more comprehensive understanding of the strongly localized surface plasmon modes at the singularities. This transformation approach can be generalized to investigate some other complicated plasmonic systems containing singularities and, therefore, may lead to a large number of practical consequences, such as single molecular detection, new generations of sensors, efficient light harvesting, and surface enhanced Raman scattering.

## 5.2.

### Plasmonic Nanostructures with Blunt Edges/Corners

In real-world applications, singularities or perfectly sharp boundaries in those structures are unlikely to be realized due to limitations in fabrication techniques and the surface tension of the metal. Therefore, the possibility of quantitatively examining how the edge rounding at the sharp boundary will alter the optical responses has great significance on both theoretical and practical levels. TO enables a systematic investigation of a general class of blunt nanostructures by applying conformal mappings to the truncated metallodielectric system associated with the singular structures.^{83}^{,}^{84} A number of such examples are displayed in Fig. 8. The analytical studies reveal how the edge rounding of the nanostructures can be engineered to achieve peculiar optical responses including broadband or selective spectral dependence of absorbed and scattered light as well as large local electric field enhancements. Take the blunt crescent geometry shown in Fig. 8(h) as the example. The resonance condition of the localized surface plasmon modes can be deduced from the corresponding slab geometry in Fig. 8(g),^{83}

## Eq. (17)

$${\left(\frac{\epsilon -1}{\epsilon +1}\right)}^{2}{\left\{\mathrm{exp}\right[\frac{n\pi (2\pi -\theta )}{{l}_{1}+{l}_{2}}]-\mathrm{exp}(\frac{n\pi \theta}{{l}_{1}+{l}_{2}}\left)\right\}}^{2}-{\left[\mathrm{exp}\right(\frac{2{\pi}^{2}n}{{l}_{1}+{l}_{2}})-1]}^{2}=0,$$^{85}

^{–}

^{87}and motivating further synthetic work in these areas.

## 5.3.

### Nonlocal Effects in Plasmonic Nanostructures

In singular structures containing sharp asperities/corners or nearly touching particles, the nonlocal effect, i.e., spatial dispersion in the metal permittivity, plays a key role in the performance of nanodevices, where the classical macroscopic electromagnetism breaks down. An accurate description of optical properties in the subnanometer regime requires the implementation of spatially dispersive permittivities beyond the Drude free electron gas, taking into account the effect of electron–electron interactions. Incorporating nonlocal effects into the TO approach requires the transformation of the permittivity tensor with transverse and longitudinal components under the conformal inversion. The mapping $z\to {z}^{\prime}$ only modifies the longitudinal permittivity according to ${\epsilon}_{L}^{\prime}({\mathbf{k}}^{\prime},\omega )={\epsilon}_{\infty}[1-\frac{{\omega}_{p}^{2}}{\omega (\omega +i\gamma )-{\beta}^{\prime}{({z}^{\prime})}^{2}{|{\mathbf{k}}^{\prime}|}^{2}}]$, yielding a spatially dependent nonlocal parameter ${\beta}^{\prime}({z}^{\prime})=\frac{\mathrm{d}{z}^{\prime}}{\mathrm{d}z}\beta $ (here, $\beta $ is taken as the hydrodynamic parameter, which is proportional to the surface charge thickness according to $\delta =\beta /{\omega}_{p}$). Figures 10(a) and 10(b) show the inversion of a pair of nearly touching nanowires in local and nonlocal cases, respectively. The original uniform charge thickness along the nanowires boundaries maps into a periodic nonuniform width, where the corresponding nonlocal problem can be solved by implementing eikonal approximation.^{88}^{,}^{89}

A more general strategy for the inclusion of nonlocal effects in the TO frame can be done through a recently developed simplified model for nonlocality.^{90} This method treats the nonlocality by replacing the spatially dispersive metal with a composite material consisting of a thin virtual dielectric layer located on top of a local metal. The thickness, $\mathrm{\Delta}d$, and permittivity, ${\epsilon}_{NL}$, of this layer are designed so that the transmission and reflection coefficients for all incident wave vectors and at all frequencies are the same as in the spatially dispersive metal. This model has been implemented in the TO description of separated nanowires and fully-analytic closed-form expressions for the absorption cross section and field enhancement have been obtained.^{90} Figures 10(c) and 10(d) plot the absorption spectrum of the nanowire dimer in the local and nonlocal cases. As the separation between the two nanowires decreases, the surface plasmon modes in the local case redshift toward zero frequency without bound [see Fig. 10(c)]. In sharp comparison, the nonlocal surface charge screening blueshifts the plasmonic resonances, especially at small separation, setting an upper bound for the resonance shifting. This nonlocal saturation phenomenon has been experimentally confirmed in a number of different plasmonic configurations.^{91}^{–}^{93}

The nonlocal effects also set an ultimate bound for the maximum field enhancements. As shown in Fig. 11(a), the electric field along the surface of touching nanowires dramatically decreases with increased nonlocal screening,^{89} and the enhancement in the vicinity of the touching point for $\beta =2\pi c\times {10}^{-3}$ (corresponding to the realistic nonlocality in silver) is reduced by about 50 times. To maximize the field enhancement, the interplay between the nonlocal and radiative effects must be carefully considered. As shown in Fig. 11(b), small particles exhibit poor enhancement owing to strong nonlocal electron screening, whereas, large particles are affected by pronounced radiative losses and hence show reduced enhancement. The balance of these two effects leads to an optimum choice of the radius between 30 and 80 nm.

## 5.4.

### Three-Dimensional Structures and Quantum Fluctuation

TO has been extensively employed in 2-D scenarios to investigate the interaction of light with a number of metal structures^{70}^{,}^{75}^{,}^{77}^{–}^{84}^{,}^{88}^{,}^{89}^{,}^{94}^{–}^{96} processing translation symmetry (see Refs. 66, 97, and 98 for a review). However, those studies describe only a subset of the spectrum with electric fields perpendicular to the symmetric axis. The extension of the TO framework to 3-D is nontrivial because many plasmonic and quantum-fluctuation-related phenomena require a complete description of the plasmonic modes over the whole spectrum. The difficulty of the 3-D generalization lies in the fact that the permittivity of the transformed material is no longer preserved but acquires a spatial dependence according to $\epsilon (\overline{\mathbf{r}})=\left|\frac{\mathrm{d}{\overline{\mathbf{r}}}^{\prime}}{\mathrm{d}\overline{\mathbf{r}}}\right|{\epsilon}^{\prime}$ under an arbitrary 3-D conformal mapping $\overline{\mathbf{r}}\to {\overline{\mathbf{r}}}^{\prime}$.^{71}^{,}^{99} Take the 3-D inversion ${\overline{\mathbf{r}}}^{\prime}=\frac{{R}_{T}^{2}(\overline{\mathbf{r}}-{\overline{\mathbf{R}}}_{0})}{{|\overline{\mathbf{r}}-{\overline{\mathbf{R}}}_{0}|}^{2}}$ in Fig. 12 as the example. The transformed annulus geometry has an inhomogeneous permittivity given by

## Eq. (18)

$$\epsilon (\overline{\mathbf{r}})=\frac{{R}_{T}^{2}}{{|\overline{\mathbf{r}}-{\overline{\mathbf{R}}}_{0}|}^{2}}{\epsilon}^{\prime},$$^{100}

## Eq. (19)

$$\varphi (r,\theta ,\phi )=\sum _{l=0}^{\infty}\sum _{m=-l}^{l}|\overline{\mathbf{r}}-{\overline{\mathbf{R}}}_{\mathbf{0}}|[{a}_{lm}^{+}{(r/{R}_{0})}^{l}+{a}_{lm}^{-}{({R}_{0}/r)}^{l+1}]{Y}_{lm}(\theta ,\phi ),$$The factor $|\overline{\mathbf{r}}-{\overline{\mathbf{R}}}_{\mathbf{0}}|$ in Eq. (19) results from the inhomogeneous permittivity of the annulus geometry. Its presence spoils the spherical symmetry, and as a result, the total angular momentum $l$ is no longer conversed. Nevertheless, detailed calculations show that such a minor complication couples the boundary equations only to the neighboring $l\pm 1$, giving rise to a tridiagonal reflection matrix, which can be solved easily.^{100} The reflection matrix contains all the necessary information of the modes supported by the system. For instance, the resonance frequencies of the modes can be obtained by setting the determinant of the reflection matrix to zero. Results obtained in this way are plotted in Fig. 13(a), which shows that as the separation $\delta $ between the two metal spheres decreases, the surface plasmon modes factor into bonding and antibonding. The bonding modes (red-dashed lines) redshift to zero frequency with decreasing $\delta $. On the other hand, the antibonding modes (blue-dashed lines) blueshift as $\delta $ decreases and finally fall into two branches, one below and the other above the surface plasmons frequency. The branch below the surface plasmon frequency is widely discussed in the literature in the contents of plasmonic hybridization,^{102}^{–}^{105} where the electric fields of these modes are driven out of the gap because of the repulsion of surface charges. In sharp contrast, the branch of modes above the surface plasmon frequency shows a much more interesting behavior. Their mode energy is drawn into the gap despite the repulsion between the charges accumulated on the two particles. Because of their anomalous behavior, these modes are named the anomalous antibonding modes.^{100} They are rarely investigated in previous studies perhaps because they only exist at extremely small gaps, where traditional approaches, such as multipole scattering, demand too much time to get a convergent result.

As mentioned above, the description of the full set of modes enables us to go beyond the optical responses of the structure to study related phenomena in other fields, such as Casimir forces,^{106}^{–}^{108} heat transfer,^{109}^{,}^{110} and quantum friction.^{111} Figure 13(b) displays one such result, where the Casimir energy between two gold spheres is plotted against the gaps between them. Both local^{112} and nonlocal^{101} results are presented, which show that the nonlocal screening effects significantly saturate the Casimir energy and hence decrease the force at extremely small separations. The black-dot line corresponds to the calculation using Drude description of metal permittivity. The large deviations between the Drude mode and realistic metal parameters highlight the significant role of bonded electrons in Casimir interactions.^{101} Apart from the separated metal spheres, the 3-D TO framework described above can be easily extended to other complex plasmonic systems, including 3-D blunt crescents,^{99} oblate and prolate spheroids,^{95} etc.

## 5.5.

### Compact Dimension in Singular Plasmonic Metasurfaces

As the last example of this review article, we show how to use the TO approach to compress a whole spatial dimension of a 3-D system into a set of singular points of a 2-D metasurface.^{113}

In traditional optical systems, the number of characteristic $k$ vectors is normally equal to the dimensionality of the system. To design a plasmonic metasurface with compact dimensions, we start with a 3-D periodically layered metallodielectric stack in Fig. 14(a) and apply a conformal transformation ${z}^{\prime}=(d/2\pi )\mathrm{ln}[1/({e}^{z}+a)-2a]$, converting it into a 2-D singular plasmonic surface in Fig. 14(b). The theory of TO shows that the surface plasmon modes supported by the transformed geometry are characterized by three wave vectors despite the 2-D nature of the metasurface.^{113}^{,}^{114} For instance, the modes propagating along the vertical direction have two characteristic wave vectors [see Fig. 14(c)], i.e., the Bloch wave vector ${k}_{\text{Bloch}}$, which characterizes the energy flow between neighboring unit cells, and the hidden wave vector ${k}_{\text{Hidden}}$, which characterizes how the surface plasmons propagate toward the singular points. The latter is generally missing in traditional gratings or gradient metasurfaces but plays an important role in the singular plasmonic surface considered here.

To further illustrate the role of the hidden modes, we set ${k}_{\text{Bloch}}=0$. Then, for traditional gradient metasurfaces, surface plasmons have to form standing waves resulting in zero power flow everywhere within a unit cell. On the contrary, the singular metasurface is quite different because the hidden modes still have a nonzero ${k}_{\text{Hidden}}$ [see Fig. 14(d)]. These modes can propagate toward the singular points without reflection, as shown by the field plot of Fig. 14(b). In other words, the singular points resulting from the compression of an infinite hidden dimension behave like energy sinks, which completely decouple the power flows within and between unit cells.

Singular points associated with compact dimensions can be achieved by a number of ways. Apart from the geometric singularities depicted in Fig. 14(b), a flat graphene sheet with periodically doped conductivity, the values of which are singular at certain points, can also have a hidden dimension.^{115} Since the continuum of modes associated with a 3-D structure is compressed into a 2-D geometry, the singular metasurface exhibits a number of striking properties that are difficult to achieve with a conventional 2-D plasmonic system. For example, to couple free-space radiation to surface plasmons with traditional methods, the phase matching condition must be carefully considered. In general, at fixed frequencies, the surface modes can be excited only at specific incident angles, whereas at fixed incident angles, the modes can be excited only at discrete frequencies. Nevertheless, the singular metasurfaces no longer have to satisfy such rules. The interplay of Bloch and hidden modes in singular metasurfaces enables broadband wide-angle excitation of surface plasmons. To illustrate this point, Figs. 15(b) and 15(c) plot the absorbance of singular graphene metasurfaces at the normal incidence. The results show that an absorption level around 50% can be achieved over a continuous frequency band of $>10\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{THz}$.

## 6.

## Conclusions

Transformational optics offers unprecedented control over the propagation and confinement of electromagnetic fields at both supra- and subwavelength scales. It not only enables the design of a variety of electromagnetic devices such as the invisibility cloaks and illusion devices but also allows the study of plasmonic nanostructures, revealing the key features that optimize light harvesting and plasmonic field enhancement. Taking into account practical factors including geometric bluntness, nonlocal effects, material limitations, etc., brings the TO approach a step forward toward realistic applications. In fact, the range of applicability of this powerful method is beyond the manipulation of electromagnetic waves and description of the optical properties of metals but can be extended to other physical fields, such as acoustic waves,^{116}^{–}^{118} elastic waves,^{119}^{–}^{121} matter waves,^{122} and thermal fields.^{123}^{,}^{124} The recent breakthrough that shows continuum of modes of the 3-D structure can be directly mapped into the 2-D metasurface, which opens up a brand new direction for the applications of TO. We expect more exciting applications of this methodology and more experimental realizations reported in the future.

## Acknowledgments

This work was supported in part by the Singapore Ministry of Education under Grant Nos. 2017-T1-001-239 (RG91/17 (S)) and MOE2015-T2-1-145, and in part by the Gordon and Betty Moore Foundation.

## References

## Biography

**Jingjing Zhang** joined Nanyang Technological University in 2015 as a senior research fellow. Her research interests include the design and applications of metamaterials, and transformation optics applied to plasmonics.

**John B. Pendry** is a condensed matter theorist and has worked at Imperial since 1981. He has designed a series of completely novel materials with properties not found in nature subsequently the basis for new concepts with radical consequences, such as the first material with a negative refractive index and a prototype cloaking device, which have both caught the imagination of the world’s media, and laid the foundations for the new field of metamaterials.

**Yu Luo** is an assistant professor in the School of Electrical and Electronic Engineering of Nanyang Technological University. He works on a wide range of topics within the realm of metamaterials and plasmonics ranging from the design of invisibility cloaks and plasmonic light-harvesting devices to the study of nonlocal and quantum phenomena in mesoscopic plasmonic systems.