## Abstract

The Airy beams are analyzed in order to provide a cogent physical explanation to their intriguing features which include *weak diffraction, curved* propagation trajectories in free-space, and *self healing*. The asymptotically exact analysis utilizes the method of uniform geometrical optics (UGO), and it is also verified via a uniform asymptotic evaluation of the Kirchhoff-Huygens integral. Both formulations are shown to fully agree with the exact Airy beam solution in the paraxial zone where the latter is valid, but they are also valid outside this zone. Specifically it is shown that the beam along the curved propagation trajectory is not generated by contributions from the main lobe in the aperture, i.e., it is not described by a local wave-dynamics along this trajectory. Actually, this beam is identified as a caustic of rays that emerge sideways from points in the initial aperture that are located far away from the main lobe. The field of these focusing rays, described here by the UGO, fully agrees with the Airy beam solution. These observations explain that the “weak-diffraction” and the “self healing” properties are generated, in fact, by a continuum of sideways contributions to the field, and not by local self-curving dynamics. The uniform ray representation provides a systematic framework to synthesize aperture sources for other beam solutions with similar properties in uniform or in non-uniform media.

© 2010 Optical Society of America

## 1. Introduction

Recently, a class of “curving beam” solutions of the paraxial wave equation has been introduced [1–7]. These solutions were also termed “accelerating beams,” borrowing from the original work of Berry and Balazs [8] that has dealt with accelerating wavepacket solutions of the time-dependent potential-free Schrödinger equation, and the similarity between that equation and the paraxial (parabolic) wave equation. Originally, these beams were formulated in a 2D coordinate space, say (*x*,*z*), and were generated by setting an Airy function as the initial field distribution in the aperture plane *z* = 0. Later on, finite energy beams were obtained by multiplying the Airy function distribution by an exponential window or by a Gaussian window. Closed form solutions for the radiated field in these cases were derived in [2] and [3], respectively. Extensions to 3D-space have been presented in [2,6,7]. Henceforth we shell refer mainly to the “Airy beams” as representatives of this class of waves.

The Airy beams attracted a lot of attention because of their intriguing features. Their most distinctive characteristic is the propagation along *curved* trajectories in free-space [2,9]. These beams are also weakly diffractive along their trajectories, i.e., they retain their structure and remain essentially diffraction-free for distances that are much longer than the conventional diffraction length of Gaussian beams with the same beamwidth [2]. Another interesting feature, pointed out in [10], is the ability of the Airy beam to “heal” itself, i.e., regenerate itself when part of the aperture is obstructed.

Our aim here is to clarify these features viz a viz what we consider to be the basic definition of a “beam.” We regard a beam as *a local wave function whose evolution along its propagation trajectory is described by a local self-dynamics* (for a broader perspective see [11] and references therein). Thus we would like to identify the wave mechanisms that give rise to the special properties of the Airy beams.

While previous studies of the Airy beams were based on their closed form solutions, here we shall explore their properties using basic wave theory. We present an asymptotically exact wave analysis based on the method of uniform geometrical optics (UGO) [12,13], and also on the Kirchhoff-Huygens integral. Both formulations are shown to fully agree with the exact Airy beam solution in the paraxial zone where the latter is valid, but they are also valid outside this zone. The results provide a full wave explanation to the various phenomena associated with the Airy beams. For simplicity, we shall present here only the 2D case and focus on the Airy beam as a representative of this class of waves, although these concepts can easily be extended to the 3D case or to the finite energy Airy beams.

The presentation starts in Sec. 2 with a short review of the Airy beams. In Sec. 3 we consider the field generated by the main lobe of the Airy function initial distribution. It is demonstrated that this contribution decays faster than a Gaussian beam with the same width, thus confirming that the Airy beam is not generated by a self-curving dynamical evolution of this main lobe. The field radiated by the other part of the initial field distribution is calculated next in Sec. 4 using uniform geometrical optics (UGO), and the results are also verified by a uniform asymptotic evaluation of the Kirchhoff-Huygens integral. The curving Airy beam is shown to be a caustic of rays emerging sideways from points in the aperture which are located far away from the main lobe. The total field is calculated in Sec. 5, where it is shown that the UGO formulation agrees extremely well with the closed form solution of the parabolic wave equation in the range where the latter is valid. Beyond this range, the closed form solution fails whereas the UGO formulation remains valid. This ray interpretation verifies that the curving Airy beams *are not* “beam fields” in the sense discussed above. This interpretation also explains the weak diffraction and the self healing properties. It also provides a framework for synthesizing other aperture distributions that give rise to beam fields that propagate along any desired trajectory, and may also be applied to non-uniform media. These and other conclusions are discussed in Sec. 6.

