Focal shift in vector beams

We examine the three-dimensional intensity distribution of vector Bessel–Gauss beams with general polarization near the focus of a nonaperturing thin lens. Recently reported linearly and azimuthally polarized Bessel–Gauss beams are members of this family. We define the width and focal plane of such a vector beam using an encircled-energy criterion and calculate numerically that, as for scalar beams, the true focus occurs not at the geometric focus of the lens but rather somewhat closer to the lens. This focal shift depends on the mode number and system parameters and is largest for a narrow beam, long lens focal length, and large wavelength. © 1999 Optical Society of America OCIS codes: (050.1940) Diffraction; (220.2560) Focus References and links 1. W. H. Carter, “Focal shift and concept of effective Fresnel number for a Gaussian laser beam,” Appl. Opt. 21, 1989-1994 (1982). 2. Y. Li, “Dependence of the focal shift on Fresnel number and f number,” J. Opt. Soc. Am. 72, 770-775 (1982). 3. Y. Li and E. Wolf, “Three-dimensional intensity distribution near the focus in systems of different Fresnel numbers,” J. Opt. Soc. Am. A 1, 801-808 (1984). 4. 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). 5. A. E. Siegman, Lasers, (University Science Books, Mill Valley, Calif., 1986), 769-773. 6. Y. Li and E. Wolf, “Focal shifts in diffracted converging spherical waves,” Opt. Commun. 39, 211-215 (1981). 7. M. P. Givens, “Focal shifts in diffracted converging spherical waves,” Opt. Commun. 41, 145-148 (1982). 8. T.-C. Poon, “Focal shift in focused annular beams,” Opt. Commun. 65, 401-406 (1988). 9. M. Martínez-Corral, C. J. Zapata-Rodríguez, P. Andrés and M. Kowalczyk, “Analytical formula for calculating the focal shift in apodized systems,” J. Mod. Opt. 45, 1671-1679 (1998). 10. Y. Li and F. T. S. Yu, “Intensity distribution near the focus of an apertured focused Gaussian beam,” Opt. Commun. 70, 1-7 (1989). 11. Y. Li, “Focal shift formula for focused, apertured Gaussian beams,” J. Mod. Opt. 39, 1761-1764 (1992). 12. Y. Li, “Oscillations and discontinuity in the focal shift of Gaussian laser beams,” J. Opt. Soc. Am. A 3, 17611765 (1986). 13. Y. Li, “Axial intensity and the focal shift of a focused non-truncated elliptical Gaussian beam,” Optik 80, 121123 (1988). 14. W. H. Carter and M. F. Aburdene, “Focal shift in Laguerre-Gaussian beams,” J. Opt. Soc. Am. A 4, 1949-1952 (1987). 15. S. De Nicola, “On-axis focal shift effects in focused truncated J0 Bessel beams,” Pure Appl. Opt. 5, 827-831 (1996). 16. B. Lü and W. Huang, “Focal shift in unapertured Bessel-Gauss beams,” Opt. Commun. 109, 43-46 (1994). 17. B. Lü and W. Huang, “Three-dimensional intensity distribution of focused Bessel-Gauss beams,” J. Mod. Opt. 43, 509-515 (1996). 18. R. Borghi, M. Santarsiero and S. Vicalvi, “Focal shift of focused flat-topped beams,” Opt. Commun. 154, 243248 (1998). 19. F. Gori, G. Guattari and C. Padovani, “Bessel-Gauss beams,” Opt. Commun. 64, 491-495 (1987). (C) 1999 OSA 10 May 1999 / Vol. 4, No. 10 / OPTICS EXPRESS 411 #9094 $15.00 US Received March 08, 1999; Revised May 07, 1999 20. R. H. Jordan and D. G. Hall, “Free-space azimuthal paraxial wave equation: the azimuthal Bessel-Gauss beam solution,” Opt. Lett. 19, 427-429 (1994). 21. D. G. Hall, “Vector-beam solutions of Maxwell's wave equation,” Opt. Lett. 21, 9-11 (1996). 22. P. L. Greene and D. G. Hall, “Properties and diffraction of vector Bessel-Gauss beams,” J. Opt. Soc. Am. A 15, 3020-3027 (1998). 23. L. W. Casperson, D. G. Hall and A. A. Tovar, “Sinusoidal-Gaussian beams in complex optical systems,” J. Opt. Soc. Am. A 14, 3341-3348 (1997). 24. See, for example, J. W. Goodman, Introduction to Fourier Optics (McGraw-Hill, New York, 1968), pp. 59-60.


