# Enhancement of artificial magnetism via resonant bianisotropy

###### Abstract

All-dielectric ”magnetic light” nanophotonics based on high refractive index nanoparticles allows controlling magnetic component of light at nanoscale without having high dissipative losses. The artificial magnetic optical response of such nanoparticles originates from circular displacement currents excited inside those structures and strongly depends on geometry and dispersion of optical materials. Here a new approach for increasing magnetic response via resonant bianisotropy effect is proposed and analyzed. The key mechanism of enhancement is based on electric-magnetic interaction between two electrically and magnetically resonant nanoparticles of all-dielectric dimer nanoantenna. It was shown that proper geometrical arrangement of the dimer in respect to the incident illumination direction allows flexible control over all vectorial components of magnetic polarizability, tailoring the later in the dynamical range of 100 and enhancement up to 36 relative to performances of standalone spherical particles. The proposed approach provides pathways for designs of all-dielectric metamaterials and metasurfaces with strong magnetic response.

## I Introduction

Intrinsic magnetic polarizabilities of natural materials have strong frequency dependence with the fundamental cut-off in GHz range, originating from relatively low spin and orbital susceptibilities Polydoroff (1960). Recently, effective polarization currents in subwavelength structured loops, organized in ordered arrays, became sources of high-frequency artificial magnetism Smith et al. (2000). Nanostructured noble metals, supporting localized plasmon resonances could serve as building blocks for metamaterials with artificial magnetic polarizability Shalaev (2007). However, inherent material losses set severe limitations on performances of such structures Khurgin and Sun (2011). Another approach for obtaining magnetic optical response is to employ circular displacement currents in high-index dielectric nanoparticles Kuznetsov et al. (2012). This is the essence of so-called all-dielectric nanophotonics which opened the way to control magnetic component of light at nanoscale without high-dissipation, inherent for metallic nanostructures Popa and Cummer (2008); Evlyukhin et al. (2012); Miroshnichenko and Kivshar (2012); Ginn et al. (2012); Fu et al. (2013); Slovick et al. (2013a); Shcherbakov et al. (2014a). The ”magnetic light” concept found use in various applications, such as nanoantennas Krasnok et al. (2014, 2012), quantum interface for NV-centers Krasnok et al. (2015), photonic topological insulators Slobozhanyuk et al. (2015), broadband perfect reflectors Slovick et al. (2013b), waveguides Savelev et al. (2014), cloacking Mirzaei et al. (2015), harmonics generation Shcherbakov et al. (2014b), wave-front engineering, and dispersion control Decker et al. (2015).

Magnetic response of a dielectric particle strongly depends on it’s refractive index, shape, and external environment. The eigen frequencies of electric and magnetic resonances could span the entire visible range. However, the value of those multipole moments is limited by dispersion properties of optical materials.

Here a new approach for tailoring magnetic response of dielectric nanoparticles via resonant effect of bianisotropy is proposed. Microscopically, bianisotropy is the effect of magneto-electric coupling, where electric polarization induces magnetic and vice versa. The signature of this effect appears in the constitutive relations, e.g the dependence of electrical induction also on magnetic field and magnetic induction also on electric field Kriegler et al. (2010); Asadchy et al. (2015). Bianisotropy is used for achieving high values of effective polarizability in metamaterials Marques et al. (2002), unique properties of metasurfaces Asadchy et al. (2015), and designed directivity of nanoantennas Arango and Koenderink (2013). Previously, it was shown that the interaction of dielectric nanoparticles with substrates may increase induced magnetic moment due to the effect of non-resonant bianisotropy Markovich et al. (2014). The approach proposed here is based on electric-magnetic interaction between two resonant nanoparticles of dimer nanoantenna (see Fig. 1). The nanoparticles were designed in the following way: the electric dipolar resonance of the bigger sphere overlaps with magnetic response of the smaller one. In this case, the effect of bi-resonance anisotropy is achieved: the resonant electric moment of the bigger nanoparticle induces the additional magnetic moment in the smaller one, tailoring its overall response.