## 2. Preliminaries

We begin by considering the 2D paraxial wave equation in free-space, that governs the propagation of the electric field amplitude *ϕ* associated with planar optical beams

where *z* is the coordinates along the propagation axis, *x* is the transverse coordinate, and *k* = 2*π*/*λ* is the wavenumber. A time dependence exp(-*iωt*) is assumed and suppressed throughout. The Airy beam is generated by an initial field distribution in the aperture plane *z* = 0, given by

where Ai is the Airy function and *x*_{0} is an arbitrary scaling factor. Since the *L*_{2} norm of the Airy function is unbounded, one may multiply the initial distribution in Eq. (2) by a window function *F*(*x*) in order to generate a finite energy beam. Possible choices for *F*(*x*) are the exponential and Gaussian functions, resulting in the exponential Airy beam [2] or the Airy-Gauss beam [3], respectively. For simplicity, we shell study only the basic Airy beam resulting from the initial distribution in Eq. (2) but the same techniques may be used for the finite energy Airy beams as discussed in the conclusions.

The field *U* at **r** = (*x*,*z*) associated with the initial distribution in Eq. (2) is given by

where we used the normalized coordinates *x*̃ = *x*/*x*_{0} and *z*̃ = *z*/*z*_{0} with *z*_{0} = *kx*^{2}_{0}. One readily verifies that the beam envelope shifts transversely without change along a parabolic trajectory *x*̃ = z̃^{2}/4, as may also be seen in Fig. 1.

In this work we shall employ scalar diffraction theory in order to calculate the field generated by the aperture field at *z* = 0. Using the Kirchhoff-Huygens integral with the Rayleigh-Sommerfeld Green function, the radiated field is given by

where **r**′ = (*x*′, *z*′) with *z*′ = 0 denotes points in the aperture, ∂_{n′} is the derivative with respect to the normal at **r**′, *G* = (*i*/4)*H*^{(1)}_{0}(*kρ*) is the Green function of the 2D free-space wave equation, *H*^{0}^{(1)} is the zeroth order first kind Hankel function, and *ρ* = [(*x*-*x*′)^{2}+*z*^{2}]^{1/2} = ∣**r**−**r**′∣. For observation points where *kρ* ≫ 1, *G*~exp(*ikρ*+*iπ*/4)/(8*πkρ*)^{1/2}. In this case *∂ _{n′}*

*G*~ −

*ik*cos

*θG*where

*θ*= cos

^{−1}

*z*/

*ρ*is the angle between the vector from

**r**′ to

**r**and the

*z*axis.

In the sequel we consider separately the fields *U*_{1} and *U*_{2} representing, respectively, the contributions of the main lobe in the aperture and of the rest of the aperture. They are defined as

where *x*_{ep} = *α*_{1}*x*_{0} with *α*_{1} = -2.338 being the first zero of Ai(*x*). These two contributions will be calculated, separately, in Secs. 3 and 4, respectively.

## 3. The field *U*_{1} generated by the main lobe

The field *U*_{1} due to the main lobe of the Airy function initial distribution is calculated numerically via Eq. (5), and compared with the field of a Gaussian beam (GB) having the same effective half-width, defined as

Here 〈*x*〉 and *I* are, respectively, the centroid and the intensity of the aperture distribution of *U*_{1}, defined via

As a reference, we compare *U*_{1} to the field due to a Gaussian initial distribution *U*_{GB}∣_{z=0} = *I*_{0}^{1/2} exp [ − (*x* − 〈*x*〉)^{2}/2*W*^{2}_{0}], centered about *x* = 〈*x*〉, where *W*_{0} is the *e*^{−1} half width. This distribution has the same effective width *W* and intensity *I* as in Eqs. (6)–(7) if *W*_{0} = √2*W* and ${I}_{0}=I/\left(W\sqrt{2\pi}\right).$. The GB field radiating by this distribution is *U*_{GB}(**r**)=*I*_{0}^{1/2}[(−*iz _{F}*)/(

*z*−

*iz*)]

