Tailoring the focal region of abruptly autofocusing ... - OSA Publishing

0 downloads 0 Views 3MB Size Report
Sep 12, 2018 - Abruptly autofocusing beams, often referred as CABs or ring-Airy .... study of ring-Airy beams is their 1D equivalent, the colliding 1D Airy ..... and D. Sand, “Abruptly autofocusing vortex beams,” Opt. Express 20, 13302–13310.
Vol. 1, No. 1 | 15 Sep 2018 | OSA CONTINUUM 104

Tailoring the focal region of abruptly autofocusing and autodefocusing ring-Airy beams D IMITRIS M ANSOUR 1,2

AND

D IMITRIS G. PAPAZOGLOU 1,2,*

1 Institute

of Electronic Structure and Laser, Foundation for Research and Technology-Hellas, P.O. Box 1527, 71110, Heraklion, Greece 2 Department of Material Science and Technology, University of Crete, P.O. Box 2208, 71003, Heraklion, Greece * [email protected]

Abstract: We theoretically study the properties of the focal region of abruptly autofocusing and, their variant, autodefocusing ring-Airy beams, under the action of a conical phase gradient. By expanding the analysis of 1D Airy beams to cylindrically symmetric Airy beams, we derive analytic formulas for the position and dimensions of the focus. Our analysis covers in a unified way both beam types, and numerical simulations over a broad parameters range are in excellent agreement with our theoretical predictions. Our results allow the tailoring of the focal region both in size and in position by tuning the initial parameters. © 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1.

Introduction

Abruptly autofocusing beams, often referred as CABs or ring-Airy beams, are a subject of intense research effort [1–8], since they were first introduced [9, 10]. Their most exciting property is that instead of spreading due to diffraction, they abruptly autofocus [9, 10], exhibiting a more than two orders of magnitude increase of their intensity. Likewise, in the non-linear propagation regime [3], these beams preserve their abrupt autofocusing, with a minimal nonlinear focal shift. This makes them ideal candidates for applications that require controlled energy deposition, with minimal effects before the focus [10]. For example, ring-Airy beams are advantageous compared to Gaussian or Bessel beams in multiphoton polymerization [11] and in two color laser plasma THz generation [4]. Furthermore, a tighter focus, and thus further increase of the focal intensity, can be achieved by using a lens. As it has been recently shown [12], due to the twin wave nature of ring-Airy beams, focusing by a lens leads to two foci instead of one. This counterintuitive behavior is characteristic for a larger family of waves, referred to as Janus waves [12]. Although this effect is interesting and beneficial to some applications [4], it prevents the use of focusing elements in applications where the twin focus is undesirable. An alternative approach for controlling the focal properties is to use a conical phase gradient. As it has been recently shown [13] for the case of abruptly autofocusing beams this conical phase acts in a cooperative action with the autofocusing properties of the beam and provides some additional control to the focus position. On the other hand, besides the extensive work in the field, to our knowledge analytical expressions that predict the dimensions and position of the focus for abruptly autofocusing and autodefocusing beams are sparse and incomplete. The position of the focus is rather accurately predicted for autofocusing beams [3, 11] even in the case of the presence of a conical phase gradient [13], but to our knowledge no analytical description exists for their autodefocusing variant [14]. Furthermore, analytical formulas that predict the focal dimensions [11] are applicable only to abruptly autofocusing beams for a limited range of beam parameters, and without the presence of a conical phase gradient. #332560 Journal © 2018

https://doi.org/10.1364/OSAC.1.000104 Received 24 May 2018; revised 10 Aug 2018; accepted 11 Aug 2018; published 12 Sep 2018

Vol. 1, No. 1 | 15 Sep 2018 | OSA CONTINUUM 105

Fig. 1. Typical intensity profiles (normalized values) of (a) ring-Airy (b) inverse ring-Airy (r0 =4 mm, w=200 µm), (c) I(r, z = 0) intensity profiles for rotationally symmetric Airy beams of decreasing radius. b=1 (ring-Airy), b=-1 (inverse ring-Airy). Arrows show the direction which the radius is increasing.

In this work, we follow a unified theoretical analysis for both autofocusing, and autodefocusing, ring-Airy beams in the presence of a conical phase gradient. Starting from the analytical expressions that describe the propagation of 1D Airy beams [15], we derive generic analytical expressions that predict the position and dimensions of the focus for both types of beams. Our theoretical predictions are in excellent agreement with numerical simulations, performed over an extended parameters range. 2.