Introduction
The phenomenon of focal shift, wherein the point of maximum intensity of a diffracted field is located not at the geometric focus, but rather somewhat closer to the diffraction plane, is well known.In general, the magnitude of this shift depends strongly on the Fresnel number of the aperture, or on the effective Fresnel number of the field in unapertured cases [1][2][3][4].The shift is significant only when the Fresnel number of the system is sufficiently small; for most imaging applications, therefore, it may safely be neglected.However, in laser cavities and other laser systems the Fresnel number is often on the order of unity or smaller [5], in which case the focal shift becomes important.A focal shift has been observed for a variety of systems, including a converging spherical wave diffracted by a circular aperture [2,3,6,7], a uniformly illuminated aperture with a central obscuration [8], and more general diffracting screens [4,9].It has likewise been noted and thoroughly studied for Gaussian beams focused by an aperturing [10,11] or nonaperturing [1,12,13] lens with a low Fresnel number, as well as for other types of focused beams, including Laguerre-Gauss [14], Bessel [15], J 0 -Bessel-Gauss [16,17], and flattened Gaussian beams [18].However, all these diverse systems involve only fields that are linearly polarized and satisfy a scalar wave equation.
In the present work we examine the focal shift of a new family of beams, found by solving the vector wave equation in the paraxial limit.These vector beams have general polarization; recently reported linearly [19] and azimuthally [20] polarized Bessel-Gauss beams may be identified as members of the set [21].An attractive property of these more general vector Bessel-Gauss beam solutions is their correspondence to modes of important cylindrical waveguides and resonators, in particular the modes emitted by a concentric-circlegrating surface-emitting (CCGSE) semiconductor laser [22].A correspondence can also be found with the modes supported within an optical fiber, in the limit of large Gaussian beam size [23], and with the free-space modes excited by light exiting the end of such a fiber [21].
We begin with a summary of the set of general vector beam solutions to be considered, including a description of their form after passing through a nonaperturing thin lens.We then define the waist and focal plane of such beams and present images and animated figures exploring their qualitative features near the focus of a long-focal-length lens, as well as numerical solutions showing the dependence of the focal shift on the mode number and various beam parameters.We find that the focal shift depends strongly on the beam and system parameters; it is largest for a narrow beam, long focal length, and large wavelength.

