Backscattering-immune one-way surface magnetoplasmons at terahertz frequencies

Surface magnetoplasmons (SMPs) in a basic physical model for the terahertz regime, which consists of a semi-infinite magnetized semiconductor with dielectric cladding terminated by a metal slab, are theoretically investigated. The dispersion properties of such SMPs are analyzed and examined in detail. It is shown that SMPs may follow three different kinds of dispersion diagrams, depending on the applied dc magnetic field intensity. Complete one-way propagation that operates within the band gap of the semiconductor is available for SMPs, and the one-way bandwidth reaches a maximum at a certain magnetic field intensity. Regular modes guided by the dielectric layer are also analyzed. These modes may cause the (complete) SMP one-way region to be compressed or even removed, but they can be suppressed by reducing the dielectric layer thickness. Owing to the mirror effect of the metal slab, one-way propagating and backscattering-immune basic SMPs can exhibit a larger propagation length than those sustained by a single dielectric–semiconductor interface. © 2015 Optical Society of America OCIS codes: (240.6680) Surface plasmons; (230.3810) Magneto-optic systems; (230.7380) Waveguides, channelled; (040.2235) Far infrared or terahertz. References and links 1. F. D. M. Haldane and S. Raghu, “Possible realization of directional optical waveguides in photonic crystals with broken time-reversal symmetry,” Phys. Rev. Lett. 100(1), 013904 (2008). 2. S. Raghu and F. D. M. Haldane, “Analogs of quantum-Hall-effect edge states in photonic crystals,” Phys. Rev. A 78(3), 033834 (2008). 3. The Quantum Hall Effect, edited by R. E. Prange and S. M. Girvin (Springer, 1987). 4. J. D. Joannopoulos, R. D. Meade, and J. N. Winn, Photonic Crystals: Molding the Flow of Light, 2nd ed. (Princeton University, 2008). 5. Z. Wang, Y. D. Chong, J. D. Joannopoulos, and M. Soljacic, “Reflection-free one-way edge modes in a gyromagnetic photonic crystal,” Phys. Rev. Lett. 100(1), 013905 (2008). 6. X. Ao, Z. Lin, and C. T. Chan, “One-way edge mode in a magneto-optical honeycomb photonic crystal,” Phys. Rev. B 80(3), 033105 (2009). 7. Z. Wang, Y. D. Chong, J. D. Joannopoulos, and M. Soljacic, “Observation of unidirectional backscatteringimmune topological electromagnetic states,” Nature 461(7265), 772-775 (2009). 8. Z. Yu, G. Veronis, Z. Wang, and S. Fan, “One-way electromagnetic waveguide formed at the interface between a plasmonic metal under a static magnetic field and a photonic crystal,” Phys. Rev. Lett. 100(2), 23902 (2008). 9. A. B. Khanikaev, A. V. Baryshev, M. Inoue, and Y. S. Kivshar, “One-way electromagnetic Tamm states in magnetophotonic structures,” Appl. Phys. Lett. 95(1), 011101 (2009). 10. L. Feng, M. Ayache, J. Huang, Y. L. Xu, M. H. Lu, Y. F. Chen, Y. Fainman, and A. Scherer, “Nonreciprocal light propagation in a silicon photonic circuit,” Science 333(6043), 729-733 (2011). #224928 $15.00 USD Received 15 Oct 2014; accepted 23 Dec 2014; published 15 Jan 2015 (C) 2015 OSA 26 Jan 2015 | Vol. 23, No. 2 | DOI:10.1364/OE.23.000950 | OPTICS EXPRESS 950 11. K. Fang, Z. Yu, and S. Fan, “Realizing effective magnetic field for photons by controlling the phase of dynamic modulation,” Nature Photon. 6(11), 782-787 (2012). 12. A. B. Khanikaev, S. H. Mousavi, W. K. Tse, M. Kargarian, A. H. MacDonald, and G. Shvets, “Photonic topological insulators,” Nature Mater. 12(3), 233-239 (2013). 13. M. C. Rechtsman, J. M. Zeuner, Y. Plotnik, Y. Lumer, D. Podolsky, F. Dreisow, S. Nolte, M. Segev, and A. Szameit, “Photonic Floquet topological insulators,” Nature 496(7444), 196-200 (2013). 14. H. B. Zhu and C. Jiang, “Broadband unidirectional electromagnetic mode at interface of anti-parallel magnetized media,” Opt. Express 18(7), 6914-6921 (2010). 15. B. Hu, Q. J. Wang, and Y. Zhang, “Broadly tunable one-way terahertz plasmonic waveguide based on nonreciprocal surface magneto plasmons,” Opt. Lett. 37(11), 1895-1897 (2012). 16. V. Kuzmiak, S. Eyderman, and M. Vanwolleghem, “Controlling surface plasmon polaritons by a static and/or time-dependent external magnetic field,” Phys. Rev. B 86(4), 045403 (2012). 17. X. G. Zhang, W. Li, and X. Y. Jiang, “Confined one-way mode at magnetic domain wall for broadband highefficiency one-way waveguide, splitter and bender,” Appl. Phys. Lett. 100(4), 041108 (2012). 18. J. J. Brion, R. F. Wallis, A. Hartstein, and E. Burstein, “Theory of surface magnetoplasmons in semiconductors,” Phys. Rev. Lett. 28(22), 1455-1458 (1972). 19. R. F. Wallis, J. J. Brion, E. Burstein, and A. Hartstein, “Theory of surface polaritons in anisotropic dielectric media with application to surface magnetoplasmons in semiconductors,” Phys. Rev. B 9(8), 3424-3437 (1974). 20. M. S. Kushwaha and P. Halevi, “Magnetoplasmons in thin films in the Voigt configuration,” Phys. Rev. B 36(11), 5960-5967 (1987). 21. M. S. Kushwaha and B. D. Rouhani, “Theory of magnetoplasmons in semiconductor superlattices in the Voigt geometry: A Green-function approach,” Phys. Rev. B 43(11), 9021-9032 (1991). 22. T. H. Isaac, W. L. Barnes, and E. Hendry, “Determining the terahertz optical properties of subwavelength films using semiconductor surface plasmons,” Appl. Phys. Lett. 93(24), 241115 (2008). 23. Y. S. Lee, Principles of Terahertz Science and Technology (Springer, 2009). 24. L. D. Landau, L. P. Pitaevskii, and E. M. Lifshitz, Electrodynamics of ContinuousMedia, Second Edition: Volume 8 (Course of Theoretical Physics) (Butterworth-Heinemann, 1984).


