###### Abstract

We revisit the fundamental topic of light scattering by single homogenous nanoparticles from the new perspective of excitation and manipulation of toroidal dipoles. It is revealed that besides within all-dielectric particles, toroidal dipoles can also be efficiently excited within homogenous metallic nanoparticles. Moreover, we show that those toroidal dipoles excited can be spectrally tuned through adjusting the radial anisotropy parameters of the materials, which paves the way for further more flexible manipulations of the toroidal responses within photonic systems. The study into toroidal multipole excitation and tuning within nanoparticles deepens our understanding of the seminal problem of light scattering, and may incubate many scattering related fundamental researches and applications.

Efficient excitation and tuning of toroidal dipoles within individual homogenous nanoparticles

Wei Liu, Jianhua Shi, Bing Lei, 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, Centre for Ultrahigh-bandwidth Devices for Optical Systems (CUDOS), Research School of Physics and Engineering, Australian National University, Acton, ACT 0200, Australia.

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

OCIS codes: (290.5850) Scattering, particles; (290.4020) Mie theory; (160.1190) Anisotropic optical materials

## References and links

- [1] T. Kaelberer, V. Fedotov, N. Papasimakis, D. Tsai, and N. Zheludev, “Toroidal dipolar response in a metamaterial,” Science 330, 1510 (2010).
- [2] 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).
- [3] 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).
- [4] 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).
- [5] 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).
- [6] 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).
- [7] V. Savinov, V. Fedotov, and N. Zheludev, “Toroidal dipolar excitation and macroscopic electromagnetic properties of metamaterials,” Phys. Rev. B 89, 205112 (2014).
- [8] 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).
- [9] 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).
- [10] 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. (2015). arXiv:1412.0299.
- [11] 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). arXiv:1412.4931.
- [12] 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).
- [13] 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).
- [14] V. M. Dubovik and V. V. Tugushev, “Toroidmoments in electrodynamics and solid-state physics,” Phys. Rep. 187, 145 (1990).
- [15] 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).
- [16] J. D. Jackson, Classical electrodynamics (Wiley New York, 1962).
- [17] C. F. Bohren and D. R. Huffman, Absorption and Scattering of Light by Small Particles (Wiley, 1983).
- [18] L. Novotny and B. Hecht, Principles of nano-optics (Cambridge University Press, Cambridge, 2012), 2nd ed.
- [19] W. Liu, R. F. Oulton, and Y. S. Kivshar, “Geometric interpretations for resonances of plasmonic nanoparticles,” Sci. Rep. 5, 12148 (2015).
- [20] A. I. Kuznetsov, A. E. Miroshnichenko, Y. H. Fu, J. B. Zhang, and B. S. Lukyanchuk, “Magnetic light,” Sci. Rep. 2, 492 (2012).
- [21] A. B. Evlyukhin, S. M. Novikov, U. Zywietz, R. L. Eriksen, C. Reinhardt, S. I. Bozhevolnyi, and B. N. Chichkov, “Demonstration of magnetic dipole resonances of dielectric nanospheres in the visible region,” Nano Lett. 12, 3749 (2012).
- [22] 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).
- [23] W. T. Doyle, “Optical properties of a suspension of metal spheres,” Phys. Rev. B 39, 9852 (1989).
- [24] 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).
- [25] I. Fernandez-Corbaton, S. Nanz, and C. Rockstuhl, “On the dynamic toroidal multipoles,” arXiv:1507.00755 (2015).
- [26] I. Fernandez-Corbaton, S. Nanz, R. Alaee, and C. Rockstuhl, “Exact dipolar moments of a localized electric current distribution,” arXiv:1507.00752 (2015).
- [27] C. C. Lam, P. T. Leung, and K. Young, “Explicit asymptotic formulas for the positions, widths, and strengths of resonances in mie scattering,” J. Opt. Soc. Am. B 9, 1585 (1992).
- [28] W. Liu, “Ultra-directional super-scattering of homogenous spherical particles with radial anisotropy,” Opt. Express 23, 14734 (2015).
- [29] P. B. Johnson and R. W. Christy, “Optical constants of the noble metals,” Phys. Rev. B 6, 4370 (1972).
- [30] C.-W. Qiu, L.-W. Li, T.-S. Yeo, and S. Zouhdi, “Scattering by rotationally symmetric anisotropic spheres: Potential formulation and parametric studies,” Phys. Rev. E 75, 026609 (2007).
- [31] 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).
- [32] Y. X. Ni, L. Gao, A. E. Miroshnichenko, and C. W. Qiu, “Controlling light scattering and polarization by spherical particles with radial anisotropy,” Opt. Express 21, 8091 (2013).
- [33] 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).
- [34] 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).
- [35] R. F. Oulton, V. J. Sorger, T. Zentgraf, R. M. Ma, C. Gladden, L. Dai, G. Bartal, and X. Zhang, “Plasmon lasers at deep subwavelength scale,” Nature 461, 629 (2009).
- [36] M. A. Noginov, G. Zhu, A. M. Belgrave, R. Bakker, V. M. Shalaev, E. E. Narimanov, S. Stout, E. Herz, T. Suteewong, and U. Wiesner, “Demonstration of a spaser-based nanolaser,” Nature 460, 1110 (2009).
- [37] A. V. Kabashin, P. Evans, S. Pastkovsky, W. Hendren, G. A. Wurtz, R. Atkinson, R. Pollard, V. A. Podolskiy, and A. V. Zayats, “Plasmonic nanorod metamaterials for biosensing,” Nat. Mater. 8, 867 (2009).
- [38] L. Novotny and N. van Hulst, “Antennas for light,” Nat. Photon. 5, 83 (2011).
- [39] H. A. Atwater and A. Polman, “Plasmonics for improved photovoltaic devices,” Nat. Mater. 9, 865 (2010).
- [40] L. Lu, J. D. Joannopoulos, and M. Soljacic, “Topological photonics,” Nat. Photon. 8, 821 (2014).
- [41] F. Xia, H. Wang, D. Xiao, M. Dubey, and A. Ramasubramaniam, “Two-dimensional material nanophotonics,” Nat. Photon. 8, 899 (2014).
- [42] G. N. Afanasiev and Y. P. Stepanovsky, “The electromagnetic field of elementary time-dependent toroidal sources,” J. Phys. A Math. Gen. 28, 4565 (1995).
- [43] E. A. Marengo and R. W. Ziolkowski, “Nonradiating sources, the aharonov-bohm effect, and the question of measurability of electromagnetic potentials,” Radio Sci. 37, 19 (2002).

