Azimuthally polarized surface plasmons as effective terahertz waveguides

Quite recently, it was found that metal wires can effectively guide terahertz radiation. Based on the fact that the absolute values of the relative permittivities of metals in the spectral region of terahertz radiation are huge, we here analyse the properties of this kind of waveguide and explain the related experimental results. In particular, we show that the observed waveguiding is due to the propagation of an azimuthally polarized surface plasmon along the wire. Some related aspects, such as the choice of metal and the slowly decaying modal field, are also discussed. In particular, we show that, if a copper wire with a radius of 0.45 mm is used, the attenuation coefficient is smaller than 2×10 cm in the whole range of 0.1~1 THz. © 2005 Optical Society of America OCIS codes: (260.3090) Infrared, far; (240.6680) Surface plasmon, (230.7370) Waveguides; (260.3910) Metal, optics of. References and links 1. D. M. Mittleman, ed. Sensing with Terahertz Radiation (Springer, Heidelberg, 2002). 2. P. R. Smith, D. H. Auston, and M. C. Nuss, “Subpicosecond photoconducting dipole antennas,” IEEE J. Quant. Electron. 24, 255-260 (1998). 3. M. Exter and D. Grischkowsky, “Characterization of an optoelectronic terahertz beam system,” IEEE Trans. Microwave Theory Tech. 38, 1684-1691 (1990). 4. P. U. Jepsen, R. H. Jacobsen, and S. R. Keiding, “Generation and detection of terahertz pulses from biased semiconductor antennas,” J. Opt. Soc. Am. B 13, 2424-2436 (1996). 5. D. M. Mittleman, R. H. Jacobsen, and M. C. Nuss, “T-ray imaging,” IEEE J. Select.Top. Quant. Electron. 2, 679-692 (1996). 6. R. H. Jacobsen, D. M. Mittleman, and M. C. Nuss, “Chemical recognition of gases and gas mixtures with terahertz waves,” Opt. Lett. 21, 2011-2013 (1996). 7. R. M. Woodward, V. P. Wallace, D. D. Arnone, E. H. Linfield, and M. Pepper, “Terahertz pulsed imaging of skin cancer in the time and frequency domain,” J. Biol. Phys. 29, 257-261 (2003). 8. K. Kawase, Y. Ogawa, and Y. Watanabe, “Non-destructive terahertz imaging of illicit drugs using spectral fingerprints,” Opt. Express 11, 2549-2554 (2003). 9. S. Wang and X. –C. Zhang, “Pulsed terahertz tomography,” J. Phys. D 37, R1-R36 (2004). 10. R. W. McGowan, G. Gallot, and D. Grischkowsky, “Propagation of ultrawideband short pulses of THz radiation through submillimeter-diameter circular waveguides,” Opt. Lett. 24, 1431-1433 (1999). 11. G. Gallot, S. P. Jamison, R. W. McGowan, and D. Grischkowsky, “Terahertz waveguides,” J. Opt. Soc. Am. B 17, 851-863 (2000). 12. R. Mendis and D. Grischkowsky, “Plastic ribbon THz waveguides,” J. Appl. Phys. 88, 4449-4451 (2000). 13. S. P. Jamison, R. W. McGown, and D. Grischkowsky, “Single-mode waveguide propagation and reshaping of sub-ps terahertz pulses in sapphire fiber,” Appl. Phys. Lett. 76, 1987-1989 (2000). 14. R. Mendis and D. Grischkowsky, “Undistorted guided-wave propagation of subpicosecond terahertz pulses,” Opt. Lett. 26, 846-848 (2001). 15. R. Mendis and D. Grischkowsky, “THz interconnect with low loss and low group velocity dispersion,” IEEE Microwave Wireless Comp. Lett. 11, 444-446 (2001). 16. S. Coleman and D. Grischkowsky, “A THz transverse electromagnetic mode two-dimensional interconnect layer incorporating quasi-optics,” Appl. Phys. Lett. 83, 3656-3658 (2003). 17. K. Wang and D. M. Mittleman, “Metal wires for terahertz wave guiding,” Nature (London), 432, 376-379 (2004). 18. H. Raether, Surface Plasmons on Smooth and Rough Surfaces and on Gratings, (Springer, Berlin, 1988). 19. http://www.surfaceplasmonoptics.org (C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 511 #6079 $15.00 US Received 15 December 2004; revised 12 January 2005; accepted 12 January 2005 20. Q. Cao and Ph. Lalanne, “Negative role of surface plasmons in the transmission of metallic gratings with very narrow slits,” Phys. Rev. Lett. 88, 057403 (2002). 21. S. I. Bozhevolnyi, “Waveguiding in surface plasmon polariton band gap structures,” Phys. Rev. Lett. 86, 3008-3011 (2001). 22. M. I. Stockman, “Nanofocusing of optical energy in tapered plasmonic waveguides,” Phys. Rev. Lett. 93, 137404 (2004). 23. C. A. Pfeiffer, E. N. Economou, and K. L. Ngai, “Surface polaritons in a circularly cylindrical interface: surface plasmons,” Phys. Rev. B 10, 3038-3051 (1974). 24. U. Schröter and A. Dereux, “Surface plasmon polaritons on metal cylinders with dielectric core,” Phys. Rev. B 64, 125420 (2001). 25. M. Born and E. Wolf, Principles of Optics, 5th ed. (Pergamon Press, Oxford, 1975). 26. G. N. Watson, A Treatise on the Theory of Bessel Functions, 2nd ed. (Cambridge U. Press, Cambridge, UK 1966). 27. M. A. Ordal, R. J. Bell, R. W. Alexander, L. L. Long, and M. R. Querry, “Optical properties of fourteen metals in the infrared and far infrared: Al, Co, Cu, Au, Fe, Pb, Mo, Ni, Pd, Pt, Ag, Ti, V, and W,” Appl. Opt. 24, 4493-4499 (1985). 28. http://hyperphysics.phy-astr.gsu.edu/hbase/tables/magprop.html#c1 29. http://www.stainless-rebar.org/grades.htm 30. C. Weber and J. Fajans, “Saturation in “nonmagnetic” stainless steel,” Rev. Scientific Instruments 69, 36953696 (1998). 31. S. Quabis, R. Dorn, M. Eberler, O. Glöckl, and G. Leuchs, “Focusing light to a tighter spot,” Opt. Commun. 179, 1-7 (2000). 32. R. Dorn, S. Quabis, and G. Leuchs, “Sharper focus for a radially polarized light beam,” Phys. Rev. Lett. 91, 233901 (2003). 33. Q. Zhan and J. R. Leger, “Focus shaping using cylindrical vector beams,” Opt. Express 10, 324-331 (2002).


