OPEN

SUBJECT AREAS: THERMODYNAMICS APPLIED PHYSICS

Received 22 January 2015

Accepted 23 March 2015

Published 30 April 2015

Correspondence and requests for materials should be addressed to M.F. (mohamed. farhat@kaust.edu.sa)

Thermal invisibility based on scattering cancellation and mantle cloaking

M. Farhat1, P.-Y. Chen2, H. Bagci1, C. Amra3, S. Guenneau3 & A. Alu4

1Division of Computer, Electrical, and Mathematical Sciences and Engineering, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia, 2Department of Electrical and Computer Engineering, Wayne State University, Detroit, Michigan 48202, U.S.A, 3Aix-Marseille University, CNRS, Centrale Marseille, Institut Fresnel, Campus universitaire de Saint-Jerome, 13013 Marseille, France, 4Department of Electrical and Computer Engineering, The University of Texas at Austin, Austin, TX, 78712, U.S.A.

We theoretically and numerically analyze thermal invisibility based on the concept of scattering cancellation and mantle cloaking. We show that a small object can be made completely invisible to heat diffusion waves, by tailoring the heat conductivity of the spherical shell enclosing the object. This means that the thermal scattering from the object is suppressed, and the heat flow outside the object and the cloak made of these spherical shells behaves as if the object is not present. Thermal invisibility may open new vistas in hiding hot spots in infrared thermography, military furtivity, and electronics heating reduction.

The realization of electromagnetic invisibility cloaks1-4 is undoubtedly one of the most exciting and challenging applications of metamaterials5. In the previous decades, thanks to the astonishing development of micro- and nano-fabrication and 3D printing, this goal has got closer to reality. In 2005, Alu and Engheta proposed a transparency device that relies on the so-called scattering cancellation technique (SCT)6. This mechanism consists of using a low or negative electric permittivity cover to cancel the different scattering multipoles of the object to hide. This class of cloaking devices has been shown to be quite robust to changes in the geometry of objects and the frequency of operation7-9. Moreover, a recent experimental study has shown that these cloak designs can actually be realized at microwave frequencies10. Applications in furtivity, non-invasive sensing, and probing can be envisaged11,12, opening new directions in medicine, defense, and telecommunications. Recent findings also suggest that objects can be made invisible using the mantle cloaking technology, where a metasurface can produce similar effects with a simpler and thinner geometry. This is achieved by tailoring the surface current on the metasurface and consequently the phase of re-radiated fields13-16. It should also be mentioned here that other cloaking techniques have been put forward in the recent years based on various concepts such as conformal mapping1, transformation optics2,3,17, homogenization of multistructures18,19, active plasmonic cloaks20, anomalous localized resonances21, and waveguide theory22.

The concept of invisibility has been extended to other realms of physics. Cloaks capable of hiding objects from

acoustic waves23-25, surface water waves26, flexural bending waves27, seismic waves28,29, quantum matter waves30,31

and even diffusive light propagation32,33 have been developed. And more recently, after the seminal work of Guenneau et al.34, invisibility cloaks for heat waves has become another exciting venue for cloaking applica-tions35-37. Thermal cloak designs inspired by transformation optics2 have been subsequently proposed38-40 to control the flow of heat in metamaterial structures. Their experimental validation followed shortly41-44. Thermal cloaking may find interesting applications in modern electronics. It can be used to reduce the heat diffused from computers or to protect a specific nano-electronic component by re-directing the flow of heat. This technique can also be used for isolation in buildings to reduce the consumption of energy required in heating or cooling.

In this paper, we propose to use the concept of scattering cancellation to generate the invisibility effect for heat diffusion waves. The peculiarity of our cloak is that, unlike earlier designs, we consider both static and time-harmonic dependence (note that time-harmonic heat sources can be generated using pulsating lasers45). This scenario requires cancellation of two scattering orders for small objects, i.e. the monopole and dipole ones, corresponding to the specific heat capacity and the heat conductivity, respectively. Numerical simulations confirm that a scattering reduction of over 40 dB can be obtained for optimized cloak parameters. Additionally, it is shown that the proposed cloak suppresses both the near and far heat fields. We also demonstrate that coating an object with an ultra-thin layer or thermal metasurface is a viable way for scattering reduction (mantle cloaking).