Radially symmetric Airy beams

The 1D Airy beams impart their unique properties to ring-Airy beams, their rotationally symmetric variant. As shown in Fig. 1(a) the intensity profile of a ring-Airy beam [9] is comprised by concentric rings that become thinner and denser as the radius is increased. The beam propagation properties are controlled by the dimensions of the primary ring radius and width [9, 10] but can also be tailored by introducing a tilt though a linear phase gradient [13, 16] which in the case of these rotationally symmetric beams takes the form of a conical phase distribution. On the other hand, the case where the secondary rings extend inwards instead of outwards , as shown in Fig. 1(b), is an alternative way to define a rotationally symmetric Airy beam. We call this less studied [14, 17] rotationally symmetric Airy beam as inverse ring-Airy. The field distribution of both beams can be described as follows: u(r, 0) = uo Ai(bρ)ebaρ e−i sin θkr

(1)

where, Ai(·) denotes the Airy function, ρ = (r0 − r)/w, r is the radius, r0, w are respectively the radius and width parameters, a is an exponential truncation parameter, θ is the equivalent cone angle of the linear phase gradient, k is the wavenumber, and b is actually a switch parameter that can take only two values b = ±1. When b = 1 Eq. (1) describes a ring-Airy beam distribution, as shown in Fig. 1(a), while for b = −1 it describes an inverse ring-Airy beam distribution, as

Vol. 1, No. 1 | 15 Sep 2018 | OSA CONTINUUM 106

shown in Fig. 1(b). Using this unifying description, the primary radius and width, defined as full width at half maximum (FWHM), can be estimated respectively from R0  r0 + b w and FW H M  2.28w. Figure 1(c) shows a graphical representation of radial intensity distribution I(r, z = 0) of ring-Airy and inverse ring-Airy beams as the radius parameter is varied. As one can clearly observe, using the description of Eq. (1), there is a smooth transition from the ring-Airy to the inverse ring-Airy beams. Despite the unified analytical description of Eq. (1) for ring and inverse ring-Airy beams, as shown in Fig. 2 their propagation properties are quite different. Ring-Airy beams exhibit abrupt autofocusing keeping very low intensity up to the focus, as shown in Fig. 2 (a), in contrast to inverse ring-Airy beams that exhibit abrupt autodefocus. Furthermore, as shown in Fig. 2(a), 2(b) by varying the conical phase angle we can shift the focus position along the propagation axis. We have to note here that although the typical parabolic trajectory of accelerating ring-Airy beams [9–12], is shadowed by the presence of the linear gradient, the abrupt autofocusing behavior is clearly evident in the profiles of the peak intensity along z. 3.

Focal region of rotationally symmetric Airy beams

By analytically describing the focus position and spatial dimensions, one is able to tailor the focal region of any beam. For Gaussian beams this is straightforward and is well described in textbooks [18]. On the other hand, although the recently introduced [9, 10], rotationally symmetric beams have been intensively studied [1, 4, 6, 11, 19, 20], such an analytic description is partially available for ring-Airy beams [11, 13] but, to our knowledge, not for inverse ring-Airy beams. In the following we will provide analytic descriptions of the focal position and spatial dimensions of these rotationally symmetric Airy beams, that enable us to tailor their focal region. 3.1.

Colliding Airy beams

An excellent model for the study of ring-Airy beams is their 1D equivalent, the colliding 1D Airy beams [12, 21]. These beams consist of two mirror symmetric 1D Airy beams that propagate along parabolic trajectories. As they propagate these beams coherently interference leading to a focal point along the propagation axis, similar to the autofocus of the ring-Airy beams. Since analytical formulas [15, 16] can be used to describe the propagation of an 1D Airy beam we can tailor the focusing behavior of the two colliding 1D Airy beams and expand this to cylindrically symmetric Airy beams. The conical linear phase gradient introduced in Eq. (1) is in this 1D scheme equivalent to an initial tilt angle for each beam. Using the analytical formulas for the propagation of tilted 1D Airy beams [16] we can describe the field distribution of the colliding 1D Airy beams: U1D (x, z; θ) = Φ(x, z, θ) + Φ(−x, z, −θ), Φ(x, z, θ) = Ai(b ×e

ba

x 0 +x−θ z w



a z2 4k 2 w 4

x0 + x z2 az z − 2 4 + i 4 − bθ ) w w 4k w kw