## Appendix A Introduction

Stimulated by the recent experimental demonstration of toroidal dipoles (TDs) excited on the platform of metamaterials consisting of specially organized split ring resonators [1], studies into TDs and more generally toroidal multipoles (TMs) within photonic nano-systems have attached surging attention [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13]. Compared to conventional electric and magnetic multipoles, TMs can be viewed as poloidal currents flowing on the surface of a torus along its meridians, or as a series of magnetic dipoles (enclosed circulating currents) aligned along enclosed paths [14, 15]. For the near field in terms of charge-current distributions, TMs are contrastingly different from their electric and magnetic counterparts [5, 7]. In the far field however, TMs and conventional electric multipoles are indistinguishable [5, 7], which results in the situation that in various scattering theories concerning particle scattering, TMs are not considered as separate entities [16, 17]. Generally under most circumstances, the contributions of TMs have been incorrectly attributed to conventional electric multipoles, rendering TMs unfortunately as missing or hidden terms.

Following the first experimental demonstration of TDs [1], TMs have been observed in various configurations, most of which (influenced by the configuration designed in Ref. [1]) are kind of toroid-shaped composite structures [2, 3, 4, 5, 6, 9]. This is understandable: the corresponding charge-current distribution of TMs is in a toroidal shape and naturally in such composite structures TMs can be efficiently excited. But this unfortunately leads to more or less a dogma that TMs can only be observed in carefully engineered compound structures and consequently there have been very few searches for the existence of TMs in simple individual structures [8]. Only recently it is shown decisively that TMs can be excited within single scattering particles of fundamental (spherical or cylindrical) shapes, which also have been proven to play indispensable roles for some exotic scattering properties [10, 11, 13]. The incorporation of TMs extends significantly our understanding of the fundamental problem of light scattering by spherical and cylindrical particles. Moreover, considering that the scattering of spherical and cylindrical particles is the cornerstone of many scattering related fundamental studies and applications, the incorporation of TMs into scattering photonic systems can possibly spur further breakthroughs in sensing, nanoantennas, nanoscale lasing and photovoltaic devices [18].