Calculations
The system under consideration is shown schematically in Fig. 1.The beam is a member of the set of paraxial vector Bessel-Gauss beams, which are most easily expressed in cylindrical coordinates.We can write its electric field in terms of a vector amplitude U(ρ,φ,z): The exp(-iωt) harmonic time dependence will be taken as implicit throughout the following work.The field amplitude is assumed to include a Gaussian envelope g(ρ,z) to provide the necessary transverse localization in free space,  Here g(ρ,z) is the familiar Gaussian solution, where w 0 is the Gaussian beam waist, L kw . Considering only transversely polarized fields for concreteness, we seek a vector function f ρ φ , , z ( ) with only ρ and φ components, which may be described by [21,22] Here the notation {ρ, φ} denotes either the radial or the azimuthal component, m is the azimuthal mode number, β is an arbitrary constant, and the constant a m is the mode amplitude.The subscript m has been suppressed from f f , , and Θ Θ here for economy of notation.The radial component f ρ uses the upper operator (+) in Eq. ( 4) and either of the upper two trigonometric functions in Eq. ( 7), and the azimuthal component f φ uses the lower operator (-) in Eq. ( 4) and the corresponding lower function in Eq. (7).Note that the sinusoidal functions in Eq. ( 7) provide a characteristic 2m-fold azimuthal symmetry for most parameter sets [22].The familiar scalar Fresnel diffraction integral [24] can be converted to cylindrical coordinates and used to obtain a diffraction integral for each component of a general transverse vector beam affected by a circularly symmetric pupil (e.g., a lens or circular aperture).For the case of a nonaperturing thin lens, with a pupil function modeled by a quadratic phase object, the vector field amplitude following the lens is found to be where I ν is the modified Bessel function of order ν, f is the geometric focal length of the lens, and the notation is the same as that used in Eq. ( 4).We further define a Fresnel number N w for the Gaussian portion of the beam, as well as a separate Fresnel number N β for the Bessel portion: Because we are considering higher-order Bessel functions, this definition for N β differs by a factor of (m+2) 2 from that used in the discussion of linearly polarized J 0 -Bessel-Gauss beams [16,17]; note that the spot radius of J m (βρ) is approximately (m + 2)/β.It is worthwhile to note that the focal shift seen in scalar beams is embodied in this more general vector analysis.The familiar Gaussian beam, as well as other linearly polarized beams, may be found as members of the family of vector Bessel-Gauss beams [21], and thus the origins of their focal shift must be the same as for the more general beams.In particular, the focal shift appears when the distance from the wavefront at the lens to the point of observation is not approximated as the focal length of the lens [7].
In usual treatments of focal shift in linearly polarized, scalar beams, the location of the actual focus is defined to be the point at which the beam reaches its maximum on-axis intensity.However, this definition is not applicable to all types of beams.In the case of a general vector beam, the presence of the higher-order Bessel functions in Eq. ( 4) results in an on-axis null for all mode numbers except m = 1, since all Bessel functions J ν (u) except ν = 0 vanish as u goes to zero.We therefore use a more general definition for the focal plane of these beams.For a standard linearly polarized Gaussian beam, slightly more than 80% of the energy propagates within a radius w 0 of the axis, so we define the width ρ 0 of a vector beam as that radius within which 80% of the beam's power is enclosed, and the focal plane z = z 0 as the plane in which that width is at its smallest value.Note that the plane containing the 2m points of maximum intensity of a vector beam is not this focal plane, but rather is in general between the focal plane defined here and the geometric focus.
The intensity of the beam is the time average of the z component of the Poynting vector, The magnetic field H is easily found using Maxwell's equations, and the asterisk (*) denotes (C) 1999 OSA 10 May 1999 / Vol. 4, No. 10 / OPTICS EXPRESS 414 the complex conjugate.The power carried may be found by integrating this intensity over the ρ-φ plane, so we seek the propagation distance z 0 at which the smallest value of ρ 0 satisfies , , , , . .The width of a vector beam at a given propagation distance z must be calculated from Eq. ( 13) numerically, and therefore the focal shift ∆z = f -z 0 must also be found numerically.