+i(b x o +x+b3a 2k w

2w

z−

z3 12k 3 w 6

+θk(x0 +x)−b

θ z2 2k w 3

−θ

2kz 2

(2)

)

,

where x, x0 are respectively the transverse coordinate and initial displacement, w is the Airy width parameter, θ is the equivalent cone angle of the initial phase gradient exp(−i sin θ k x), k is the wavenumber, and a is the truncation coefficient. As in Eq. (1) the parameter b can take only two values. For b = 1 we get the 1D equivalent of the autofocusing ring-Airy beam, while for b = −1 we get the 1D equivalent of the abruptly autodefocusing inverse ring-Airy beam. For small values of the truncation coefficient we can predict analytically, the position of the focus: h i 1/2 1D f Ai  4bz Ai −2θ˜ ± (4θ˜2 + b x) ˜ , (3)

Vol. 1, No. 1 | 15 Sep 2018 | OSA CONTINUUM 107

Fig. 2. Intensity I(x, z) profiles of rotationally symmetric Airy beams and peak intensity along z with different conical phase gradients at the input plane (z=0) θ i=1,2,3 =30 mrad, 20 mrad, 10 mrad (r0 =1.0 mm and w=100 µm, u0 =1) (a) ring-Airy (b) inverse ring-Airy. Insets show the position and the dimensions of the focus.

where z Ai ≡ kw 2 /2 , x˜ ≡ xo /w + b is a shape factor and θ˜ ≡ z Ai θ/w is a normalized cone angle. The two solutions in Eq. (3) are a result of the Janus nature of these waves [12], and correspond to the foci of the "real" and "virtual" waves. Furthermore, by properly manipulating the initial shift and linear gradient we can tailor the distribution of the intensity maximum along propagation. This focus tailoring technique can then be generalized to cylindrically symmetric Airy beams by introducing a conical phase, equivalent to that introduced by an axicon [22]. As shown in Fig. 3 there is a strong analogy between the autofocusing colliding 1D Airy beams Fig. 3(i) and the ring-Airy Fig. 3(ii), as well as between the autodefocusing colliding 1D Airy beams Fig. 3(iv), and the inverse ring-Airy beams Fig. 3(v). Besides the position of the focus, the intensity distribution along the propagation axis, especially in the focal region, is practically identical as shown in Figs. 3(iii), 3(vi). On the other hand, the I(x, z) intensity profiles, as shown in Figs. 3(iii), 3(iv) insets are quite different especially along x direction. This is expected since in one case (colliding Airy) we have the interference of two beams and in other (radially symmetric Airy) the interference of a conical wavefront. Using the analytic predictions of Eq. (2) that refer to the propagation of colliding 1D Airy beams, and taking into account the above-mentioned similarities and differences, we retrieve approximate analytic expressions for the focal properties of cylindrically symmetric Airy beams.

Vol. 1, No. 1 | 15 Sep 2018 | OSA CONTINUUM 108

Fig. 3. Comparison of intensity I(x, z) along propagation for autofocusing colliding 1D Airy beams (1st row) with with radially symmetric Airy beams (2nd row). 3d row: Normalized intensities along the z propagation axis I(0, z)/Imax . Insets: I(x, z) intensity profiles (normalized, false colors), (1D) colliding 1D Airy beams (2D) radially symmetric Airy beams. (Numerical simulation parameters: λ=0.8 µm, α = 0.05 (i) x0 = −2.0 mm, w=100 µm, θ=10 mrad (ii) r0 =2.0 mm, w=100 µm, θ=10 mrad (iv) x0 =2.0 mm, w=100 µm, θ=10 mrad (v) r0 =2.0 mm, w=100 µm, θ=10 mrad).

3.2.

Cylindrically symmetric Airy beams