Introduction
The terahertz (THz) region of the electromagnetic spectrum is located between microwave and optical frequencies, normally defined as the range from 0.1 to 10 THz (or correspondingly, from 30 µm to 3 mm in wavelength).In recent years, THz radiation has attracted a lot of interests [1][2][3][4][5][6][7][8][9], because it offers significant scientific and technological potential in many fields, such as in sensing, in imaging, and in spectroscopy.However, waveguiding in this intermediate spectral region is still a challenge [10][11][12][13][14][15][16].Both of the conventional metal waveguides for microwave radiation and the dielectric fibers for visible radiation cannot be used to effectively guide THz radiation.The obstacles come from the high loss from the finite conductivity of metals or the high absorption coefficient of dielectric materials in this spectral range.Quite recently, Wang and Mittleman [17] found that a simple metal wire can be used as effective THz waveguide.This finding paves the way for a wide range of new applications for THz sensing and imaging.In Ref. [17], the authors focused on the report of experimental observations, however, without a deep theoretical explanation.
In this paper, we present the needed theory and explain the related experimental results.In particular, we show that the observed waveguiding effect is due to the propagation of an azimuthally polarized surface plasmon (APSP) along the wire.In addition, we also discuss some related problems.In particular, we suggest the use of copper wires instead of stainless steel wires as more effective THz waveguides.

