## Abstract

We report the generation of a zero-order Bessel beam of continuously variable spot size using a simple optical setup. We have used a pair of metal axicon mirrors to generate a hollow beam of variable dark diameter. This beam was subsequently focused by a convex lens to get a Bessel beam of variable spot size. We also studied the effect of a hollow-beam ring width on nondiffracting propagation range of the generated beam.

© 2012 Optical Society of America

## 1. Introduction

Bessel beams have received significant attention since their conception and realization by Durnin *et al.* [1]. The size of the central spot as well as the peak intensity of a Bessel beam does not change significantly over a propagation distance. This unique and remarkable property, called nondiffracting property, of Bessel beams has made them useful for various technological applications such as atom guiding [2,3], optical tweezers [4], laser ablation, and laser machining [5–7]. Besides nondiffracting property, Bessel beams also show self-healing property which can be exploited for some purposes [8]. While generation of high-order Bessel beams is possible [3], we consider, in the present work, generation of only zero-order Bessel beams. The cross-section of a zero-order Bessel beam has a bright spot in the center surrounded by concentric rings, with overall intensity distribution given as square of the zero-order cylindrical Bessel function (${J}_{0}$). The intensity of side-rings decreases rapidly as one goes away from the beam center, however power in any ring is equal to that in the central lobe in intensity profile of the beam. These beams also show interesting electromagnetic field structure when spot size of the central lobe is comparable to the wavelength [9].

The electromagnetic field, having transverse variation as ${J}_{0}$-Bessel function, results from the superposition of plane waves with wave vectors distributed over the surface of a cone. The spot size, which is referred to here as full width at half-maximum (FWHM), of the central peak in the transverse intensity profile, is equal to $2.27/\alpha $ for an ideal ${J}_{0}$-Bessel beam; where $\alpha =k\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta $, $k$ is wave-vector magnitude and $\theta $ is half of the vertex angle of the cone formed by the wave vectors in the superposition. It is well known that generation of an ideal ${J}_{0}$-Bessel beam is not possible as it requires plane waves having infinite width as its constituent waves. However, beams having an intensity profile close to ${J}_{0}$-Bessel beam intensity profile have been experimentally generated over finite transverse extent as well as with finite nondiffracting propagation range. Such a beam is referred to here as a generated ${J}_{0}$-Bessel beam (or Bessel beam).

Various optical elements and setup have been used to generate Bessel beams. These include a ring aperture placed in the focal plane of a convex lens [1], binary-coded holograms [10], axicon lenses [11], a converging lens with spherical aberrations [12], and reflective axicon mirrors [13–18]. The reflective axicon mirrors are easy to fabricate and have no aberrations (as compared to refractive elements such as axicon lenses) due to propagation in bulk of refractive medium. The inhomogeneity in bulk medium can be a source of such aberrations.

Generally, the spot size of a Bessel beam is fixed for a given experimental setup. To change the Bessel-beam spot size, one needs to replace the components of the setup with those having modified specifications. However, some active optical devices such as spatial light modulators (SLMs) [19] and tunable acoustic gradient lens [20,21] have been used to generate Bessel beams with variable spot parameters. But, SLMs are more expensive whereas acoustic gradient lenses have a complex design as compared to optical axicons. In an interesting paper, Milne *et al.* [22] have reported use of a fluidic axicon lens to obtain a Bessel beam with variable spot size. This method required refilling of fluid having a different value of refractive index to obtain the changed spot size of the Bessel beam. While changing fluid each time is inconvenient, this approach may also suffer from limitations on continuous variation of spot size due to discrete values of the refractive index of various fluids. Recently, Brousseau *et al.* [23] and Brunne *et al.* [24] have used axicon mirrors whose conical surface could be manipulated by applied magnetic field and applied voltage, respectively, to obtain variable spot parameters of the generated Bessel beam. Earlier in a work by Vaičaitis and Paulikas [25], the generation of a Bessel beam of variable spot size has also been studied by using passive optical elements. In this paper we report the generation of a ${J}_{0}$-Bessel beam of continuously variable spot size using a simple and convenient method with a low-cost setup. We used a setup consisting of a pair of axicon mirrors and a few other standard optical elements. The generated hollow beam (HB) was focused using a convex lens to obtain the ${J}_{0}$-Bessel beam. The Bessel-beam spot size could be varied just by varying the separation between axicon mirrors in the setup. The range over which spot size remains variable can be altered by changing the focal length of the lens. A ${J}_{0}$-Bessel beam with such a continuous variation possible in its spot size can be useful in several applications such as atom guiding, optical tweezers, and laser machining.

