OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 20, Iss. 24 — Nov. 19, 2012
  • pp: 26852–26867
« Show journal navigation

Generalizing higher-order Bessel-Gauss beams: analytical description and demonstration

Damian N. Schimpf, Jan Schulte, William P. Putnam, and Franz X. Kärtner  »View Author Affiliations


Optics Express, Vol. 20, Issue 24, pp. 26852-26867 (2012)
http://dx.doi.org/10.1364/OE.20.026852


View Full Text Article

Acrobat PDF (2932 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We report on a novel class of higher-order Bessel-Gauss beams in which the well-known Bessel-Gauss beam is the fundamental mode and the azimuthally symmetric Laguerre-Gaussian beams are special cases. We find these higher-order Bessel-Gauss beams by superimposing decentered Hermite-Gaussian beams. We show analytically and experimentally that these higher-order Bessel-Gauss beams resemble higher-order eigenmodes of optical resonators consisting of aspheric mirrors. This work is relevant for the many applications of Bessel-Gauss beams in particular the more recently proposed high-intensity Bessel-Gauss enhancement cavities for strong-field physics applications.

© 2012 OSA

1. Introduction

Bessel-Gauss beams offer properties that Gaussian beams do not possess. For instance, they exhibit an extended depth of field, and away from a strong on-axis intensity region they exhibit an annular intensity distribution in the far field [1

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

, 2

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

]. These attributes have made these beams interesting for applications in biomedical imaging, material processing, particle trapping, or laser-based acceleration of particles [3

3. F. O. Fahrbach, P. Simon, and A. Rohrbach, “Microscopy with self-reconstructing beams,” Nat. Photonics 4, 780–785 (2010). [CrossRef]

7

7. D. Li and K. Imasaki, “Vacuum laser-driven acceleration by two slits-truncated Bessel beams,” Appl. Phys. Lett. 87, 091106 (2005). [CrossRef]

].

Laser resonators supporting Bessel-Gauss beams have also been designed and implemented [8

8. J. C. Gutiérrez-Vega, R. Rodríguez-Masegosa, and S. Chávez-Cerda, “Bessel-Gauss resonator with spherical output mirror: geometrical- and wave-optics analysis,” J. Opt. Soc. Am. A 20, 2113–2122 (2003). [CrossRef]

10

10. B. Ma, F. Wu, W. Lu, and J. Pu, “Nanosecond zero-order pulsed Bessel beam generated from unstable resonator based on an axicon,” Opt. Laser Technol. 42, 941–944 (2010). [CrossRef]

]. More recently, enhancement cavities based on Bessel-Gauss-type modes have been proposed for strong-field applications, such as cavity-enhanced high-harmonic generation (HHG). Bessel-Gauss modes offer the benefits of near-perfect out-coupling of the generated intra-cavity high-harmonics and an increased intensity ratio from the focus to the cavity mirror surfaces [11

11. 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]

]. In these applications the well-known Bessel-Gauss beam has been exclusively considered as the operating mode. Having observed higher-order azimuthally symmetric mode solutions with numerical mode-solvers, the most pertinent questions we ask are: Do higher-order Bessel-Gauss beams exist? How can they be described analytically and demonstrated experimentally?

The existence of higher-order Bessel-Gauss beams showing azimuthal phase dependence (i.e. carry optical angular momentum) is known [12

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

, 13

13. R. Vasilyeu, A. Dudley, N. Khilo, and A. Forbes, “Generating superpositions of higherorder Bessel beams,” Opt. Express 17, 23389–23395 (2009). [CrossRef]

]. However, azimuthally symmetric higher-order Bessel-Gauss beams are hitherto not reported in the literature, to the best of our knowledge.

In the present paper, we first analytically describe higher-order Bessel-Gauss beams. Based on these findings, we then design an optical resonator and show numerically that it supports modes resembling azimuthally symmetric higher-order Bessel-Gauss beams. Finally, we demonstrate these higher-order Bessel-Gauss modes experimentally.

2. Generalizing higher-order Bessel-Gauss beams

Conventional Bessel-Gauss beams can be built up by a superposition of decentered Gaussian beams, whose centers are positioned on a circle and whose beam directions are pointing to the apex of a cone [14

14. C. Palma, “Decentered Gaussian beams, ray bundles, and Bessel-Gauss beams,” Appl. Opt. 36, 1116–1120 (1997). [CrossRef] [PubMed]

, 15

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

]. In the following we will show that higher-order Bessel-Gauss beams can be formed in a similar way by using the more general decentered Hermite-Gaussian beams as component beams for the superposition on a single cone.

