Mimicking electromagnetically induced transparency in the magneto-optical activity of magnetoplasmonic nanoresonators

We show that the interaction between a plasmonic and a magnetoplasmonic metallic nanodisk leads to the appearance of magnetooptical activity in the purely plasmonic disk induced by the magnetoplasmonic one. Moreover, at specific wavelengths the interaction cancels the net electromagnetic field at the magnetoplasmonic component, strongly reducing the magneto-optical activity of the whole system. The MO activity has a characteristic Fano spectral shape, and the resulting MO inhibition constitutes the magneto-optical counterpart of the electromagnetic induced transparency. ©2013 Optical Society of America OCIS codes: (250.5403) Plasmonics; (240.6680) Surface plasmons; (230.3810) Magneto-optic systems; (160.3820) Magneto-optical materials; (260.3910) Metal optics; (350.4238) Nanophotonics and photonic crystals. References and Links 1. E. Prodan, C. Radloff, N. J. Halas, and P. Nordlander, “A hybridization model for the plasmon response of complex nanostructures,” Science 302(5644), 419–422 (2003). 2. P. Mühlschlegel, H.-J. Eisler, O. J. F. Martin, B. Hecht, and D. W. Pohl, “Resonant optical antennas,” Science 308(5728), 1607–1609 (2005). 3. D. P. Fromm, A. Sundaramurthy, P. J. Schuck, G. Kino, and W. E. Moerner, “Gap-dependent optical coupling of single «bowtie» nanoantennas resonant in the visible,” Nano Lett. 4(5), 957–961 (2004). 4. P. Ghenuche, S. Cherukulappurath, T. H. Taminiau, N. F. van Hulst, and R. Quidant, “Spectroscopic mode mapping of resonant plasmon nanoantennas,” Phys. Rev. Lett. 101(11), 116805 (2008). 5. K.-H. Su, Q.-H. Wei, X. Zhang, J. J. Mock, D. R. Smith, and S. Schultz, “Interparticle coupling effects on plasmon resonances of nanogold particles,” Nano Lett. 3(8), 1087–1090 (2003). 6. W. Rechberger, A. Hohenau, A. Leitner, J. R. Krenn, B. Lamprecht, and F. R. Aussenegg, “Optical properties of two interacting gold nanoparticles,” Opt. Commun. 220(1-3), 137–141 (2003). 7. A. Dmitriev, T. Pakizeh, M. Käll, and D. S. Sutherland, “Gold-silica-gold nanosandwiches: tunable bimodal plasmonic resonators,” Small 3(2), 294–299 (2007). 8. L. V. Brown, H. Sobhani, J. B. Lassiter, P. Nordlander, and N. J. Halas, “Heterodimers: plasmonic properties of mismatched nanoparticle pairs,” ACS Nano 4(2), 819–832 (2010). 9. C. Wadell, T. J. Antosiewicz, and C. Langhammer, “Optical absorption engineering in stacked plasmonic AuSiO2-Pd nanoantennas,” Nano Lett. 12(9), 4784–4790 (2012). 10. M. Hentschel, M. Saliba, R. Vogelgesang, H. Giessen, A. P. Alivisatos, and N. Liu, “Transition from isolated to collective modes in plasmonic oligomers,” Nano Lett. 10(7), 2721–2726 (2010). 11. J. B. Lassiter, H. Sobhani, J. A. Fan, J. Kundu, F. Capasso, P. Nordlander, and N. J. Halas, “Fano resonances in plasmonic nanoclusters: geometrical and chemical tunability,” Nano Lett. 10(8), 3184–3189 (2010). 12. S. Zou, N. Janel, and G. C. Schatz, “Silver nanoparticle array structures that produce remarkably narrow plasmon lineshapes,” J. Chem. Phys. 120(23), 10871–10875 (2004). 13. B. Auguié and W. L. Barnes, “Collective resonances in gold nanoparticle arrays,” Phys. Rev. Lett. 101(14), 143902 (2008). 14. J. Aizpurua, G. W. Bryant, L. J. Richter, F. J. G. de Abajo, B. K. Kelley, and T. Mallouk, “Optical properties of coupled metallic nanorods for field-enhanced spectroscopy,” Phys. Rev. B 71(23), 235420 (2005). #194818 $15.00 USD Received 30 Jul 2013; revised 23 Sep 2013; accepted 24 Sep 2013; published 4 Nov 2013 (C) 2013 OSA 4 November 2013 | Vol. 21, No. 22 | DOI:10.1364/OE.21.027356 | OPTICS EXPRESS 27356 15. A. Sundaramurthy, K. Crozier, G. Kino, D. Fromm, P. Schuck, and W. Moerner, “Field enhancement and gapdependent resonance in a system of two opposing tip-to-tip Au nanotriangles,” Phys. Rev. B 72(16), 165409 (2005). 16. G. Dolling, C. Enkrich, M. Wegener, J. F. Zhou, C. M. Soukoulis, and S. Linden, “Cut-wire pairs and plate pairs as magnetic atoms for optical metamaterials,” Opt. Lett. 30(23), 3198–3200 (2005). 17. N. Verellen, Y. Sonnefraud, H. Sobhani, F. Hao, V. V. Moshchalkov, P. Van Dorpe, P. Nordlander, and S. A. Maier, “Fano resonances in individual coherent plasmonic nanocavities,” Nano Lett. 9(4), 1663–1667 (2009). 18. Z.-J. Yang, Z.-S. Zhang, L.-H. Zhang, Q.-Q. Li, Z.-H. Hao, and Q.-Q. Wang, “Fano resonances in dipolequadrupole plasmon coupling nanorod dimers,” Opt. Lett. 36(9), 1542–1544 (2011). 19. B. Luk’yanchuk, N. I. Zheludev, S. A. Maier, N. J. Halas, P. Nordlander, H. Giessen, and C. T. Chong, “The Fano resonance in plasmonic nanostructures and metamaterials,” Nat. Mater. 9(9), 707–715 (2010). 20. Y. Francescato, V. Giannini, and S. A. Maier, “Plasmonic systems unveiled by Fano resonances,” ACS Nano 6(2), 1830–1838 (2012). 21. B. Gallinet and O. J. F. Martin, “Influence of electromagnetic interactions on the line shape of plasmonic Fano resonances,” ACS Nano 5(11), 8999–9008 (2011). 22. A. E. Cetin and H. Altug, “Fano resonant ring/disk plasmonic nanocavities on conducting substrates for advanced biosensing,” ACS Nano 6(11), 9989–9995 (2012). 23. E. Cubukcu, Yu. Nanfang, E. J. Smythe, L. Diehl, K. B. Crozier, and F. Capasso, “Plasmonic laser antennas and related devices,” IEEE J. Sel. Top. Quantum Electron. 14(6), 1448–1461 (2008). 24. J. N. Anker, W. P. Hall, O. Lyandres, N. C. Shah, J. Zhao, and R. P. Van Duyne, “Biosensing with plasmonic nanosensors,” Nat. Mater. 7(6), 442–453 (2008). 25. J. A. Schuller, E. S. Barnard, W. Cai, Y. C. Jun, J. S. White, and M. L. Brongersma, “Plasmonics for extreme light concentration and manipulation,” Nat. Mater. 9(3), 193–204 (2010). 26. V. V. Temnov, G. Armelles, U. Woggon, D. Guzatov, A. Cebollada, A. Garcia-Martin, J. M. Garcia-Martin, T. Thomay, A. Leitenstorfer, and R. Bratschitsch, “Active magneto-plasmonics in hybrid metal-ferromagnet structures,” Nat. Photonics 4(2), 107–111 (2010). 27. V. I. Belotelov, I. A. Akimov, M. Pohl, V. A. Kotov, S. Kasture, A. S. Vengurlekar, A. V. Gopal, D. R. Yakovlev, A. K. Zvezdin, and M. Bayer, “Enhanced magneto-optical effects in magnetoplasmonic crystals,” Nat. Nanotechnol. 6(6), 370–376 (2011). 28. J. C. Banthí, D. Meneses-Rodríguez, F. García, M. U. González, A. García-Martín, A. Cebollada, and G. Armelles, “High magneto-optical activity and low optical losses in metal-dielectric Au/Co/Au-SiO2 magnetoplasmonic nanodisks,” Adv. Mater. 24(10), OP36–OP41 (2012). 29. G. Armelles, A. Cebollada, A. García-Martín, and M. U. González, “Magnetoplasmonics: combining magnetic and plasmonic functionalities,” Adv. Opt. Mater. 1(1), 10–35 (2013). 30. J. Y. Chin, T. Steinle, T. Wehlus, D. Dregely, T. Weiss, V. I. Belotelov, B. Stritzker, and H. Giessen, “Nonreciprocal plasmonics enables giant enhancement of thin-film Faraday rotation,” Nat Commun 4, 1599 (2013). 31. B. Sepúlveda, J. B. González-Díaz, A. García-Martín, L. M. Lechuga, and G. Armelles, “Plasmon-induced magneto-optical activity in nanosized gold disks,” Phys. Rev. Lett. 104(14), 147401 (2010). 32. S. E. Harris, J. E. Field, and A. Imamoğlu, “Nonlinear optical processes using electromagnetically induced transparency,” Phys. Rev. Lett. 64(10), 1107–1110 (1990). 33. Q. Xu, S. Sandhu, M. L. Povinelli, J. Shakya, S. Fan, and M. Lipson, “Experimental realization of an on-chip alloptical analogue to electromagnetically induced transparency,” Phys. Rev. Lett. 96(12), 123901 (2006). 34. S. Zhang, D. A. Genov, Y. Wang, M. Liu, and X. Zhang, “Plasmon-induced transparency in metamaterials,” Phys. Rev. Lett. 101(4), 047401 (2008). 35. N. Papasimakis, V. A. Fedotov, N. I. Zheludev, and S. L. Prosvirnin, “Metamaterial analog of electromagnetically induced transparency,” Phys. Rev. Lett. 101(25), 253903 (2008). 36. N. Liu, L. Langguth, T. Weiss, J. Kästel, M. Fleischhauer, T. Pfau, and H. Giessen, “Plasmonic analogue of electromagnetically induced transparency at the Drude damping limit,” Nat. Mater. 8(9), 758–762 (2009). 37. T. Pakizeh, A. Dmitriev, M. S. Abrishamian, N. Granpayeh, and M. Käll, “Structural asymmetry and induced optical magnetism in plasmonic nanosandwiches,” J. Opt. Soc. Am. B 25(4), 659–667 (2008). 38. Y. Ekinci, A. Christ, M. Agio, O. J. F. Martin, H. H. Solak, and J. F. Löffler, “Electric and magnetic resonances in arrays of coupled gold nanoparticle in-tandem pairs,” Opt. Express 16(17), 13287–13295 (2008), http://www.opticsexpress.org/abstract.cfm?URI=oe-16-17-13287. 39. H. Fredriksson, Y. Alaverdyan, A. Dmitriev, C. Langhammer, D. S. Sutherland, M. Zäch, and B. Kasemo, “Hole–mask colloidal lithography,” Adv. Mater. 19(23), 4297–4302 (2007). 40. J. B. González-Díaz, A. García-Martín, J. M. García-Martín, A. Cebollada, G. Armelles, B. Sepúlveda, Y. Alaverdyan, and M. Käll, “Plasmonic Au/Co/Au nanosandwiches with enhanced magneto-optical activity,” Small 4(2), 202–205 (2008). 41. A. García-Martín, G. Armelles, and S. Pereira, “Light transport in photonic crystals composed of magnetooptically active materials,” Phys. Rev. B 71(20), 205116 (2005). 42. B. Caballero, A. García-Martín, and J. C. Cuevas, “Generalized scattering-matrix approach for magneto-optics in periodically patterned multilayer systems,” Phys. Rev. B 85(24), 245103 (2012). 43. FDTD Solutions” software from Lumerical Inc. (www.lumerical.com). 44. COMSOL Multiphysics®. #194818 $15.00 USD Received 30 Jul 2013; revised 23 Sep 2013; accepted 24 Sep 2013; published 4 Nov 2013 (C) 2013 OSA 4 November 2013 | Vol. 21, No. 22 | DOI:10.1364/OE.21.027356 | OPTICS EXPRESS 27357 45. A. A. Yanik, A. E. Cetin, M. Huang, A. Artar, S. H. Mousavi, A. Khanikaev, J. H. Connor, G. Shvets, and H. Altug, “Seeing protein monolayers with naked eye through plasmonic Fano resonances,” Proc. Natl. Acad. Sci. U.S.A. 108(29), 11784–11789 (2011). 46. C. Wu, A. B. Khanikaev, R. Adato, N. Arju, A. A. Yanik, H. Altug, and G. Shvets, “Fano-resonant asymmetric metamaterials for ultrasensitive spectroscopy and identification of molecular monolayers,” Nat. Mater. 11(1), 69–75 (2011). 47. P. Offermans, M. C. Schaafsma, S. R. K. Rodriguez, Y. Zhang, M. C


