Linear to Radial Polarization Conversion in the Thz Domain Using a Passive System

Linear to radial polarization conversion in the THz domain using a passive system

T. Grosjean1, F. Baida1, R. Adam2, J-P. Guillet2, L. Billot1, P. Nouvel2, J. Torres2, A. Penarier2, D. Charraut1, L. Chusseau2 1Institut FEMTO-ST, Université de Franche-Comté, UMR 6174 CNRS, Département d’Optique P.M. Dufﬁeux, 16 route de Gray, 25030 Besançon cedex, France.

2Institut d’Électronique du Sud, Université Montpellier 2, UMR 5214 CNRS, Place E. Bataillon, 34095 Montpellier, France. [email protected]

Abstract: This paper addresses a passive system capable of converting a linearly polarized THz beam into a radially polarized one. This is obtained by extending to THz frequencies and waveguides an already proven concept based on mode selection in optical ﬁbers. The approach is validated at 0.1 THz owing to the realization of a prototype involving a circular waveguide and two tapers that exhibits a radially polarized beam at its output. By a simple homothetic size reduction, the system can be easily adapted to higher THz frequencies. © 2008 Optical Society of America

OCIS codes: (220.4830) Systems design; (230.5440) Polarization-selective devices; (230.7370) Waveguides; (260.5430) Polarization

References and links 1. F. C. De Lucia, Sensing with Terahertz radiation, chap. Spectroscopy in the Terahertz spectral region, pp. 39–115 (Springer Series of Optical Science, Springer, Berlin, 2003). 2. A. Markelz, A. Roitberg, and E. Heilweil, “Pulsed terahertz spectroscopy of DNA, bovine serum albumin and collagen between 0.1 and 2 THz,” Chem. Phys. Lett. 320, 42–48 (2000). 3. S. Quabis, R. Dorn, M. Eberler, O. Glöckl, and G. Leuchs, “Focusing Light to a Tighter Spot,” Opt. Commun. 179, 1–7 (2000). 4. K. Youngworth and T. Brown, “Focusing of high numerical aperture cylindrical-vector beams,” Opt. Express 7, 77–87 (2000). 5. L. Helseth, “Roles of Polarization, Phase and Amplitude in Solid Immersion Lens Systems,” Opt. Commun. 191, 161–172 (2001). 6. C. Sheppard and A. Choudhury, “Annular pupils, radial polarization, and superresolution,” Appl. Opt. 43, 4322– 4327 (2004). 7. T. Grosjean and D. Courjon, “Smallest focal spots,” Opt. Commun. 272, 314–319 (2007). 8. L. Novotny, M. Beversluis, K. Youngworth, and T. Brown, “Longitudinal field modes probed by single molecules,” Phys. Rev. Lett. 86, 5251 (2001). 9. A. Bouhelier, J. Renger, M. Beversluis, and L. Novotny, “Plasmon–coupled tip–enhanced near–field optical microscopy,” J. Microsc. 210, 220–224 (2002). 10. E. Descrovi, L. Vaccaro, W. Nakagawa, L. Aeschimann, U. Staufer, and H. Herzig, “Collection of transverse and longitudinal fields by means of apertureless nanoprobes with different metal coating characteristics,” Appl. Phys. Lett. 85, 5340–5342 (2004). 11. V. Niziev and A. Nesterov, “Influence of beam polarization on laser cutting efficiency,” J. Phys. D: Appl. Phys. 32, 1455–1561 (1999).

