###### Abstract

Copyright © 2015 Wei Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

###### Abstract

For homogenous isotropic dielectric nanospheres with incident plane waves, Cartesian electric and toroidal dipoles can be tunned to cancel each other in terms of far-field scattering, leading to the effective anopole excitation. At the same time however, other multipoles such as magnetic dipoles with comparable scattered power are simultanesouly excited, mixing with the anopole and leading to a non-negligible total scattering cross section. Here we show that for homogenous dielectric nanospheres, radial anisotropy can be employed to significantly suppress the other multipole excitation, which at the same time does not compromise the property of complete scattering cancallation between Cartesian electric and toroidal dipoles. This enables an elusive pure anopole excitation within radially anisotropic dielectric nanospheres, which may shed new light to many scattering related fundamental researches and applications.

Elusive pure anapole excitation in homogenous spherical nanoparticles with radial anisotropy

Wei Liu, Bing Lei, Jianhua Shi, Haojun Hu, and Andrey E. Miroshnichenko

College of Optoelectronic Science and Engineering, National University of Defense Technology, Changsha, Hunan 410073, P. R. China.

Nonlinear Physics Centre, Research School of Physics and Engineering, Australian National University, Acton, ACT 0200, Australia.

Corresponding author: wei.liu.pku@gmail.com

## 1 Introduction

Since the discovery by J. J. Thomson that electrons are constituents of atoms [1], there had been a lot of searches for charge-current distributions in which no electromagnetic energy are radiated, which can be applied to explain the stability of atoms [2]. This is due to the existing dogma since the establishment of Maxwell equations that accelerating charges would lose energy continuously through electromagnetic radiation, and consequently to prevent the collapse of atoms, the charges within them must be arranged into specific configurations. Such special charge-current configurations are also termed as non-radiating sources [2, 3]. Those non-radiating sources are deeply related to the concepts of anapoles in quantum field theory, which have their roots in the law of parity unconservation [4], and can be used as possible explanations for the existence and features of neutrinos and cold dark matter [5, 6].

For time-varying oscillating charge-current distributions, the dynamic version of anapoles is highly related to the so called toroidal multipoles, which recently have attracted enormous attention [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20]. As independent terms in the multipole expansion of electromagnetic fields besides conventional cartesian magnetic and electric multipoles [3], toroidal multipoles correspond to significantly different near-field charge-current distributions and will inevitably become significant when the operating spectral regime goes beyond the electrostatic region [21, 20, 22]. In the far field however, the toroidal multipoles have identical scattering patterns compared with their conventional Cartesian electric and magnetic counterparts, and that is to say, they are indistinguishable in terms of far-field scattering [21, 7, 23, 13]. This at the same time also means that when properly engineered, the contributions of conventional and toroidal multipoles to the scattered field can be tuned to the same amplitude and out of phase, which leads to the far-field scattering cancellation and thus generate the so called electrodynamic anapoles [24, 21, 11, 15, 16, 17, 19]. The simplest case of such anapoles is induced by the destructive cancellation of Cartesian electric dipole and toroidal dipole, which has been demonstrated firstly in composite carefully engineered structures consisting of particle clusters [11, 15]. Recently it has also been theoretically proposed and experimentally demonstrated in surprise that the anapole can be effectively excited even in individual homogenous dielectric nanoparticles [16]. However, the problem of the above demonstration in simple structures is that besides the anapole excitation, other multipoles (such as magnetic dipoles) of similar magnitudes in terms of scattered power are also excited and mixed with the excited electric and toroidal dipoles [16, 17]. Such problems can be conquered within single core-shell nanoparticles of fundamental shapes (spherical or cylindrical) where the extra multipole excitations can be significantly suppressed, enabling elusive pure anapole excitations [17, 19].

