28 January 2019 Transformation optics from macroscopic to nanoscale regimes: a review
Author Affiliations +
Transformation optics is a mathematical method that is based on the geometric interpretation of Maxwell’s equations. This technique enables a direct link between a desired electromagnetic (EM) phenomenon and the material response required for its occurrence, providing a powerful and intuitive design tool for the control of EM fields on all length scales. With the unprecedented design flexibility offered by transformation optics (TO), researchers have demonstrated a host of interesting devices, such as invisibility cloaks, field concentrators, and optical illusion devices. Recently, the applications of TO have been extended to the subwavelength scale to study surface plasmon-assisted phenomena, where a general strategy has been suggested to design and study analytically various plasmonic devices and investigate the associated phenomena, such as nonlocal effects, Casimir interactions, and compact dimensions. We review the basic concept of TO and its advances from macroscopic to the nanoscale regimes.



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.12.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, making them an enabling path for transformation optical structures.


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 ε and μ 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 ε and μ 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).

Fig. 1

A simple coordinate transformation that compresses a space along the x axis. As a result, light follows a distorted trajectory, as shown by the red solid line, but emerges from the compressed region traveling in exactly the same direction with the same phase as before. We can predict the metamaterial properties in the brown region that would realize this trajectory for a light ray. Figure reprinted with permission: Ref. 9, © 2010 by the Imperial College London.


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



Note that ε and μ 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 εyμx=εyμx=1. Hence, we can deduce



To conclude, if we compress a coordinate system along a certain axis, both ε and μ are decreased by the compression factor in the direction of distortion. On the other hand, perpendicular to the direction of distortion, ε and μ 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


where ε¯¯(μ¯¯) and ε¯¯(μ¯¯) are the permittivity (permeability) tensors in the original and transformed space, respectively, and
is the Jacobian matrix of the transformation relating differential elements between the primed and unprimed frames. Therefore, we can rewrite Maxwell’s equations in the transformed frame and derive the physical quantities in terms of Λjj in the transformed space accordingly. In the new coordinate system, we must use modified values of the permittivity and permeability to ensure that Maxwell’s equations are satisfied. Table 1 summarizes the transformation of different physical quantities in terms of the transformation matrix Λjj. Maxwell’s equations thus have the same general form as expected so long as the constitutive parameters and source terms are multiplied by the appropriate factors of the Jacobian matrix and its determinant.

Table 1

Summary of transformations of different physical quantities.

Physical quantitiesBefore transformationAfter transformation
Scalar potentialϕϕ=ϕ
Charge densityρρ=ρ/det(Λ¯¯)
Electric fieldE¯[E¯H¯k¯]=(Λ¯¯T)1·[E¯H¯k¯]
Magnetic fieldH¯
Wave vectork¯
Magnetic flux densityB¯[B¯D¯j¯S¯]=Λ¯¯·[B¯D¯j¯S¯]/det(Λ¯¯)
Electric displacementD¯
Current densityj¯
Poynting vectorS¯

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 R1 by directing the rays around this region while confining the distorted rays within a larger sphere radius R2 [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)].

Fig. 2

(a) The undistorted coordinate system, where a ray of light in free space travels in a straight line. (b) The coordinates are transformed to exclude the cloaked region. Trajectories of rays are pinned to the coordinate mesh and therefore avoid the cloaked region, returning to their original path after passing through the cloak. (c) The coordinates are transformed to fold the space into the annulus region. (d) The field distribution for a cloak under the Gaussian beam illumination. (e) The field distribution for a concentration under the Gaussian beam illumination.


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



This transformation compresses the space in the radial direction. Therefore, along the angular coordinates the values of ε and μ are increased by the inverse of the compression factor, whereas along the radial direction the values are reduced. The corresponding cloak parameters are



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 ε and μ 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.1112.13.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.1617.

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


f(r)={αr,when  r<R1R2αR1R2R1(rR1)+αR1,when  R1<r<R2,
where the constant α stands for the magnification/shrink factor related to the compression/expansion of space. This design gives εt=μt=f(r), εr=μr=ε0f2(r)r2f(r). Note that the cloak can be taken as the extreme case where α0. The incident power is suppressed at the core when 0<α<1 and enhanced when α1. In particular, when α>R2/R1, the negative slope of the function gives a folded region where ε and μ are negative. In this case, the power that flows through the inner region (core) is larger than that which flows through the whole concentrator, due to the energy circulation between the coating and the inner media, as shown in Fig. 2(e). The device enhances the scattering cross section of an object so that it looks like a scatterer larger than the scale of the whole device.2829.30 Space transformation with folded region can also lead to perfect imaging with negative refractive index superlens,3132.33.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.3637.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.


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=αγ, where α is the compression/expansion factor. The intermediate domain is then conformally mapped to the physical domain, giving



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 n2 ranges from 0.68 to 1.96.40 Thus, the anisotropy of the carpet cloak can be neglected. Figure 3(b) shows the Ez 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

Fig. 3

A carpet cloak designed with quasi-conformal mapping. (a) n2 profile of the cloak. (b) Ez distribution with the cloak located within the rectangle in dashed line when a Gaussian beam is launched at 45 deg toward the ground plane. (c) The reflected ray from the bump of the carpet cloak has been shifted toward the incident point. Figure reprinted with permission: (a) and (b) Ref. 40, © 2008 by the American Physical Society (APS); (c) Ref. 41, © 2010 by APS.


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,4243.44 near-infrared,4546.47.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.


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=α1x+β1y+γ, y=α2x+β2y+γ2, z=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


which is homogeneous and anisotropic.51,52

Fig. 4

A linear transformation that transforms an arbitrary triangular region to another one in the physical space.


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)],