#97435 - $15.00 USD Received 19 Jun 2008; revised 15 Sep 2008; accepted 25 Sep 2008; published 3 Nov 2008 (C) 2008 OSA 10 November 2008 / Vol. 16, No. 23 / OPTICS EXPRESS 18895 12. Q. Zhan, “Trapping metallic Rayleigh particles with radial polarization,” Opt. Express 12, 3377–3382 (2004). 13. J. Fontana and R. Pantell, “A high-energy, laser accelerator for electrons using the inverse Cherenkov effect,” J. Appl. Phys. 54, 4285–4288 (1983). 14. M. Roth, E. Wyss, H. Glur, and H. Weber, “Generation of radially polarized beams in a Nd:YAG laser with self-adaptive overcompensation of the thermal lens,” Opt. Lett. 30, 1665–1667 (2005). 15. I. Moshe, S. Jackel, and A. Meir, “Production of of radially or azimuthally polarized beams in solid-state lasers and the elimination of thermally induced birefringence effects,” Opt. Lett. 28, 807–809 (2003). 16. R. Dorn, S. Quabis, and G. Leuchs, “Sharper focus for a radially polarized light beam,” Phys. Rev. Lett. 91, 233,901 (2003). 17. K. Wang and D. Mittelman, “Metal wires for terahertz waveguiding,” Nature 432, 373–379 (2004). 18. J. Deibel, K. Wang, M. Escarra, and D. Mittelman, “Enhanced coupling of terahertz radiation to cylindrical wire waveguides,” Opt. Express 14, 279–290 (2006). 19. Q. Cao and J. Jahns, “Azimuthally polarized surface plasmons as effective terahertz waveguides,” Opt. Express 13, 511–518 (2005). 20. G. Chang, C. Divin, C.-H. Liu, S. Williamson, A. Galvanauskas, and T. Norris, “Generation of radially polarized terahertz pulses via velocity-mismatched optical rectification,” Opt. Lett. 32, 433–435 (2007). 21. B. Knoll and F. Keilmann, “Near-field probing of vibrational absorption for chemical microscopy,” Nature 399, 134–137 (1999). 22. K. Wang, D. Mittleman, N. van der Valk, and P. Planken, “Antenna effects in terahertz apertureless near-field optical microscopy,” Appl. Phys. Lett. 85, 2715–2717 (2004). 23. N. van der Valk and P. Planken, “Electro-optic detection of subwavelength terahertz spot sizes in the near-field of a metal tip,” Appl. Phys. Lett. 81, 1558–1560 (2002). 24. H.-T. Chen, R. Kersting, and G. Cho, “Terahertz imaging with nanometer resolution,” Appl. Phys. Lett. 83, 3009–3011 (2003). 25. H.-T. Chen, S. Kraatz, G. C. Cho, and R. Kersting, “Identification of a resonant imaging process in apertureless near-field microscopy,” Phys. Rev. Lett. 93, 267,401 (2004). 26. R. Kersting, H.-T. Chen, N. Karpowicz, and G. C. Cho, “Terahertz microscopy with submicrometre resolution,” J. Opt. A: Pure and Applied Optics 7, S184–S189 (2005). 27. R. Lecaque, S. Grésillon, N. Barbey, R. Perreti, J.-C. Rivoal, and A.-C. Boccara, “THz near-field optical imaging by a local source,” Opt. Commun. 262, 125–128 (2006). 28. Y. Mushiake, K. Matsumura, and N. Nakajima, “Generation of radially polarized optical beam mode by laser oscillation,” Proc. IEEE 60, 1107–1109 (1972). 29. J.-L. Li, K.-I. Ueda, M. Musha, A. Shirakawa, and L.-X. Zhong, “Generation of radially polarized mode in Yb fiber laser by using dual conical prism,” Opt. Lett. 31, 2969–2971 (2006). 30. Y. Kozawa and S. Sato, “Generation of a radially polarized laser beam by use of a conical Brewster prism,” Opt. Lett. 30, 3063–3065 (2005). 31. R. Oron, S. Blit, N. Davidson, A. Friesem, Z. Bomzon, and E. Hasman, “The Formation of Laser Beams with Pure Azimuthal or Radial Polarization,” Appl. Phys. Lett. 77, 3322–3324 (2000). 32. D. Pohl, “Operation of a Ruby Laser in the Purely Transverse Electric Mode TE01,” Appl. Phys. Lett. 20, 266– 267 (1972). 33. K. Yonezawa, Y. Kozawa, and S. Sato, “Generation of a radially polarized laser beam by use of the birefringence of a c-cut Nd:YVO4 crystal,” Opt. Lett. 31, 2151–2153 (2006). 34. A. Nesterov, V. Niziev, and V. Yakunin, “Generation of High-Power Radially Polarized Beam,” J. Phys. D: Appl. Phys. 32, 2871–2875 (1999). 35. T. Moser, H. Glur, V. Romano, M. Ahmed, F. Pigeon, O. Parriaux, and T. Graf, “Polarization-selective grating mirrors used in the generation of radial polarization,” Appl. Phys. B 80, 707–713 (2005). 36. S. Tidwell, D. Ford, and W. Kimura, “Generating radially polarized beams interferometrically,” Appl. Opt. 29, 2234–2239 (1990). 37. S. Tidwell, G. Kim, and W. Kimura, “Efficient radially polarized laser beam generation with a double interfer- ometer,” Appl. Opt. 32, 5222–5229 (1993). 38. N. Passilly, D. de Saint Denis, K. Aït-Ameur, F. Treussart, R. Hierle, and J.-F. Roch, “Simple interferometric technique for generation of a radially polarized light beam,” J. Opt. Soc. Am. A 22, 984–991 (2005). 39. E. Churin, J. Hosfeld, and T. Tschudi, “Polarization Configurations with Singular Point Formed by Computer Generated Holograms,” Opt. Commun. 99, 13–17 (1993). 40. M. Stalder and M. Schadt, “Linearly Polarized Light with Axial Symmetry Generated by Liquid-Crystal Polar- ization Converters,” Opt. Lett. 21, 1948–1950 (1996). 41. M. Neil, F. Massoumian, R. Juskaitis, and T. Wilson, “Method for the generation of arbitrary complex vector wave front,” Opt. Lett. 27, 1929–1931 (2002). 42. R. Hongwen, L. Yi-Hsin, and W. Shin-Tson, “Linear to axial or radial polarization conversion using a liquid crystal gel,” Appl. Phys. Lett. 89, 051,114 (2006). 43. K. Moh, X.-C. Yuan, D. Tang, W. Cheong, and L. Zhang, “Generation of femtosecond optical vortices using a single refractive optical element,” Appl. Phys. Lett. 88, 091,103 (2006).