_{F}^{1/2}exp[

*ik*(

*z*+(

*x*−〈

*x*〉)

^{2}/2(

*z*−

*iz*))] where

_{F}*z*=

_{F}*kW*

^{2}

_{0}is the Fresnel (or Rayleigh) length.

Figure 2 compares the intensity of *U*_{1} calculated via Eq. (5) to the intensity of the field *U*_{GB} mentioned above and also to the Airy beam of Eq. (3), at different ranges. Here and throughout *λ* = 0.5*μm* and *x*_{0} = 200*λ*, giving *W* = 139.4*λ* and *z _{F}* = 2

*kW*

^{2}= 24.410

^{4}

*λ*. In order to present the results using physically meaningful scales, we normalize the

*z*and

*x*axes with respect to

*z*and

_{F}*W*, respectively. One readily observes that

*U*

_{1}does not generate the weakly diffractive Airy beam but it rather behaves like the matched

*U*

_{GB}above. In fact, it diffracts a little faster: At

*z*= √3

*z*, the effective width of

_{F}*U*

_{GB}and of

*U*

_{1}increase by the factors of 2 and 2.3, respectively, while at

*z*= √8

*z*, the factors are 3 and 3.6.

_{F}To conclude this section we note that the main curved lobe of the Airy beam is not generated by the first main lobe of the aperture field, but by some other part, which is identified in the next section. In other words, the curved lobe of the Airy beam is not described by a local self-curving wave dynamics, as was explain in the Introduction.

## 4. The field *U*_{2} generated by the aperture beyond the main lobe

We now turn to calculate the field *U*_{2} in Eq. (5). The contribution of this part will be described in terms of rays that emerge from the aperture in directions specified by the local phase gradient. Special care should be given to the field near the curving main beam, where the rays focus to a caustic. To this end, the ray field is expressed in Sec. 4.1 using uniform geometrical optics (UGO). This provides a complete and accurate representation for the field of the Airy beam.

The artificial separation of the aperture field distribution into two parts at *x* ≶ *x*_{ep} gives rise to a spurious end-point diffraction at *x*_{ep} which should also be added to the asymptotic representation of *U*_{2}. This contribution is derived in Sec. 4.2 using the uniform asymptotic theory (UAT), which is the uniform variant of the geometrical theory of diffraction (GTD). Clearly, when the total field *U*_{1} + *U*_{2} is calculated, this end-point contribution is canceled out by the truncation effect which is included implicitly in the numerical calculation of *U*_{1}, leaving the total field described only by the UGO contribution, as discussed in Sec. 5.

The field representation above is obtained by following the rules of uniform ray theory. For completeness, they are verified rigorously in Sec. 4.3 using uniform asymptotic evaluation of the kirchhoff integral in Eq. (5).

Readers who are not interested in the analytic details may skip Secs. 4.2–4.3 and proceed directly to Sec. 5 which includes a summary of the main results and of the numerical results.

#### 4.1. The uniform geometrical optics (UGO) field

In order to facilitate ray analysis, we use the asymptotic expression of Ai(*x*) for *x*≪−1 and decompose the aperture field at *x* ≤ *x*_{ep} into a sum of two local plane wave constituents, viz,

$$\phantom{\rule{2.em}{0ex}}={A}_{0}^{+}\mathrm{exp}\left(ik{\psi}_{0}^{+}\right)+{A}_{0}^{-}\mathrm{exp}\left(ik{\psi}_{0}^{-}\right)$$

$$\phantom{\rule{1.8em}{0ex}}\equiv {\varphi}_{0}^{+}\left(x\right)+{\varphi}_{0}^{-}\left(x\right),$$

where *ψ*^{±}_{0} and *A*^{±}_{0}, the initial phase and amplitude functions, are given by

The initial fields *ϕ*_{0}^{+} and *ϕ*_{0}^{−} in Eq. (8) give rise to rays that emerge from points *x*′ ≤ *x*_{ep} in the aperture at local angles (measured with respect to the *z*-axis)

The rays follow the trajectories