Here we propose an alternative approach and demonstrate that elusive pure anapole can be observed in homogenous radially anisotropic dielectric spherical nanoparticles under incident plane waves. It is shown that when the radial anisotropy of material permittivity is employed, the Cartesian electric and toroidal dipoles can be tuned to the spectral positions where other multipoles are significantly suppressed, without compromising the feature of complete destructive interference between them. Since Cartesian electric and toroidal dipoles can totally cancel each other in terms of far-field scattering and at the same time the excitation of other multipoles is negligible, it is proved that effectively pure anapole has been obtained. The anapole excitation is highly related to the exotic scattering phenomena of invisibility and possible applications of non-invasive sensing and detections, which may pave the way for many further scattering related fundamental researches and applications.

## 2 Methods and Expressions

### A Multipole Expansions with toroidal components

For any charge-current distribution, when the geometric size of the source is far smaller than the effective wavelength (the retardation effect within the source can be neglected), the quasi-static approximation can be applied and then the radiated fields of the source can be expanded into two sets of conventional Cartesian electric and magnetic multipoles [3, 21]. The lowest order terms of the two sets of multipoles are electric and magnetic dipoles, which can be expressed respectively as [throughout the paper we adopt the notation for electromagnetic waves]:

(1) |

where is the speed of light, is the angular frequency and J is the current. When the retardation effect within the source cannot be neglected (the size of the source gets larger and/or the effective wavelength gets smaller), the two sets of multipoles are not complete and thus other multipoles have to be included as correcting terms [21, 20, 22]. The lowest order correcting multipoles are toroidal multipoles, with the toroidal dipole expressed as as [24, 21, 7, 23]:

(2) |

The radiated power of P, T and M are respectively:

(3) | |||

(4) | |||

(5) |

where is the vacuum permeability and is the angular wave number. We emphasize that the expressions above are for individual multipoles only and a simple sum of them does not lead to the total scattered power as there is interference between them [21, 7, 23] [see also Eq. (10) below]. It is also worth mentioning that the incorporation of toroidal multipoles does not make the expansion complete [21, 20, 22]. When the operating spectral region is far beyond the quasi-static regime, besides toroidal multipoles, other higher order correcting multipoles have to be also included [21, 20, 22].

For scattering particles, based on the charge conservation relation, the current (displacement current) distribution is

(6) |

where is the vacuum permittivity; is the relative permittivity, which is a scalar for isotropic particles or tensor for anisotropic particles. With the current distribution shown in Eq.(6), all the Cartesian multipoles including P, T and M and their scattered power can be directly calculated [see Eqs.(1)-(5)].

### B Electric and magnetic multipoles deduced from vector spherical harmonics of Mie scattering particles

For spherical (both isotropic and radially anisotropic) particles with incident plane waves of electric field and magnetic field , both the fields inside the particle and the scattered fields can be expressed in transverse electric and transverse magnetic (in terms of electric field distribution)vector spherical harmonics [25, 25, 26, 27, 28, 20]. On one hand, according to Eqs.(1)-(6), all the Cartesian multipoles including P, T and M and their scattered power can be analytically expressed and numerically calculated. On the other hand, based on the far-field equivalence, the scattered fields can be viewed as the radiated fields of a combination of conventional Cartesian electric and magnetic multipoles (which correspond to transverse magnetic and transverse electric vector spherical harmonics respectively). As a result the far-field deduced electric dipole and magnetic dipole can be expressed respectively as [25, 29, 16, 30, 31, 17, 32, 28, 20]:

(7) |

where and are the Mie scattering coefficients of the lowest order. The total scattered power is [25, 26, 27, 28]:

(8) |

where and are the Mie scattering coefficients of order and . It is clear that the scattered power of and is respectively:

(9) |

### C Anapole excitation induced by destructive interference between P and T

As is discussed above, shown in Eq.( 7) has been obtained based on the complete vector spherical harmonic expansion and thus has included all the transverse magnetic dipolar scattering components [29, 16, 30, 31, 17, 32, 28, 20]. While the conventional Cartesian electric dipole P shown in Eq. (1) is calculated through current integration and is only the lowest order expansion of [21, 20, 22]. For operating spectral region beyond the quasi-static regime, the next order correcting term, which is the toroidal dipole T shown in Eq.( 7), should also be included. When higher order correcting terms beyond toroidal terms can be neglected, the scattered dipolar transverse magnetic field can be expressed as as [21, 23, 16, 17]:

