Switchable plasmonic routers controlled by external magnetic fields by using magneto-plasmonic waveguides

We analytically and numerically investigate magneto-plasmons in metal films surrounded by a ferromagnetic dielectric. In such waveguide using a metal film with a thickness exceeding the Skin depth, an external magnetic field in the transverse direction can induce a significant spatial asymmetry of mode distribution. Superposition of the odd and the even asymmetric modes over a distance leads to a concentration of the energy on one interface which is switched to the other interface by the magnetic field reversal. The requested magnitude of magnetization is exponentially reduced with the increase of the metal film thickness. Based on this phenomenon, we propose a waveguide-integrated magnetically controlled switchable plasmonic routers with 99-%-high contrast within the optical bandwidth of tens of THz. This configuration can also operate as a magneto-plasmonic modulator.

however this can be increased to 12% 13 by adding a higher refractive index dielectric layer on top of the interferometer and by optimizing the ferromagnet layer.
In this paper, we present a new conception for a switchable plasmonic router based on the magnetically-induced spatial asymmetry of planar SPP modes in an external magnetic field. Asymmetric guiding modes are generally supported by waveguides with asymmetric cross sections. The asymmetry can also be achieved through active approaches, for example, by using power-dependent nonlinear dielectrics 26,27 or by introducing an external magnetic field [28][29][30] . We analytically and numerically study the magnetically-induced mode asymmetry in a metal film surrounded by a ferromagnetic dielectric. An experimentally achievable magnetization can induce a significant mode asymmetry for a metal film with thickness exceeding the Skin depth. For a certain degree of mode asymmetry, a superposition of the even and the odd modes over a distance leads to a high-contrast-asymmetric splitting within a broad optical bandwidth. This configuration can operate as a waveguide-integrated switchable magneto-plasmonic router with high contrast and broad optical bandwidth controlled by an external magnetic field.