where *x*′ is the exit point and *ρ* measures the distance along the ray (see Figs. 1 and 3). The ray species *ϕ*^{+}_{0} emerges sideways, converging to a caustic along the curving main beam trajectory *x*̃ = *z*̃^{2}/4 (this expression is valid within the parabolic approximation; a more accurate
expression is discussed below). Each observation point on the lit side of the caustic (region II in Fig. 1) is illuminated by two rays of species *ϕ*^{+}_{0}, one which has touched the caustic before reaching the observer, and another which has not. They are denoted, respectively as rays 2 and 1, referring to the fact that the phase of ray 1 is smaller than that of ray 2 (implying that in the time domain, ray 1 reaches the observer first).

Due to the artificial truncation of the aperture at *x*_{ep}, the rays leaving this point form artificial shadow boundaries for the *ϕ*^{+}_{0} and *ϕ*^{−}_{0} species denoted, respectively, as SB^{+} and SB^{−} (see Fig. 1). Referring first to species *ϕ*^{+}_{0}, this implies that the caustic is truncated at the point where it touches SB^{+} (see Fig. 3). The shadow side of SB^{+} (region III) is illuminated only by ray 1 of this species.

Finally, we note that the lit side of SB^{−} corresponding to species *ϕ*^{−}_{0} is far from the main beam region and therefore this contribution will not be considered below.

The geometrical optics (GO) field at points **r** along a ray of the *ϕ*^{+}_{0} species emerging from a given point *x*′ [see Eq. (11)] is given by

The phase along this ray is

where the initial phase *ψ*_{0}^{+} is given in Eq. (9) and *ρ* is the distance along this ray. The ray amplitude is

where *J* = ∂(*x*,*z*)/∂(*ρ*,*x*′) is the Jacobian for the transformation (*ρ*,*x*′) → (*x*,*z*) evaluated at the observation point **r**, and *J*_{0} is the Jacobian at the exit point. We also identify *ρ _{c}* = − cos

^{2}

*θ*/∂

^{2}

_{x′}

*ψ*

^{+}

_{0}as the distance from

*x*′ to the caustic (note that ∂

^{2}

_{x′}

*ψ*

^{+}

_{0}< 0 so that

*ρ*> 0). Eq. (14) fails near the caustic at

_{c}*ρ*≈

*ρ*. For

_{c}*ρ*>

*ρ*, the square root in Eq. (14) becomes −

_{c}*i*∣1 −

*ρ*/

*ρ*∣−

_{c}^{1/2}. The caustic coordinates are obtained by setting

*ρ*→

*ρ*in Eq. (11) for each exit point

_{c}*x*′.

The exit points *x*′ of the rays that are associated with a given observation point **r** are found via

For observation points in the lit region of the caustic (region II in Fig. 1), this equation has two solutions, denoted as *x*′_{1,2} in accord with rays 1 and 2 discussed above (see Fig. 3). On the shadow side of SB^{+} (region III) only the *x*′_{1} solution exists, while the other solution *x*′_{2} is beyond the end point (i.e., *x*′_{2} > *x*_{ep}). In the paraxial region one may replace in Eq. (15) *ρ* by *z*, leading to a simple algebraic expression for *x*′_{1,2}.

As an illustrative example we depict in Fig. 4 the intensity of the initial field distribution at *z* = 0 and the exit points corresponding to certain observation points at the observation plane *z* = 10*z _{F}*. The observation points are tagged by their location in the 1st, 2nd, or 3rd lobes of the field at this observation plane. One observes, for example, that the exit points

*x*′

_{1}and

*x*′

_{2}corresponding to an observation point in the first lobe at

*z*= 10

*z*are located at the 44th and the 13th lobes of the initial field distribution, respectively.

_{F}The field on the lit side of the caustic is the sum of the two ray fields *U*_{GO1}^{+} + *U*_{GO2}^{+} [recall that ray 2 is beyond the caustic hence its amplitude function *A*^{+}_{2} contains a factor of -*i* as discussed after Eq. (14)]. The representation fails near the caustic at *ρ* ≈ *ρ _{c}* where

*A*

^{+}

_{1,2}in Eq. (14)explode. A uniform expression for the field, valid through the caustic, is obtained by replacing the GO fields, derived above, by the uniform geometrical optics (UGO) field [12,13]

