Mantle cloaking due to ideal magnetic dipole scattering

One of the most exciting applications of metaparticles and metasurfaces consists in the magnetic light excitation. However, the principal limitation is due to parasitic extra multipoles of electric family excited in magnetic dipole meta-particles characterized by a radiating nature and corresponding radiating losses. In this paper, we propose the “ideal magnetic dipole” with suppressed additional multipoles except of magnetic dipole moment in the scattered field from a cylindrical object by using mantle cloaking based on metasurface and on anapole concept. The considered metasurface consists of a periodic width modulated microstrip line, with a sinusoidally shaped profile unit cell printed on a dielectric substrate.

Mantle cloaking due to ideal magnetic dipole scattering Barbara cappello 1 , Anar K. ospanova 2,3 , Ladislau Matekovits 1* & Alexey A. Basharin 2,3 one of the most exciting applications of metaparticles and metasurfaces consists in the magnetic light excitation. However, the principal limitation is due to parasitic extra multipoles of electric family excited in magnetic dipole meta-particles characterized by a radiating nature and corresponding radiating losses. in this paper, we propose the "ideal magnetic dipole" with suppressed additional multipoles except of magnetic dipole moment in the scattered field from a cylindrical object by using mantle cloaking based on metasurface and on anapole concept. the considered metasurface consists of a periodic width modulated microstrip line, with a sinusoidally shaped profile unit cell printed on a dielectric substrate.
Despite of the progress in electromagnetics, the light-matter interaction is still associated especially with electric component of light, while magnetic component is suppressed down to negligible level, particularly at more anticipated optical frequencies. On the other hand, magnetic light-matter interaction is promising playground for unusual effects like negative refraction [1][2][3][4] , magnetoinductive waves 5 , fluorescent microscopy 6,7 , nanoscale imaging and others 8,9 . Therefore, magnetic response of subwavelength particles contributes to optical magnetism leading to effects such as magnetic light 10 , magnetic nanoantennas [11][12][13][14][15] and magnetic Purcell effect [16][17][18][19] . For this sake, recent works propose various techniques for artificial magnetism that enable strong magnetic field localization even for nonferromagnetic particles that is usually several times weaker than the electric field component counterpart [20][21][22][23][24] . The first realization of artificial magnetic dipole (MD) has been demonstrated in metal split-ring resonator (SRR), thereafter becoming a basic element of metamaterials 3,25,26 . Overall, strongly concentrated magnetic field can be excited by a circular current oscillating within the SRR and mimicking a MD. The idea of SRR is still of high demand; however, its application is restricted by the impossibility of implementation in visible frequency range due to intrinsic metallic particles with high Joule losses 26 . A second limitation is due to the excitation of parasitic extra multipoles of electric family characterized by a radiating nature and corresponding radiating losses. Thus, researchers are approaching to achieve an ideal magnetic dipole without additional multipoles (except of magnetic dipole moment) in the system.
The idea of the ideal MD scatterer can be elegantly described by multipole decomposition theory of the scattering from resonant particles. This approach allows studying the radiation properties by description of the scattering cross-section as a sum of electric a E (l, m) and magnetic a M (l, m) scattering coefficients. Indeed, the resulting scattering efficiency Q sca is given by the superposition of electric and magnetic scattered multipoles: where l is the number of spherical harmonics defining multipoles order. In this way, we can calculate the scattering intensity of electric multipoles described by the first term and of magnetic multipoles, expressed as second term in Eq. (1), directly through the induced currents inside the particles 27 . Inspired by the ideal MD scatterer approach, we should find the current distribution with electric multipoles a E (l, m) scattering intensity tending to zero. In this case, Q sca is almost described by the second term in Eq. (1), i.e., a M (l, m), given by the magnetic dipole mode. One of the elegant solutions is anapole mode excitation, defined as electric type scattering elimination by toroidal dipole moment at the same frequency 28,29 . Consequently, the spectral overlapping of magnetic and anapole modes provides the pure MD scattering in far-field zone. This concept has recently been introduced by Feng et al. 30 for complex high refractive index core-shell nanoparticles.
In this work, we propose the metal-dielectric hybrid design for ideal MD scatterer with eliminating electric type scattering due to anapole mode excitation. Therefore, the operating frequency corresponds to the minimum radar cross section (RCS) and, additionally, it demonstrates collateral mantle cloaking effect.
The anapole mode interpretation requires the introduction of the toroidal multipoles -a separate family of multipoles of complex electromagnetic configuration 31,32 . Let us imagine the poloidal currents j flowing along the torus meridians and generating a magnetic field m in metaparticles of toroidal shape. Respectively, the toroidal dipole (TD) moment appears as an oscillation along the torus main axis and manifests a strong field localization even within the point-like source, possessing far-field radiation pattern similar to the electric dipole (ED) mode. Namely, the interference between ED and TD modes under the condition P = ikT, where P stands for ED moment and T is the TD moment, mutually cancels their far-field radiations 29 . This condition, also called anapole mode, leads to nonradiating states intended for development of nonradiating sources, noninvasive sensors, cloaking devices, etc. [33][34][35][36][37][38][39][40][41][42][43][44][45][46][47] .
The multipole decomposition of far-field intensities that includes toroidal multipoles contribution is presented in the work by Vaman and Radescu 48 and are expressed as: Thus, the system characterized by anapole (P = ikT) state accompanied by the presence of magnetic dipole M will scatter as an ideal MD, since only that component is present on the right hand side in the equation above.
Results the metasurface design. For the realization of the ideal MD scatterer concept, here we propose the design of a structure consisting of an infinitely elongated Perfectly Electric Conductor (PEC) cylinder core and a hybrid PEC metasurface coating with a dielectric layer embedded between them. Providing that the PEC core is a strong electromagnetic scatterer, the dielectric coating with imprinted metasurface of sinusoidal-like pattern is assumed to considerably reduce its scattering properties. From electromagnetic point of view, a metasurface can be considered as a mantle cloaking device [49][50][51] , since cloaking can be realized by the use of a thin layer usually composed of periodic arrangement of unit cells in 2D. Thus, we treat this structure based on its scattering properties and estimate the field distribution properties [52][53][54] .
The characteristics of the metasurface can be equivalently defined in terms of surface impedance. This latter modifies the boundary condition of the object connecting the incident tangential electric field with the induced surface currents 51,55 .
Considering a PEC cylinder, the total field, in cylindrical coordinates, can be expressed as an infinite sum of harmonics in the form: where J m and H m (2) are Bessel and Hankel functions that describe the incident and the scattered field, respectively, c m are the scattering coefficients, k b is the wavenumber in the background medium that here is considered as vacuum.
Scattering cancellation can be therefore achieved by the annulment of coefficients c m . This leads to the formu- , being a and b the (external) radius of the metallic cylinder without and with coating, respectively. In the previous expression "prime" denotes the derivative (with respect to the radial coordinate ρ) and γ . In this way, it can be proven that a single value of the surface impedance is directly associated to one harmonic index m and therefore to the annulment of a single scattering coefficient (usually the dominant one, if it does exist).
In this paper we consider a non-electrically small cylinder. In this case, a larger number of harmonics contribute to the scattered field 49 ; therefore, a homogeneous surface impedance is not sufficient to cloak the cylinder. For this reason, in the following we propose a modulated cell profile, in order to provide a non-homogeneous value of surface impedance and therefore to enlarge the cloaking bandwidth. The expected bandwidth enlargement can also be explained by the fact that, similar to biconical structure, the modulated profile incorporates different resonant lengths of closely located frequencies that give rise to a wide-band resonance.
In particular, the structure consists of a PEC cylinder of radius r PEC = a = 20 mm wrapped into a hybrid coating composed of a dielectric substrate of relative permittivity ε r = 3, for an easier practical realization, and an infinitely thin metasurface layer (let us suppose PEC layer) (Fig. 1a). The thickness of the dielectric substrate is t = 2.92 mm and its height is D v = 32 mm.
For convenient field distribution estimation, we introduce the cell profile described by the analytical formulation in 54 : where the parameter α is related to the cell modulation that ensures a variation of the dielectric constant and the offered surface impedance inside the unit cell leading, therefore, to an enlargement of the operation bandwidth [52][53][54] . The unit cell period is determined by the length D u = 4 mm, while the minimum and maximum height of metallic pattern are W min = 0.2 mm and W max = 27.3 mm, respectively. The one-round wrapping contains N = 36 unit cells (Fig. 1b). During the numerical simulations, performed with commercial solver Microwave Studio by CST, periodic boundary conditions along the ±z directions have been set, imitating an infinitely long cylinder. In such a way, no edge diffraction effects are present in the results that could act as noise in the analysis described below.
Variations of the unit cell parameters enable to control the operation frequency of the scatterer. In particular, the working frequency is inversely proportional to W max 54 .