Results and Discussion
Magnetically-induced mode asymmetry. We study analytically and numerically by solving the Maxwell equations the magnetically-induced mode asymmetry in a metal film surrounded by a ferromagnetic dielectric under an external magnetic field inducing a magnetization My M = , where ŷ is the unit vector in y-direction, as shown in Fig. 1. The higher absorption in ferromagnetic metals results in shorter propagation lengths of SPPs and makes these materials unfavorable for many applications. Ferromagnetic dielectric materials offer the advantage of much lower loss which makes a noble metal film surrounded by a ferromagnetic dielectric attractive for diverse applications. Excitation of SPP modes in diverse metal film structures can be achieved by the Kretschman configuration or a grating (see e.g. ref. 31 ). For a sufficiently large thickness of the surrounding dielectric layer, the SPP modes are not significantly influenced by the substrate material. We assume that the metal layer is surrounded by semi-infinite dielectric layers.
The permittivity tensor ε of the ferromagnetic dielectric can be expressed as follows.
where a is the magneto-optical susceptibility. For the fundamental TM mode, we can get the following dispersion relation from the Maxwell equation and boundary conditions.
, t is the thickness of the metal film, β is the wavenumber, and ε m and ε d are the permittivity of the metal and the ferromagnetic dielectric, respectively. In the absence of the external magnetic field g = 0, the dispersion relation Eq.
The superscripts o and e represent the odd and the even modes, respectively. We note that Δβ g is in accordance with Eq. (11) of ref. 33 and Eq. (4) of ref. 18 which describes the magnetically-induced wavenumber change by a single interface between a metal and a ferromagnetic dielectric. If the thickness t is infinite, the dispersion relation Eq. 18 .
We define the degree of the mode asymmetry as follows.
In the absence of an external magnetic field, the intensity distributions for both modes are symmetric, m o = 1 and m e = −1 but an external magnetic field induces an asymmetry of the mode distribution. By using Eq. (4), in the first order of |(β − β 0 )/β 0 |, Eq. (7) is approximated as follows.
here g 0 characterize a gyration requested to induce a significant asymmetry of Eq. (9) shows that the characteristic gyration g 0 is exponentially reduced with the increase of the metal film thickness t. For a metal film thickness exceeding the Skin depth, a significant asymmetry can be achieved by introducing a small value of g or a corresponding external magnetic field. Figure 2 shows the asymmetry m for the odd (red curves) and the even (blue curves) modes in the magneto-plasmonic waveguide, as shown in Fig. 1, with t = 100 nm at λ = 800 nm in dependence on g. The diagonal permittivity of the ferromagnetic dielectric is assumed to be ε d = 6.25. The experimental data for the permittivity of silver 34  High-contrast magneto-plasmonic modulation. Based on the above results of the magneticallyinduced mode asymmetry, a configuration for a high-contrast magnetically controlled switchable plasmonic router is suggested. A coupled plasmonic waveguide without external magnetic field is added to the left side of the magneto-plasmonic waveguide structure, as shown in Fig. 3. The long-range SPP mode (odd mode), which is incident to the left port, is presented by a superposition of the odd and the even modes in the magneto-plasmonic waveguide. The coefficients c 1 and c 2 can be obtained from the continuity of field distribution at the interface between the non-magnetic and the magneto-plasmonic waveguides, shown by the vertical white-dashed-line in Fig. 3. While the magneto-plasmon propagates along the waveguide, the field distribution during propagation varies due to the difference in wavenumber between the odd and the even modes, β e − β o . By using Eqs (7) and (10), the intensity ratio   From Eq. (14), one can see that the ratio R becomes extremely high at the characteristic gyration g = ±g 0 . The magnitude of magnetization for a high contrast exponentially decreases with the increase of the metal film thickness t.
Let's investigate the half period L 0 for the case of a high contrast. At the high-contrast points g = ±g 0 , we have  in the magneto-plasmonic waveguide. The left and the right parts to the white dashed line are in the absence and the presence of the external magnetic fields inducing g = ±g 0 , respectively. One can see that the field intensity is concentrated on the bottom interface, which is switched to the top interface by the reversal of the magnetic field direction. The distance between the point 1 and the point 2 is coincident with L 0 = 16 μm calculated by Eq. (16). The field distribution is recovered at the point 3 except a power attenuation due to the ohmic loss in the metal.
A configuration for the magnetically controlled switchable plasmon router is shown in Fig. 4. The value of g 0 can be reduced by increasing the metal thickness to a value of 110 nm. g 0 = 0.026 and L 0 = 26 μm have been calculated by Eqs (9) and (16). We note that the experimentally observed gyration values of 0.03 or 0.06 for Bi-substituted iron garnet (BIG) have been reported in ref. 35 , in the Table 1 of ref. 36 and in the Methods of ref. 37 . We assume a bending radius of 300 nm for which the numerical simulation shows a small value of the bending loss of 1.2%. For a more realistic simulation we take into account the absorption of 300 dB/cm 38 in the ferromagnetic dielectric. The calculation using Eq. (3) predicts Imβ ≈ 5 × 10 4 /m including a contribution of the absorption in the ferromagnetic dielectric, 1×10 4 /m, which is not ignorable compared to β e − β o ≈ 12 × 10 4 /m. The absorption in the metal and the ferromagnetic dielectric results in an insertion loss, however does not influence the high contrast because the intensity reduction due the absorption occur equally in the both channels.
In Fig. 5 numerical results for the configuration show a high-contrast magneto-plasmonic modulation and channel switching by the magnetic field reversal.
We estimate the modulation contrast A 1 = P ch1 /(P ch1 + P ch2 ) and A 2 = P ch2 /(P ch1 + P ch2 ) defined by the ratio of powers through the two channels in the configuration in dependence on g (Fig. 6a) and the wavelength λ (Fig. 6b) for L 0 = 26 μm. The numerical results (red crosses) well agree with the analytically predicted results by Eq. (12) (blue curves). The power splitting ratio is tunable by controlling the external magnetic field (Fig. 6a). In the vicinity of g = ±g 0 = ±0.026 the power splitting ratio reaches a value larger than 20 dB, which allows a high-contrast magneto-plasmonic modulation. The 99-%-high contrast is maintained in a broad wavelength range from 770 nm to 820 nm (Fig. 6b).
The broad bandwidth shown in Fig. 6b allows a high-contrast modulation of femtosecond pulses. The temporal power of the SPP pulses in both channels are shown in Fig. 7. The full-width at half maximum (FWHM) of the pulse is 20 fs and the central wavelength is 800 nm. By the magnetic field reversal, the ultrashort pulse is amplitude-shift keying (ASK) modulated with 99-%-high contrast.
In conclusion, we derived an analytical formula for the magnetically-induced mode asymmetry in metal films, with a thickness exceeding the Skin depth, surrounded by ferromagnetic dielectrics. In such waveguide, a significant spatial asymmetry of the mode distribution can be induced by an external magnetic field in the transverse direction which is exponentially reduced with the increase of the metal film thickness. Superposition of the odd and the even asymmetric modes over a distance presents complete energy concentration on one interface which is switched to the other interface by the magnetic field reversal. Based on this phenomenon, we proposed a Here g 0 = 0.026 and L 0 = 26 μm are calculated by Eqs (9) and (16). The metal thickness has been increased to a value of 110 nm. The other parameters are the same as in Fig. 3. Figure 6. Modulation contrast A 1 = P ch1 /(P ch1 + P ch2 ) and A 2 = P ch2 /(P ch1 + P ch2 ) defined by the ratio of powers through the two channels in the configuration in dependence on the external magnetic field g (a) and the wavelength λ (b). Other parameters are the same as in Fig. 5. The red crosses and the blue circles are calculated numerically and the red and the blue curves from the analytical formula. waveguide-integrated magnetically controlled switchable plasmonic router with 99-%-high contrast within an optical bandwidth of tens of THz. This configuration can also operate as a magneto-plasmonic modulator.

Methods
Analytical formulas (4) and (8)  The temporal response to the input pulse was calculated by the superposition of the responses of the Fourier components of the time-depending pulse.
The experimental data for the permittivity spectra of silver 34 have been used as the permittivity ε m of the metal film in the structure of Fig. 1.  Figure 7. Temporal power characteristics of the SPPs in the configuration (Fig. 4a). The full-width at half maximum (FWHM) of the pulse is 20 fs and the central wavelength is 800 nm. The red crosses represent the power through channel 1 while the blue circles the power through channel 2. Other parameters are same as in Fig. 5.