#97435 - $15.00 USD Received 19 Jun 2008; revised 15 Sep 2008; accepted 25 Sep 2008; published 3 Nov 2008 (C) 2008 OSA 10 November 2008 / Vol. 16, No. 23 / OPTICS EXPRESS 18896 44. T. Grosjean, D. Courjon, and M. Spajer, “An All–Fiber Device for Generating Radially and Other Polarized Light Beams,” Opt. Commun. 203, 1–5 (2002). 45. G. Volpe and D. Petrov, “Generation of cylindrical vector beams with few-mode fibers by Laguerre-Gaussian beams,” Opt. Commun. 237, 89–95 (2004). 46. T. Grosjean, A. Sabac, and D. Courjon, “A versatile and stable device allowing the efficient generation of beams with radial, azimuthal or hybrid polarizations,” Opt. Commun. 252, 12–21 (2005). 47. M. Ordal, R. Bell, R. Alexander, L. Newquist, and M. Querry, “Optical properties of Al, Fe, Ti, Ta, W, and Mo at submillimeter wavelengths,” Appl. Opt. 27, 1203–1209 (1988). 48. M. Ordal, R. J. Bell, R. Alexander, L. Long, and M. R. Querry, “Optical properties of Au, Ni, and Pb at submillimeter wavelengths,” Appl. Opt. 26, 744–752 (1987). 49. A. Azad, Y. Zhao, W. Zhang, and M. He, “Effect of dielectric properties of metals on terahertz transmission through subwavelength hole arrays,” Opt. Lett. 31, 2637–2639 (2006). 50. D. Marcuse, Light Transmission Optics, chap. 8, pp. 293–294 (New York: Van Nostrand Reinhold, 1972). 51. B. Davidson and R. Ziolkowski, “Body–of–revolution finite–difference time–domain Modeling of space–time Focusing by a three–dimensional Lens,” J. Opt. Soc. Am. A 11, 1471–1490 (1994). 52. F. Baida, D. Labeke, and Y. Pagani, “Body-of-Revolution FDTD Simulations of Improved Tip Performance for Scanning Near-Field Optical Microscopes,” Opt. Commun. 255, 241–252 (2003). 53. R. Oron, Y. Danziger, N. Davidson, A. Friesem, and E. Hasman, “Discontinuous phase elements for transverse mode selection in laser resonators,” Appl. Phys. Lett. 74, 1373–1375 (1999). 54. F. Baida, D. Labeke, G. Granet, A. Moreau, and A. Belkhir, “Origin of the super-enhanced light transmission through a 2-D metallic annular aperture array: a study of photonic bands,” Appl. Phys. B 79, 1–8 (2004). 55. F. Baida, A. Belkhir, D. Labeke, and O. Lamrous, “Subwavelength metallic coaxial waveguides in the optical range: Role of the plasmonic modes,” Phys. Rev. B 74, 205,419 (2006). 56. N. Klein, P. Lahl, U. Poppe, F. Kadlec, and P. Kuzel, “A metal-dielectric antenna for terahertz near-field imaging,” J. Appl. Phys. 98, 014910 (2005).

1. Introduction For twenty years, considerable improvements have been brought to THz systems and the field has seen an astonishing development of its experimental tools in the inspection of materials and the spectroscopy of chemical species [1] and biological samples [2]. As experienced in the optical domain, the generation and use of radially polarized beams should have a significant impact in the evolution of research in the THz domain, especially in near-field microscopy (and derived techniques) and THz plasmonics. At optical frequencies, these highly symmetric fields have brought an important added value in high numerical aperture focusing and conventional microscopy [3-7], single molecule probing [8], near-field optical microscopy and spectroscopy [9, 10], laser cutting [11], particle trapping and accelerating [12, 13], high-power laser emission “stability” [14, 15]. All these studies take benefit from the total field symmetry provided by radial polarization. Such a property can lead, for example, to focused fields with either a strong longitudinal component [3, 4, 16] or a p-polarization state with symmetry of revolution. So far, only a very limited number of systems allows the generation of radially polarized THz beams. Most of them are based on the excitation of the mode sustained by metallic wires whose electric field distribution is radially distributed around the wire [17-19]. The use of pho- toconductive antennas with radial symmetry seems to solve the problem of coupling efficiency between THz radiation and the wire mode [18]. However, the radiation efficiency of the wire mode in free space is questionable. Recently, it has been shown that optical rectification gives an alternative solution for the generation of THz radially polarized pulses [20]. However, both systems require the design of a proper active THz source. In this paper, we propose a new concept of passive radial polarizer at THz frequencies. Such a system is aimed at working with continuous waves and it is specially designed for radiating in free space a radially polarized doughnut mode. It has the advantage to be adaptable to any kind of existing THz sources. In a first section, the principle of the proposed THz radial polarizer is described. Then, the concept is validated with a first prototype.