Introduction
In complex optical structures, the electromagnetic interaction between the constituent elements strongly determines the response of the whole system.Plasmonics offer numerous examples of this, with a variety of complex architectures including core-shell particles [1], nanoantennas [2][3][4], dimers [5][6][7][8][9], oligomers [10,11], periodic arrays [12,13], etc.The diverse direct consequences of these interactions carry both fundamental and technological implications.The strong enhancement of the electromagnetic field in subwavelength spatial regions [14,15], the building up of bright and dark modes [16][17][18] or the Fano resonances exhibited in the optical signal [19][20][21], to name a few, open research and application areas in the development of nanoantennas, higher sensitivity sensors, novel telecommunication architectures/devices, etc. [22][23][24][25].In order to control and even boost at will such optical responses, one can make use of the geometry and the spatial arrangement of the constituting elements.These morphological and configurational factors determine how the whole entity interacts with an electromagnetic field.In this scenario, the incorporation of magnetic active character into the plasmonic system by adding a ferromagnetic component, to form the socalled magnetoplasmonic structures, introduces an additional degree of freedom.Examples of the potentiality of magnetoplasmonics are the development of active plasmonic configurations whose properties can be externally tuned by the application of a magnetic field or structures with enhanced magneto-optical activity [26][27][28][29][30].
Focusing our attention on complex interacting plasmonic systems with magneto-optical activity, it is convenient to consider a simple picture by factorizing the component elements in simple units, those excited by the electromagnetic field, i.e. point-like dipoles.A point dipole in the presence of an external, steady, magnetic field, experiences a modification of its dipole moment that depends on the relative dipole-magnetic field orientation.In particular, if they are perpendicularly oriented, the dipole rotates due to the Lorentz force [31], and this rotation can be seen as a new degree of freedom for the interactions in plasmonic structures.Moreover, the Lorentz force depends on the material, being stronger for ferromagnetic components.This makes it possible to design structures where this force is spatially different by selecting different materials in its interior, therefore enriching the interaction pattern.In this work we will explore this phenomenology using nanodisks dimers consisting of a magnetoplasmonic and a plasmonic unit.We will show that, due to the electromagnetic interaction, a pronounced dip in the spectral dependence of the magneto-optical activity is observed, exhibiting a characteristic Fano shape.This reduction of the magneto-optical activity is similar to the electromagnetically induced transparency effects observed in numerous physical systems [32][33][34][35][36].

