Theoretical Introduction and Generation Method of a Novel Nondiffracting Waves : Olver Beams

In this paper, we introduce a new class of scalar nondiffracting Helmholtz-equation solution. We demonstrate that this novel wave-equation solution has some specific orders; among these ordinary Airy beams which are regarded as the zeroth order. Moreover, a general expression of these novel beams, which are named Olver Beams and referred to OBs, is developed. The zeroth and the first high orders of the incident OBs are presented theoretically and numerically in this paper. Yet, based on a computer generated holograms method, the generation’s masks of the Finite OBs in first orders are given in this work. Also, the incident transverse intensity distribution in 1-D and 2-D of the first orders of OBs is performed.


Introduction
Olver beams which are the subject of this work, are proved to be a new class of nondiffracting beams and therefore as solution of scalar nondiffracting Helmholtz-equation.This new beams family is defined for the non-negative integer index beam's order m.The zeroth order of this class of beams has described the ordinary Airy beams, which were introduced firstly by Berry et al. [1] in 1979, in the context of quantum mechanics.In this last work, the authors were interested to explain all the propagation properties of the Airy wave packet.The latter is a solution of Schrödinger equation which is equivalent to the equation of the linear propagation under the paraxial approximation of diffraction.The Airy beams class had attracted interest since their introduction theoretically and experimentally by Siviloglu et al. [2] [3] in 2007.Recently, for physical realization, this beams family has been modulated by an exponential term to obtain the finite Airy beam which is expressed as a product of the classical Airy mode and the linear exponential function.For many reasons, researchers have tried to properly study these laser beams and to be interested in their most remarkable properties.These beams are already considered as a single type of nondiffracting beams that are currently described in literature as a one-dimensional system [2]- [10].The concept of Airy beams has been deeply explored on optical transforms for light beams; Airy transform can be applied in optics laser in many contexts by using Fourier Transform, Hankel Transform and Hilbert Transform.These transforms are all useful in beam shaping, beam analysis, image processing and others [11]- [22].
Since 2013, our research laboratory group has started its researches over the Airy class beams through different optical systems.At first, we found a theoretical formalism for the study of the propagation of these nondiffracting beams family through an ABCD optical system with a rectangular annular aperture [23].Another study has demonstrated the propagation characteristics of Airy-Gaussian beams passing through a misaligned ABCD optical system with aperture [24].Thereafter, our studies have gone further than the previous work to do the same study of the propagation through an optical ABCD system, but this time with a new form of Airy-Hermite-Gaussian beam [25].This series of studies continue to give an interesting work that explores another new form of Finite Airy-Gaussian-Hermite-Gaussian Beams propagating through a paraxial ABCD optical system [26], which is a special case of the Olver beams.
The present work will be an extension of the study of nondiffracting Airy beams.We will give a general analytical expression of a novel beams family, called Olver beams and is a solution of Schrödinger equation.This solution can be expressed in terms of Olver functions obtained directly by solving Olver's differential equation given by [27]: where ( ) f t is a real continuous function equal to a constant factor that times a power of t .The analysis of this function, which represents the solution of Equation ( 1), is the subject of the Olver beams that will be used in a large number of more general applications concerning their distributions through a different optical system.Therefore, the form of these solutions changes depending on the sign of t which will be discussed in the second part of this paper besides reporting the first idea of Olver optical beams.In the third part, we will give all the integral representations and analytical expressions of Olver functions noted as

( )
. m O of different orders m.Then we will introduce the finite Olver electrical field formulas and numerical representations of their corresponding intensity distributions for different orders m = 1, 2 and 3. Finally, in the fourth section of this report and by using a computer-generated holograms technique, we present the masks used in the experimental generation of the finite Olver beams.Also, in the same section we will draw all the masks of the finite Olver beams for each order m.Thanks to their alternative properties, below, we separate between four classes of Olver beams corresponding to different orders.

Olver's Equations and Their Solutions
The Airy function which is a special case of generalized Olver function at m = 0 is a solution of the scalar Helmholtz-equation and their generalization is also a solution of differential wave-equation.So, the Olver functions evolve according to the scalar wave-equation in terms of the variables ( ) whose solution is, if only if we consider that s is a purely imaginary complex, ( ) is the Olver function of order m, with 0,1, 2, m =  which would be calculated analytically for different orders.Equation ( 3) is verified directly by using the integral representation of the Olver function ex-pressed as follows (Appendix A & B): Our analysis suggests also a potential relationship which links the solutions of Olver's equation with the Olver functions.It is the starting point for the discovery of Olver beams family using a well-known form of Olver's ordinary differential equation in the unknown function ( ) W t , given by Equation (1), when ( ) where m denotes a positive integer as 2 m ≥ .Equation ( 1) adopts as one standard solution according to the sign of t denoted ( ) m V t .These solutions can be written in terms of Olver functions with the following equation [27]: Depending on the sign of the parameter t , three essentially cases arise.The classification of all the solutions is given by: Case 1: 0, t ∈ +∞ , the Olver's function, mathematically, is given by: ( ) ( ) . O t given by Equation ( 6), we will find: ( ) .