Numerical results
The three-dimensional structure of these vector beams in the region of the geometric focus is complicated, as shown in Fig. 2. In the animation associated with each of these snapshots the viewer is traveling along the direction of propagation of a vector beam, from z = 0.9 f to z = 1.1 f.Animations (a)-(d) show modes m = 0, 1, 2, and 7, respectively, for λ = 632.8nm, f = 80 mm, w 0 = 1 mm, and β = 0 00004 .µm -1 ≈ 0.006 µm -1 .These values correspond to N w ≈ 20 and N β ≈ 0 5 2 2 .( ) m + .As described in previous work, this relatively large value of the inverse Bessel width β ensures that the beams show the characteristic 2m-fold azimuthal symmetry and that the m = 0 and m = 2 beams differ [22].The white lines in these snapshots have been added to aid in the interpretation of Fig. 3.The most noteworthy features are the change in the shape of the beam and the contrast reversal that occur in the vicinity of the geometric focus; the distinction between the plane at which each beam is narrowest and the plane containing its highest intensity is also seen in these sequences.The structure of the beams is shown further in Fig. 3, which presents the intensity of each mode in a longitudinal cross-section containing the beam axis.Images (a)-(d) show m = 0, 1, 2, and 7, respectively.A specific choice of azimuthal coordinate is necessary to fully describe the ρ-z plane; in each case we have selected φ = π/(2m), indicated by a white line added to the snapshots of Fig. 2. The beam parameters are the same as those in Fig. 2. In images (b) and (c), note that the beam does not vanish entirely in the region surrounding z = f, but because of the contrast reversal noted above the energy is no longer found at the particular azimuth φ corresponding to this plane.The exact focal shift is not evident from Fig. 3, but a distinct asymmetry in the intensity of the field on either side of the geometric focus may be seen, indicating that the energy is somewhat more concentrated just before z = f, and therefore that the beam is likely to be narrower there.
This focal shift is confirmed in Fig. 4, a plot of the beam width in micrometers, calculated from Eq. ( 13), against the propagation distance, for modes m = 0 through 7, with the same beam parameters as used in Fig. 2. The asymmetry about z = f shown in Fig. 3  also apparent here.Each of the modes reaches a minimum waist just before the geometric focal length, showing a small focal shift ∆z toward the lens that decreases (moves toward the geometric focus) slowly with mode number for m ≥ 1.The shift depends more strongly on the other beam parameters, however, as demonstrated in the following two figures.Fig. 5 shows its dependence on the Gaussian Fresnel number N w for the modes m = 0 and m = 1.Both beams have N β ≈ 2, corresponding to β ≈ 0.006 µm -1 and β ≈ 0.0095 µm -1 , respectively, and to change N w the Gaussian width w 0 is varied from 300 µm to 1900 µm.Similarly, Fig. 6 plots the focal shift against the Bessel Fresnel number N β for these two modes, with N w ≈ 20 and the change in N β accomplished by varying β from 0.001 µm -1 (0.0015 µm -1 ) to 0.010 µm -1 for m = 0 (m = 1).As these two figures show, the focal shift decreases as either N w or N β increases, and is slightly larger for m = 1 than for m = 0, an effect which becomes more evident for smaller N w or N β .Note that differences between the two modes in Fig. 6 arise (C) 1999 OSA 10 May 1999 / Vol. 4, No. 10 / OPTICS EXPRESS 417 primarily from the ( ) m + 2 2 factor in N β , as defined in Eq. ( 11); when the focal shift is plotted against β itself, the lines are nearly identical.

Fig. 1 .
Fig. 1.Schematic diagram of the geometry under consideration, showing the notation used.The blue dashed lines represent the geometric focus of the lens, and the red curves the profile of the focused beam.

Fig. 2 .
Fig. 2. (1 MB each) Snapshots from animations showing vector beam intensity patterns as the viewer travels along the direction of propagation, near the geometric focus.Images (a)-(d) show modes m = 0, 1, 2, and 7, respectively.The beam parameters f = 80 mm and β ≈ 0.006 µm -1 correspond to Fresnel numbers Nw ≈ 20 and N β ≈ 0.5 (m + 2) 2 .Each mode has been independently normalized, but the individual images within each animation are not separately scaled.The white lines have been added to aid in the interpretation of Fig. 3.

7 Fig. 4 .
Fig. 4. Beam width in micrometers plotted against normalized propagation distance z/f for several mode numbers m.The beam parameters are as in Fig. 2. The minimum width for each mode occurs slightly before z/f = 1, indicating a positive focal shift ∆z that decreases slowly with mode number for m ≥ 1.

1 Fig. 6 .
Fig.6.Relative focal shift ∆z/f plotted against the Bessel Fresnel number N β , varied by changing β, for modes m = 0 (blue line) and 1 (red line).The Gaussian Fresnel number is Nw ≈ 20.Differences between the two modes arise primarily from the (m + 2) 2 factor in N β , as defined in Eq. (11); when the focal shift is plotted against β itself, the lines are nearly identical.The interaction between the Bessel and Gaussian Fresnel numbers is shown in Fig.7, which plots the relative focal shift against N w for m = 0 and several values of N β .The largest