The manuscript is organized as follows: first, the optical properties of isolated spherical particles are briefly discussed in the context of the resonance tuning. Next, the analytical coupled dipoles formulation of the problem is developed and verified by numerical modeling. The expression describing the magnetic moment of the nanoparticle, considering the bianisotropy, effect is derived. Furthermore, it was shown that proper geometrical arrangement of the dimer in respect to the incident illumination direction allows achieving additional vectorial component of magnetic polarization.

## Ii Coupled dipoles theory and numerical results

In order to obtain the dimer design, properties of isolated components will be briefly discussed. First, isolated silicon sphere of radius nm, having magnetic dipolar resonance at wavelength 477 nm, i.e. in the visible range is considered. The material dispersion of crystalline silicone (c-Si) is taken from Palik (1997). The scattering cross-section of the nanoparticle have been calculated using CST Microwave Studio, and corresponding results as the function of wavelength are presented in Fig. 2(a) (blue curve). Electric dipolar resonance is blue shifted in respect to the magnetic one and appears on the spectrum at 418 nm wavelength. This Mie resonances hierarchy is red shifted when the radius of nanoparticle is increased. It is worth noting, that high-order quadrupole resonances, also contributing to the scattering cross-section, are suppressed in this region due to perceptible losses of silicon. The bigger sphere with radius of nm exhibits its electric resonance at the similar spectral position where 52-nm sphere exhibits resonant magnetic response (Fig. 2a, red curve). The electric and magnetic polarizabilities which are associated with induced electric and magnetic moments have been also calculated (Fig. 2b). Results of numerical simulations (scatter plots) are in good agreement with analytical Mie theory (lines plot) verifying the validity of the numerical tool. Values of those electric and magnetic moments will be subsequently used in the analytical model, based on discrete dipoles approximation.

Next, the scattering of a plane wave on all-dielectric dimer nanoantenna will be analyzed. The electromagnetic scattering problem could be solved by employing Coupled Electric and Magnetic Dipole Approximation Hopkins et al. (2015) (CEMDA). In this method, complex nanostructures are represented by converging series of point electric and magnetic dipoles Novotny and Hecht (2012), while the problem of two spheres could be approximated by only two. This approximation is particularly accurate, if the gap between the spheres is bigger than their radii Evlyukhin et al. (2010). The goal of the subsequent analytical analysis is obtaining a simple formulation, underlining the interference phenomena, affecting the magnetic dipolar polarizability of the smaller particle and revealing the bianisotropic nature of the interaction. The high-order modes and their interplay have been taken into account using the full-wave numerical calculations using the CST Microwave Studio. Following the CEMDA method, total electromagnetic fields [both electric () and magnetic ()] were decomposed into incident and scattered components:

(1) |

where indices denote first (bigger) and second (smaller) nanoparticles, respectively, and are electric and magnetic fields of the incident plane wave, and are the full electric and magnetic fields at the position of the first (electric) nanoparticle (with radius-vector ), and are the electric and magnetic fields scattered by the second (smaller) nanoparticle at the first particle center, and analogously for the second particle. Both particles are polarized by the incident field, as well as by the scattered one:

(2) |

where and are the electric and magnetic polarizabilities of single particle :

(3) |

and is the dielectric constant, is the imaginary unit. The coefficients and are called the first order Mie scattering amplitudes, and they can be written in simplified form which is suitable for analysis: , where the following notation is introduced:

(4) |

where , , - radius of the sphere , - the permittivity of the material (silicon, here), is the wave number in the surrounding space, and is a wavelength.

The values of single particle polarizabilities, calculated using Eqs. (3) and (II), are in a perfect agreement with our numerical calculations [see Fig.2(b)]. Thus, the scattered fields of the dipoles can be obtained through the Green function of a point dipole in free-space Novotny and Hecht (2012), following the method CEMDA:

(5) |

where and are the wavenumber and speed of light in vacuum. The Green’s function of a point dipole in free-space is well known:

(6) |

where is the radius vector, connecting the center of the first dipole (coordinate origin, or center of the bigger particle) with the second one (center of the smaller particle). The tensor and vector have the following symmetry of indices permutation: and .