Case 2:
On the contrary, when the variable , the Olver's function yields to : where ( ) cot . is the cotangent function,

Case 3:
When the variable t is nul ( ) With the use of Equation ( 6), one obtains ( )

Introduction of Olver Beams as Nondiffracting Beams Family
In this part, we will include the exact analytical and illustrative graphs of Olver functions ( ) for different orders m .It is shown that these functions have the following integral representation: For the reason to obtain the accurate analytical expression of this integral established in the Equation (13) which will be calculated according to the parity of γ .So, four cases are considered.
Case 1: When the parameter γ is multiple of an integer ( ) , in this case, the Olver order m can take the values 1,5,9, m =  In this case,

( )
m O x will be given by the following integral representation After some calculations, we find giving: ) After developing this last equation, we find we find: ) In this case, it is pointed that the zeroth-order of Olver beams corresponds to the pur Airy beam [1].

Special Cases Representation
Based on the above formulas elaborated in Equations ( 15), ( 17), ( 19) and ( 21), the closed-forms analytical expressions of the first orders of Olver's nondiffracting beams are given by: ( ) In Figure 1, we present the one-dimensional distribution of Olver's function for different orders m, as m = 0, 1, 2 and 3 determined by Equation (22).
A comparison is conducted between the ordinary Airy beam representation [28] and our zeroth-Olver function calculated at m = 0 by using Equation (22).This result validates the theoretical concept used to calculate Olver functions ( ) m O x for any m order.Whereas in Figure 2, we plot the two-dimensional intensity distribution of the Olver beams in comparison with finite Olver beams.
We notice the appearance of some secondary lobes in the left of the central lobe for the even orders of the Olver beams.While we observe these sidelobes at the right of the central lobe for the odd orders.See Figures .  (a-b)

Olver's Holograms Masks
In this section, we will introduce how to generate the Olver beams in all orders by using a computer-generated holograms technique which is defined as a recording of the interference pattern between two light beams, making this technique an easily use full tool.The Fourier Transform (TF) provides a good understanding of Franhaufer diffraction for holograms.Hence, the result of this Franhaufer diffraction is proportional to the Fourier transform of the product of the input field (Gaussian beam) (E inp ) and of the Hologram transmission function noted H.By producing such holograms, the transformation of the interference patterns is written as [29] [30]: This transformation will be used to simulate the holographic masks illustrated in Figure 3, and expressed analytically for the Finite Olver beams for the first special orders as: For the experimental realisation, one can use the proposed masks to elaborate the desired Finite Olver beams.Note that the first mask corresponds to the one used in [2] [3].

Conclusion
In conclusion, a general analytical expression of the Olver beams (OBs) was devoted for various orders as a novel family class of conical beams and was considered as a generalization of the ordinary Airy beams.We had   demonstrated that the OBs were the solution of the scalar Helmholtz equation.Yet, in this work, the finite Olver beams are reported as pseudo-nondiffracting beams related to the OBs.We demonstrated the ability to generate and to produce this last class of optical beams using a computer-generated holograms technique.The obtained results provide more ideas for studying the propagation of the Olver beams and for examining their tuning properties through different optical systems.
Analytically, for m = 0 and by using the following form of integration given by [31]- [34]: ( ) we find ( ) For m = 0, the integral given by Equation (A.1) is theoretically valid.For its numerical evaluation, we first calculate numerically the integral given by Equation (A.1) using the numerical Trapeze integration method.Then we simulate numerically the second part given as follows: ( ) ( ) ( ) The integral given by Equation ( 1) is also validated numerically as shown in Figure A1.We note that the calculation parameters chosen are given by: 1 3 a = − , and s = 0.05.

Case 2:
For m = 1, the numerical evaluation of the integral expression given by Equation (A.1) is obtained by the same way as the first case, where m = 0. We calculate numerically the integral given by Equation (A.1) by using the numerical Trapeze integration method, and then we simulate numerically the second part of the integral given as follows: Figure A2 shows that our result is validated numerically by the following calculation parameters:

Case 3:
For m = 5, as it has been done before for m = 0 and 1, we calculate numerically the integral expression presented in Equation (A.1) by using the numerical Trapeze integration method.Then, we simulate numerically the second part of the integral given as follows: ( ) ( ) ( ) Figure A3 shows that the result of 5 I is validated numerically using the following calculation parameters: In this appendix A, we tried to validate the integral given by Equation (A.1) with a numerical calculation.This will be useful to approach many theoretical optical problems concerning Olver expressions.We are confident to say that our results are satisfying.

Appendix B
We shall prove in this Appendix that the integral representation of Olver function is a solution of Schrödinger equation for a particle with mass m, namely: and for different values of the beam order, the probability density which is sketched in Figure B1, is given by: Bessel function of a second kind.Using the relationship between function.Using Equation (6), we find: Olver order m stands to m = 2, 6, 10,••• .The integral representation of ( ) m O x in this case is:

Figure 1 .
Figure 1.Distribution of Olver's functions for various orders m.
Figure A2shows that our result is validated numerically by the following calculation parameters: 1 4 a = − and s = 0.05.

Figure B1 .
Figure B1.Probability density for the Olver wave packet with