where Ai′ is the derivative of the Ai function with respect to its argument. The parameters in Eq. (16) are functions of **r** and defined in terms of the GO parameters via

They are denoted, respectively, as the average and differential phases and the average and differential amplitudes. In the lit region, *σ* is real and positive and *σ*^{1/4} is taken to be the real and positive root. In the shadow region near the caustic, *σ* is real and negative and *σ*^{1/4} is taken such that its argument is *π*/4.

Equation (16) describes the field on both sides of the caustic. Far from the caustic on the lit side it reduces uniformly to *U*_{GO1}^{+} + *U*_{GO2}^{+} while on the shadow side it describes the decaying field there (see numerical results in Sec. 5). Note also that the Ai′ term can be neglected on the caustic.

#### 4.2. The uniform diffracted field

As mentioned earlier, the artificial end point *x* = *x*_{ep} gives rise to edge diffraction fields. They have the GTD form

where the ± sign corresponds to *ϕ*_{0}^{±}, *ρ _{d}*(

**r**) = [(

*x*−

*x*

_{ep})

^{2}+

*z*

^{2}]

^{1/2}is the distance from the end point at

*x*=

*x*

_{ep}to the observation point

**r**(see Fig. 3),

*ψ*

_{0d}

^{±}and

*A*

_{0d}

^{±}are the initial phase and amplitude functions at the end point, obtained by substituting

*x*′ =

*x*

_{ep}in Eq. (9).

*G*(

*ρ*) ~ exp(

_{d}*ikρ*+

_{d}*iπ*/4)/(8

*πkρ*)

_{d}^{1/2}is the asymptotic expression for the Green function, and D

^{±}are the diffraction coefficients

where *θ _{d}*(

**r**) = sin

^{−1}[(

*x*−

*x*

_{ep})/

*ρ*] is the exit angle of the ray from

_{d}*x*

_{ep}to

**r**and

*θ*

_{SB}

^{±}=

*θ*

^{±}(

*x*′) are the exit angles of rays SB

^{±}from

*x*

_{ep}(see Eq. (10) and Fig. 3; note that

*θ*

_{SB}

^{−}= −

*θ*

_{SB}

^{+}). The expression in Eq. (20) is the GTD diffraction coefficient under the physical optics approximation for a Neumann boundary condition at the edge (the Neumann boundary condition ∂

_{n}

*U*= 0 is implied by the Rayleigh-Sommerfeld formulation used in the the Kirchhoff integral in Eq. 4).

The edge diffraction fields in Eq. (19) explode near the shadow boundaries where *θ _{d}* ≈

*θ*

_{SB}

^{±}. A uniform expression for the fields, valid near the shadow boundary, is obtained via the uniform asymptotic theory (UAT) [14,15]. Henceforth we consider only the field near SB

^{+}since, as noted earlier, SB

^{−}is far from the region of interest (see Figs. 1 and 3). Recalling the discussion following Eq. (11), the ray that is shadowed by SB

^{+}is ray 2 in species

*ϕ*

_{0}

^{+}, hence the UAT replaces the fields

*U*

_{d}^{+}+

*U*

_{GO2}

^{+}. It is given by [14,15]

where *U*_{GO2}^{+} is given in Eqs. (12)-(14), *U _{d}*

^{+}is given in Eqs. (19)-(20), and

*F* is the Fresnel integral, *F*̂ is the asymptotic end point contribution in the Fresnel integral, sgn(*x*) = ±1 for *x* ≷ 0, *θ*_{SB}^{+} is defined after Eq. (20), *θ*_{2} is the exit angle of ray 2 as defined in Eq. (15), *ψ _{d}*

^{+}=

*ρ*+

_{d}*ψ*

_{0d}

^{+}is the phase of the diffracted ray [see Eq. (19)]. Note that

*ψ*

_{2}

^{+}>

*ψ*

_{d}^{+}hence the square root in Eq. (23) is taken to be positive.

One may show that the singularity of *U _{d}*

^{+}and of

*F*̂ near SB

^{+}(i.e.,

*η*≈ 0) are canceled, leaving the field there described by the

*F*term. For observation points far enough from SB

^{+}, the large

*η*asymptotics

*F*(

*η*) ~

*H*(-

*η*)+

*F*̂(

*η*), where

*H*(

*x*) = 0 or 1 for