## 2. Experimental Setup

Figure 1 shows the schematic of our experimental setup. We used two homemade metal axicon mirrors [${\mathrm{AX}}_{1}$ (convex) and ${\mathrm{AX}}_{2}$ (concave)] of equal angle $\gamma =1\xb0$ for conical surface. These gold-coated diamond-turned copper axicon mirrors have high reflectivity ($\sim 93\%$) at 633 nm laser wavelength which we used in the experiments. A linearly polarized (vertical-direction) laser beam having a nearly Gaussian spatial profile with FWHM $\sim 0.4\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ and power $\sim 2.5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mW}$ from an H–Ne laser system was used as an input laser beam. Appropriate polarization-selective optical elements such as a polarizing-cube beam splitter (PBS) and quarter-wave ($\lambda /4$) plates (Fig. 1) were used for coupling the input laser beam to axicon mirrors. The collimated HB generated from the setup was focused using a convex spherical lens ($L$) of $\sim 1000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ focal length to obtain the ${J}_{0}$-Bessel beam. A digital CCD camera (Pixelfly, PCO, Germany) having $1392\times 1024$ pixels (pixel-size $6.45\times 6.45\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{\mu m}$), in conjunction with peripheral component interconnect interface-board on a computer, was used to record the beam images at different positions. We have recorded an image of a HB by taking an image of the screen on which the HB was incident. Direct imaging of the HB was not possible because its size was larger than the size of the CCD element. Therefore surface quality of the screen used may also affect the intensity profiles of the recorded images. On the other hand, for a Bessel beam, images were recorded keeping CCD directly in the beam path. The transverse intensity profiles of the beams were obtained from the recorded images. The overall conversion efficiency from input laser-beam power to ${J}_{0}$-Bessel-beam power in our experiments was $\sim 73\%$. This efficiency can be further increased by using high-reflectivity dielectric coating on the metal axicon mirrors.

After the lens (L), the generated HB gets converged such that wave vectors of different rays form a conical surface with cone axis the same as optical axis ($z$-axis) of the lens (Fig. 1). The region of this conical superposition of waves [shown by shaded part in Fig. 1] is the region of formation of the Bessel beam with average direction of propagation along the $z$-axis. The transverse intensity variation in this region is approximately a square of ${J}_{0}$-Bessel function. Here we note that angle $\theta $, which governs the Bessel-beam spot size, is dependent on the dark diameter of the HB at lens (L). The dark diameter can be varied by varying the separation ($d$) between axicon mirrors in the setup. This provided the opportunity to vary continuously the spot size of the generated Bessel beam just by varying the mirror separation ($d$).

## 3. Results and Discussion

#### A. Variation of FWHM with Axicon Mirror Separation ($\mathsf{d}$)

We recorded images of generated Bessel beam at different values of position $z$ and for different values of mirror separation $d$. Corresponding to different values of $d$, the images of the HB were also recorded at lens plane ($z=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$). Figures 2 and 3 show the CCD images and intensity profiles of the generated HB (at the lens plane $z=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$) and Bessel beam (at $z=1230\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$) respectively when we kept mirror separation of $d=90\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ in our setup. The graphs shown in Figs. 2 and 3 are the plot of intensity profiles of hollow and generated Bessel beams across the horizontal diameter in the beam cross-section. Nearly identical profiles were observed across the vertical diameter of the beam, indicating a good circular symmetry in the generated beam. In Fig. 3(b), the curve having variation as square of ${J}_{0}$-Bessel function is also shown for comparison with experimentally observed profile.

