1. Introduction
The analysis of fractional differential equations involving Hadamard fractional derivatives has increased interest in mathematical analysis, as proved, for example, by the publication of the recent monograph [
1]. On the other hand, few results regarding the applications of Hadamard fractional differential equations in mathematical physics (see, for example, [
2] and probability (see [
3]) exist. In [
4], some analytical results regarding Hadamard fractional equations with time-varying coefficients have been pointed out. In particular, the following
-Mittag–Leffler function was introduced:
The
-Mittag–Leffler function is an entire function of the complex variable
z if the parameters are such that
,
and
(see [
5]). This function was independently introduced and studied by Gerhold in [
6]. In [
4], we called this special function
-Mittag–Leffler function, since it includes for
the well-known two-parameters Mittag–Leffler function
widely used in the theory of fractional differential equations (see the recent monograph [
7] and the references therein). Moreover, the
-Mittag–Leffler function is a generalization of the so-called Le Roy function [
8]
In the more recent paper [
5], the authors studied the asymptotic behavior and numerical simulation of this new class of special functions.
We observe that the Le Roy functions are used in probability in the context of the studies of COM-Poisson distributions [
9], which are special classes of weighted Poisson distributions (see for example [
10]). We show in
Section 3, as a first probabilistic application, that the
-Mittag-Leffler functions can also be used in the construction of a new generalization of the COM-Poisson distribution that can be interesting for statistical applications and in physics in the context of generalized coherent states [
11].
The aim of this paper is to study some particular classes of Hadamard fractional integral equations and differential equations whose solutions can be written in terms of the
-Mittag–Leffler function (
1) and are somehow related to a fractional-type generalization of the COM-Poisson distribution.
A second application in probability is contained in
Section 4. We observe that Imoto [
12] has recently introduced the following generalization of the COM-Poisson distribution,
involving the normalizing function
with
,
and
or
,
and
. The special function (
5) is somehow related to the generalization of the Le Roy function, while the distribution (
4) includes the COM-Poisson for
and
,
. In
Section 4, we show that this function is related to integral equations with a time-varying coefficient involving Hadamard integrals.
In
Section 5, we present other results concerning the relation between Le Roy-type functions and Hadamard fractional differential equations. We introduce a wide class of integro-differential equations extending the hyper-Bessel equations. This is a new interesting approach in the context of the mathematical studies of fractional Bessel equations (we refer for example to the recent paper [
13] and the references therein).
In conclusion, the main aim of this paper is to establish a connection between some generalizations of the COM-Poisson distributions and integro-differential equations with time-varying coefficients involving Hadamard integrals or derivatives. As a by-product we suggest a possible application of the -Mittag–Leffler function to build a generalized COM-Poisson distribution that in future should be investigated in more detail.
2. Preliminaries About Fractional Hadamard Derivatives and Integrals
Starting from the seminal paper by Hadamard [
14], many papers have been devoted to the analysis of fractional operators with logarithmic kernels (we refer in particular to [
15]). In this section, we briefly recall the definitions and main properties of Hadamard fractional integrals and derivatives and their Caputo-like regularizations recently introduced in the literature.
Definition 1. Let and . The Hadamard fractional integral of order α, applied to the function , , , for , is defined as Before constructing the corresponding differential operator we must define the following space of functions.
Definition 2. Let be a finite interval such that and let be the space of absolutely continuous functions on . Let us denote and define the space Clearly for .
Definition 3. Let , and , where is the integer part of α. The Hadamard fractional derivative of order α applied to the function , , is defined as It has been proved (see e.g., Theorem 4.8 in [
16]) that in the
space,
,
, the Hadamard fractional derivative is the left-inverse operator to the Hadamard fractional integral, i.e.,
Analogously to the Caputo fractional calculus, the regularized Caputo-type Hadamard fractional derivative is defined in terms of the Hadamard fractional integral in the following way (see, for example, [
17])
where
,
and
, with
. In this paper, we will use the symbol
for the Caputo-type derivative in order to distinguish it from the Riemann-Liouville type definition (
8) and also to underline the fact that essentially it coincides with the fractional power of the operator
. Moreover, by definition, when
,
. The relationship between the Hadamard derivative (
8) and the regularized Caputo-type derivative is given by ([
17], Equation (
12))
In the sequel we will use the following useful equalities (that can be checked by simple calculations),
for
and
. It is immediate to see that
Finally, we observe that formally the relationship between Hadamard-type derivatives and Riemann-Liouville derivatives is given by the change of variable , leading to the logarithmic kernel. According to this change of variable, the Hadamard derivative represents the counterpart of the fractional power of the operator , i.e., the fractional power of the operator .
3. Generalized COM-Poisson Processes Involving the -Mittag–Leffler Function and the Related Hadamard Equations
Here we show a possible application of the
-Mittag–Leffler function (
1) in the context of fractional-type Poisson statistics and we discuss the relation with Hadamard equations.
We first recall that weighted Poisson distributions have been widely studied in statistics in order to consider over or under-dispersion with respect to the homogeneous Poisson process. The probability mass function of a weighted Poisson process
,
, is given by (see [
10])
where
N is a Poisson distributed random variable
is a non-negative weight function with non-zero, finite expectation, i.e.,
The so-called COM-Poisson distribution introduced by Conway and Maxwell [
9] in a queueing model belongs to this wide class of distributions. A random variable
is said to have a COM-Poisson distribution if
where
is the Le-Roy function.
This distribution can be viewed as a particular weighted Poisson distribution with
. The Conway-Maxwell-Poisson distribution is widely used in statistics in order to take into account over-dispersion (for
) or under-dispersion (for
) in data analysis. Moreover it represents a useful generalization of the Poisson distribution (that is obtained for
) for many applications (see for example [
18] and the references therein). Observe that when
and
it reduces to the geometric distribution, while for
it converges to the Bernoulli distribution.
Another class of interesting weighted Poisson distributions are the fractional Poisson-type distributions
where
is the Mittag–Leffler function (see the recent monograph [
7] for more details). In this case, the weight function is given by
. There is a recent wide literature about fractional generalizations of Poisson processes (see e.g., [
19,
20,
21,
22,
23,
24,
25,
26,
27]) and we should explain in which sense Equation (
19) can be regarded as a fractional-type Poisson distribution. The relationship of this kind of distributions with fractional calculus, first considered by Beghin and Orsingher in [
28], is given by the fact that the probability generating function
satisfies the fractional differential equation
where
denotes the Caputo fractional derivative w.r.t. the variable
u (see [
16] for the definition). Essentially, the connection between the distribution (
19) and fractional calculus is therefore given by the normalizing function, since the Mittag–Leffler plays the role of the exponential function in the theory of fractional differential equations. In this context, we discuss a new generalization of the COM-Poisson distribution by using the generalized
-Mittag–Leffler function Equation (
1), independently introduced by Gerhold in [
6] and by Garra and Polito in [
4]. This new fractional distribution includes a wide class of special distributions already studied in the recent literature: the fractional Poisson distribution Equation (
19) considered by Beghin and Orsingher [
28], the classical COM-Poisson distribution and so on.
In this paper, we introduce the following generalization, namely the fractional COM-Poisson distribution
where
is the
-Mittag–Leffler function, considered by Gerhold [
6], Garra and Polito [
4]. The distribution Equation (
23) is a weighted Poisson distribution with weights
. Obviously the probability generating function is here given by a ratio of
-Mittag–Leffler functions, i.e.,
With the next theorem we explain the reason why we consider the distribution Equation (
23), a kind of fractional generalization of the COM-Poisson classical distribution related to Hadamard fractional equations with varying coefficients.
Theorem 1. For , the functionsatisfies the equationwhere is the Caputo-Hadamard fractional derivative of order ν. Proof. By direct calculations we obtain that
where we changed the order of summation with fractional differentiation since the
-Mittag–Leffler function is, for
, an entire function. ☐
A more detailed analysis about the statistical properties and possible applications of this new class of distributions is not the object of this paper and will be considered in a future work.
We finally observe that the fractional COM-Poisson distribution (
23) includes as a special case, for
the fractional Poisson distribution (
19) previously introduced by Beghin and Orsingher in [
28] and recently treated by Chakraborty and Ong [
29], Herrmann [
21] and Porwall and Dixit [
30]. Asymptotic results for the multivariate version of this distribution have been recently analyzed by Beghin and Macci in [
19].
4. Hadamard Fractional Equations Related to the GCOM-Poisson Distribution
In a recent paper, Imoto [
12] studied a different generalization of the Conway-Maxwell-Poisson distribution (in the following GCOM Poisson distribution), depending on two real parameters
r and
. According to [
12], a random variable
has a GCOM Poisson distribution if
where
is the normalizing function which converges for
,
and
(by applying the ratio criterion) or
,
and
. In the context of weighted Poisson distributions, it corresponds to the choice
.
This distribution includes as particular cases the COM-Poisson distribution for and , the geometric for , and the homogeneous Poisson distribution for . The new parameter introduced in the distribution plays the role of controlling the length of tails within the framework of queueing processes. Indeed it was proved that, it displays over-dispersion for and under-dispersion for . Moreover, from the ratios of successive probabilities, when and it emerges to be heavy tailed, while for light-tailed.
We also underline that for
, with
, the function (
29) becomes a sort of generalization of the Wright function which can be written as
Recall that, in the general form, the Wright function is defined as follows (we refer to [
31] for a good survey):
We show now the relation between the normalizing function Equation (
29) and an Hadamard fractional differential equation with time-varying coefficients. Therefore, we can consider in some sense also the GCOM Poisson distribution as a fractional-type modification of the homogeneous distribution, in the sense that its probability generating function satisfies a fractional differential equation.
Proposition 1. The functionsatisfies the integro-differential equationinvolving an Hadamard fractional integral of order r. Proof.
where we used Equation (
13). ☐
Corollary 1. The probability generating function of the GCOM Poisson distribution Equation (28)satisfies the integro-differential equation 5. Further Connections between Modified Mittag–Leffler Functions and Hadamard Fractional Equations
The analysis of fractional differential equations with non-constant coefficients is an interesting and non-trivial topic. In particular, the analysis of equations involving fractional-type Bessel operators (i.e., the fractional counterpart of singular linear differential operators of arbitrary order) has attracted the interest of many researchers (see e.g., [
13,
32] and references therein). In [
4], some results about the connection between Le-Roy functions and equations with space-varying coefficients involving Hadamard derivatives and Laguerre derivatives have been obtained. Here we go further in the direction started in [
4], showing other interesting applications of Hadamard fractional equations in the theory of hyper-Bessel functions. With the next theorem, we find the equation interpolating classical Bessel equations of arbitrary order possessing exact solution in terms of Le Roy functions. We remark, indeed, that Le Roy functions include as special cases hyper-Bessel functions.
Theorem 2. The Le Roy function satisfies the equation Remark 1. For , we have that the function is a solution ofFor , the function satisfiesand so forth for any integer value of α. From this point of view, the solution of Hadamard Equation (37), for non integer values of α, leads to an interpolation between successive hyper-Bessel functions. Let us consider the operator
where
H stands for an Hyper-Bessel type operator involving Caputo-type Hadamard derivatives. We have the following general Theorem.
Theorem 3. A solution of the equationis given by To conclude this section, we restate Theorem 3.3 of [
4], in view of the comments of Turmetov [
33] in the case in which Caputo-type Hadamard derivatives appear in the governing equations.
Theorem 4. The function , with , , is an eigenfunction of the operatorwhere is the integer part of β and denotes the Caputo-type regularized Hadamard derivative of order . The difference w.r.t. the previous version is simply given by the fact that, by using the regularized Caputo-type Hadamard derivative , that is necessary for the correctness of the result.