A unique feature of ring-Airy beams is that they abruptly autofocus along propagation, while their peak intensity remains very low up to their abrupt focus. The basic parameters that describe the focus of such a beam are depicted in the insets of Fig. 2(a), 2(b). These include the effective focal length f Ai [10], which denotes the position of the abrupt focus, the focal spot length ∆ f Ai which denotes the longitudinal (FWHM) of the intensity focal distribution, and focal spot width w Ai which denotes the transverse dimensions of the focus (FWHM). Recently approximate analytic relations [11] were introduced for the estimation of abruptly autofocusing ring-Airy focal parameters f Ai , ∆ f Ai , w Ai as a function of the initial radius r0 and width w parameters, valid only for r0 /w  1. Using the equivalence to 1D+1 colliding Airy beams, we expanded this analysis to cover abruptly autofocusing and autodefocusing ring-Airy beams, under the action of a conical phase gradient. Furthermore, our unified analytical description now covers an extended initial parameters range. In more detail, using the analytic prediction of the trajectory of the peak of the primary intensity lobe of Eq. (3) to radially symmetric Airy beams shaped by a conical phase we get for the focus position: h i f Ai = 4b −2θ˜ ± (4θ˜2 + bs)1/2 , (4) z Ai where z Ai ≡ kw 2 /2 , s ≡ ro /w + b is a shape factor and θ˜ ≡ z Ai θ/w is a normalized cone angle.

Vol. 1, No. 1 | 15 Sep 2018 | OSA CONTINUUM 109

Fig. 4. Focal distance f Ai , of abruptly autofocusing ring-Airy beams normalized over z Ai as a function, (a) of the shape factor s (θ˜ = 0), (b) of the normalized cone angle θ˜ (s=21).

As in the case of 1D colliding Airy beams, Eq. (3) the Janus nature [12] of these waves leads to two solutions for the position of the focus. For the case of θ˜ = 0 we get that the two foci are symmetrically positioned relative to z = 0. The next crucial focal parameter is the focal spot length ∆ f Ai . Taking into account the trajectory of the 1D toy model of Eq. (3), we can estimate the projection of the primary lobe on the axis at the focus. This enables us to reach to an approximate analytic estimation for the focal spot length:     1/2  −1 ∆ f Ai 42b − 1  ˜2 2 ˜ ˜ 4θ + bs − 2θ 4θ + bs  . (5) f Ai 50 The case of the focal spot width w Ai is more complicated. The cylindrically symmetric nature of the problem enables us to locally treat the focal spot as a result of the interference of a conical wavefront. This approach makes it possible to use the well-known analytic estimations for the spot size of Bessel beams [23] where the wavefront cone angle at focus is predicted by the parabolic trajectory of the 1D toy model as described in Eq. (3). Thus, we reach to an approximate analytic estimation of the focal spot width for both ring-Airy and inverse ring-Airy Eq. (6) beams: ! 1/2   w Ai 1 4 w2  C 1− + , (6) w 2 − 2bθ˜2 + bs 4θ˜2 + bs z 2Ai where C ≡ (76 + 22b)/(65 + 15b) is a constant. 3.2.1.

Abruptly autofocusing ring-Airy beams

To confirm the validity of our analytic predictions as described in Eqs. (4) – (6) we used numerical simulations of the paraxial wave equation, for a broad range of parameters. More specifically we study the evolution of the normalized focal parameters f Ai /z Ai , ∆ f Ai / f Ai , w Ai /w as a function ˜ In our of two parameters; the normalized shape factor s and the normalized conical angle θ. simulations we either vary the normalized shape factor s keeping a constant angle θ˜ = 0 or we vary θ˜ keeping s constant. In the first case, since θ˜ = 0, our results are comparable to the estimations presented in [11]. We have to note though that Eqs. (4) – (6) are valid for a much larger initial parameters range from r0 /w > 1, θ˜ > 0 up to r0 /w ≈ 40, θ˜ ≈ 15 rad .

Vol. 1, No. 1 | 15 Sep 2018 | OSA CONTINUUM 110

Fig. 5. Ring-Airy focal spot length normalized over the focal distance f Ai as a function (a) of the shape factor s (θ˜ = 0), (b) of normalized cone angle θ˜ (s=21).

Fig. 6. Focal spot width w Ai of ring-Airy, normalized over the width parameter w as a function (a) of the shape factor s (θ˜ = 0) and (b) of the normalized cone angle θ˜ (s=21).