Results

Heat diffusion waves and their dispersion relation. Using the first principle of thermodynamics in a closed system46, one can show that in the absence of radiation and convection, the temperature of a physical system obeys the Fourier relation V ■ jQ+pc8T/8t_Q. Here, Jq , p and T represent the density of heat flux (heat flow per unit surface per unit time), the density of the fluid, and the temperature field, respectively. c is the specific heat capacity and Q denotes the heat energy generated per unit volume per unit time (Fig. 1). Using the Fourier law, i.e. the linear and instantaneous relation Jq_—kVT, where k is the heat conductivity of the medium, one can derive,

pc~V-(kVT) + Q.

Cloak ■ Obstacle

Figure 1 | Thermal scattering problem. (a) Cross-sectional view of the heat transfer scenario, with the object to cloak in the middle. (b) Cross-sectional view of the cloaked object.

values of the shell conductivity and the specific heat capacity. First, a spherical object centered at the origin of a spherical coordinate system is considered. Two parallel plates set at different temperatures T1 < T2, generate a heat flux (plane heat diffusion wave) that impinges on the scattering object [Fig. 1(a)]. In this first section, the case of static (steady-state) regime is considered, i.e. 8T/8t_0. So Eq. (1) is simplified to AT+Q/k_0. The scalar temperature field T in the different regions of space can be expressed in spherical coordinates (Fig. 1) as,

T(r, h) = AlrlPl(cos h), 0 < r < a1,

T (r, h) = J2 [Blrl + C/r-(/+1)] Pl(cos h), ai < r < a2, (4)

For a constant conductivity and in the absence of heat sources, Eq. (1) simplifies to 8T/8t_K/(pc)AT. To solve this equation, one can assume that T(r,t) _ <(T(r,t)), with T(r,t) _ T0eikr—ivt, where k is the wave number of the pseudo diffusion plane wave and v its angular frequency. This ansatz is valid, only because Eq. (1) is a linear equation, meaning that T is a solution, if and only if, T is a solution. The dispersion relation of heat diffusion waves is thus ;®_k2K/pc. If one assumes that v is real, then k_+(1 + ;)/d, with d _ \J2k/(pc®). The general solutions are thus attenuated diffusing plane waves. Now, under the assumption of time-harmonic dependence e —ivt, generated for instance by a pulsating laser, and constant conductivity, Eq. (1) simplifies to

AT + k2T _—Q/k ■ (2)

For the structure in Fig. 1, Eq. (1) is supplied with two boundary conditions that should be satisfied at the surface of both spherical object and the cloak. Across the boundaries r_a1 and r_a2, we have the continuity of the temperature and the density of heat flux,

Le. TUa—2 _T Ir-aV and (K8T/8n)|r_a" _ (k8T/8n) jr_a+2,

where the signs + and — refer respectively to the inner and outer regions, and 8=8n denotes the normal derivative, which only depends upon the radial coordinate in the case of circular objects. Here, a1 and a2 are the inner and outer radii of the shell. Moreover, [k;-, p;-, c;-, k;]i = 012 represent the conductivity, fluid density, heat capacity, and wave number in the background medium (r>a2), object (r<a1), and shell (a1<r<a2), respectively.

Scattering cancellation technique for heat diffusion waves: static regime. The aim of this study is to show that scattering from various spherical objects can be reduced drastically by carefully choosing the

T (r, h) = Yj W + Dlr{(l+1)] Pl(cos h), r > a2,

where P/(-) represents the Legendre polynomial of order For r?œ, T (r > a2 ) ~ — (Q/k0)r cos h, therefore Ej ~ — Q/k0 and all the other coefficients £/= j are zero. The remaining coefficients are obtained by solving the linear system

1 1 1/a1 0

ki k2 -2k2/a3 0

0 1 1/a2 1 =a2

B0 k2 -2k2/a2 2k0/a2

Bi Ci \dJ

Q/ko Q

which is obtained by applying the continuity conditions at the boundaries r_a1 and r_a2. The scattering cancellation condition is obtained by enforcing that the first scattering coefficient D1 is zero,

D _ _Q (K0 — K2)(k1 + 2K2) — c3(K1 — K2)(k0 + 2K2) _0 7 1 2k0 c3(K0 —K2)(k1 —K2) —(2k0 + K2XK1+2K2) '