First, we define decentered Hermite-Gaussian beams and describe their propagation through optical systems. After that we superpose these component beams to the higher-order Bessel-Gauss beams and also discuss their propagation properties. Finally, we show that these beams can collapse to the azimuthally symmetric Laguerre-Gaussian beams.

2.1. Decentered Hermite-Gaussian beams

Fig. 1 Geometry for the superposition of generalized decentered Hermite-Gaussian beams.

2.2. Propagation of decentered Hermite-Gaussian beams through optical systems

2.3. Superposing decentered Hermite-Gaussian beams to higher-order Bessel-Gauss beams

We form the higher-order Bessel-Gauss beams by superimposing the component beams with the following attributes: the Hermite-Gaussian component is modulated only in the radial direction. The vectors rd0 and ε0 are collinear with rd0= const. and ε0= const. This corresponds to a superposition of beams whose mean centers are placed on a circle with radius rd0 in the z=0 plane and whose mean beam directions point to the apex of a single cone with a semi-aperture angle ε0, as shown in Fig. 1.

For the superposition, we express each component beam in an auxilliary coordinate system, which is rotated by an angle γ with regard to a fixed reference coordinate system (r,θ,z) so that its abscissa points in the direction of the vector rd0. In this auxiliary system, the decentered Hermite-Gaussian beam consists solely of the Hm-term (i.e. Hn=0 = 1) and the abscissa coordinate is given by r·cos(θγ). The propagated decentered Hermite-Gaussian beam is expressed in the cylindrical coordinates (r,θ,z) as
vγ(r,θ,L)=1wexp(ikL)exp(iϕ(m+1))Hm(2(rcos(θγ)rd)w)×exp(ik2q(r2+rd22rrdcos(θγ)))exp(ikεrcos(θγ))exp(iφ).
(7)
Because of the parallelism of rd0 and ε0, the following beam-propagation transformations hold
rd=xd2+yd2=Ard0+Bε0
(8)
ε=εx2+εy2=Crd0+Dε0
(9)
φ=k2[Crd0rd+Bε0ε]
(10)
in addition to the relations for the Gaussian beam parameters, Eqs. (5) and (6). The formulation of higher-order Bessel-Gauss beams in terms of ABCD-matrix parameters allows us to study their transformation under the impact of optical components. We will expand upon this feature in section 3.1 describing optical resonators for these beam solutions.

The generalization of the Bessel-Gauss beam solution is found by integrating over the γ angle from 0 to 2π:
uml(r,θ,L)=2πwexp(ikL)exp(iϕ(m+1))exp(iφ)exp[ik2q(r2+rd2)]ml.
(11)
The term ℐml is given by
ml=12π02πdγHm(acos(θγ)b)exp(iαcos(θγ))exp(ilγ),
(12)
where we substitute
α=kr(ε(rd/q));a=2r/w;b=2rd/w.
(13)
To generalize the analysis, we have introduced an azimuthal phase-variation while superposing the component beams. This phase is characterized by the index l. The integral ℐml can be solved analytically (see Appendix), and the first four expressions are given by
0l=eilθilJl(α)
(14)
1l=eilθ2il[2ia((1)lJ1l(α)+J1+l(α))4bJl(α)]
(15)
2l=eilθ2il[2a2((1)lJ2l(α)+J2+l(α))i8ab((1)lJ1l(α)+J1+l(α))
(16)
+(4a2+8b24)Jl(α)]
(17)
3l=eilθ2il[i2a3((1)lJ3l(α)+J3+l(α))+12a2b((1)lJ2l(α)J2+l(α))+i(6a3+24ab212a)((1)lJ1l(α)+J1+l(α))+2(8b312a2b+12b)Jl(α)]
(18)
If m = 0 and l = 0, the well-known Bessel-Gauss beam is obtained [14

14. C. Palma, “Decentered Gaussian beams, ray bundles, and Bessel-Gauss beams,” Appl. Opt. 36, 1116–1120 (1997). [CrossRef] [PubMed]

]. Furthermore, the solutions comprising the terms ℐ0l (i.e. m = 0) are identical to previously defined higher-order Bessel-Gauss beams exhibiting an azimuthal phase-variation [12

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

, 13

13. R. Vasilyeu, A. Dudley, N. Khilo, and A. Forbes, “Generating superpositions of higherorder Bessel beams,” Opt. Express 17, 23389–23395 (2009). [CrossRef]

].

