Focal shift in vector Mathieu-Gauss beams

doi:10.1364/OE.16.005838

« Show journal navigation

Focal shift in vector Mathieu-Gauss beams

Raul I. Hernandez-Aranda and Julio C. Gutiérrez-Vega »View Author Affiliations

Optics Express, Vol. 16, Issue 8, pp. 5838-5848 (2008)
http://dx.doi.org/10.1364/OE.16.005838

View Full Text Article

Acrobat PDF (481 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Vector Mathieu-Gauss beams
3. Numerical results
4. Conclusions
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Multimedia
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (33)
Cited By
Figures (5)
Multimedia (8)
Metrics

Abstract

The three dimensional distribution of focused vector Mathieu-Gauss beams (vMG) is studied in the vicinity of the geometrical focus of an unapertured thin lens. We adopt two different intensity based criteria for defining the actual focus. Our analysis confirms the existence of a focal shift towards the lens for this type of beams. The dependence of the focal shift on the different parameters of the beams is discussed in detail. Beams with different states of polarization are studied as well, and it is shown that the focal shift is independent of the polarization.

1. Introduction

Over the past decades, considerable attention has been paid to the phenomenon of focal shift, in which the maximum on-axis intensity for a monochromatic converging wave does not occur exactly at the geometric focus, but it is slightly displaced towards the aperture plane [1

E. Collet and E. Wolf, “Symmetry properties of focused fields,” Opt. Lett. 5, 264–266 (1980). [CrossRef]

, 2

J. J. Stamnes and B. Spjelkavik, “Focusing at small angular apertures in the Debye and Kirchhoff approximations,” Opt. Commun. 40, 81–85 (1981). [CrossRef]

, 3

Y. Li and E. Wolf, “Focal shifts in diffracted converging spherical waves,” Opt. Commun. 39, 211–215 (1981). [CrossRef]

, 4

Y. Li, “Dependence of the focal shift on Fresnel number and f number,” J. Opt. Soc. Am. 72, 770–774 (1982). [CrossRef]

, 5

W. H. Carter, “Focal shift and concept of effective Fresnel number of a Gaussian laser beam,” Appl. Opt. 21, 1989–1994 (1982). [CrossRef] [PubMed]

, 6

V. N. Mahajan, “Axial irradiance and optimum focusing of laser beams,” Appl. Opt. 19, 3042–3053 (1983). [CrossRef]

, 7

M. Martínez-Corral, C. J. Zapata Rodríguez, P. Andrés, and E. Silvestre, “Effective Fresnel-Number concept for evaluating the relative focal shift in focused beams,” J. Opt. Soc. Am. A. 15, 449–455 (1998). [CrossRef]

, 8

C. J. R. Sheppard and P. Török, “Dependence of focal shift on Fresnel number and angular aperture,” Opt. Lett. 23, 1803–1804 (1998). [CrossRef]

, 9

S. De Nicola, D. Anderson, and M. Lisak, “Focal shift effects in diffracted focused beams,” Pure Appl. Opt. 7, 1249–1259 (1998). [CrossRef]

]. It has been shown by different authors that this shift depends on the Fresnel number of the Gaussian apodization, for unapertured beams, and on the effective Fresnel number associated with the aperture, for apertured systems [1

E. Collet and E. Wolf, “Symmetry properties of focused fields,” Opt. Lett. 5, 264–266 (1980). [CrossRef]

, 5

W. H. Carter, “Focal shift and concept of effective Fresnel number of a Gaussian laser beam,” Appl. Opt. 21, 1989–1994 (1982). [CrossRef] [PubMed]

, 7

In imaging applications where the Fresnel number is usually large, the focal shift is negligible, however in other physical systems with Fresnel numbers closer to unity, or even smaller, the effect cannot be neglected. This phenomenon has been extensively studied in the context of coherent scalar beams, for instance it has been shown to exist in systems with an obscured pupil [6

V. N. Mahajan, “Axial irradiance and optimum focusing of laser beams,” Appl. Opt. 19, 3042–3053 (1983). [CrossRef]

, 10

T. C. Poon, “Focal shift in focused annular beams,” Opt. Commun. 65, 401–406 (1988). [CrossRef]

], in Laguerre-Gaussian beams [11

W. H. Carter and M. F. Aburdene, “Focal shift in Laguerre-Gaussian beams,” J. Opt. Soc. Am. A 10, 1949–1952 (1987). [CrossRef]

], in unapertured Bessel-Gauss beams [12

B. Lü and W. Huang, “Focal shift in unapertured Bessel-Gauss beams,” Opt. Commun. 109, 43–46 (1994). [CrossRef]

] and apertured Bessel beams [13

S. De Nicola, “On axis focal shift effects in focused truncated J ₀ Bessel beams,” Pure Appl. Opt. 5, 827–831 (1996). [CrossRef]

]. More recently, Hricha et al. [14

Z. Hricha, L. Dalil-Essakali, and A. Belafhal, “Axial intensity distribution and focal shifts of focused partially coherent conical Bessel-Gauss beams,” Opt. Quantum. Electron. 35, 101–110 (2003). [CrossRef]

] and Dong and Pu [15

M. Dong and J. Pu, “Effective Fresnel number and the focal shifts of focused partially coherent beams,” J. Opt. Soc. Am. A 24, 192–196 (2007). [CrossRef]

] have studied this effect in partially coherent beams. However, to our knowledge very few references have addressed this effect for vector beams [16

P. L. Greene and D. G. Hall, “Focal shift in vector beams,” Opt. Express 4, 411–419 (1999). [CrossRef] [PubMed]

, 17

J. Pu and B. Lü, “Focal shifts in focused nonuniformly polarized beams,” J. Opt. Soc. Am. A 18, 2760–2766 (2001). [CrossRef]

On the other hand, Mathieu beams, which constitute a family of nondiffracting optical beams [18

J. C. Gutiérrez-Vega, M. D. Iturbe-Castillo, and S. Chávez-Cerda, “Alternative formulation for invariant optical fields,” Opt. Lett. 25, 1493–1495 (2000). [CrossRef]

], have been employed in different areas of optics, for instance in photonic lattices to study soliton behavior in periodic media [19

Y. V. Kartashov, A. A. Egorov, V. A. Vysloukh, and L. Torner, “Shaping soliton properties in Mathieu lattices,” Opt. Lett. 31, 238–240 (2006). [CrossRef] [PubMed]

], in optical tweezers for the transfer of orbital angular momentum [20

C. López-Mariscal, J. C. Gutiérrez-Vega, G. Milne, and K. Dholakia, “Orbital angular momentum transfer in helical Mathieu beams,” Opt. Express. 14, 4182–4187 (2006). [CrossRef] [PubMed]

], and in the study of localized X waves [21

C. A. Dartora and H. E. Hernández-Figueroa, “Properties of a localized Mathieu pulse,” J. Opt. Soc. Am. A 21, 662–667 (2004). [CrossRef]

, 22

J. Davila-Rodriguez and J. C. Gutiérrez-Vega, “Helical Mathieu and parabolic localized pulses,” J. Opt. Soc. Am. A 24, 3449–3455 (2007). [CrossRef]

]. In this paper, we study the focal shift in vector Mathieu-Gauss beams (vMG), which constitute a particular family of the more general vector Helmholtz-Gauss beams (vHzG) [23

M. A. Bandres and J. C. Gutiérrez-Vega, “Vector Helmholtz-Gauss and vector Laplace-Gauss beams,” Opt. Lett. 30, 2155–2157 (2005). [CrossRef] [PubMed]

]. We start our discussion by introducing the vMG beams and setting up the mathematical formulae. Next, we discuss about the criteria for defining the focal distance of a vMG beam, and confirm the existence of a focal shift for these beams based on this criteria. The numerical results are presented and the dependence of the focal shift on the beam parameters is discussed in detail.

2. Vector Mathieu-Gauss beams

Vector Mathieu-Gauss (vMG) beams constitute a particular family of the more general vector Helmholtz-Gauss beams, introduced by Bandres and Gutiérrez-Vega [23

M. A. Bandres and J. C. Gutiérrez-Vega, “Vector Helmholtz-Gauss and vector Laplace-Gauss beams,” Opt. Lett. 30, 2155–2157 (2005). [CrossRef] [PubMed]

], and correspond to localized solutions of the paraxial wave equation in elliptical coordinates. These solutions represent a monochromatic field propagating along the positive z axis, and have a transverse distribution whose amplitude is modulated by a Gaussian factor. Recently we have explored the propagation dynamics of focused vector Mathieu-Gauss beams in free space [24

R. I. Hernández-Aranda, M. A. Bandres, and J. C. Gutiérrez Vega, “Propagation dynamics of vector Mathieu-Gauss beams,” in Laser Beam Shaping VII, Fred M. Dickey and David L. Shealy, eds., Proc. SPIE 6290, 629011 (2006). [CrossRef]

], and through paraxial ABCD systems [25

R. I. Hernandez-Aranda, J. C. Gutiérrez-Vega, M. Guizar-Sicairos, and M. A. Bandres, “Propagation of generalized vector Helmholtz-Gauss beams through paraxial optical systems,” Opt. Express. 14, 8974–8988 (2006). [CrossRef] [PubMed]

], giving closed form expressions for their evaluation. Chafiq et al. [26

A. Chafiq, Z. Hricha, and A. Belafhal, “Propagation properties of vector Mathieu-Gauss beams,” Opt. Commun. 275, 165–169 (2007). [CrossRef]

] have also studied the propagation properties of vMG beams.

2.1. General description of a vector Mathieu-Gauss beam

Consider a monochromatic electromagnetic field propagating along the positive z axis. The electric and magnetic components of this field can be expressed in terms of their transverse and longitudinal components as E=(E _t+ẑ E_z )exp(ikz) and H=(H _t+ẑ H_z )exp(ikz) respectively, where

k = ω {(μ_{0} ε_{0})}^{\frac{1}{2}} = \sqrt{k_{t}^{2} + k_{z}^{2}}

is the wave number, and the subscripts t and z refer to transverse and longitudinal components. Throughout the rest of the paper a time dependence factor exp(-iωt) will be assumed, and it will be omitted from the equations for the sake of simplicity. For a vHzG beam the transverse component of the electric field E _t is of the form [23

M. A. Bandres and J. C. Gutiérrez-Vega, “Vector Helmholtz-Gauss and vector Laplace-Gauss beams,” Opt. Lett. 30, 2155–2157 (2005). [CrossRef] [PubMed]

]

E_{t} = Z (ζ) G (r, ζ) U (X, Y),

(1)

where

r = \sqrt{x^{2} + y^{2}}

is the usual radial variable in polar coordinates, (X,Y)=(x/ζ,y/ζ) are scaled Cartesian coordinates, and ζ(z)=1+iz/z_R , with z_R being the Rayleigh distance. The function Z(ζ)=exp[k ² _t w ² ₀(ζ ^-1-1)/4] depends only on the propagation coordinate, with w ₀ being the Gaussian beam waist at plane z=0. The factor G(r,ζ) represents a Gaussian amplitude modulation given by

G (r, ζ) = ζ^{- 1} \exp (\frac{{- r}^{2}}{w_{0}^{2} ζ}) .

(2)

Finally, the function U is known to be a solution to the vector Helmholtz equation in 2D, and following Stratton [23

M. A. Bandres and J. C. Gutiérrez-Vega, “Vector Helmholtz-Gauss and vector Laplace-Gauss beams,” Opt. Lett. 30, 2155–2157 (2005). [CrossRef] [PubMed]

, 27

J. A. Stratton, Electromagnetic Theory (McGraw-Hill, New York, 1941).

], it represents two independent solutions corresponding to the TM and TE polarizations. These solutions are of the form:

U^{(1)} = \nabla_{T} W (X, Y), U^{(2)} = \hat{z} \times U^{(1)} .

(3)

It is important to mention that the gradient ∇_T is taken over the transverse scaled coordinates. Once the electric field in Eq. (1) is determined, its corresponding magnetic field can be computed as H _t=(ε ₀/µ ₀)^1/2 ẑ×E _t, which is a direct consequence of the transversality of the fields from Maxwell’s equations.

Inserting the solution U ⁽¹⁾ from Eq. (3) into Eq. (1) we obtain the following expressions for the electric and magnetic fields of the TM polarization

E_{t} = Z (ζ) G (r, ζ) \nabla_{T} W (X, Y),

(4)

H_{t} = \sqrt{\frac{ε_{0}}{μ_{0}}} Z (ζ) G (r) (\hat{z} {\times \nabla}_{T} W (X, Y)) .

(5)

A similar procedure but now employing U ⁽²⁾ will give the expressions corresponding to TE polarization. Based on the Lax series expansion of Maxwell’s equations [28

M. Lax, W. H. Louisell, and W. B. McKnight, “From Maxwell to paraxial wave optics,” Phys. Rev. A 11, 1365–1370 (1975). [CrossRef]

], E _t and H _t would correspond to the zeroth order terms, and for the next-order of correction the presence of small longitudinal field components is required. These longitudinal terms can be computed from the previous knowledge of E _t and H _t as follows

E_{z} = \frac{i}{k} \nabla_{T} \cdot E_{t}, H_{z} = \frac{i}{k} \nabla_{T} \cdot H_{t} .

(6)

It is important to remark that these solutions hold for the paraxial approximation, where the z component of the field is small compared to the transverse component. The longitudinal components of vMG beams include a term of the order (k_t /k), and therefore to satisfy the paraxial condition it is required that k≫k_t , i.e. the transverse spatial oscillations are many wavelengths wide.

For vMG beams, the function W(X,Y) in Eq. (3) is a scalar solution to the 2D Helmholtz equation in elliptical coordinates, we shall refer to it as the seed function, and may be expressed in terms of the m-th order radial and angular Mathieu functions [18

J. C. Gutiérrez-Vega, M. D. Iturbe-Castillo, and S. Chávez-Cerda, “Alternative formulation for invariant optical fields,” Opt. Lett. 25, 1493–1495 (2000). [CrossRef]

, 29

J. C. Gutiérrez-Vega and M. A. Bandres, “Helmholtz-Gauss waves,” J. Opt. Soc. Am. A 22, 289–298 (2005). [CrossRef]

]

W_{m}^{e} (ξ, η; q) = {Je}_{m} (ξ, q) {ce}_{m} (η, q),

(7)

W_{m}^{o} (ξ, η; q) = {Jo}_{m} (ξ, q) {se}_{m} (η, q),

(8)

where Je_m(·),Jo_m(·) represent even and odd radial Mathieu functions, and ce_m(·), se_m(·) correspond to the even and odd angular Mathieu functions respectively. The (e,o) superscripts stand for even an odd solutions, m is a non-negative integer indicating the order of the functions, and q is a free parameter called the ellipticity parameter. The scaled variables (ξ,η) in the elliptical coordinates system are related to the scaled Cartesian coordinates by the relation (X,Y)=(h coshξ cos η,h sinh ξ sin η), where h is the semifocal distance of the elliptical coordinates at the plane z=0, and is also related to the transverse component of the wave vector k_t , and to the ellipticity parameter q by h=2√q/k_t .

Since the transverse structure of vMG beams is primarily determined by the gradient of the function W, as can be seen from Eq. (3), we provide an expression for the gradient in the scaled elliptic coordinates

\nabla_{T} W (ξ, η) = (\frac{1}{h_{ξ}} \frac{\partial}{\partial ξ} \hat{ξ} + \frac{1}{h_{η}} \frac{\partial}{\partial η} \hat{η}) W (ξ, η),

(9)

where the scaling factors are given by

h_{ξ} = h_{η} = h \sqrt{\sinh^{2} ξ + \sin^{2} η}

. This gradient can also be expressed in Cartesian coordinates by an appropriate transformation

\nabla_{T} W (ξ, η) = \frac{\partial}{\partial X} W (ξ, η) \hat{x} \frac{\partial}{\partial Y} W (ξ, η) \hat{y},

(10)

\frac{\partial}{\partial X} W (ξ, η) = \frac{h \sinh ξ \cos η}{h_{ξ}^{2}} \frac{\partial}{\partial ξ} W (ξ, η) - \frac{h \cos ξ \sin η}{h_{η}^{2}} \frac{\partial}{\partial η} W (ξ, η),

(11)

\frac{\partial}{\partial Y} W (ξ, η) = \frac{h \cosh ξ \sin η}{h_{ξ}^{2}} \frac{\partial}{\partial ξ} W (ξ, η) + \frac{h \sinh ξ \cos η}{h_{η}^{2}} \frac{\partial}{\partial η} W (ξ, η),

(12)

doing so is helpful for numerical purposes, since it directly provides expressions for the x̂ and ŷ components.

2.2. Focusing of vector Mathieu-Gauss beams

We now consider the propagation of a vMG beam through an unapertured thin lens of focal distance f. Figure 1 shows the system under consideration, with the lens located at plane z=0, which implies ζ(0)=1, and therefore the scaled coordinates (X,Y)=(x,y). Without loss of generality we can take a TM type vMG beam as the field incident on the lens, and applying the Huygens-Fresnel diffraction integral we obtain a closed form expression for the field amplitude, this analysis has been presented in a previous work [24

], and it was shown that the propagating field is described by

Fig. 1. Focusing system under consideration. The lens is assumed to be thin and unapertured, f defines the lens focal distance, i.e. geometrical focus, z ₀ the actual focus, and Δz the focal shift.

E_{t} (\tilde{ξ}, \tilde{η}, z) = \frac{k}{i 2 F z} \exp [ik (z + \frac{r^{2}}{2 z})] \exp [- \frac{1}{4 F} (k_{t}^{2} + \frac{k^{2} r^{2}}{z^{2}})] {\tilde{\nabla}}_{t} W (\tilde{ξ}, \tilde{η}),

(13)

where the transverse gradient ∇̃ _t is taken over the complex elliptic variables (ξ̃,η̃). The factor F is defined as

F = \frac{1}{w_{0}^{2}} - \frac{ik}{2} (\frac{1}{z} - \frac{1}{f}),

(14)

and the complex elliptic variables (ξ̃,η̃) are determined by the relations

x = \frac{2 iFzh}{k} \cosh \tilde{ξ} \cos \tilde{η}, y = \frac{2 iFzh}{k} \sinh \tilde{ξ} \sin \tilde{η} .

(15)

The corresponding amplitude for the magnetic field is obtained by applying the relation H _t =(ε ₀/µ ₀)^1/2 ẑ×E _t . It should be mentioned that although the expressions for the electric and magnetic fields have been derived for a TM type vMG beam, the scalar factors to the left of the gradient operator in Eq. (13) remain unaltered for TE or any other state of polarization. We point out that the following analysis is done under the paraxial approximation, however it should be mentioned that for strongly focused beams, i.e. highly non-paraxial beams, the existence of large longitudinal components has been demonstrated [31

L. Novotny, M. R. Beversluis, K. S. Youngworth, and T. G. Brown, “Longitudinal Field Modes Probed by Single Molecules,” Phys. Rev. Lett. 86, 5251–5254 (2001). [CrossRef] [PubMed]

, 32

R. Dorn, S. Quabis, and G. Leuchs, “Sharper Focus for a Radially Polarized Light Beam,” Phys. Rev. Lett. 91, 233901 (2003). [CrossRef] [PubMed]

], and in that case they cannot be neglected.

2.3. Focal shift in vector Mathieu-Gauss beams

A common criteria employed to define the actual focus of a beam is to consider the point where the maximum on axis intensity occurs, and the difference between this point and the geometric focus corresponds to the focal shift. However, for vMG beams, and in general for beams with an on-axis null in their intensity profile, this criterion is not applicable. To deal with this situation we adopt two different approaches, the first one was also employed by Greene and Hall [16

P. L. Greene and D. G. Hall, “Focal shift in vector beams,” Opt. Express 4, 411–419 (1999). [CrossRef] [PubMed]

], and Lü and Peng [30

B. Lü and R. Peng, “Focal shift in Hermite-Gaussian beams based on the encircled-power criterion,” Opt. Laser Technol. 35, 435–440 (2003). [CrossRef]

], which defines the beam’s width as the radius r_w within which 80% of the beam’s power is encircled, defining the focal plane z=z ₀ as the plane where the width is minimum. The second approach is to compute the beam’s mean square radius r_m of the diffracted beam [9

S. De Nicola, D. Anderson, and M. Lisak, “Focal shift effects in diffracted focused beams,” Pure Appl. Opt. 7, 1249–1259 (1998). [CrossRef]

], and define the focal distance as the value z=z ₀ at which r_m achieves its minimum.

2.3.1. Encircled energy criterion

The power carried by the beam is calculated by means of the Poynting vector defined as 〈S〉=Re(E×H*)/2, which can be decomposed into transverse and longitudinal components as 〈S〉=〈s〉+〈s_z 〉ẑ. The time average of the longitudinal component 〈s_z 〉 gives the intensity distribution, whereas the total power carried by the beam along the propagation coordinate is calculated by integrating the transverse intensity over the whole plane (x,y), namely

P = \iint 〈 S 〉 \cdot \hat{z} dA = \iint 〈 S_{z} 〉 dA .

(16)

For a vMG beam, the s_z component of the Poynting vector is proportional to the squared amplitude of the transverse electric field, this is

〈 s_{z} 〉 = \frac{1}{2} {(\frac{ε_{0}}{μ_{0}})}^{\frac{1}{2}} (E_{t} \cdot E_{t}^{*}),

(17)

where E _t is defined in Eq. (13). We look for the smallest value of r_w satisfying the condition

\int_{0}^{2 π} \int_{0}^{r_{w}} 〈 s_{z} 〉 r d r d θ = 0.8 P,

(18)

the actual focus of the beam is the propagation distance z ₀ at which this value occurs, to find this value Eq. (18) has to be solved numerically. Once z ₀ is found, a straightforward calculation allows the determination of the focal shift Δz=f-z ₀ using this criterion.

2.3.2. Beam’s mean squared radius (MSR)

The computation of the beam’s mean squared radius, or which is equivalent its waist with respect to the radial coordinate, is accomplished through the calculation of the second order moment for the intensity distribution [9

S. De Nicola, D. Anderson, and M. Lisak, “Focal shift effects in diffracted focused beams,” Pure Appl. Opt. 7, 1249–1259 (1998). [CrossRef]

, 33

J. Alda, “Laser and Gaussian beam propagation and transformation,” in Encyclopedia of Optical Engineering , (Marcel Dekker, New York, 2003).

r_{m}^{2} (z) = \frac{1}{P} [\iint r^{2} 〈 s_{z} 〉 dA - m_{r}],

(19)

where r is the usual radial coordinate, r_m is the mean squared radius of the beam, and m_r =∬r〈s_z 〉dA its centroid, also known as the beam’s “center of mass”. The value of the centroid is always zero for a non-decentered distribution, which means that it is always located along the propagation axis. Similarly to the encircled energy criterion, the value z ₀ where r_m achieves its minimum will be defined as the focal distance, and the focal shift is thus calculated as the difference between the geometrical focus and this value, namely Δz=f-z ₀.

3. Numerical results

An important parameter of vMG beams is the ellipticity parameter q, when this parameter is equal to zero, vMG beams reduce to the special case of vector Bessel-Gauss beams studied in [16

P. L. Greene and D. G. Hall, “Focal shift in vector beams,” Opt. Express 4, 411–419 (1999). [CrossRef] [PubMed]

], however as q increases the elliptic features of the beam become more evident. The transverse intensity patterns and evolution of focusedvMG beams for different orders are shown in Fig. 2.

The top row images in Fig. 2 correspond to the initial profiles at plane z=0, whereas the bottom row consists of the initial snapshots from animations showing the evolution of the beam near the focal region. Note that we have used normalized units of length for all images. The videos show the three dimensional field structure in the focal region, a relevant feature is the on-axis null of all beams close to the geometrical focus. Unless stated otherwise, all following values for the simulation parameters remain unchanged for the rest of the paper, λ=632.8 nm, lens focal distance f=80 mm, waist of Gaussian envelope w ₀=1 mm, transverse wave number k_t =0.006 µm^-1, and ellipticity parameter q=5. Three different cases are presented, for instance all distributions in Fig. 2 were obtained using a TM polarized vMG beam, whose seed function is given by Eq. (7). Plots in Fig. 3(a) and (b) also correspond to a TM type solution, however the seed function in this case is known as an helical mode [29

J. C. Gutiérrez-Vega and M. A. Bandres, “Helmholtz-Gauss waves,” J. Opt. Soc. Am. A 22, 289–298 (2005). [CrossRef]

], and is given by a linear combination between Eqs. (7) and (8), namely W(ξ,η;q)=W^e _m (ξ,η;q)+iW^o _m (ξ,η;q), a distinctive feature of helical modes is that they carry orbital angular momentum (OAM) as they propagate [20

C. López-Mariscal, J. C. Gutiérrez-Vega, G. Milne, and K. Dholakia, “Orbital angular momentum transfer in helical Mathieu beams,” Opt. Express. 14, 4182–4187 (2006). [CrossRef] [PubMed]

], a behavior clearly observed in the corresponding videos of Fig. 3(a) and (b). Finally, in Fig. 3(c) and (d) we present distributions for a right circularly polarized beam, constructed from the superposition of TM and TE type solutions given in Eq. (3). It is a straightforward calculation to show that right (+) and left (-) circularly polarized solutions can be constructed as E^± _t=Z(ζ)G(r,ζ)[U ⁽¹⁾±i U ⁽²⁾].

Fig. 2. Transverse intensity profiles of a TM type focused vMG beam with seed function W=W^e _m (ξ,η;q), q=5 and beam orders (a)m=0, (b) m=1, (c) m=2, and (d) m=7. Plots in the top row correspond to initial profiles at plane z=0, bottom row images are snapshots for animations showing the three dimensional evolution in the focal region from z/f=0.95 to 1.05. (Movie files: (a) 1.68 MB, (b) 2.07 MB, (c) 1.94 MB, and (d) 2.31 MB). [Media 1][Media 2][Media 3][Media 4]

Fig. 3. Transverse intensity profiles of a TM type focused vMG beam with seed function W=W^e _m (ξ,η;q)+iW^o _m (ξ,η;q), q=5 and beam orders (a) m=4, and (b) m=7. Plots in (c) and (d) correspond to a right-circularly polarized vMG beam with seed function W=W^e _m (ξ,η;q), q=5, and beam orders m=0 and m=7 respectively. Plots in the top row show initial profiles at plane z=0, bottom row images are snapshots for animations showing the three dimensional evolution in the focal region from z/f=0.95 to 1.05. (Movie files: (a) 2.12 MB, (b) 2.33 MB, (c) 1.65 MB, and (d).2.25 MB) [Media 5][Media 6][Media 7][Media 8]

Propagation in the vicinity of the geometric focus goes from z/f=0.95 to 1.05. Apart from the change in the shape of the beam, and the appearance of an on-axis null, one can appreciate the existence of a plane at which the beam is narrowest, and a plane where it achieves its maximum intensity.

The existence of a focal shift in vMG beams was confirmed by means of the two approaches previously presented in Section 2.3. For instance, Fig. 4(a) shows the plots of the beam width according to the encircled energy criterion for different values of q, and mode orders m=0 through 7. It can be readily seen from these curves that the minimum width of the beam is achieved at a plane z<f, therefore a focal shift Δz towards the lens exists. For the special case when q=0, our findings are in agreement to those reported in [16

P. L. Greene and D. G. Hall, “Focal shift in vector beams,” Opt. Express 4, 411–419 (1999). [CrossRef] [PubMed]

], and we observe a decreasing value for Δz as the mode order is increased, however for larger values of q this behavior is lost, and some higher order modes may present larger shifts than lower order modes, but our results suggest that no specific dependence exist between the magnitude of the focal shift and q.

It is also observed that, in general, as the value of q increases the beam width decreases, i.e. the energy is confined to a smaller radius and the beam spreading is also reduced. Similar results were obtained using the second approach, which calculates the beam waist through the second order moment of the intensity with respect to the radial coordinate. These results are presented in Fig. 4(b), and again only for the special case q=0 the focal shift decreases with an increasing mode order. It is noteworthy that the magnitude of Δz calculated in Fig. 4(b) is larger than the one for Fig. 4(a), therefore depending on the specific feature of the beam in which we are interested, we may adopt either of the two criteria to calculate the focal shift.

The dependence of the shift on other beam parameters such as its polarization state and on the Gaussian-Fresnel number is of importance. Our simulations reveal the fact that the focal shift of a vMG beam is independent of its polarization state. This can be understood from the fact that the two criteria employed to define the actual focus of the beam are based on intensity. It is seen from Eqs. (1)–(4) that the intensities for a TM and TE vMG beams are equal in magnitude, thus the magnitude of a circularly polarizedvMG beam can be computed as |E ^± _t|²=|E ^TM _t |²+|E ^TE _t |², which indeed is equivalent to two times the magnitude of either polarization. The latter result is a consequence of the orthogonality between TM and TE polarized vMG beams.

The Gaussian Fresnel number of the beam is defined as

N_{w} = \frac{w_{0}^{2}}{λ f},

(20)

where f is the lens focal distance. Figure 5 shows the variation of the focal shift with respect to N_w , the simulation parameters are q=5, m=0, 1, and 7. The rest of the parameters are the same as in Fig. 2. A change in N_w was achieved by varying the waist of the Gaussian envelope w ₀ from 600 µm to 2000 µm. As expected, the focal shift of the beam for a given mode order approaches zero as N_w increases, in other words, for large values of w ₀ the Gaussian amplitude modulation tends to unity, it means that in free space we would have wider beams that do not spread significantly, however smaller values of w ₀ imply narrower beams, which spread rapidly on propagation. It should also be noted from these curves that higher order modes present smaller values of Δz, as we have mentioned in previous paragraphs, but the difference between the focal shifts for different orders is significantly reduced for large values of N_w .

Fig. 4. Beam waist plotted against the normalized propagation distance z/f, for q=[0,5,25] and mode orders m=0 through 7. (a) Calculation of beam waist using the encircled energy criteria and (b) using the mean squared radius of the intensity distribution. The minimum beam waist is achieved just before the geometric focus, i.e. z/f=1, indicating a focal shift Δz towards the lens.

Fig. 5. Relative focal shift versus Gaussian-Fresnel number N_w of the beam, for q=5 and beam orders m=[0,1,7]. A change in N_w was achieved by changing the Gaussian waist w ₀ from 600 µm to 2000 µm. All simulation parameters are the same as in Fig. 2. Notice how narrower beams present larger shifts.

4. Conclusions

In this paper we have studied the behavior of vMG beams when they are focused by an unapertured thin lens. The existence of a focal shift for vMG beams was confirmed by means of two different approaches to define the actual focus of a beam, the two are based on intensity measurements, one computes the encircled energy of the beam within a radius of interest, and the other is based on the calculation of the beam waist through the second order moment of the intensity distribution. The focal shift was found to be closer to the lens, slightly displaced from the geometric focus. The magnitude of this shift is highly dependent on the beam parameters and criterion used, for the encircled energy criterion typical values are a small percentage of the lens focal distance, in the range of 0.1% to 0.25%.

We also explored the dependence of this shift on other beam parameters such as its Gaussian Fresnel number N_w and polarization. Small values of N_w correspond to narrower beams which present larger focal shifts, in contrast to wider beams with large N_w that possess smaller focal shifts. Our findings also revealed that for vMG beams the focal shift is independent of its polarization state. The latter is a direct consequence of the orthogonality between the different solutions of the vector paraxial wave equation given by Eqs. (1) and (3).

Acknowledgments

This research was supported by Consejo Nacional de Ciencia y Tecnología (grant 42808) and by Tecnológico de Monterrey (grant CAT-007).

References and links

1.	E. Collet and E. Wolf, “Symmetry properties of focused fields,” Opt. Lett. 5, 264–266 (1980). [CrossRef]
2.	J. J. Stamnes and B. Spjelkavik, “Focusing at small angular apertures in the Debye and Kirchhoff approximations,” Opt. Commun. 40, 81–85 (1981). [CrossRef]
3.	Y. Li and E. Wolf, “Focal shifts in diffracted converging spherical waves,” Opt. Commun. 39, 211–215 (1981). [CrossRef]
4.	Y. Li, “Dependence of the focal shift on Fresnel number and f number,” J. Opt. Soc. Am. 72, 770–774 (1982). [CrossRef]
5.	W. H. Carter, “Focal shift and concept of effective Fresnel number of a Gaussian laser beam,” Appl. Opt. 21, 1989–1994 (1982). [CrossRef] [PubMed]
6.	V. N. Mahajan, “Axial irradiance and optimum focusing of laser beams,” Appl. Opt. 19, 3042–3053 (1983). [CrossRef]
7.	M. Martínez-Corral, C. J. Zapata Rodríguez, P. Andrés, and E. Silvestre, “Effective Fresnel-Number concept for evaluating the relative focal shift in focused beams,” J. Opt. Soc. Am. A. 15, 449–455 (1998). [CrossRef]
8.	C. J. R. Sheppard and P. Török, “Dependence of focal shift on Fresnel number and angular aperture,” Opt. Lett. 23, 1803–1804 (1998). [CrossRef]
9.	S. De Nicola, D. Anderson, and M. Lisak, “Focal shift effects in diffracted focused beams,” Pure Appl. Opt. 7, 1249–1259 (1998). [CrossRef]
10.	T. C. Poon, “Focal shift in focused annular beams,” Opt. Commun. 65, 401–406 (1988). [CrossRef]
11.	W. H. Carter and M. F. Aburdene, “Focal shift in Laguerre-Gaussian beams,” J. Opt. Soc. Am. A 10, 1949–1952 (1987). [CrossRef]
12.	B. Lü and W. Huang, “Focal shift in unapertured Bessel-Gauss beams,” Opt. Commun. 109, 43–46 (1994). [CrossRef]
13.	S. De Nicola, “On axis focal shift effects in focused truncated J ₀ Bessel beams,” Pure Appl. Opt. 5, 827–831 (1996). [CrossRef]
14.	Z. Hricha, L. Dalil-Essakali, and A. Belafhal, “Axial intensity distribution and focal shifts of focused partially coherent conical Bessel-Gauss beams,” Opt. Quantum. Electron. 35, 101–110 (2003). [CrossRef]
15.	M. Dong and J. Pu, “Effective Fresnel number and the focal shifts of focused partially coherent beams,” J. Opt. Soc. Am. A 24, 192–196 (2007). [CrossRef]
16.	P. L. Greene and D. G. Hall, “Focal shift in vector beams,” Opt. Express 4, 411–419 (1999). [CrossRef] [PubMed]
17.	J. Pu and B. Lü, “Focal shifts in focused nonuniformly polarized beams,” J. Opt. Soc. Am. A 18, 2760–2766 (2001). [CrossRef]
18.	J. C. Gutiérrez-Vega, M. D. Iturbe-Castillo, and S. Chávez-Cerda, “Alternative formulation for invariant optical fields,” Opt. Lett. 25, 1493–1495 (2000). [CrossRef]
19.	Y. V. Kartashov, A. A. Egorov, V. A. Vysloukh, and L. Torner, “Shaping soliton properties in Mathieu lattices,” Opt. Lett. 31, 238–240 (2006). [CrossRef] [PubMed]
20.	C. López-Mariscal, J. C. Gutiérrez-Vega, G. Milne, and K. Dholakia, “Orbital angular momentum transfer in helical Mathieu beams,” Opt. Express. 14, 4182–4187 (2006). [CrossRef] [PubMed]
21.	C. A. Dartora and H. E. Hernández-Figueroa, “Properties of a localized Mathieu pulse,” J. Opt. Soc. Am. A 21, 662–667 (2004). [CrossRef]
22.	J. Davila-Rodriguez and J. C. Gutiérrez-Vega, “Helical Mathieu and parabolic localized pulses,” J. Opt. Soc. Am. A 24, 3449–3455 (2007). [CrossRef]
23.	M. A. Bandres and J. C. Gutiérrez-Vega, “Vector Helmholtz-Gauss and vector Laplace-Gauss beams,” Opt. Lett. 30, 2155–2157 (2005). [CrossRef] [PubMed]
24.	R. I. Hernández-Aranda, M. A. Bandres, and J. C. Gutiérrez Vega, “Propagation dynamics of vector Mathieu-Gauss beams,” in Laser Beam Shaping VII, Fred M. Dickey and David L. Shealy, eds., Proc. SPIE 6290, 629011 (2006). [CrossRef]
25.	R. I. Hernandez-Aranda, J. C. Gutiérrez-Vega, M. Guizar-Sicairos, and M. A. Bandres, “Propagation of generalized vector Helmholtz-Gauss beams through paraxial optical systems,” Opt. Express. 14, 8974–8988 (2006). [CrossRef] [PubMed]
26.	A. Chafiq, Z. Hricha, and A. Belafhal, “Propagation properties of vector Mathieu-Gauss beams,” Opt. Commun. 275, 165–169 (2007). [CrossRef]
27.	J. A. Stratton, Electromagnetic Theory (McGraw-Hill, New York, 1941).
28.	M. Lax, W. H. Louisell, and W. B. McKnight, “From Maxwell to paraxial wave optics,” Phys. Rev. A 11, 1365–1370 (1975). [CrossRef]
29.	J. C. Gutiérrez-Vega and M. A. Bandres, “Helmholtz-Gauss waves,” J. Opt. Soc. Am. A 22, 289–298 (2005). [CrossRef]
30.	B. Lü and R. Peng, “Focal shift in Hermite-Gaussian beams based on the encircled-power criterion,” Opt. Laser Technol. 35, 435–440 (2003). [CrossRef]
31.	L. Novotny, M. R. Beversluis, K. S. Youngworth, and T. G. Brown, “Longitudinal Field Modes Probed by Single Molecules,” Phys. Rev. Lett. 86, 5251–5254 (2001). [CrossRef] [PubMed]
32.	R. Dorn, S. Quabis, and G. Leuchs, “Sharper Focus for a Radially Polarized Light Beam,” Phys. Rev. Lett. 91, 233901 (2003). [CrossRef] [PubMed]
33.	J. Alda, “Laser and Gaussian beam propagation and transformation,” in Encyclopedia of Optical Engineering , (Marcel Dekker, New York, 2003).

OCIS Codes
(050.1940) Diffraction and gratings : Diffraction
(140.3300) Lasers and laser optics : Laser beam shaping
(220.2560) Optical design and fabrication : Propagating methods
(260.1960) Physical optics : Diffraction theory
(350.5500) Other areas of optics : Propagation

ToC Category:
Physical Optics

History
Original Manuscript: January 18, 2008
Revised Manuscript: April 4, 2008
Manuscript Accepted: April 8, 2008
Published: April 11, 2008

Citation
Raul I. Hernandez-Aranda and Julio C. Gutiérrez-Vega, "Focal shift in vector Mathieu-Gauss beams," Opt. Express 16, 5838-5848 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-8-5838

Sort: Year | Journal | Reset

References

E. Collet and E. Wolf, "Symmetry properties of focused fields," Opt. Lett. 5, 264-266 (1980). [CrossRef]
J. J. Stamnes and B. Spjelkavik, "Focusing at small angular apertures in the Debye and Kirchhoff approximations," Opt. Commun. 40, 81-85 (1981). [CrossRef]
Y. Li and E. Wolf, "Focal shifts in diffracted converging spherical waves," Opt. Commun. 39, 211-215 (1981). [CrossRef]
Y. Li, "Dependence of the focal shift on Fresnel number and f number," J. Opt. Soc. Am. 72, 770-774 (1982). [CrossRef]
W. H. Carter, "Focal shift and concept of effective Fresnel number of a Gaussian laser beam," Appl. Opt. 21, 1989-1994 (1982). [CrossRef] [PubMed]
V. N. Mahajan, "Axial irradiance and optimum focusing of laser beams," Appl. Opt. 19, 3042-3053 (1983). [CrossRef]
M. Martinez-Corral, C. J. Zapata Rodriguez, P. Andres, and E. Silvestre, "Effective Fresnel-Number concept for evaluating the relative focal shift in focused beams," J. Opt. Soc. Am. A. 15, 449-455 (1998). [CrossRef]
C. J. R. Sheppard and P. Torok, "Dependence of focal shift on Fresnel number and angular aperture," Opt. Lett. 23, 1803-1804 (1998). [CrossRef]
S. De Nicola, D. Anderson, and M. Lisak, "Focal shift effects in diffracted focused beams," Pure Appl. Opt. 7, 1249-1259 (1998). [CrossRef]
T. C. Poon, "Focal shift in focused annular beams," Opt. Commun. 65, 401-406 (1988). [CrossRef]
W. H. Carter and M. F. Aburdene, "Focal shift in Laguerre-Gaussian beams," J. Opt. Soc. Am. A 10, 1949-1952 (1987). [CrossRef]
B. Lu and W. Huang, "Focal shift in unapertured Bessel-Gauss beams," Opt. Commun. 109, 43-46 (1994). [CrossRef]
S. De Nicola, "On axis focal shift effects in focused truncated J0 Bessel beams," Pure Appl. Opt. 5, 827-831 (1996). [CrossRef]
Z. Hricha, L. Dalil-Essakali, and A. Belafhal, "Axial intensity distribution and focal shifts of focused partially coherent conical Bessel-Gauss beams," Opt. Quantum. Electron. 35, 101-110 (2003). [CrossRef]
M. Dong and J. Pu, "Effective Fresnel number and the focal shifts of focused partially coherent beams," J. Opt. Soc. Am. A 24, 192-196 (2007). [CrossRef]
P. L. Greene and D. G. Hall, "Focal shift in vector beams," Opt. Express 4, 411-419 (1999). [CrossRef] [PubMed]
J. Pu and B. Lu, "Focal shifts in focused nonuniformly polarized beams," J. Opt. Soc. Am. A 18, 2760-2766 (2001). [CrossRef]
J. C. Gutierrez-Vega, M. D. Iturbe-Castillo, and S. Chavez-Cerda, "Alternative formulation for invariant optical fields," Opt. Lett. 25, 1493-1495 (2000). [CrossRef]
Y. V. Kartashov, A. A. Egorov, V. A. Vysloukh, and L. Torner, "Shaping soliton properties in Mathieu lattices," Opt. Lett. 31, 238-240 (2006). [CrossRef] [PubMed]
C. Lopez-Mariscal, J. C. Gutierrez-Vega, G. Milne, and K. Dholakia, "Orbital angular momentum transfer in helical Mathieu beams," Opt. Express. 14, 4182-4187 (2006). [CrossRef] [PubMed]
C. A. Dartora and H. E. Hernandez-Figueroa, "Properties of a localized Mathieu pulse," J. Opt. Soc. Am. A 21, 662-667 (2004). [CrossRef]
J. Davila-Rodriguez and J. C. Gutierrez-Vega, "Helical Mathieu and parabolic localized pulses," J. Opt. Soc. Am. A 24, 3449-3455 (2007). [CrossRef]
M. A. Bandres and J. C. Guti’errez-Vega, "Vector Helmholtz-Gauss and vector Laplace-Gauss beams," Opt. Lett. 30, 2155-2157 (2005). [CrossRef] [PubMed]
R. I. Hernandez-Aranda, M. A. Bandres, and J. C. Gutierrez Vega, "Propagation dynamics of vector Mathieu- Gauss beams," in Laser Beam Shaping VII, Fred M. Dickey and David L. Shealy, eds., Proc. SPIE 6290, 629011 (2006). [CrossRef]
R. I. Hernandez-Aranda, J. C. Gutierrez-Vega, M. Guizar-Sicairos, and M. A. Bandres, "Propagation of generalized vector Helmholtz-Gauss beams through paraxial optical systems," Opt. Express. 14, 8974-8988 (2006). [CrossRef] [PubMed]
A. Chafiq, Z. Hricha, and A. Belafhal, "Propagation properties of vector Mathieu-Gauss beams," Opt. Commun. 275, 165-169 (2007). [CrossRef]
J. A. Stratton, Electromagnetic Theory (McGraw-Hill, New York, 1941).
M. Lax, W. H. Louisell, and W. B. McKnight, "From Maxwell to paraxial wave optics," Phys. Rev. A 11, 1365- 1370 (1975). [CrossRef]
J. C. Gutierrez-Vega and M. A. Bandres, "Helmholtz-Gauss waves," J. Opt. Soc. Am. A 22, 289-298 (2005). [CrossRef]
B. Lu and R. Peng, "Focal shift in Hermite-Gaussian beams based on the encircled-power criterion," Opt. Laser Technol. 35, 435-440 (2003). [CrossRef]
L. Novotny, M. R. Beversluis, K. S. Youngworth, and T. G. Brown, "Longitudinal Field Modes Probed by Single Molecules," Phys. Rev. Lett. 86, 5251-5254 (2001). [CrossRef] [PubMed]
R. Dorn, S. Quabis, and G. Leuchs, "Sharper Focus for a Radially Polarized Light Beam," Phys. Rev. Lett. 91, 233901 (2003). [CrossRef] [PubMed]
J. Alda, "Laser and Gaussian beam propagation and transformation," in Encyclopedia of Optical Engineering, (Marcel Dekker, New York, 2003).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper

Figures


Fig. 1.	Fig. 2.	Fig. 3.


Fig. 4.	Fig. 5.

Multimedia

Multimedia Files	Recommended Software
» Media 1: AVI (1725 KB)	QuickTime
» Media 2: AVI (2125 KB)	QuickTime
» Media 3: AVI (1992 KB)	QuickTime
» Media 4: AVI (2370 KB)	QuickTime
» Media 5: AVI (2175 KB)	QuickTime
» Media 6: AVI (2388 KB)	QuickTime
» Media 7: AVI (1695 KB)	QuickTime
» Media 8: AVI (2312 KB)	QuickTime

« Previous Article | Next Article »

OSA is a member of CrossRef.

Select a Conference
Session or Paper #
Year

OpticsInfoBase

Optics Express

Focal shift in vector Mathieu-Gauss beams

Article Tools

Article Navigation

Abstract

1. Introduction

2. Vector Mathieu-Gauss beams

2.1. General description of a vector Mathieu-Gauss beam

2.2. Focusing of vector Mathieu-Gauss beams

2.3. Focal shift in vector Mathieu-Gauss beams

2.3.1. Encircled energy criterion

2.3.2. Beam’s mean squared radius (MSR)

3. Numerical results

4. Conclusions

Acknowledgments

References and links

References

Appl. Opt.

J. Opt. Soc. Am.

J. Opt. Soc. Am. A

J. Opt. Soc. Am. A.

Opt. Commun.

Opt. Express

Opt. Express.

Opt. Laser Technol.

Opt. Lett.

Opt. Quantum. Electron.

Phys. Rev. A

Phys. Rev. Lett.

Pure Appl. Opt.

Other

Cited By

Figures

Multimedia