## Abstract

We derive solutions for radially polarized Bessel-Gauss beams in free-space by superimposing decentered Gaussian beams with differing polarization states. We numerically show that the analytical result is applicable even for large semi-aperture angles, and we experimentally confirm the analytical expression by employing a fiber-based mode-converter.

© 2013 Optical Society of America

## 1. Introduction

Since their introduction several decades ago, Bessel-Gauss beams have sparked tremendous interest. As an approximation to the non-diffracting Bessel beam, Bessel-Gauss beams offer an extended depth of focus as well as other intriguing features such as an annular intensity distribution in the far-field [1–3]. Such features have fueled wide-ranging interest in potential applications ranging from materials processing [4] and micro-manipulation [5] to strong-field physics [6].

Interest has also grown in Bessel-Gauss beams with vector polarization states, in particular cylindrical polarization, i.e. radial or azimuthal polarization. Such cylindrical Bessel-Gauss beams have been explored in the context of laser electron acceleration [7] as well as in studies concerning tight-focusing of optical beams for microscopy and lithography applications [8–10]. Additionally, emission patterns from surface-emitting semiconductor lasers have been described with radially and azimuthally polarized Bessel-Gauss beam expressions formulated via paraxial vector wave equations [11–13].

In this contribution, we derive expressions for radially polarized Bessel-Gauss beams through a novel scalar-wave approach. The approach involves superposing multiple decentered scalar Gaussian beams with differing polarization states. The expressions we derive for radially polarized Bessel-Gauss beams provide new insight into the behavior of such beams in first-order, ABCD, optical systems. Moreover, they allow for straightforward numerical calculations. Additionally, we confirm our radially polarized Bessel-Gauss beam expressions through an experiment based on a fiber mode-converter.

The outline of the paper is as follows: In section 2, we define linearly polarized decentered Gaussian beams, and we superpose these component beams to formulate expressions for cylindrical, i.e. radially and azimuthally, polarized Bessel-Gauss beams in a way hitherto not demonstrated in the literature, to the best or our knowledge. In section 3, we derive the electromagnetic fields of the cylindrical Bessel-Gauss beams from the vector potential description developed in the preceding section, and we numerically confirm our analytical expressions through non-paraxial numerical simulations. In section 4, we describe an experimental verification of our radially polarized Bessel-Gauss expression. In section 5, we conclude this work.

## 2. Cylindrical Bessel-Gauss beams as superpositions of decentered Gaussian beams

In the following section, we first define the decentered Gaussian beam and propagate this beam in free-space using an adapted paraxial approximation. Next, we superpose many such decentered Gaussian component beams to generate cylindrical, i.e. radially and azimuthally polarized, Bessel-Gauss beams.

#### 2.1. Decentered Gaussian beams

Throughout this section and the following ones, we will work with the vector potential **A** in the Lorenz gauge. With this gauge condition, the vector potential satisfies the wave equation, and we can write **A** in a familiar form as a superposition of plane-wave solutions:

**n̂**represents a cartesian vector (or any linear combination thereof),

*â*(

*k*,

_{x}*k*) denotes the distribution of plane-wave amplitudes, and the familiar condition ${\left(\omega /c\right)}^{2}={k}^{2}={k}_{x}^{2}+{k}_{y}^{2}+{k}_{z}^{2}$ is satisfied. In the following, we consider, but are not limited to, monochromatic light with angular frequency

_{y}*ω*

_{0}.

Let us consider a Gaussian distribution of plane-wave amplitudes in the vicinity of a principal wavevector of fixed orientation given by **K** = (*K _{x}*,

*K*,

_{y}*K*). We denote deviations from this principal orientation in the

_{z}*x*and

*y*direction as

*α*and

*β*, respectively, such that:

*a*(

*α*,

*β*) ≡

*â*(

*K*+

_{x}*α*,

*K*+

_{y}*β*) is given by:

*w*

_{0}describes the width of the Gaussian distribution (and will describe the waist of the resulting beam), and

*x*

_{d}_{0}and

*y*

_{d}_{0}are decenter parameters whose importance will soon be elucidated. Inserting Eq. (4) into Eq. (1), setting

*z*= 0, and integrating, we find an expression for the decentered Gaussian beam at its focus (

*z*= 0) [14, 15]:

**A**

*is the vector potential of the decentered Gaussian beam,*

_{dG}*ε*=

_{x}*K*and

_{x}/K_{z}*ε*=

_{y}*K*. Additionally,

_{y}/K_{z}*q*

_{0}= −

*iz*, where ${z}_{R}={K}_{z}{w}_{0}^{2}/2$. Inspecting Eq. (5), we see that the expression for

_{R}**A**

