Introduction

During the past ten years, considerable attention has been devoted to the study and application of Transformation Optics (TO)1, 2. Based on geometric transformations, the idea provide a powerful theoretical framework to manipulate ad libitum the propagation of electromagnetic waves. Generally, the applied change of coordinates amounts to replacing the initial materials by equivalent inhomogeneous anisotropic materials3. From a practical point of view, those exotic properties can be mimicked by metamaterials with suitably designed unit cells4. However, most of the TO-based media require conducting or plasmonic inclusions, but the propagation loss represent a major limitation for practical applications, particularly at visible wavelengths.

Certainly the most studied application of TO rely on the idea that an object surrounded by a coating consisting of a material with specially tailored properties becomes invisible to electromagnetic waves5,6,7,8,9,10,11,12,13,14,15,16,17,18. The device that implements this phenomenon has been called an “invisibility cloak”, in reference to Harry Potter, the popular character of J. K. Rowling. But the young wizard has other spells beside his famous cloak. One of the most impressive is the “polyjuice potion” that is able to turn somebody into anybody else’s appearance. This concept as been translated and verified theoretically and experimentally by applying TO to create illusions19,20,21,22,23, i.e. to design a device that gives an arbitrary predefined electromagnetic response to any object placed inside it.

However, as mentioned before, the complexity and limitations of the material properties derived from TO have called for other approaches for practical applications. One route employed to engineering specific wave behaviour is topology optimization24, 25. It has been successfully applied to design devices such as invisibility cloaks26,27,28,29,30,31,32,33, photonic crystals based devices34,35,36, metamaterials37,38,39,40, frequency selective surfaces41 or antennas42,43,44,45, using various strategies such as local gradient-based approaches or global stochastic algorithms.

Recent advances in additive manufacturing, whereby a component is built up layer by layer, open up the design domain significantly and relax fabrication constraints. Techniques such as 3D-printing allows to realize components with intricate complexities, enabling the production of optimal microwave materials and devices with improved performances41, 46,47,48,49.

Herein we report the numerical simulation, optimization, manufacturing and experimental validation of an illusion device at microwave frequencies. The key idea is that we want to make a metallic cylinder (or any object placed inside a perfectly conducting shell with the same outer diameter) behave like a predefined, arbitrary shaped dielectric object (the reference) for a given excitation. A 2D finite element model is employed to simulate the electromagnetic response of the considered structures excited by a Transverse Electric (TE) polarized line source. We use the polymer acrylonitrile butadiene styrene (ABS), a standard thermoplastic used in commercial 3D printers as our low index dielectric material. A gradient based topology optimization algorithm is used to find the layout of a dielectric shell that minimizes the difference between the fields scattered by the illusion and reference structures. The obtained design and reference object were then 3D-printed by a stereolithography based process. The electric field maps of the various configurations are measured in a waveguide configuration and show a good agreement compared too simulated results, demonstrating the illusion effect. We believe this work may serve for designing improved microwave devices such as antennas or sensors where one can control the near and far fields in a predefined manner, while incorporating fabrication and material constraints necessary for real life applications.

Results

Optimized illusion device

As to prove the versatility of the design process, our reference object has an arbitrary complex shape described by Equation(3) and is depicted on Fig. 1(a). It is made of ABS, the same material that we use for the illusion cover, although this is not mandatory nor a constraint of the method. We run the optimization process for various uniform initial densities (see Methods) and obtain different designs. Because of the nature of our optimization algorithm, starting with different material distribution yields different solutions as the problem exhibits many local minima. We retain the design illustrated on Fig. 1(a) as it is a single piece contiguous domain for obvious practical reasons. The resulting 2D objects are extruded by 14 mm in the z direction and then 3D printed (see Methods). A photograph of the samples is shown on Fig. 1(c), proving the ability to built complex structures using additive manufacturing techniques.

Figure 1
figure 1

Optimization results. (a) Reference dielectric object. (b) Optimized illusion device. (c) Photograph of the 3D printed structures.

Near field

