Abstract
The evanescent wave of the cylindrical vector field is analyzed using the vector angular spectrum of the electromagnetic beam. Comparison between the contributions of the TE and TM terms of both the propagating and the evanescent waves associated with the cylindrical vector field in free space is demonstrated. The physical pictures of the evanescent wave and the propagating wave are well illustrated from the vectorial structure, which provides a new approach to manipulating laser beams by choosing the states of polarization in the cross-section of the field.
© 2013 Optical Society of America
1. Introduction
Recently the subject of vector optical field has attracted much interest in linear and nonlinear optics realms due to its unique and novel properties and potential applications [1–3]. In contrast to most of the past research based on spatially homogeneous polarization states such as linear, elliptical, and circular polarization, the vector optical fields deal with inhomogeneous states of polarization in the field cross section. The vector optical fields have been studied intensively. For instance, focusing a radially polarized beam can lead to the smallest possible spot size [4, 5], and the angular optical momentum from the vector optical field has been demonstrated [6]. Evanescent waves, also known as inhomogeneous waves [7], have a long history and have been the subject of many studies and appear in many branches of physics [8]. The evanescent waves of nonparaxial beams have been extensively studied [8–13], due to the novel properties and potential applications in many areas, such as nanophotonics, guided optics, internal reflection spectroscopy, and optical microscopy [8]. Generation and applications of evanescent wave with different systems and instruments have been reported [14–16]. R. Merlin and his colleagues described the contribution of pure evanescent waves to focused fields for generating a small focal spot beyond diffraction limit [17, 18]. More recently, the evanescent waves of radially polarized and azimuthally polarized optical fields have been demonstrated [19–21].
In this paper, we study the evanescent wave of an arbitrary vector optical field with inhomogeneous states of polarization in the cross section of the field. The TE and TM terms (i.e. transverse electric mode (TE) term and transverse magnetic mode (TM) term in term of its vector angular spectrum) of the evanescent wave and the propagating wave of the cylindrical vector optical field are derived by the vector angular spectrum method. Propagation of the wavefield in z direction in free space is described by the propagation of the angular spectrum in the near field as well as in the far field since the angular spectrum propagation method is a solution of the Helmholtz equation [22, 23]. Two other ways to characterize such propagation are in terms of the Fresnel and the Fraunhofer approximation, but they are valid only under certain constraints. Recently, the focused azimuthal-variant vector beams have been studied under the nonparaxial and paraxial approximations by using the vectorial Rayleigh-Sommerfeld formulae [24] which can also be derived from the angular method [23]. Here in this paper, the ratio of the evanescent wave and the propagating wave of a cylindrical vector optical field as a function of the propagation distance in the near field is investigated. The intensity (squared modulus) distributions of the TE and TM terms of the propagating and the evanescent waves are analyzed in detail. These results, therefore, provide useful information on how to spatially manipulate the evanescent waves of a vector cylindrical optical field in the near field by choosing appropriate states of polarization in the cross-section of the field.
2. Theoretical formulation
In the Cartesian coordinate system, the z-axis is taken to be the direction of wave propagation. A cylindrical optical vector field can be described as [1–3]
where and are the polar radius and the azimuthal angle in the polar coordinate system, respectively. is the unit vector describing the distribution of the states of polarization of the vector field. m is the topological charge, and θ0 is the initial phase. and are the unit vectors in x and y directions respectively in the Cartesian coordinate system. A(r) represents the amplitude distribution in the cross-section of the cylindrical vector field. The time dependent factor is omitted in Eq. (1) ( is the circular frequency). When m = 1 with θ0 = 0 and π/2, the vector fields correspond to the radially and azimuthally polarized vector fields, respectively. When m = 0, Eq. (1) degenerates to the horizontal (for θ0 = 0) and vertical (for θ0 = π/2) linearly polarized fields, respectively.The transverse components of the vector angular spectrum and are given by the Fourier transform of the input field:
where k is the wavenumber. ρ and ϕ are related to the transverse Fourier-transform variables u, ν and ξ, and we haveFor Gaussian distributionwith an arbitrary polarized electromagnetic field as Eq. (1) in the source plane z = 0, where w indicates the waist of Gaussian beam, the angular spectrum iswhere I(.) are the Bessel functions of the second kind and is the unit vector in z direction. The electric field component of the vector cylindrical optical field in the z plane can be represented as:The electromagnetic field component of the vector cylindrical optical field in the z plane can be decomposed accordingly into the TE and TM terms [25–30]:Both the TE and TM components are perpendicular to the propagation axis. With reference to [13] and [20], they can be written in the following forms: And the corresponding evanescent wave components can be obtained according to the calculations similar in Refs [19, 21]: where J(.) are the Bessel functions of the first kind.In order to compare the contributions of the propagating and the evanescent waves associated with the cylindrical vector field, the integrated intensity (squared modulus) of the propagating and the evanescent fields, Ipr and Iev, are calculated respectively [19, 21]:
with Substituting Eqs. (11)-(13) into Eqs. (9) and (10), and performing integrations over ø from zero to 2π yield a single integration of ρ. Then, the corresponding results can be obtained by performing numerical integration over ρ. The ratio δ = (Ipr - Iev)/Ipr provides the direct information about the propagating and the evanescent components of the field.The ratio δ = (Ipr - Iev)/Ipr for w = 0.1λ (highly nonparaxial case) and w = 0.5λ cases as a function of different distances z from the initial plane z = 0 are shown in Fig. 1. It can be clearly seen that the evanescent field dominates near the source plane and the relative weight of Iev would drastically decrease with the increase of the propagation distance. It becomes negligible in the propagation distance z = 0.5λ as shown in Fig. 1. Comparing Figs. 1(a) with 1(b), one can recognize that the relative weight of Iev would reduce when the waist size w increases, i.e., the relative weight of the evanescent wave component is increasing with increasing nonparaxial feature as the waist size w decreases. As it is expected, bigger values of the waist size w leads to reduced nonparaxial behavior. The transverse and the longitudinal components of the evanescent wave become negligible with larger waist size w and hence the relative weight of Iev would be smaller. Furthermore, it can also be found from Fig. 1 that the relative weight of Iev would increase with the increase of the topological charge. It can be explained by the fact that the field distribution will increasingly diverge and extend from the center of beam with the increasing topological charge. As a result, the relative weight of Iev would increase with the increasing topological charge under the same conditions and beam parameters. The calculation results indicate that the evanescent wave component with different initial phase θ0 for any topological charge (except m = 1) is the same as shown in Fig. 1. For m = 1 case with different initial phase θ0, the relative weight of Iev of the azimuthal polarization (m = 1, θ = π/2) is the maximum and the relative weight of Iev of the radial polarization (m = 1, θ = 0) is the minimum as shown in Fig. 1. The physical explanation for the exception m = 1 is that the values of the z-components of either the propagation wave or the evanescent wave are different for different initial phase θ0, and especially the values of the z-components of either the propagation wave or the evanescent wave are the maximum for θ0 = π/2 (i.e., radial polarization) whereas the z-components of either propagation wave or evanescent wave are zero for θ0 = 0 (i.e., azimuthal polarization), which can also be recognized from Eqs. (5)-(8).
Once we have discussed the appropriate ranges for the beam waist w and the propagation distance z, we will intuitively show the contribution of the evanescent and the propagation waves to the optical field. Examples of the intensity distributions of the TE and TM components for the evanescent and the propagation waves are depicted in Figs. 2, 3, and 4. The value of w is set to be 0.1λ for the highly nonparaxial case. The plane z = 0.2λ is selected as the reference plane. θ is considered to be π/4 for m = 0, 1, and 2, respectively. Apparently, the TE and TM terms are orthogonal to each other in the field. The intensity distribution profiles of the TE and TM components of the evanescent and the propagation terms are similar, but the magnitude of intensity for the propagation term is smaller than that of the evanescent term. This verifies that the evanescent field dominates near the source plane and the contribution of the evanescent wave to the cylindrical vector field is considerable in magnitude. From Fig. 2 and Fig. 4, for m = 0 and m = 2 with w = 0.1λ and propagation distance z = 0.2λ, it can be found that the intensity distribution profiles of both the evanescent and the propagation waves are composed of a pattern with two peaks and some pairs of side lobes. Moreover, the intensities of the side lobes are much lower than those of the peaks. In these cases, the magnitude of the TM term is greater than that of the TE term. As analyzed above, the magnitude of the intensities of the evanescent wave and its two terms (TE and TM) decrease with the increase of the waist width w and the propagation distance z. The two peaks will superimpose as a Gaussian distribution. The numerical results indicate that the two peaks will superimpose as a Gaussian distribution for m = 0 and the propagation distance z = 0.2λ when w = 0.5λ as described in Ref [30]. For m = 2, it needs farther propagation distance or greater waist width (e.g., w = 3λ for z = 0.2λ). Figure 3 illustrates the intensity distributions of the propagating and the evanescent components of the cylindrical vector beam when m = 1, θ0 = π/4, w = 0.1λ, z = 0.2λ. These patterns resemble the doughnut shape with a dark spot in the center and some bright rings around it. Such a field has important applications in optical trapping and manipulating of nano-objects [3]. We also present the contributions of the evanescent and the propagation waves to the radially and the azimuthally polarized optical fields in Fig. 5. The TE components of both the evanescent and the propagation waves of the radial polarization optical field are zero as the TM components of the azimuthal polarization.
3. Conclusion
By using the angular spectrum method, the Maxwell’s equation is resolved. The TE and TM terms of the propagating wave and the evanescent wave components of nonparaxial Gaussian beam with inhomogeneous polarization in the cross-section of the field are obtained in the integral form. The contribution of the evanescent wave to the vector optical field is analyzed with the numerical results. The relative weight of the evanescent wave would increase with the increasing topological charge but would reduce with the increasing waist size w and drastically decrease with the increasing propagation distance. The intensity distribution profiles of the evanescent term are similar to those of the propagating term, but the magnitude of the intensities of the propagating term is different from that of the evanescent term. The intensity distribution is different with the increasing waist width w or propagation distance z.
Acknowledgments
G. Li would like to thank the support from the OSU Ophthalmology and ECE departments, the ElectroScience Laboratory CERF fund, National Institutes of Health National Eye Institute (through grant R01 EY020641), National Institute of Biomedical Imaging and Bioengineering (through grant R21 EB008857), National Institute of General Medical Sciences (through grant R21 RR026254/R21 GM103439), and Wallace H. Coulter Foundation Career Award (through grant WCF0086TN). G. Li also thanks the reviewers for their constructive comments and M. Yun for some editing work.
References and links
1. Q. W. Zhan, “Cylindrical vector beams: from mathematical concepts to applications,” Adv. Opt. Photon. 1(1), 1–57 (2009). [CrossRef]
2. X. L. Wang, Y. N. Li, J. Chen, C. S. Guo, J. P. Ding, and H. T. Wang, “A new type of vector fields with hybrid states of polarization,” Opt. Express 18(10), 10786–10795 (2010). [CrossRef] [PubMed]
3. B. Gu, F. Ye, K. Lou, Y. Li, J. Chen, and H. T. Wang, “Vectorial self-diffraction effect in optically Kerr medium,” Opt. Express 20(1), 149–157 (2012). [CrossRef] [PubMed]
4. Q. Zhan and J. Leger, “Focus shaping using cylindrical vector beams,” Opt. Express 10(7), 324–331 (2002). [CrossRef] [PubMed]
5. R. Dorn, S. Quabis, and G. Leuchs, “Sharper Focus for a Radially Polarized Light Beam,” Phys. Rev. Lett. 91(23), 233901 (2003). [CrossRef] [PubMed]
6. X. L. Wang, J. Chen, Y. Li, J. Ding, C. S. Guo, and H. T. Wang, “Optical orbital angular momentum from the curl of polarization,” Phys. Rev. Lett. 105(25), 253602 (2010). [CrossRef] [PubMed]
7. G. C. Sherman, J. J. Stamnes, and E. Lalor, “Asymptotic approximations to angular-spectrum representations,” J. Math. Phys. 17(5), 760–776 (1976). [CrossRef]
8. F. De Fornel, Evanescent Waves: From Newtonian Optics to Atomic Optics (Springer. 2001).
9. G. P. Agrawal and M. Lax, “Free-space wave propagation beyond the paraxial approximation,” Phys. Rev. A 27(3), 1693–1695 (1983). [CrossRef]
10. C. J. R. Sheppard and S. Saghafi, “Electromagnetic Gaussian beams beyond the paraxial approximation,” J. Opt. Soc. Am. A 16(6), 1381–1386 (1999). [CrossRef]
11. A. Doicu and T. Wriedt, “Plane wave spectrum of electromagnetic beams,” Opt. Commun. 136(1-2), 114–124 (1997). [CrossRef]
12. D. G. Hall, “Vector-beam solutions of Maxwell’s wave equation,” Opt. Lett. 21(1), 9–11 (1996). [CrossRef] [PubMed]
13. A. V. Shchegrov and P. S. Carney, “Far-field contribution to the electromagnetic Green’s tensor from evanescent modes,” J. Opt. Soc. Am. A 16, 2583–2584 (1999). [CrossRef]
14. M. F. Imani and A. Grbic, “Generating evanescent Bessel beams using near-field plates,” IEEE Trans. Antenn. Propag. 60(7), 3155–3164 (2012). [CrossRef]
15. A. Grbic, R. Merlin, E. M. Thomas, and M. F. Imani, “Near-field plates: metamaterial surfaces/arrays for subwavelength focusing and probing,” Proc. IEEE 99(10), 1806–1815 (2011). [CrossRef]
16. M. F. Imani and A. Grbic, “Tailoring near-field patterns with concentrically corrugated plates,” Appl. Phys. Lett. 95(11), 111107 (2009). [CrossRef]
17. R. Merlin, “Radiationless electromagnetic interference: evanescent-field lenses and perfect focusing,” Science 317(5840), 927–929 (2007). [CrossRef] [PubMed]
18. A. Grbic, L. Jiang, and R. Merlin, “Near-field plates: subdiffraction focusing with patterned surfaces,” Science 320(5875), 511–513 (2008). [CrossRef] [PubMed]
19. L. E. Helseth, “Radiationless electromagnetic interference shaping of evanescent cylindrical vector waves,” Phys. Rev. A 78(1), 013819 (2008). [CrossRef]
20. R. Martínez-Herrero, P. M. Mejías, I. Juvells, and A. Carnicer, “Transverse and longitudinal components of the propagating and evanescent waves associated to radially polarized nonparaxial fields,” Appl. Phys. B 106(1), 151–159 (2012). [CrossRef]
21. R. Martínez-Herrero, P. M. Mejías, and A. Carnicer, “Evanescent field of vectorial highly non-paraxial beams,” Opt. Express 16(5), 2845–2858 (2008). [CrossRef] [PubMed]
22. J. W. Goodman, Introduction to Fourier Optics. (Greenwood Village: Roberts and Company, 2004).
23. K. E. Okan, Diffraction, Fourier Optics and Imaging. (Wiley & Sons, 2007).
24. B. Gu and Y. Cui, “Nonparaxial and paraxial focusing of azimuthal-variant vector beams,” Opt. Express 20(16), 17684–17694 (2012). [CrossRef] [PubMed]
25. R. Martínez-Herrero, P. M. Mejías, S. Bosch, and A. Carnicer, “Vectorial structure of nonparaxial electromagnetic beams,” J. Opt. Soc. Am. A 18(7), 1678–1680 (2001). [CrossRef] [PubMed]
26. H. Guo, J. Chen, and S. Zhuang, “Vector plane wave spectrum of an arbitrary polarized electromagnetic wave,” Opt. Express 14(6), 2095–2100 (2006). [CrossRef] [PubMed]
27. G. Zhou, “Analytical vectorial structure of Laguerre-Gaussian beam in the far field,” Opt. Lett. 31(17), 2616–2618 (2006). [CrossRef] [PubMed]
28. K. Duan and B. Lü, “Polarization properties of vectorial nonparaxial Gaussian beams in the far field,” Opt. Lett. 30(3), 308–310 (2005). [CrossRef] [PubMed]
29. R. P. Chen and K. H. Chew, “Far field properties of a vortex Airy beam,” Laser and Particle Beams 31(01), 9–15 (2013). [CrossRef]
30. G. Zhou, “The analytical vectorial structure of a nonparaxial Gaussian beam close to the source,” Opt. Express 16(6), 3504–3514 (2008). [CrossRef] [PubMed]