*(*

_{dG}*x*,

*y*,

*z*= 0,

*t*) resembles that of a conventional Gaussian beam at its focus with a decentered focal spot and a tilted propagation direction. The decenterd focal spot is described by the parameters

*x*

_{d}_{0}and

*y*

_{d}_{0}. They represent the coordinates of the beam focus in the

*x*–

*y*plane, which may also be written in terms of the vector

**r**

_{d}_{0}= (

*x*

_{d}_{0},

*y*

_{d}_{0}). The tilted propagation direction is characterized by the parameters

*ε*and

_{x}*ε*which we may likewise write in vector form as

_{y}**= (**

*ε**ε*,

_{x}*ε*). The vector

_{y}**is the projection of the principal wavevector**

*ε***K**onto the

*x*–

*y*plane (divided by

*K*). Thus, the tilt angle with respect to the

_{z}*z*direction amounts to

*φ*= tan

^{−1}

*ε*where

*ε*= |

**|. The geometry is illustrated in Fig. 1(a).**

*ε*Having established the form of the decentered Gaussian beam at the focal plane (*z* = 0), we can find the full form of the beam by invoking a paraxial approximation and integrating for all *z* in Eq. (1). We perform the paraxial approximation as follows

Since decentered Gaussian beams are tilted beams, the ramifications of the tilt on the validity of the paraxial approximation is of particular interest. Specifically, the approximation error is given by *δ* = (1 − *χ*)^{(1/2)} − (1 − *χ*/2), where *χ* stands for the term in the brackets in Eq. (7). Using
$\alpha =\beta \approx \sqrt{2}/{w}_{0}$[17] and
$\left({K}_{x}+{K}_{y}\right)\approx \sqrt{2}{k}_{0}\text{sin}\phi $, we can make a conservative estimate for the approximation error as *χ* ≤ 4tan*φ*/(cos*φ* · (2*πw*_{0}/*λ*))+4/(cos*φ* ·(2*πw*_{0}/*λ*))^{2} where *φ* is the tilt angle and *w*_{0}/*λ* is the Gaussian waist parameter (in units of wavelength). As an example, consider a decentered beam with a tilt angle *φ* = 5°, a waist *w*_{0} ≈ 10*λ*, and a wavelength *λ* = 1*μm*. For these parameters, *δ* = 10^{−5}, and we find that the residual phase *K _{z}δz* reaches

*π*at

*z*≈ 50 mm. Since the above assessment is conservative, the preceding analysis helps identify regimes in which the paraxial approximation is strictly valid, however, it is not necessarily limited to these regimes.

Inserting our paraxial approximation, i.e. Eq. (8), into the plane-wave decomposition of the vector potential, i.e. Eq. (1), and integrating, we find the complete propagated form of the decentered Gaussian beam:

*q*=

*q*

_{0}+

*z*is the Gaussian beam

*q*-parameter at

*z*, and $w={w}_{0}\sqrt{1+{\left(z/{z}_{R}\right)}^{2}}$ is the corresponding Gaussian waist. The Gouy phase is given by

*ϕ*= tan

^{−1}(

*z/z*). The beam center coordinates in the

_{R}*x*–

*y*plane at propagation distance

*z*are expressed as

*x*=

_{d}*x*

_{d}_{0}+

*ε*and

_{x}z*y*=

_{d}*y*

_{d}_{0}+

*ε*which may be combined into the vector

_{y}z**r**

*= (*

_{d}*x*,

_{d}*y*) =

_{d}**r**

_{d}_{0}+

*ε**z*. We see that the full decentered Gaussian beam given in Eq. (9) takes on a very intuitive form. It resembles a conventional Gaussian beam with a decentered focal spot and a propagation direction tilted with respect to the

*z*-axis as illustrated in Fig. 1(a).

Before proceeding, let us make a brief note concerning polarization. We have described the decentered Gaussian beam with vector potential **A*** _{dG}* pointing in the direction of a cartesian unit vector. If

**A**

*points in some direction in the*

_{dG}*x*–

*y*plane, then we expect the corresponding electromagnetic fields to resemble a linearly polarized decentered Gaussian beam by analogy with the familiar Gaussian beam [16, 17]. However with

**A**

*pointing in the*

_{dG}*z*-direction, we expect a radially polarized decentered Gaussian beam now by analogy with the radially polarized Gaussian beam [18].

#### 2.2. Superposition to cylindrical Bessel-Gauss beams

Bessel-Gauss beams are formed by superposing decentered Gaussian beams [15]. As shown in Fig. 1(a), the centers of the decentered Gaussian beams are placed on a circle with a constant radius *r _{d}*

_{0}= |

**r**

_{d}_{0}| in the

*x*–

*y*plane, and the tilt angles of the beam directions point to the apex of a single cone. Specifically, this corresponds to a superposition of beams whose vectors

**r**

_{d}_{0}and

**are collinear and where the parameters**

*ε**r*

_{d}_{0}and

*ε*are constant.

To carry out the superposition, we introduce cylindrical coordinates (*r*, *θ*, *z*). The focus of each decentered Gaussian beam in the superposition is given by coordinates (*r _{d}*

_{0},

*γ*, 0). The angle between the vectors

**r**and

**r**

_{d}_{0}is (

*θ*−

*γ*). Consequently, we can rewrite Eq. (9) in cylindrical coordinates as

**A**

*is the vector potential of the decentered Gaussian beam with its focus at coordinates (*

_{dG,γ}*r*

_{d}_{0},

*γ*, 0), and we should recall that

**r**

*= (*

_{d}*x*,

_{d}*y*) =

_{d}**r**

_{d}_{0}+

*ε**z*. Note that one difference between Eq. (10) and our preceding expression for the decentered Gaussian beam, Eq. (9), is that in Eq. (10) we allow the direction of the vector potential

**n̂**be a function of the beam focal position, specifically the parameter

*γ*. A variation of

**n̂**throughout the superposition enables us to generate Bessel-Gauss beams with cylindrical polarization states.

Finally, let us carry out the superposition of the component beams as defined in Eq. (10). The superposition constitutes an integration of *γ* from 0 to 2*π* and yields a solution for the vector potential of the cylindrical Bessel-Gauss beam:

**n̂**oriented along the radial or azimuthal direction. As we will see in section 3, this will produce electromagnetic fields representing radially and azimuthally polarized Bessel-Gauss beams. Next, we consider

**n̂**oriented along a cartesian unit vector. Here, we generate expressions for the well-known linearly polarized Bessel-Gauss beam as well as an alternative formulation for a radially polarized Bessel-Gauss beam. For these integrations, we apply the definition of the Bessel-function, e.g. as given in [3]:

*J*

_{−}

*(*

_{k}*τ*) = (−1)

*(*

^{k}J_{k}*τ*).

### 2.2.1. **n̂** with radial or azimuthal orientation direction

Consider the decentered Gaussian component beam with its focal spot given by coordinates (*r _{d}*

_{0},

*γ*, 0). In the superposition

*r*

_{d}_{0}is fixed and

*γ*is variable. With local linear polarization

**n̂**in the radial or azimuthal directions,

**n̂**= cos

*γ*

**x̂**+ sin

*γ*

**ŷ**or

**n̂**= −sin

*γ*

**x̂**+ cos

*γ*

**ŷ**, respectively. Expanding the trigonometric functions into exponential functions and inserting these expressions into the cylindrical Bessel-Gauss beam integral, Eq. (11), the integration can easily be performed by applying Eq. (12). The superposition is illustrated in Fig. 1(b). The result is given by:

**n̂**=

**r̂**or

**, and the vector potential amplitude**

*θ̂**ψ̃*(

*r*,

*z*) is given by:

*J*

_{1}-Bessel function. Similar expressions have previously been reported for azimuthal and radial polarization [11, 12]; however, one advantage of superposing decentered Gaussian beams to radially and azimuthally polarized Bessel-Gauss beams is that the propagation of these component beams through first-order optical systems, i.e. ABCD optical systems, is known. Thus, a solution for radially or azimuthally polarized Bessel-Gauss beams in ABCD optical systems can be derived as well. This will be discussed in the appendix.

Figs. 2 (a) and (b) show this beam’s vector potential amplitude in the focus and far-field, respectively. The parameters are *w*_{0} = 1 mm, *ε* = 0.45°, *r _{d}*

_{0}= 0, and wavelength of

*λ*= 1050 nm (corresponding to the experimental configuration described in section 4). The propagation of the amplitude through the focus is shown in Fig. 2 (c).

### 2.2.2. **n̂** along the Cartesian-unit vector

Consider the case for which **n̂** is orientated along any cartesian unit vector (**x̂**, **ŷ**, **z ̂**) or any linear combination thereof. For this case, the superposition is illustrated in Fig. 1(c), and the integral of Eq. (11) is given by:

**n̂**=

**x̂**or

**ŷ**this represents the well-known linearly polarized Bessel-Gauss beam solution. For

**n̂**=

**ẑ**, we will see in section 3 that this expression represents a radially polarized Bessel-Gauss beam in the same form as that described in the preceding subsection. Note that the cartesian unit vector solutions are determined by

*J*

_{0}instead of

*J*

_{1}. Figs. 2 (d) and (e) show this beam’s vector potential amplitude in the focus and far-field, respectively. Fig. 2 (f) displays the complete beam propagation through the focus. In the focus this beam solution shows an on-axis amplitude. The beam parameters are the same as those for the preceding subsection.

Before closing this section, it is worth noting that this concept of a ‘ray’ bundle of decentered Gaussian beams, i.e. our superposition procedure, can be extended to describe more general beam-shapes. For example, the component beams do not need to be continuously placed over the circle but can be situated discretely [19]. In this way, the output of recently demonstrated of fiber-laser arrays can be described [20].

## 3. Electromagnetic fields of cylindrical Bessel-Gauss beams

In this section we obtain the electromagnetic fields of cylindrical Bessel-Gauss beams from the vector potential description of the previous section. We highlight the utility of the *z*-oriented solution in describing radially polarized Bessel-Gauss beams. Finally, we compare the fields derived from our analytical results to non-paraxial numerical simulations.

#### 3.1. Electromagnetic fields and the vector potential

The electric and magnetic field are obtained from the vector potential **A**(*x,y,z,t*) and scalar potential *φ*(*x,y,z,t*) by:

**A**need not be zero (unlike ∇ ·

**E**). If we choose the Lorenz gauge condition: it simply defines the scalar potential. Considering monochromatic light with angular frequency

*ω*

_{0}, we rewrite the fields as

**E**(

*x,y,z,t*) is expressed by a term proportional to the vector potential, as well as a gradient term. The latter is required to ensure that the divergence of the electric field vanishes.

#### 3.2. Fields of cylindrical beams

In section 2, we formulated cylindrical Bessel-Gauss beams in terms of the vector potential. Specifically, the vector potential descriptions took the form:

with*ψ*=

*ψ*

^{(z)}for

**n̂**=

**ẑ**and

*ψ*=

*ψ*

^{(r)}for

**n̂**=

**r̂**,

**being the cases of interest for this work as these lead to radial or azimuthal polarization. Using the field expressions in Eq. (19), we can tabulate expressions for the field components in these cases:**

*θ̂*Inspecting the expressions in Tables 1 and 2, we see the two cases correspond to *transverse electric* and *transverse magnetic* beams. Considering the azimuthal orientation ** θ̂** of the vector potential, we see that the resulting electric field in Table 1 is purely azimuthally polarized while the magnetic field possesses transverse as well as longitudinal components. Such beams are thus referred to as

*transverse electric*beams. Orientation of the vector potential in the radial

**r̂**and

**ẑ**directions produces radially polarized electric fields, as shown in Table 2. For this case, the magnetic field does not have a longitudinal component but the electric field does. This case corresponds to

*transverse magnetic*beams.

Considering further the field expressions in Table 2, it is intriguing to note that there are two cylindrical Bessel-Gauss solutions that offer radial polarized fields. Specifically, the radial polarization can either arise from a radially oriented vector potential with amplitude *ψ*^{(r)} or it can be derived from a *z*-oriented vector potential with amplitude *ψ*^{(z)}. In section 2, we have derived two different explicit paraxial expressions for *ψ*^{(r)} and *ψ*^{(z)}. It is then reasonable to inquire how similar or different the resulting radially polarized electromagnetic fields are. Considering the magnetic field of the beam arising from the radially oriented vector potential, i.e. the first row of Table 2, and utilizing *ψ*^{(r)} in the form of Eq. (14), we find *B _{θ}* ≈

*iK*

_{z}ψ^{(r)}(the primary

*z*variation in

*ψ*

^{(r)}follows from the fast-varying exp(

*iK*) term). For the second radially polarized beam solution arising from the

_{z}z*z*-oriented vector potential, we employ

*ψ*

^{(z)}in the form of Eq. (16), and noting that the primary radial variation of

*ψ*

^{(z)}should come from the Bessel function term, it is straightforward to show that near the focus i.e.

*z*≪

*z*and with small decenter parameter i.e.

_{R}*r*

_{d}_{0}≪

*z*the magnetic field takes on the form

_{R}*B*≈

_{θ}*K*

_{z}ε*ψ*

^{(r)}. So, for this region of interest these two radially polarized solutions appear virtually identical. Although, these two solutions look nearly identical in this region, we have not thoroughly explored their properties outside this region. Such an investigation merits further study.

We should also emphasize that the amplitudes *ψ*^{(r)} and *ψ*^{(z)}, which are used for the general field expressions in Tables 1 and 2, obey two different wave equations. As mentioned, in the Lorenz gauge and in a current-free medium, the vector potential **A**(*x,y,z,t*) satisfies the vector wave equation:

*z*-oriented vector potential, i.e. Eq. (20) with

*ψ*=

*ψ*

^{(z)}and

**n̂**=

**ẑ**, the vector wave equation reduces simply to:

*x*–

*y*plane; thus, this familiar wave equation is also used to describe linearly polarized beams. However, for the radially oriented vector potential, i.e. Eq. (20) with

*ψ*=

*ψ*

^{(r)}and

**n̂**=

**r̂**or

**, the vector wave equation reduces to:**

*θ̂**ψ*

^{(r)}/

*r*

^{2}term follows from the vector Laplacian in the vector wave equation. Eq. (23) has been used previously to study azimuthally and radially polarized Bessel-Gauss beams [11, 12]; however, it does not share the familiar form of Eq. (22) and does not allow straightforward plane-wave solutions. This is a key point; we have formulated expressions for radially polarized Bessel-Gauss beams by superposing component beams following the familiar wave equation for Cartesian unit vector orientation. In this way, we have indirectly solved for the more complicated wave equation in the form of Eq. (23). The developed radially polarized Bessel-Gauss beam expressions (at

*z*= 0) serve as an input for non-paraxial, numerical propagation through Fourier transform techniques, as described in the following.

#### 3.3. Non-paraxial numerical propagation

In this section we show that the vector potential can also be applied to numerically calculate the propagation of radially polarized beams. The *z*-oriented vector potential is advantageous for numerical calculations as the corresponding wave equation, Eq. (22), can be solved with a superposition of plane waves as given by Eq. (1). The electromagnetic fields are calculated from the numerical solution for the vector potential by using Eq. (19).

In a numerical implementation, we do not need to invoke a paraxial approximation, i.e. *k _{z}* does not need to be approximated by a parabolic function, as performed in Eq. (8). This approximation was required to obtain an analytical solution for the Bessel-Gauss beam in form of Eq. (15) and Eq. (16). Thus, the numerical solution is directly calculated by using Eq. (1), which we can rewrite in terms of Fourier-transforms (

*FT*) and inverse Fourier-transforms (

*IFT*):

*ψ*

^{(z)}again refers to the amplitude of the

*z*-oriented vector potential. We only need to know the initial amplitude

*ψ*

^{(z)}(

*x,y,z*= 0,

*t*= 0) in order to obtain the complete vector potential amplitude

*ψ*

^{(z)}(

*x,y,z,t*), and consequently, the electromagnetic fields by way of Eq. (19). Additionally, we should note that since the cylindrical beams are azimuthally symmetric, the Fourier-transform in Eq. (24) is advantageously calculated by using a quasi-discrete Hankel-transform [21]. For radially polarized light, this is a fast alternative to the commonly employed integral representation of the field in the focus [22].

For the example, we simulate the evolution of a radially polarized Bessel-Gauss beam with semi-aperture angle *φ* = 60° (= 1.05 rad), a Gaussian waist parameter of *w*_{0} = 10 *μ*m, and a wavelength of *λ* = 1 *μ*m. Fig. 3 (a) shows the initial state, and the numerical simulation of the propagation is shown in Fig. 3 (b). The radial and *z*-component of the electric fields are calculated by applying Eq. (19). Close-ups of the instantaneous fields are shown in Figs. 3(c) and (d). It can be seen that the radial component shows zero amplitude in the beam center, while the *z*-component is the strongest there. For the same beam parameters, in Fig. 4(a) we compare the numerically calculated *z*-component of the electric field to the analytically derived field at the beam center (*r* = 0), and in Fig. 4(b) we show a comparison of the corresponding radial electric fields at the off-axis maximum (at about *r* = 0.75 *μ*m). Note that the *z*-component of the electric field exceeds the radial component for the present semi-aperture angle of *φ* = 60°. This effect has previously been described [8–10].

We additionally numerically calculate the ratio of the maximum *z*-component of the electric field to the maximum radial component of the electric field for radially polarized Bessel-Gauss beams as the semi-aperture angle of the beam is varied. The results of this computation are given in Fig. 4(c). We see that the *z*-component exceeds the radial component for large semi-aperture angles *φ*. The beam parameters are *w*_{0} = 10 *μ*m and *λ* = 1 *μ*m for this simulation; however, the results look virtually identical for other beam waists. Furthermore, we find that the results from the analytical solution parallel the numerical results. This demonstrates that the analytical model based on Eq. (16), which is a result of the approximation Eq. (8), holds even for large semi-aperture angles. Consequently, the error of the specially developed paraxial approximation, Eq. (8), is minimal.

## 4. Experiment

For experimental investigation, we have designed an optical system that permits the production of radially polarized Bessel-Gauss beams. The beam generation is based on a specialty optical fiber in which a long period grating efficiently converts a linearly polarized waveguide mode to a radially polarized higher-order waveguide mode [23]. Compared to free-space optics approaches, such as liquid crystal polarization converters or space-variant waveplates [24,25], the main advantage of the fiber-based approach is that it permits a fully integrated solution to the transformation from a linearly polarized beam to a radially polarized donut beam. Beam pointing stability and remote delivery of the generated beams are additional beneficial attributes.

A schematic of the setup is shown in Fig. 5(a). We use a tunable grating-stabilized external cavity diode laser with a wavelength of 1050 nm and a line width of 100 kHz. The laser is attached to a single-mode, polarization-maintaining fiber (PM980-XP, mode-field diameter of approximately 6.6 *μ**m*); at the fiber output, we obtain 30 mW of power. The light is then free-space coupled (using lenses L1 and L2 with focal lengths f = 4 mm) through a half-wave plate (HWP) and a quarter-wave plate (QWP) for polarization control and into a fiber module comprising a single mode fiber (SM980) spliced to a specialty optical fiber [23]. The refractive index profile of this fiber is shown in Fig. 5(b). The linearly polarized fundamental HE_{11} fiber mode is mainly located in the inner core of the fiber, and a radially polarized TM_{01} fiber mode is dominantly located in the outer index ring. The calculated effective area of this TM_{01} mode at a wavelength of 1050 nm is A* _{eff}* = 51

*μ*m

^{2}. The long period mircobend grating (LPG) converts the modes by matching the effective indices of the fundamental HE

_{11}fiber mode to the radially polarized TM

_{01}fiber mode [23]. The LPG is inscribed in this specialty fiber by using a LDS3 fusion splicer [26]. The LPG has a grating period of 762

*μ*m and a length of about 23 mm (30 periods). After the mode conversion using the LPG, the power ratio

*P*

_{HE11}/

*P*

_{TM01}is measured to be −19 dB at a wavelength of 1050 nm for the employed device. Thus, only approximately 1% of the power is left in the fundamental HE

_{11}fiber mode after conversion. These values are obtained with an independent measurement of the LPG spectrum, as shown in Fig. 5(c). Specifically, this measurement is obtained by splicing the output of the specialty optical fiber to a single mode optical fiber.

We magnify the TM_{01} mode coming out of the fiber by about 60 times with a single lens L3 with f = 4 mm. Fig. 6(a) shows an image of the magnified mode that has been recorded in plane P1. The mode clearly shows an annular mode shape with some unconverted, residual HE_{11} power still visible in the center. With a polarizer situated directly before plane P1, the parts of the mode aligned with the polarizer are transmitted. This projection leads to a double crescent-shaped intensity profile rotating with the orientation of the polarizer as illustrated in Fig. 6(b). This confirms the radial polarization of the TM_{01} mode. Removing the polarizer, the mode is further magnified 10 times with a telescope (lenses L4 and L5 with focal lengths of 10 mm and 100 mm, respectively). The magnified mode is then sent through an axicon to produce a radially polarized Bessel-Gauss beam. Fig. 7(a) shows the partial beam profile at the position of the axicon, plane P2. Note that the whole beam cannot be recorded simultaneously due to limited size of the camera chip. The beam exhibits a ring radius of *r _{d}* = 2350

*μ*m and a Gaussian waist parameter of

*w*≈ 1050

*μ*m. These parameters are determined by fitting the intensity profile to a Bessel-Gauss beam in the far-field [1]. The imposed semi-aperture angle by the axicon is given by 0.45°. The on-axis crossing point of the Bessel-Gauss beam occurs at a distance

*r*/ tan(0.45°) ≈ 300 mm behind the axicon (plane P3). Fig. 8(a) shows the image acquired in plane P3, and Fig. 8(b) shows a comparison between the cross-section through the experimental data and our analytical expressions for radially polarized Bessel-Gauss beams. The analytical form is based on the electric field of the

_{d}*z*-oriented vector potential Bessel-Gauss solution, where the fields are calculated with Eq. (16) and Eq. (19). The analytical expression agrees very well with the measured beam profile.

As a final note concerning our experimental work, we should point out that axicons are a common means to produce Bessel-Gauss beams; however, fabrication tolerances limit the quality of axicon-produced Bessel-Gauss beams [27]. We have measured the relative optical path length of the axicon with a wavefront sensor (imagine optics) permitting a *λ*/100 rms absolute measurement accuracy while having a high dynamic range of hundreds of wavelengths. It can be seen in Fig. 7(b) that the axicon shows a round tip instead of a conic apex due to fabrication tolerances of diamond turning. The optical path impressed by the real axicon can be well approximated by a hyperbolic function
$\left(n-1\right)\sqrt{{\left(\delta \cdot r\right)}^{2}+{\left(10.5\mu m\right)}^{2}}$, where *n* is the refractive index of axicon substrate (fused silica, *n* = 1.45) and *δ* is the base angle of the axicon, *δ* = 1°. In contrast to commonly used Gaussian illumination of axicons, the beam here possesses a ring shape, and thus, it avoids the round center of the axicon. This is why our measured beam profiles so precisely match our analytical Bessel-Gauss expressions. This technique may be useful for producing cleaner Bessel-Gauss beams in the future.

## 5. Conclusion

In conclusion, we have obtained expressions for radially and azimuthally polarized Bessel-Gauss beams by superposing decentered Gaussian beams with differing polarization states. This description can be applied to ABCD optical systems, which will be useful for the design of resonators for such beams. Additionally, we have developed a nonparaxial, fully numerical beam propagator for radially polarized, azimuthally symmetric beams. We have employed this to investigate radially polarized Bessel-Gauss beams and demonstrate the accuracy of our analytical expressions even for large semi-aperture angles. Finally, the form of the analytical solution has been experimentally verified by employing a fiber-based mode-converter generating an annular mode to produce a radially polarized Bessel-Gauss beam. Our radially polarized Bessel-Gauss beam generation scheme avoids the imperfect center of the axicon, and thus, results in a high-quality, radially polarized Bessel-Gauss beam.

## Appendix

## Propagation of cylindrical Bessel-Gauss beams in ABCD optical systems

In the following, the solutions for radially and azimuthally polarized Bessel-Gauss beams are generalized for propagation in ABCD optical systems. The analytical description of decentered Gaussian beams in ABCD optical systems is known [14,15], and their essential form is identical to Eq. (9). A superposition of these beams to radially and azimuthally polarized Bessel-Gauss beams is then analogous to section 2.2. For an orientation of **n̂** along a Cartesian unit vector, the results agree with previously derived expressions [15, 28]. If the orientation of **n̂** is in the radial or azimuthal direction, we find the following expression:

**n̂**stands for either radial or azimuthal unit vector, and the vector potential amplitude is given by:

*q*,

*w*and

*ϕ*transform like that of conventional on-axis Gausian beams, as demonstrated by Eq. (27). The formulation of cylindrical Bessel-Gauss beams in terms of ABCD-matrix parameters allows us to study their transformation under the influence of optical components.

It is worth noting that by comparing Eq. (26) with Eq. (14) one sees that *K _{z}* is replaced by

*k*

_{0}. This change is required for consistency with the Collins integral [29]. Additionally,

*K*depends on the inclination angle which changes during propagation through an ABCD optical system. In fact, for propagation in a uniform medium, the phase-term Σ in Eq. (26) can be regarded as a correction of

_{z}*k*

_{0}to approximate

*K*. However, for such propagation the description based on

_{z}*K*in form of Eq. (16) is more accurate.

_{z}## Acknowledgments

The authors acknowledge the contributions of Gilberto Abram to the implementation of the quasi-discrete Hankel-transform, as well as Jan Schulte for helpful discussions. We thank Imagine Optic, Inc. for the loan of the wavefront sensor (model HASO3). This work is supported by AFOSR grant FA9550-12-1-0499 and the Center for Free-Electron Laser Science. William P. Putnam acknowledges support from the NSF graduate research fellowships program.

## References and links

**1. **C. J. R. Sheppard and T. Wilson, “Gaussian-beam theory of lenses with annular aperture,” IEE J. Microw. Opt. Acoust. **2**, 105–112 (1978) [CrossRef] .

**2. **F. Gori, G. Guattari, and C. Padovani, “Bessel-Gauss beams,” Opt. Commun. **64**, 491–495 (1987) [CrossRef] .

**3. **V. Bagini, F. Frezza, M. Santarsiero, G. Schettini, and G. Spagnolo Schirripa, “Generalized Bessel-Gauss beams,” J. Mod. Opt. **43**: 6, 1155–1166 (1996).

**4. **M. Duocastella and C. B. Arnold, “Bessel and annular beams for materials processing,” Laser Photonics Rev. **6**, 607–621 (2012) [CrossRef] .

**5. **V. Garcés-Chávez, D. McGloin, H. Melville, W. Sibbett, and K. Dholakia, “Simultaneous micromanipulation in multiple planes using a self-reconstructing light beam,” Nature **419**, 145–147 (2002) [CrossRef] [PubMed] .

**6. **W. P. Putnam, D. N. Schimpf, G. Abram, and F. X. Kärtner, “Bessel-Gauss beam enhancement cavities for high-intensity applications,” Opt. Express **20**, 24429–24443 (2012) [CrossRef] [PubMed]

**7. **R. D. Romea and W. D. Kimura, “Modeling of inverse Cerenkov laser acceleration with axicon laser-beam focusing,” Phy. Rev. D **42**, 1807–1818 (1990) [CrossRef] .

**8. **S. Quabis, R. Dorn, M. Eberler, O. Glöckl, and G. Leuchs, “Focusing light to a tighter spot,” Opt. Commun. **179**, 1–7 (2000) [CrossRef] .

**9. **K. Youngworth and T. Brown, “Focusing of high numerical aperture cylindrical-vector beams,” Opt. Express **7**, 77–87 (2000) [CrossRef] [PubMed] .

**10. **E. Y. S. Yew and C. J. R. Sheppard, “Tight focusing of radially polarized Gaussian and Bessel-Gauss beams,” Opt. Lett. **32**, 3417–3419 (2007) [CrossRef] [PubMed] .

**11. **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) [CrossRef] [PubMed] .

**12. **D. G. Hall, “Vector-beam solutions of Maxwell’s wave equation,” Opt. Lett. **21**, 9–11 (1996) [CrossRef] [PubMed] .

**13. **P. L. Greene and D. G. Hall, “Properties and diffraction of vector Bessel-Gauss beams, ” J. Opt. Soc. Am. A **15**, 3020–3027 (1998) [CrossRef] .

**14. **A. R. Al-Rashed and B. E. A. Saleh, “Decentered Gaussian beams,” Appl. Opt. **34**, 6819–6825 (1995) [CrossRef] [PubMed] .

**15. **D. N. Schimpf, J. Schulte, W. P. Putnam, and F. X. Kärtner, “Generalizing higher-order Bessel-Gauss beams: analytical description and demonstration,” Opt. Express **20**, 26852–26867 (2012) [CrossRef] [PubMed] .

**16. **L. W. Davis, “Theory of electromagnetic beams,” Phy. Rev. A **19**, 1177–1179 (1979) [CrossRef] .

**17. **H. A. Haus, *Waves and Fields in Optoelectronics* (CBLS, 2004).

**18. **L. W. Davis, “TM and TE electromagnetic beams in free space,” Opt. Lett. **6**, 22–23 (1981) [CrossRef] [PubMed] .

**19. **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] .