#97435 - $15.00 USD Received 19 Jun 2008; revised 15 Sep 2008; accepted 25 Sep 2008; published 3 Nov 2008 (C) 2008 OSA 10 November 2008 / Vol. 16, No. 23 / OPTICS EXPRESS 18897 2. Principle of the passive THz radial polarizer 2.1. From Optics to terahertz The wide development of the terahertz techniques is partly due to the successful transposition of some optical concepts to the THz domain. For example, the extension to THz frequencies of the well-known near-field optical microscopy has opened new perspectives in the local observation of physical and biological materials with submicrometer resolution [21-27]. We suggest here that the generation of THz radially polarized beams can also get benefits from the know-how of the optics community. The various techniques proposed for more than 15 years for generating radially polarized light can be divided into active and passive systems. Active systems are generally conventional lasers whose cavity has been modified for radiating a radially polarized mode. These modifi- cations consist of inserting optical components such as axicons [28-30], phase-step elements [31], birefringent crystals [32, 33] or diffractive elements [34, 35] into the cavity in order to se- lect the desired mode. Despite good results in terms of polarization purities and/or efficiencies, these solutions do not appear to be well-adapted to the current THz sources. Passive systems are designed for converting linearly or circularly polarized laser beams into doughnut beams that are radially polarized. They are set outside the laser cavity and do not produce by them- selves THz waves. Among the passive systems developed so far, one can find Mach-Zehnder interferometers which combine the orthogonal TEM01 and TEM10 modes [36-38], computer generated holograms [39], liquid crystal polarization converters [40-42], spiral phase system [43] and the mode selection inside few-mode optical fibers [44-46]. All these concepts seem to be easier to extend to the THz domain as they are not limited to a specific radiation source. The challenge is here to extend to THz frequencies the components specifically developed for the optical range.

2.2. Modes sustained by circular metallic waveguides It turns out that the optical fibers used in Refs. [44-46] can find in metallic circular waveguides a straightforward counterpart for the generation of radial polarization in the THz domain. Figure 1 displays the classification of the allowed propagating modes of a perfectly conducting metallic cylindrical waveguide with respect to its cutoff diameter.

TE11 Fundamental TE01 mode TM01 TE21 TM11 TE31 Modes

0.59 0.77 0.97 1.22 1.34 Cutoff diameter Waveguide (*l) diameter

Fig. 1. Classiﬁcation of the ﬁrst modes carried by a hollow circular waveguide of perfect metal (λ is the wavelength).

Following in wavelength, scale after the linearly polarized fundamental TE11 mode, are the radially polarized TM01 and the hybrid polarized TE21! and TE21!! modes. Figure 2 shows the intensity distribution and electric ﬁeld orientation of these ﬁrst modes into a 7 λ wide circular metallic waveguide (λ is the wavelength). The primes on the TE21 mode indicate that one is rotated by 45° with respect to the other. The twofold degeneracy of the TE21 mode is a consequence of cylindrical symmetry.

Fig. 2. Intensity distribution and electric ﬁeld orientation (pointed out by arrows) of the TE11, TM01 and TE21 modes sustained by a 7 λ wide metallic circular waveguide.

2.3. Coupling efficiency between the waveguide modes and an incoming Gaussian beam The generation of radially polarized beams with metallic cylindrical waveguides requires first, the selective excitation of the TM01 mode and second, the efficient emission of this mode in free space. Since the concept of radial polarizer has to be independent of the configurations of THz continuous-wave sources and wavelength, it is assumed in the following that the incoming beam is a linearly polarized Gaussian beam whose beam-waist W is about 10 λ. Such field confinements can be obviously achieved for example with parabolic mirrors. We also as- sume that the waveguides are made with infinitely conducting perfect metal. Perfect metal is a good approximation of the real metals of interest over the THz spectral range. For example, in the case of aluminum, Drude model leads to complex permittivity comprised between values around 3.62104 + 7.32106i at 0.1 THz and 1.62104 + 5.85104i at 10 THz. In this portion − − of the spectrum, the skin depth does not exceed λ/1230. Therefore, the penetration of the THz radiation into the metal is weak enough so that the transmission loss of the guided modes can be neglected over the distance of propagation of a few tens wavelengths considered here. Let us note that this approximation is also valid for other metals as iron, gold and silver which exhibit skin depths smaller than the ones of aluminum in the THz domain [47-49]. Generally, the energy transfer from the incident beam and the various waveguide modes is investigated by means of the coupling coefficients Cm defined as

Pm Cm = , (1) Pi

where Pm and Pi are the powers carried by the modes and the input beam, respectively. When the incoming beams and the waveguides are much larger than the wavelength, the longitudinal components of the various fields are negligible. In that case, the powers Pm and Pi can be seen as the Poynting vector flow through the waveguide cross section 1 P = ℜ r dr dθ(E H∗) e , (2) i 2 i × i · z !! 1 P = ℜ r dr dθ(a E b∗ H∗ ) e . (3) m 2 m m × m m · z !! Unit vector ez defines the direction of the waveguide axis. Constants am and bm are the components of respectively the input electric and magnetic field distributions onto the m–mode. Their expressions are based on the overlap integrals between the incoming and the mode field