Explanations for THz metal wire waveguides
It is well known that there exists an electromagnetic bound state at a flat metal-dielectric interface.This bound state can only exist for the TM polarization and is called surface plasmon (SP) [18,19].SPs can be excited by periodic structures like metal gratings [20], and can propagate along flat metal-dielectric interfaces [21].Similar to THz radiation, SPs have also attracted much attention in recent years [19].
It is relatively less well known that SP can also exist at a cylindrical metal-dielectric interface [22][23][24].For an SP at a flat metal-dielectric interface, there are one magnetic field component H y , and two electrical field components E x and E z .The only transverse magnetic field component H y indicates the TM polarization.For an SP at a cylindrical metal-dielectric interface, however, there are one magnetic field component H φ , and two electric field components E r and E z .The only transverse magnetic field component H φ indicates the TM polarization.The latter kind of SP can be reasonably called APSP, because the only magnetic field component H φ is angular.
Consider the eigenproblem of a cylindrical metal wire surrounded by air.We denote by ε and µ the relative permittivity and the relative magnetic permeability, respectively.For the metal, we use ε m and µ m .And for air, ε=µ=1.We denote by λ 0 , c, and R the wavelength in free space, the light speed in free space, and the radius of the metal wire, respectively.We assume a temporal factor exp(-jωt), where ω=2πc/λ 0 , is the angular frequency.
We are interested in axially symmetrical eigenmodes.For them, the relations ∂E/∂φ=0 and ∂H/∂φ=0 hold in the cylindrical coordinates.By use of these relations, one can separate the electromagnetic field into two families.One is the TE polarization, and the other is the TM polarization.For the former, there are three components E φ , H r and H z .For the latter, there are another three components H φ , E r and E z .The two families of TE and TM polarizations are decoupling in this case.We now focus on the TM polarization.By substituting the relations ∂E/∂φ=0 and ∂H/∂φ=0 into Maxwell's equations [25], one obtains All of H φ , E r and E z are functions of the two variables r and z.For the formulation of the eigenproblem, they can be written as H φ (r,z)= H φ (r)exp(jk 0 n eff z), E r (r,z)= E r (r)exp(jk 0 n eff z), and E z (r,z)= E z (r)exp(jk 0 n eff z), respectively, where n eff is the effective index of the eigenmode.
Then the operator ∂/∂z can be replaced by jk 0 n eff .In the remainder of this paper, we refer to H φ , E r and E z specifically as H φ (r), E r (r) and E z (r), respectively.By using the relation ∂/∂z=jk 0 n eff in Eq. (3), one gets Substituting Eq. (4) into Eq.(1), and using the relation ∂/∂z=jk 0 n eff again, one can obtain Further substituting Eq. ( 5) into Eq.( 2), one can obtain where ρ=k 0 [(n eff ) 2 −µε] 1/2 r.The suitable solution for E z has the form I 0 (k 0 κ m r) in the metal and K 0 (k 0 κ a r) in air, because this solution approximately decays exponentially from the interface, where I 0 (.) and K 0 (.) are the generalized Bessel functions [26], . Substituting this kind of solution into Eq.( 5), and using the relations [26] dI 0 (k 0 κ m r)/dr=k 0 κ m I 1 (k 0 κ m r) and dK 0 (k 0 κ a r)/dr=−k 0 κ a K 1 (k 0 κ a r), one can determine the solution form of H φ .By using the continuities of E z and H φ at the interface r=R, one can obtain the following eigenvalue equation The existence of the solution of Eq. ( 7) results from the negative real part of ε m .This negative value is due to the electron plasma contribution when the frequency is lower than the plasma frequency.Accordingly, this eigenmode is traditionally called SP.In particular, when R→∞, the relations I 1 (k 0 κ m R)/I 0 (k 0 κ m R)=1 and K 1 (k 0 κ a R)/K 0 (k 0 κ a R)=1 hold.By use of these properties, one can obtain, for very large R, Eq. ( 8) is exactly valid for an SP at a flat metal-air interface.It can be expected that Eq. ( 8) is also approximately valid for an APSP at the interface of a very thick metal wire.In particular, for nonmagnetic metals with µ m =1, Eq. ( 8) reduces to n eff =[ε m /(1+ε m )] 1/2 .Eq. ( 8) is already enough for some qualitative analyses, though the exact n eff value of the APSP should be obtained by numerically solving Eq. ( 7).We now briefly discuss the dimensions of the APSP modal field.For simplicity, we discuss the only magnetic field component H φ (r).The H φ (r) field is proportional to K 1 (k 0 κ a r)/K 0 (k 0 κ a R) outside the metal.As we shall point out below, n eff is approximately equal to 1.As a consequence, κ a is very small.Because κ a is very small, then the field distribution K 1 (k 0 κ a r)/K 0 (k 0 κ a R) decays very slowly with the radial coordinate r.Typically, the modal field has a width of about several tens times of the radius R outside the metal.In the metal, the H φ (r) field is proportional to I 1 (k 0 κ m r)/I 0 (k 0 κ m R).As we shall see below, |ε m | is far larger than 1 in the THz spectral region.Using the relations |ε m |>>1 and n eff ≈1, one can get the relation κ m ≈(-ε m ) 1/2 .As a consequence, the real part of κ m is quite large.Typically, it has the order of magnitude of 1.0×10 3 or larger.This property leads to the fast decay of the modal field in the metal.Typically, the modal field can exist in the metal for only about 1 µm or less.Therefore, the modal field of an APSP has a very large beam width but with a hollow center.To get an intuitive impression, we calculate the H φ (r) field distribution of an APSP of a copper wire with a radius of 0.45 mm.The frequency is 0.5 THz (the corresponding wavelength is 0.6 mm).The related calaulation method used here will be stated in Section 3.
From Fig. 1(a), one can see that the modal field outside the metal has a long tail of about 20 mm.In contrast, the modal field decays very fast in the metal.From Fig. 1(b), one can see that the penetration depth in the metal is smaller than 0.5 µm.
We now attribute the observed waveguiding of a metal wire [17] to the propagation of an APSP along the wire.In the following, we try to use this model to explain the related experimental results of Ref. [17].
1).It is well known that the absolute values of ε m of metals are huge in the spectral region of THz radiation.This property can be deduced from the Drude model [25,27].Typically [27], −Re(ε m ) is larger than 10 4 and Im(ε m ) is larger than 10 5 , where Re(.) and Im(.) indicate the real and the imaginary parts, respectively.By using the properties |ε m |>>1 and |ε m |>>µ m in Eq. ( 8), one can find that n eff ≈1.This result explains the low loss of metal wire waveguides because Im(n eff ) is very small.
2).The relation n eff ≈1 holds for a very wide bandwidth, because |ε m | is huge in the whole range of THz radiation.This property explains the very low dispersion of metal wire waveguides.
3).It can be proved that the ration |E z /E r | is in the order of magnitude of |κ a /n eff | in air.By using the relation n eff ≈1, one can find that |E z |<<|E r |.This property implies that the longitudinal field component E z is negligible in air.As a result, the modal field is an approximate TEM field in air.Note that it is no longer true in the metal.In the metal, |E r |<<|E z |.It is worth mentioning that, it is not the radial field component E r but the longitudinal field component E z that leads to the attenuation of the eigenmode.4).Because the longitudinal field E z is negligible, the remaining only electric field component in air is E r .Then it looks like that the electric field is radially polarized.5).The electrical field of the light source is linearly polarized, but the electric field of the APSP is approximately radially polarized.The polarization mismatch between them leads to a very low coupling efficiency, though some additional setup was used to help the coupling.
6).As we state above, the modal field is quite large.As a result, the guided propagation can be easily coupled between two waveguides that are in contact with each other.7).As stated above, the modal field is slowly decaying outside the metal.As a result, almost all of the energy is propagated in air whose refractive index is 1.Also, n eff ≈1.These two properties imply that the Fresnel reflectivity at the end of metal wire is negligible.Therefore, the mode can propagate off the end of the waveguide.8).The metal wire is axially symmetric.For axially symmetric electromagnetic fields, TM and TE polarizations are decoupling.Thus, the radial polarization of the electric field can be maintained during propagation, even after propagating off the end of the metal wire.9).The beam width of the modal field is quite large because it decays slowly outside the metal.As a result, the waveguide can be only slightly curved.Otherwise, emission due to waveguide bending may happen.
Clearly, these excellent explanations confirm that the waveguide effect of a metal wire is indeed due to the propagation of an APSP.It is worth mentioning that, besides the zero-order SP (i.e., APSP) mentioned above, there also exist higher-order SPs [23].However, they are not really guided modes.In fact, those higher-order modes were called virtual radiative SPs in Ref. [23], in contrast to APSP, which is called real nonradiative SP in Ref. [23].Therefore, the existence of those virtual radiative SPs does not affect our analyses.