Here in this article we restudy the seminal problem of plane wave light scattering by homogenous spherical nanoparticles, but from a new angle of excitation and manipulations of TDs. It is demonstrated that, similar to all-dielectric particles, even a homogenous metallic spherical nanoparticle can surprisingly support TDs in the spectral region beyond the quasi-static regime. We further show that radial anisotropy can be employed to tune spectrally the TDs excited, providing extra freedom for TMs manipulations, which can possibly find significant applications in TMs related topics, especially those where near-field distributions of fields and currents play a crucial role.

## Appendix B Methods and Expressions

### b.1 Conventional Multipole Expansions and Emerging Toroidal Multipoles

To obtain the radiated fields of an arbitrarily distributed charge-current source, the general approach is to firstly calculate the retarded vector potentials and then get the associated electromagnetic fields from those potentials [16]. When the retardation effect within the charge-current source can be neglected (in other words the geometrical dimension of the source is far smaller than the effective wavelength and thus the quasi-static approximation can be applied), the vector potentials can be Taylor-expanded into multiple terms of different orders, from which the conventional electric and magnetic multipoles can be deduced [16]. For example, the zeroth order expansion term corresponds to the conventional electric dipole P [the notation has been adopted where k and denote wave-vector and angular frequency respectively]:

(1) |

and the first order expansion term includes both conventional electric quadrupole and magnetic dipole M:

(2) |

where is the speed of light. Since M corresponds to the higher order term of the Taylor expansion, it results in more or less the mind-sets that: (1) P should be stronger than the M and (2) like most higher order modes [19], M should be supported at higher frequencies than P. Nevertheless, we should keep in mind that the Taylor-expansion of the vector potentials [16] is based on the quasi-static approximation, and thus when the precondition of this approximation does not exist, the above conclusion is not valid anymore. An outstanding example of this is the recent demonstration of optically-induced magnetic response of higher permittivity dielectric particles [20, 21, 22], where M is supported at lower frequencies and can be stronger than the electric dipole.

The other consequence of the breakdown of the quasi-static approximation is that the conventional electric and magnetic dipoles [Eqs. (1)-(2)] and higher order multipoles do not form a complete set and thus other multipoles have to be included as correcting terms [14, 15, 7]. For dipolar responses, besides M and P, the lowest order correcting term would be the toroidal dipole:

(3) |

The scattered power of P and T are respectively:

(4) |

where is the vacuum permeability.

### b.2 Electric and magnetic multipoles associated with spherical harmonics of Mie scattering particles

For spherical (both homogenous and multi-layered) particles with incident plane waves of electric field , the scattered electric field can be expressed as [17]£º

(5) |

where and are the Mie scattering coefficients; and are vector spherical harmonics; and . In the far field the scattered fields [Eq.(5) for ] associated with [] and [] are exactly the same as the far fields of conventional electric and magnetic multipoles of order , with conventional electric and magnetic dipoles P and M () defined in Eqs.(1)-(2). For this far-field equivalence, it is natural to define another spherical harmonics based electric dipole as [17, 23, 10, 11, 24]:

(6) |

where is the vacuum permittivity, and the associated scattering power of is [15, 10, 11]:

(7) |

We note here that has been obtained based on the far scattered field expansion and thus has included all the dipolar transverse magnetic components [17, 11]. In contrast, the conventional electric dipole P shown in Eq.(1) has been obtained based on the near-field current integration, and here the current can be expressed as:

(8) |

where is a scalar if the particle is isotropic and a tensor if the particle is radially anisotropic. The conventional electric dipole P is only the lowest order expansion term of shown in Eq.(6). Similarly, the toroidal dipole T shown in Eq.(3) is the next higher order expansion of . So basically can be expressed as [14, 15, 10, 11]:

(9) |

where denote higher order correcting terms of with respect to [14, 15, 25, 26]. As is shown in Eq.(9) the multipole expansion of the far scattered fields based on spherical harmonics is complete [Eqs.(5)-(6)] while the conventional multipole expansion based on the integration of the near-field current is incomplete [Eqs.(1)-(3)]. As the normalized particle radius () is getting larger and larger (that is to say the particle radius R is getting larger and larger or the functioning wavelength is getting smaller or smaller), higher order correcting terms including the toroidal dipole T would arise. That is to say, multipole expansions even including toroidal multipoles are still incomplete. For sufficiently larger , other higher order correcting terms must also be incorporated [14, 15, 10, 11, 13, 25, 26].