From the recorded CCD images, we measured the dark diameter of the generated HB for various values of separation $d$. In the measurements, we have taken dark diameter as separation between nearest half-maximum intensity points (on the intensity profile) which are diametrically opposite about the beam center. We found that measured dark diameter was nearly equal to the geometrically estimated dark diameter ($\varphi $) given by $\varphi =2d\mathrm{tan}(2\gamma )$, and its variation with mirror separation $d$ was linear as shown by hollow circles in Fig. 4. Similar results have been observed by us earlier also [18]. However, this striking dependence of hollow-beam diameter on mirror separation has been exploited in the present work to generate the Bessel beam of variable spot size after focusing the HB using the lens. Owing to dependence of Bessel-beam spot size on hollow-beam diameter in our geometry (Fig. 1), we could obtain a continuous variation in the spot size of the generated Bessel beam by varying the mirror separation $d$, as shown in Fig. 4. In this figure, data show the measured variation in hollow-beam diameter (hollow circles) at the lens plane (position $z=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$) and corresponding variation in Bessel beam spot size (filled circles) measured at $z=1150\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ with separation $d$. It is clearly evident from the data in this figure that a larger mirror separation resulted in a larger hollow-beam diameter and a smaller Bessel-beam spot size. Though data shown in Fig. 4 are for some selected values of $d$, using our setup we could vary the Bessel-beam spot size continuously over a range by varying the mirror separation $d$. We note that for our geometry, Bessel-beam spot size is given as $\mathrm{FWHM}=2.27/(k\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta )\sim 2.27/(k\hspace{0.17em}\mathrm{tan}\text{\hspace{0.17em}}\theta )\sim 2.27/(k\varphi /2{z}_{\mathrm{BB}})\sim 4.54f/(k\varphi )$, where ${z}_{\mathrm{BB}}$ ($\sim f$, i.e. focal length of lens) is the distance of the Bessel-beam formation region from the lens and $\varphi =2d\mathrm{tan}(2\gamma )$ is the geometrically estimated dark diameter of the HB. Hence FWHM of a Bessel beam is expected to vary inversely with mirror separation $d$ as well as with diameter $\varphi $. Our data presented in Fig. 5 show similar behavior. A small variation in $d$ can result in variation in Bessel-beam spot size which is given approximately as $\mathrm{\Delta}(\mathrm{FWHM})\sim \mathrm{FWHM}(\mathrm{\Delta}d/d)$, where $\mathrm{\Delta}$ refers the variation. Thus a variation $\mathrm{\Delta}d=1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ (at $d=70\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$) is expected to change FWHM of Bessel beam by value $\sim 1.5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{\mu m}$. However, such a small change is difficult to sense in our setup as resolution of our CCD is $\sim 6.45\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{\mu m}$ (size of a pixel). Nevertheless, such a small variation in Bessel-beam spot size can be measured by other more accurate methods [16].

We note that for our geometry, Bessel-beam spot size varies with focal length $f$ and dark diameter $\varphi $ which is given as $\mathrm{FWHM}\sim 4.54f/(k\varphi )$. For a given value of focal length $f$, the range over which FWHM can be varied is governed by the cross-section area of optical components (i.e. size of optics) to accommodate the generated HB and the minimum achievable separation $d$ in the setup. In our setup using $f=1000\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$, we could vary spot size from $\sim 50\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{\mu m}$ to $\sim 100\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{\mu m}$ by varying mirror separation $d$ from 130 to 70 mm. As FWHM is varying linearly with focal length $f$, the spot size of the generated Bessel beam can be further reduced if a lens of smaller focal length is used. We have earlier demonstrated the generation of a Bessel beam of FWHM smaller than the size of a CCD pixel by using a shorter focal length lens [16]. Thus the range over which the spot size can be varied by varying separation $d$ can be further extended by replacing the lens in the setup.