*x*≶ 0, respectively, reduces Eq. (21) to

*U*

_{UAT}

^{+}~

*U*

_{GO2}

^{+}+

*U*

_{d}^{+}or to

*U*

_{d}^{+}in the lit or shadow sides of SB

^{+}, respectively.

#### 4.3. Asymptotic evaluation of the Kirchhoff integral

The results presented in the previous sections have been derived directly via the rules of the UGO and the UAT. Here they will be re-derived by asymptotic evaluation of the kirchhoff integrals in Eq. (5). We present only an outline of the analysis since the final expressions have been given above. For further analytical details the readers are referred to [16, Ch. 4].

Substituting the decomposition *ϕ*_{0} = *ϕ*_{0}^{+}+*ϕ*_{0}^{−} of Eq. (8) into the kirchhoff integral for *U*_{2} in Eq. (5), it is expressed as a sum of two generalized ray integrals *U*_{2} = *U*_{2}^{+}+*U*_{2}^{−}, with

where we have used the asymptotic expression for ∂_{n′}
*G* presented after Eq. (4). *A*^{±}_{0} and *ψ*_{0}^{±} are defined in Eq. (9).

We start with the evaluation of *U*_{2}^{+} which, as discussed above, describes the main contribution near the caustic. We utilize the complex steepest descent path (SDP) method. The staddle points are found by solving the equation ∂_{x′}
*ψ*^{+} (*x*′) = 0. In view of Eq. (25), this equation is identical to Eq. (15). There are, therefore, two real saddle points *x*′_{s1,2} which are the same as the exit points *x*′_{s1,2} found in Sec. 4.1. The saddle point contributions are given by

Noting that ∂_{x′}^{2}*ψ*^{+}(*x*′) ≷ 0 at *x*′_{1,2}, respectively, it follows that the sign of the square root in Eq. (26) is exp(±*iπ*/4), so that in the total amplitude in Eq. (26), *A*_{0}^{+} is multiplied by a real positive factor for *U*_{s1}^{+}, and by a negative imaginary factor for *U*_{s2}^{+}.

The expressions in Eq. (26) are identical with *U*_{GO1,2}^{+} discussed in Eq. (12)–(15). They fail near the caustic where the two saddle points coincide to a 2nd order saddle point with ∂_{x′}^{2}*ψ*^{+} = 0. A uniform expression, valid near the caustic, is derived via the Chester-Friedman-Ursell method [16, Sec. 4.5], [17]. The method utilizes a mapping of *U*^{+} to a canonical integral whose phase is a 3rd order polynomial that can be matched to both saddle points *x*_{s1,2}, simultaneously. The canonical integral can be expressed in terms of the Airy function and its derivative and it depends on the phases and amplitudes at the stationary points. Applying the procedure of [16, Sec. 4.5] to the integral in Eq. (24), using the the phases and amplitudes from Eq. (26), one obtains the UGO result in Eq. (16).

The asymptotic evaluation of *U*^{±} in Eq. (24) also yields the end-point contributions

These expressions are identical to the GTD fields *U _{d}*

^{±}in Eq. (19).

*U*

^{+}

_{ep}fails near the shadow boundary SB

^{+}where ∂

_{x′}

*ψ*

^{+}→0. A uniform asymptotic expression valid through this region is obtained by following the procedure in [16, Sec. 4.6] for a stationary point near an end point. The result can be expressed in the UAT form of Eq. (21).

## 5. Numerical results

Summarizing the results of the previous sections, the total radiated field is given by

where it is assumed here that the observation point is in region II of Fig. 1 (lit side of the caustic but not too close to the caustic or to SB^{+}). Near the caustic *U*_{GO1}^{+}+*U*_{GO2} are replaced by *U*_{UGO}^{+} as in Eq. (16). Near SB^{+}, *U*_{GO2}^{+}+*U _{d}*

^{+}are replaced by

*U*

_{UAT}

^{+}as in Eq. (21).

The separate contribution of *U*_{1} and *U*_{2} as well as their sum are depicted in Fig. (5), together with the paraxial closed form solution. The sum *U*_{1} + *U*_{2} agrees well with the closed form solution except for a small 1% error around the peak of the main lobe. This error is due to a small error in *U _{d}*

^{±}arising from the approximation of the initial field in Eq. (8) near the end point