## Appendix C Results and Discussions

In Fig. 1(a) the schematic of the scattering configuration is shown: the pane wave is scattered by a spherical particle, which could be homogenous and isotropic (metal or dielectric) with refractive index or radially anisotropic with radial permittivity and transverse (perpendicular to the radial direction) permittivity . The anisotropy parameter is defined as . The plane wave is polarized along direction and propagating along direction. In Fig. 1(b) we show the excited toroidal dipole T within the particle by the magnetic field H and current distribution J.

### c.1 TD excitation within homogenous and isotropic dielectric particles

As a first step we study the plane wave scattering by a homogenous and isotropic dielectric sphere of nm and . The scattered power spectra [in terms of , P and T] are shown in Fig. 2(a). Four points A-D are selected and indicated by black dots in Fig. 2(a) and the corresponding near-field distributions are shown in Fig. 2(A-D): the vector-plots correspond to the displacement field on the plane [only the field within the particle is plotted, as is also the case for Fig. (2)-Fig. (3)] and the color-plots correspond to the normalized magnetic fields along direction ().

At point A of , the quasi-static approximation can be applied, and thus the total dipolar field can be approximated as a conventional electric dipole: . This is due to the fact that in the quasi-static spectral regime, T and other higher order correcting terms in Eq. (9) can be neglected. To further confirm this, we show the near fields at point A in Fig. 2(A), where a typical electric dipole field is exhibited.

At Points B and C of and respectively, which is beyond the quasi-static spectral regime, the quasi-static approximation breaks down and thus the contributions of TD have to be included [see Eq. (9)]. At those points, P been totally suppressed and there are only efficient TD excitation: . The corresponding near-field distributions are shown respectively in Fig. 2(B-C). It is worth mentioning that the field distribution at points B in Fig. 2(B) corresponds exactly to the filed and currents shown in Fig. 1(b), which is a typical and fundamental representation of TD [the magnetic field is confined within the particle and there are two field intensity peaks (thus two oppositely circulating current loops on plane) across]. The field distribution at point C shown in Fig. 2(C), though different from what is shown in Fig. 1(B), also corresponds to a TD, which however is of a higher mode number with more field lobes [27, 19, 28] [the magnetic field is also confined within the particle while there are four field intensity peaks across].

At Points D of which is far beyond the quasi-static spectral regime, even the incorporation of TD will not account for all the dipolar scattering, which can be confirmed in Fig. 2(a): though at point D the overall dipolar scattering is negligible [, dashed black curve], TD response T is dominant over P and this means that they cannot interfere totally in a destructive way, thus catcalling the scattering of each other [11, 13]; As a result, besides P and T, there must be higher order correcting terms for as shown in Eq. (9). This further proves that the multipole expansions based on near-field integration is still incomplete even of the toroidal responses have been considered. Besides conventional electric, magnetic and toroidal multipoles, there are other terms which will become more and more significant as the particle dimension gets larger and larger compared to the functioning wavelength [14, 15, 25, 26]. The corresponding near-field distribution of point D is shown in Fig. 2(D).

### c.2 TD excitation within homogenous and isotropic silver particles

As a next step we study the plane wave scattering by a homogenous and isotropic Ag sphere of nm and the permittivity of silver is adopted from the experimental data in Ref. [29]. The scattered power spectra [in terms of , P and T] are shown in Fig. 3(a). Two points A and B are selected and indicated by black dots in Fig. 3(a) and the corresponding near-field distributions are shown in Fig. 3(A-B). Similar to the dielectric sphere discussed in the subsection above: At point A of in the quasi-static regime, P is the only present term [, see Fig. 3(a)] and the near field exhibits a conventional electric dipole configuration [see Fig. 3(A)]; At point B of beyond the quasi-static spectral regime, the TD term arises [see Fig. 3(a)] and the corresponding near field shown in Fig. 3(B) is similar to that shown in Fig. 2(B). The difference is that in Fig. 2(B) there is only TD excitation while in Fig. 3(B) there are also P and higher order correcting terms, similar to point D in Fig. 2(D). In a word, here we show that for homogenous metallic particles, T can be effectively excited beyond the quasi-static spectral regime, which can be stronger over P [Fig. 3(a)].

