Abstract
We present a systematic, theoretical investigation of the polar magneto-optical (MO) Kerr effects of a single Ni nanorod in the Mie regime. The MO Kerr rotation, ellipticity, amplitude ratio, and phase shift are calculated as a function of the length and width of the nanorod. The electric field amplitude ratio of the MO Kerr effect is locally maximized when the nanorod supports a plasmonic resonance in the polarization state orthogonal to the incident light. The plasmonic resonances directly induced by the incident light do not enhance the amplitude ratio. In the Mie regime, multiple local maxima of the MO Kerr activity are supported by the resonant modes with different modal characteristics. From the viewpoint of first-order perturbation analysis, the spatial overlap between the incident-light-induced electric field and the Green function determines the local maxima.
© 2016 Optical Society of America
1. Introduction
Magneto-optic (MO) effects make it possible to engineer the polarization and phase of electromagnetic waves by spin–orbit interactions in magnetized media [1] and have been used for magnetic detection [2], optical isolation [3], information recording [4], and optical modulation [5]. On the basis of their plasmonic resonances, metal nanostructures have recently been employed to modify or boost the MO effects [6,7]. Periodic metal gratings combined with a thin magnetic film yielded enhancement of the Faraday rotation [8], MO Kerr effect [9], and MO transparency [10]. Metal–ferromagnet–metal waveguides can modulate the wave vector of guided surface plasmon polaritons (SPPs) and can be used to construct an active SPP interferometer [11]. In such studies, magnetic media have been employed just as an assistant to add MO functions to plasmonic nanostructures. On the other hand, magnetic metal materials such as Ni, Fe, and Co can support plasmonic resonances and MO activity simultaneously, and intermediate their direct interplay [12–14]. Recently, the magneto-plasmonic responses of Ni nanodisks have received wide attention [15,16], and their periodic or random arrangements have been applied to demonstrate resonant MO Kerr effects [17] and phase-sensitive molecule-level sensing [18].
To obtain general guidelines for engineering the MO effects of magnetic metal nanostructures, a systematic understanding of the coupling between the MO activity and the plasmonic resonances is required. However, most theoretical analyses are based on simplified models using quasi-static analysis or the effective medium theory [16–21]. The quasi-static model assumes uniformly distributed internal fields and fails to fully include the modal characteristics of nanostructures with a nonspherical, complex geometry or a size exceeding the Rayleigh criterion (~λ/10). Modified analytical methods, such as the use of the polarizability tensor considering the retarded depolarization to predict the MO properties of ellipsoidal particles [22], have only been suggested. In this work, we systematically investigate the polar MO Kerr effect of a single Ni nanorod in the Mie regime depending on its geometry using the finite-difference time-domain (FDTD) method, and present a universal explanation of the resonant MO Kerr effects from the viewpoint of first-order perturbation analysis.
2. Schematics
In the presence of a magnetic field, either external or internal, two orthogonal electric polarizations in the magnetized medium are coupled by spin–orbit interactions. Because the MO activity changes the permittivity, typically by on the order of 1% or less [23], the MO permittivity can be approximated by adding the perturbation of an off-diagonal tensor to the nonmagnetic permittivity . When the direction of the magnetization is along the z axis, the permittivity of the Ni medium can be formulated as a Hermitian tensor,
Here and Q are the nonmagnetic permittivity and Voigt parameter, respectively. In this case, first-order perturbation theory can successfully model the MO effects of the Ni nanostructure [24]. The first-order MO perturbation in the steady state with a frequency is obtained by is the electric field without the MO activity and V is the volume of the Ni medium. is the Green function of the electromagnetic wave equation involving only the nonmagnetic element of the permittivity: . The Green function acts as the propagator and gives the probability amplitude of coupling of the initial-state mode induced by the incident light to the final-state mode with an orthogonal polarization state. The coupling probability amplitude increases with increasing local density of optical states (LDOS) of the final-state mode. For a spherical nanoparticle, the Green function and first-order perturbation integral can be simply obtained in terms of the spherical harmonics and spherical Bessel functions. We employed Eq. (2) for the qualitative interpretation of the simulation results from the viewpoint of first-order perturbation analysis. For a structure with complicated boundary conditions, numerical simulations are required to investigate the MO activity, but first-order perturbation theory is still useful for interpreting numerical simulation results based on modal analysis.Figures 1(a) and 1(b) show schematics of the MO activity in a single Ni nanorod under a normally incident planewave in the quasi-static and Mie regimes, respectively. The incident light is linearly polarized along the y axis, and a static external magnetic field is applied along the z axis. For simplicity, the surrounding medium is considered to be air. In the quasi-static regime, electric fields are uniformly induced inside the Ni nanorod, and the MO perturbations to the electric fields, ΔEx, must also be uniform. On the other hand, in the Mie regime, the electric fields and MO perturbations are not uniform due to spatial retardation effects. The incident light polarized along the nanorod length excites SPPs of the transverse-magnetic (TM) mode, and the MO activity interfaces the TM mode to the transverse-electric (TE) mode. We define the TM and TE polarization states according to the length direction of the Ni nanorod. The two ends of the nanorod reflect propagating SPPs and support longitudinal Fabry–Pérot (FP) resonances of the TM and TE modes. Because of their different dispersion relations, the TM and TE modes exhibit FP resonances at different nanorod lengths for a given wavelength of light.
In the FDTD simulations, we suppose that, to support the polar MO Kerr effects, the magnetization of the Ni medium is saturated along the + z direction and a 660-nm-wavelength planewave is normally incident from the top. At the wavelength of 660 nm, the nonmagnetic permittivity of Ni is −10.726 + 15.842i [25], and the Voigt parameter Q in the saturation condition is 0.00588 + 0.0039i [26]. To apply the Hermitian permittivity tensor, we used the grid attribute technique of the FDTD calculator (Lumerical Solutions, Inc.). In the magnetized Ni medium, the electric fields are first calculated in the unitarily transformed coordinates by diagonalizing the Hermitian permittivity tensor and then converted to those in the Cartesian coordinates again. Because the field components in the Yee grid are spatially offset by a half-step and the unitary coordinate transformation is valid only in the magnetized medium, a crucial error occurs at the interface between the magnetized and nonmagnetized media. To avoid this interface error, we suppose that the outermost grids at the boundary of the Ni medium are not magnetized. The artificial effects of the nonmagnetized outer layer can be suppressed sufficiently by reducing the grid size in the Ni medium. We set the spatial grid size inside the Ni medium to 2 nm and gradually increase the grid size in the other region from 2 to 16 nm.
3. Modal properties of nonmagnetized Ni nanorods
We first investigate the modal properties of the nonmagnetized Ni nanorod by calculating the extinction and scattering efficiencies. The extinction (scattering) efficiency is defined as the ratio between the extinction (scattering) cross section and the physical cross section of the nanorod. Figures 2(a) and 2(b) show two-dimensional (2D) maps of the extinction and scattering efficiencies, respectively, as a function of the length and width of the nanorod under a normally incident TM-polarized planewave. The length and width vary from 80 to 1120 nm, and the thickness is fixed at 40 nm. The FP resonances of the guided modes locally maximize the extinction and scattering efficiencies of single nanorod structures [27,28]. For the 80-nm-wide nanorod, the extinction efficiency is locally maximized at lengths of ~200 and ~720 nm by the first- and third-order longitudinal FP resonances of the TM mode, with extinction efficiencies of ~5.35 and ~2.58, respectively. The even-order FP resonances cannot be excited by the normally incident planewave because of their antisymmetric field distributions. The electric field amplitude and surface charge density distributions in Figs. 2(c)–2(f) clearly show the spatial properties of the first- and third-order FP resonances of the TM mode. The electric field and charge density distributions of the first-order resonance have the features of typical dipolar oscillation. As the width of the nanorod decreases into the quasi-static regime, the absorption efficiency becomes larger than the scattering efficiency, and the extinction efficiency gradually increases. As the nanorod width increases, the phase pickup for reflecting SPPs decreases, and FP resonances appear for slightly longer nanorods. In addition, we note that, even though they are not excited by the TM-polarized incident light, the longitudinal FP resonances of the TE mode also exist and play an important role in the MO Kerr effect, which will be discussed later.
4. Polar magneto-optic Kerr effects of magnetized Ni nanorods
To investigate the polar MO Kerr effect of the Ni nanorod, we calculated the Kerr rotation and Kerr ellipticity of the far-field radiation scattered back to the vertical direction ( + z direction):
The scattered electric field in the far-field region can be obtained by the near-to-far-field transformation [29]. When the phase difference between the Ex and Ey fields, , is small enough, the Kerr rotation and ellipticity can be approximated by the amplitude ratio and the product of the phase difference and amplitude ratio, respectively. Figures 3(a)–3(d) show 2D maps of the Kerr rotation, Kerr ellipticity, amplitude ratio, and phase difference as a function of the length and width of the nanorod. Showing support for four local maxima, all the 2D maps differ from that of the extinction efficiency in Fig. 2(a). The width and length supporting the local maxima of the amplitude ratio absolutely do not correlate with those supporting the maxima of the extinction efficiency for TM-polarized incidence. To enhance the amplitude ratio and allow a large Kerr rotation, a high LDOS in the TE polarization is required; the local maxima of the amplitude ratio depend entirely on the resonances of the TE mode.The MO perturbation ΔEx profile provides a better understanding of the MO Kerr effect of the Ni nanorod. Figures 3(e)–3(h) show the ΔEx profiles in the xy plane for the four selected nanorods that support the local maxima of the amplitude ratio in Fig. 3(c). The one- or two-dimensional standing wave pattern in the ΔEx profiles verifies the resonance of the plasmonic TE mode, the polarization state of which is orthogonal to the incident-light-induced TM mode. The TE mode resonances associated with each local maximum of the amplitude ratio/Kerr rotation have different modal characteristics. First, the enhancement of the Kerr effect for short nanorods with a length of 80 nm originates from the transverse FP resonances of the TE mode along the width direction of the nanorod. Rotated by 90°, the transverse FP resonances of the TE mode are identical to the longitudinal FP resonances of the TM mode. As shown in Figs. 3(e) and 3(f), the 200-nm- and 720-nm-wide nanorods support the first- and third-order FP resonances of SPPs propagating along the x axis, respectively. On the other hand, the local maximum of the amplitude ratio in the 680-nm-long and 200-nm-wide nanorod originates from the longitudinal FP resonance of the TE mode, the modal characteristics of which are completely different from those of the longitudinal (transverse) FP resonance of the TM (TE) mode.
The ΔEx profile in Fig. 3(g) clearly shows the longitudinal FP resonance of the TE mode propagating along the length direction. Three equally spaced lobes exhibit the typical characteristics of the third-order longitudinal FP resonance. We calculated the propagation wavelength of the TE-polarized SPPs and the reflection pickup at the end facet for the Ni nanorod with a 40-nm-thick and 200-nm-wide cross section. From the calculated propagation wavelength and reflection phase pickup of ~650 nm and ~166°, respectively, the nanorod length supporting the third-order longitudinal FP resonance is estimated to be ~676 nm, in excellent agreement with the length maximizing the Kerr effect. In contrast to the TM mode resonances, the TE mode resonances cannot be effectively excited by the normally incident planewave because of their low radiation efficiencies in the vertical direction. However, the internal electric fields induced inside the Ni nanorod cause the MO electric polarizations acting as point sources to excite the TE mode resonances. Two-dimensional SPP resonances can also be allowed with increasing length and width of the nanorod. Figure 3(h) shows the 2D standing wave pattern of the ΔEx field in the 680-nm-long and 720-nm-wide nanorod. The locally maximized MO Kerr rotation values of the third-order longitudinal FP resonance of the TE mode and the 2D SPP resonance are ~2.20 and ~1.56 mrad, respectively. The Kerr rotations maximized by the first- and third-order transverse FP resonances of the TE mode are ~3.18 and ~1.48 mrad, respectively.
To support a local maximum of the first-order perturbation integral of the MO Kerr effect, it is necessary to optimize the spatial overlap and phase matching between the incident-light-induced electric field and the Green function, in particular, in the resonance condition of the final-state mode, where the Green function has large amplitudes and is accompanied by a high LDOS. In the numerical simulations, the optimization condition can be identified by examining whether the amplitude and phase distributions of the incident-light-induced field inside the MO medium can promote the resonance condition of the final-state mode with an orthogonal polarization (Fig. 4(a)). This mechanism is analogous to that of the Purcell effect, in which the emission rate enhancement is proportional to the LDOS at the position of the emitter. In Fig. 4, we calculate the amplitude and relative phase distributions of the y component of the incident-light-induced field Ey inside the Ni nanorod and the x component of the resulting MO perturbation ΔEx. Here, the lengths of the nanorods are 480, 680, and 760 nm, respectively, and the widths are fixed at 200 nm.
First, we focus on the 680-nm-long Ni nanorod showing the local maximum of the amplitude ratio of the MO Kerr effect. The 680-nm-long nanorod does not support a resonance of the TM mode, whereas the 480-nm- and 760-nm-long nanorods support the second- and third-order longitudinal FP resonances of the TM mode, respectively. However, the 680-nm-long nanorod supports the third-order longitudinal FP resonance of the TE mode, which is the final-state mode in the first-order perturbation integral of the MO Kerr effect. As shown in Fig. 4(c), the incident-light-induced Ey field has three evenly spaced amplitude antinodes along the longitudinal direction. Neighboring antinodes with comparable magnitudes oscillate with opposite phases to each other. In the transverse direction, the amplitude antinodes of the Ey field are positioned at the sides of the nanorod, where the Ex field of the TE mode is maximized. This profile of the Ey field causes the MO perturbation ΔEx to seamlessly couple to the third-order FP resonance of the TE mode and thus the electric field amplitude ratio of the MO Kerr effect is locally maximized.
On the other hand, the amplitude and phase distributions of the Ey field of the 760-nm- and 480-nm-long nanorods are far from those of a standing wave pattern (Figs. 4(b) and 4(d)). Despite its length, the 760-nm-long nanorod has only two antinodes of the Ey field along the longitudinal direction. Further, the antinodes oscillate in the same phase. The Ey field profile of the 480-nm-long nanorod has a single antinode in the middle along the longitudinal direction, but the nanorod is too long to support the first-order FP resonance of the TE mode. The ΔEx profiles also show features of the transitional stages between the FP resonance conditions. Indeed, the 480-nm- and 760-nm-long nanorods cannot support a high LDOS in the TE mode or enhance MO Kerr activity, even though they support the FP resonances for TM-polarized incidence.
5. Summary
We presented a systematic study of the coupling between the polar MO Kerr activity and the plasmonic resonances of a single Ni nanorod in the Mie regime. The electric field amplitude ratio of the MO Kerr effect is locally maximized when the Ni nanorod supports a plasmonic resonance in the polarization state orthogonal to the incident light. The plasmonic resonances directly induced by the incident light do not contribute to the maximization of the amplitude ratio. From the viewpoint of first-order perturbation analysis, the spatial overlap and phase matching between the incident-light-induced electric field and the Green function determine the MO Kerr effect. Our understanding of the MO effects of optical nanostructures will be useful for various applications such as MO sensors, optical logic devices, and MO metamaterials.
Funding
National Research Foundation of Korea (NRF) (2013R1A2A2A01014224, 2014M3C1A3052537, and 2014M3A6B3063709).
References and Links
1. P. N. Argyres, “Theory of the Faraday and Kerr Effects in Ferromagnetics,” Phys. Rev. 97(2), 334–345 (1955). [CrossRef]
2. Z. Q. Qiu and S. D. Bader, “Surface magneto-optic Kerr effect,” Rev. Sci. Instrum. 71(3), 1243 (2000). [CrossRef]
3. L. Bi, J. Hu, P. Jiang, D. H. Kim, G. F. Dionne, L. C. Kimerling, and C. A. Ross, “On-chip optical isolation in monolithically integrated non-reciprocal optical resonators,” Nat. Photonics 5(12), 758–762 (2011). [CrossRef]
4. S. Tsunashima, “Magneto-optical recording,” J. Phys. D. 34(17), R87–R102 (2001). [CrossRef]
5. L. A. Williamson, Y. H. Chen, and J. J. Longdell, “Magneto-Optic Modulator with Unit Quantum Efficiency,” Phys. Rev. Lett. 113(20), 203601 (2014). [CrossRef] [PubMed]
6. 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). [CrossRef] [PubMed]
7. 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). [CrossRef]
8. 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). [CrossRef] [PubMed]
9. 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). [CrossRef] [PubMed]
10. V. I. Belotelov, L. E. Kreilkamp, I. A. Akimov, A. N. Kalish, D. A. Bykov, S. Kasture, V. J. Yallapragada, A. Venu Gopal, A. M. Grishin, S. I. Khartsev, M. Nur-E-Alam, M. Vasiliev, L. L. Doskolovich, D. R. Yakovlev, K. Alameh, A. K. Zvezdin, and M. Bayer, “Plasmon-mediated magneto-optical transparency,” Nat. Commun. 4, 2128 (2013). [CrossRef] [PubMed]
11. V. V. Temnov, G. Armelles, U. Woggon, D. Guzatov, A. Cebollada, A. García-Martín, J. M. García-Martín, T. Thomay, A. Leitenstorfer, and R. Bratschitsch, “Active magneto-plasmonics in hybrid metal–ferromagnet structures,” Nat. Photonics 4(2), 107–111 (2010). [CrossRef]
12. E. Th. Papaioannou, V. Kapaklis, E. Melander, B. Hjörvarsson, S. D. Pappas, P. Patoka, M. Giersig, P. Fumagalli, A. Garcia-Martin, and G. Ctistis, “Surface plasmons and magneto-optic activity in hexagonal Ni anti-dot arrays,” Opt. Express 19(24), 23867–23877 (2011). [CrossRef] [PubMed]
13. E. Th. Papaioannou, V. Kapaklis, P. Patoka, M. Giersig, P. Fumagalli, A. García-Martín, E. Ferreiro-Vila, and G. Ctistis, “Magneto-optic enhancement and magnetic properties in Fe antidot films with hexagonal symmetry,” Phys. Rev. B 81(5), 054424 (2010). [CrossRef]
14. G. Ctistis, E. Papaioannou, P. Patoka, J. Gutek, P. Fumagalli, and M. Giersig, “Optical and Magnetic Properties of Hexagonal Arrays of Subwavelength Holes in Optically Thin Cobalt Films,” Nano Lett. 9(1), 1–6 (2009). [CrossRef] [PubMed]
15. V. Bonanni, S. Bonetti, T. Pakizeh, Z. Pirzadeh, J. Chen, J. Nogués, P. Vavassori, R. Hillenbrand, J. Åkerman, and A. Dmitriev, “Designer magnetoplasmonics with nickel nanoferromagnets,” Nano Lett. 11(12), 5333–5338 (2011). [CrossRef] [PubMed]
16. N. Maccaferri, A. Berger, S. Bonetti, V. Bonanni, M. Kataja, Q. H. Qin, S. van Dijken, Z. Pirzadeh, A. Dmitriev, J. Nogués, J. Åkerman, and P. Vavassori, “Tuning the magneto-optical response of nanosize ferromagnetic Ni disks using the phase of localized plasmons,” Phys. Rev. Lett. 111(16), 167401 (2013). [CrossRef] [PubMed]
17. M. Kataja, T. K. Hakala, A. Julku, M. J. Huttunen, S. van Dijken, and P. Törmä, “Surface lattice resonances and magneto-optical response in magnetic nanoparticle arrays,” Nat. Commun. 6, 7072 (2015). [CrossRef] [PubMed]
18. N. Maccaferri, K. E. Gregorczyk, T. V. A. G. de Oliveira, M. Kataja, S. van Dijken, Z. Pirzadeh, A. Dmitriev, J. Åkerman, M. Knez, and P. Vavassori, “Ultrasensitive and label-free molecular-level detection enabled by light phase control in magnetoplasmonic nanoantennas,” Nat. Commun. 6, 6150 (2015). [CrossRef] [PubMed]
19. N. Maccaferri, M. Kataja, V. Bonanni, S. Bonetti, Z. Pirzadeh, A. Dmitriev, S. Dijken, J. Åkerman, and P. Vavassori, “Effects of a non-absorbing substrate on the magneto-optical Kerr response of plasmonic ferromagnetic nanodisks,” Phys. Status Solidi., A Appl. Mater. Sci. 211(5), 1067–1075 (2014). [CrossRef]
20. K. Lodewijks, N. Maccaferri, T. Pakizeh, R. K. Dumas, I. Zubritskaya, J. Åkerman, P. Vavassori, and A. Dmitriev, “Magnetoplasmonic design rules for active magneto-optics,” Nano Lett. 14(12), 7207–7214 (2014). [CrossRef] [PubMed]
21. I. Zubritskaya, K. Lodewijks, N. Maccaferri, A. Mekonnen, R. K. Dumas, J. Åkerman, P. Vavassori, and A. Dmitriev, “Active magnetoplasmonic ruler,” Nano Lett. 15(5), 3204–3211 (2015). [CrossRef] [PubMed]
22. N. Maccaferri, J. B. González-Díaz, S. Bonetti, A. Berger, M. Kataja, S. van Dijken, J. Nogués, V. Bonanni, Z. Pirzadeh, A. Dmitriev, J. Åkerman, and P. Vavassori, “Polarizability and magnetoplasmonic properties of magnetic general nanoellipsoids,” Opt. Express 21(8), 9875–9889 (2013). [CrossRef] [PubMed]
23. A. K. Zvezdin and V. A. Kotov, Modern Magneto-optics and Magneto-optical Materials (informa, 1997).
24. R.-J. Tarento, K.-H. Bennemann, P. Joyes, and J. Van de Walle, “Mie scattering of magnetic spheres,” Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 69(2), 026606 (2004). [CrossRef] [PubMed]
25. E. D. Palik, Handbook of Optical Constants of Solids I (Academic, 1985).
26. Š. Višňovský, V. Pařízek, M. Nývlt, P. Kielar, V. Prosser, and R. Krishnan, “Magneto-optical Kerr spectra of nickel,” J. Magn. Magn. Mater. 127(1–2), 135–139 (1993). [CrossRef]
27. E. S. Barnard, J. S. White, A. Chandran, and M. L. Brongersma, “Spectral properties of plasmonic resonator antennas,” Opt. Express 16(21), 16529–16537 (2008). [CrossRef] [PubMed]
28. H.-S. Ee, J.-H. Kang, M. L. Brongersma, and M.-K. Seo, “Shape-dependent light scattering properties of subwavelength silicon nanoblocks,” Nano Lett. 15(3), 1759–1765 (2015). [CrossRef] [PubMed]
29. K. C. Y. Huang, Y. C. Jun, M.-K. Seo, and M. L. Brongersma, “Power flow from a dipole emitter near an optical antenna,” Opt. Express 19(20), 19084–19092 (2011). [CrossRef] [PubMed]