Arxiv:1402.5106V1 [Physics.Optics] 20 Feb 2014 Vr Sntteol Osbedrcint Xlr Fone Beams

Home , Bessel beam

arXiv:1402.5106v1 [physics.optics] 20 Feb 2014 vr sntteol osbedrcint xlr fone beams. if optical explore of to classes direction new build possible to only wants the how- not This, is coordinate conditions. ever, boundary determinate suitable a with or paraxial in system the equation of solutions nonparaxial) represent (or all com- a they have re- feature: however, been examples, mon these also All have modified [19,20]. beams ported suitably Throcloidal as in known pattern cavities experimental field the On unconventional [18]. side, variables spatial so- the polynomial in beams and lutions [17] [16], spectra function) pendulum (Pearcey quantum with caustics form diffraction cusp the paraxial the of on in pattern the based of proposed beams propagation-invariant solutions were of of equations kind name nonparaxial to new and just Recently, beams, few. [15] a Laplace-Gauss Helmholtz- and [13], [10], [14] beams Gauss accelerating modes [11,12], Hypergeometric beams Airy [9], [6–8], beams beams generalized Ince-Gaussian [5], Laguerre-Gaussian beams and Bessel Hermite-Gaussian as scalar such the many regimes in however, vector both decade, and investigated past been have the solutions During other laser. emitted most field electromagnetic by the the represent of structure and transverse and respectively, Cartesian in coordinates are equation Cylindrical wave which paraxial [4], the Hermite- of beams solutions so-called Laguerre-Gaussian the and are Gaussian configurations plane the famous Perhaps, the propagates. most across distribution and intensity direction distributed this propagation how the be to spa- there- respect can their with (and transverse energy) intensity of the the study how the fore on namely to structure, dedicated amount tial was large fairly literature the a optical investigated, Among of an been [1–3]. of have al. freedom that of et beam degrees Maiman and the by properties of 1960s different development the very and in been invention Laser has the the form, either since study paraxial obey flourishing its the and or direction equation i.e., propagate preferred Helmholtz that a beams, along configurations optical mainly field of electromagnetic field of research The 3 nttt o pis nomto n htnc,Universi Photonics, and Information Optics, for Institute 2 a lnkIsiuefrteSineo ih,G Light, of Science the for Institute Planck Max 1 nttt fApidPyis rerc-cilrUnivers Friedrich-Schiller Physics, Applied of Institute pia oit fAmerica of Society Optical em.W netgt h omo hs e em swl sth as well as beams new consider zone these We paraxial of beam. the form in Bessel the has, investigate a We beams of beams. generalized spectrum such angular of the spectrum of form the erpr nanwcaso xc ouin ftesaa Helmh scalar the of solutions exact of class new a on report We CScodes: OCIS eeaie eslbaswt w indices two with beams Bessel Generalized 2021)EetoantcOtc;(6.02 igurO Singluar (260.6042) Optics; Electromagnetic (260.2110) ac Ornigotti Marco *[email protected] opldAgs ,2018 2, August Compiled u ¨ te-caosySrse1Bu4 15 ragn Germa Erlangen, 91058 1/Bau24, nther-Scharowsky-Strasse Germany 1 ∗ yo ragnNeneg tuttas /2 15 Erlan 91058 7/B2, Staudtstrasse Erlangen-Nuernberg, of ty n nraAiello Andrea and 1 t,Jn,MxWe lt ,073Jn,Germany Jena, 07743 1, Platz Max-Wien Jena, ity, where where sthe is h upsso hsLte,i scnein ocalculate distribution field to [21]: this convenient as to is associated it spectrum angular Letter, the this of purposes the pcrmrpeetto faBse em as beam, Bessel angular well-known a the of obtain representation we spectrum (2) into (1) Eq. tuting where h ria nua oetm(A)o h beam, the of (OAM) momentum R angular orbital the ino clrmncrmtcBse emcharacter- aperture beam angle cone Bessel a monochromatic by ized scalar a of tion the of maxima radial beam. the LG of corresponding position the by determined h oiin fterda aiao Gba.A a with As beams beam. Bessel apertures LG of angular a superposition different of a maxima obtains with general one radial function more result, the Delta a the of an- of with positions zeros the the beam the in Bessel connects function which a delta one of Dirac spectrum the replace gular To we beams. (LG) paraxial this, Laguerre-Gaussian the do of (in one equivalent the to fully a region) is possessing that beams show structure Bessel we radial generalized particular, an- In construct the beam. to Bessel of how a engineering of the spectrum completely approach: gular a different of and means new by obtained equation Helmholtz ewitni h olwn form: following the in written be ebgnoraayi ybifl ealn h defini- the recalling briefly by analysis our begin We nti etr epeetanwcaso ouin fthe of solutions of class new a present we Letter, this In = h aerda tutr sLaguerre-Gaussian as structure radial same the , lzeuto bandb aeul engineering carefully by obtained equation oltz p i euirpoaainproperties. propagation peculiar eir K K l l t re eslfnto ftefis id For kind. first the of function Bessel order -th x npriua h aei hc h angular the which in case the particular in A 2 sa nee ubrta xstevleof value the fixes that number integer an is = = E + ( k q k A ( x y x ,y z y, x, ( k , k x 2 ˆ k tc;(6.10 Interference (260.3160) ptics; tan , x 2 x + y + k , = ) 2 k = ) y k , y 3 = ) y 2 θ y ˆ 2 n tan and 1 = π J and l i ( Z l y/x K K 1 d E ϑ 0 0 2 R 0 , δ ( K ϑ hs au r uniquely are value whose , φ ,y, x, ( ) K 0 K e hs lcrcfil can field electric whose , = = ilθ 0 − = k e 0) dk iz y K e /k x k √ i

0 dk 0 K c ) k x sin e · 0 2 ti important is It . R y − ilφ 2018 pnsubsti- Upon . d K ϑ 2 , 0 2 0 K, , and ny gen, J l ( (2) (3) (1) x ) of a LG-beam results in generating Bessel beams with 2 indices, one (l) related to the amount of OAM carried by the beam itself, and the other one (p) related to their radial structure. Notice that we implicitly assumed that both the LG- and Bessel-beams have the same azimuthal index l. Although this is not required by our constrain, without loss of generality we can limit our analysis to this case. To understand how to apply this constrain, let us ﬁrst analyze the case p = 0. The position of the radial maxima of a l-th order Bessel beams are found as the points in which its derivative vanishes, namely

Jl−1(ξ) Jl+1(ξ)=0, (6) Fig. 1. Comparison between the transverse intensity dis- − tribution of (a) Bessel-like beam E2(x, y, 0) and (b) the where ξ = K0R. Since we intend to establish a direct 2 LG-beam LG0(x, y). Note how the central lobe of the connection between the structure of the two beams in Bessel -like beam has the same radial structure of the the paraxial region, we only consider the smallest root LG-beam. The axes of both graphs span the interval of Eq. (6) and we name it ξl. The position of the radial [ 3, 3] in units of the beams waist w of the LG-beam. maximum of a LG-beam with p = 0 is instead given by − 0 R = w0 l /2; if we define η = √2R/w0, being w0 the waist ofp the| | beam, we can introduce the dimensionless parameter η = l . This result is immediate to verify to remark that the diffractionless nature of a Bessel beam l p| | |l| is a direct consequence of its singular angular spectrum. since for p = 0 we have L0 (η) = 1 and the radial part of a LG-beam reduces to η|l| exp( η2/2). The values of Suppose now to operate the following replacement: − Kn to put in Eq. (5) are then found by requiring that N 1 the position of the maxima given by Rl = ξl/Kn for the δ(K K0) Cnδ(K Kn), (4) generalized beam in Eq. (5) and R = w l /2 of the K0 − → X − l 0 n=1 LG-beam coincide. This constrain brings top the| | following where Cn are arbitrary constants and the quantities Kn result: are in general arbitrary subjected to the only constrain ξl sin ϑl = θ0, (7) that K = k sin ϑ < k . By inverting Eq. (2), the elec- 2 l n 0 n 0 | | tric field distribution associated with the general angular p where θ0 = 2/(k0w0) is the beam divergence associated spectrum (4) assumes the form: to the LG-beam [4]. In this case, as we choose the az- N imuthal index of the two beams to be the same, we have 2 2 iz√k0 −Kn ilθ El(x,y,z)= CnJl(KnR)e e . (5) a single index l and the summation in Eq. (5) contains nX=1 only one term. The result for p = 0 is therefore just a single Bessel beam whose angular aperture is determined Equation (5) is still an exact solution of Helmholtz equa- as a function of the amount of OAM carried by the beam tion, and represents a whole new family of Bessel-like itself, namely beams. Note that since Kn = k0 sin ϑn, this solution cor- ik0z cos ϑl responds in general to the superposition of N different El(x,y,z)= ClJl(k0 sin ϑlR)e , (8) Bessel beams, each one characterized by its own angular aperture ϑn. The form and properties of these new beams where we have choses the arbitrary constant Cl = |l| 2 are determined by specifying the explicit form of Kn in ηl exp( ηl /2) to be the value of the LG-beam at its Eq. (4). In the present work, we construct a new class radial maximum.− The reason of this choice will be clear of exact solutions of the Helmholtz equation, which we later, when we will discuss the general case with p =0.A call two-index Bessel beams, by requiring that the radial comparison of this solution with the correspondent6 LG- maxima in the paraxial region of these new beams are beam is shown in Fig. 1. As it can be seen, this operation matched with the radial maxima of a LG-beam. This for p = 0 only results in a rescaling (upon the beam waist requirement translates then in a condition that deter- w0 of the correspondent LG-beam) of a normal Bessel mines the quantities Cn and Kn in Eq. (4). As it is well beam. Before turning our attention to the case p = 0, it known [4], LG-beams are characterized by two indices: a is instructive to analyze the case l =0= p. In this6 case, radial index p and an azimuthal index l, the latter quan- in fact, since the smallest root of Eq. (6) is ξ0 = 0, the tifying the OAM carried by the beam. Bessel beams, on condition expressed by Eq. (7) has no sense anymore. In the other hand, only posses the azimuthal index l, as this particular case, however, we can still rescale the 0- can be seen from Eq. (1). However, the position of the th order Bessel beam with respect to the LG-beam waist radial zeros of a LG-beam depend on both the p and l w0 by simply requiring that sin ϑ0 = θ0. index, and therefore applying as constrain that the ra- We now turn our attention to the most interesting dial maxima of Eq. (5) correspond to the radial maxima case of p = 0. In general, a LG-beam with radial index p 6

2 possesses p+1 radial maxima [22]. In this case, therefore, we will obtain N = p + 1 different expressions of Kn, whose value depend upon the choice of both indices l and p. Each term in the summation (5) therefore represents a Bessel beam whose first radial maximum coincides with the n-th radial maximum of the correspondent LG-beam. The consequence of this is that the beam described by Eq. (5) depends upon two indices: the index l that defines its OAM (as it is done for a standard Bessel beam), and the index p that determines the radial structure of the beam (in the same way that happens for a LG-beam). We call this new solution of the Helmholtz equation a two-index Bessel beam. To find the values of Kn, we recall that the position of the maxima of a LG-beam with p = 0 are found by finding the roots of the following equation:6

|l| 2 2 |l|+1 2 |l| 2 l Lp (η ) η 2Lp−1 (η )+ Lp (η ) =0, (9) | | − h i where as before η = √2R/w0. The condition given by Eq. (7) now generalizes to

ξl sin ϑn = θ0, (10) Fig. 2. Comparison between the transverse intensity dis- √2ηn l tribution of a two-index Bessel beam Ep(x, y, 0) (left l where we indicated with ηn the n-th root of Eq. (9). column) and the correspondent LG-beam LGp(x, y, 0) Inserting this result into Eq. (5) brings to (right column) in z = 0 for. The values of the index l and p are reported n top of each figure. As can be noted, p+1 n the number of rings of the two-index Bessel beams in the E (x,y,z)= C J (k sin ϑ R)eik0z cos ϑ eilθ. (11) pl n l 0 n paraxial regime equals the number of rings of the corre- nX=1 spondent LG-beams. The axes of all the graphs span the It is useful at this point to discuss the form of the interval [ 3, 3] in units of the beams waist w of the − 0 weight coefficients Cn. In order to make the correspon- LG-beam. dence complete, in fact, it is not sufficient to require that the first maximum of each Bessel beam in the summation above corresponds to the n-th maximum of the features and patterns. If we define k = k cos ϑ as the correspondent LG-beam, but it is also necessary that zn 0 n propagation constant of each Bessel beam component, the relative weight between the maxima is also repro- κ = w K and we introduce the dimensionless radial duced. If this condition is fulfilled, in fact, each term of n 0 n coordinate ρ = R/w , the total intensity of a two-index the summation will contribute with the proper weight 0 Bessel beam can be written as to the construction of the generalized beam. We there- p+1 fore require that the coefficients Cn are given by Cn = 2 2 2 |l| (−ηn/2) |l| 2 Ipl(x,y,z)= Epl(x,y,z) = C J (κnρ) ηn exp Lp (ηn). The transverse intensity distribu- | | j l l Xj=1 tion of Ep(x, y, 0) for some values of the two indices p and l, together with their comparison with the corre- p+1 p+1 +2 Λ cos [(k k )z] , (12) spondent LG-beams are depicted in Fig. 2. This is the j,n zj − zn main result of our Letter. Two-index Bessel beams con- nX=1 j=Xn+1 stitute a novel class of exact solutions of the Helmholtz where Λ = C C J (κ ρ) J (κ ρ). Let us analyze three equation, and they can be represented as a superposition j,n j n l n l j distinct cases. If p = 0, the second summation in Eq. (12) of conventional Bessel beam whose cone angle ϑ is cho- n is zero and Eq. (12) reduces to sen (in this particular case) in such a way to give to the two-index Bessel beam the same paraxial structure of a 2 2 I0l(x,y,z)= Cl Jl (κnρ) , (13) LG-beam. We conclude our analysis of this new class of beams by investigating their propagation features. Since and the corresponding beams are diffraction-less. As it the two-index Bessel beams are written in Eq. (5) as a was pointed out before, the case p = 0 reduces the two- superposition of conventional Bessel beams, it is interest- index Bessel beam to the conventional Bessel beam with ing to study the evolution of the intensity of such beams rescaled coordinates with respect to the scaling parame- along the z-direction, as interference between the vari- ter w0. If now we put l = 0 and p arbitrary, we have that ous components of the beam can give rise to interesting κ =2= κ and k k = 0, therefore the intensity n j z,j − z,n

3 References

1. T. H. Maiman, Nature 187, 493 (1960). 2. J. P. Gordon, H. J. Zeiger and C. H. Townes, Phys. Rev. 99, 1264 (1955). 3. N. G. Basov and A. M. Prokhorov, Zh. Ekzperim. i Teor. Fiz., 27 431 (1954). 4. O. Svelto, Principles of Lasers, Springer (2010). 5. J. Durnin, J. Opt. Soc. Am. A 4, 651 (1987). 6. M. A. Banders, and J. C. Gutierrez-Vega, Opt. Lett. 32, Fig. 3. Section (at y = 0) of the evolution of the normal- 3459 (2007). ized intensity of the two-indices Bessel beam Epl(ξ, 0, ζ) 7. M. A. Banders, and J. C. Gutierrez-Vega , Opt. Lett. 33 as a function the normalized coordinate ξ = x/w0 and , 177 (2008). 2 propagation direction ζ = z/zR (being zR = k0w0/2 8. M. A. Banders, and J. C. Gutierrez-Vega, Opt. Express the Rayleigh range of the beam), for (a) (p,l) = (1, 2) 16, 21087 (2008). and (b) (p,l)=(3, 4). The beams are propagating along 9. M. A. Banders, and J. C. Gutierrez-Vega , Opt. Lett. 29 ζ with a periodic pattern, whose period depends essen- , 144 (2004). tially on p. As can be seen from the figures, in fact, as p 10. E. Karimi, G. Zito, B. Piccirillo, L. Marrucci and E. 32 increases, the periodicity of the beam along the propa- Santamato, Opt. Lett. , 3053 (2007). 47 gation direction increases accordingly. 11. M. V. Berry, and N. L. Balazs, Am. J. Phys. , 264 (1979). 12. G. A. Sivioglou, and D. N. Christodoulides, Opt. Lett. 32, 979 (2007). 13. J. Zhao, P. Zhang, D. Deng, J. Liu, Y. Gao, I. D. Chrem- is given by nos, N. K. Efremedis, D. N Christodoulides, and Z. Chen 2 Opt. Lett. 38, 498 (2013). Ip0(x,y,z)=ΓJ (2ρ), (14) 0 14. J. C. Gutierrez-Vega, and M. A. Banders, J. Opt. Soc. Am. A 22, 289 (2005). p+1 2 p+1 p+1 where Γ = j=1 Cj +2 n=1 j=n+1 CnCj . 15. M. A. Banders, and J. C. Gutierrez-Vega, Opt. Lett. 30, Two-indexP Bessel beamsP withP l = 0 are also 2155 (2003). diffraction-less. The last case to consider is the general 16. J. D. Ring, J. Lindberg, A. Mourka, M. Malizu, K. Dho- case where l = 0 and p = 0. In this case, the intensity lakia, and M. R. Dennis, Opt. Express 20, 18955 (2012). 6 6 distribution is periodic with respect to z, with a period 17. M. R. Dennis, and J. D. Ring, Opt. Lett. 38, 3325 that is essentially determined by the index p through the (2013). term kzj kzn inside the cosine function and through 18. M. R. Dennis, J. B. Götte, R. P. King, M. A. Morgan, the fact that− p defines also the number of such cosine and M. A. Alonso, Opt. Lett. 36, 4452 (2011). terms that have to be summed together. However, al- 19. T. H. Lu et al., Phys. Rev. Lett. 101, 233901 (2008). 96 though formally this solution is not diffraction-less any- 20. Y. F. Chen et al., Phys. Rev. Lett. , 213902 (2006). more in the traditional sense (i.e. the intensity distribu- 21. L. Mandel and E. Wolf, Optical Coherence and Quantum tion is not anymore z-independent), Fig. 3 shows that Optics, Cambridge University Press (1995). during propagation, the intensity distribution given by 22. National Institute of Standard and Technology, Dig- ital Library of mathematical Functions (NIST 2010), Eq. (12) reproduces itself periodically, thus exhibiting a http://dlmf.nist.gov quasi-diffraction-free behavior. In conclusion, we have introduced a new class of solutions of the scalar Helmholtz equation based on a gen- eralization of the angular spectrum of a Bessel beam according to Eq. (4). We have show that by choosing Kn in such a way that the maxima of Eq. (5) correspond to the maxima of a LG-beam it is possible to obtain a new class of solutions having two indices that we named two- index Bessel beams. We have shown that these beams are parametrized upon the parameter w0 that in this case is the beam waist of the corresponding LG-beam and that their propagation properties depend upon the choice of the index couple (p,l), and that this class of beams contains full diffraction-less beams as well as periodically self-reproducing beams.

4 Informational Fifth Page References 1. T. H. Maiman, “Stimulated optical radiation in ruby”, Nature 187, 493 (1960). 2. J. P. Gordon, H. J. Zeiger and C. H. Townes, “The Maser - new type of microwave amplifier, frequency standard and spectrometer”, Phys. Rev. 99, 1264 (1955). 3. N. G. Basov and A. M. Prokhorov, Zh. Ekzperim. i Teor. Fiz., 27 431 (1954). 4. O. Svelto, Principles of Lasers, Springer (2010). 5. J. Durnin, “Exact solutions for nondiffracting beams I. The scalar theory”, J. Opt. Soc. Am. A 4, 651 (1987). 6. M. A. Banders, and J. C. Gutierrez-Vega, “Cartesian beams”, Opt. Lett. 32, 3459 (2007). 7. M. A. Banders, and J. C. Gutierrez-Vega , “Circular beams”, Opt. Lett. 33, 177 (2008). 8. M. A. Banders, and J. C. Gutierrez-Vega, “Elliptical beams”, Opt. Express 16, 21087 (2008). 9. M. A. Banders, and J. C. Gutierrez-Vega, “Ince- Gaussian beams”, Opt. Lett. 29, 144 (2004). 10. E. Karimi, G. Zito, B. Piccirillo, L. Marrucci and E. San- tamato, “Hypergeometric-Gaussian modes”, Opt. Lett. 32, 3053 (2007). 11. M. V. Berry, and N. L. Balazs, “Nonspreading wave packets”, Am. J. Phys. 47, 264 (1979). 12. G. A. Sivioglou, and D. N. Christodoulides, “Accelerat- ing finite energy Airy beams”, Opt. Lett. 32, 979 (2007). 13. J. Zhao, P. Zhang, D. Deng, J. Liu, Y. Gao, I. D. Chrem- nos, N. K. Efremedis, D. N Christodoulides, and Z. Chen, “Observation of self-accelerating Bessel-like optical beams along arbitrary trajectories”, Opt. Lett. 38, 498 (2013). 14. J. C. Gutierrez-Vega, and M. A. Banders, “Helmholtz- Gauss waves”, J. Opt. Soc. Am. A 22, 289 (2005). 15. M. A. Banders, and J. C. Gutierrez-Vega, “Vector Helmholtz-Gauss and vector Laplace-Gauss beams”, Opt. Lett. 30, 2155 (2003). 16. J. D. Ring, J. Lindberg, A. Mourka, M. Malizu, K. Dho- lakia, and M. R. Dennis, “Auto-focusing and self-healing of Pearcey beams”, Opt. Express 20, 18955 (2012). 17. M. R. Dennis, and J. D. Ring, “Propagation-invariant beams with quantum pendulum spectra: from Bessel beams to Gaussian beam-beams”, Opt. Lett. 38, 3325 (2013). 18. M. R. Dennis, J. B. Götte, R. P. King, M. A. Morgan, and M. A. Alonso, “Paraxial and nonparaxial polynomial beams and the analytic approach to propagation”, Opt. Lett. 36, 4452 (2011). 19. T. H. Lu, Y. C. Lin, Y. F. Chen and K. F. Huang, “Three-dimensional coherent optical waves localized on trochloidal parametric surfaces”, Phys. Rev. Lett. 101, 233901 (2008). 20. Y. F. Chen, T. H. Lu, K. W. Su and K. F. Huang, “Devil’s staircase in three-dimensional coherent waves localized on Lissajous parametric surfaces”, Phys. Rev. Lett. 96, 213902 (2006). 21. L. Mandel and E. Wolf, Optical Coherence and Quantum Optics, Cambridge University Press (1995). 22. National Institute of Standard and Technology, Dig- ital Library of mathematical Functions (NIST 2010), http://dlmf.nist.gov