### c.3 Tuning the TDs excited within homogenous dielectric particle through radial anisotropy

At the end, we study the scattering of radially anisotropic dielectric spheres of radius nm. The anisotropy parameter is and the transverse permittivity is fixed at . For radially anisotropic particles, both the Mie scattering coefficients and the near-fields can be calculated analytically in a similar manner to that of isotropic particles through modifying the expressions of the orders of related spherical Bessel functions to include the anisotropy parameters [30, 31, 32, 28]. As a result, current J [as is shown in Eq. (8) where now is a tensor], P, T, and their scattered power for radially anisotropic spherical particles can also be calculated directly.

The scattered power spectra are shown in Fig. 4(a) of and in Fig. 4(b) of with contributions from P, T, and . Four points where only TD are excited [] with both suppressed P and higher order correcting terms are indicated by black dots (A: , B: , C: , D: ). The corresponding near fields are shown respectively in Fig. 4(A)-(D). It is known that through employing radial anisotropy, can be tuned to larger/smaller wavelengths (smaller/larger ) through employing smaller/larger [28]. Similarly, TD excited within dielectric spheres can also be tuned by adjusting the radial anisotropy parameters, which is clear according to Fig. 4(a) and (b). Compared to the isotropic case of shown in Fig. 2(a), the TD responses [e.g. characterized by the positions where only TD is excited, as indicated by points B and C in Fig. 2(a) and points A-D in Fig. 4(a-b)] can be tuned to larger wavelengths (smaller ) through making [see Fig. 4(a)] and to smaller wavelengths (larger ) through rendering [Fig. 4(b)]. For the near-field distributions, both TD with low mode number [see Fig. 4(A) and (C), which is similar to Fig. 2(B)] and higher mode number [see Fig. 4(B) and (D), which is similar to Fig. 2(C)] can be effectively tuned within radially anisotropic dielectric particles.

## Appendix D Conclusions and Outlook

To conclude, we revisit the fundamental problem of Mie scattering of spherical particles from the new perspective of TD excitation and tuning. Besides the demonstration of TDs within dielectric particles, we further demonstrate that TDs can be effectively excited even within homogenous metallic particles as long as it is beyond the quasi-static spectral regime. It is further shown that the TDs excited can be spectrally tuned by adjusting the radial anisotropy parameters, which paves the way for further flexible manipulations of toroidal responses.

It is worth mentioning that here we have confined our study to the electric dipolar responses (transverse magnetic dipolar responses [17, 11]) and similar studies can certainly conducted for the magnetic dipolar responses (transverse electric dipolar responses [17, 11]), or both transverse magnetic and transverse electric responses of higher orders. Such multipole expansions including toroidal responses can also be implemented for other non-spherical structures where the fields and currents, and thus related multipoles can be calculated through numerical methods. As for the expansion of the dipolar scattering shown in Eq. (9), we have only discussed in detail the correcting terms up to toroidal dipoles, and similarly higher order correcting terms can also be further investigated [14, 15, 10, 11, 13, 25, 26]. For the tuning of TDs excited, Here in this paper we employ only the radial anisotropy, and it is natural to expect that the other kinds of anisotropy, nonlinear and/or gain effects can also be used for further more flexible manipulations of toroidal multipole responses.

We believe that the incorporation of toroidal multipoles into scattering nanoparticles and the explorations into their tunability have not only expanded our understanding of the fundamental Mie scattering physics, but also would be of great significance for both far-field scattering pattern shaping [22, 33, 34] and near-field manipulations, which could possibly stimulate many scattering related fundamental researches and applications such as nano-lasing [35, 36], sensing [37], nanoantennas [38], and photovoltaic devices [39]. Moreover, investigations into light-matter interactions involving toroidal multipoles can be possibly conducted on the recently emerging topological [40] and/or two-dimensional [41] photonic platforms, and probably can lay a solid foundation for the decisive solution of the dynamical Aharonov-Bohm effect controversy concerning quantum information processing [42, 43].

## 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.