Here, we introduce higher-order Bessel-Gauss beams with a radial index m > 0. In Fig. 2 we illustrate the effect of the radial index by plotting the azimuthally symmetric (l = 0) higher-order Bessel-Gauss beams for different orders m in the near and far-field. The parameters are w0 = 200μm, rd0 = 0, ε = 0.3°, and λ = 1040nm. In the far-field (here, calculated at z= 4 · zR, where zR=kw02/2), all beams show an annular beam with variations in the radial direction such that the number of nodes corresponds to m. If l > 0 these solutions exhibit azimuthal phase-variations, as displayed in Fig. 3. To highlight the impact of the azimuthal index l on the solution, we consider the superposition (um,l + um,−l) resulting in a cosine-variation over the angular coordinate.

Fig. 2 Absolute value of the azimuthally symmetric Bessel-Gauss beam at z=0 for l=0 and (a): m=0, (b): m=1, (c): m=2, and in the far-field (z= 4 · zR); (d): m=0, (e): m=1, (f): m=2.
Fig. 3 Absolute value of the superposition of higher-order Bessel-Gauss beams in the far-field, complementary to Fig. 2 (at z= 4 · zR) (a): (u0,1 + u0,−1); (b): (u1,1 + u1,−1); (c): (u2,1 + u2,−1); (d): (u0,2 + u0,−2); (e): (u1,2 + u1,−2); (f): (u2,2 + u2,−2).

How does the Bessel-Gauss beam solution depend on the parameters rd0 and ε? - A classification is known for the fundamental Bessel-Gauss beam [12

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

]: Generalized Bessel-Gauss beam, ordinary Bessel-Gauss beam and modified Bessel-Gauss beam denote the cases (rd0 ≠ 0, ε0 ≠ 0); (rd0 = 0, ε0 ≠ 0), and (rd0 ≠ 0, ε0 = 0), respectively. We extend this terminology to the higher-order Bessel-Gauss beams. For the example of a Bessel-Gauss beam with m=2, we plot the different cases in Fig. 4. On the left side, the radial profile is plotted over the propagation distance. On the right side, we display the transverse beam pattern at the position of minimum Gaussian waist (here at z=0). Figures 4(a) and 4(b) show the absolute value of the amplitude for a generalized higher-order (m=2) Bessel-Gauss beam for w0 = 200μm, rd0 = 0.8mm, ε = 0.3°, and λ = 1040nm. The vertical blue line highlights the (z=0) position. For this beam type the intersection with the symmetry axis is shifted from the position where the minimum Gaussian waist occurs. Figures 4(c) and 4(d) displays the ordinary beam type. The intersection coincidences with the position of minimum waist. Figures 4(e) and 4(f) picture the modified beam type. Furthermore, if both rd0 and ε0 are zero, a azimuthally symmetric Laguerre-Gaussian beam is recovered, as displayed in Fig. 4(g) and 4(h).

Fig. 4 Beam propagation and transverse pattern (at z=0) for: (a) and (b) generalized higher-order (m=2, l=0) Bessel-Gauss beam; (c) and (d) ordinary higher-order (m=2, l=0) Bessel-Gauss beam; (e) and (f) modified higher-order (m=2, l=0) Bessel-Gauss beam; (g) and (h) azimuthally symmetric Laguerre-Gaussian beam of order p=1, respectively.

2.4. Collapse of the solution to Laguerre-Gaussian beams for ε = 0 and rd0 = 0

To show the collapse of the azimuthally symmetric higher-order Bessel-Gauss beams to Laguerre-Gaussian beams, we write the integral in Eq. (12) for ε = 0 and rd0 = 0. For uneven m the integral vanishes due to symmetry, and for even m (= 2 · p) we find:
12π02πdγH2p(2(r/w)cos(θγ))=(1)p(2p)!p!Lp(2r2/w2),
(19)
This integration is performed by substituting t = cos(θγ) and applying [18

18. G. Ryshik, Tables of Series, Products and Integrals (Verlag Harri Deutsch, 1981).

]:
11dt(1t2)α12H2p(ut)=(1)pπ(2p)!Γ(α+12)Γ(p+α+1)Lpα(u)
(20)
which holds for Im[α] > −1/2.

