Fox H -Functions in Self-Consistent Description of a Free-Electron Laser

: A fractional calculus concept is considered in the framework of a Volterra type integro-differential equation, which is employed for the self-consistent description of the high-gain free-electron laser (FEL). It is shown that the Fox H -function is the Laplace image of the kernel of the integro-differential equation, which is also known as a fractional FEL equation with Caputo–Fabrizio type fractional derivative. Asymptotic solutions of the equation are analyzed as well.


Introduction
In this paper, we discuss a fractional calculus concept for the classical electrodynamics of free-electron lasers (FEL)s. It is well known that a self-consistent description of FELs is presented in the framework of an integro-differential equation. The latter can also be considered as a specific form of a fractional integro-differential equation. We study the kernel of this equation and show that it can be presented in the form of the Fox H-function. The Fox H-functions is widely used in fractional calculus, and it plays an essential role in a variety of applications of fractional calculus [1,2]. Between many of these examples, a new one is in the field of fractional electrodynamics of the FEL, recently considered in ref. [3]. Experimental implementation and theoretical description of the FEL is a long-lasting problem that started in the seventies of the last century. This extensively studied phenomenon is well described and reviewed [4][5][6][7], to mention a few. Contemporary studies are also reflected in recent publications and related to both classical and quantum descriptions [8,9]. In a paradoxical way, the classical electrodynamics of electrons interacting with electromagnetic fields explains this quantum lasing phenomenon [4][5][6][7]10]. In particular, the self-consistent evolution in the small-signal slow-varying amplitude approximation (of the electromagnetic field) is accounted by the Volterra-type integro-differential equation [3].
Here τ = (t + z/c)/∆t is a dimensionless gain time-scale, where ∆t is the interaction time, z is the longitudinal coordinates, t is the current time, and c is the light speed. The dimensionless parameters of the system include the resonance parameter v, which is linked to the laser frequency and scaled by ∆t, the small-signal gain coefficient g 0 , and the coefficient µ, which relates to a parameter regulating the effects of the gain reduction due to the electron's energy distribution [3].
In the absence of the attenuation of the gained signal, µ = 0, this integro-differential Equation (1) has a solution in the superposition form E(τ) = ∑ j E j exp(iΛ j τ), where E j and Λ j are related to the roots of a cubic equation [10]. Performing the Laplace transform of Equation (1) we have for µ = 0Ẽ where E 0 is the initial condition for Equation (1). The poles s j are the roots of the cubic equation determined by the Cardano rule. This equation defines both Λ j and E j . In the Volterra Equation (1), studied in ref. [3], the integro-differential operator in the r.h.s. was considered by analogy to the Caputo-Fabrizio fractional derivative [11], which is a fractional derivative without a singular kernel [3,11]. This fractional derivative/operator has been treated in the form of an iteration technique, based on an expansion employing a family of two variable Hermite polynomials that eventually leads to the analytical solution [3]. Following this fractional calculus concept of refs. [3,11], it is reasonable to suggest an alternative approach for the kernel of the integral operator, presenting it in the form of the Fox H-functions. In this case, the Laplace transformation of Equation (1) becomes feasible for µ = 0.

Fox H-Function in Laplace Space
Performing the Laplace transformation of Equation (1), one obtains whereG(s) is defined by the integral The way of introducing the Fox H-function inside the integrand, based on the representation of the exponential function in the form of the Fox H-function by means of the Mellin-Barnes integration see, e.g., [2,12,13]. Here Γ(ξ) is a gamma function: Γ(ξ + 1) = ξΓ(ξ). Performing the variable change y = τ 2 in the integral (5), we havẽ where s v = s + iv is used for brevitys sake. Then taking Z = s v y 1/2 and substituting the Mellin-Barnes integral (6) inside integration (7), we obtain the chain of transformations as follows Eventually, we arrive at the definition of the Fox H-function, which is presented in terms of the Mellin-Barnes integral. Its general definition reads [2] H m,n p,q (Z ) = H m,n p,q Z (a 1 , where with 0 ≤ n ≤ p, 1 ≤ m ≤ q and a i , b j ∈ C, while A i , B j ∈ R + , for i = 1 , . . . , p, and j = 1 , . . . , q. The contour Ω starting at c − i∞ and ending at c + i∞, separates the poles of the functions Γ(b j + B j ξ), j = 1 , . . . , m from those of the function Γ(1 − a i − A i ξ), i = 1 , . . . , n.

Series Expansion and Asymptotics
The Fox H-function can be presented in the form of series expansion [2]. Then we have the l.h.s. of Equation (14) Let us consider asymptotic behavior, when |s/v| 1 and |vµ −1/2 | 1. Then to obtain a gained signal at least the first four terms in the expansion should be accounted for. Then we have where Then, Equation (12) reduces to the cubic equation with roots s j = s j (µ, v).
The coefficients of the expansion in the form of the Fox H-functions of large arguments behave as follows [2] Then C r = µ/v 2 for r = 0, 1, 2, 3, and we have Thus, by analogy with Equation (12), the poles ofẼ(s) are determined by a cubic equation, which now reads Note that this expression is independent of µ, which results from the asymptotic consideration for both |s/v| 1 and |vµ −1/2 | 1. In any case of the cubic equation, the solution is the superposition of three waves where s j are roots of the cubic equation, which is an approximation of Equation (12) that also determine E j with the initial condition ∑ 3 j=1 E j = E 0 .

Asymptotics of Small τ and Series Expansion
Disregarding the fractional concept, related to the Fox H-function and looking for the asymptotic solution of Equation (1) for the small gain time-scale τ 1, a simplified consideration can be suggested as follows. Noting that the Laplace imageG(s) in Equation (5) is a table integral [14], which reads Erfc s + iv Taking into account the asymptotic behavior of the Erfc-function for the large values of |s + iv| 1 [15], we havẽ Therefore, the initial gain solution of Equation (1) is G(τ) = ∑ 3 j=1 E j e s j τ , where s j is the root of the cubic equation.

Conclusions
A fractional calculus concept was considered in the framework of a Volterra type integro-differential equation, which is known is employed for the self-consistent description of the high-gain free-electron laser (FEL). We have shown that the Fox H-function can be employed for the Laplace image of the kernel of the integro-differential equation. The analysis was performed in the framework of the Laplace transformation with respect to the gain time-scale τ. Note that the FEL geometry can be chosen in such a way that τ > 0. This approach makes it possible to obtain an exact analytical expression for the Laplace image of the gained signalẼ(s), and its singular behavior is determined by the roots of the transcendent Equation (12). Further analytical analysis is possible (and presented) in the asymptotic approximation for both large τ 1 and small τ 1. In either case, these solutions are described by cubic equations with coefficients depending on µ. It is worth mentioning that an alternating approach to the electron energy distribution has been considered as well [16] in the framework of the fractional generalization of the FEL Equation (1) with µ = 0.
Discussing a mathematical aspects related to fractional calculus, it should be admitted that, the r.h.s. of Equation (1) can be considered as an FEL fractional derivative by analogy with the Caputo-Fabrizio fractional derivative [11,17], which reads as follows where m(α) is a normalization term with constant α. Therefore the FEL fractional derivative reads Therefore, the Laplace image of the kernel of the FEL fractional derivative,G(s) obtained in Equation (11) is the Fox H-function in Laplace space. We also note in passing that D FEL in Equation (26) differs from those introduced in ref. [3] by the term (τ − ξ) f (ξ) → ξ f (ξ).

Conflicts of Interest:
The author declares no conflict of interest.