Coupled oscillators: mechanical and dipolar models
Let us start by using a simplified model system that considers the plasmonic components, from a classical point of view, as masses coupled by springs (coupled spring resonators) [9].In this case, the simple plasmonic structure without external magnetic field would be two point masses coupled in one dimension [Fig.1(a), left side].This system is characterized by two resonant eigenmodes, symmetric and anti-symmetric in character.Now, still classically, if one of the masses is charged and a magnetic field is present (applied along the z-direction), that mass would experience a Lorentz force along the y-direction [Fig.1(a), right side] (Media 1) and therefore it will also move along the y-axis.Due to the interaction, this movement is transferred to the uncharged particle (see Appendix A for the full development of equations governing this system).Going now to the real plasmonic case, simple structures that resemble this two coupled spring resonators system are metallic disks dimers spatially separated by a dielectric.The interaction between the electric dipoles gives rise to bright and dark modes (or more precisely, superradiant and subradiant modes), whose spectral position and character can be tuned by the thickness of the separating dielectric [7,37,38].In the simple model presented in Fig. 1(b), left side, the bright mode corresponds to the configuration in which the electric dipoles are oriented parallel to each other (symmetric mode), while in the dark case the dipoles are

(a) (b)
Aceptado para su publicación en Optics Express DOI: 10.1364/OE.21.027356 oriented antiparallel (anti-symmetric mode).In both cases the electric dipoles are along x-axis (parallel to the electric field of the incoming light).If one of the disks has ferromagnetic character, the application of a magnetic field along z direction will induce a rotation of its associated dipole about this axis [Fig. 1 Extrapolating what it is obtained in the spring-mass case, the transverse oscillation induced in the uncharged mass will manifest here as an induced rotation of the dipole of the non-ferromagnetic metallic disk.The combined rotation of both dipoles will cause polarization conversion effects.In this sense, magnetooptics (MO) is the appropriate tool to explore this new degree of freedom, since it provides a direct measurement of the magnetic field induced polarization conversion.

Optical and magneto-optical response of nanoresonators
The analyzed structures consist of a pure Au nanodisk separated by a SiO 2 spacer from a 2nmCo/4nmAu multilayer nanodisk.The presence of Au in the bottom magnetoplasmonic disk reduces its optical losses as compared to a pure ferromagnetic metal disk [29].The disks were fabricated using hole-mask colloidal lithography and evaporation (electron beam evaporation for Co and SiO 2 layers and thermal evaporation for Au layers) [28,39].In Fig. 2(a) we show a sketch of the internal structure of each disk with the characteristic individual layer thickness.The Au/Co multilayer disk exhibits perpendicular magnetic anisotropy, reducing the magnetic field required to achieve saturation in polar configuration [Fig.2(b)].Fig. 2(c) presents a cross section Scanning Electron Microscopy image of the same structure showing the truncated nanocone shape of the obtained nanoresonators, whose typical lower and upper disk diameters are 110-120nm and 70-90nm respectively.In this cross section image, both the metallic and dielectric components of the nanodisks are clearly distinguishable.Additionally, in Fig. 2   In the left column of Fig. 3 we show the measured extinction spectra as a function of the SiO 2 spacer thickness.As it can be seen, the spectral dependence of the extinction strongly varies as a function of the dielectric thickness, i.e. with the interdisk interaction.For the extreme case of thick SiO 2 (50 nm), both metallic disks interact weakly, and as a consequence two peaks corresponding to the individual disk resonances are observed.Thus, the highenergy peak can be identified with the resonance of the top metallic disk (smaller diameter) and the low energy one with that of the bottom one (larger diameter).As the SiO 2 spacer gets thinner, the position and relative intensity of both resonances changes.In this situation, both modes become hybridized [1,7] and we can no longer talk of individual modes but of complex modes of symmetric and anti-symmetric character.The symmetric mode, occurring at higher energies, strongly couples to the incident light and as a consequence has a larger extinction (bright or superradiant mode).On the other hand, the anti-symmetric mode, occurring at lower energies, couples weakly to the light, resulting in a lower extinction peak (dark or subradiant mode).Qualitatively, the symmetric mode is clearly observed for all the structures, while the anti-symmetric one gradually decreases in intensity, shifting to lower energies as the spacer becomes thinner, becoming practically unobservable for SiO 2 thickness below 15 nm.
The right hand column of Fig. 3 shows the spectral dependence of the MO activity of the same nanostructures, corresponding to the modulus of the complex Polar Kerr rotation, Φ (where i θ ϕ Φ = + , being θ and ϕ the Polar Kerr rotation and ellipticity respectively [29]), measured at normal incidence and magnetic saturation.Unlike simple fully metallic Au/Co/Au nanodisks where MO activity and extinction spectra are directly related [40], the spectra of our nanodisk dimers exhibit a strong dependence on the spacer thickness but with distinctive differences with respect to the extinction spectra.For example, for thick spacer layers only one peak is observed, contrary to the double peaked spectral dependence of the extinction.This peak corresponds to the MO activity of the bottom disk, which actually is the only one with a ferromagnetic constituent.The upper disk, made up of pure Au, does not contribute to the MO activity of the structure.However, as the SiO 2 becomes thinner and the electromagnetic (EM) interaction between both metallic disks increases, a new peak appears at higher energy, whose position hardly depends on the spacer thickness.On the other hand, the low energy MO peak gradually red shifts and loses intensity as the spacer thickness is reduced.Additionally, and contrary to what was observed in the extinction measurements, this low-energy peak is observed all the way down to the thinnest dielectric spacer in the MO case.Special mention deserves the specific MO peak shapes for intermediate spacer thickness (between 15 and 30 nm), where the two MO peaks are energetically close.In this regime the MO spectra are clearly not the result of the convolution of two peaks that would originate from independent modes, but rather exhibit a typical Fano resonance shape, which on the other hand is absent in the extinction spectra.
To understand this phenomenology we have performed rigorous simulations for both the extinction and MO response, based on scattering matrix techniques [41,42] and FDTD [43] as well as FEM [44] codes.The three methods provide equivalent results.The evolution of the spectra as a function of the SiO 2 thickness is depicted as 3D graphs on top of the corresponding columns in Fig. 3, exhibiting a very good agreement with the experimental results, both regarding the intensity, spectral shape, and peak position evolution.In addition to the mentioned extinction and MO activity, the performed simulations allow also obtaining the EM field distribution in regions around the nanodisks.As an example we present in Fig. 4 the near field distribution of the z-component of the electric field, for the sample with 20 nm SiO 2 spacer, in two planes at ten nanometres above the top and below the bottom disks respectively, both in the absence of an external magnetic field and as the difference for magnetic saturations along opposite directions.These distributions have been calculated for the three wavelengths corresponding to the two maxima and the minimum of the MO activity (labelled with red, green and blue diamonds in Fig. 3).The E z component without external magnetic field (extinction measurements situation) is presented in the left panel of Fig. 4. There, for the low energy position (red framed) the EM distribution for each metallic disk resembles that of a point dipole oscillating along the xdirection, with similar E z intensity for both upper and lower disk, and with the individual dipoles (represented by the black arrows in Fig. 4) oriented antiparallel.(Note that the E z field of a dipole changes sign when it is seen from above or from below).This situation therefore corresponds to the aforementioned anti-symmetric configuration.On the other hand, for the high energy position (blue framed) the EM distribution for each metallic disk corresponds also to a dipole-like distribution along the x-direction but now the individual dipoles corresponding to the disks are oriented parallel to each other (symmetric mode).As it can be seen, the intensity of the E z component for both disks is, again, similar, but larger than that for the low energy case.This is consistent with the higher extinction value obtained for the

(a) (b)
symmetric mode respect to the anti-symmetric one.Finally, at the intermediate energy (green framed), the distribution is still dipolar-like and would correspond to a symmetric-like mode, but the intensity distribution is not equally distributed, with a larger value of the E z component in the upper disk than for the lower one.Next, in the right panel of Fig. 4 we show the difference between E z components at magnetic saturation along opposite directions, E z (H)-E z (-H) (MO measurements situation) which reflects the magnetic field effect on the field distributions of the system eliminating the purely optical contribution.For the three energies, a dipolar-like distribution is still observed, but the resulting balance between magnetic saturation along opposite directions is an E z component corresponding to a point dipole oscillating along the y-direction.In other words, the effect of the magnetic field is to induce a "dipole" along the y-direction that is two orders of magnitude smaller in intensity than the dipole generated along the x-direction.
On the other hand, regarding the symmetry character of the resulting modes, it is maintained for high and low energies with respect to the case of no magnetic field applied: anti-symmetric and symmetric for low and high energy respectively.Regarding the intensities of the E z fields, it is similar for top and bottom disks and for both energies.The induction by the magnetic field of a dipole along y-direction is pretty intuitive for the bottom disk, since it is simply due to the presence of ferromagnetic material in it.What is in principle not as intuitive is the presence of the same component in the upper disk, since it has no ferromagnetic nature.This y-component in the upper disk is induced by the "dipole" of the bottom disk, reaching a sizeable magnitude.More interestingly, the situation for the intermediate energy, corresponding to the minimum of the MO activity, is quite different.In this case both lower and upper regardless the presence of a ferromagnetic component in its interior, exhibit almost zero MO-induced "dipoles".

Discussion
From the observation of the resulting E z components, one can conclude that the application of a magnetic field has three main effects: First, it generates a "dipolar" component along the ydirection not present without magnetic field; second, mediated by the interaction between the disks, the bottom disk, which has a ferromagnetic component, induces a "dipole" along the ydirection in the upper disk, in other words it generates a MO activity in a disk which has no ferromagnetic nature; third, there is a spectral region were both "dipoles" vanish almost completely, and therefore, the MO activity is strongly reduced.
To understand the physical origin of this phenomenology, let us consider the field distributions obtained for the intermediate energy region shown in the left panel of Fig. 4. In that region, the intensity of the field at the bottom disk clearly decreases with respect to the other two energies considered.This suggests that the electric field induced by the upper disk destructively interferes with the incoming electric field.This can be understood from the simple model based on two interacting point dipoles.In this model the different components of the dipoles can be expressed as: where ( ) is the amplitude of the incident electric field at point dipole p 1 , δ is the phase difference of the incident field at the two dipole positions, α i,jk are the corresponding components of the polarizability tensor of each dipole, and  is the field propagator (the complete derivation of these expressions is given in Appendix B).
As it can be seen, there is a common term [ ] for both x-and y-components of the dipole with MO activity (p 2 ) and the y component of the dipole with no MO activity (p 1 ).This term is related to the total electric field along the x-direction at the MO active dipole, and results from the addition of the incident field and the field generated at this point by the dipole with no MO activity.If this term vanishes the three components die out, and therefore the MO activity is suppressed.On the other hand, as the x component of the dipole with no MO activity does not depend on this term, it is not cancelled being the only one that contributes to the optical response of the whole system.From the point of view of the mechanical model this situation corresponds to the transition region between the in-phase and out-of-phase oscillations of the masses: the charged mass decreases its oscillation and therefore the Lorentz force, proportional to the velocity, is strongly attenuated.As a consequence, the movement of this mass along the y direction and the induced y movement of the uncharged mass are also reduced (see also Fig. 6 and Media2 in Appendix A).This phenomenon can be viewed as the MO counterpart of the electromagnetic induced transparency.For a specific spectral region, the EM interaction between the two disks gives rise to a situation where the total EM field at the MO active element is minimized, therefore strongly reducing the MO activity of the whole system (see Eq. 1).This is in Fig. 5(a), where we present the experimental MO activities for two different situations corresponding to very weak interaction (the 50 nm SiO 2 structure) and a situation where the two disks clearly interact (the 20 nm SiO 2 structure).For the very weak interacting system only a broad peak is observed, corresponding to the MO activity of the magnetoplasmonic disk.However, for the system where the two metallic disks interact, a narrow dip with a clear Fano-like shape is obtained.The difference between both spectra is also shown (inset), presenting a narrow peak.Fano resonances are observed whenever a broad and a narrow state interfere destructively, giving rise to a sharp spectral feature.In our specific case, the broad state corresponds to the EM field at the MO disk in the absence of interaction, whereas the narrow one is due to the electromagnetic field induced by the non-MO disk at the MO disk.Due to its intrinsic interaction origin, this Fano shape strongly depends on the SiO 2 spacer thickness.As it can be seen in Fig. 3, for thick enough SiO 2 spacers, the MO activity is simply that of an individual disk (the bottom disk) and, as the spacer thickness is gradually reduced, the interaction between the disks becomes more relevant and, especially for intermediate spacer thickness, the Fano-like shape is clear.This interaction dependence appears explicitly in the common term of Eq. 1 as the field propagator  depends on the distance between the two dipoles.In Fig. 5(b) we also present the calculated MO activity using a transfer matrix formalism or a simple dipolar model (inset) for two cases corresponding to a structure without the upper disk (no interaction) and a structure with both disks separated by a 20 nm thick SiO 2 layer.As it can be observed the calculated results follow the same trends as the experimental ones.

Summary
In conclusion, MO activity in non-MO elements can be induced via electromagnetic interaction with MO active elements.This effect has been shown in magnetoplasmonic resonators consisting of a Au plasmonic disk separated by a SiO 2 layer from a MO active magnetoplasmonic (Au/Co) disk.When a magnetic field is applied perpendicular to the disks, it produces an electric dipole along the perpendicular direction of the electric field of the incoming light in the MO active disk.This electric dipole induces a MO activity in the nonmagnetic Au disk.Additionally, we have shown the existence of an electromagnetically induced magneto-optical transparency in these systems.This is due to the interference at the MO disk between the incidence electromagnetic field and that generated by the non-MO one, leading to a Fano-like spectral dependence of the MO activity and to a strong reduction of the MO activity.This may find applications in sensing architectures based on MO-Fano resonances.Optical Fano resonances have already demonstrated their potential applicability [45][46][47].Additionally, plasmon excitation enhancement of the MO activity has also shown improvement of the performance of standard plasmonic sensors [48,49].Combining both effects may allow the generation of novel sensing platforms based on MO-Fano resonances.Since the observed phenomenon is driven by the resonance of the upper disk, modifications of the dielectric environment may strongly affect it and, as a consequence, the overall observed Fano resonance characteristics.

Appendix A -Spring-mass model
A physically simple and intuitive model for studying interacting, resonant systems is the damped driven coupled oscillators.This model has been already employed to give a simple picture for systems with two interacting plasmonic modes [9,11,50,51].Going one step further, we can introduce in this mechanical model effects like those of a MO activity in the materials by considering charged masses in a static magnetic field.The presence of the magnetic field and a charged element give rise to the appearance of a Lorentz force l F q B = × r  .In the case considered here, the static magnetic field is applied perpendicularly to the driving, external force, therefore the Lorentz force will be perpendicular to both an thus define a plane of movement.Let's say that the external force is in the x-direction and the magnetic field in the z-direction, then the Lorentz force will be in the y-direction and the plane of movement of the masses will be the x-y plane.
The equations of motion of this model read: ) where the particle positions i r are restricted to the x-y plane and referred to its equilibrium state, i k are spring constants, 12 k is the interaction between the particles, i m are their masses and, more importantly, the damping terms i  are tensors: arising from the Lorentz force: Aceptado para su publicación en Optics Express DOI: 10.1364/OE.21.027356 It is easy to see that in the absence of a magnetic field, and considering only the friction force in a homogeneous medium,  can be regarded as a scalar term (diagonal with equal eigenvalues) b.Now, without losing generality, we further assume that 1 2 m m m = ≡ , and assume that the external driving force has an harmonic temporal dependence i t e ω − .Thus, using the following definitions: the equations of motion in frequency space become: where i  is the corresponding Fourier component of the external force, / i i F m ≡  .In matrix form, this reads as: , where  is the 2x2 unit matrix, 2 σ is the 2x2 Pauli matrix ( 2 0 0 In the case we are interested in, only one of the particles has MO activity, which is the same as to say that either ,1 c ω or ,2 c ω is zero (we have chosen ,1 c ω to be zero).In order to mimic the photonic case, the external driving force is also special, and is applied to both particles along the x-direction only.The solution is then: where In Fig. 6 (Media 2) we show the results obtained from Eq. ( 8), using Aceptado para su publicación en Optics Express DOI: 10.1364/OE.21.027356Notice that the above exposed problem to solve can be generally written as: Therefore the solution reads as: where ( ) . This in the end gives rise to:

Appendix B -Point dipole model
In the previous appendix we have described with detail the case of the classical oscillator model for charged spheres and in presence of a magnetic field.The problem we face is of optical nature, thus the model must be different, but we will use some of the equations from above.The simplest model one can use is to substitute the disks by point dipoles [see Fig. y-amplitude (F x /ω 2 1(b)].These dipoles must reflect the characteristics introduced by the geometrical aspect ratio, thus the polarizability is given by that of an oblate spheroid (see e.g.Chapter 5 in Ref. [52]), and we take into account that the dielectric tensor is not diagonal, and consider radiative corrections [53], ( ) where V is the volume of the spheroid,  is a diagonal matrix that takes care of geometrical aspect, and the environment is vacuum.
It is convenient to work with dimensionless quantities The interaction between the two dipoles is given by the Green tensor, which in its dimensionless form is: where R is the vector connecting r and r'.
With those definitions the electric field incident on a particle of an ensemble of N particles is given by ( ) 0, , .
and the Green tensor that describes the interaction of the two particles (placed along z axis) is diagonal given by: ( ) ( ) ( ) ( ) being d the distance between the particles.One has to take into account that with the geometry described above the z-direction plays no role and can be excluded.The equations to deal with are: where i α is the 2x2 x-y part of the tensor described in Eq. ( 15), and δ = -kd is the phase difference on the incoming wave due to the separation distance d (normally 1 δ  ).
Aceptado para su publicación en Optics Express DOI: 10.1364/OE.21.027356 With that it is easy to see that we end with the following set of equations: which has the same structure as Eq. ( 9) for only one MO dipole (dipole 2).The solution is now straightforward: As it can be seen in Fig. 5, the agreement of the point dipole model with the experimental results and with more elaborated models is very good.Fig. 7 shows the results obtained from Eq. ( 19) for intermediate interaction.It is worth noticing that for the interaction considered the dipole in the y-direction is two orders of magnitude smaller than that in the x-direction, following the results obtained for the fields using the rigorous calculation [Fig.4].Also the relative intensities of the dipoles follow that of the fields.

Fig. 1 .
Fig. 1.(a) Spring model representing a two coupled masses system excited by an harmonic force, F(t), along x axis.Left side: uncharged masses.Right side-(Media 1): one of the masses (blue) is charged (q) and a static magnetic field (B) is applied along the z direction, inducing a Lorentz force, F l (t), along the y direction.The y-movement is transferred to the other mass through the coupling.(b) Two interacting electric dipoles, representing two metallic disks, excited by an incident beam polarized along the x axis.Left side: No disk has magneto-optical activity and the reflected (E r ) and transmitted (E t ) light have the same polarization direction than the incident (E i ) light.Right side: one of the disks (blue) has magneto-optical activity and a static magnetic field (B) applied along the z direction induces a rotation of its electric dipole, which is transferred to the other dipole through the interaction.The rotation modifies the polarization direction of the reflected (E r ) and transmitted (E t ) light.
(d) a top view image of a representative structure, showing the homogeneous distribution of the disks over large areas, is presented.

Fig. 2 .
Fig. 2. (a) Schematic drawing of the nanoresonators composed of a purely plasmonic Au disk and a magnetoplasmonic Au/Co superlattice disk separated by a dielectric spacer.(b) Polar Kerr loop of a characteristic sample.The presence of multiple Co/Au interfaces reduces the value of the magnetic field needed to saturate the nanodisks in the direction perpendicular to the sample plane.(c), (d) Cross section and planar view SEM pictures, respectively, of a representative sample.The images show the homogeneous and random distribution of nanoresonators and their truncated conical shape and internal structure.

Fig. 3 .
Fig. 3. (a) Extinction and (b) MO activity spectra of the nanoresonators as a function of SiO 2 thickness.The different dashed horizontal lines indicate the zero value for each spectrum above the line.The 3D graphs on top of each panel show the results obtained from theoretical calculations of the same structures.The blue, green and red diamonds appearing in the experimental MO response for 20 nm SiO 2 correspond to the spectral positions where the E z field distributions are calculated [see Fig. 4].

Fig. 4 .
Fig. 4. (a) Calculated near field intensity of the E z component for a nanoresonator with 20 nm SiO spacer in two planes above the top disk and below the bottom one, with the incident field polarized along the x direction and in the absence of an external magnetic field.This distribution reflects the excitation of two dipoles along the x direction.The insets indicate the corresponding charges and dipole orientations (indicated by the black arrows) according to the E z distribution for the different cases.The red (blue) arrows in the insets represent the positive (negative) values of the Ez field component.(b)Difference of the E z components for magnetic saturation along opposite directions in the same planes and for the same structure.This difference accounts for the effect of the applied magnetic field: The appearance of a dipole along the y direction for both the magnetoplasmonic (intrinsic dipole) and the plasmonic (induced dipole) disks.In both cases, the components for three different wavelengths labelled as diamonds in Fig.3are shown.

Fig. 5 .
Fig. 5. Fano resonance in the magneto-optical activity.(a) Experimental spectra of the MO activity for the structures with 20 nm and 50 nm SiO 2 spacer, corresponding to interacting and non interacting situations, respectively.For the interacting situation a clear Fano resonance shape is observed.The inset shows the difference between the two spectra, resulting in a narrow peak.(b) Theoretical MO activity spectra obtained for a structure with 20nm SiO 2 spacer and for the same structure but without the upper metallic disk.The inset shows the MO spectra obtained using the simple point dipole model.

Fig. 6 .
Fig. 6.Oscillation amplitude as a function of the frequency of masses 1 (uncharged) and 2 (charged) along x-(left panel) and y-(right panel) direction by using a simple spring-mass model.Inset: Frame of a video showing schematically the evolution with the frequency of the masses' oscillation for frequency values around the minimum in x-amplitude for the charged mass 2 (in the video, the y-amplitude has been multiplied by a factor 10, and the frequency has been scaled as ω vid = (1+(ω-0.93)×50)ω 2 for clarity in the visualization) (Media 2).

2 |
p  the usual electric dipole.It is important to notice that with this definition the optical theorem for a scalar absorptionless particle reads as polarized along the x-axis and its wavevector is parallel to the z-axis (short axis of the particles), the polarization of the particles lies in the x-y plane,

Fig. 7 .
Fig. 7. Wavelength dependence of the normalized dipole magnitudes along the x-(left panel) and y-(right panel) directions obtained using a point dipole model.