r dr dθ(Ei Hm∗ ) ez am = × · , (4) r dr dθ(E H ) e "" m × m∗ · z "" r dr dθ(Em∗ Hi) ez bm = × · . (5) r dr dθ(E H ) e "" m∗ × m · z "" The expression of the field distribution (Em,Hm) can be easily calculated by adapting the procedure given in Ref. [50] to the hollow cylindrical waveguide made in perfect metal. In the following, the axis of symmetry of the input Gaussian beam is supposed to be aligned and centered with respect to the cylindrical waveguide. Assuming that (Ei,Hi) defines the input beam waist, we have 1 E (x,y,0,t)= exp( iωt) e (u,v)G(u,v)exp[i(ux + vy)] dudv, (6) i 4 2 − i π !! 1 H (x,y,0,t)= exp( iωt) h (u,v)G(u,v)exp[i(ux + vy)] dudv. (7) i 4 2 − i π !! 2 (u2+v2)W 2 where G(u,v) = πW exp − , and h =(k e )/(ωµ ) with k =(u,v,w) being the 4 i × i 0 wave vector, ω the wave# angular frequency$ and µ0 the permeability of the vacuum. Expres- sions ei(u,v)G(u,v) and hi(u,v)G(u,v) define the plane wave spectrum of the input electric and magnetic fields, respectively. Figure 3(a) reports the coupling between the incident linearly polarized Gaussian incoming beam and the cylindrical waveguide modes as a function of the waveguide diameter. In this basic configuration, only the fundamental mode TE11 that is linearly polarized is excited with a maximum efficiency of 0.87 for a waveguide diameter of 13 λ.

2.4. Selection of the radial mode of a cylindrical metallic waveguide

Our approach for selectively exciting the upper mode TM01 radially polarized is divided in two steps. The first one, borrowed from Ref. [46], consists of inverting the polarization direction of the input beam (configuration 1) or the fundamental waveguide mode (configuration 2) over half their cross sections. Then, the two halves are in phase opposition. Note that in configuration 1, this polarization inversion is carried out before the injection of the input beam inside the waveguide, whereas in configuration 2, it is realized inside the waveguide on the linearly polarized fundamental mode. Figures 3(b) and 3(c) show the coupling efficiencies between the Gaussian beam and the waveguide modes in configurations 1 and 2, respectively. In both cases, the energy of the fundamental mode that was solely excited without field reversal as shown in Figure 3(a), is now transferred into modes TM01 and TE21 (see Figs. 3(b) and 3(c)). The maximum coupling efficiency between the incident gaussian beam and the waveguide mode TM01 reaches 44% in configuration 2 (at a waveguide diameter of 13 λ) whereas it is limited to 28% in configuration 1 (at a waveguide diameter of 10.3 λ). The second step of the TM01 mode selection consists of filtering out the TE21 mode. This task can be fulfilled by channeling the TM01 and the TE21 modes excited in configurations 1 or 2 into a smaller hollow cylinder whose diameter is comprised between the cutoff diameters of the two modes (cf Fig. 1). In that case, only TM01 (with the lower cutoff diameter) is transmitted through the small waveguide, TE21 being completely reflected. This operation can be carried out by tapering the large waveguide so that its output aperture matches the entrance aperture of the small one (see inset of Fig. 4). Figure 4 displays the channeling efficiency of mode TM01 as a function of the taper angle θ. This channeling efficiency is defined as the ratio between the energies carried by the TM01 mode inside the small and the large waveguides. BOR-FDTD (Body-Of-Revolution Finite Difference

#97435 - $15.00 USD Received 19 Jun 2008; revised 15 Sep 2008; accepted 25 Sep 2008; published 3 Nov 2008 (C) 2008 OSA 10 November 2008 / Vol. 16, No. 23 / OPTICS EXPRESS 18900 0.9 Incident gaussian beam TE11 0.85

0.8

C 0.75

0.7

Other modes : C=0 0.65 8 l 10 l 12 l 14 l 16 l 18 l 20 l waveguide diameter

(a)

0.35 0.5 Incident gaussian Fundamental beam mode 0.45 0.3 0.4 TM TM01 01

0.25 0.35

C 0.3 0.2 TE21

0.25 TE21 0.15 0.2

Other modes : C=0 Other modes : C=0 0.1 0.15 8 l 10 l 12 l 14 l 16 l 18 l 20 l 8 l 10 l 12 l 14 l 16 l 18 l 20 l waveguide diameter waveguide diameter

(b) (c)

Fig. 3. (a) Coupling efficiency C between a linearly polarized Gaussian beam and the various modes sustained by a cylindrical hollow waveguide of perfect metal, as a function of the waveguide diameter. (b) Configuration 1, half the incoming Gaussian beam cross- section has been phase retarded by π (see inset). (c) Configuration 2, half the fundamental mode cross section initially excited has been phase retarded by π (as shown in the inset).

Time Domain) algorithm [51, 52] with spatial uniform mesh has been used to simulate the propagation of the mode inside the tapered waveguide. The transmission has been calculated as the normalized transmitted Poynting vector flux by the incident one. The flux in each case is determined by integrating the Poynting vector over the cylinder sections. In the case reported here, the diameters of the large and the small hollow cylinders are 7 λ and 0.9 λ, respectively. As shown in Fig. 1, the small circular waveguide fulfills the condition required for transmitting TM01 and filters out TE21. The transmission efficiency of the TM01 mode through the system keeps lower than 5% for taper angles smaller than 62 degrees but grows rapidly for larger angles and almost reaches 90% for θ = 83◦ which is considered as a quasi adiabatic regime of the taper. ≈