Equation (20) describes the remarkable property that the azimuthally symmetric Laguerre-Gaussian beams are a subset of the presented novel beam class. Thus, the presented solutions are a more general description for azimuthally symmetric beams.

3. Comparing the Bessel-Gauss beam solutions to modes of optical resonators

In the following, we would like to compare the analytical expressions for the fundamental and azimuthally symmetric higher-order Bessel-Gauss beams with numerical results. Specifically, we simulate a resonant optical cavity that supports a conventional generalized Bessel-Gauss beam for the fundamental mode, and should exhibit higher-order Bessel-Gauss beams as higher-order modes.

3.1. Design of the optical resonator

Figure 5 shows a schematic of the optical resonator, it consists of a spherical mirror with a radius of curvature R and an axicon mirror with a base angle ε.

Fig. 5 Schematic of the optical resonator.

Ray matrices can be employed for the design of the optical resonator. One roundtrip is described by
(ABCD)=(1L01)(102/R1)(1L01).
(21)
Geometrical optics determines the radius of the displacement and tilt of Bessel-Gauss beams during propagation, as expressed by Eqs. (8) and (9). At the axicon mirror, these parameters are described by (rd0, ε) and (rd0, −ε) before and after one roundtrip, respectively. Thus, the resonator length is given by
L=Rrd0/ε.
(22)
Gaussian beam optics determines the Gaussian beam parameter q, which must repeat after one round trip. At the axicon mirror, one finds
q=iLRL2.
(23)
The resonator is stable for L < R. The minimum waist for the Gaussian component of the Bessel-Gauss beam occurs at the axicon mirror and can be determined from q = −izR.

To summarize, the resonator length and the waist w0 of the Gaussian component of the Bessel-Gauss beam solution are fully governed by the radius of curvature of the spherical resonator mirror, the base angle of the axicon mirror, as well as the radius of the ring of the generalized Bessel-Gauss beam on the axicon, rd0. For the following calculations, we chose a radius of curvature R = 250 mm, ε = 0.5° and rd0= 1.5 mm, which result in L ≈ 78 mm and w0 ≈ 196μm for a light wavelength of 1040 nm. These particular parameters were selected since they approximate the experimental configuration which, in turn, was determined by the availability of base angles for the axicon mirror and radii of curvatures for the spherical mirror. It results in a ring radius of the Bessel-Gauss beam that avoids the round tip on the axicon mirror and results in a not too large ring radius on the curved mirror. This situation is advantageous as it reduces effects of surface variations on the resonator mode.

The spherical mirror serves as the coupler to the resonator. For the coupling to the resonator with an excitation beam, it is important to know the beam parameters at the plane backside of this mirror. Geometrical optics of the meridional ray, characterizing the propagation of the Bessel-Gauss mode, gives a semi-aperture angle ε′ ≈ nm · ε, where nm is the refractive index of the mirror substrate, and a ring radius rdrd + Δ · ε, where the ring-radius at the front surface is given by rd = ε · R and Δ is the geometrical thickness of the mirror, which is given by Δ = Δouts(rout), where Δout is the mirror thickness at the outer radius rout and the surface function is given by s(r)=RR2r2.

3.2. Numerical mode solver

Diffraction can be described in terms of a Hankel transform, which we numerically implement by employing a quasi-discrete Hankel transform (i.e. a Fourier-Bessel series expansion) [19

19. H. F. Johnson, “An improved method for computing a discrete hankel transform,” Comp. Phys. Comm. 43, 181–202 (1987). [CrossRef]

21

21. M. Guizar-Sicairos and J. C. Gutiérrez-Vega, “Computation of quasi-discrete Hankel transforms of integer order for propagating optical wave fields,” J. Opt. Soc. Am. A 21, 53–58 (2004). [CrossRef]

]. In the following, we solve for the azimuthally symmetric modes.

3.3. Comparison of the numerical mode-solution to the Bessel-Gauss beam solutions

