Extensive studies of self-accelerating beams have been made recently. The basic such beam is the Airy solution
governed by the paraxial equation
in free space. Here, and are, respectively, dimensionless transverse and longitudinal variables, defined in terms of the original variables X and Z, the wavenumber k, and a scaling factor with units of meters. The complexification , where is a positive parameter, ensures the square integrability (finite energy) of the input function and, hence, of for . The finite-energy version of the solution given in Equation (1) was first formulated analytically by Siviloglou and Christodoulides  and subsequently demonstrated experimentally by Siviloglou et al.  . Their work was motivated by the infinite-energy (nonspreading) accelerating Airy solution to the Schrödinger equation introduced by Berry and Balazs  in the context of quantum mechanics. A full wave theoretical analysis of the Airy beam has been undertaken by Kaganovsky and Heyman  . An Airy beam is slowly diffracting; it can retain its intensity over several diffraction lengths while bending laterally along a parabolic path even though its centroid is constant. Another feature, which has been demonstrated both analytically and experimentally, is that an Airy beam propagating in free space can perform ballistic dynamics akin to those of projectiles moving under the action of gravity.
Both bending Airy beams and accelerating Airy wavepackets are characterized by self-healing properties; they tend to reform in spite of the severity of imposed perturbations; this is due to the reinforcement of the main lobe by the side lobes. The robustness of such beams has been studied in the presence of material dispersion  , deterministic inhomogeneities (see  and references therein), and turbulent environments   . These exotic properties suggest various physical applications, such as Airy beam-mediated particle cleaning and vacuum electron acceleration. A recent review of the theory, generation and applications of Airy beams has been published by Hu et al.  .
2. Critical Comments on Accelerating Olver Beams
Belafhal et al.  have reported a new class of “accelerating” (more precisely, nonlinearly self-bending) nondiffracting solutions of the free-space parabolic equation [cf. Equation (2)] in terms of the Olver functions
which, in turn, are related to the solutions of the ordinary differential equation
Specifically, the authors state that a novel class of nondiffracting solutions to Equation (2) is given as follows:
The first equality above is correct. The integral expression does, indeed, satisfy the parabolic Equation (2) for all values of m, as the authors have shown in Appendix B of Ref.  . However, the second equality above is incorrect. It holds only for , which yields the well-known Airy beam solution given in Equation (1). This affects the statements in Equations (3) and (4), Equation (A-1) and Equations (B-6) and (B-7) in Ref.  . A general proof of our assertion can be provided by using the ansatz
into the parabolic Equation (2). As a result, one obtains the ordinary differential equation
a solution of which is the Airy function; specifically, . This is in contradiction to the second equality in Equation (5) in which the Olver functions are related to the solutions to Equation (4) for .
The authors give in Equation (24) of Ref.  the spectra of exponentially apodized initial conditions as follows:
for m = 0, 1 and 2. Here, denotes the Fourier transform of the Olver function . For , one obtains, then,
The wave function is the well-known Airy beam and contains finite energy for a positive parameter a. The Fourier representations of the wavefunctions given by the integral expressions above do not yield analytical solutions for . The expressions for m = 1 and 2 could justifiably be called Olver beams (but not necessarily accelerating Olver beams) because they are associated with apodized fourth-order diffusion (super diffusion) and fifth-order Airy (hyper Airy) functions, respectively, at . However, a careful examination is needed in order to establish whether the exponential apodization is sufficient for ensuring finite energy. It turns out that the Airy beam contains finite energy only for . The choice of appropriate apodization for requires careful examination.
Due to their complexity, the wave functions for m = 1 and 2 can be examined only numerically. The authors allude to such a program in Sec. 4, where they propose hologram masks. However, it is not clear whether the integrations given in Equation (9) above, which are required for , have been carried out, and what is the basis for the plots in Figure 3 of Ref.  . Why do the plots for the fourth-order and sixth-order super diffusion, characterized essentially by single lobes, appear multi-lobed under finite-energy (apodised) conditions?
Consider the second expression in Equation (5), viz.,
It has been shown in this note that this does not represent a solution of the parabolic Equation (2), unless . Is it possible, however, that there exists another type of equation the solutions of which can be expressed as in Equation (10) for all values of m? Such an equation does exist and it is given as follows:
Introducing the ansatz given in Equation (6) into Equation (11), one obtains the ordinary differential equation
which has been studied extensively by Olver  . The solutions are related to the Olver functions; for example,
Since Equation (12) above is identical to Equation (1) in Ref.  , with given in Equation (6), the specific result in Equation (13) above indicates that Equation (6) in  is invalid.
3. Critical Comments on Accelerating Olver-Gauss Beams
Hennani et al.  have used the Fresnel diffraction formula for a paraxial ABCD optical system, viz.,
in order to determine the output due to the input Olver-Gauss function
Here, is a normalization parameter with units of length, and are positive dimensionless parameters used to ensure finite energy. The main result is given in Equation (7) of Ref.  , viz.,
Unfortunately, this result is incorrect, in general, because it is based on a variation of the expression given in Equation (5); specifically,
This formula and, as consequence, Equation (16), is valid only for . In this case, the solution in Equation (16) corresponds to a finite-energy accelerating Airy-Gauss beam which has been studied previously (see, e.g., Ez-Zariy et al.  ). All the other cases examined analytically and numerically for different ABCD parameters in Ref.  are incorrect.
The free-space version of Equation (16) corresponding to the optical ABCD matrix
with z the direction of propagation, has been used by Hennani et al.  in order to study the action of the radiation forces produced by a highly focused finite-energy Olver-Gaussian beam on a Rayleigh dielectric sphere. For the reasons detailed above, the numerical results presented by the authors are valid only for the case .
4. Concluding Remarks
In the seminal work of Berry and Balazs  it was determined that the infinite-energy Airy wave packet is the only self-accelerating nondiffracting solution to the (1 + 1)D Schrödinger equation. Within the framework of optics, Siviloglou and Christodoulides  rendered the Airy beam given in Equation (1) physically realizable by appropriate apodization of the initial wave function, viz., , with a positive parameter .
The introduction by Belafhal et al.  of the Olver beams as a novel class of self-accelerating nondiffracting solutions to the paraxial Equation (2) is a direct contradiction to the work of Berry and Balazs. It has been shown in this note that the work in  is valid only for , which coincides with the Airy solution. Unfortunately, the mistakes in  have been carried over to the work by Hennani et al.  on the propagation of Olver-Gauss functions through an optical ABCD system, and by Hennani et al.  on the action of the radiation forces produced by a highly focused finite-energy Olver-Gaussian beam on a Rayleigh dielectric sphere.
In closing, it should be pointed out that the Olver functions
can be considered as “incomplete” elementary catastrophe integrals; the latter are defined as
is precisely equal to and corresponds to the fold catastrophe. for n = 4, 5 and 6 are the cusp, swallowtail and butterfly catastrophes, respectively. The dynamics of the cusp, swallowtail and the butterfly optical catastrophes have been studied in terms of solutions to the 3D paraxial equation    . There exist more complicated catastrophe integrals than those in Equation (20); for example, the hyperbolic umbilic. Self-accelerating solutions of the 3D paraxial equation associated with this catastrophe have been determined recently   .
 Siviloglou, G.A., Broky, J., Dogariu, A. and Christodoulides, D.N. (2007) Observation of Accelerating Airy Beams. Physical Review Letters, 99, 213901.
 Besieris, I.M. and Shaarawi, A.M. (2008) Accelerating Airy Wave Packets in the Presence of Quadratic and Cubic Dispersion. Physical Review E, 78, 046005.
 Besieris, I.M., Shaarawi, A.M. and Zamboni-Rached, M. (2016) Accelerating Airy Beams in the Presence of Inhomogeneities. Optics Communications, 369, 56-64.
 Hu, Y., Siviloglou, G.A., Zhang, P., Efremidis, N.K., Christodoulides, D.N. and Chen, Z. (2012) Self-Accelerating Airy Beams: Generation, Control and Applications. In: Chen, Z. and Marandotti, R., Eds., Nonlinear Photonics and Novel Optical Phenomena, Springer Series in Optical Sciences 170, New York, 1-44.
 Belafhal, A., Ez-Zariy, L., Hennani, S. and Nebdi, H. (2015) Theoretical Introduction and Generation Method of Novel Nondiffracting Waves: Olver Beams. Optics and Photonics Journal, 5, 234-246.
 Hennani, S., Ez-Zariy, L. and Belafhal, A. (2015) Propagation of Finite Olver-Gaussian Beams Passing Through a Paraxial ABCD Optical System. Optics and Photonics Journal, 5, 273-294.
 Ez-Zariy, L., Hennani, S., Nebdi, H. and Belafhal, A. (2014) Airy-Gaussian Beams Passing Through a Misaligned Optical System with Finite Aperture. Optics and Photonics Journal, 4, 325-336.
 Hennani, S., Ez-Zariy, L. and Belafhal, A. (2015) Radiation Forces on a Dielectric Sphere Produced by Finite Olver-Gaussian Beams. Optics and Photonics Journal, 5, 344-353.
 Ren, Z.J., Ying, C.F., Jin, H.Z. and Chen, B. (2015) Generation of a Family of Pearcy Beams Based on Fresnel Diffraction Catastrophes. Journal of Optics, 17, 105608.
 Zannotti, A., Diebel, F., Bogulawski, M. and Denz, C. (2017) Optical Catastrophes and the Swallowtail and Butterfly Beams. New Journal of Physics, 19, 053004.