Theory of generalized Bessel functions

Generalized revival and splitting of an arbitrary optical field in GRIN media

Assuming a non-paraxial propagation operator, we study the propagation of an electromagnetic field with an arbitrary initial condition in a quadratic GRIN medium. We show analytically that at certain specific periodic distances, the propagated field is given by the fractional Fourier transform of a superposition of the initial field and of a reflected version of it. We also prove that for particular wavelengths, there is a revival and a splitting of the initial field. We apply this results, first to an initial field given by a Bessel function and show that it splits into two generalized Bessel functions, and second, to an Airy function. In both cases our results are compared with the numerical ones.

Discrete-like diffraction dynamics in free space

We introduce a new class of paraxial optical beams exhibiting discrete-like diffraction patterns reminiscent to those observed in periodic evanescently coupled waveguide lattices. It is demonstrated that such paraxial beams are analytically described in terms of generalized Bessel functions. Such effects are elucidated via pertinent examples.

Recursive algorithm for arrays of generalized Bessel functions: Numerical access to Dirac-Volkov solutions

In the relativistic and the nonrelativistic theoretical treatment of moderate and high-power laser-matter interaction, the generalized Bessel function occurs naturally when a Schrödinger-Volkov and Dirac-Volkov solution is expanded into plane waves. For the evaluation of cross sections of quantum electrodynamic processes in a linearly polarized laser field, it is often necessary to evaluate large arrays of generalized Bessel functions, of arbitrary index but with fixed arguments. We show that the generalized Bessel function can be evaluated, in a numerically stable way, by utilizing a recurrence relation and a normalization condition only, without having to compute any initial value. We demonstrate the utility of the method by illustrating the quantum-classical correspondence of the Dirac-Volkov solutions via numerical calculations.