Figure 4(a) shows the position of the focus f Ai as a function of the primary ring parameters. Numerical simulations are in both cases in excellent agreement with the theoretical predictions of Eq. (4). As the √ shape factor s is increased, the normalized focal distance is increased proportionally to s as predicted by Eq. (4) for θ˜ = 0. It should be noted here that increasing s by decreasing the ring width w parameter also changes z Ai . On the other hand, Fig. 4(b) depicts the effect of the conical phase gradient on the focus position. As expected, increasing the ˜ conical angle θ˜ decreases the normalized focal distance. Clearly, by varying the conical angle θ, instead of the shape factor s, we achieve a similar dynamic range in controlling the focus position, although f Ai saturates for θ˜ >6 rad. Figure 5 shows the focal spot length ∆ f Ai as a function of the primary ring parameters. Again the numerical simulation results are in very good agreement

Vol. 1, No. 1 | 15 Sep 2018 | OSA CONTINUUM 111

to the analytical predictions of Eq. (5). As the shape factor s is increased the normalized focal spot length is decreased. Combining this with the fact that the focal distance f Ai is increased, as shown in Fig. 4(a), we get that as the focus is shifted further away it becomes shorter compared to the focal distance f Ai . On the other hand, as shown in Fig. 5(b), when the conical angle is increased the the normalized focal length quickly saturates after θ˜ > 3 rad. Likewise, Fig. 6 shows the dependence of the focal spot width w Ai , normalized over the width parameter w as function of the initial parameters. The analytical predictions of Eq. (6) are in very good agreement with the numerical simulation results. As the shape factor s is increased, the normalized focus width monotonically decreases. A similar behavior is observed in Fig. 6(b), where the cone angle θ˜ is increased. The shape of the focal region is well described by it’s aspect ratio, defined as the ratio of the focal length to the focal width AR = ∆ f Ai /w Ai . Figure 7 depicts the analytically estimated values of the aspect ratio AR for autofocusing ring-Airy beams as a function of the radius ro and width w parameters. The results in Fig. 7(a) are retrieved for a conical angle θ = 0 while in Fig. 7(b) for a conical angle θ=5 mrad. The contour lines show that a constant aspect ratio, i.e. a preservation of the focus shape, can be achieved by controlling the shape parameter through the radius parameter ro while keeping the width parameter w constant. This interesting behavior was first demonstrated by [11] for the case of θ = 0. Here we expand this result for the case of a conical angle θ > 0. As it can be seen from Fig. 7 a typical characteristic of these beams is the needle-like shape of the focal region since the aspect ratio is quite high (AR > 100). This from the practical point of view, and in respect of the shape of the focal region, sets autofocusing ring-Airy beams between Bessel beams [23] where the aspect ratio is much higher and typical Gaussian beams where the aspect ratio can be quite lower.

0 90

700

300



500

 100

UR PP









(a)



Z ȝP

0 90

700

300



500

 100

UR PP









(b)



Z ȝP $5 







  

Fig. 7. Aspect ratio ∆ f Ai /w of the focal region of a ring-Airy as a function of the radius and width parameters (ro, w). (a) θ = 0, (b) θ = 5 mrad.

Vol. 1, No. 1 | 15 Sep 2018 | OSA CONTINUUM 112

Fig. 8. Position of the focus f Ai of an autodefocusing inverse ring-Airy beam, normalized over z Ai . (a) as a function of the shape factor s (θ˜ = 7.854 rad), (b) as a function of the normalized cone angle θ˜ (s=19).

3.2.2.

Abruptly autodefocusing inverse ring-Airy beams