where c_a1/a2. Solving Eq. (7) for k2 yields the value of the shell conductivity, which ensures that there is no temperature perturbation with a uniform temperature gradient, as if the object does not exist,

— [(K0 — 2k1)c3 + (2K0 — K1)] + J [(K0 — 2k1)c3 + (2k0 — K1)]2 + 8k0K1(y3 — 1)2

K2 _---■ (8)

4(c3 — 1)

Scattering cancellation technique for heat diffusion waves: time-harmonic regime. The scattering coefficients relate the scattered fields to the incident ones, and depend on the geometry of the object and the frequency ®. Moreover, for a given size a of the object, only contributions up to a given order 10 are relevant, since the amplitude of the scattering coefficients changes as o(k)a)2iz1. The incident heat excitation is an oblique plane diffusion wave, of incidence angle h, and is of the form eik°r cos h. In a spherical coordinate system, it can be expressed as

Tinc(r, h) ~ To £ ¿l(2/ + 1)ji(kor)Pi(cos h),

where j/ denotes the fth spherical Bessel function and T0 is the amplitude of the incident temperature field. The scattered field (r>a2) can be expressed in a spherical coordinate system as

Tscat(r, h) = T0 jr il(2l + 1)5/h|1)(k0r)P/(cos h) ,

where s, are the complex scattering coefficients and are spherical Hankel functions of the first kind. Therefore, the temperature field can be expressed in the different regions of space as

T(r, 8)~ToY^ aiji(ki r)Pi(cos h), 0 < r < a , (11)

T(r, d)~ToJ2 [biji (k2r) + clyl (k2r)\Pl (cos 8),ai<r<a2, (12)

T(r, h)~Tinc + Ts'

= To^2 ii(2i +1) \ji (kor) + slh<l1)(kor) P, (cos h), r > a2,

with ai and (bi, ci) complex coefficients of the temperature field inside the object and the shell, respectively. Applying the continuity conditions at the boundaries r~ai and r~a2 yields the different coefficients. In particular, si ~ — Ui/(Ui + iVi). Here Ui and Vi are given by the determinants

ji(hai) yi (Mi) —ji (kiai) 0

ji(k2a2) yi (k2a2) 0 ji(koa2)

K2k2j[ (k2ai ) K2k2y\(k2ai) — Kikij'l(kiai) 0

Kjkjkjaj) K2k2y'l(k2a2)

ji(k2a1 ) ji(k2a2 )

yi (k2 a1) yi (k2 a2)

—ji(k1ai) 0

K2k2j'(k2a1) K2k2y'(k2a1) -K1k1j'(k1a1) K2k2?''(k2a2) K2k2y'(k2a2) 0

Kokoji' (koa2) 0

yi(koa2)

Kokoy'(koa2)

fscat = dß\g(Ö)\2.

Here, dV is the incremental solid angle, in spherical coordinates, dV~2n sin hdh, and g(ff) is expressed as

g(h)~- fg (2i+1)siPi(cos h). k

Inserting Eq. (17) into Eq. (16) yields 4p ?

fscat=TnïE (2i+1)\si\2.

\ko\ i~o

(\so\2 + 3\s1\2).

Consequently, canceling these two modes, i.e. s0 ~0 and sj~0, will ensure that fscat<0, and the thermal scattering from the object will be suppressed. Namely, the SCT conditions on the parameters of the cloaking shell k2, p2c2, and a2 are

r2C2-roCo ~ ~y3, forso~o,

P2C2-P1C1 \a2y

(Ko-K2)(k1 +2K2) (K1 -K2)(Ko+2K2)

~y3, fors1 ~o.

The monopole SCT condition in Eq. (20), depends only on the product of the density and the specific heat of the shell, and the ratio of radii of the object and the shell y. Similarly, the condition in Eq. (21) depends only on the conductivity of the shell and y. By enforcing these two conditions, the total scattering from the spherical object can be suppressed in the quasistatic limit.