To demonstrate the accuracy of the analytical expressions, the Bessel-Gauss beam solutions of different radial orders m are evaluated and compared to numerical results at the axicon mirror, and at two z-positions outside the cavity. The mode field is propagated in free space from the axicon mirror to the on-axis crossing point at a distance −(RL) via an intermediate position −(0.9 · RL). In Fig. 5, these positions are denoted as P1, P2 and P3. Figure 6 shows a comparison between numerical results and the analytical solutions. It can be seen that the azimuthally symmetric higher-order Bessel-Gauss modes represent an accurate analytical model for higher-order modes in resonators supporting the (m=0, l=0) Bessel-Gauss beam as the fundamental mode. Conversely, the numerical results validate the existence of azimuthally symmetric higher-order Bessel-Gauss beams, and constitute a reference for the hitherto unknown analytical beam solution.

Fig. 6 (a–c) Absolute value of the amplitude of the fundamental (m=0) Bessel-Gauss beam, (d–f) m=1 higher-order Bessel-Gauss mode, and (d–f) m=2 higher-order Bessel-Gauss mode at distance 0, −(0.9 · RL), and −(RL) from the axicon.

Fig. 7 (a) Absolute value of the amplitude of the generalized Bessel-Gauss beam vs distance from the minimum Gaussian waist point, (b) and (c): patterns for the (m=1) and (m=2) azimuthally symmetric higher-order Bessel-Gauss beams, respectively.

3.4. Experimental setup

To observe the higher-order Bessel-Gauss beams experimentally, we have designed an enhancement cavity according to the guidelines, which have been presented in section 3.1. As shown in Fig. 5, the resonator is formed by an axicon mirror (AM) and a spherical mirror (SM). In the experiment, the conic side of the axicon mirror has a base angle of 0.5° and the coating has a reflectance of 0.99. The spherical mirror has a radius of curvature of R = 250 mm and a reflectance of 0.9. The 1-inch diameter spherical mirror (with a 1/4-inch thickness of the BK-7 substrate) serves as the input coupler to the resonator. The back sides of both mirrors are anti-reflection coated. In the design configuration, the separation distance between the two mirrors is L = 75 mm, and the ring radius of the Bessel-Gauss mode on the axicon is approximately 1.5 mm.

A schematic of the excitation path to the resonator is shown in Fig. 8. We use a tunable grating-stabilized external cavity diode laser, which is coupled to a single-mode fiber (PM980-XP, mode-field diameter of approximately 6.6 μm). The laser is operated at a wavelength of 1040 nm and it emits up to 30 mW of power at the fiber output with a line width of 100 kHz. The laser frequency is varied by tuning the grating angle of the external cavity with a piezoelectric element. This allows us to lock the laser frequency to the resonance of the enhancement cavity by dither locking. For this, the light that is reflected from the in-coupling mirror SM is detected at the side port of an optical isolator with a photodiode and is used as a feedback signal for the locking loop. The leakage through the axicon mirror is used to record the resonator modes with a camera, which is placed behind AM.

Fig. 8 Schematic of the experimental setup.

To match the laser beam to the resonator mode, we first magnify the beam emanating from the fiber by a factor of 450 using the lenses L1, L2, and L3 with focal lengths of 4 mm, 10 mm, 100 mm, respectively. An axicon (A) generates an approximate Bessel-Gauss beam with a semi-aperture angle of (n − 1) · δ = 0.45°, where n is the refractive index of fused silica (n=1.45) and δ the base angle of the axicon, δ = 1°. The waist of the beam illuminating the axicon is 1.5 mm. Due to the round, and thus imperfect, apex of the axicon, the resulting beam in the far field is not a pure ring-like beam but still contains unwanted light in the central region. The quality of the beam is increased by employing Fourier filtering with an opaque disk of 1 mm radius in a 4-f system (150 mm focal-length lenses L4 and L5) [24

24. 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]

]. Behind the 4-f system, a 300 mm focal-length lens (L6) and a 200 mm focal-length lens (L7) match the excitation beam to the beam parameters required at the plane side of SM.

To determine the positions of the lenses L6 and L7, at first, we consider a configuration that assumes an ideal Bessel-Gauss beam being produced by the axicon. After the 4-f system the excitation beam is collimated by the lens L6 and is then focused into the resonator by the lens L7. This results in a semi-aperture angle of 3/2· 0.45° = 0.68°. The ring radius of the beam on SM determines the distance between L7 and SM. To better match the beam to the cavity mode, the excitation path is simulated numerically. We account for an imperfect axicon showing a parabolic shape that deviates by 8 μm from the ideal conic shape at the tip. As a consequence, we iteratively improve the distances between L5 and L6, L6 and L7, as well as L7 and SM for best excitation. The resulting values are given in parentheses in Fig. 8, and are used as a starting point for the experiment, in which we also adjust the resonator mirrors and set the final positions of L6 and L7 so that the photodiode signal shows the most pronounced resonance feature.