(10) |

where n is the unit vector along r. The existence of the anapole lies in which can be satisfied when:

(11) |

It is worth mentioning that when Eq.( 11) is satisfied, only the transverse magnetic dipolar scattering will be vanishing and consequently the anapole excited (consisting of destructively interfering P and T) can be mixed with other excited multipoles, such as magnetic dipoles, electric quadrupoles and so on [16, 17, 19].

## 3 Results and Discussions

In Fig. 1 we provide a schematic illustration of anopole excitation in a homogenous particle with plane wave incidence. The plane wave is polarized along direction (in terms of electric field) and propagating along direction. The radius of the particle is R, and the particle is isotropic with refractive index , or is radially anisotropic with transverse (perpendicular to the radial direction) permittivity and radial permittivity , with the anisotropy parameter defined as: . Both T and P are excited, out of phase and of the same magnitude with regard to the scattered power, indicating excitation of the anapole [Eqs. (10)-(11)]. Also the typical current J and magnetic field H distributions for T only have been shown.

### A Mixed excitation of anapole and other multipoles within homogeneous isotropic dielectric nanospheres

At the beginning, we investigate the scattering of homogeneous isotropic () dielectric nanospheres with plane wave incidence with nm and [see the inset of Fig. 2(a)]. The scattered power spectra are in Fig. 2(a): the contributions from far field deduced electric and magnetic dipoles and , conventional Cartesian electric and toroidal dipoles P and T, and the total scattered power are shown. The shaded region in Fig. 2(a) is the area of interest, and thus is zoomed in and shown in Fig. 2(b), where two points A and B are pinpointed and indicated by black dots. At point A of nm, for the transverse magnetic dipolar components, there is only toroidal dipole T excitation [], and the conventional Cartesian electric dipole P is completely suppressed; At point B of nm, there is equal out of phase T and P excitation in terms of scattering power, leading to an effective anapole excitation defined in Eqs. (10)-(11) [16, 17, 19, 20]. The corresponding near-field distributions are shown in Fig. 2(A-B): the vector-plots correspond to the displacement field on the plane and the color-plots correspond to the normalized magnetic fields along direction (). We note here that for better visibility we plot only transverse magnetic dipolar fields and as is shown: (1) At point A shown in Fig. 2(A), the near field exhibits a typical toroidal dipolar distribution, as is schematically shown in Fig. 1 (the magnetic field is confined within the particle and there are accompanying oppositely circulating current loops); (2) Point B shown Fig. 2(B) corresponds to an anapole excitation, and at this point the fields outside the particle are vanishing, confirming the negligible overall dipolar scattering and thus effective anapole excitation.

According to Fig. 2(b), at Point B besides the anapole excitation, there is the coexistence of other multipoles, of which the far field deduced magnetic dipole is the most significant one. At the same time, the scattered power of is close to that of both P and T, rendering the anapole excitation of a hybrid and mixed nature. Basically to obtain the sole excitation of a pure anapole, the suppression of is required.

### B Overlapping the minimum scattering points of and within radially anisotropic nanospheres

It is recently demonstrated that radial anisotropy can be employed to tune spectrally the far field deduced electric dipoles and thus the anapoles excited [28, 20]. In contrast when the transverse permittivity is fixed, the spectra of far field deduced magnetic dipoles are independent of the anisotropy parameter due to its transverse electric nature [28]. In Fig. 3(a) we show the scattering spectra of radially anisotropic dielectric nanosphere (transverse permittivity is fixed at and nm). Only the contributions from and are shown. It is clear that would blueshift/redshift with larger/smaller , while is stable and independent of . We indicate the first (counted from large wavelength) minimum scattering points of and by circular dots (which corresponds to the anapole excitation [20]) and respectively in Fig. 3(a), and their spectra dependence on is shown Fig. 3(b). According to Fig. 3(b), besides the confirmation of the spectral tunability of and stability of , it is further shown that at the minimum scattering points of both and overlap. This means the anapole excitation is achieved with simultaneous magnetic dipole suppression.