As described in Methods, we map experimentally the electric field in a waveguide to illustrate the illusion effect. Our measurement results are reported on Fig. 2(e–h) and compared with simulations (a–d). We first consider the empty waveguide and obtain a cylindrical propagating wave with excellent agreement with simulations (the free space 2D Green’s function). This indicates that the line source we used in simulations is a fairly good approximation of the waveguide horn of our experiments. We then map the field in the case of the PEC cylinder, which shows a strong scattering and a shadow region in the forward direction as expected from theory. Due to its asymmetrical shape, the reference dielectric object diffracts the incident wave in a characteristic pattern. The measured values agree really well with simulations results and prove our ability to map complex diffracted fields with sub-wavelength resolution. Finally, the illusion device wrapped around the metallic cylinder shows a scattering in free space very similar to the reference object. The measured electric field matches the reference field as expected from the numerical optimization, more particularly the shadow caused by the PEC cylinder is reduced and the characteristic asymmetric pattern is retrieved.

Figure 2
figure 2

Field maps. Real part of the dielectric field at the design frequency f = 10 GHz. Top row (ad) FEM simulations. Bottom row (eh) measurements. Various case are considered: free space (a,e), PEC cylinder (b,f), reference dielectric object (c,g) and the same PEC cylinder with the proposed illusion device (d,h).

To illustrate further these results, we plot cuts of those field maps on Fig. 3 at a constant radius r = 1.5R ill. The amplitude normalized to the incident field and phase are displayed for measurements (thick solid lines) and simulations (thin dashed lines) and are in fairly good agreement. As mentioned before, the PEC cylinder provoke a strong attenuation in the forward direction (θ = 0 degree), and backscattering. As for the reference and illusion cases, their amplitude patterns show strong similarities, particularly in the forward direction where we can observe a symmetric angular response in both cases. Reference and illusion phase patterns are in good agreement and mostly differ from the PEC case in the forward direction, albeit some noise in the measured data for the illusion case.

Figure 3
figure 3

Amplitude and phase. Electric field magnitude normalized to the free space value (in dB, left column) and phase (in degree, right column) for the PEC cylinder (a,b), the reference dielectric object (c,d) and the illusion device (e,f). The measured (thick solid lines) and computed (thin dashed lines) values are plotted as a function of the observation angle θ (in degree) for a fixed radius r = 1.5 R ill.

Performances

In order to confirm the previous qualitative observations, we calculate two figures of merit to asses the performances of the illusion device. The first one is the partial objective function ϕ 1 (defined in Equation 5) which is essentially the relative error between the electric fields for the reference and illusion cases. As can be seen on Fig. 4(a), ϕ 1 exhibits a dip in frequency at 10 GHz, showing that the difference is minimal between the reference and illusion fields. While the simulated value is 0.01, the measured one is around 0.6, which is mainly due to experimental errors and to the local nature of this figure of merit.

Figure 4
figure 4

Performances of the illusion device. Measured (thick solid lines) and computed (thin dashed lines) figures of merit as a function of frequency f: partial objective function ϕ 1 representing the relative error between the reference and illusion electric fields (a), and correlation coefficient between those fields for the real (b) and imaginary (c) parts.

Another parameter suitable to quantify the performance of the illusion device is the correlation coefficient50 defined as:

$$\rho (X,Y)=\frac{{\rm{E}}(XY)-{\rm{E}}(X){\rm{E}}(Y)}{\sqrt{{\rm{E}}({X}^{2})-{\rm{E}}{(X)}^{2}}\sqrt{{\rm{E}}({Y}^{2})-{\rm{E}}{(Y)}^{2}}}$$
(1)

where X = E ref is the field produced by the reference dielectric object and Y = E ill is the field scattered by the PEC cylinder covered by the illusion device. We calculated this coefficient for both real and imaginary parts of the fields and plot it as a function of frequency on Fig. 4(b) and (c) respectively. In both cases, it reaches its maximum at the design frequency of 10 GHz, with an experimental value around 0.9, while the theoretical simulated value is 0.99 (a coefficient of 1 means that both fields are perfectly correlated). This again prove that our device reproduces the desired scattering with good fidelity.

Discussion