3.5. Experimental results: comparison of the experimental modes to the beam solutions

Figure 9 shows the photodiode signal during a scan of the laser frequency. The resonances of the higher-order modes are much less pronounced compared to the one of the fundamental Bessel-Gauss mode. This is due to the lower excitation efficiency of the higher-order beam with the beam coupled into the cavity. However, the resonance can be still recognized in the photodiode signal and images of the resonator modes can be taken, as shown in Fig. 10(a)–(d). To obtain these pictures, we limit the range of the frequency scan to the resonance region of the mode of interest. The images are recorded at a scan rate of the piezo that is fast compared to the integration time of the camera. In this way, the camera averages several scans. These images were taken at a distance z= 25 mm behind the flat surface of the axicon mirror. The evolution of the beam profile as a function of other distances from the axicon mirror will be similar to the one shown in Fig. 7. However, it must be noted that in Fig. 7 we have neglected refraction effects by the axicon-mirror substrate. Behind the axicon the semi-aperture angle of the higher-order Bessel-Gauss mode is changed from 0.5° to 0.725° due to refraction. Compared to Fig. 7, this will move the on-axis crossing-point of the beam closer to the axicon mirror. In Fig. 10 azimuthal variations of the intensity on the annulus can be seen, which are due to surface imperfections of the resonator mirrors [11

11. 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]

]. Faint central rings can be observed in the experimental data for the higher-order modes, Fig. 10 (b), 10(c) and 10(d). These artifacts are due to transmitted excitation light. The main reason for the show up are the weak resonances of the higher order modes compared to the fundamental mode, as shown in Fig. 9. Figures 10(e)–(h) shows numerical simulations corresponding to the experimental data (taking the refraction of the beam in the axicon mirror substrate into account).

Fig. 9 Signal of the photodiode. The labels correspond to the mode images in Fig. 10.
Fig. 10 (a)–(d): Experimentally obtained images of higher-order modes m= 0, 1, 2, 3, respectively. (e)–(h): corresponding numerical simulations.

4. Conclusion

In this work, Bessel-Gauss beams have been generalized to include a set of higher-order solutions that are analytically derived and experimentally demonstrated. Expressions for these beams have been obtained by superposing decentered Hermite-Gaussian beams. The well-known Bessel-Gauss beam is the fundamental solution of this class. Moreover, we have rigorously demonstrated that the azimuthally symmetric Laguerre-Gaussian beams are special cases of these higher-order Bessel-Gauss beams. Thus, these formerly distinct beam types i.e. Bessel-Gauss beams and Laguerre-Gauss beams have been unified. Furthermore, we have shown that the generalized higher-order Bessel-Gauss beams are present as higher-order modes in optical resonators consisting of aspheric mirrors.

Enhancement cavities based on Bessel-Gauss beams are seen to be an increasingly important way to employ strong-field physics in practical applications. Higher-order Bessel-Gauss beams can be used to tailor the intensity distribution in such devices. Other fields of application for these beams are imaging or control and manipulation of particles.

Appendix

To solve the integral of Eq. (12), we write the Hermite-polynomials Hm in their explicit form. For this we use their recursion relation: H0(t) = 1, H1(t) = 2t and Hm+1(t) = 2tHm(t) − 2mHm−1(t). This results in an expansions in powers of cos-functions. The n-th power of the cos-functions can be written in terms of the harmonics of its argument by using De Moivre’s formula, e.g.:
cos2(β)=(1+cos(2β))/2cos3(β)=(3cos(β)+cos(3β))/4cos4(β)=(3+4cos(2β)+cos(4β))/8
Then, we make use of the following relation
eilθ2π02πdβeiαcosβcos(nβ)eilβ=eilθin+l2[(1)lJnl(α)+Jn+1(α)],
(24)
which can be easily shown by writing the cosine function in terms of exponentials and applying the definition of the Bessel-function, e.g. as given in [12

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

]:
Jk(t)=12πik02πdβexp{i(tcosβkβ)}.
(25)
Moreover, we use the relation Jk(t) = (−1)k · Jk(t).