2.5. Efficiency of the linear to radial polarization conversion From the study detailed above, the efficiency of the linear to radial polarization conversion can be evaluated by multiplying the efficiencies calculated previously for the two steps of the TM01 mode selection. Figure 5 shows the total efficiency as a function of the large waveguide diameter, for a width of the small waveguide of 0.9 λ and a taper angle of 80◦. For this angle,

0.8

TM +TE 0.7 01 21 _ TM q 01 0.6

0.5

0.4

0.3

Power transmission 0.2

0.1

0 0 10 20 30 40 50 60 70 80 90 Taper angle q (degree)

Fig. 4. Channeling efﬁciency of the TM01 mode as a function of the taper angle θ. The diameters of the large and small waveguides are 7 λ and 0.9 λ, respectively.

the channeling efﬁciency, calculated for a diameter of 7 λ, is about 77%. This value has been kept constant for calculating the efﬁciency of the radial polarizer.

0.4

0.35

0.3

0.25

0.2

Efficiency

0.15

0.1

0.05 8 l 10 l 12 l 14 l 16 l 18 l 20 l waveguide diameter

Fig. 5. Efficiency of the TM01 mode selection as a function of the diameter of the first (large) waveguide; dashed line: configuration 1; solid line: configuration 2, configurations are detailed in §2.4.

The second configuration of polarization reversal leads to a maximum efficiency of linear to radial polarization conversion larger than 34% for a large waveguide diameter of 13 λ whereas the maximum efficiency of the first one is limited to 22% at large waveguide diameter of 10.3 λ. Therefore, the configuration of radial polarizer which involves a polarization reversal inside the first (large) waveguide seems to be better adapted for radially polarizing THz waves.

2.6. Far-ﬁeld emission

Because the TM01 mode frequency is close to the cutoff inside the small waveguide, the free space radiation efﬁciency of the polarizing device is not maximum. This problem can be solved by placing a circular horn antenna at the end of the small waveguide for adapting the electro- magnetic impedance of the output TM01 guided mode to vacuum.

#97435 - $15.00 USD Received 19 Jun 2008; revised 15 Sep 2008; accepted 25 Sep 2008; published 3 Nov 2008 (C) 2008 OSA 10 November 2008 / Vol. 16, No. 23 / OPTICS EXPRESS 18902 2.7. Prototype design and fabrication Figure 6 depicts the preliminary system aimed at validating the above presented concept of radial polarizer at a frequency of 0.1 THz. It consists of a tapered waveguide coupled to a discontinuous phase element (DPE). The latter is a dielectric plate whose output interface exhibits a step. The step height h = λ/2(n 1) induces a π-phase retardation for the thicker part of the plate with respect to the thinner one− when a radiation goes through the system at normal inci- dence. When the step is orthogonal to the direction of the incident linear polarization, the DPE induces the desired phase conversion for the generation of radial polarization (Fig. 7). Such a dielectric system is widely used in optics for the generation of radially polarized beams [53, 31, 38, 46]. DPE Waveguide system, Al Teflon n f=100Ghz

40mm l/(2(n-1))

Fig. 6. Schema of the ﬁrst radial polarizer prototype. It is composed of a DPE and a focusing

waveguide system. _

DPE

_ TE11

TM01+TE21

Fig. 7. Effect of the DPE onto the incident free space Gaussian beam. The waveguide modes excited with and without DPE are also indicated.

The DPE is fabricated by micromachining a 3.5 mm high step onto one of the two faces of a polytetrafluoroethylene (PTFE) plate whose index of refraction n is close to 1.43 at 0.1 THz. The DPE is set outside the waveguide system in order to finely center the phase step with respect to the waveguide entrance aperture thanks to a precision translation stage. As discussed in §2.4, the fine positioning of the DPE outside the waveguide system is made at the expense of some loss of efficiency. The focusing waveguide system has been fabricated by micromachining an aluminum rod. At its heart, a 6 mm long hollow cylindrical waveguide of 2.7 mm (0.9 λ) diameter supports the

#97435 - $15.00 USD Received 19 Jun 2008; revised 15 Sep 2008; accepted 25 Sep 2008; published 3 Nov 2008 (C) 2008 OSA 10 November 2008 / Vol. 16, No. 23 / OPTICS EXPRESS 18903 simultaneous propagation of both the TE11 and TM01 modes. Input is tapered with a circular horn up to an external aperture of diameter 40 mm for channelling an incoming 30 mm ≈ (10 λ) wide Gaussian beam to the small waveguide. Adiabatic conversion of the free-space propagating mode is ensured by a taper angle θ higher than 80◦ of resulting attenuation 1.1 dB as a result of calculations of Fig. 4. The output of the waveguide is also tapered but to≤ a limited diameter of 4 mm. This circular horn antenna ensures a sufﬁcient power density at the output aperture with an enhanced free space radiation of the system. Relatively high intensities at the system output are required for the detection of the radially polarized beam with low efﬁcient micro-probes.