We designed an illusion device made of a relatively thin (1.5λ) low index dielectric cover with an optimized topology. When wrapped around a given PEC cylinder its function is to emulate the scattering of a predefined, arbitrarily shaped dielectric object so that an observer probing the electromagnetic field outside of the device is tricked into thinking that he detects the said reference object. Note that if the PEC cylinder is replaced by a thin PEC surface, the device will produce the illusion for any object places inside as the field does not penetrate. This device has been fabricated using 3D printing using a standard stereolithography process and common ABS thermoplastic material. The illusion effect is present in the near field as evidences experimental by spatially resolved waveguide measurements. The performances of the device have been assessed by two figures of merit and are in reasonable accordance with the numerical solutions. Although the illusion effect has been proven here experimentally at microwave frequencies, the design methodology is general and suitable for the optical domain. Furthermore, it only requires a low index dielectric and doe not necessitate anisotropic or plasmonic materials. We believe that the optimization route used here might be of interest for the design of other devices and materials with enhanced characteristics able to control the propagation of electromagnetic waves, not only in the far field but in the near field. The development and refinement of additive manufacturing process, including the ability to use different materials, allows the design of topologies difficult to realize with traditional fabrication techniques and claim for innovative engineering approaches. Advanced numerical techniques, efficient global optimization algorithms and innovative strategies need to be developed in order to improve the design process and the designed structures, by including multiple objectives, constraints, multiple materials and physics. It is believed that this could benefits applications such as telecommunications, antennas, sensors and energy harvesting to name a few.

Methods

Numerical simulations: electromagnetic problem and optimization procedure

The problem at stake is governed by Maxwell’s equation in time-harmonic regime (with convention ext (−iωt) discarded hereafter). The material parameters do not depend on z and this medium is excited by a line source located at r′ = (x′, y′)T and directed along the z axis. This is the so called TE propagation with the total electric field E = E z satisfying the following scalar wave equation:

$$\nabla \cdot ({\mu }^{-1}\nabla E)+{k}^{2}\varepsilon E=\delta ({\bf{r}}{\boldsymbol{^{\prime} }}),$$
(2)

where ε and μ are the relative permittivity and permeability respectively, δ is the Dirac delta distribution, r = (x, y)T is the spatial position vector and k = ω/c is the free space wavenumber.

To solve this problem numerically, we use a Finite Element Method51 (FEM) commercial package (COMSOL Multiphysics). A schematics of the problem under study is given in Fig. 5(b). The line source is located at x′ = −260 mm and y′ = 0 mm. The background domain is truncated by cartesian Perfectly Matched Layers52, 53 (PML) in order to damp propagating waves, with thickness t pml = λ 0, constant stretching coefficients along x and y (s x  = s y  = 1 − i), and homogeneous Neumann condition on their outward boundaries.

Figure 5
figure 5

Numerical and experimental methods. (a) Optimization algorithm workflow. (b) Computational setup for the FEM optimization problem. (c) Measurement setup.

For this work, we used ABS thermoplastic commonly used for 3D printing as our dielectric material, with a permittivity ε ABS = 2.6980 − 0.0068i as measured at the operating frequency f 0 = 10 GHz (see Methods). Furthermore, all materials are considered non magnetic (μ = 1).

The aim of the proposed illusion device is to design a dielectric cover that will make a metallic cylinder of radius R cyl = 35.7 mm look like an arbitrarily shaped dielectric. The shape of this reference object is given in polar coordinates (R, θ) as a Fourier series:

$$R(\theta )={R}_{0}+\sum _{n\mathrm{=1}}^{5}{A}_{n}\,\cos (n\theta )+{B}_{n}\,\sin (n\theta ),$$
(3)

with R 0 = 15 mm, A = (−0.83, −4.50, 4.02, 4.44, −0.09) mm and B = (−0.11, −1.62, 4.00, −1.31, −3.89) mm.

The objective is to minimize the difference between the field E ill scattered by the PEC cylinder covered by the illusion layer and the field E ref produced by the reference object in the background domain Ω b. Defining the integral quantities

$$I={\int }_{{{\rm{\Omega }}}_{{\rm{b}}}}{|{E}_{{\rm{ill}}}-{E}_{{\rm{ref}}}|}^{2}{\rm{d}}{\bf{r}},{\rm{and}}\quad {I}_{{\rm{ref}}}={\int }_{{{\rm{\Omega }}}_{{\rm{b}}}}{|{E}_{{\rm{ref}}}|}^{2}{\rm{d}}{\bf{r}},$$
(4)

the objective function to minimize is

$${\varphi }_{1}({\varepsilon }_{{\rm{des}}})=I/{I}_{{\rm{ref}}}$$
(5)

where ε des is the permittivity in the design domain (an annular layer of thickness t des = 1.5 λ 0 = 45 mm), subject to constraints ε min ≤ ε des ≤ ε max, with ε min = ε air = 1 and ε max = ε ABS = 2.6980 − 0.0068i.

In addition, we set additional constraint on the gradient of permittivity to avoid chessboard patterns and very small features by defining a second objective function