Acknowledgments

This work is supported by AFOSR grant FA9550-10-1-0063 and the Center for Free-Electron Laser Science. William P. Putnam acknowledges support by a NSF Graduate Fellowship.

References and links

1.

C. J. R. Sheppard and T. Wilson, “Gaussian-beam theory of lenses with annular aperture,” IEEE 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.

F. O. Fahrbach, P. Simon, and A. Rohrbach, “Microscopy with self-reconstructing beams,” Nat. Photonics 4, 780–785 (2010). [CrossRef]

4.

T. A Planchon, L. Gao, D. E Milkie, M. W Davidson, J. A Galbraith, C. G Galbraith, and Eric Betzig, “Rapid three-dimensional isotropic imaging of living cells using Bessel beam plane illumination,” Nat. Methods 8, 417–423 (2011). [CrossRef] [PubMed]

5.

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

6.

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]

7.

D. Li and K. Imasaki, “Vacuum laser-driven acceleration by two slits-truncated Bessel beams,” Appl. Phys. Lett. 87, 091106 (2005). [CrossRef]

8.

J. C. Gutiérrez-Vega, R. Rodríguez-Masegosa, and S. Chávez-Cerda, “Bessel-Gauss resonator with spherical output mirror: geometrical- and wave-optics analysis,” J. Opt. Soc. Am. A 20, 2113–2122 (2003). [CrossRef]

9.

A. N. Khilo, E. G. Katranji, and A. A. Ryzhevich, “Axicon-based Bessel resonator: analytical description and experiment,” J. Opt. Soc. Am. A 18, 1986–1992 (2001). [CrossRef]

10.

B. Ma, F. Wu, W. Lu, and J. Pu, “Nanosecond zero-order pulsed Bessel beam generated from unstable resonator based on an axicon,” Opt. Laser Technol. 42, 941–944 (2010). [CrossRef]

11.

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]

12.

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

13.

R. Vasilyeu, A. Dudley, N. Khilo, and A. Forbes, “Generating superpositions of higherorder Bessel beams,” Opt. Express 17, 23389–23395 (2009). [CrossRef]

14.

C. Palma, “Decentered Gaussian beams, ray bundles, and Bessel-Gauss beams,” Appl. Opt. 36, 1116–1120 (1997). [CrossRef] [PubMed]

15.

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

16.

A. R. Al-Rashed, “Spatial and temporal modes of resonators with dispersive phase-conjugate mirrors,” PhD Thesis (1997).

17.

S. A. Collins Jr., “Lens-System Diffraction Integral Written in Terms of Matrix Optics,” J. Opt. Soc. Am. 60, 1168–1177 (1970). [CrossRef]

18.

G. Ryshik, Tables of Series, Products and Integrals (Verlag Harri Deutsch, 1981).

19.

H. F. Johnson, “An improved method for computing a discrete hankel transform,” Comp. Phys. Comm. 43, 181–202 (1987). [CrossRef]

20.

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

21.

M. Guizar-Sicairos and J. C. Gutiérrez-Vega, “Computation of quasi-discrete Hankel transforms of integer order for propagating optical wave fields,” J. Opt. Soc. Am. A 21, 53–58 (2004). [CrossRef]

22.

A. Fox and T. Li, “Computation of optical resonator modes by the method of resonance excitation,” IEEE J. Quantum Electron. 4, 460–465 (1968). [CrossRef]

23.

A. E. Siegman, Lasers (University Science Books, 1986).

24.

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]

OCIS Codes
(070.2580) Fourier optics and signal processing : Paraxial wave optics
(140.3300) Lasers and laser optics : Laser beam shaping
(140.3410) Lasers and laser optics : Laser resonators

ToC Category:
Physical Optics

History
Original Manuscript: October 1, 2012
Revised Manuscript: November 4, 2012
Manuscript Accepted: November 7, 2012
Published: November 14, 2012