Figures 2(a) and 2(b) illustrate numerical solutions to Eqs. (20) and (21), where the variation of the relative specific heat capacity r2c2/p0c0 and the relative heat conductivity k2/k0 are plotted versus y and pjCj/poCo and Kj/k0, respectively. From the solution of Eq. (20), given in Fig. 2(a), one can see that the relative specific heat capacity of the shell r2c2/r0c0, given here in logarithmic scale, takes positive and negative values, depending on y and the heat capacity of the object. The red line represents the curve obeying the equation y3pjci/p0c0 ~ 1 implying r2c2/r0c0~0. The specific heat capacity takes negative (positive near-zero) values above (below) this curve. From the solution ofEq. (21), given in Fig. 2(b), it can be seen that the required relative heat conductivity of the shell k2/k0 needs to be almost always negative, for varying y and k1/k0. However, for an object with small heat conductivity and small radius [lower part of Fig. 2(b), in blue color], the required shell conductivity is close to zero. In fact, from Eq. (21), one can derive that for the negative solution ofEq. (21)

K2/Ko ~ -

The scattering cross-section (SCS) fscat is a measure of the overall visibility of the object to external observers. It is obtained by integrating the scattering amplitude g(h), defined such that Tscat/T0<g(h)/reik°r, for r??,

-^/a2 + 8K1/Ko(1-y3) 4(1 -y3)

In the quasistatic limit (long diffusion length k0a = 1), only few scattering orders contribute to the overall scattering cross-section, namely the first two orders (l~0 for the monopole, and l ~ 1 for the dipole mode, unlike in the electrodynamic case, where the first dominant mode is the dipole one). In this scenario, one has

where a ~ 2 + y3 — (1 + 2y3)k1 /k0. It can be clearly seen that for positive conductivities of the object, the condition k2/k0 <0 has to be satisfied to achieve the optimal heat cloaking effect.

Let us move now to the analysis of a specific scenario, where the heat scattering of a spherical object is characterized. The relative specific heat capacity of the object is p1c1/p0c0 ~ 1.25 and its relative conductivity is k1/k0 ~ 0.5. The radius of the object a1 ~ 1, and the wave numbers are normalized to a1. The free space wave number is chosen as k0a1 ~0.5. This object is coated with a shell of outer radius a2~1.1a1. fscat of the total object-shell structure, defined in Eqs. (16)-(18), is normalized to the SCS of the bare object, and plotted against varying values of p2 c2 /p„ c0 and k2 / k0. The result is shown in Fig. 3(a) in logarithmic scale. The blue regions correspond to significant scattering reduction, whereas red regions correspond to enhanced scattering from the structure. It can be noticed that ranges of k2/k0 between 1 and 4, and p2c2/p0c0 between 0.05 and 0.5, are best for thermal scattering cancellation (now using the positive solution ofEq. (21), for practical realizations). The white dot has coordinates (3.1, 0.15) that correspond to the theoretical SCT condition obtained from Eqs. (20) and (21). It is also interesting to note that numerical simulations taking into account many scattering orders, give scattering reduction of 40 dB, sensibly around the same point.

These results show the importance of taking into account both the shell conductivity and specific heat capacity, in contrast to previous studies that only considered the effect of conductivity through the static analysis. This can be better understood from Figs. 3(b) and 3(c), where the normalized SCS is plotted versus k2 / k0 for various values of p2c2/p0c0, and versus p2c2/p0c0 for various values of k2/k0, respectively. The sensitivity to variations in k2/k0 is more evident from these figures, since a small variation from the optimum value

Figure 2 | Optimal cloaking parameters. (a) Relative specific heat capacity of the shell p2 c2 /p0c0 in logarithmic scale, versus the ratio c _ a1 /a2 and the relative specific heat capacity of the object p1 c1 / p0 c0. The color bar denotes the plot of log|p2 c2/p0 c01. (b) Relative heat conductivity of the shell k2/k0 in logarithmic scale, versus the ratio c_a1 /a2 and the relative heat conductivity of the object K1 /k0 . The color bar denotes the plot of log |k2/k0 |and the dashed black line represents log |k2/k0 | _ 0.

results in fast deterioration of the scattering reduction: when k2 / k0 is equal to 1 or 5, there is no dip in the SCS and the scattering is high, as can be seen from Fig. 3(c). The sensitivity to variations in p2c2 /p0 c0 is less important, as can be seen from Fig. 3(b), but it is important to choose values around those predicted by Eqs. (20) and (21). On the other hand, when V/ _ 0, peaks corresponding to modal resonances start appearing in the scattering cross-section (related to Fano-like response of the system due to interference between dark and bright scattering modes)47.