#### B. Variation in Nondiffracting Propagation Range with Axicon Mirror Separation ($\mathsf{d}$)

The variation in peak intensity of a Bessel beam with propagation distance is another important issue. The range over which peak intensity as well as spot size remain unchanged is called nondiffracting propagation range (${Z}_{\mathrm{PR}}$). For its experimental determination, we define ${Z}_{\mathrm{PR}}$ as a range along $z$-axis over which variation in both the parameters (i.e., intensity and spot size) is within 10% from their values at the peak intensity position. For our geometry, ${Z}_{\mathrm{PR}}$ is expected to scale as the length of Bessel beam formation region (and is given as ${Z}_{\mathrm{PR}}\sim 2w/\mathrm{sin}\text{\hspace{0.17em}}\theta \sim 2w/\mathrm{tan}\text{\hspace{0.17em}}\theta \sim 4wf/\varphi $, for $2w\ll \varphi $). Here $w$ is ring-width parameter defined as $\sim 0.85\text{\hspace{0.17em}}$ times the FWHM of the Gaussian intensity profile of the ring of the HB [expressed mathematically later in Eq. (3)]. Figure 6 shows the measured variation in peak intensity of the Bessel beam (hollow circles) generated in our experiments with distance $z$, for some selected different values of separation $d$. A faster variation in peak intensity with $z$ for the larger $d$ value indicates the smaller value of ${Z}_{\mathrm{PR}}$ for larger $\varphi $. The experimentally measured ${Z}_{\mathrm{PR}}$ ($\sim 200\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ for $d=70\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$) for the generated Bessel beam was much larger (approximately five times) than the Rayleigh range for a Gaussian beam of equal FWHM. It can be pointed out here that since the nondiffracting propagation range is large, for a sample having a length much smaller than ${Z}_{\mathrm{PR}}$ value, a careful positioning of the sample in the beam path would not be required for unchanged spot-size and intensity values on a sample during continuous variation in spot size over a certain range. This may be useful in several applications.

Further, the nondiffracting propagation range ${Z}_{\mathrm{PR}}$ can be connected to FWHM as ${Z}_{\mathrm{PR}}\sim 2w/\mathrm{sin}\text{\hspace{0.17em}}\theta \sim 2kw/\alpha \sim 2kw(\mathrm{FWHM})/2.27$ [as $\mathrm{FWHM}=2.27/\alpha $]. Thus the propagation-invariant range is larger for a larger FWHM spot size of the Bessel beam. This is reflected from our data in Fig. 6 where the nondiffracting range is reduced with a decrease in FWHM [see graphs in Figs. 6(a)–6(d)].

As can be noted from Fig. 6, the FWHM of the generated beam (filled circles) also varies with distance $z$. These observations are in agreement with our simulation results presented afterwards in this paper. Since the variation in FWHM spot size, as well as peak intensity, is slow with $z$, the beam seems to evidently preserve its nondiffracting property.

#### C. Effect of HB Ring-Width on Bessel-Beam Generation