*x*

_{ep}. Note that the first lobes are generated mainly by

*U*

_{2}and that

*U*

_{1}is negligible there.

As noted in Sec. 4, the artificial separation of the aperture field distribution into two parts at *x* ≶ *x*_{ep} gives rise to a spurious end-point effect in the numerical calculation of *U*_{1}. This contribution is canceled by the asymptotic end point diffractions *U _{d}*

^{+}+

*U*

_{d}^{−}in Eq. (28). Since the contribution of the main lobe in

*U*

_{1}is otherwise negligible, we obtain from Eq. (28) only the contribution of

*U*

_{GO1}

^{+}and of

*U*

_{GO2}

^{+}, which, as mentioned above, may be replaced by the UGO field, i.e.

Figure 6 compares the paraxial closed form solution of Eq. (3) to the UGO solution. An excellent agreement between the two is seen.

## 6. Discussion and conclusions

The Airy beam belongs to a new class of *curving beam* solutions with intriguing features such as *curved* propagation trajectories in free-space, *weak diffraction*, and *self healing*. Our goal in this paper has been to provide a cogent physical description to the Airy beam solutions viz a viz what we consider to be the basic definition of a “beam,” namely a local *wave function, whose evolution along its propagation trajectory is described by a local wave-dynamics*.

In order to establish the wave mechanism of the Airy beam we explored its properties within the rigorous framework of the Kirchhoff-Huygens integral, and also provided an asymptotically exact analysis utilizing the method of uniform geometrical optics (UGO). In Sec. 3 we calculated the field *U*_{1} radiated only by the main lobe of the Airy function initial field distribution, and demonstrated that it diffracts faster than a Gaussian beam with the same effective width (Fig. 2). The field *U*_{2} radiated by the other part of the initial field distribution has been studied in Sec. 4 using a uniform asymptotic evaluation of the corresponding Kirchhoff-Huygens integral. It was found that the Airy beam is described by sideways radiating rays that focus onto a caustic along the curved beam trajectory, and thereafter diverge away from it. It has been shown in Fig. 6 that the field of these focusing rays, expressed via the UGO, fully agrees with the Airy beam solution, thus providing a simple physical interpretation for its properties. It has thus been established that the evolution of the main lobe of the Airy beam along the curved trajectory is not described by a local wave dynamics, and it cannot be regarded as a “beam field” in the sense discussed above.

In addition to these focusing rays, *U*_{2} also consists of spurious contributions of end point diffraction, due to the “artificial” aperture truncation after it has been separated into two parts: the main lobe part and the rest of the aperture. These spurious contributions which are calculated asymptotically via a uniform version of the GTD, are canceled out by the numerical truncation effect in *U*_{1}, as has been discussed in Sec. 5, leaving the field described only by the UGO field.

The ray representation readily explains the “self-healing” property of the Airy beams reported in [10], i.e., its regeneration when the main lobe is obstructed. Noting that the main lobe field is not generated by a local field evolution but rather by sideways ray contributions, it follows that obstructing the main lobe at some point will have no effect on the forming of the main lobe at some range further on.

The ray interpretation can also be used to synthesize curving beam solutions that follow any prescribed (convex) trajectory, by synthesizing the aperture field distribution that gives rise to a caustic along this trajectory (for a recent physical realization see [18]). This is done by backtracking rays that emerge tangentially from that caustic to the aperture plane *z* = 0. The angle of these rays with respect to the aperture plane define the local gradient of the phase distribution there (see Eq. (10)). This concept can readily be generalized to the synthesis of curving beams in non-uniform media.

The concepts above may also be used to synthesize finite-energy Airy beam solutions. As discussed after Eq. (3), such beam fields are obtained by multiplying the Airy function aperture field distribution in Eq. (2) by a proper window function *F*(*x*). In view of the ray interpretation above, *F*(*x*) must be wide enough on the scale of *x*_{0} to incorporate many lobes of the Airy function aperture field distribution, and hence it can be regarded as a slowly varying amplitude function in the asymptotic evaluation of the Kirchhoff-Huygens integral. The resulting field is the same as the UGO result in Eq (16) except that *A*^{+}_{1,2} in Eq. (14) and in Eq. (25) are now multiplied by *F*(*x*′_{1,2}), the value of *F* at the ray exit points *x*′_{1,2}.

