Faraday-Rotation in Iron Garnet Films Beyond Elemental Substitutions
Miguel Levy, Olga V. Borovkova, Colin Sheidler, Brandon Blasiola, Dolendra Karki, François Jomard, Mikhail A. Kozhaev, Elena Popova, Niels Keller, Vladimir I. Belotelov
1 Faraday-Rotation in Iron Garnet Films Beyond Elemental Substitutions
Miguel Levy , Olga V. Borovkova , Colin Sheidler , Brandon Blasiola , Dolendra Karki , François Jomard , Mikhail A. Kozhaev , Elena Popova , Niels Keller , Vladimir I. Belotelov Physics Department, Michigan Technological University, Houghton, Michigan 49931, USA Russian Quantum Center, Skolkovo, Moscow Region 143025, Russia Groupe d’Etude de la Matière Condensée (GEMaC),
CNRS-UVSQ, Université Paris-Saclay, 78035 Versailles, France Prokhorov General Physics Institute RAS, Moscow 119991, Russia Lomonosov Moscow State University, Moscow 119991, Russia
Abstract
In previous decades significant efforts have been devoted to increasing the magneto-optical efficiency of iron garnet materials for the miniaturization of nonreciprocal devices such as isolators and circulators. Elemental substitutions or proper nano-structuring to benefit from optical resonances have been pursued. However, all these approaches still require film thicknesses of at least several tens of microns to deliver useful device applications, and suffer from narrow bandwidths in the case of optical resonance effects. This Letter reports on a newly discovered enhancement of the Faraday Effect observed experimentally in nanoscale bismuth-substituted iron garnet films. It is shown here that this enhancement is not due to elemental substitution or compositional variations. Nor is it due to photon trapping or resonance effects. Comprehensive experimental and theoretical analysis of the Faraday rotation reveals a dramatic seven-fold amplification in the magneto-optic gyrotropy within only 2 nm of the air-surface interface, corresponding to just a couple of atomic monolayers as a result of symmetry-breaking at the air-film interface. This finding opens up an avenue to the application of monolayer magnetic garnets for the control of light.
Rare-earth iron garnet films are technologically important since they are among the main materials used in magneto-photonics, ultrafast magnetism and spintronics [1, 2]. In particular, they are used extensively in optical circuits for the fabrication of nonreciprocal optical devices. Isolators and circulators, operating as optical filters rely on the non-reciprocal functionality of iron garnets to protect laser sources from reflected downstream light that might otherwise enter the cavity and destabilize the signal [2-4]. Magneto-optical light modulators and magnetic sensors utilize iron garnet films for efficient magneto-optical response and outstanding performance [3-5]. On the other hand, iron garnet thin films are being successfully used in magnonics in experiments on the spin-orbit torque transfer and nonlinear spin-waves dynamics [6-9]. Since the development of yttrium iron garnet in the 1950s, subsequent bismuth and cerium substitutions in the dodecahedral sub-lattice of these materials were found to enhance their Faraday rotation response [1, 5, 10]. Such substitutions amplify the spin-orbit splitting of the 3d Fe excited states, impacting the refractive index difference between photons of opposite helicities [1]. The Faraday rotation at wavelengths longer than the absorption band is affected as a result. These discoveries impacted the efficiency of nonreciprocal devices, and made it possible to reduce their footprint and hence help drive their miniaturization. However, no other elemental substitutions with similar effect have been found since the 1990s. The need for on-chip isolator integration for the manufacturing of optical circuits, especially for telecommunication purposes, has led to extensive research on and prototyping of these devices over the years [11-15]. Smart (or engineered/designed) magnetophotonic crystal and magnetoplasmonic structures incorporating bismuth-substituted iron-garnet films were demonstrated [16-19]. However, better and stronger non-reciprocal efficiency than at present is still needed by industry to further miniaturize these systems down to nanometer size. This Letter reports on newly discovered enhancement effects observed experimentally in ultrathin rare-earth substituted iron garnet films [20-22]. It is shown here that these effects are not due to elemental substitution or compositional variations. Nor are they due to photon trapping or resonance effects. Theoretical modelling presented in this Letter and verified experimentally, shows that surface reflections do impact the magneto-optic response of these films but cannot account for the large Faraday rotation enhancement observed below 50nm-film-thickness. There are four physical mechanisms that are capable of producing enhancement in the Faraday rotation in magnetic garnets: a) bismuth or cerium substitutions into the dodecahedral sub-lattice of these crystals, b) optical resonances due to reflections, c) surface symmetry-breaking, and d) quantum confinement. This Letter provides strong theoretical and experimental evidence that the observed effects are not due to rare-earth substitutions nor classical electrodynamic interference, thus pointing to surface symmetry breaking effects as the alternative mechanism of significance in the physical response of ultra-thin iron garnet magneto-optics. Modelling of these experimental results is consistent with a surface enlargement in the magneto-optic gyrotropy parameter responsible for the Faraday effect in these films by a factor of seven within only 2nm of the surface . Two different sets of liquid-phase-epitaxy (LPE)-grown films with Bi Gd Lu Fe O , and Bi Gd Lu Fe Ga O per-formula-unit (pfu) were studied. These films were fabricated on gadolinium gallium garnet (GGG) (100)-oriented substrates. The films for any set of measurements were all taken from the same wafer and wet-etch-thinned-down sequentially in ortho-phosphoric acid, thus avoiding possible composition differences due to growth conditions. Light from a CW laser source was used to probe the Faraday rotation at 532nm. Measurements normal to the film surface were conducted in a rotating-polarizer configuration, with an accuracy of 0.001 [3]. Faraday-rotation hysteresis loops were recorded for each sample and the Faraday rotation at saturation was chosen to characterize the response of each film. The paramagnetic signal of the GGG substrate was subtracted out from the overall Faraday rotation signal. The experimental setup is described in [20]. Film thickness was measured via ellipsometry, and Faraday rotation per unit length is plotted in Fig. 1 for both sets of samples as a function of thickness. For the ferrimagnetic films thinner than 50 nm the specific Faraday rotation (SFR) grows dramatically with a decrease in film thickness. One can see that the SFR reaches m in the 20-nm-thick film, almost 2.5 times greater than SFR in the same thick film. For the thin films we observe oscillations due to interference of the emerging waves due to multiple reflections inside the thin film. Fig. 1. Specific Faraday rotation at 532 nm versus film thickness for Bi Gd Lu Fe O (Sample 1) and Bi Gd Lu Fe Ga O (Samples 2 and 3), showing a dramatic increase below 50 nm. Inset: Schematic depiction of the samples under consideration. In this section and in Methods we describe the theoretical model, based on classical electrodynamics, used to analyze the measured Faraday rotation, and the results of this analysis. Linearly polarized light impinges on the magnetic film at normal incidence. The magnetization is directed along the propagation direction, perpendicular to the film surface. Refractive indices and film thickness are as shown in the inset to Fig. 1. The magneto-optic gyration in the ferrimagnetic film is initially assumed to be uniform inside the film and parameterized by the gyrotropy parameter g . Multiple reflections are included in the analysis, giving rise to an oscillatory behavior in FR as a function of film thickness h . The exact expression for Faraday rotation is derived in Methods [see Eq. (M3)]. It has two asymptotes, for h , ultra-thin films, and for h , bulk ferrimagnetic material, with specific Faraday rotations (SFR) given by g if hnSFR gh if hn (1) From (1), the ratio of the SFR at small h with respect to the SFR at large h is, thus, hh SFR nSFR n (2) The refractive indices at wavelength λ = 532nm are given by 𝑛 = 1.980 + 0.003𝑖 in the GGG substrate, and 𝑛 = 2.610 + 0.056𝑖 in the Bi-substituted iron garnet films for all samples. Therefore, at λ = 532nm this ratio is predicted to be about 1.75, but the experimentally measured ratios are always larger than 2. In particular, hh SFRSFR , for Sample 1, and hh SFRSFR , for Samples 2 and 3. Thus, the analytical expression derived from classical electrodynamic cannot account for all the effects operating in this system. Besides the preceding analysis we also performed an electromagnetic modelling based on a rigorous coupled-wave analysis (RCWA) [23, 24]. The result of the corresponding numerical simulation is given in Fig. 2 by the red line (Model 1). One can see that this model does not reproduce the pronounced increase in the SFR obtained experimentally. Therefore, the assumption of constant gyration inside the ferrimagnetic film does not describe the observed abrupt growth in the SFR in ultrathin ferromagnetic films shown in Fig. 1. Fig. 2. (a) Theoretical fits to Sample-2 data based on uniform-gyrotropy-parameter in the film (Model 1), gyrotropy-parameter proportional to Bi content in the film (Model 2), and gyrotropy-parameter-enhancement at the top surface (Model 3). A seven-fold g -value magnification is predicted within 2 nm of the surface (Model 3), and four-fold magnification over the next 2 nm, as compared to bulk g -value. (b) Schematic depiction of the gyrotropy parameter distribution across the film in the three models. BiREIG stands for Bi-rare-earth-substituted iron garnet. This abrupt growth in SFR in ultrathin films might possibly have been explained by composition gradients in the film/substrate interface, such as an increase in Bi content [1, 5, 10]. However, detailed cross-sectional composition analysis, discussed below, shows that there is no such increase in Bi content, and that the only composition gradient is manifested as a decrease in Bi, Fe, and Lu, content in the transient layer (Fig. 3). The cross-sectional film composition was measured by two different methods, via secondary-ion-mass-spectroscopy (SIMS), Fig. 3(a) and cross-sectional scanning transmission electron microscopy (S-TEM), Fig. 3(b) (Methods). Fig. 3. (a) Bi and Ga concentrations in Bi Gd Lu Fe O (Sample 1), arbitrary units, determined by SIMS on a 60nm-thick film. Bi concentration is shown multiplied by factor 10 for the sake of clarity. (b) Elemental concentrations (pfu) for Bi, Lu, Gd, Ga and Fe in Bi Gd Lu Fe Ga O (Samples 2 and 3) determined by cross-sectional S-TEM EDX. Both show uniform concentrations as a function of position above (film) and below (substrate) the transient layer. The zero of horizontal coordinates is set at the film-substrate interface. Assuming a gyration proportional to the Bi-content, comprising also its decline in the transient layer (Fig. 2(b), blue line), the numerical simulation produces an even poorer coincidence with the experimental data than Model 1, constant gyration. This is shown by the blue line (Model 2) in Fig. 2(a). Notice that Model 1 in effect already assumes a significant increase in magneto-optic gyration in the interfacial region, if one considers the decrease in Bi content in the transient layer. Yet the calculated SFR in that case still remains well below the experimental data. To achieve coincidence with the experimental data, we allow the gyrotropy parameter g to experience an increase near the air/film interface. Best fits to the experimental data are produced by a rise in g value taking place over a very thin layer adjacent to the surface, which is h’ = 4nm thick. In Model 3 the gyration evinces a near-surface step-like growth, as depicted in Fig. 2(b). One can see from the Fig. 2 (Model 3) that this assumption provides very good agreement with the SFR, including for sub-50nm-thick ferrimagnetic films. These fits predict that the magneto-optic gyrotropy parameter exhibits a seven-fold amplification in magnitude within 2 nm of the surface, and a four-fold amplification in the next 2 nm, over its bulk value (Methods). To summarise, detailed Faraday rotation data and electrodynamic analysis of ultra-thin Bi-substituted iron garnet films predict a several-fold increase in near-surface magneto-optic gyrotropy within 2 nm of the surface. Electron dispersive X-ray spectroscopy and SIMS confirm that the measured Faraday rotation enhancement in ultra-thin films is not due to compositional non-uniformities near the film-substrate interface, thus pointing to surface symmetry breaking effects as the origin of these results. This dramatic enhancement in the gyrotropy parameter occurs within a few monolayers of the films, as the unit-cell in these iron garnets is ~1.2 nm. These results hold a significant potential for the nanofabrication of nonreciprocal devices, given the technological role of these materials in optical isolator and circulator components. Acknowledgements
This study was supported by Russian Foundation for Basic Research (project no. 16-32-60135 mol_a_dk), and by Foundation for the advancement of theoretical physics BASIS. ML and DK acknowledge support from the Henes Center for Quantum Phenomena. The authors also thank Ashim Chakravarty for valuable discussions and Pinaki Mukherjee for assistance with the cross-sectional EDX S-TEM measurements of film composition. Methods Compositional analysis
SIMS measurements were performed with O bombardment and positive ion detection mode using an IMS-7f (CAMECA) micro-analyzer. The primary beam was rastered over an area of 125 x 125 µm and the secondary ions were collected from the central part of this area (diameter 33 µm). In order to avoid charging effect, the sample was coated with a layer of gold (50 nm) and electron flooding was employed using a normal incidence e-gun. Electron micrographs and Energy dispersive X-ray spectroscopy (EDX) maps were obtained on Bi Gd Lu Fe Ga O (Samples 2 and 3) using an FEI Titan Themis aberration-corrected S-TEM operated at 200kV. The point resolution in this aberration-corrected mode is 0.08nm. The microscope is fitted with a Super-XTM X-ray detector, which is a combination of 4 detectors for fast X-ray mapping in S-TEM mode. For the present experiment, one-nm resolution EDX maps were taken with an average beam current of 100pA. The size of the maps were 512 X 512 pixels and 50µs/pixel dwell time was used for collecting the signal. All maps are generated by summing over 10 frames. Drift correction during data collection and subsequent analysis were performed using Velox software. These results are shown in Fig. 4(a), with a standard deviation of pfu for Bi, and pfu for Fe. Magnetic samples
Two different sets of liquid-phase-epitaxy (LPE)-grown films with Bi Gd Lu Fe O , and Bi Gd Lu Fe Ga O per-formula-unit (pfu) were studied. These films were grown on (100) gadolinium gallium garnet (GGG) substrates. The bulk gyration, obtained from measurements in m − thick films, is different for the two different composition. For Bi Gd Lu Fe Ga O pfu films (Samples 2 and 3), the bulk 𝑔 = 0.029 + 0.003𝑖 , and for Bi Gd Lu Fe O (Sample 1), it is 𝑔 = 0.038 + 0.004𝑖 . Theoretical analysis
The transmission coefficients for right- and left-circularly polarized light are given by
21 cosh sinh i Et E nn ik n h n ik n hn , (M1) where denotes the optical polarization, i E , E are the amplitudes of the incident and transmitted to the substrate light circularly polarized components, respectively, n is the refractive index of the substrate, and where the refractive index of the ferrimagnetic film can be written as
22 2 n n g , k is the wave number in vacuum, and λ is the wavelength [25]. From Eq. (M1) and the fact that sinh sin ix i x and cosh cos ix x , one can find an expression for the transmitted light in terms of the component of the incident light
22 21 cos sin i EE nn h n g i n g h n gn g . (M2) Consequently, the Faraday rotation (FR), taking into account surface reflections and absorption losses, can be expressed as EEnn k gh i g k ghgnn k gh i g k ghg , (M3) where
22 2 n . References G. F. Dionne, Magnetic Oxides (Springer) 355-380 (2009). 2.
B. J. H. Stadler, T. Mizumoto, “Integrated Magneto-Optical Materials and Isolators: A Review,” IEEE Photonics Journal 6, 0600215 (2014). T. Mizumoto, R. Takei, and Y. Shoji, “Waveguide Optical Isolators for Integrated Optics,” IEEE J. of Quantum Electronics 48 (2), 252-260 (2012). 4.
M. Levy, “The On-Chip Integration of Magnetooptic Waveguide Isolators,” IEEE J. of Selected Topics in Quantum Electronics 8 (6), 1300-1306 (2002). 5.
S. Geller, M. A. Gilleo U.S. Patent 3,281,363, “Bismuth-Containing Garnets and Their Preparation,” (1963). 6.
A. Hamadeh, O. d’Allivy Kelly, C. Hahn, H. Meley, R. Bernard, A.H. Molpeceres, V.V. Naletov, M. Viret, A. Anane, V. Cros, S.O. Demokritov, J.L. Prieto, M. Muñoz, G. de Loubens, O. Klein, “Full Control of the Spin-Wave Damping in a Magnetic Insulator Using Spin-Orbit Torque”, Phys. Rev. Lett. 113(19), 197203 (2014). 7.
B. Heinrich, C. Burrowes, E. Montoya, B. Kardasz, E. Girt, Y.Y. Song, Y. Sun, M. Wu, “Spin Pumping at the Magnetic Insulator (YIG)/Normal Metal (Au) Interfaces”, Phys. Rev. Lett. 107, 066604 (2011). 8.
Z. Wang, M. Cherkasskii, B.A. Kalinikos, M. Wu, “Observation of spin-wave dark soliton pairs in yttrium iron garnet thin films”, Phys. Rev. B 91, 174418 (2015). 9.
I.V. Savochkin, M. Jäckl, V.I. Belotelov, I.A. Akimov, M.A. Kozhaev, D.A. Sylgacheva, A.I. Chernov, A.N. Shaposhnikov, A.R. Prokopov, V.N. Berzhansky, D.R. Yakovlev, A.K. Zvezdin, M. Bayer, “Generation of spin waves by a train of fs-laser pulses: a novel approach for tuning magnon wavelength”, Scientific Reports 7, 5668 (2017). 10.
M. Gomi, H. Furuyama, and M. Abe, "Strong magneto-optical enhancement in highly Ce-substituted iron garnet films prepared by sputtering," J. Appl. Phys. 70 (11), 7065-7067 (1991). 11.
R. Wolfe, R. A. Lieberman, V. J. Fratello, R. E. Scotti, and N. Kopylov, “Etch-tuned ridged waveguide magneto-optic isolator,” Appl. Phys. Lett. 56 (5), 426–428 (1990). 12.
M. Levy, R. M. Osgood, Jr., H. Hegde, F. J. Cadieu, R. Wolfe, and V. J. Fratello, “Integrated Optical Isolators with Sputter-Deposited Thin-Film Magnets,” IEEE Phot. Technol. Lett. 8 (7), 903-905 (1996). 13.
J. Fujita, M. Levy, R. M. Osgood, Jr., L. Wilkens, and H. Dötsch, “Waveguide Optical Isolator Based on Mach-Zehnder Interferometer,” Appl. Phys. Lett. 76 (16), 2158-2160 (2000). 14.
Y. Shoji, T. Mizumoto, H. Yokoi, I-Wei Hsieh, and R. M. Osgood, Jr., “Magneto-optical isolator with silicon waveguides fabricated by direct bonding,” Appl. Phys. Lett. 92, 071117 (2008). 15.
L. Magdenko, E. Popova, M. Vanwolleghem, C. Pang, F. Fortuna, T. Maroutian, P. Beauvillain, N. Keller, B. Dagens, “Wafer-scale fabrication of magneto-photonic structures in Bismuth Iron Garnet thin film”, Microelectronic Engineering 87, 2437 (2010). 16.
V.I. Belotelov, A.K. Zvezdin, “Magnetooptics and extraordinary transmission of the perforated metallic films magnetized in polar geometry”, JMMM 300, 1, e260-e263 (2006). A.N. Kalish, D.O. Ignatyeva, V.I. Belotelov, L. Kreilkamp, I.A. Akimov, A.V. Gopal, M. Bayer, and A.P. Sukhorukov, “Transformation of mode polarization in gyrotropic plasmonic waveguides.” Laser Physics, 24(9), p.094006 (2014). 18.
A.R. Prokopov, P.M. Vetoshko, A.G. Shumilov, A.N. Shaposhnikov, A.N. Kuz'michev, N.N. Koshlyakova, V.N. Berzhansky, A.K. Zvezdin, V.I. Belotelov, “Epitaxial BiGdSc iron-garnet films for magnetophotonic applications”, Journal of Alloys and Compounds 671 403-407 (2016). 19.
A.N. Kalish, R.S. Komarov, M.A. Kozhaev, V.G. Achanta, S.A. Dagesyan, A.N. Shaposhnikov, A.R. Prokopov, V.N. Berzhansky, A.K. Zvezdin, V.I. Belotelov. Magnetoplasmonic Quasicrystals: an Approach for Multiband Magneto-optical Response. Optica 5, issue 5, 617-623 (2018). 20.
M. Levy, A. Chakravarty, H. –C. Huang, R. M. Osgood, Jr., “Large magneto-optic enhancement in ultra-thin liquid-phase-epitaxy iron garnet films,” Appl. Phys. Lett. 107, 011104 (2015). 21.
O.V. Borovkova, H. Hashim, М.A. Kozhaev, S.A. Dagesyan, A. Chakravarty, M. Levy, and V.I. Belotelov, “TMOKE as efficient tool for the magneto-optic analysis of ultra-thin magnetic films,” Appl. Phys. Lett. 112, 063101 (2018). 22.
E. Popova, M. Deb, L. Bocher, A. Gloter, O. Stéphan, B. Warot-Fonrose, B. Berini, Y. Dumont, N. Keller, “Interplay between epitaxial strain and low dimensionality effects in a ferrimagnetic oxide”, J. Appl. Phys. 121, 115304 (2017). 23.
M. G. Moharam, E. B. Grann, D. A. Pommet, and T. K. Gaylord, J. Opt. Soc. Am. A 12, 1068 (1995). 24.
L. Li, J. Opt. A 5, 345 (2003). 25.
A. Zvezdin and V. Kotov, Modern Magnetooptics and Magnetooptical Materials,