Both nanoparticles, being isolated, have 3-fold degenerated (magnetic and electric) dipolar resonances, oriented along the unit vectors of a Cartesian coordinate system. The excitation of an isolated sphere is solely defined by the polarization of the incident wave – for example, linearly polarized beam will excite only one of 3 components of the dipolar mode. However, the geometry of the coupled dipoles together with the excitation (not necessarily coinciding with one of the symmetry axis of the system) will break the degeneracy and, as the result, all three vectorial components must be taken into account. The obtained set of equations can be solved analytically by means of the matrix inversion or numerically in the same fashion. In order to verify the validity of the proposed theoretical model, we consider a particular case, where the system of Eqs. (II) has a simple and intuitive solution. Arranging the nanoparticles along the x-axis and exciting the system with linearly polarized plane wave along y-axis with angles being and , the set of 12 coupled equations was reduced just for 4, since the symmetry considerations allow the moment components to be induced only along z-axis (magnetic dipole) and y-axis (electric dipole). Similar configuration was studied in Noskov et al. (2012) for hybrid metal-dielectric nanoantennas with non-resonant dielectric nanoparticles. This set of coupled equations has particularly simple solution, which agrees well with the numerical calculations. The values of the magnetic polarizability enhancement, calculated both analytically and numerically as the function of the distance between nanoparticles are shown in Fig. 3. It is clearly seen, that the analytical and numerical models agree with each other and showing 15 enhancement of the magnetic moment of the single nanoparticle, as the result of the bianisotropic coupling. For clear understanding of oscillatory behaviour of the magnetic moments, the simplest model, when the polarizabilities and are supposed to be zero, was considered. In this case the simple formula can be obtained in form showing the bianisotropic nature of the effect:

(7) |

The qualitative analysis of magnetic moment enables to observe the clear interference phenomena – the direct excitation of magnetic dipole by the plane wave [ term in Eq. (7)] and the contribution of the scattered field through the term . Thus, oscillation behavior of magnetic moment enhancement is the result of the interference phenomenon (see Fig. 3). For longer distances the coupling between the particles becomes weaker and converges to the value of the isolated particle. An almost perfect fit of full numerical modeling with the CEMDA method enables to use the later for analysis of more complex structures without involving full wave simulations. This approach will be subsequently employed.

## Iii Vectorial structure of magnetic moments

In the subsequent studies the illumination was chosen to propagate along -axis and being polarized along . There are three geometrical parameters, affecting the magnetic polarizability of the smaller particle: the distance between the spheres’ centres , and the angles and . The angular dependence of the induced magnetic moment components of the smaller particle for 3 characteristic separation distances , , and nm (gap between the nanoparticles’ surfaces) will be studied next. In the case of noninteracting nanoparticles there is the only one non-zero component of nanoparticle’s magnetic moment – codirected with the magnetic field. In the case of bianisotropic coupling there are also - and - non-zero components of the magnetic moment (, ). Fig. 4 shows the three-dimensional angular dependencies of normalized magnetic moment components (, ) of the smaller nanoparticle, where is the polarizability of the single smaller nanoparticle. These results have been obtained by exact solving Eqs. (II) taking into account the electric and magnetic responses of both nanoparticles. As it could be seen, the induced components are smaller than the main moment. Those additional vectorial components of magnetic moments increase when the distance between the particles goes down and may being up to 50% for nm. But it should be noted that for such small distances between the particles accuracy of CEMDA method is low, sow this value is rough Evlyukhin et al. (2010). When the distance between the particles goes up, secondary vectorial components of magnetic moment decrease. Therefore, the overall variation of magnetic moments is mainly determined by its -component. Diagrams of -component of the magnetic moment are asymmetrical and their forms are dissimilar for different separation distances between the particles (). For nm enhancement of magnetic moment is weak, because of the vanishing coupling between the particles. For 2 m this diagram is symmetrical with good accuracy, as it replicates the performance of the isolated particle. For smaller distances the maximum of magnetic moment corresponds to the cases (, and ), and for bigger distances maximum of magnetic moment corresponds to the case (, ). It should be noted that case (, and ) corresponds to magnetic dipoles coupling only, while the case (, ) corresponds to electric-magnetic dipoles coupling only. Additional analysis shows that the maximum of relative magnetic moment enhancement is nearly and it is achieved for nm for the case of (, ). The performance of the dimer with this separation distance will be investigated in details hereafter. Both induced magnetic and electric moments of the smaller nanoparticle were calculated using the theoretical model CEMDA. Enhancement of the magnetic polarizability of smaller nanoparticle for different angular arrangements (the distance =168 nm is kept constant) is shown in the Fig. 5. These dependences show that the effect has a resonance character. Moreover, the maximum value is achieved at the wavelength of 480 nm i.e. at the electric dipole resonance of bigger nanoparticle and magnetic dipole resonance of smaller one underlining the impact of the resonant nature.