Finally, it should be noted that the closed form Airy function solution in Eq. (3) is valid only within the paraxial zone where the parabolic wave equation in Eq. (1) is valid. The UGO solution, on the other hand, is not restricted by the paraxial approximation and hence it is valid beyond the paraxial zone. One may show that the error in the paraxial solution in Eq. (3) is mainly due to the phase, causing a transversal shift of the solution. An estimate of the region of validity of the paraxial approximation can be obtained, therefore, by calculating the transversal shift of the paraxial caustic *x*̃ = *z*̃^{2}/4 with respect to the exact caustic, obtained via the ray procedure outlined in Sec. 4.1.

## Acknowledgement

This work is supported in part by the Israeli Science Foundation under Grant No. 674/07.

## References and links

**1. **G. A. Siviloglou, J. Brokly, A. Dogariu, and D. N. Christodoulides, “Observation of accelerating Airy beams,” Phys. Rev. Lett. **99** (2007). [CrossRef]

**2. **G. A. Siviloglou and D. N. Christodoulides, “Accelerating finite energy Airy beams,” Opt. Lett. **32**, 979–981 (2007). [CrossRef] [PubMed]

**3. **M. A. Bandres and J. C. Gutierrez-Vega, “Airy-Gauss beams and their transformation by paraxial optical systems,” Opt. Express **15**, 16719–16728 (2007). [CrossRef] [PubMed]

**4. **I. M. Besieris and A. M. Shaarawi, “A note on an accelerating finite energy Airy beam,” Opt. Lett. **32**, 2447–2449 (2007). [CrossRef] [PubMed]

**5. **M. A. Bandres, “Observation of accelerating parabolic beams,” Opt. Express **16**, 12866–12871 (2008). [CrossRef] [PubMed]

**6. **M. A. Bandres, “Accelarating parabolic beams,” Opt. Lett. **33**, 1678–1680 (2008). [CrossRef] [PubMed]

**7. **M. A. Bandres, “Accelerating beams,” Opt. Lett. **34**, 3791–3793 (2009). [CrossRef] [PubMed]

**8. **M. V. Berry and N. L. Balazs, “Nonspreading wave packets,” Am. J. Phys. **47**, 264–267 (1979). [CrossRef]

**9. **G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Ballistic dynamics of Airy beams,” Opt. Lett. **33**, 207–209 (2008). [CrossRef] [PubMed]

**10. **J. Broky, G. A. Siviloglou, A. Dogariu, and D. N. Christodoulides, “Self-healing properties of optical Airy beams,” Opt. Express **16**, 12880–12891 (2008). [CrossRef] [PubMed]

**11. **E. Heyman and L. B. Felsen, “Gaussian beam and pulsed beam dynamics: complex-source and complex-spectrum formulations within and beyond paraxial asymptotics,” J. Opt. Soc. Am. **18**1588–1611 (2001). [CrossRef]

**12. **D. Ludwig, “Wave propagation near a smooth caustic,” Bull. Amer. Math. Soc. **71**, 776–779 (1965). [CrossRef]

**13. **Yu. A. Kravtsov, “A modification of the geometrical optics method,” Izv. VUZ Radiofiz. Engl. transl., Radiophys. Quantum Electron. **7**, 664–673 (1964).

**14. **R. M. Lewis and J. Boersma, “Uniform asymptotic theory of edge diffraction,” Math. Phys. **10**, 2291–2305 (1969). [CrossRef]

**15. **S. W. Lee and G.A. Deschamps, “A uniform asymptotic theory of electromagnetic diffraction by a curved wedge,” IEEE Trans. Antennas Propag. **AP-24**, 25–34 (1976).

**16. **L. B. Felsen and N. Marcuvitz, *Radiation and Scattering of Waves* (Prentice-Hall, New Jersey, 1973).

**17. **C. Chester, B. Friedman, and F. Ursell, “An extension of the method of steepest descenets,” Proc. of Cambridge Phi.Soc. **53**, 599–611 (1957). [CrossRef]

**18. **T. Ellenbogen, N. Voloch-Bloch, A. Ganany-Padowicz, and A. Arie “Nonlinear generation and manipulation of Airy beams,” Nature Photonics **3**, 395–398 (2009). [CrossRef]