2.8. Simulation of the prototype A preliminary numerical study of the ability of this prototype to polarize radially THz radiations is conducted by using a home-made BOR-FDTD code. This latter is based on the discretiza- tion of the Maxwell equations when they are expressed in cylindrical coordinates. By this way, the axis-symmetry of the structure is fully analytically treated to reduce the dimension of the studied problem by explicitly expressing the azimuthal dependence in Maxwell equations. This leads to consider a 2D meshing in the (r,z) space instead of a 3D one for the three cylindrical coordinates (r,z,φ). Because aluminum shows very high permittivities in the THz domain [47], transmission losses inside the metallic structure are negligible and the prototype has been modelled with perfect metal for convenience. In Fig. 8, the transmission through the waveguide system is evaluated for two incoming free space eigenmodes with (Fig. 8(a)) radial and (Fig. 8(c)) hybrid polarizations. These two modes are supposed to be coupled to the TM01 and TE21 modes inside the metallic structure, respectively. In both cases, the injection plane is located at 20 mm in front of the entrance facet of the device. The field distributions of the input beams in the injection plane are the ones of the TM01 (Fig. 8(a)) and TE21 (Fig. 8(c)) guided mode of a cylindrical waveguide at the desired frequency. These initial conditions are calculated using a N-order FDTD code [54] especially elaborated for the determination of the eigenfrequencies of an axis-symmetrical structure [55]. Figures (8(d) and (f)) show that the radially polarized incident beam (coupled to mode TM01) is transmitted through the system (Fig. 8(d)) whereas the incoming beam with hybrid polarization (which couples to mode TE21) is stopped at the entrance of the hollow cylinder waveguide (Fig. 8(f)). For a better view of this phenomenon, the fifth root of the electric intensities are represented in Fig. 8(d,e,f). Note that the length of the cylinder is larger than the evanescent tail of TE21 inside the waveguide. This explains the negligible output far-field re-emission of this mode (Fig. 8(i)). The transmission of the degenerated mode that is usually produced by the DPE (Figs. 8(b,e,h)) can be simulated simply by adding the field distributions calculated in the two last cases. We see that the two lobe incident field distribution (Fig. 8(b)) is converted into a doughnut output beam that is related to the selective transmission of the TM01 mode radially polarized through the metallic structure (Fig. 8(h)). An efficiency of 11% has been calculated for this preliminary configuration of waveguide system. Note that this value do not take into account the reflections of the input beam onto the interfaces of the DPE.

3. Experimental setup The experimental setup of generation and characterization of the radially polarized beam is sketched in Fig. 9. A continuous-wave linearly polarized radiation is generated at 0.1 THz. The source is an electronic synthesizer followed by a sextupler from Spacek Labs Inc. Available operating bandwidth of this commercial system is 75-110 GHz with a nearly constant 0 dBm output. A WR10 horn antenna is connected to the source for improving the free space radiation efﬁciency. After their collection and collimation with a parabolic mirror, waves propagate under

Input intensity 31.25 mm 31.25 mm 31.25 mm +

20 mm TM01 20 mm TE21

20 mm

Sagittal sections (intensity) (d) (e) (f)

(g) (h) (i) -5 X10

Output intensity 15.2 mm 15.2 mm 15.2 mm

Fig. 8. FDTD simulation of the focusing waveguide system. The real case is simulated in the middle column whereas the right and left columns show the projection of the ﬁeld distribution in a basis of eigenmodes. We see that the degenerated space mode (b) that is produced with the DPE is the result of the combination of (a) a radially polarized mode and (c) a four-spot mode with hybrid polarization. (d,e,f) show the ﬁfth root of the electric intensities in a longitudinal cross-section of the device for the three input modes. (g,h,i) exhibit the electric intensities in a lateral plane located at 15 mm from the output side of the device.

#97435 - $15.00 USD Received 19 Jun 2008; revised 15 Sep 2008; accepted 25 Sep 2008; published 3 Nov 2008 (C) 2008 OSA 10 November 2008 / Vol. 16, No. 23 / OPTICS EXPRESS 18905 a form close to a Gaussian beam whose waist is about 30 mm (10 λ), as shown in Fig. 10(a).

Fig. 9. Scheme of the experimental setup.

The polarization converter is placed into the THz beam so that the DPE is set close to the beam waist and the waveguide system is positioned in contact to the DPE. The waveguide is carefully centered and aligned with respect to the incoming beam propagation axis. The DPE can be translated transversally with a 1D translation stage that ensures an accurate centering of the DPE phase step with respect to the waveguide input aperture. Detection involves a Schottky diode coupled to a PTFE pyramidal probe of principle similar to the one described in Ref. [56]. The probe was built from a WR10 waveguide (inner sizes 1.27 mm 2.54 mm) ﬁlled with a PTFE parallelepiped ended by a pyramid of 3 mm height protruding from× the waveguide end. The two opposite largest facets of the tip were metal coated using Ni pulverization. The sensitive end facet area exhibits dimensions of 20 40 µm. The detection provided by this tip is thus specially engineered to be polarization sensitive× since this collection system picks up only one transverse component of the electric ﬁeld. The distance between the polarizer end facet and the probe is kept lower than 0.2 mm (λ/15) during image acquisitions by raster-scanning. The probe system is mounted onto a 2D motorized translation stage, with optimal resolution of 100 nm, for scanning the end aperture of our device.