Scattering properties and field distribution.
To study the ideal magnetic dipole scatterer, we consider the scattering properties by computing the object RCS and the field maps of the electric and magnetic field components intensities. Moreover, we estimate the evolution of the chosen parameters and we study the electromagnetic response of scatterers with similar pattern of a different ellipse height W max , namely, W 1 = 27.3 mm, W 2 = 23.4 mm, W 3 = 19.5 mm, W 4 = 17.5 mm and W 5 = 15.6 mm, while maintaining D v and W min constant.
The structure is illuminated by a normally incident plane wave of TM polarization and it shows minimum RCS (~1000 mm 2 ) for W 1 = 27.3 mm at f = 2.96 GHz (Fig. 2a). As for the other patterns of the metasurface, the RCS is higher up to 1200 mm 2 and shifted toward higher frequencies (Fig. 2b-e). Moreover, strong oscillations at higher than the resonance frequencies correspond to a resonance behavior of scattering of the system. Furthermore, we study the electric and magnetic field distribution at the minimum scattering frequency f = 2.96 GHz for cylinder coated with PEC metasurface (Fig. 3c,d) and for bare PEC cylinder at the same frequency (Fig. 3a,b). In the case of the bare cylinder, electric and magnetic field amplitude distributions evidently resemble the scattering from simple electric scatterer. On the other hand, it can be observed that the field properties of the metasurface coated PEC cylinder exhibit a cloaking effect. Indeed, the electric and magnetic field perturbations are strongly reduced, and field lines replicate the plane wave propagation.
Moreover, we compare far-field scattering of the bare and the cloaked cylinder. The results are reported in Fig. 4, showing a reduction of the RCS for all directions around the cylinder. Further details on the cloaking effect of the metasurface such as the dependence on the materials properties or the scalability of the structure are reported in the Supplementary material.
The multipolar description of near-field interaction. The origin of the suppressed RCS and accompanied mantle cloaking effect can be accurately explained in terms of multipole decomposition approach. We have studied the near-field properties of the interaction between incident TM-polarized wave and ideal MD scatterer. In this instance, multipolar decomposition provides an accurate description of the near-field properties due to charge-current redistribution inside the object 48 and interpreted far-field scattering fields. The lowest www.nature.com/scientificreports www.nature.com/scientificreports/ order multipoles, such as dipoles and quadrupoles, are considered as the strongest ones and are related to electric, magnetic and toroidal families. This selection of multipoles correctly characterizes the pseudotorus topological electromagnetic excitations 29,31,47 .
For explanation of the origin of the encountered cloaking effect, we have carried out multipolar decomposition near the minimum scattering value at f = 2.96 GHz (Fig. 2f). The scattering intensity of the strongest multipoles is given in terms of electric P, magnetic M and toroidal T dipoles, electric Qe, magnetic Qm and toroidal Qt www.nature.com/scientificreports www.nature.com/scientificreports/ quadrupoles. The key role in the scattering characterization is played by the magnetic dipole mode M which dominates all other standard multipoles (Fig. 2f). In the proximity of the resonance frequency, the power scattered by the magnetic dipole moment is more than 2 times higher than that of the electric and toroidal dipole moments and corresponding quadrupoles. However, this pronounced scattering resonance is due to the suppression of the far-field scattering defined by electric-type of multipoles -we remind here that electric and toroidal radiations exhibit identical radiation patterns. More significantly, minimal scattering at f = 2.96 GHz is accompanied by the overlapping of electric P and toroidal T dipole moments resulting in anapole mode excitation, more appreciated  www.nature.com/scientificreports www.nature.com/scientificreports/ for total elimination of the electric scattering due to the appearance of the P = ikT condition. That said, the scattering origin is interpreted as a magnetic response of higher magnetic dipole intensity associated with the fully eliminated electric component of scattering.
We also note that the initial response of PEC cylinder is due to electric dipole mode with the suppression of other multipoles contribution, including magnetic response. On the other hand, the metasurface plays the role of cloak that could be described as a compensating element, having influence only on electric types of multipoles. Thus, the cloaking technique accompanied by a RCS minimum is also determined by a magnetic scattering response, hence the cancelation of electric scattering coefficient in Eq. (1). In this view, we consider the different parameters of the metasurface. Even if the magnetic dipole moment has the maximum intensities for W max = 23.4 mm and W max = 19.5 mm, anapole mode cannot be established due to unfulfillment of the P = ikT relation, i.e., the toroidal mode is shifted in frequency from the resonance for these cases. Thus, the RCS achieves 1193 mm 2 and 1233 mm 2 for W max = 23.4 mm and W max = 19.5 mm, respectively. However, the toroidal quadrupole Qt plays the key role at resonance appearing in cases of W max = 17.5 mm and W max = 15.6 mm.
It should be stressed that the RCS damping by anapole mode is crucial to realize a magnetic dipole scattering by the cylindrical particle for a TM-polarized wave excitation. Our results demonstrate that the proposed nonmagnetic metaparticle can be used for magnetization without complex systems providing currents circulation as in SRR and hybrid nanoparticles. Therefore, we have demonstrated a novel approach for realization of a mantle cloaking of simplified design for the generation of magnetic light due to ideal magnetic dipole scattering, previously proposed only for hybrid metal-dielectric nanoparticles.