## Iv Outlook and conclusions

Coupled particles approach for controlling magnetic polarizabilities of nanoscale spheres was proposed. While standalone nanoparticles allow obtaining dipolar and high-multipolar resonant responses at the desired wavelength in the visible range specifying their radii and materials, a system of two coupled nanoparticles possesses more degrees of freedom. Altering the radii of both nanoparticles allows to investigate the impact of all the combinations of multipolar coupling effects on the properties of the system. Naturally, the amplitude of the effects is dependent on their mutual displacement, and it decreases at longer distances between the particles, as it was proved both analytically and numerically in the case of dipolar magnetic-electric coupling. Symmetry considerations allow exciting only one dominant induced electric and magnetic moment component when a plane wave is incident upon a single nanoparticle. It was shown, that for a nanoparticle system, altering the spherical angles and distance between nanoparticles allows to excite all vectorial components of moments simultaneously. The secondary components being up to 50 of the amplitude value of dominant ones. Furthermore, these parameters also define the spectral position of deeps and peaks of the dominant electric moment components and amplitude values, allowing full on-demand control of the electromagnetic properties of the system of coupled nanoparticles. The proposed approach can find use in designs of more complex structures such as all-dielectric metamaterials and metasurfaces with strong magnetic responses.

## V Acknowledgments

This work was financially supported by the Ministry of Education and Science of the Russian Federation (14.584.21.0009 10).

## References