Fig. 5

(a) The linear transformation for the design of a carpet cloak. (b) Full-wave simulation of Hz profiles with an incident beam at 45 deg to the homogeneous carpet cloak. Figure reprinted with permission: (b) Ref. 51, © 2009 by IEEE.


In the transformed geometry, the permittivity tensor of the nonmagnetic carpet cloak takes the form of



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 crystals5354.55 or metamaterials with effective uniaxial material profiles, such as dielectric grating structures.5657.58.59 The experimental demonstration of this kind of carpet cloak has been done in visible spectrum, with the calcite crystal53,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 adapters6263.64.65 and magnifying lens.34


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 ϕ(z) in the complex plane with z=x+iy, it must satisfy the Laplace’s equation:



If we make a coordinate transformation z=x+iy=f(z), ϕ(z) also satisfies the Laplace’s equation in the new coordinate system:


provided that the transformation f(z) is analytic everywhere in the region under consideration. Hence, in electrostatics, a conformal transformation preserves the potential in each coordinate system:



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.6667.68.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.


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.

Fig. 6

Schematic of the conformal transformation that maps canonical plasmonic systems to singular structures. (a) A thin metal slab that couples to a 2-D line dipole is transformed to a crescent-shaped nanocylinder illuminated by a uniform electric field. (b) Two semi-infinite metal slabs separated by a thin dielectric film that are excited by a 2-D dipole source are transformed to two touching metallic nanowires illuminated by a uniform electric field.


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:


for the crescent;75


for two touching nanowires.76 Here, ε is the permittivity of the metal, Di and Do denote the inner and outer diameters of the crescent; D1 and D2 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.

Fig. 7

Absorption cross section as a fraction of the physical cross section as a function of frequency for different shapes of (a) crescents and (b) touching nanowires. The absorption spectrum of one individual cylinder is also shown for comparison. (c) The normalized electric field Ex/E0 plotted along the surface of the (a) crescents and (b) touching nanowires at a frequency of 664 THz (where εm=7.02+0.25i). Here the angle θ is defined in the inset.


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,7778. 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.


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


which indicates that the positions of the plasmonic resonances and hence the optical responses of the structure can be controlled by the three geometrical parameters l1, l2, and θ. Properly altering these parameters can tune the frequency and the linewidth of the plasmonic resonances to any predesigned values. Figure 9 shows how the quantized plasmonic modes merge into a continuum when the crescent cusp becomes sharper and sharper. The ability to engineer the physical parameters of metallic nanostructures and thus to maneuver their optical responses provides great promise for the development of new devices and applications. Furthermore, theoretical models can predict unique plasmonic responses of structures that are yet to be fabricated, providing important guidelines for the experiments8586.87 and motivating further synthetic work in these areas.

Fig. 8

The original plasmonic systems are truncated periodic metallo-dielectric structures depicted in (a), (c), (e), and (g), where the EM source is an array of line dipoles (red arrows), aligned along the y axis, with a pitch of 2π. Under the transformation, the structures shown in (a), (c), (e), and (g) are mapped to a (b) metal groove, (d) metal wedge, (f) pair of overlapping nanowires, and (h) a crescent-shaped nanocylinder, respectively. The source of the line dipole array is transformed into a uniform electric field.


Fig. 9