It is known that a Fourier transform of a ring function is a ${J}_{0}$-Bessel function, where the radius of the ring determines the transverse frequency of the ${J}_{0}$-Bessel function variation. Hence, focusing of an annular intensity profile with a very small ring width should result in a Bessel-beam field, and Bessel-beam FWHM spot size should vary with hollow-beam diameter as we experimentally observed (Fig. 5). These conditions were satisfied in our experiments described above (with $w=0.57\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$), and a good quality Bessel beam was obtained as shown in Fig. 3. In order to investigate the effect of $w$ on the quality of a beam generated in the focal region of the lens, we measured the transverse intensity profile of a generated beam for different values of ring-width parameter $w$ of the HB. In these measurements, the input beam was expanded to increase the ring width of a HB at the lens plane for a fixed value of mirror separation ($d=90\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$). When input beam spot size was changed, we noted that, in addition to change in $w$, the peak intensity radius ${r}_{0}$ of HB was also changed slightly. This may be due to diffraction effects involved in the propagation of a ring from source position to the lens plane. The measured profiles of a generated Bessel beam are shown in Fig. 7 for different values of ring-width parameter $w$ of a HB. As can be noted from this figure, the profile of a generated beam gets modified in a qualitative way with an increase in a ring-width parameter $w$. The higher central-peak intensity (with relatively weaker side peaks) for larger values of $w$ shows that a generated beam was degraded and became much different from a ${J}_{0}$-Bessel beam at ring-width parameter $w=1.29\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm}$ [Fig. 7(c)]. We have also measured the variation in peak intensity and FWHM of the beam with axial distance $z$, for different values of ring-width parameter $w$. Results of these measurements are shown in Fig. 8. The results show reduction in nondiffracting propagation range ${Z}_{\mathrm{PR}}$ with increase in $w$. This is in agreement with the nature of spatial profiles shown in Fig. 7. These experimental results are also qualitatively in agreement with the results of simulations presented in next subsection.

Thus, contrary to the geometrically estimated behavior as discussed before that nondiffracting propagation range (${Z}_{\mathrm{PR}}\sim 2w/\mathrm{sin}\text{\hspace{0.17em}}\theta $) should scale linearly with ring width $w$, we observed reduction in the nondiffracting propagation range with increase in $w$. This was due to deviation of generated beam profile from a ${J}_{0}$-Bessel beam profile for larger values of $w$.

#### D. Results of Simulations

Earlier in a work by Pu *et al.* [12], it was shown that a focusing of a Gaussian beam (with its central part blocked) by using an aberrated converging lens can result in zero-order Bessel function type transverse field distribution. To show theoretically for our setup that focusing a HB using a converging spherical lens can result in a field distribution close to that of a ${J}_{0}$-Bessel beam, we numerically solved the Fresnel integral for our geometry to find field distribution after the lens. We used formulism similar to that used in earlier works [12,26]. When a lens is illuminated by a HB (Fig. 9), we can write the field distribution at point $P$ ($\rho $, $z$) after the lens as

In the above definition of $t(r)$, $A(r)$ is the function describing the variation of the HB field with $r$ at lens plane. Here $\beta $ is the lens aberration coefficient and ${f}^{\prime}$ is paraxial focal length of the lens. For a uniform intensity HB, $A(r)=1$. For a HB having a spatial intensity profile of its ring, $A(r)$ is a function of position $r$. To model our observed hollow-beam intensity profile, we have assumed $A(r)$ to be of form

where ${r}_{0}$ is peak intensity radius of HB and $w$ (separation between peak position ${r}_{0}$ and the point having $1/{e}^{2}$ times the peak intensity) is the ring-width parameter of the HB. In our simulations, to obtain field distribution after the lens by evaluating the integral of Eq. (1), we have set parameters ${r}_{0}$ and $w$ by comparing Eq. (3) with the experimentally measured intensity profile of the HB at the lens plane.Figure 10 shows the transverse intensity profiles, obtained after evaluating Eq. (1), for different values of axial position $z$ for two different values of ring-width parameter ($w$). These profiles were evaluated for different hollow-beam parameters which were set after measuring the hollow-beam profiles for a fixed $d(=90\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{mm})$ but varying the expansion factor in the input beam (as done for taking data shown in Figs. 7 and 8 for different values of $w$). While changing the spot size of the input beam, we noted that, in addition to change in $w$, there was a small change in the hollow-beam peak radius ${r}_{0}$ also. This may be because of diffraction effects. We estimated the lens parameters $\beta $ and ${f}^{\prime}$ by measuring the focal positions of HBs of different radii and applying the formulism similar to given in [12].