Further discussions
In Section 2 we have explained the experimental results of Ref. [17].In this Section we discuss some further problems.
As we state above, the low loss and the very low dispersion of an APSP come from the huge absolute values |ε m |.In the spectral region of THz radiation, different metals still have different |ε m | [27] though those values are all huge.For example, the value |ε m | of copper is more than 10 times higher than that of iron [27].Generally speaking, the larger the value |ε m |, the better the waveguide.In addition, the relative magnetic permeability µ m also plays a role.By use of a Taylor expansion of Eq. ( 8) for huge ε m , one can prove that n eff ≈1−µ m /(2ε m ).From this expression, one can deduce that, for the same ε m , a magnetic metal will increase the loss by µ m times.Therefore, the use of ferromagnetic metal [25,28], such as iron and nickel, should be avoided.We suggest the use of nonmagnetic metals with very large |ε m |.Copper, silver and gold are good candidates.Copper, of course, would be preferrable because it is much cheaper than silver and gold.We notice that stainless steel wires are used in the experiments of Ref. [17].Stainless steels can be magnetic or "nonmagnetic" [29, 30], depending on the kinds.In Ref. [17], there is no report about the µ m value of the used stainless steel.Therefore, it is difficult to concretely calculate the waveguide parameters of the stainless steel wires used there [17].However, the µ m value should be not smaller than 1 [30].Also, the |ε m | value of a stainless steel should be smaller than that of copper, because the |ε m | value of iron, which is the main chemical composition of stainless steel, is at least 10 times smaller than that of copper [27].Therefore, a copper wire waveguide is superior to a stainless steel wire waveguide.
As we state above, the exact n eff value can only be obtained by numerically solving the eigenvalue equation ( 7), though Eq. ( 8) is very useful for qualititative analyses.It is known that both the functions I 1 (.) and I 0 (.) approximately increase exponentially with the increase of the real part of the argument.Unfortunately, the real part of the argument k 0 κ m R is quite large in the investigated case.For example, in the example of Fig. 1, the argument k 0 κ m R is as large as 6.21×10 3 −j4.95×10 3 .Because the argument k 0 κ m R is quite large, numerical overflow will happen if one directly calculates the functions I 1 (k 0 κ m R) and I 0 (k 0 κ m R) that appear in Eq. ( 7).To solve this problem, we use the asymptotic expression [26] where ν is the order and u is the argument.By using this relation in the ratio I 1 (k 0 κ m R)/I 0 (k 0 κ m R), the diverging factor exp(k 0 κ m R) is removed and the overflow is overcome.Then the eigenvalue equation ( 7) can be solved by complex-root search methods.
Local search methods would be effective.We use the Newton method to solve the eigenvalue problem.Once obtaining the eigenvalue n eff , one can immediately obtain the complete modal field.To see the influence of the finite radius of the metal wire, we calculate the effective indices n eff of the APSPs at the interfaces of copper wires with different R. The frequency is chosen to be 0.5 THz (i.e., λ 0 =0.6 mm).The ε m value is determined to be ε m =−6.3×10 5 +j2.77×10 6 according to a fitted Drude formula for copper [27].The relation µ m =1 is used.The calculated n eff values, as a function of R, are shown in Fig. 2. One can see that, just as expected, n eff approach to the value determined by Eq. ( 8) for very large R.One can also see that, both the values of Im(n eff ) and Re(n eff )−1 increase with the decrease of the radius R. Thick metal wires are sometimes inconvenient for certain applications, but with lower attenuation.Therefore, one should balance between the size and the attenuation when designing a concrete metal wire waveguide.
To see the dependency on the frequency, we calculate the n eff values, as a function of frequency.The metal is chosen to be copper, and the radius R is chosen to be 0.45 mm.The frequency-dependent ε m of copper is obtained from the corresponding Drude formula of Ref. [27].The calculated n eff , as a function of frequency is shown in Fig. 3(a).By use of the relation α=k 0 Im(n eff ), one can further obtain the attenuate coefficient α.The calculated α value is shown in Fig. 3(b).One can see that the attenuation is smaller than 2×10 -3 cm -1 in the range of 0.1~1THz.This result explicitly shows the superiority of a copper wire waveguide.One may find that the attenuations shown in Fig. 3(b) here are smaller than the measured values of Ref. [17].Also, in our calculations, the attenuation increases with increasing frequency, but an opposite trend was observed in Ref. [17].These differences indicate that the dominant loss does not result from the attenuation of the eigenmode (i.e., APSP) itself but from other sources.We guess some other metal components that are not far away from the metal wire are probably the loss sources.As we state above, the modal field of an APSP is quite large.If some other metal components (especially magnetic metal components, note that optical tables are normally magnetic) are in the range of the modal field, then the loss will increase.In addition, through numerical calculations, we found that, for a certain metal wire, the width of the modal field increases with the decrease of the frequency.Then, the APSPs for low frequencies will have higher losses if some other metal components are not far from the metal wire, because larger modal fields have higher losses in this case.
As we state above, the beam width of an APSP is quite large, though the metal wire is thin.This property is important for the focusing application.One should be aware of the large beam width when choosing a focusing lens.If the lens is not large enough, the truncation of the modal field may happen and lead to the loss of a lot of energy.It is worth mentioning that APSPs are related to radially polarized ligh beams [31][32][33].Radially polarized light beams can be used for many applications, such as microscopy, optical tweezer, and optical data storage.Clearly, a metal wire waveguide also represents a viable method to generate radially polarized light beams with high purity.The drwaback is the low efficiency.

Conclusions
We have attributed the observed THz metal wire waveguide to the propagation of an APSP, and explained the related experimental results of Ref. [17].Some further problems, such as the choice of metal and the influence of the finite radius of metal wire, have been also discussed.In particular, we have predicted that the attenuation coefficient is smaller than 2×10 -3 cm -1 in the whole range of 0.1~1 THz if a copper wire with a radius of 0.45 mm is used.The work presented in this paper not only provides a theoretical basis for THz metal wire waveguides, but also hopefully further paves the way for a wide range of new applications for THz sensing and imaging.

Fig. 1 .
Fig. 1.Normalized modal field of an APSP of a copper wire.The red curve is the modal field outside the metal, and the blue curves are the modal field in the metal.(a) The total profile.(b)The detailed distribution of the very small penetration of the modal field in the metal.

Fig. 2 .
Fig. 2. Change of effective index with the radius R of metal wire.The red curves are the calculated results, and the black curves are the values given by Eq. (8).The dashed curves are Im(n eff ), and the solid curves are Re(n eff )−1.

Fig. 3 .
Fig. 3. (a) Change of effective index with the frequency.The red curve is Im(n eff ), and the black curve is Re(n eff )−1.(b) Change of attenuation with the frequency.