In order to reveal their abrupt autodefocusing behavior inverse ring-Airy beams, and in contrast to autofocusing ring-Airy beams, require the presence of a linear conical phase. This a conical phase gradient is represented by conical angle θ in Eq. (1). As shown in Fig. 2(b) when this conical angle exceeds a value of e θ > (s/4)1/2 abrupt autodefocusing behavior is observed . As in the case of autofocusing ring-Airy beams, we used numerical simulations of the paraxial wave equation, to confirm the validity of our analysis for a broad range of parameters, in the s < 4θ˜2 regime. Figure 8(a) shows the position of the focus f Ai of an autodefocusing inverse ring-Airy beam as a function of the scale parameter s.The simulations are performed for a conical angle θ = 10 mrad (θ˜ = 7.854 rad). Numerical simulations are in excellent agreement with the theoretical predictions of Eq. (4). The normalized focus position shows a practically linear dependence on the shape parameter s. Indeed from Eq. (4) and by setting b = −1 we get f Ai /z Ai  s/θ˜ for s  4θ˜2 . Likewise, Fig. 8(b) shows the effect of the conical phase gradient on the focus position of the abruptly autodefocusing ring-Airy. Again, the numerical simulations are in excellent agreement to the predictions of Eq. (4). The normalized focus position exhibits a monotonic decrease as the conical angle is increased, in similar fashion the the case of the abruptly autofocusing ring-Airy beams shown in Fig. 4. Figure 9(a) depicts the dependence of focal spot length ∆ f Ai of an autodefocusing inverse ring-Airy beam as a function of the scale parameter s. The numerical simulation results are in excellent agreement to the analytical predictions of Eq. (5). Similarly to the case of the ring-Airy beams shown in Fig. 5(a) the normalized focal length shows a monotonic decrease as the s parameter is increased. On the other hand, as shown in Fig. 9(b), although there is a very good agreement between the numerical simulations and the analytical predictions for θ˜ > 7 rad, Eq.(5) underestimates ∆ f Ai / f Ai by < 5% when θ˜ < 6 rad. This is due to the fact that for low values of the cone angle θ˜ we are approaching the validity limit of our analysis which in holds for 4θ˜2 > s. Figure 10(a) depicts the focal width w Ai of an autodefocusing inverse ring-Airy beam, normalized over the width parameter w, as a function of the shape parameter s. The analytic

Vol. 1, No. 1 | 15 Sep 2018 | OSA CONTINUUM 113

Fig. 9. Inverse ring-Airy focal spot length ∆ f Ai normalized over f Ai as a function (a) of the shape factor s (θ˜ =7.854 rad). (b) of the normalized cone angle θ˜ (s=29).

Fig. 10. Focus width w Ai of an abruptly autodefocusing inverse ring-Airy, normalized over the width parameter w as a function (a) of the shape factor s (θ˜ =10.21 rad). (b) of the normalized cone angle θ˜ (s=19).

predictions are again in excellent agreement to the numerical simulations. Although the normalized focal width w Ai exhibits a linear dependence on the shape parameter s, this holds only for if the width parameter w remains constant, as we can deduce from Eq. (6) for the case of s  4θ˜2 . The dependence of w Ai /w on the conical phase angle θ˜ is shown in Fig. 10(b). While the fit between analytical predictions and numerical simulation results is excellent for θ˜ >8 rad, Eq. (6) overestimates by 10% for θ˜ w and ro /w < 2(1 + θ˜2 ) as shown by the respective curves. The contour lines mark the areas of constant aspect ratio. As in the case of abruptly autofocusing ring-Airy beams AR , and thus the focal voxel shape, is preserved when the radius parameter ro is increased while keeping the width parameter w constant. Likewise, the aspect ratio takes large values (AR > 100), indicating a needle-like focal region. 4.

Conclusions

We have studied the properties of the focal region of abruptly autofocusing ring-Airy, and their variant autodefocusing inverse ring-Airy beams, under the action of a conical phase gradient. Exploiting the similarities between the propagation of 1D Airy beams to cylindrically symmetric Airy beams, we present, using a unified formulation for both beam types, analytic approximations for the position, length and width of the focus. Our results enable the tailoring of the focal region of such abruptly autofocusing and autodefocusing beams by simply tuning the initial radius ro and width w parameters. Furthermore, We have confirmed the validity of our analytical results by numerical simulations over a broad parameters range. We expect that our results will have a significant impact to applications that exploit the unique properties of cylindrically symmetric Airy beams like for instance, materials processing, optical trapping, and deposition of high-laser powers at remote locations.

Vol. 1, No. 1 | 15 Sep 2018 | OSA CONTINUUM 115