As can be noted from Fig. 10, for larger value of $w$ [Fig. 10(b)] the calculated intensity profile shows more deviation from an ideal ${J}_{0}$-Bessel beam profile. These calculated profiles show the trend similar to the experimentally observed profiles shown in Fig. 7. Similarly, it can be seen that the peak intensity (Fig. 11) and FWHM spot size (Fig. 12) obtained from numerical simulations follow the similar trend in variation with axial distance $z$ as observed experimentally for different values of $w$ (Fig. 8). Thus our simulations show a qualitative agreement with the experimentally observed results. From these findings, it is clear that smaller ring width of the HB is required for generation of a good-quality Bessel beam.

## 4. Conclusion

Using a simple setup we have demonstrated the generation of a zero-order Bessel beam of continuously variable central spot size. This was achieved by focusing a variable-diameter collimated HB which was generated after passing a laser beam through a setup consisting of a pair of convex and concave metal axicon mirrors. In the experimental setup, by varying the mirror separation, we could obtain the variation in the spot size of the generated Bessel beam. Such a variable spot-size ${J}_{0}$-Bessel beam may find various potential applications in atom guiding, optical tweezers, and laser machining. We note that range over which Bessel-beam spot size can be varied in our setup remains limited due to limits on separation between mirrors and size of the optics used. This range can be extended by using a lens of different focal length.

We have also investigated the effect of hollow-beam ring width on spatial profile and non-diffracting propagation range of Bessel beams generated in our setup. Our results show that smaller ring width of the HB gives the better quality of generated Bessel beams.

We thank R. Balasubramaniam (Machine Dynamics Division, BARC, Mumbai), C. Mukharjee (MOSS, RRCAT, Indore), and L. Sudarsanam (Laser Workshop, RRCAT, Indore) for helping in fabrication of axicon mirrors; and H. S. Vora (LESD, RRCAT, Indore) for providing image processing software. We acknowledge the contribution of Optical Workshop, RRCAT, Indore for lapping the copper disks used in making axicon mirrors.

## References

**1. **J. Durnin, J. J. Miceli Jr., and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. **58**, 1499–1501 (1987). [CrossRef]

**2. **S. Schmid, G. Thalhammer, K. Winkler, F. Lang, and J. H. Denschlag, “Long distance transport of ultracold atoms using a 1D optical lattice,” New J. Phys. **8**, 159 (2006). [CrossRef]

**3. **J. Arlt and K. Dholakia, “Generation of high-order Bessel beams by use of an axicon,” Opt. Commun. **177**, 297–301 (2000). [CrossRef]

**4. **J. Arlt, V. Garcés-Chávez, W. Sibbett, and K. Dholakia, “Optical micromanipulation using a Bessel light beam,” Opt. Commun. **197**, 239–245 (2001). [CrossRef]

**5. **M. Yoshihiko and H. Makoto, “Micro grooving of metallic material using a Bessel beam,” Rev. Laser Eng. **34**, 842–847 (2006).

**6. **Y. Matsuoka, Y. Kizuka, and T. Inoue, “The characteristics of laser micro drilling using a Bessel beam,” Appl. Phys. A **84**, 423–430 (2006). [CrossRef]

**7. **M. Riox, R. Tremblay, and P. A. Belanger, “Linear, annular and radial focusing with axicons and applications to laser machining,” Appl. Opt. **17**, 1532–1536 (1978). [CrossRef]

**8. **X. Tsampoula, V. Garcés-Chávez, M. Comrie, D. J. Stevenson, B. Agate, C. T. A. Brown, F. Gunn-Moore, and K. Dholakia, “Femtosecond cellular transfection using a nondiffracting light beam,” Appl. Phys. Lett. **91**, 053902 (2007). [CrossRef]

**9. **S. R. Mishra, “A vector wave analysis of a Bessel beam,” Opt. Commun. **85**, 159–161 (1991). [CrossRef]

**10. **J. Turunen, A. Vasara, and A. T. Friberg, “Realization of general nondiffracting beams with computer-generated holograms,” J. Opt. Soc. Am. A **6**, 1748–1754 (1989). [CrossRef]