To better illustrate the efficiency of the proposed cloak, the far-field scattering patterns, i.e. the heat scattering amplitude |g(h)| in polar coordinates, in the x— y plane, are shown in Figs. 4(a) and 4(b). These figures demonstrate that the object is almost undetectable at all angles with scattering amplitude orders of magnitude lower than that of the bare object. As a result, there is no temperature perturbation around the object immersed in the thermal fields. To further demonstrate the functionality of the cloak, Figs. 4(c) and 4(d) plot the amplitude distribution of the scattered thermal field when the heat from the infinite sheet of oscillating heat source is impinging from left to right on the

structure, without and with the cloaking shell, respectively. When the object is cloaked, the field amplitude is constant everywhere in space in contrast to the case of the object without the cloak.

Discussion

Mantle cloaking for heat diffusion waves. As stated in the introduction, recent findings suggest that objects can be made invisible using the surface cloaking technology, where a metasurface may produce similar cloaking effects in a simpler and thinner geometry14-16. The ultrathin mantle cloak with an averaged surface reactance metasurface13 reduces the scattering from the hidden object, comparable to bulk metamaterial cloaks. The setup of the problem is similar to the previous section, except for the fact that scattering cancellation is achieved by a surface, instead of a shell. This is illustrated in the inset of Fig. 5(a). The impedance boundary condition results in jumps in the radial component of the density of heat flux, on the interface between the two media.

In what follows it is shown that the scattering from various spherical objects can be drastically reduced by choosing the appropriate

0.2 0.4 0.6 0.8 Relative specific heat capacity

Figure 3 | Thermal scattering reduction. (a) Normalized (analytical) SCS fscat in logarithmic scale, versus the relative heat conductivity k2/k0 and the relative specific heat capacity p2 c2/p0 c0. The white dot represents the position of optimized scattering reduction, with a value of 40 dB. The color bar denotes the plot of 10 log( f]cat / fscat), where the subscripts 1 and 2 refer to the scattering cross-section of the obstacle and cloaked structure, respectively. (b) Normalized SCS versus the relative heat conductivity for various values of the specific heat capacity p2 c2 / p0 c0. (c) Normalized SCS versus the relative specific heat capacity p2 c2/p0 c0 for various values of the relative heat conductivity k2 /k0.

■ yt

—> fK-

IL-.«-«-**.

_. ___>» / j.

I 3 ' v-

V- ^ f'

Figure 4 | Near and far-field characterization. Analytical scattering amplitude |g(h)|, given by Eq. (17), in polar coordinates, and in logarithmic scale (a) for the bare object with ki /k0 = 0.5 and (b) for the cloaked object, with k2/k0 = 3. i. Amplitude of the oscillating temperature in the near-field of (c) the same bare object of Fig. 4(a) and (d) the same cloaked object of Fig. 4(b) for k0 ai =0.5. Arrows show the direction of+T and the color bar denotes the plot of It/tinc|.

surface impedance, and thus their visibility to heat diffusion waves can be suppressed.

To design a mantle cloak, we keep the boundary conditions at r = ai same as those in the previous sections, while we replace the boundary conditions at r~a2 with

It should be noted that for the mantle cloak design considered here, k2 ~ k0 and p2c2 ~ p0c0. In Eq. (25), the dimensionless function y is defined as

T |r = a— T |r = az T1 r = a

ÔT Ln

= Zs—i T|r = £

/ ivZ— ' i

y=-K0kT = K0k0

For k0 a2 = i, the spherical Bessel functions take a simpler polynomial form, and the approximate cloaking condition in this limit can be written as

Eq. (24) is a surface impedance condition that implies a jump in the density of heat flux. Here, Zs ~Rs + iXs is the averaged surface impedance that relates the temperature to the density of heat flux on the surface.

Following the procedure described in the previous section, with these new boundary conditions, one can show that the lth spherical scattering harmonic can be suppressed, provided that the following determinant is canceled,

3y3vai

K0 Z 2ki Ki — K0

ji (k0ai) k0ji (k0 ai ) ji (ha2) ji (hai) + yji (k0 a2)

yi (k0ai) hy[ (k0 ai )/p0 yi(k0a2) yi (k0a2) + yyi (k0a2)

—ji (ki ai ) Ki/K0kiji (ki ai ) 0

ji (k0 a2 ) ji (k0 a2 )

. (25)

This clearly shows that by properly choosing the thermal surface reactance (expressed in units of J=(m2K)), it is possible to suppress the dominant multipolar scattering in the quasistatic limit.

Figure 5(a) plots the SCS versus Xs for cloaked objects with various y. The SCS of a bare object is plotted for comparison. For Xs ? ?,we notice that the metasurface does not reduce the scattering, consistent with the limit of no-surface. For specific values of Xs, however, a relevant scattering reduction is achieved, and this may be obtained for different values of a2, even in the limit of a cloak winding con-formal to the object (a2 ~ a1; y ~ 1).

Figure 5(b) plots the SCS versus the frequency for cloaked objects with a2 ~ a1 (conformal) and a2 ~ 1.1a1. We suppose here that the surface reactance does not vary with frequency and is given with Xs ~ 1.36 x 10 —7 for y = 1 and Xs ~1.62x10 —7 for y~0.9. The

Figure 5 | Thermal mantle cloaking. Normalized (analytical) SCS fscat versus Xs 10 — 7 for the object with normalized radius k0a1 _ 0-5 and relative heat conductivity K1 /k0 _0-1 for various values of the ratio y. (b) Normalized SCS fscat versus normalized frequency for the same object with Xs _L36x10 —7 and Xs _ 1 -62 x 10 — 7 for values of the ratio y _ 1 and y _ 0-9. The inset of Fig. 5(a) illustrates the object coated with a thermal metasurface, projected in the x— y plane. Amplitude of the temperature field on the x—y plane for the same object (c) with mantle cloak with Xs _ 1-62 x 10 — 7 and (d) without the cloak. The color bar denotes the plot of | T/ Tinc |.

SCS of uncloaked objects with radius r _ a and r _ 1 ■ 1 a are plotted for comparison. It is evident that excellent scattering reduction may be achieved over a large range of frequencies for both cases.

Figures 5(c) and 5(d) plot the amplitude of the temperature field scattered by a cloaked and uncloaked object, on the x—y plane at a time instant, respectively. When the object is cloaked, both forward and backward scattering almost vanish. This reduction of scattering is achieved due to the proper choice of the surface impedance, which restores almost uniform amplitude all around the cloak.

Summary. In conclusion, we have proposed an original route towards designing thermal cloaks based on the scattering cancellation technique. This technique is inspired by the plasmonic cloaking, which makes use of shells with induced negative polarization to suppress scattered electromagnetic fields. And contrary to invisibility cloaks based on transformation optics, SCT offers simple cloaking designs (without the need of anisotropy and inhomogeneity of the physical parameters).

One may envision that using this design may further make the thermal cloaking closer to its practical and feasible realization. We believe that such a structured cloak could be manufactured within current technology, having in mind some potential applications in invisibility, sensing and thermography. The range of industrial applications is vast, and our proof of concept should foster research efforts in this emerging area of thermal cloaks and metamaterials.

Methods

Analytical methods based on scattering Mie theory of spherical thermal scatterers are used to obtain the results presented in Figs. 2, 3, 4(a), 4(b), and 5. In the quasistatic limit, where the size of the object is much smaller than the wavelength and only the lowest-order Mie coefficients are kept, analytical formulas are obtained [Eqs. (20) and

(21)]. Those give results similar to the ones obtained from full Mie series solutions [Fig. 3(a)]. The results given in Figs. 4(c) and 4(d) are obtained using COMSOL Multiphysics software, which solves Eq. (2) with proper boundary conditions using a finite element scheme.

1. Leonhardt, U. Optical conformal mapping. Science 312, 1777-1780 (2006).

2. Pendry, J. B., Schurig, D., and Smith, D. R. Controlling electromagnetic fields. science 312, 1780-1782 (2006).

3. Schurig, D. eta/. Metamaterial electromagnetic cloak at microwave frequencies. Science 314, 977-980 (2006).

4. Guenneau, S. eta/. The colours of cloaks. J. Opt. 13, 024014 (2011).

5. Sihvola, A. Metamaterials in electromagnetics. Metamateria/s 1, 2-11 (2007).

6. Alu, A. and Engheta, N. Achieving transparency with plasmonic and metamaterial coatings. Phys. Rev. E 72, 016623 (2005).

7. Aluu, A. and Engheta, N. Plasmonic materials in transparency and cloaking problems: mechanism, robustness, and physical insights. Opt. Express 15, 3318-3332 (2007).

8. Muhlig, S., Farhat, M., Rockstuhl, C., and Lederer, F. Cloaking dielectric spherical objects by a shell of metallic nanoparticles. Phys. Rev. B 83, 195116 (2011).

9. Muhlig, S. eta/. A self-assembled three-dimensional cloak in the visible. Sci. Rep. 3, 2328 (2013).

10. Rainwater, D. eta/. Experimental verification of three-dimensional plasmonic cloaking in free-space. New J. Phys. 14, 013054 (2012).

11. Alu, A. and Engheta, N. Cloaking a sensor. Phys. Rev. Lett. 102, 233901 (2009).

12. Alu, A. and Engheta, N. Cloaked near-field scanning optical microscope tip for noninvasive near-field imaging. Phys. Rev. Lett. 105, 263906 (2010).

13. Munk, B. A. Frequency Selective Surfaces: Theory andDesign. (John Wiley & Sons, 2005).

14. Alu, A. Mantle cloak: Invisibility induced by a surface. Phys. Rev. B 80, 245115 (2009).

15. Chen, P. -Y., Farhat, M., Guenneau, S., Enoch, S., and Alu, A. Acoustic scattering cancellation via ultrathin pseudo-surface. App/. Phys. Lett. 99, 191913 (2011).

16. Farhat, M., Chen, P. -Y., Guenneau, S., Enoch, S., and Alu, A. Frequency-selective surface acoustic invisibility for three-dimensional immersed objects. Phys. Rev. B 86, 174303 (2012).

www.nature.com/scientificreports \ > _

17. Cai, W., Chettiar, U. K., Kildishev, A. V., and Shalaev, V. M. Optical cloaking with metamaterials. Nature photon. 1, 224-227 (2007).

18. Farhat, M., Guenneau, S., Movchan, A., and Enoch, S. Achieving invisibility over a finite range of frequencies. Opt. Express 16, 5656-5661 (2008).

19. Farhat, M. etaZ. A homogenization route towards square cylindrical acoustic cloaks. New J. Phys. 10, 115030 (2008).

20. Vasquez, F. G., Milton, G. W., and Onofrei, D. Active exterior cloaking for the 2d laplace and helmholtz equations. Phy. Rev. Left. 103, 073901 (2009).

21. Nicorovici, N., Milton, G. W., McPhedran, R. C., and Botten, L. C. Quasistatic cloaking of two-dimensional polarizable discrete systems by anomalous resonance. Opt. Express 15, 6314-6323 (2007).

22. Smolyaninov, 1.1., Smolyaninova, V. N., Kildishev, A. V., and Shalaev, V. M. Anisotropic metamaterials emulated by tapered waveguides: application to optical cloaking. Phy. Rev. Lett. 102, 213901 (2009).

23. Torrent, D. and Sánchez-Dehesa, J. Acoustic cloaking in two dimensions: a feasible approach. New J. Phys. 10, 063015 (2008).

24. Cummer, S. A. et aZ. Scattering theory derivation of a 3d acoustic cloaking shell. Phy. Rev. Lett. 100, 024301 (2008).

25. Norris, A. N. Acoustic cloaking theory. Proc. R. Soc. A-Math. Phys. Eng. Sci. 464, 2411-2434 (2008).

26. Farhat, M., Enoch, S., Guenneau, S., and Movchan, A. Broadband cylindrical acoustic cloak for linear surface waves in a fluid. Phy. Rev. Lett. 101, 134501 (2008).

27. Farhat, M., Guenneau, S., and Enoch, S. Ultrabroadband elastic cloaking in thin plates. Phy. Rev. Lett. 103, 024301 (2009).

28. Brun, M., Guenneau, S., and Movchan, A. B. Achieving control of in-plane elastic waves. AppZ. Phy. Lett. 94, 061903 (2009).

29. Brûlé, S., Javelaud, E., Enoch, S., and Guenneau, S. Experiments on seismic metamaterials: Molding surface waves. Phys. Rev. Lett. 112, 133901 (2014).

30. Zhang, S., Genov, D. A., Sun, C., and Zhang, X. Cloaking of matter waves. Phy. Rev. Lett. 100, 123002 (2008).

31. Greenleaf, A., Kurylev, Y., Lassas, M., and Uhlmann, G. Approximate quantum cloaking and almost-trapped states. Phy. Rev. Lett. 101, 220404 (2008).

32. Schittny, R., Kadic, M., Bueckmann, T., and Wegener, M. Invisibility cloaking in a diffusive light scattering medium. Science 345, 427-429 (2014).

33. Schittny, R. et aZ. Transient behavior of invisibility cloaks for diffusive light propagation. Optica 2, 84-87 (2015).

34. Guenneau, S., Amra, C., and Veynante, D. Transformation thermodynamics: cloaking and concentrating heat flux. Opt. Express 20, 8207-8218 (2012).

35. Alù, A. Thermal cloaks get hot. Physics 7, 12, Feb (2014).

36. Maldovan, M. Sound and heat revolutions in phononics. Nature 503, 209-217 (2013).

37. Leonhardt, U. Applied physics: Cloaking of heat. Nature 498, 440-441 (2013).

38. Chen, T., Weng, C. -N., and Chen, J. -S. Cloak for curvilinearly anisotropic media in conduction. AppZ. Phys. Lett. 93, 114103 (2008).

39. Guenneau, S. and Amra, C. Anisotropic conductivity rotates heat fluxes in transient regimes. Opt. Express 21, 6578-6583 (2013).

40. Moccia, M., Castaldi, G., Savo, S., Sato, Y., and Galdi, V. Independent manipulation of heat and electrical current via bifunctional metamaterials. Phys. Rev. X 4, 021025 (2014).

41. Narayana, S. and Sato, Y. Heat flux manipulation with engineered thermal materials. Phys. Rev. Lett. 108, 214303 (2012).

42. Schittny, R., Kadic, M., Guenneau, S., and Wegener, M. Experiments on transformation thermodynamics: molding the flow of heat. Phys. Rev. Lett. 110, 195901 (2013).

43. Han, T. et al. Experimental demonstration of a bilayer thermal cloak. Phys. Rev. Lett. 112, 054302 (2014).

44. Xu, H., Shi, X., Gao, F., Sun, H., and Zhang, B. Ultrathin three-dimensional thermal cloak. Phys. Rev. Lett. 112, 054301 (2014).

45. Shendeleva, M. L. Thermal wave reflection and refraction at a plane interface: Two-dimensional geometry. Phys. Rev. B 65, 134209 (2002).

46. Reif, F. Fundamentals of statistical and thermal physics. (Waveland Press, 2009).

47. Argyropoulos, C., Chen, P. -Y., Monticone, F., D'Aguanno, G., and Alù, A. Nonlinear plasmonic cloaks to realize giant all-optical scattering switching. Phys. Rev. Lett. 108, 263905 (2012).

Acknowledgments

This work is partially funded by King Abdulaziz City for Science and Technology (KACST) TIC (Technology Innovation Center) for Solid-state Lighting at KAUST. P.-Y.C. would like to acknowledge fruitful discussion with David Piech. S.G. would like to acknowledge a funding of the European Research Council through ERC grant ANAMORPHISM.

Author contributions

M.F. and P.-Y.C conceived the idea of this study. M.F. performed numerical simulations and wrote the manuscript. P.Y.C., H.B., C.A., S.G., and A.A. contributed to the analysis of the results and reviewed the manuscript. S.G. and A.A. supervised the project.

Additional information

Competing financial interests: The authors declare no competing financial interests.

How to cite this article: Farhat, M. et a/. Thermal invisibility based on scattering cancellation and mantle cloaking. Sci. Rep. 5, 9876; DOI:10.1038/srep09876 (2015).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder in order to reproduce the material. To view a copy ofthis license, visit http://creativecommons.org/licenses/by/4.0/