Discussion
Summarizing, we have presented a cloaking technique appreciated for the number of applications from widely known camouflage and stealth techniques for military purposes to modern medical and biological probes. The interest in cloaking techniques arises due to strong reduction of fields scattered from the cloaked object in all directions around the scatterer position. This leads to undisturbed field interaction between an incident light so that the cloaked object becomes invisible for external observers.
Modern literature reports on various cloaking techniques available due to their simple practical implementation. Apart from fundamentally important transformation optics (TO) technique, the plasmonic cloaking (PC) and mantle cloaking (MC) [49][50][51][60][61][62]65,66 are more common from a practical point of view. In literature, different approaches to control the scattering of light are found [67][68][69][70][71][72][73][74][75][76][77] . In particular, MC differs from others since it does not require bulk anisotropic materials and is realized by using a thin layer metasurface usually consisting of a periodic arrangement of unit cells. Moreover, we considered the cloaked device in terms of magnetic response, discussed in previous literature to enhance the Purcell effect 78 , and anapole mode 79,80 . In our analysis, the cloaking effect is explained by the suppression of scattering multipoles of the same polarization as of the incident wave (electric and toroidal dipole moment). Additionally, magnetic dipole-like response has been demonstrated that is due to the excess part of the anapole effect of electric type.

conclusion
In this paper, we considered a new type of ideal magnetic metal-dielectric hybrid scatterer based on well-pronounced magnetic dipole moment with simultaneously suppressed electric response leading to minimization of total scattering. Therefore, this strong scattering suppression has a similar consequence as mantle cloaking effect; at the specific resonance frequency, one could describe the cloaking effect due to multipoles interplay. Moreover, we proposed a properly design metasurface for practical realization of ideal magnetic dipole effect. Its performances have been numerically proved.

Methods
Simulations. The electromagnetic properties of metal-dielectric hybrid design is computed by commercial Maxwell's equation solver, namely, CST Microwave Studio, using the standard transient modeling approach. The simulations provide values of RCS, giving data on scattering properties of the structure. Furthermore, the field maps of electric and magnetic field distributions in near-field of the structure are obtained. The current densities induced in the structure are used to calculate the powers radiated by conventional multipoles, including those of toroidal dipoles.