**11. **R. M. Herman and T. A. Wiggins, “Production and uses of diffractionless beams,” J. Opt. Soc. Am. A **8**, 932–942 (1991). [CrossRef]

**12. **J. Pu, H. Jhang, and S. Nemoto, “Lens axicons illuminated by Gaussian beams for generation of uniform-axial intensity Bessel fields,” Opt. Eng. **39**, 803–807 (2000). [CrossRef]

**13. **J. F. Fortin, G. Rousseau, N. McCarthy, and M. Piche, “Generation of quasi-Bessel beams and femtosecond optical $X$-waves with conical mirrors,” Proc. SPIE **4833**, 876–884 (2003). [CrossRef]

**14. **S. R. Mishra, S. K. Tiwari, S. P. Ram, and S. C. Mehendale, “Generation of hollow conic beams using a metal axicon mirror,” Opt. Eng. **46**, 084002 (2007). [CrossRef]

**15. **K. B. Kuntz, B. Braverman, S. H. Youn, M. Lobino, E. M. Pessina, and A. I. Lvovsky, “Spatial and temporal characterization of a Bessel beam produced using a conical mirror,” Phys. Rev. A **79**, 043802 (2009). [CrossRef]

**16. **S. K. Tiwari, S. P. Ram, J. Jayabalan, and S. R. Mishra, “Measuring a narrow Bessel beam spot by scanning charge-coupled device (CCD) pixel,” Meas. Sci. Technol. **21**, 025308 (2010). [CrossRef]

**17. **E. E. Ushakova and S. N. Kurilkina, “Formation of Bessel light pulses by means of a conical mirror,” J. Appl. Spectrosc. **77**, 827–831 (2011). [CrossRef]

**18. **S. K. Tiwari, S. R. Mishra, and S. P. Ram, “Generation of a variable diameter collimated hollow laser beam using metal axicon mirrors,” Opt. Eng. **50**, 014001 (2011). [CrossRef]

**19. **N. Chattrapiban, E. A. Rogers, D. Cofield, W. T. Hill III, and R. Roy, “Generation of nondiffracting Bessel beams by use of a spatial light modulator,” Opt. Lett. **28**, 2183–2185 (2003). [CrossRef]

**20. **E. McLeod, A. B. Hopkins, and C. B. Arnold, “Multiscale Bessel beams generated by a tunable acoustic gradient index of refraction lens,” Opt. Lett. **31**, 3155–3157 (2006). [CrossRef]

**21. **E. McLeod and C. B. Arnold, “Optical analysis of time-averaged multiscale Bessel beams generated by a tunable acoustic gradient index of refraction lens,” Appl. Opt. **47**, 3609–3618 (2008). [CrossRef]

**22. **G. Milne, G. D. Jeffries, and D. T. Chiu, “Tunable generation of Bessel beams with a fluidic axicon,” Appl. Phys. Lett. **92**, 261101 (2008). [CrossRef]

**23. **D. Brousseau, J. Drapeau, M. Piche, and E. F. Borra, “Generation of Bessel beams using a magnetic liquid deformable mirror,” Appl. Opt. **50**, 4005–4010 (2011). [CrossRef]

**24. **J. Brunne, M. Bock, A. Treffer, U. Wallrabe, and R. Grunwald, “Adaptive generation of Bessel-like beams by reflective multi-electrode piezo-axicons,” in Conference on Lasers and Electro-Optics—European Quantum Electronics Conference, OSA Technical Digest Series (Optical Society of America, 2011), paper CF_P14.

**25. **V. Vaičaitis and Š. Paulikas, “Formation of Bessel beams with continuously variable cone angle,” Opt. Quantum Electron. **35**, 1065–1071 (2003). [CrossRef]

**26. **Z. Jaroszewicz and J. Morales, “Lens axicons: systems composed of a diverging aberrated lens and a perfect converging lens,” J. Opt. Soc. Am. A **15**, 2383–2390 (1998). [CrossRef]