4. Results and discussion Figure 10 shows the properties of the field before and after its transmission by the waveguide system when the DPE is removed. Images in Fig. 10(b) reports the intensity of the field components parallel and perpendicular to the incident polarization direction as observed at the waveguide system end. Those images have been taken by two successive acquisitions, with a scan step of 200 µm and by rotating the probing device by 90◦ between the two scans. As a major consequence,≈ the upper part of the Fig. 10(b) displays a bright spot whereas the lower part shows an intensity distribution which just exceeds the detection background. These images demonstrate that the output field distribution is linearly polarized along the polarization direction of the incoming beam. Such a polarization property is due to the excitation of the TE11 mode inside the waveguide structure. The acquisitions realized with the DPE in front of the focusing system are reported in Fig. 11. Images obtained when the axis of the polarizing probe stage is set parallel and perpendicular to the incoming polarization direction are displayed in Figs. 11(a) and 11(b), respectively. These orientations are indicated by white arrows on each image. A numerical reconstruction of the outgoing beam from Figs. 11(a) and 11(b) is provided in Fig. 11(c) whereas cross-sections of Fig. 11(a) along the horizontal direction (solid curve), and Fig. 11(b) along the vertical direction

(b)

Fig. 10. (a) Measured intensity before the DPE (dots) compared with a theoretical gaussian beam (solid line). (b) Images of the transmitted intensity obtained without DPE when the polarizing probe axis is parallel (upper part) and perpendicular (lower part) to the incident polarization direction. Intensities are normalized to the same maximum value for both images.

(dashed curve) are plotted in Fig. 11(d). Two-grain structures jump out in Figs. 11(a) and 11(b). As expected, their directions follow the prescribed axis of the polarizing micro-detection system. The numerical combination of these two orthogonal patterns leads to an annular shape intensity distribution. The null intensity at the beam center (see Fig. 11(d)) evidences that the fundamental TE11 mode, whose maximum is expected at the center, has been totally rejected by the structure. Moreover, the visibility of the two-spots pattern remains unchanged when the probe axis is rotated (Figs. 11(a) and 11(b)). This is another evidence that all higher modes, except TM01, are reflected by the system. These observations validate our concept of THz radial polarizer. Note that the efficiency of this prototype cannot be measured with precision by means of the detection tools that we used. This is partly due to the fact that there is no direct coupling between a radially polarized beam and the fundamental mode of a rectangular waveguide. From the discussion of §2.4, it can be enhanced by inserting the DPE inside the focusing structure. This can be achieved, for example, by adapting a cylindrical waveguide to the entrance of the focusing taper. This will be made in a close future for fabricating the first THz experiments involving radial polarization.

(d)

Fig. 11. Acquisition results of the ﬁeld distribution transmitted by the prototype over a scan of 7 7 mm2. (a) and (b) are images acquired for two orthogonal axis of the polarizing × detection probe (axis indicated by arrows). (c) Numerical combination of (a) and (b). (d) Horizontal cross-section of (a) (solid curve) and vertical one of (b) (dashed curve).

#97435 - $15.00 USD Received 19 Jun 2008; revised 15 Sep 2008; accepted 25 Sep 2008; published 3 Nov 2008 (C) 2008 OSA 10 November 2008 / Vol. 16, No. 23 / OPTICS EXPRESS 18908 5. Conclusion A new THz radial polarizer suitable to convert the usual linearly polarized signal from a millimeter or sub-millimeter source is described and demonstrated at 100 GHz. Its principle of operation involves an adequate mode filtering in a metallic cylindrical≈ waveguide that supports only TE11 and TM01 propagating modes. The correct system operation arises with the selection of the radially polarized TM01 mode that is ensured by means of a discontinuous phase element placed at the entrance of the system. A nearly optimum design of such a passive radial polarizer has been given owing to numerical FDTD calculations of the coupling efficiency of a Gaussian propagating linearly polarized beam to the TM01 mode of a large circular waveguide. Aperture size of this input cylindrical waveguide as well as the design of the taper that follows have been optimized. The experimental realization has been conducted at 0.1 THz. The built polarizer has shown an operation in excellent agreement with the theoretical design. Doughnut modes have been observed at polarizer output with a very high rejection of the fundamental TE11 mode. Provided that micromechanical machining difficulties can be overcome, the proposed design is straightforwardly scalable to much higher frequencies in the whole THz domain. In the future, this polarizer coupled to an axicon is aimed at generating very small focal spots with enhanced longitudinal electric field for THz near field imaging purposes.

Acknowledgments Authors are indebted to the Agence Nationale de la Recherche for supporting this research under contract #ANR06-BLAN-0073. They are grateful to D. Courjon for helpful discussions. One of us (RA) also thanks jointly the CNRS and the Région Languedoc-Roussillon for doc- toral fellowship.