Absorption spectrum as a function of the frequency and the bluntness. Figure reprinted with permission: Ref. 83, © 2012 by APS.



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 zz only modifies the longitudinal permittivity according to εL(k,ω)=ε[1ωp2ω(ω+iγ)β(z)2|k|2], yielding a spatially dependent nonlocal parameter β(z)=dzdzβ (here, β is taken as the hydrodynamic parameter, which is proportional to the surface charge thickness according to δ=β/ω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

Fig. 10

Conformal transformation of a metal-vacuum-metal geometry into a nanowire dimer with (a) local82 and (b) nonlocal88 treatments. Absorption cross-section of 10-nm radius Ag nanowire dimers versus frequency and gap size in (c) local and (d) nonlocal cases.


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, Δd, and permittivity, ε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.9192.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 β=2πc×103 (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.

Fig. 11

(a) Electric field enhancements in the vicinity of the touching point at different degrees of nonlocality β=2πc×104 (blue line), β=2πc×103 (green line), and β=2πc×102 (red line). Solid and dot lines correspond to the analytical and numerical calculations. Gray line plots the field enhancement under local description. (b) Electric field enhancement at the touching point of two nanowires versus frequency and nanowire radius. Figure reprinted with permission: Ref. 89, © 2012 by APS.



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 structures70,75,7778.,88,89,9495.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 ε(r¯)=|dr¯dr¯|ε under an arbitrary 3-D conformal mapping r¯r¯.71,99 Take the 3-D inversion r¯=RT2(r¯R¯0)|r¯R¯0|2 in Fig. 12 as the example. The transformed annulus geometry has an inhomogeneous permittivity given by


where ε is the permittivity of the metal sphere or background in the physical space. Despite this minor complication, the simplification of geometry still enables us to find a closed-form solution of the electrical potential in the inhomogeneous system,100


where Ylm(θ,φ) represents spherical harmonics; alm+ and alm are expansion coefficients, which can be determined through boundary conditions.

Fig. 12

The 3-D inversion that maps a metal-dielectric-metal annulus geometry into a pair of nanospheres separated by a small gap.


The factor |r¯R¯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±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 δ 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 δ. On the other hand, the antibonding modes (blue-dashed lines) blueshift as δ 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,102103.104.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.

Fig. 13

Resonance frequencies of (a) plasmonic modes and (b) Casimir energy versus the separation between two 5-nm-radius gold spheres. Figure reprinted with permission: Ref. 101, © 2013 by the National Academy of Sciences of the United States of America.


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,106107.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 local112 and nonlocal101 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.


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=(d/2π)ln[1/(ez+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 kBloch, which characterizes the energy flow between neighboring unit cells, and the hidden wave vector kHidden, 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.

Fig. 14

A conformal transformation compacts (a) the x dimension of a metallodielectric structure into (b) periodic singularities of a plasmonic metasurface. Dispersion of the singular plasmonic metasurface with respect to (c) both the hidden and Bloch wave vectors, (d) the hidden vector alone. Figure reprinted with permission: Ref. 113, © 2017 by the American Association for the Advancement of Science.


To further illustrate the role of the hidden modes, we set kBloch=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 kHidden [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  THz.

Fig. 15

(a) Schematic of singular graphene metasurface with periodical conductivity. (b) and (c) Absorption spectra of two singular graphene metasurfaces showing how the plasmonic resonances merge into a continuum with increasing dissipation losses. Figure reprinted with permission: Ref. 115, © 2018 by the American Chemical Society.




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,116117.118 elastic waves,119120.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.


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.



Ward A. J., Pendry J. B., “Refraction and geometry in Maxwell’s Equations,” J. Mod. Opt., 43 (4), 773 –793 (1996). https://doi.org/10.1080/09500349608232782 JMOPEW 0950-0340 Google Scholar


Pendry J. B., Schurig D., Smith D. R., “Controlling electromagnetic fields,” Science, 312 (5781), 1780 –1782 (2006). https://doi.org/10.1126/science.1125907 SCIEAS 0036-8075 Google Scholar


Kundtz N. B., Smith D. R., Pendry J. B., “Electromagnetic design with transformation optics,” Proc. IEEE, 99 (10), 1622 –1633 (2011). https://doi.org/10.1109/JPROC.2010.2089664 IEEPAD 0018-9219 Google Scholar


Pendry J. B. et al., “Extremely low frequency plasmons in metallic mesostructures,” Phys. Rev. Lett., 76 (25), 4773 –4776 (1996). https://doi.org/10.1103/PhysRevLett.76.4773 PRLTAO 0031-9007 Google Scholar


Pendry J. B. et al., “Magnetism from conductors and enhanced nonlinear phenomena,” IEEE Trans. Microwave Theory Tech., 47 (11), 2075 –2084 (1999). https://doi.org/10.1109/22.798002 IETMAB 0018-9480 Google Scholar


Smith D. R. et al., “Composite medium with simultaneously negative permeability and permittivity,” Phys. Rev. Lett., 84 (18), 4184 –4187 (2000). https://doi.org/10.1103/PhysRevLett.84.4184 PRLTAO 0031-9007 Google Scholar


Smith D. R., Pendry J. B., Wiltshire M. C. K., “Metamaterials and negative refractive index,” Science, 305 (5685), 788 –792 (2004). https://doi.org/10.1126/science.1096796 SCIEAS 0036-8075 Google Scholar


Ramakrishna S. A., “Physics of negative refractive index materials,” Rep. Prog. Phys., 68 (2), 449 –521 (2005). https://doi.org/10.1088/0034-4885/68/2/R06 RPPHAG 0034-4885 Google Scholar


Luo Y., “Transformation optics applied to plasmonics,” Imperial College London, London, (2012). Google Scholar


Leonhardt U., “Optical conformal mapping,” Science, 312 (5781), 1777 –1780 (2006). https://doi.org/10.1126/science.1126493 SCIEAS 0036-8075 Google Scholar


Cummer S. A. et al., “Full-wave simulations of electromagnetic cloaking structures,” Phys. Rev. E, 74 (3), 036621 (2006). https://doi.org/10.1103/PhysRevE.74.036621 Google Scholar


Cai W. S. et al., “Optical cloaking with metamaterials,” Nat. Photonics, 1 (4), 224 –227 (2007). https://doi.org/10.1038/nphoton.2007.28 NPAHBY 1749-4885 Google Scholar


Zhang J. J., Luo Y., Mortensen N. A., “Minimizing the scattering of a nonmagnetic cloak,” Appl. Phys. Lett., 96 (11), 113511 (2010). https://doi.org/10.1063/1.3366726 APPLAB 0003-6951 Google Scholar


Cai W. et al., “Nonmagnetic cloak with minimized scattering,” Appl. Phys. Lett., 91 (11), 111105 (2007). https://doi.org/10.1063/1.2783266 APPLAB 0003-6951 Google Scholar


Schurig D. et al., “Metamaterial electromagnetic cloak at microwave frequencies,” Science, 314 (5801), 977 –980 (2006). https://doi.org/10.1126/science.1133628 SCIEAS 0036-8075 Google Scholar


Tretyakov S. et al., “Broadband electromagnetic cloaking of long cylindrical objects,” Phys. Rev. Lett., 103 (10), 103905 (2009). https://doi.org/10.1103/PhysRevLett.103.103905 PRLTAO 0031-9007 Google Scholar


Semouchkina E. et al., “An infrared invisibility cloak composed of glass,” Appl. Phys. Lett., 96 (23), 233503 (2010). https://doi.org/10.1063/1.3447794 APPLAB 0003-6951 Google Scholar


Xu S. et al., “Experimental demonstration of a free-space cylindrical cloak without superluminal propagation,” Phys. Rev. Lett., 109 (22), 223903 (2012). https://doi.org/10.1103/PhysRevLett.109.223903 PRLTAO 0031-9007 Google Scholar


Li C., Liu X., Li F., “Experimental observation of invisibility to a broadband electromagnetic pulse by a cloak using transformation media based on inductor-capacitor networks,” Phys. Rev. B, 81 (11), 115133 (2010). https://doi.org/10.1103/PhysRevB.81.115133 Google Scholar


Gomory F. et al., “Experimental realization of a magnetic cloak,” Science, 335 (6075), 1466 –1468 (2012). https://doi.org/10.1126/science.1218316 SCIEAS 0036-8075 Google Scholar


Narayana S., Sato Y., “DC magnetic cloak,” Adv. Mater., 24 (1), 71 –74 (2012). https://doi.org/10.1002/adma.201104012 ADVMEW 0935-9648 Google Scholar


Zhu J. F. et al., “Three-dimensional magnetic cloak working from d.c. to 250 kHz,” Nat. Commun., 6 8931 (2015). https://doi.org/10.1038/ncomms9931 NCAOBW 2041-1723 Google Scholar


Luo Y. et al., “Interaction of an electromagnetic wave with a cone-shaped invisibility cloak and polarization rotator,” Phys. Rev. B, 78 (12), 125108 (2008). https://doi.org/10.1103/PhysRevB.78.125108 Google Scholar


Yang F. et al., “Dc electric invisibility cloak,” Phys. Rev. Lett., 109 (5), 053902 (2012). https://doi.org/10.1103/PhysRevLett.109.053902 PRLTAO 0031-9007 Google Scholar


Han T. C. et al., “Manipulating DC currents with bilayer bulk natural materials,” Adv. Mater., 26 (21), 3478 –3483 (2014). https://doi.org/10.1002/adma.201305586 ADVMEW 0935-9648 Google Scholar


Rahm M. et al., “Design of electromagnetic cloaks and concentrators using form-invariant coordinate transformations of Maxwell’s equations,” Photonics Nanostruct., 6 (1), 87 –95 (2008). https://doi.org/10.1016/j.photonics.2007.07.013 Google Scholar


Luo Y. et al., “Design and analytical full-wave validation of the invisibility cloaks, concentrators, and field rotators created with a general class of transformations,” Phys. Rev. B, 77 (12), 125127 (2008). https://doi.org/10.1103/PhysRevB.77.125127 Google Scholar


Yang T. et al., “Superscatterer: enhancement of scattering with complementary media,” Opt. Express, 16 (22), 18545 –18550 (2008). https://doi.org/10.1364/OE.16.018545 OPEXFF 1094-4087 Google Scholar


Wee W. H., Pendry J. B., “Shrinking optical devices,” New J. Phys., 11 (7), 073033 (2009). https://doi.org/10.1088/1367-2630/11/7/073033 NJOPFM 1367-2630 Google Scholar


Luo Y. et al., “High-directivity antenna with small antenna aperture,” Appl. Phys. Lett., 95 (19), 193506 (2009). https://doi.org/10.1063/1.3264085 APPLAB 0003-6951 Google Scholar


Pendry J. B., “Negative refraction makes a perfect lens,” Phys. Rev. Lett., 85 (18), 3966 –3969 (2000). https://doi.org/10.1103/PhysRevLett.85.3966 PRLTAO 0031-9007 Google Scholar


Yan M., Yan W., Qiu M., “Cylindrical superlens by a coordinate transformation,” Phys. Rev. B, 78 (12), 125113 (2008). https://doi.org/10.1103/PhysRevB.78.125113 Google Scholar


Chen H. Y., Chan C. T., Sheng P., “Transformation optics and metamaterials,” Nat. Mater., 9 (5), 387 –396 (2010). https://doi.org/10.1038/nmat2743 NMAACR 1476-1122 Google Scholar


Wee W. H., Ye Y. J., Luo Y., “Towards a practical compact magnifying superlens-a simple simplical design,” J. Opt. UK, 18 (4), 044011 (2016). https://doi.org/10.1088/2040-8978/18/4/044011 Google Scholar


Lai Y. et al., “Illusion optics: the optical transformation of an object into another object,” Phys. Rev. Lett., 102 (25), 253902 (2009). https://doi.org/10.1103/PhysRevLett.102.253902 PRLTAO 0031-9007 Google Scholar


Jiang W. X. et al., “Broadband all-dielectric magnifying lens for far-field high-resolution imaging,” Adv. Mater., 25 (48), 6963 –6968 (2013). https://doi.org/10.1002/adma.v25.48 ADVMEW 0935-9648 Google Scholar


Li C. et al., “Experimental realization of a circuit-based broadband illusion-optics analogue,” Phys. Rev. Lett., 105 (23), 233906 (2010). https://doi.org/10.1103/PhysRevLett.105.233906 PRLTAO 0031-9007 Google Scholar


Jiang W. X., Cui T. J., “Radar illusion via metamaterials,” Phys. Rev. E, 83 (2), 026601 (2011). https://doi.org/10.1103/PhysRevE.83.026601 Google Scholar


Chang Z. et al., “Design method for quasi-isotropic transformation materials based on inverse Laplace’s equation with sliding boundaries,” Opt. Express, 18 (6), 6089 –6096 (2010). https://doi.org/10.1364/OE.18.006089 OPEXFF 1094-4087 Google Scholar


Li J., Pendry J. B., “Hiding under the carpet: a new strategy for cloaking,” Phys. Rev. Lett., 101 (20), 203901 (2008). https://doi.org/10.1103/PhysRevLett.101.203901 PRLTAO 0031-9007 Google Scholar


Zhang B. L., Chan T., Wu B. I., “Lateral shift makes a ground-plane cloak detectable,” Phys. Rev. Lett., 104 (23), 233903 (2010). https://doi.org/10.1103/PhysRevLett.104.233903 PRLTAO 0031-9007 Google Scholar


Liu R. et al., “Broadband ground-plane cloak,” Science, 323 (5912), 366 –369 (2009). https://doi.org/10.1126/science.1166949 SCIEAS 0036-8075 Google Scholar


Ma H. F., Cui T. J., “Three-dimensional broadband ground-plane cloak made of metamaterials,” Nat. Commun., 1 21 (2010). https://doi.org/10.1038/ncomms1023 NCAOBW 2041-1723 Google Scholar


Ma H. F. et al., “Compact-sized and broadband carpet cloak and free-space cloak,” Opt. Express, 17 (22), 19947 –19959 (2009). https://doi.org/10.1364/OE.17.019947 OPEXFF 1094-4087 Google Scholar


Valentine J. et al., “An optical cloak made of dielectrics,” Nat. Mater., 8 (7), 568 –571 (2009). https://doi.org/10.1038/nmat2461 NMAACR 1476-1122 Google Scholar


Gabrielli L. H. et al., “Silicon nanostructure cloak operating at optical frequencies,” Nat. Photonics, 3 (8), 461 –463 (2009). https://doi.org/10.1038/nphoton.2009.117 NPAHBY 1749-4885 Google Scholar


Ergin T. et al., “Three-dimensional invisibility cloak at optical wavelengths,” Science, 328 (5976), 337 –339 (2010). https://doi.org/10.1126/science.1186351 SCIEAS 0036-8075 Google Scholar


Tamma V. A. et al., “Dispersion characteristics of silicon nanorod based carpet cloaks,” Opt. Express, 18 (25), 25746 –25756 (2010). https://doi.org/10.1364/OE.18.025746 OPEXFF 1094-4087 Google Scholar


Gharghi M. et al., “A carpet cloak for visible light,” Nano Lett., 11 (7), 2825 –2828 (2011). https://doi.org/10.1021/nl201189z NALEFD 1530-6984 Google Scholar


Ergin T., Fischer J., Wegener M., “Optical phase cloaking of 700 nm light waves in the far field by a three-dimensional carpet cloak,” Phys. Rev. Lett., 107 (17), 173901 (2011). https://doi.org/10.1103/PhysRevLett.107.173901 PRLTAO 0031-9007 Google Scholar


Luo Y. et al., “A rigorous analysis of plane-transformed invisibility cloaks,” IEEE Trans. Antennas Propag., 57 (12), 3926 –3933 (2009). https://doi.org/10.1109/TAP.2009.2027824 IETPAK 0018-926X Google Scholar


Xi S. et al., “One-directional perfect cloak created with homogeneous material,” IEEE Microwave Wireless Compon. Lett., 19 (3), 131 –133 (2009). https://doi.org/10.1109/LMWC.2009.2013677 Google Scholar


Chen X. Z. et al., “Macroscopic invisibility cloaking of visible light,” Nat. Commun., 2 176 (2011). https://doi.org/10.1038/ncomms1176 NCAOBW 2041-1723 Google Scholar


Zhang B. L. et al., “Macroscopic invisibility cloak for visible light,” Phys. Rev. Lett., 106 (3), 033901 (2011). https://doi.org/10.1103/PhysRevLett.106.033901 PRLTAO 0031-9007 Google Scholar


Liang D. C. et al., “Robust large dimension terahertz cloaking,” Adv. Mater., 24 (7), 916 –921 (2012). https://doi.org/10.1002/adma.201103890 ADVMEW 0935-9648 Google Scholar


Zhang J. et al., “Homogeneous optical cloak constructed with uniform layered structures,” Opt Express, 19 (9), 8625 –8631 (2011). https://doi.org/10.1364/OE.19.008625 Google Scholar


Xu X. F. et al., “Broad band invisibility cloak made of normal dielectric multilayer,” Appl. Phys. Lett., 99 (15), 154104 (2011). https://doi.org/10.1063/1.3648116 APPLAB 0003-6951 Google Scholar


Popa B. I., Zigoneanu L., Cummer S. A., “Experimental acoustic ground cloak in air,” Phys. Rev. Lett., 106 (25), 253901 (2011). https://doi.org/10.1103/PhysRevLett.106.253901 PRLTAO 0031-9007 Google Scholar


Zigoneanu L., Popa B. I., Cummer S. A., “Three-dimensional broadband omnidirectional acoustic ground cloak,” Nat. Mater., 13 (4), 352 –355 (2014). https://doi.org/10.1038/nmat3901 NMAACR 1476-1122 Google Scholar


Landy N., Smith D. R., “A full-parameter unidirectional metamaterial cloak for microwaves,” Nat. Mater., 12 (1), 25 –28 (2013). https://doi.org/10.1038/nmat3476 NMAACR 1476-1122 Google Scholar


Li W. et al., “A near-perfect invisibility cloak constructed with homogeneous materials,” Opt. Express, 17 (26), 23410 –23416 (2009). https://doi.org/10.1364/OE.17.023410 OPEXFF 1094-4087 Google Scholar


Zhang J. J., Luo Y., Mortensen N. A., “Transmission of electromagnetic waves through sub-wavelength channels,” Opt. Express, 18 (4), 3864 –3870 (2010). https://doi.org/10.1364/OE.18.003864 OPEXFF 1094-4087 Google Scholar


Wang Y. K. et al., “Waveguide devices with homogeneous complementary media,” Opt. Lett., 36 (19), 3855 –3857 (2011). https://doi.org/10.1364/OL.36.003855 OPLEDP 0146-9592 Google Scholar


Xu H. Y. et al., “Dielectric waveguide bending adapter with ideal transmission: practical design strategy of area-preserving affine transformation optics,” J. Opt. Soc. Am. B, 29 (6), 1287 –1290 (2012). https://doi.org/10.1364/JOSAB.29.001287 JOBPDE 0740-3224 Google Scholar


Zhang J. J. et al., “Surface plasmon wave adapter designed with transformation optics,” ACS Nano, 5 (6), 4359 –4364 (2011). https://doi.org/10.1021/nn200516r ANCAC3 1936-0851 Google Scholar


Pendry J. B. et al., “Transformation optics and subwavelength control of light,” Science, 337 (6094), 549 –552 (2012). https://doi.org/10.1126/science.1220600 SCIEAS 0036-8075 Google Scholar


Luo Y. et al., “Harvesting light with transformation optics,” Sci. China Inf. Sci., 56 (12), 1 –13 (2013). https://doi.org/10.1007/s11432-013-5031-2 Google Scholar


Zhao R. K., Luo Y., Pendry J. B., “Transformation optics applied to van der Waals interactions,” Sci. Bull., 61 (1), 59 –67 (2016). https://doi.org/10.1007/s11434-015-0958-x Google Scholar


Pendry J. B., Luo Y., Zhao R. K., “Transforming the optical landscape,” Science, 348 (6234), 521 –524 (2015). https://doi.org/10.1126/science.1261244 SCIEAS 0036-8075 Google Scholar


Luo Y., Pendry J. B., Aubry A., “Surface plasmons and singularities,” Nano Lett., 10 (10), 4186 –4191 (2010). https://doi.org/10.1021/nl102498s NALEFD 1530-6984 Google Scholar


Fernandez-Dominguez A. I., Maier S. A., Pendry J. B., “Collection and concentration of light by touching spheres: a transformation optics approach,” Phys. Rev. Lett., 105 (26), 266807 (2010). https://doi.org/10.1103/PhysRevLett.105.266807 PRLTAO 0031-9007 Google Scholar


Zhang J. J., Zayats A., “Multiple fano resonances in single-layer nonconcentric core-shell nanostructures,” Opt. Express, 21 (7), 8426 –8436 (2013). https://doi.org/10.1364/OE.21.008426 OPEXFF 1094-4087 Google Scholar


Mcphedran R. C., Perrins W. T., “Electrostatic and optical resonances of cylinder pairs,” Appl. Phys., 24 (4), 311 –318 (1981). https://doi.org/10.1007/BF00899728 APPYEK 1080-9198 Google Scholar


Mcphedran R. C., Milton G. W., “Transport-properties of touching cylinder pairs and of the square array of touching cylinders,” Proc. R. Soc. London, Ser. A, 411 (1841), 313 –326 (1987). https://doi.org/10.1098/rspa.1987.0069 Google Scholar


Aubry A. et al., “Broadband plasmonic device concentrating the energy at the nanoscale: the crescent-shaped cylinder,” Phys. Rev. B, 82 (12), 125430 (2010). https://doi.org/10.1103/PhysRevB.82.125430 Google Scholar


Lei D. Y. et al., “Broadband nano-focusing of light using kissing nanowires,” New J. Phys., 12 (9), 093030 (2010). https://doi.org/10.1088/1367-2630/12/9/093030 NJOPFM 1367-2630 Google Scholar


Aubry A. et al., “Plasmonic hybridization between nanowires and a metallic surface: a transformation optics approach,” ACS Nano, 5 (4), 3293 –3308 (2011). https://doi.org/10.1021/nn200438e ANCAC3 1936-0851 Google Scholar


Lei D. Y. et al., “Plasmonic interaction between overlapping nanowires,” ACS Nano, 5 (1), 597 –607 (2011). https://doi.org/10.1021/nn102819p ANCAC3 1936-0851 Google Scholar


Luo Y., Aubry A., Pendry J. B., “Electromagnetic contribution to surface-enhanced Raman scattering from rough metal surfaces: a transformation optics approach,” Phys. Rev. B, 83 (15), 155422 (2011). https://doi.org/10.1103/PhysRevB.83.155422 Google Scholar


Aubry A. et al., “Plasmonic light-harvesting devices over the whole visible spectrum,” Nano Lett., 10 (7), 2574 –2579 (2010). https://doi.org/10.1021/nl101235d NALEFD 1530-6984 Google Scholar


Aubry A. et al., “Conformal transformation applied to plasmonics beyond the quasistatic limit,” Phys. Rev. B, 82 (20), 205109 (2010). https://doi.org/10.1103/PhysRevB.82.205109 Google Scholar


Aubry A. et al., “Interaction between plasmonic nanoparticles revisited with transformation optics,” Phys. Rev. Lett., 105 (23), 233901 (2010). https://doi.org/10.1103/PhysRevLett.105.233901 PRLTAO 0031-9007 Google Scholar


Luo Y. et al., “Broadband light harvesting nanostructures robust to edge bluntness,” Phys. Rev. Lett., 108 (2), 023901 (2012). https://doi.org/10.1103/PhysRevLett.108.023901 PRLTAO 0031-9007 Google Scholar


Luo Y. et al., “Transformation-optics description of plasmonic nanostructures containing blunt edges/corners: from symmetric to asymmetric edge rounding,” ACS Nano, 6 (7), 6492 –6506 (2012). https://doi.org/10.1021/nn3022684 ANCAC3 1936-0851 Google Scholar


Lei D. Y. et al., “Revealing plasmonic gap modes in particle-on-film systems using dark-field spectroscopy,” ACS Nano, 6 (2), 1380 –1386 (2012). https://doi.org/10.1021/nn204190e ANCAC3 1936-0851 Google Scholar


Hanham S. M. et al., “Broadband terahertz plasmonic response of touching InSb disks,” Adv. Mater., 24 (35), OP226 –OP230 (2012). https://doi.org/10.1002/adma.201202003 ADVMEW 0935-9648 Google Scholar


Li G. C. et al., “Metal-substrate-mediated plasmon hybridization in a nanoparticle dimer for photoluminescence line-width shrinking and intensity enhancement,” ACS Nano, 11 (3), 3067 –3080 (2017). https://doi.org/10.1021/acsnano.7b00048 ANCAC3 1936-0851 Google Scholar


Fernandez-Dominguez A. I. et al., “Transformation-optics insight into nonlocal effects in separated nanowires,” Phys. Rev. B, 86 (24), 241110 (2012). https://doi.org/10.1103/PhysRevB.86.241110 Google Scholar


Fernandez-Dominguez A. I. et al., “Transformation-optics description of nonlocal effects in plasmonic nanostructures,” Phys. Rev. Lett., 108 (10), 023901 (2012). https://doi.org/10.1103/PhysRevLett.108.106802 PRLTAO 0031-9007 Google Scholar


Luo Y. et al., “Surface plasmons and nonlocality: a simple model,” Phys. Rev. Lett., 111 (9), 093901 (2013). https://doi.org/10.1103/PhysRevLett.111.093901 PRLTAO 0031-9007 Google Scholar


Ciraci C. et al., “Probing the ultimate limits of plasmonic enhancement,” Science, 337 (6098), 1072 –1074 (2012). https://doi.org/10.1126/science.1224823 SCIEAS 0036-8075 Google Scholar


Wiener A. et al., “Electron-energy loss study of nonlocal effects in connected plasmonic nanoprisms,” ACS Nano, 7 (7), 6287 –6296 (2013). https://doi.org/10.1021/nn402323t ANCAC3 1936-0851 Google Scholar


Zhao Y. D. et al., “Effects of surface roughness of Ag thin films on surface-enhanced Raman spectroscopy of graphene: spatial nonlocality and physisorption strain,” Nanoscale, 6 (3), 1311 –1317 (2014). https://doi.org/10.1039/C3NR05303B NANOHL 2040-3364 Google Scholar


Kraft M. et al., “Transformation optics and hidden symmetries,” Phys. Rev. B, 89 (24), 245125 (2014). https://doi.org/10.1103/PhysRevB.89.245125 Google Scholar


Huidobro P. A. et al., “Hidden symmetries in plasmonic gratings,” Phys. Rev. B, 95 (15), 155401 (2017). https://doi.org/10.1103/PhysRevB.95.155401 Google Scholar


Kraft M. et al., “Designing plasmonic gratings with transformation optics,” Phys. Rev. X, 5 (3), 031029 (2015). https://doi.org/10.1103/PhysRevX.5.031029 PRXHAE 2160-3308 Google Scholar


Schuller J. A. et al., “Plasmonics for extreme light concentration and manipulation,” Nat. Mater., 9 (3), 193 –204 (2010). https://doi.org/10.1038/nmat2630 NMAACR 1476-1122 Google Scholar


Xu L., Chen H. Y., “Conformal transformation optics,” Nat. Photonics, 9 (1), 15 –23 (2015). https://doi.org/10.1038/nphoton.2014.307 NPAHBY 1749-4885 Google Scholar


Fernandez-Dominguez A. I. et al., “Theory of three-dimensional nanocrescent light harvesters,” Nano Lett., 12 (11), 5946 –5953 (2012). https://doi.org/10.1021/nl303377g NALEFD 1530-6984 Google Scholar


Pendry J. B. et al., “Capturing photons with transformation optics,” Nat. Phys., 9 (8), 518 –522 (2013). https://doi.org/10.1038/nphys2667 NPAHAX 1745-2473 Google Scholar


Luo Y., Zhao R. K., Pendry J. B., “Van der Waals interactions at the nanoscale: the effects of nonlocality,” Proc. Natl. Acad. Sci. U. S. A., 111 (52), 18422 –18427 (2014). https://doi.org/10.1073/pnas.1420551111 Google Scholar


Prodan E. et al., “A hybridization model for the plasmon response of complex nanostructures,” Science, 302 (5644), 419 –422 (2003). https://doi.org/10.1126/science.1089171 SCIEAS 0036-8075 Google Scholar


Nordlander P. et al., “Plasmon hybridization in nanoparticle dimers,” Nano Lett., 4 (5), 899 –903 (2004). https://doi.org/10.1021/nl049681c NALEFD 1530-6984 Google Scholar


Prodan E., Nordlander P., “Plasmon hybridization in spherical nanoparticles,” J. Chem. Phys., 120 (11), 5444 –5454 (2004). https://doi.org/10.1063/1.1647518 JCPSA6 0021-9606 Google Scholar


Willingham B., Brandl D. W., Nordlander P., “Plasmon hybridization in nanorod dimers,” Appl. Phys. B-Lasers Opt., 93 (1), 209 –216 (2008). https://doi.org/10.1007/s00340-008-3157-5 Google Scholar


Rodriguez A. W., Capasso F., Johnson S. G., “The Casimir effect in microstructured geometries,” Nat. Photonics, 5 (4), 211 –221 (2011). https://doi.org/10.1038/nphoton.2011.39 NPAHBY 1749-4885 Google Scholar


Emig T. et al., “Casimir forces between arbitrary compact objects,” Phys. Rev. Lett., 99 (17), 170403 (2007). https://doi.org/10.1103/PhysRevLett.99.170403 PRLTAO 0031-9007 Google Scholar


Bimonte G., Emig T., “Exact results for classical casimir interactions: dirichlet and drude model in the sphere-sphere and sphere-plane geometry,” Phys. Rev. Lett., 109 (16), 160403 (2012). https://doi.org/10.1103/PhysRevLett.109.160403 PRLTAO 0031-9007 Google Scholar


Kim K. et al., “Radiative heat transfer in the extreme near field,” Nature, 528 (7582), 387 –391 (2015). https://doi.org/10.1038/nature16070 Google Scholar


Song B. et al., “Enhancement of near-field radiative heat transfer using polar dielectric thin films,” Nat. Nanotechnol., 10 (3), 253 –258 (2015). https://doi.org/10.1038/nnano.2015.6 NNAABX 1748-3387 Google Scholar


Pendry J. B., “Shearing the vacuum—quantum friction,” J. Phys. Condens. Matter, 9 (47), 10301 –10320 (1997). https://doi.org/10.1088/0953-8984/9/47/001 JCOMEL 0953-8984 Google Scholar


Zhao R. et al., “Description of van der Waals interactions using transformation optics,” Phys. Rev. Lett., 111 (3), 033602 (2013). https://doi.org/10.1103/PhysRevLett.111.033602 PRLTAO 0031-9007 Google Scholar


Pendry J. B. et al., “Compacted dimensions and singular plasmonic surfaces,” Science, 358 (6365), 915 –917 (2017). https://doi.org/10.1126/science.aap7939 SCIEAS 0036-8075 Google Scholar


Yang F., Huidobro P. A., Pendry J. B., “Transformation optics approach to singular metasurfaces,” Phys. Rev. B, 98 (12), 125409 (2018). https://doi.org/10.1103/PhysRevB.98.125409 Google Scholar


Galiffi E., Pendry J. B., Huidobro P. A., “Broadband tunable THz absorption with singular graphene metasurfaces,” ACS Nano, 12 (18), 1006 –1013 (2018). https://doi.org/10.1021/acsnano.7b07951 ANCAC3 1936-0851 Google Scholar


Chen H. Y., Chan C. T., “Acoustic cloaking in three dimensions using acoustic metamaterials,” Appl. Phys. Lett., 91 183518 (2007). https://doi.org/10.1063/1.2803315 APPLAB 0003-6951 Google Scholar


Cummer S. A. et al., “Scattering theory derivation of a 3D acoustic cloaking shell,” Phys. Rev. Lett., 100 (2), 024301 (2008). https://doi.org/10.1103/PhysRevLett.100.024301 PRLTAO 0031-9007 Google Scholar


Zhang S., Xia C. G., Fang N., “Broadband acoustic cloak for ultrasound waves,” Phys. Rev. Lett., 106 (2), 024301 (2011). https://doi.org/10.1103/PhysRevLett.106.024301 PRLTAO 0031-9007 Google Scholar


Farhat M., Guenneau S., Enoch S., “Ultrabroadband elastic cloaking in thin plates,” Phys. Rev. Lett., 103 (2), 024301 (2009). https://doi.org/10.1103/PhysRevLett.103.024301 PRLTAO 0031-9007 Google Scholar


Stenger N., Wilhelm M., Wegener M., “Experiments on elastic cloaking in thin plates,” Phys. Rev. Lett., 108 (1), 014301 (2012). https://doi.org/10.1103/PhysRevLett.108.014301 PRLTAO 0031-9007 Google Scholar


Buckmann T. et al., “Mechanical cloak design by direct lattice transformation,” Proc. Natl. Acad. Sci. U. S. A., 112 (16), 4930 –4934 (2015). https://doi.org/10.1073/pnas.1501240112 Google Scholar


Zhang S. et al., “Cloaking of matter waves,” Phys. Rev. Lett., 100 (12), 123002 (2008). https://doi.org/10.1103/PhysRevLett.100.123002 PRLTAO 0031-9007 Google Scholar


Schittny R. et al., “Experiments on transformation thermodynamics: molding the flow of heat,” Phys. Rev. Lett., 110 (19), 195901 (2013). https://doi.org/10.1103/PhysRevLett.110.195901 PRLTAO 0031-9007 Google Scholar


Han T. C. et al., “Full control and manipulation of heat signatures: cloaking, camouflage and thermal metamaterials,” Adv. Mater., 26 (11), 1731 –1734 (2014). https://doi.org/10.1002/adma.201304448 ADVMEW 0935-9648 Google Scholar


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.

© The Authors. Published by SPIE and CLP under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Jingjing Zhang, John B. Pendry, and Yu Luo "Transformation optics from macroscopic to nanoscale regimes: a review," Advanced Photonics 1(1), 014001 (28 January 2019). https://doi.org/10.1117/1.AP.1.1.014001
Received: 11 November 2018; Accepted: 6 December 2018; Published: 28 January 2019

Back to Top