Funding General Secretariat for Research and Technology (GSRT); Hellenic Foundation for Research and Innovation (HFRI) (HRFI-844). Acknowledgments D.M. acknowledges financial support from the General Secretariat for Research and Technology (GSRT) and the Hellenic Foundation for Research and Innovation (HFRI) (HRFI-844). References 1. I. Chremmos, P. Zhang, J. Prakash, N. K. Efremidis, D. N. Christodoulides, and Z. Chen, “Fourier-space generation of abruptly autofocusing beams and optical bottle beams,” Opt. Lett. 36, 3675–3677 (2011). 2. I. Chremmos, N. K. Efremidis, and D. N. Christodoulides, “Pre-engineered abruptly autofocusing beams,” Opt. Lett. 36, 1890–1892 (2011). 3. P. Panagiotopoulos, D. G. Papazoglou, a. Couairon, and S. Tzortzakis, “Sharply autofocused ring-Airy beams transforming into non-linear intense light bullets,” Nat. Commun. 4, 2622 (2013). 4. K. Liu, A. D. Koulouklidis, D. G. Papazoglou, S. Tzortzakis, and X.-C. Zhang, “Enhanced terahertz wave emission from air-plasma tailored by abruptly autofocusing laser beams,” Optica 3, 605–608 (2016). 5. H. Deng, Y. Yuan, and L. Yuan, “Annular arrayed-waveguide fiber for autofocusing Airy-like beams,” Opt. Lett. 41, 824–827 (2016). 6. Y. Jiang, K. Huang, and X. Lu, “Propagation dynamics of abruptly autofocusing Airy beams with optical vortices,” Opt. Express 20, 18579–18584 (2012). 7. J. A. Davis, D. M. Cottrell, and D. Sand, “Abruptly autofocusing vortex beams,” Opt. Express 20, 13302–13310 (2012). 8. H. Chi-Young, K. Kyoung-Youm, and L. Byoungho, “Dynamic control of circular Airy beams with linear optical potentials,” IEEE Photon. J. 4, 174–180 (2012). 9. N. K. Efremidis and D. N. Christodoulides, “Abruptly autofocusing waves,” Opt. Lett. 35, 4045–4047 (2010). 10. D. G. Papazoglou, N. K. Efremidis, D. N. Christodoulides, and S. Tzortzakis, “Observation of abruptly autofocusing waves,” Opt. Lett. 36, 1842–1844 (2011). 11. M. Manousidaki, D. G. Papazoglou, M. Farsari, and S. Tzortzakis, “Abruptly autofocusing beams enable advanced multiscale photo-polymerization,” Optica 3, 525–530 (2016). 12. D. G. Papazoglou, V. Y. Fedorov, and S. Tzortzakis, “Janus waves,” Opt. Lett. 41, 4656–4659 (2016). 13. J. Zhang, Y. Li, Z. Tian, and D. Lei, “Controllable autofocusing properties of conical circular Airy beams,” Opt. Commun. 391, 116–120 (2017). 14. C.-Y. Hwang, K.-Y. Kim, and B. Lee, “Bessel-like beam generation by superposing multiple Airy beams,” Opt. Express 19, 7356–7364 (2011). 15. G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Observation of Accelerating Airy Beams,” Phys. Rev. Lett. 99, 213901 (2007). 16. G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Ballistic dynamics of Airy beams,” Opt. Lett. 33, 207–209 (2008). 17. P. Zhang, J. Prakash, Z. Zhang, M. S. Mills, N. K. Efremidis, D. N. Christodoulides, and Z. Chen, “Trapping and guiding microparticles with morphing autofocusing Airy beams,” Opt. Lett. 36, 2883–2885 (2011). 18. M. Born and E. Wolf, Principles of Optics (Elsevier, 1980), 6th ed. 19. I. D. Chremmos, Z. Chen, D. N. Christodoulides, and N. K. Efremidis, “Abruptly autofocusing and autodefocusing optical beams with arbitrary caustics,” Phys. Rev. A 85, 023828 (2012). 20. R.-S. Penciu, K. G. Makris, and N. K. Efremidis, “Nonparaxial abruptly autofocusing beams,” Opt. Lett. 41, 1042–1045 (2016). 21. Y. Zhang, M. R. Belić, H. Zheng, H. Chen, C. Li, Y. Li, and Y. Zhang, “Interactions of Airy beams, nonlinear accelerating beams, and induced solitons in Kerr and saturable nonlinear media,” Opt. Express 22, 7160–7171 (2014). 22. J. H. McLeod, “The Axicon: A new type of optical element,” J. Opt. Soc. Am. 44, 592–597 (1954). 23. J. Durnin, J. J. Miceli, and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58, 1499–1501 (1987).