**20. **R. S. Kurti, K. Halterman, R. K. Shori, and M. J. Wardlaw, “Discrete cylindrical vector beam generation from an array of optical fibers,” Opt. Express **17**, 13982–13988 (2009) [CrossRef] [PubMed] .

**21. **Li Yu, M. Huang, M. Chen, W. Chen, W. Huang, and Z. Zhu, “Quasi-discrete Hankel transform,” Opt. Lett. **23**, 409–411 (1998) [CrossRef] .

**22. **B. Richards and E. Wolf, “Electromagnetic diffraction in optical systems. II. Structure of the image field in an aplanatic system,” Proceedings of the Royal Societyof London. Series A, Mathematical and Physical Sciences **253**, 358–379 (1959) [CrossRef] .

**23. **S. Ramachandran, P. Kristensen, and M. F. Yan, “Generation and propagation of radially polarized beams in optical fibers,” Opt. Lett. **34**, 2525–2527 (2009) [CrossRef] [PubMed] .

**24. **M. Stalder and M. Schadt, “Linearly polarized light with axial symmetry generated by liquid-crystal polarization converters,” Opt. Lett. **21**, 1948–1950 (1996) [CrossRef] [PubMed] .

**25. **G. Machavariani, Y. Lumer, I. Moshe, A. Meir, and S. Jackel, “Efficient extracavity generation of radially and azimuthally polarized beams,” Opt. Lett. **32**, 1468–1470 (2007) [CrossRef] [PubMed]

**26. **I. K. Hwang, S. H. Yun, and B. Y. Kim, “Long-period fiber gratings based on periodic microbends,” Opt. Lett. **24**, 1263–1265 (1999) [CrossRef] .

**27. **O. Brzobohaty, T. Cizmar, and P. Zemanek, “High quality quasi-Bessel beam generated by round-tip axicon,” Opt. Express **16**, 12688–12700 (2008) [CrossRef] [PubMed] .

**28. **M. Santarsiero, “Propagation of generalized Bessel-Gauss beams through ABCD optical systems,” Opt. Commun. **132**, 1–7 (1996) [CrossRef] .

**29. **S. A. Collins Jr., “Lens-system diffraction integral written in terms of matrix optics,” J. Opt. Soc. Am. **60**, 1168–1177 (1970) [CrossRef] .