Citation
Damian N. Schimpf, Jan Schulte, William P. Putnam, and Franz X. Kärtner, "Generalizing higher-order Bessel-Gauss beams: analytical description and demonstration," Opt. Express 20, 26852-26867 (2012)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-20-24-26852


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. C. J. R. Sheppard and T. Wilson, “Gaussian-beam theory of lenses with annular aperture,” IEEE 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. F. O. Fahrbach, P. Simon, and A. Rohrbach, “Microscopy with self-reconstructing beams,” Nat. Photonics4, 780–785 (2010). [CrossRef]
  4. T. A Planchon, L. Gao, D. E Milkie, M. W Davidson, J. A Galbraith, C. G Galbraith, and Eric Betzig, “Rapid three-dimensional isotropic imaging of living cells using Bessel beam plane illumination,” Nat. Methods8, 417–423 (2011). [CrossRef] [PubMed]
  5. M. Duocastella and C. B. Arnold, “Bessel and annular beams for materials processing,” Laser Photon. Rev.6, 607–621 (2012). [CrossRef]
  6. 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,” Nature419, 145–147 (2002). [CrossRef] [PubMed]
  7. D. Li and K. Imasaki, “Vacuum laser-driven acceleration by two slits-truncated Bessel beams,” Appl. Phys. Lett.87, 091106 (2005). [CrossRef]
  8. J. C. Gutiérrez-Vega, R. Rodríguez-Masegosa, and S. Chávez-Cerda, “Bessel-Gauss resonator with spherical output mirror: geometrical- and wave-optics analysis,” J. Opt. Soc. Am. A20, 2113–2122 (2003). [CrossRef]
  9. A. N. Khilo, E. G. Katranji, and A. A. Ryzhevich, “Axicon-based Bessel resonator: analytical description and experiment,” J. Opt. Soc. Am. A18, 1986–1992 (2001). [CrossRef]
  10. B. Ma, F. Wu, W. Lu, and J. Pu, “Nanosecond zero-order pulsed Bessel beam generated from unstable resonator based on an axicon,” Opt. Laser Technol.42, 941–944 (2010). [CrossRef]
  11. W. P. Putnam, D. N. Schimpf, G. Abram, and F. X. Kärtner, “Bessel-Gauss beam enhancement cavities for high-intensity applications,” Opt. Express20, 24429–24443 (2012). [CrossRef]
  12. V. Bagini, F. Frezza, M. Santarsiero, G. Schettini, and G. Spagnolo Schirripa, “Generalized Bessel-Gauss beams,” J. Mod. Opt.43(6), 1155–1166 (1996).
  13. R. Vasilyeu, A. Dudley, N. Khilo, and A. Forbes, “Generating superpositions of higherorder Bessel beams,” Opt. Express17, 23389–23395 (2009). [CrossRef]
  14. C. Palma, “Decentered Gaussian beams, ray bundles, and Bessel-Gauss beams,” Appl. Opt.36, 1116–1120 (1997). [CrossRef] [PubMed]
  15. A. R. Al-Rashed and B. E. A. Saleh, “Decentered Gaussian beams,” Appl. Opt.34, 6819–6825 (1995). [CrossRef] [PubMed]
  16. A. R. Al-Rashed, “Spatial and temporal modes of resonators with dispersive phase-conjugate mirrors,” PhD Thesis (1997).
  17. S. A. Collins, “Lens-System Diffraction Integral Written in Terms of Matrix Optics,” J. Opt. Soc. Am.60, 1168–1177 (1970). [CrossRef]
  18. G. Ryshik, Tables of Series, Products and Integrals (Verlag Harri Deutsch, 1981).
  19. H. F. Johnson, “An improved method for computing a discrete hankel transform,” Comp. Phys. Comm.43, 181–202 (1987). [CrossRef]
  20. L. Yu, M. Huang, M. Chen, W. Chen, W. Huang, and Z. Zhu, “Quasi-discrete Hankel transform,” Opt. Lett.23, 409–411 (1998). [CrossRef]
  21. M. Guizar-Sicairos and J. C. Gutiérrez-Vega, “Computation of quasi-discrete Hankel transforms of integer order for propagating optical wave fields,” J. Opt. Soc. Am. A21, 53–58 (2004). [CrossRef]
  22. A. Fox and T. Li, “Computation of optical resonator modes by the method of resonance excitation,” IEEE J. Quantum Electron.4, 460–465 (1968). [CrossRef]
  23. A. E. Siegman, Lasers (University Science Books, 1986).
  24. O. Brzobohaty, T. Cizmar, and P. Zemanek, “High quality quasi-Bessel beam generated by round-tip axicon,” Opt. Express16, 12688–12700 (2008). [CrossRef] [PubMed]

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.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited