Manipulating the Loss in Electromagnetic Cloaks for Perfect Wave Absorption

We examine several ways to manipulate the loss in electromagnetic cloaks, based on transformation electromagnetics. It is found that, by utilizing inherent electric and magnetic losses of metamaterials, perfect wave absorption can be achieved based on several popular designs of electromagnetic cloaks. A practical implementation of the absorber, consisting of ten discrete layers of metamaterials, is proposed. The new devices demonstrate super-absorptivity over a moderate wideband range, suitable for both microwave and optical applications. It is corroborated that the device is functional with a subwavelength thickness and, hence, advantageous compared to the conventional absorbers.

Metamaterials can be generally defined as a class of "artificial" media that exhibit extraordinary electromagnetic properties not found in nature. For example, materials with negative refractive index can be designed [1], which can theoretically achieve infinite subwavelength resolution [2]. Recently, a dispersive electromagnetic cloak was experimentally demonstrated [3] by achieving the required material anisotropy [4]. Such structures, typically constructed from periodically-placed subwavelength unit cells based on highly conductive metals over dielectric substrates, can be analyzed with the effective medium theory [5]. They are defined by complex electromagnetic parameters: the frequency dependent permittivity ε(ω) = ε 1 (ω) + jε 2 (ω) and permeability µ(ω) = µ 1 (ω) + jµ 2 (ω). So far, metamaterial research has mainly focused on the real parts (ε 1 , µ 1 ) of the material parameters in order to design negative-index and cloaking devices. However, the imaginary parts (ε 2 , µ 2 ) of the material parameters (which characterize the losses of the medium) can also have interesting potential applications, such as the design of more efficient thermal imagers and novel absorbers. In this letter, we take advantage of the overlooked lossy properties of metamaterials, which are widely regarded as a drawback in the design of these exotic devices [6].
Three absorber designs are proposed, by manipulating both electric and magnetic losses of subwavelength cylindrical electromagnetic cloaking coatings. The first two designs utilize the material parameter sets of the matched reduced [7] and ideal [8] electromagnetic cloaks. The third design consists of a more realistic discrete ten-layer structure, based on an approach presented in [3], where only one material parameter is radius-dependent and the others have constant and non-dispersive values. The experimentally verified structure in [3] is generally regarded as an imperfect cloaking device; however, as we demonstrate below, with a proper manipulation of material losses a perfect absorber can be constructed. Specifically, by calculating the field patterns and analyzing the scattering coefficients, we demonstrate that these new devices can achieve perfect wave absorption over a narrow bandwidth and moderate performance over a broad bandwidth. Additionally, they are matched to free space without further parameter tuning and they are found to perform significantly better when compared to conventional dielectric absorbers. The effects are analyzed using a radially-dependent dispersive Finite-Difference Time-Domain (FDTD) method and scrutinized with existing analytical solutions of the cloaks [9,10].
The proposed absorbers are based on a different concept from those used in the recently-demonstrated perfect microwave metamaterial [11] and terahertz metamaterial [12,13] absorbers. The latter devices operate correctly only when there is strong coupling at the appropriate resonant mode of each metamaterial unit cell. As a result, careful tuning of the complex parameters in the device is essential, in order to achieve zero backscattering and maximum absorptivity. However, the proposed absorbers are independent from resonant modes [8] and, hence, can operate over a relatively wide frequency range.
For the two-dimensional (2-D) FDTD simulations presented here, a perfect electric conductor (PEC) cylinder is surrounded by a sub-wavelength absorption coating. Without loss of generality, a TE plane wave is incident and only three field components are non-zero: E x , E y and H z . The computational domain of the infinite (towards z-direction) cylindrical "cloaking" absorber can be seen in Fig. 1(a).
The performance of an absorber can be characterized by its absorptivity A = 1 − |S 11 | 2 − |S 21 | 2 , where S 11 and S 21 are the reflection and transmission coefficients of the device, respectively. An ideal absorber is characterized by an absorptivity equal to unity, which directly leads to no reflection (S 11 = 0) and no transmission (S 21 = 0) through the device. In addition, the scattering coefficient σ S of an absorber can be calculated with reference to free space, with no obstacles present. It is given by the formula where H z is the complex magnetic field distribution in the "cloaking" absorber and H f r z is the complex magnetic field distribution in the free space. When σ S = 0, the surrounding field of the absorber is equal to the field in free space, i.e. the structure is totally reflectionless. When σ S = 1, the field is entirely dissipated inside the device (H z = 0) and the radiation is not transmitted through the absorber. This last condition that S 21 → 0 in the transmitted region of the absorber (y > 0 in Fig. 1(a)), combined with the reflectionless property S 11 → 0 of the cloaking material in the reflected re- gion (y < 0 in Fig. 1(a)), is the ideal condition in order to achieve a perfect absorber with A → 100%. The "cloaking" absorbing structures proposed in this Letter can easily achieve A > 80%, which is rapidly enhanced by increasing the loss factor in some suitable way. The first design (reduced "cloaking" absorber) utilizes the material parameters of the matched (to free space) reduced cloak [7]: where R 1 and R 2 are the inner and outer radii of the absorber, respectively, with R 1 < r < R 2 . The permittivity value is gradually changing with the radius r of the device and it can be obtained that 0 ≤ ε r < 1 and ε φ , µ z > 1 for all values of r. The ε r parameter is mapped with the Drude dis- where ω p is the plasma frequency and γ is the collision frequency characterizing the losses of the dispersive material. Furthermore, ε φ and µ z are simulated with a conventional dielectric/magnetic model: ε(ω) = ε + σ ω , where σ measures the conductive/magnetic losses. A thorough description of the radially-dependent dispersive FDTD algorithm employed to simulate the proposed absorber can be found in [14,15].
In the case of Drude model mapping of the material parameter (ε r < 1), the lossy parameter can be presented in an alternative way:ε r = ε r (1 −  tan δ), where the radiusdependent parameter ε r is given from the reduced parameter set and tan δ is the loss tangent of the lossy material. If the previous formula is substituted in the Drude model and tan δ is assumed constant, the radially-dependent plasma frequency becomes ω p (r) = (1 − ε r )ω 2 + ε r ωγ tan δ and the collision frequency becomes γ(r) = εrω tan δ (1−εr ) . Similarly, the conductivity of the conventional dielectric/magnetic model is given by (using ε φ as an example) σ = ε φ ω tan δ, which is constant because ε φ and µ z have a fixed value. Note that, the concept of artificial loss functions was also introduced in [16], but for a different coordinate transformation function. It was used to create an absorbing boundary condition (like perfectly matched layers) in computational electromagnetics, whereas the aim of our proposed device is to be implemented in practice for engineering applications.
The device is tested at a frequency of 2 GHz, which is used throughout the Letter. Other operating frequencies can also be chosen by simply adjusting the material parameters due to the frequency-independent nature of coordinate transformation functions [4]. The dimensions of the device, in terms of the free-space wavelength λ, are R 1 = λ and R 2 = 4λ 3 . The results of the normalized real part of the magnetic field amplitude distribution when the steady-state is reached are shown in Fig. 1(b), for a device with tan δ = 0.5. The computed scattering coefficients of the device as a function of the loss tangent tan δ can be seen in Fig. 1(c), where the scattering coefficient of a bare PEC cylinder is also shown. It is observed that the backscattering of the "cloaking" absorber -i.e. angle equal to 180 • -goes to zero for all loss tangents; simultaneously, the reflection coefficients approach zero (S 11 → 0) as tan δ increases. It should also be mentioned that the backscattering coefficient of the reduced "cloaking" absorber utilized here is at least 20 dB lower, when compared to an absorber based on the simple reduced set used in [3] (not shown).
We further verify this absorber concept by introducing a second design which is based on the material parameter set of an ideal electromagnetic cloak [8] (ideal "cloaking" absorber), without altering the geometry. All the parameters are now radius-dependent: r . The scattering coefficients of the device, varying with the losses, can be seen in Fig. 1(d). The backscattering of the ideal set device has slightly improved values, compared with the reduced "cloaking" absorber in Fig. 1(c). Note that, when the losses are increased, the exact backscattering is still low, but there are weak backward reflec-tions at the sides of both absorbers, which will be discussed later in this Letter.
A big shadow which tends to increase with the loss tangent is cast at the back of both the proposed devices (especially between the angles −30 • to 30 • ), as can be seen in Figs. 1(c), (d), which is the ideal scattering pattern for an efficient absorber [17]. It is interesting that, for the low loss tangent of tan δ = 0.05 the absorption of the reduced "cloaking" absorber ( Fig. 1(c)) is better than the absorption of the ideal "cloaking" absorber ( Fig. 1(d)). For higher losses (tan δ = 0.5), it is observed that the scattering coefficients tend to one, especially inside the window −30 • to 30 • .
In the previous design both conductivities and collision frequencies are required to be continuously radius-dependent, which makes a practical implementation of this metamaterial structure challenging. In order to check the performance of a more realistic absorber, the device is designed with a discrete ten layered structure (discrete reduced "cloaking" absorber), similar to the cloak constructed in [3]. Again the reduced parameter set is used [7] and only the radius-dependent parameter is discretized to ten different values, one for each layer. The other parameters are kept constant, as before, and the loss tangent is chosen tan δ = 0.1. The scattering pattern of the device is then computed, it is given in Fig. 1(e) and is compared with the continuous reduced and ideal parameter set "cloaking" absorbers. It is seen that the patterns of the three absorbers are following a similar trend, especially at the angles between 0 • -90 • (shadow). The larger difference is the backscattering coefficient (angle of 180 • ), which is less in the case of ideal "cloaking" absorber. For tan δ = 0.1, it is calculated that the device achieves an absorptivity of at least A = 80% within a frequency range of 400 MHz. If the losses are increased, absorptivity values of A > 90% can be easily obtained, which suggests that the proposed structure can be utilized as a perfect metamaterial absorber.
Next, we demonstrate the absorber's performance for different coating thicknesses, with R 1 = 2λ 3 and R 2 = 4λ 3 . The loss parameters for these structures are graphically depicted in Figs. 2(a),(b), assuming tan δ = 0.1. Note that at the inner surface of both absorbers (r = R 1 ), the electric collision frequency γ e = 0, because ε r = 0 at r = R 1 . It can be seen that there are no infinite loss values for the reduced "cloaking" absorbers, which makes it easier to be practically implemented. Moreover, the conductivities of this particular absorber are constant, which leads directly to the usage of only one type of dielectric material as the substrate to fabricate split-ring resonators (SRRs), in order to achieve the required losses. The bandwidth response of the ideal and reduced "cloaking" absorbers can be seen in Fig. 2(c). The lowest backscattering is obtained for both sets at the central frequency of 2 GHz, while the backscattering is stronger at other frequencies due to the dispersive nature of the absorber. Based on this behavior, the proposed devices can also have moderate wideband applications, e.g. reducing the radar cross section (RCS) of an object.
Next, the "cloaking" absorber coating ( Fig. 1(a)) is re- placed with conventional absorbing materials, in order to evaluate their individual performances. A conventional dielectric absorber with ε = 2 and a matched conventional absorber with ε = 1 are compared to the discrete reduced "cloaking" absorber. All the devices have identical loss tangents of tan δ = 0.1 and the same dimensions (R 1 = 2λ 3 , R 2 = 4λ 3 ). The obtained scattering coefficient patterns of these absorption coatings are shown in Fig. 2(d), where it is observed that the proposed metamaterial device performs better when compared to the conventional ones. Firstly, minimal backscattering is only achieved with the proposed "cloaking" structure, which is not the case even for the matched conventional absorber. In addition, the absorption (shadow) of the "cloaking" absorber is more uniform and more solid than the other two cases. Furthermore, when compared with the absorption of the matched dielectric absorber, the shadow is notably larger in a broader angle range (between −90 • to 90 • ).
One of the main physical reasons that explains the superior performance of the metamaterial absorber designs is the bending of the field wavefronts within the absorption coating. Due to the anisotropy of the material parameters, the impinging field energy is guided over longer distances around the cloaked object, thus, dissipating gradually inside the thin sub-wavelength coating. The proposed 'cloaking" absorbers inherit this unique property from the electromagnetic cloaks, an effect that is not possible in conventional absorbers of the same thickness, since the field energy is limited to straight line propagation. Similar reasoning explains the slightly improved scattering coefficient of the ideal 'cloaking" absorber compared to the reduced ones ( Fig. 1(e)). The reduced parameter sets allow imperfect field bending and, thus, some energy scatters off the inner core of the absorption coating [18].
The proposed absorber designs also inherit the reflectionless property of the cloaks with both ideal and reduced parameter sets. The characteristic impedance of the absorbers at the outer boundary is ε0 , which is the free space wave impedance . As a result, the absorbers are matched to the surrounding free space and there is minimal backscattering. This is indeed observed in Fig.  1(c), (d), where the field at the backscattering angle of 180 • remains low as the losses are being increased, which is advantageous for the efficient use of the devices.
Finally, we would like to address the weak omnidirectional scattering that was previously observed in Figs. 1(c), (d) as the losses are increased, especially between the angles of 90 • to 120 • and 240 • to 270 • . In order to explain this particular response of the device, the cylindrical wave expansion is used [9]. The total magnetic field inside and outside of the cylindrical 2-D "cloaking" absorber is given, respectively, by H in H l are the l th -order Bessel and Hankel functions of the first kind, respectively. In order to computer the fields, the expansion coefficients a inc l , a sc l , a 1 l , a 2 l need to be calculated by applying the proper boundary conditions at the inner and outer surfaces of the absorption coating. It was obtained in [9] that a 2 l = 0 for both the ideal and reduced parameter sets, which also holds true in the proposed absorbers.
Nevertheless, for the scattering expansion coefficients in the current "cloaking" absorber a sc l = a 2 l = 0. This is a direct result from the differences existing between the phase variations in the arguments of the Bessel and Hankel functions in the above expansions. More precisely, at the outer surface r = R 2 of the device, the phase variation of the total incident wave is k 0 R 2 , whilst the phase variation of the total field inside the absorber is equal to k 1 (R 2 − R 1 ). However, unlike the case of the ideal and reduced cylindrical lossless cloaks, these two quantities are not equal when losses are introduced. The phase variation of the field inside the absorber is: , which is directly derived from the ideal and reduced parameters. Meanwhile, the phase variation of the field outside the absorber is k 0 R 2 = √ ε 0 µ 0 R 2 . This phase mismatch which occurs as the loss tangent is increased from zero, gives rise to the weak scattering observed at the sides of the absorber in Figs. 1(c), (d). Nevertheless, this scattering is compensated by the bigger shadow achieved when more losses are introduced and the overall performance of the "cloaking" absorber is mostly unaffected by this slight imper-fection.
The proposed discrete reduced set absorbing device can be constructed with ten concentric arrays of split-ring resonators or electric ring resonators printed on a dielectric substrate, as was firstly demonstrated in [3]. The varying collision frequency values of the "cloaking" absorber can be controlled with different practical techniques. For the microwave frequency regime, where the dielectric losses are dominant [19], different substrate dielectric materials can be chosen, or lossy nanoparticles can be introduced into the substrate. For higher frequencies, such as infrared, the ohmic losses are more significant [19], which leads to the solution of changing the metallic material of the resonator particles (for example, gold is more lossy than silver at infrared frequencies) or varying their dimensions.
To conclude, a set of novel metamaterial absorbers is proposed based on the manipulation of losses inside electromagnetic cloaks, yielding absorptivity values larger than 80%. It is explained that a discrete reduced "cloaking" absorber can be constructed with existing metamaterial technologies. These devices could be used as narrow-band bolometers, due to their dispersive nature; they also exhibit a moderate wideband response, which is ideal for RCS reduction. Another useful application is the total isolation of an antenna placed inside the absorber from low-frequency surrounding noise. It is demonstrated that the devices are reflectionless and efficient in terms of absorption. Moreover, they have subwavelength thickness, which is highly desirable for absorbers operating in VHF or UHF bands. Finally, such absorbers consist one more useful addition to the overlooked impact of metamaterials to absorbing technology.