Introduction
One-way propagating electromagnetic (EM) modes, which were proposed by Raghu and Haldan [1,2] as analogues of quantum Hall edge states [3] in photonic crystals (PhCs) [4], are not only interesting for fundamental science but also important for practical applications. One-way EM modes are such that they are allowed to propagate in only one direction, and because of the absence of a back-propagating mode in the waveguide system, they are immune to backscattering at imperfections or bends [2,5], which is of particular interest for slow light systems. It was predicted that one-way EM modes can be confined at the edges of certain two-dimensional (2D) PhCs made of magnetic-optical (MO) materials, where time-reversal symmetry is broken by applying a dc magnetic field [1,2,5,6]. The existence of such modes was first experimentally demonstrated by Wang and his colleagues using MO PhCs in the microwave regime [7], and since then one-way EM modes have received increasing attention. Different schemes for realizing EM one-way propagation have also been proposed [8][9][10][11][12][13]; among them, the one based on surface plasmons seems to be most attractive owing to its robust mechanism and simple configuration [8,[14][15][16][17].
Surface plasmons (SPs) sustained by the interface between a dielectric and a metal will become nonreciprocal under an external dc magnetic field [8]. It has been shown that, when the electron cyclotron frequency is comparable to the plasma frequency of the metal, the SP asymptotic frequency will differ significantly in the forward and backward directions, causing SPs to propagate in only one direction within the interval between the two different asymptotic frequencies. To make one-way SPs immune to backscattering, Yu et al. introduced a PhC structure into the dielectric to remove the coupling between SP and bulk modes, where the latter is completely suppressed by the PhC's band gap [8]. Though such one-way SPs seem unrealistic because of the required magnetic field strength (B 0 ∼ 10 4 T), the scheme can be extended to the terahertz regime and they may become realizable with the use of semiconductors. A semiconductor generally has a plasma frequency in the terahertz regime, and meanwhile its electron effective mass is far less than the electron mass. As a result, the electron cyclotron frequency is often comparable to the plasma frequency for a semiconductor under normal magnetic field intensities. It is interesting to note that the study of SPs in magnetized semiconductors, referred to as surface magnetoplasmons (SMPs), began forty years ago [18][19][20][21]. We refer the reader to the work of Brion et al. for unidirectional propagation of SMPs in certain frequency regions. However, the SMPs studied in [18] are not robust against imperfections on the semiconductor surface, as the dielectric cladding is semi-infinite and bulk modes are allowed to propagate within it. Certainly, the bulk modes can be suppressed by introducing a PhC structure into the dielectric as proposed in [8]. For the terahertz regime, because metal resembles a perfect electric conductor, bulk mode suppression can be achieved by simpler means, i.e., by terminating the dielectric with a metal slab. One-way SMPs in such a basic physical model have recently been studied by Hu et al., and some interesting results were reported [15].
It is well known that SMP dispersion is characterized by its asymptotic frequencies, at which the propagation constant approaches infinity while the depth of field penetration into the dielectric approaches zero. Accordingly, in the basic model for the terahertz regime, the influence of the metal slab on the SMPs vanishes at the asymptotic frequencies. However, the SMP dispersion behaviour near one asymptotic frequency in Hu's work [15] is found to be quite different from that in Brion's [18]. Moreover, the dielectric layer sandwiched by a metal and a semiconductor (with a band gap) in the basic SMP model may support regular modes, which are guided based on total internal reflection. Evidently, regular modes are generally allowed to propagate in both the forward and backward directions, and their existence may make SMPs lose their robustness of one-way propagation. Unfortunately, attention has yet to be paid to the regular modes in the basic model. Furthermore, the loss of the semiconductor cannot be neglected generally for terahertz frequencies. Therefore, for basic one-way SMPs, it is necessary to examine whether the introduction of the metal slab will significantly increase the EM energy fraction in the semiconductor and cause serious propagation loss. In this paper, we will carefully investigate the characteristics of SMPs in the basic one-way propagation model and then physically clarify all the questions mentioned above.

Basic physical model of one-way SMPs
The basic physical model of SMPs consists of a semi-infinite semiconductor (x ≤ 0) with dielectric cladding (x > 0) terminated by a metal slab (see the upper panel of Fig. 1). The metal slab is used to suppress bulk modes in the dielectric with relative permittivity ε r , and the thickness of the dielectric layer is denoted by d. The semiconductor is uniformly magnetized by an external dc magnetic field B 0 in the −y direction, and SMPs propagate along the z direction. Gyroelectric anisotropy is induced in the semiconductor, with the permittivity tensor taking the form [18] where ω is the angular frequency, ω p is the plasma frequency of the semiconductor, ω c = eB 0 /m * (where e and m * are, respectively, the charge and effective mass of the electron) is the electron cyclotron frequency, and ε ∞ is the high-frequency (relative) permittivity of the semiconductor. The semiconductor is assumed to be lossless herein, and its loss will be later taken into account when we consider the SMP propagation loss. Moreover, we assume the metal slab to be a perfect electric conductor, which is a good approximation for the terahertz regime.
It should be noted that the metal slab in this system also has a mirror effect, making the system equivalent to a more complicated system, as shown in the lower panel of Fig. 1, where the latter requires the application of two oppositely directed magnetic fields. Owing to the mirror effect of the metal slab, this waveguide system is equivalent to that in the lower panel, where two opposite magnetic fields are applied.
SMPs are transverse-electric-polarized in the (2D) basic model, and the nonzero component of the magnetic field (H) can be written as in the dielectric layer for 0 < x ≤ d and in the semiconductor for x ≤ 0, where k is the propagation constant, α d = k 2 − ε r k 2 0 with k 0 = ω/c (where c is the light speed in vacuum), and α = k 2 − ε v k 2 0 with ε v = ε 1 − ε 2 2 /ε 1 being the Voigt permittivity. The nonzero components (E x and E z ) of the electric field (E) can be obtained straightforwardly from H y . The tangential component E z must vanish at the metal boundary x = d, yielding The boundary conditions further require the field components E x and H y to be continuous at the interface x = 0; based on this we can obtain which is the SMP dispersion relation. The second term in the dispersion equation is linear with respect to k, meaning that the SMP dispersion relation is asymmetric about k = 0. It should be pointed out that the SMP dispersion relation in the guiding system as shown in the lower panel of Fig. 1 is completely identical to Eq. (4).

SMP dispersion characteristics
SMP dispersion is characterized by its asymptotic frequencies, at which k → ±∞. Owing to the nonreciprocity of SMPs, we analyze their dispersion behaviour for the cases of k > 0 and k < 0 separately. We first consider the case of k > 0. As k → +∞, α ≈ k → +∞ and α d ≈ k → +∞, and from Eq. (4) two solutions can be found for ω: The second solution, Eq. (6), is surely the SMP asymptotic frequency for k > 0, as it reduces to ω + sp = ε ∞ /(ε ∞ + ε r )ω p for ω c = 0, which is just the result for SPs without magnetization. The first solution, Eq. (5), is found to correspond to a zero point of ε v . Thus, for ω = ω a , α = k since ε v = 0, and the dispersion equation (4) reduces to (1 + ε 2 /ε 1 )k = 0. This equation seems to hold for any (positive) k value, since the factor (1 + ε 2 /ε 1 ) vanishes in this case. To obtain a physical solution for ω = ω a , Eq. (4) should be solved in a proper way. Let us consider the limit case of ω → ω a , for which ε v → 0 and α ≈ k − ε v k 2 0 /(2k). In this situation, by deleting the factor (1 + ε 2 /ε 1 ), Eq. (4) can be rewritten as which is the actual equation that can be used to directly determine k for the special case of ω = ω a . Obviously, k is generally finite at ω = ω a ; thus, ω a is not an asymptotic frequency.
Because SMPs are a nonradiative mode, their dispersion curve lies outside the zone of bulk modes in the semiconductor, whose boundaries are described by ω = c|k|/ √ ε v . For a magnetized semiconductor, there exist two bulk-mode zones where ε v ≥ 0. The lower zone is located in the frequency region of ω a ≤ ω ≤ ω s , where ω s = ω 2 c + ω 2 p , and ε v increases from 0 at ω = ω a to +∞ at ω = ω s . The upper zone is located in the region of ω ≥ ω b , where ω b is the second zero point of ε v and is given by Correspondingly, there exist two band gaps for the semiconductor: one for ω < ω a and the other for ω s < ω < ω b . By using the bulk-mode zones as the position reference, the SMP dispersion diagram can be classified into three types for k > 0, depending on the applied magnetic field intensity, and we refer to them as dispersion diagrams I, II, and III, for which ω + sp is, respectively, located in the regions ω < ω a , ω a ≤ ω ≤ ω s , and ω > ω s . Dispersion diagram I with ω + sp < ω a pertains when ω c < ω In this case, the SMP dispersion relation has only a single band and lies completely within the lower semiconductor band gap. The dispersion band starts from the origin, rises to the right of the light line in the dielectric, which is described by ω = c|k|/ √ ε r , and finally approaches ω + sp as k → +∞. Dispersion diagram II with ω a ≤ ω + sp ≤ ω r pertains when ω Similar to dispersion diagram I, the SMP dispersion relation in this case has only a single dispersion band, but the dispersion band extends outside the lower semiconductor band gap. If the applied magnetic field is so large that ω c > ω (2) c , the SMPs follow dispersion diagram III with ω + sp > ω s . In this case, the dispersion relation is split into two branches by the lower bulk-mode zone of the semiconductor. The lower branch is just terminated by the boundary of the bulk-mode zone. The upper branch starts from ω s at k = k + c f and approaches ω (+) sp as k → +∞. As ω → ω s , ε 1 → 0 and ε v → ∞, and from Eq. (4) we derive where If the upper branch starts from the light line of the dielectric, it can be derived from Eq. (4) that α + k + c f ε 2 /ε 1 = 0. This equation cannot hold for positive k + c f , since both ε 1 and ε 2 are positive for ω > ω s , so the upper branch cannot start from the light line. In the case when ω c is close to ω Clearly, the value of k + c f approaches +∞ as ω c → ω (2) c (as expected), and it decreases when ω c increases from ω (2) c . We now consider the case of k < 0. Similar to the case of k > 0, as k → −∞, two solutions are found from Eq. (4). One solution is ω = ω b and the other is Evidently, ω − sp is really the SMP asymptotic frequency for k < 0. The first solution ω b is just the second zero of ε v , at which Eq. (4) holds for any negative k value, similar to the case of k > 0 at ω = ω a . For ω = ω b , we can also similarly derive which is the actual equation that can be used to directly determine k for this special frequency. For the case of k < 0, ω − sp is smaller than ω a for any ω c , so there always exists a complete SP dispersion band, which starts from the origin, rises, and approaches ω − sp as k → −∞. However, for k < 0, an additional dispersion branch can always be found in the region where ω > ω s . This branch starts from the light line of the dielectric and ends at the boundary of the upper bulk-mode zone of the semiconductor. At the starting point of this additional branch, α d = 0; thus, from Eq. (4) we have α + kε 2 /ε 1 = 0 with k = √ ε r k 0 , from which we find that ε 1 is equal to ε r and the cutoff frequency is for the additional branch. Evidently, the existence of the additional branch only requires ε ∞ > ε r , and this condition is generally met by transparent dielectrics in the terahertz regime. To validate our analysis above, we numerically calculate the SMP dispersion relation with Eq. (4) for various ω c values. As an example, the basic parameters of the system are taken as follows: ε ∞ = 15.6 (InSb) [22], ε r = 2.28 (polymer) [23], and d = 0.16λ p (where λ p is the vacuum wavelength for the plasma frequency ω p ). From now on, the semiconductor is assumed to be InSb in all numerical calculations. For the present system, we find that ω (1) c = 0.066ω p and ω (2) c = 0.148ω p . Figures 2(a)-2(c) display the SMP dispersion relations for the cases of ω c = 0.05ω p , 0.1ω p , and 0.2ω p , respectively. It is obvious that, as expected, they correspond to dispersion diagrams I, II and III, respectively. As can be seen from Figs. 2(a)-2(c), SMPs always follow two dispersion branches for k < 0, and the upper one (i.e., the additional branch) occurs in a narrow range compared with the lower one. Figure 2(d) shows a magnified view of the upper branch (with k < 0) in Fig. 2(a); it clearly indicates that this branch starts from the light line and ends at the boundary of the bulk-mode zone. For SMPs with k > 0, the second dispersion branch is observed only in Fig. 2(c). Our numerical analysis also shows that the second dispersion branch can be found only when ω c > 0.148ω p . All the numerical results agree well with our analytical analysis.

Complete SMP one-way propagation
It is evident that, for each kind of dispersion diagram, there always exist certain frequency ranges in which SMPs are allowed to propagate unidirectionally. For dispersion diagrams I and II, as shown in Figs. 2(a) and 2(b), SMPs only propagate forward in the frequency range from ω − sp to ω + sp , while its upper branch with ω ≥ ω − c f allow only backward propagation. In dispersion diagram III, there may exist three unidirectional-propagation ranges, as shown in Fig. 2(c); two correspond to the forward direction and the third to the backward direction. However, only the portion of the unidirectional-propagation range located within the semiconductor band gap is of interest, and we refer to it as the complete one-way region (or range). In the complete one-way region, SMPs are immune to backscattering resulting from imperfections in the system, as their coupling with bulk modes in the semiconductor is completely eliminated. Because a larger complete one-way range is expected for a normal (external) magnetic field, in what follows, we restrict ourselves to the complete one-way region within the lower semiconductor band gap (ω < ω a ). This complete one-way region ranges from ω − sp to ω + sp for dispersion diagram I (ω c ≤ ω (1) c ), and it ranges from ω − sp to ω a for dispersion diagrams II and III (ω c > ω (1) c ). Correspondingly, the complete one-way bandwidth (Δω ow ) has the form for ω c > ω  The dependence of Δω ow on ω c is plotted in Fig. 3(a). Three different dielectrics are analyzed: air, polymer (ε r = 2.28), and silica (ε r = 11.68), all of which are transparent media commonly used in the terahertz regime [23]. As seen from Fig. 3(a), Δω ow is independent of ε r when ω c < 0.031ω p , which just corresponds to the value of ω (1) c for air. Therefore, under a small magnetic field, it is impossible to increase Δω ow by choosing the dielectric. However, the maximal value of Δω ow , which is equal to the value of ω (1) c and occurs at ω c = ω (1) c , closely depends on ε r . The larger ε r is, the larger the maximal Δω ow is. It is found that the maximal values of Δω ow are 0.031ω p , 0.066ω p , and 0.24ω p for air, polymer, and silica, respectively. These values also represent the peak positions ω (1) c for the corresponding dielectrics. Figure  3(b) shows the complete one-way range for the case of air. The upper limit of the one-way range increases with ω c when ω c < 0.031ω p , and it decreases with ω c when ω c > 0.031ω p , while the lower limit always decreases with ω c . Finally, we should note that the semiconductor is inherently lossy as it is dispersive; consequently, SMPs propagating in it are lossy. In general, the SMP loss will become very serious when k √ ε r k 0 . To reduce the one-way SMP loss to an acceptable level, ω c should be chosen to be a little larger than ω (1) c . By doing so, the SMP dispersion curve in the complete one-way region does not deviate too far from the light line, while the one-way bandwidth is still large.

Regular mode guided by the dielectric layer
In the basic one-way SMP model, the dielectric layer can still support (transverse-electricpolarised) regular modes, which are guided by a means of zigzag reflections at the surfaces of the metal and the semiconductor (with a band gap). As regular modes are generally allowed to propagate in both forward and backward directions, their existence causes SMPs to lose their robustness of one-way propagation against imperfections. The dispersion equation for regular modes can be obtained from Eq. (4) by replacing α d with −ip, where p = ε r k 2 0 − k 2 is realvalued, and it has the form Obviously, the dispersion relation for regular modes is also asymmetric about k = 0. Because regular modes have a tangential electric field component (E z ) that vanishes at the metal surface, they should exhibit transverse resonances and related cutoff frequencies. These cutoffs increase as d decreases; hence, it is possible to remove the regular modes from the complete one-way region by reducing the dielectric layer thickness.
Owing to the nonreciprocity, the cutoff point of a regular mode may deviate from k = 0. However, this deviation should be very small under a normal magnetic field, as ω c is generally much smaller than ω p (e.g., for InSb at room temperature, ω c is equal to 0.1ω p for B 0 = 0.1 T). If we let k ≈ 0, then p ≈ √ ε r k 0 , and from Eq.
In the region below ω a , ε v is always negative, so Eq. (17) has no solution if √ ε r k 0 d < π/2. This means that even the fundamental regular mode with the lowest cutoff will vanish at this frequency. By setting k 0 = ω a /c, the critical thickness (d c ) of the dielectric layer is found to be below which no regular mode appears in the region for ω < ω a ; i.e., the complete one-way region is preserved for SMPs. Obviously, d c is inversely proportional to √ ε r . For ω c = 0.1ω p , we find that d c = 0.26λ p for ε r = 1, d c = 0.17λ p for ε r = 2.28, and d c = 0.077λ p for ε r = 11.68. To illustrate the regular mode guided by the dielectric layer, we numerically calculate its dispersion relations using Eq. (16) for various d values, and the results are plotted in Fig. 4, where ε r = 2.28 and ω c = 0.1ω p (corresponding to d c = 0.17λ p ) for all cases. Figures 4(a)-4(d) show the dispersion relations (lines with filled circles) of the (fundamental) regular mode for d = d c , 0.2λ p , 0.25λ p , and 0.29λ p , respectively. As expected, in the case of d = d c , the fundamental regular mode has a cutoff of ω a , and the complete one-way SMP region is preserved. When d increases from d c , the cutoff of the regular mode shifts downward. As a result, the complete one-way region is compressed, and its upper limit is just determined by the cutoff frequency. Figure 4(d) indicates that, when d = 0.29λ p , the cutoff drops to ω − sp and consequently the complete one-way region vanishes completely. Using the finite element method (FEM), we performed a simulation to demonstrate the regular mode in the basic SMP model. For this purpose, the parameters of the system are taken to be the same as in Fig. 4(c). In the simulation, a magnetic current line source is used to excite SMPs and two operation frequencies of f A = 0.905 f p and f B = 0.94 f p (where f p = ω p /2π) are considered, as marked with points A and B in Fig. 4(c). The plasma frequency is set to be f p = 2 THz, so f A = 1.81 THz and f B = 1.88 THz. To examine the robustness of one-way propagating SMPs, we consider both cases with and without a metal rod obstacle. The metal rod with a radius of 8 μm is inserted at the polymer-semiconductor interface at a horizontal distance of 150 μm from the source. All simulated results are displayed in Fig. 5. The left panels show the SMP transmission at f A = 1.81 THz, where the source is placed at the polymer-semiconductor interface. It can be clearly observed from the upper left panel that SMPs propagate only in the forward direction. Furthermore, SMPs at this frequency can go around the obstacle without any power loss, as no backward-propagating wave emanates from the left side of the source (see the lower left panel). For f B = 1.88 THz, to avoid strong direct excitation of the regular mode, the source is placed inside the semiconductor at 30 μm below the surface. Nevertheless, a weak backward-propagating wave, which corresponds to the regular mode and is directly excited by the source, is still observed in the upper right panel of Fig. 5. As shown in the lower right panel, the backward-propagating wave is enhanced in the presence of the obstacle. Evidently, in this case, coupling between SMPs and the backward-propagating regular mode occurs via scattering at the obstacle. The results shown in the right panels indicate that SMPs at f B = 1.88 THz are no longer immune to backscattering owing to the existence of the regular mode in the dielectric layer. The simulated results agree well with our previous analysis.

Mirror effect of the metal slab
As a dispersive medium, the semiconductor is inherently lossy, and, consequently, SMPs suffer a propagation loss in it. For a lossy semiconductor, the expressions for ε 1 and ε 2 in Eq. (1) become where ν is the electron scattering frequency. In the lossy case, the propagation constant k in Eq. (4) becomes complex, i.e., k = k r + ik i (where both k r and k i are real-valued), and the SMP propagation length is determined by L = 1/(2k i ). Evidently, in the case of ν ω p , the SMP loss is proportional to the fraction of EM energy in the semiconductor. Correspondingly, the larger the EM energy fraction of SMPs is distributed within the dielectric layer, the larger the propagation length is. In the basic one-way SMP model, the metal slab has two opposing effects on the SMP energy distribution: 1. As the metal slab compresses the dielectric region, the energy fraction in the dielectric layer is decreased, and 2. the metal slab introduces the mirror effect, which will enhance the EM energy density and thus increase that energy fraction. Owing to the mirror effect, a SMP in the present model may be viewed as a symmetric combination of a SMP associated with a single (dielectric-semiconductor) interface and its image produced by the metal slab; the field in the dielectric layer is therefore enhanced compared with the case of a single interface. However, the mirror effect is effective only when d ≤ 1/α d , which is the field penetration depth into the dielectric. Thus, for robust one-way SMPs, an interesting question arises as to whether the propagation loss is aggravated or lessened (compared with SMPs at a single interface) by the metal slab or whether the mirror effect of the metal slab is the dominant effect on the SMP energy distribution.
To clarify this question, we numerically calculate the propagation length L using Eq. (4) for a lossy semiconductor. As a numerical example, we take ν = 0.0025ω p for the semiconductor, and other parameters of the system are as follows: ε ∞ = 15.6, ε r = 2.28, and ω c = 0.1ω p . Figure 6(a) shows the calculated L value as a function of d for the frequency ω = 0.918ω p . This frequency lies at the centre of the complete one-way region for d ≤ d c , i.e., ω = (ω − sp + ω a )/2. With the above parameters, the critical thickness is d c = 0.17λ p . For comparison, we also calculate the propagation length (L 0 ) for SMPs at a single interface, and we find that L 0 = 7.7λ (λ = 2πc/ω) for ω = 0.918ω p , which is indicated in Fig. 6(a) by the horizontal dashed line. From Fig. 6(a), it is seen that L < L 0 for small d, meaning that the metal slab aggravates the SMP loss owing to the compression of the dielectric region. L rapidly increases with d, becoming larger than L 0 at d = d c , meaning that the mirror effect becomes the dominant effect on the energy distribution. Our calculation gives L = 8.3λ for d = d c (L is inversely proportional to ν, and for this case L = 20.75λ if ν = 0.001ω p ). L reaches a maximal value of 8.4λ at d = 0.20λ p , but the cutoff of the fundamental regular mode in this case is found to be below ω. The inset of Fig. 6(a) shows the distribution of L in the complete one-way region for d = d c , and the corresponding results for L 0 are also plotted with a dashed line for comparison. Clearly, owing to the mirror effect of the metal slab, SMPs in the present system can possess a larger propagation length than SMPs at a single interface for the entire one-way region.  To gain a deeper insight into the mirror effect, we further calculate the SMP EM energy in the dielectric layer and semiconductor in the lossless case (ν = 0). The other parameters are kept the same as in Fig. 6(a). As the semiconductor is strongly dispersive, the EM energy density in it can be expressed as [24] The SMP EM energy fraction in the dielectric layer is where U d = (1/4)(ε 0 ε r |E| 2 + μ 0 |H| 2 ), which is the energy density in the dielectric layer. The calculated results of η d for various d values are plotted in Fig. 6(b). As expected, the dependence of η d on d is similar to that of L, as shown in Fig. 6(a). η d first increases with d and then reaches a maximum at d ≈ 0.2λ p . As d continues to grow, it decreases and finally approaches a limit, which corresponds to that for SMPs at a single dielectric-semiconductor interface. The behaviour of η d clearly indicates that the mirror effect of the metal slab enhances the EM energy density in the dielectric layer and thus may increase the energy fraction.

Conclusions
In summary, we have theoretically studied terahertz SMPs in a basic physical model comprising a semi-infinite magnetized semiconductor with dielectric cladding terminated by a metal slab. It has been shown that, depending on the applied magnetic field intensity, such SMPs may follow three different kinds of dispersion diagrams. For each kind of dispersion diagram, there always exists a complete one-way SMP region, which lies within the lower band gap of the semiconductor. If the magnetic field is strong enough, another complete one-way band appears in the upper band gap of the semiconductor, but it generally has a bandwidth significantly smaller than that in the lower band gap. The bandwidth of the (complete) one-way propagation closely depends on the magnetic field intensity, and it reaches a maximum at a finite magnetic field. Regular modes guided by the dielectric layer have also been investigated. The regular modes are allowed to propagate in both the forward and backward directions, and when they are present, the complete one-way region is often compressed or even removed. However, regular modes can be suppressed by reducing the dielectric layer thickness. Moreover, the mirror effect of the metal slab in this basic one-way SMP model has been analyzed. It has been shown that, because the mirror effect can significantly increase the EM energy fraction in the dielectric layer, the SMPs immune to backscattering in the basic physical model can exhibit a larger propagation length than those at a single dielectric-semiconductor interface.