$${\varphi }_{2}({\varepsilon }_{{\rm{des}}})=\frac{{h}_{{\rm{mesh}}}^{2}}{S}{|\nabla {\varepsilon }_{{\rm{des}}}|}^{2},$$
(6)

where S = πt des (t des + 2R cyl) is the area of the design domain and h mesh is the maximum size of a mesh element.

Following the so-called solid isotropic material interpolation (SIMP)54, we define a density function ξ such that

$${\varepsilon }_{{\rm{des}}}(\xi )={\varepsilon }_{{\rm{\min }}}+{\xi }^{p}({\varepsilon }_{{\rm{\max }}}-{\varepsilon }_{{\rm{\min }}})$$

where p is a penalty factor (p = 1 here).

In order to enforce a binary design, we apply a threshold projection55 at each global iteration:

$$\xi (s)=\frac{\tan \,{\rm{h}}(\beta \nu )+\,\tan \,{\rm{h}}(\beta (s-\nu ))}{\tan \,{\rm{h}}(\beta \nu )+\,\tan \,{\rm{h}}(\beta \mathrm{(1}-\nu ))},$$

with ν = 1/2 and β = 2N where N < N max is the number of optimization global iteration (we set N max = 8).

The topology optimization problem is thus to minimize the functional ϕ of the final design variable s

$$\varphi (s)=\gamma \,{\varphi }_{1}(s)+\mathrm{(1}-\gamma )\,{\varphi }_{2}(s),$$
(7)

subject to constraints 0 ≤ s ≤ 1 and where γ is a weight parameter between 0 and 1 (γ = 0.75 here). The gradient-based optimization code uses the Quasi-Newton method56 to update the design iteratively, starting from an initial homogeneous density s 0 = 0.7. The sensitivity of the objective function with respect to s is calculated by the adjoint method36. The problem is discretized by the FEM using second order Lagrange elements and the solutions are computed with a direct solver57. An overview of the optimization algorithm is given on Fig. 5(a).

Fabrication and Material Characterisation

The samples were 3D printed using an Objet30 Prime from Stratasys. The printer uses a polyjet technique to build smooth and accurate 3D models. The building process combines the techniques of inkjet and stereolythography (SLA) printing methods. The tiny droplets of liquid photopolymer are cured with a UV laser beam which instantly solidify the photopolymer creating the layers of the building model. Each of these layers has an accuracy down to 100 microns. For this particular experiment, VeroWhitePlus (RGD835) material was used. Three different models of 60 × 60 × 3 mm were manufactured using a high-resolution and a glossy finish. The printing environment was configured so the support material is built around the sample. Once the printing is finished, the gel-like support is easily removed using a high pressure water flow. The prototypes were left to dry for a couple of days thus water effects do not cause any type of uncertainties on the material characterisation.

The dielectric properties were measured using an 85072A Split-Cylinder from Keysight Technologies. The instrument operates at 10 GHz with a high Q, producing good loss tangent resolution thanks to the adjustable electrical coupling into the cavity. The latter resonator is connected to a 5230C PNA-L Network Analyser, from the same supplier, which runs the 85071E Materials Measurement Software. Multiple measurements were taken for each sample with an IF bandwidth of 100 Hz. The average results from 5 independent measurements yield a 2.698 for the real part of the complex permittivity and 0.0068 for the imaginary part.

Experimental setup

A near-field scanning system operating at frequencies between 6 and 12 GHz is used to map the electric field (see Fig. 5(c)). It is composed of two parallel conducting plates, each with a diameter of 1 m and spaced 15 mm apart. This ensures that the only propagating mode is the fundamental TEM mode for frequencies bellow 10 GHz. The source is an X-band waveguide (single-mode operational frequency range: 8.2–12.4 GHz; cutoff frequency: 6.557 GHz) fixed to the rotating bottom plate. It produces a cylindrical wave from its end, similar to a line source excitation. Microwave absorbers are added at the boundaries of the system in order to reduce scattering and reflections from the edges. Holes have been drilled along the radius of the fixed top plate with a spatial resolution of 5 mm. A monopole probe connected to the first port of a vector network analyzer (VNA) measures the amplitude and phase of the S-parameters with respect to the source waveguide connected to the other port. By varying the position of the probe along the radius and rotating the bottom plate, we obtain a 2D polar map of the near field. A computer controls and synchronizes the motor and the VNA through a LabView code.