### C Pure anapole excitation within radially anisotropic nanospheres

To further show the significance of the overlapping of the minimum scattering points of and and the sole anapole excitation, we show in Fig. 4(a) the spectra at for the radially anisotropic nanosphere with and nm [see the inset of Fig. 4(a)]. Similar to Fig. 2(a), contributions from , , P, T together with the total scattered power are shown. Two points C ( nm) and D ( nm) are indicated by black dots [the corresponding near fields are shown in Fig. 4(C)-Fig. 4(D) respectively], whose counterparts (which are similar in terms of both scattering spectra and near-field distributions) are respectively points A and B in Fig. 2(a): at point C there is dominant T excitation [with magnetic field confined within the particle and oppositely circulating currents, see Fig. 4(C)] and point D corresponds to the excited anapole [with almost zero fields everywhere outside the particle, see Fig. 4(D)]. By comparing point B [see Fig. 2] and D [see Fig. 4], the significance of introducing radial anisotropy is clear: (1) The anapole is excited at the minimum scattering point of and at the same time besides the significant suppression of , other multipoles including the electric and magnetic quadrupoles are also negligible (not shown here but is clear according to the total scattering), leading to vanishing total scattering; (2) Though completely cancel each other, the scattering of P and T is an order of magnitude larger than all the other multipoles, showing a very efficient anapole excitation. The two points above are sufficient to justify the claim of elusive pure anapole excitation within radially anisotropic nanospheres.

## 4 Conclusions and Outlook

To conclude, in this work we show the elusive excitation of pure anapoles in homogenous nanoparticles, which has been realized through employing radial anisotropy of the material permittivity to tune the anapole excited to overlap spectrally with the minimum scattering position of magnetic dipoles. At the overlapping position, Cartesian electric and toroidal dipoles can totally cancel the scattering of each other in the far field and at the same time the excitation of other multipoles (including but not limited to magnetic dipoles) is negligible, proving that the pure anapole has been achieved.

We emphasize that the anapole excited has been induced by the interferences between Cartesian electric and toroidal dipoles, which is featured by null scattering in terms of only transverse magnetic dipolar scattering. In a similar way, anapoles characterized by null transverse electric dipolar scattering, or both transverse electric and magnetic null scattering of higher orders can be achieved. The investigations into scattering nanoparticles from the new perspective of anapole and toroidal multipole excitations can also be conducted for other non-spherical structures where numerical methods can be applied to calculate the fields, currents, and thus multipoles excited. Also in this work we have confined our study to correcting terms up to toroidal multipoles [24, 21, 33, 20]. When the operating spectral region is far beyond the quasi-static regime, higher order correcting terms than toroidal multipoles (in terms of ) will arise inevitably [24, 21, 33, 20], and as a result it is expected that anapoles induced by the inferencing between more than two multipoles can be obtained.

Here in this work we show only the demonstration and tuning of anapoles based on radial anisotropy, and other effects, such as nonlinear and gain effects, can also be applied for both near and far-field scattering manipulations based on anapole control. The investigations conducted into the fundamental Mie scattering problem from the new research angle of anapole excitation and tuning is of great significance for not only deepening our understanding of the basic scattering physics, but also for quite a few scattering related applications, such as far-field scattering pattern shaping [34, 30, 35, 28], non-invasive detection and sensing [36] and so on. The anapole mode exhibits a number of nontrivial topological properties, which may render new opportunities for manipulating electromagnetic waves in an unprecedented way when combined with the emerging fields of topological and low dimensional photonics. Moreover, the investigations into anapoles should not be confined to the field of electrodynamics, and similar studies can certainly be conducted in other fields such as acoustics, fluid dynamics, quantum field theory and so on.

## Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

## Acknowledgments

We thank Jianfa Zhang and Yuri S. Kivshar for valuable discussions, and acknowledge the financial support from the National Natural Science Foundation of China (Grant numbers: and ), the Australian Research Council (Grant number: FT110100037) and the Basic Research Scheme of College of Optoelectronic Science and Engineering, National University of Defense Technology.

## References

- [1] J. Z. Buchwald and A. Warwick, Histories of the electron: the birth of microphysics (MIT Press, 2004).
- [2] G. Gbur, “Nonradiating sources and other invisible objects,” Prog. Optics 45, 273 (2003).
- [3] J. D. Jackson, Classical electrodynamics (Wiley New York, 1962).
- [4] I. B. Zel’Dovich, “Electromagnetic interaction with parity violation,” Sov. Phys. JETP 6, 1184 (1958).
- [5] V. M. Dubovik and V. E. Kuznetsov, “The toroid dipole moment of the neutrino,” Int. J. Mod. Phys. A 13, 5257 (1998).
- [6] C. M. Ho and R. J. Scherrer, “Anapole dark matter,” Phys. Lett. B 722, 341 (2013).
- [7] T. Kaelberer, V. Fedotov, N. Papasimakis, D. Tsai, and N. Zheludev, “Toroidal dipolar response in a metamaterial,” Science 330, 1510 (2010).
- [8] Z. Dong, P. Ni, J. Zhu, X. Yin, and X. Zhang, “Toroidal dipole response in a multifold double-ring metamaterial,” Opt. Express 20, 13065 (2012).
- [9] Y.-W. Huang, W. T. Chen, P. C. Wu, V. Fedotov, V. Savinov, Y. Z. Ho, Y.-F. Chau, N. I. Zheludev, and D. P. Tsai, “Design of plasmonic toroidal metamaterials at optical frequencies,” Opt. Express 20, 1760 (2012).
- [10] B. Ogut, N. Talebi, R. Vogelgesang, W. Sigle, and P. A. van Aken, “Toroidal plasmonic eigenmodes in oligomer nanocavities for the visible,” Nano Lett. 12, 5239 (2012).
- [11] V. A. Fedotov, A. Rogacheva, V. Savinov, D. Tsai, and N. I. Zheludev, “Resonant transparency and non-trivial non-radiating excitations in toroidal metamaterials,” Sci. Rep. 3, 2967 (2013).
- [12] Y. Fan, Z. Wei, H. Li, H. Chen, and C. M. Soukoulis, “Low-loss and high-q planar metamaterial with toroidal moment,” Phys. Rev. B 87, 115417 (2013).
- [13] V. Savinov, V. Fedotov, and N. Zheludev, “Toroidal dipolar excitation and macroscopic electromagnetic properties of metamaterials,” Phys. Rev. B 89, 205112 (2014).
- [14] Q. Zhang, J. J. Xiao, X. M. Zhang, D. Han, and L. Gao, “Core-shell-structured dielectric-metal circular nanodisk antenna: Gap plasmon assisted magnetic toroid-like cavity modes,” ACS Photon. 2, 60 (2014).
- [15] A. A. Basharin, M. Kafesaki, E. N. Economou, C. M. Soukoulis, V. A. Fedotov, V. Savinov, and N. I. Zheludev, “Dielectric metamaterials with toroidal dipolar response,” Phys. Rev. X 5, 011036 (2015).
- [16] A. E. Miroshnichenko, A. B. Evlyukhin, Y. F. Yu, R. M. Bakker, A. Chipouline, A. I. Kuznetsov, B. Lukyanchuk, B. N. Chichkov, and Y. S. Kivshar, “Seeing the unseen: observation of an anapole with dielectric nanoparticles,” Nat. Commun. 6, 8069 (2015).
- [17] W. Liu, J. Zhang, and A. E. Miroshnichenko, “Toroidal dipole induced transparency for core-shell nanoparticles,” Laser Photon. Rev.,doi:10.1002/lpor.201500102 (2015).
- [18] R. M. Bakker, D. Permyakov, Y. F. Yu, D. Markovich, R. Paniagua-Dom¨ªnguez, L. Gonzaga, A. Samusev, Y. Kivshar, B. Luk¡¯yanchuk, and A. I. Kuznetsov, “Magnetic and electric hotspots with silicon nanodimers,” Nano Lett. 15, 2137 (2015).
- [19] W. Liu, J. Zhang, B. Lei, H. Hu, and A. E. Miroshnichenko, “Invisible nanowires with interfering electric and toroidal dipoles,” Opt. Lett. 40, 2293 (2015).
- [20] W. Liu, J. Shi, B. Lei, H. Hu, and A. E. Miroshnichenko, “Efficient excitation and tuning of toroidal dipoles within individual homogenous nanoparticles,” arXiv:1508.02520 (2015).
- [21] E. E. Radescu and G. Vaman, “Exact calculation of the angular momentum loss, recoil force, and radiation intensity for an arbitrary source in terms of electric, magnetic, and toroid multipoles,” Phys. Rev. E 65, 046609 (2002).
- [22] I. Fernandez-Corbaton, S. Nanz, R. Alaee, and C. Rockstuhl, “Exact dipolar moments of a localized electric current distribution,” arXiv:1507.00752 (2015).
- [23] J. Chen, J. Ng, Z. Lin, and C. Chan, “Optical pulling force,” Nat. Photon. 5, 531 (2011).
- [24] V. M. Dubovik and V. V. Tugushev, “Toroidmoments in electrodynamics and solid-state physics,” Phys. Rep. 187, 145 (1990).
- [25] C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983).
- [26] C. W. Qiu and B. Luk’yanchuk, “Peculiarities in light scattering by spherical particles with radial anisotropy,” J. Opt. Soc. Am. A 25, 1623 (2008).
- [27] C. Qiu, L. Gao, J. D. Joannopoulos, and M. Soljačić, “Light scattering from anisotropic particles: propagation, localization, and nonlinearity,” Laser Photonics Rev. 4, 268 (2010).
- [28] W. Liu, “Ultra-directional super-scattering of homogenous spherical particles with radial anisotropy,” Opt. Express 23, 14734 (2015).
- [29] W. T. Doyle, “Optical properties of a suspension of metal spheres,” Phys. Rev. B 39, 9852 (1989).
- [30] W. Liu, A. E. Miroshnichenko, D. N. Neshev, and Y. S. Kivshar, “Broadband unidirectional scattering by magneto-electric core-shell nanoparticles,” ACS Nano 6, 5489 (2012).
- [31] W. Liu, A. E. Miroshnichenko, D. N. Neshev, and Y. S. Kivshar, ‘‘Polarization-independent fano resonances in arrays of core-shell nanoparticles,” Phys. Rev. B 86, 081407 (2012).
- [32] M. S. Wheeler, J. S. Aitchison, and M. Mojahedi, “Coated nonmagnetic spheres with a negative index of refraction at infrared frequencies,” Phys. Rev. B 73, 045105 (2006).
- [33] I. Fernandez-Corbaton, S. Nanz, and C. Rockstuhl, “On the dynamic toroidal multipoles,” arXiv:1507.00755 (2015).
- [34] W. Liu, A. E. Miroshnichenko, and Y. S. Kivshar, “Control of light scattering by nanoparticles with optically-induced magnetic responses,” Chin. Phys. B 23, 047806 (2014).
- [35] W. Liu, J. Zhang, B. Lei, H. Ma, W. Xie, and H. Hu, “Ultra-directional forward scattering by individual core-shell nanoparticles,” Opt. Express 22, 16178 (2014).
- [36] A. Alu and N. Engheta, “Cloaking a sensor,” Phys. Rev. Lett. 102, 233901 (2009).