- Polydoroff (1960) W. J. Polydoroff, High-Frequency Magnetic Materials: Their Characteristics and Principal Applications (Wiley, 1960).
- Smith et al. (2000) D. R. Smith, W. J. Padilla, D. C. Vier, S. C. Nemat-Nasser, and S. Schultz, Phys. Rev. Lett. 84, 4184 (2000).
- Shalaev (2007) V. M. Shalaev, Nature Photon. 1, 41 (2007).
- Khurgin and Sun (2011) J. B. Khurgin and G. Sun, Appl. Phys. Lett. 99, 211106 (2011).
- Kuznetsov et al. (2012) A. I. Kuznetsov, A. E. Miroshnichenko, Y. H. Fu, J. Zhang, and B. Luk’yanchuk, Sci. Rep. 2, 492 (2012).
- Popa and Cummer (2008) B. I. Popa and S. A. Cummer, Phys. Rev. Letters 100, 207401 (2008).
- Evlyukhin et al. (2012) A. B. Evlyukhin, S. M. Novikov, U. Zywietz, R. L. Eriksen, C. Reinhardt, S. I. Bozhevolnyi, and B. N. Chichkov, Nano Lett. 12, 3749 (2012).
- Miroshnichenko and Kivshar (2012) A. E. Miroshnichenko and S. Y. Kivshar, Nano Lett. 12, 6459 (2012).
- Ginn et al. (2012) J. C. Ginn, I. Brenner, D. W. Peters, J. R. Wendt, J. O. Stevens, P. F. Hines, L. I. Basilio, L. K. Warne, J. F. Ihlefeld, P. G. Clem, and M. B. Sinclair, Phy. Rev. Lett. 108, 097402 (2012).
- Fu et al. (2013) Y. H. Fu, A. I. Kuznetsov, A. E. Miroshnichenko, Y. F. Yu, and B. Luk’yanchuk, Nat. Comms. 4, 1527 (2013).
- Slovick et al. (2013a) B. Slovick, Z. G. Yu, M. Berding, and S. Krishnamurthy, Phys. Rev. B 88, 165116 (2013a).
- Shcherbakov et al. (2014a) M. R. Shcherbakov, D. N. Neshev, B. Hopkins, A. S. Shorokhov, I. Staude, E. V. Melik-Gaykazyan, M. Decker, A. A. Ezhov, A. E. Miroshnichenko, I. Brener, A. A. Fedyanin, and Y. S. Kivshar, Nano Lett. 14, 6488 (2014a).
- Krasnok et al. (2014) A. E. Krasnok, C. R. Simovski, P. A. Belov, and Y. S. Kivshar, Nanoscale 6, 7354 (2014).
- Krasnok et al. (2012) A. E. Krasnok, A. E. Miroshnichenko, P. A. Belov, and Y. S. Kivshar, Opt. Express 20, 20599 (2012).
- Krasnok et al. (2015) A. E. Krasnok, A. Maloshtan, D. N. Chigrin, Y. S. Kivshar, and P. A. Belov, Laser Photon. Rev. DOI: 10.1002/lpor.201400453 (2015).
- Slobozhanyuk et al. (2015) A. P. Slobozhanyuk, A. N. Poddubny, A. E. Miroshnichenko, P. A. Belov, and Y. S. Kivshar, Phys. Rev. Lett. 114, 123901 (2015).
- Slovick et al. (2013b) B. Slovick, Z. G. Yu, M. Berding, and S. Krishnamurthy, Phys. Rev. B 88, 165116 (2013b).
- Savelev et al. (2014) R. S. Savelev, D. S. Filonov, P. V. Kapitanova, A. E. Krasnok, A. E. Miroshnichenko, P. A. Belov, and Y. S. Kivshar, Applied Physics Letters 105, 181116 (2014).
- Mirzaei et al. (2015) A. Mirzaei, A. E. Miroshnichenko, I. V. Shadrivov, and Y. S. Kivshar, Sci. Reports 5, 9574 (2015).
- Shcherbakov et al. (2014b) M. R. Shcherbakov, D. N. Neshev, B. Hopkins, A. S. Shorokhov, I. Staude, E. V. Melik-Gaykazyan, M. Decker, A. A. Ezhov, A. E. Miroshnichenko, I. Brener, A. A. Fedyanin, and Y. S. Kivshar, Nano Lett. 14, 6488 (2014b).
- Decker et al. (2015) M. Decker, I. Staude, M. Falkner, J. Dominguez, D. N. Neshev, I. Brener, T. Pertsch, and Y. S. Kivshar, Adv. Opt. Mater. DOI: 10.1002/adom.201400584 (2015).
- Kriegler et al. (2010) C. E. Kriegler, M. S. Rill, S. Linden, and M. Wegener, IEEE journal of selected topics in quantum electronics 16, 367 (2010).
- Asadchy et al. (2015) V. S. Asadchy, Y. Radi, J. Vehmas, and S. A. Tretyakov, Phy. Rev. Lett. 114, 095503 (2015).
- Marques et al. (2002) R. Marques, F. Medina, and R. Rafii-El-Idrissi, Phy. Rev. B 65, 144440 (2002).
- Arango and Koenderink (2013) F. B. Arango and A. F. Koenderink, New J. Phys. 15, 073023 (2013).
- Markovich et al. (2014) D. L. Markovich, P. Ginzburg, A. K. Samusev, P. A. Belov, and A. V. Zayats, Opt. Express 22, 10693 (2014).
- Palik (1997) E. D. Palik, Handbook of Optical Constants of Solids (Elseveir Inc., 1997).
- Hopkins et al. (2015) B. Hopkins, D. S. Filonov, A. E. Miroshnichenko, F. Monticone, A. Alu, and Y. S. Kivshar, ACS Photonics 2, 724 (2015).
- Novotny and Hecht (2012) L. Novotny and B. Hecht, Principles of Nano-Optics (Cambridge University Press, 2012).
- Evlyukhin et al. (2010) A. B. Evlyukhin, C. Reinhardt, A. Seidel, B. S. Luk’yanchuk, and B. N. Chichkov, Phys. Rev. B 82, 045404 (2010).
- Noskov et al. (2012) R. E. Noskov, A. E. Krasnok, and Y. S. Kivshar, New J. Phys. 14, 093005 (2012).