Solutions to Fractional q-Kinetic Equations Involving Quantum Extensions of Generalized Hyper Mittag-Leffler Functions

Abstract

This manuscript focuses on new generalizations of q-Mittag-Leffler functions, called generalized hyper q-Mittag-Leffler functions, and discusses their regions of convergence and various fractional q operators. Moreover, the solutions to the q-fractional kinetic equations in terms of the investigated generalized hyper q-Mittag-Leffler functions are obtained by applying the q-Sumudu integral transform. Furthermore, we present solutions obtained as numerical graphs using the MATLAB 2018 program.

1. Introduction

Fractional differential equations (FDEs) have received tremendous interest in recent decades. They play a significant role in describing real-world phenomena such as energy [1], financial systems [2], hydrologic modeling [3], electrochemistry [4], mathematical biology [5], and signal processing [6], amongst others [7,8]. Later on, fractional kinetic equations can be seen as an essential part of fractional differential equations, where a wide range of analytical and numerical methods can be used to obtain their solutions. Also, several ways to study the extension and development of fractional kinetic equations involving various special functions have been presented. For example, we refer to contemporary works by Agarwal et al. [9], Saxena et al. [10], Akel et al. [11], Hidan et al. [12], Almalkia and Abdalla [13], Kolokoltsov and Troeva [14], Habenom et al. [15] and Abdalla and Akel [16].

On the other hand, fractional q-differential equations are an extension of the fractional differential equations in quantum calculus (or q-calculus) frameworks. Nowadays, the solution to the quantum analogue (or q-analogue) of the fractional kinetic equation with the Riemann–Liouville fractional q-integral operator is acquired using the technique of a q-analogue of the Laplace transform, which has been discussed by Garg and Chanchlani [17]. Meanwhile, Purohit and Faruk Uçar investigated an application of the q-Sumudu transform for fractional q-kinetic equations in [18]. Recently, Bairwa et al. [19] obtained the solution to the fractional q-kinetic equation pertaining to the generalized q-Mittag-Leffler function via the q-Sumudu transform. In a similar vein, Abujarad al. [20] introduced and investigated the solutions to q-fractional kinetic equations, including the generalized hyper Bessel function, using the q-Shehu transform.

Motivated by the previous works, we define new quantum extensions of the generalized hyper Mittag-Leffler functions in this article and establish their properties and applications in generalized fractional q-kinetic equations. This work’s organization is as follows. Section 2 presents the basic concepts of q-calculus and reviews the versions of Mittag-Leffler functions that are required in the following sections. In Section 3, we introduce the two q-analogues of the generalized hyper Mittag-Leffler functions and study the convergence properties of these functions for  $\left|q\right|<1$. Some of the fractional q-operators of the generalized hyper q-Mittag-Leffler functions using the Caputo fractional q-derivatives, the Hilfer fractional q-derivatives, the Kober q-integrals, and the Riemann-Liouville q-integrals are proposed in Section 4. In Section 5, we find the solutions to the q-fractional kinetic equations involving the generalized hyper q-Mittag-Leffler functions by using the q-Sumudu transform. In Section 6, we point out the graphical representations of the primary obtained solutions in Section 4 using the MATLAB program. Eventually, we exhibit some concluding remarks in Section 7.

2. Mathematical Preliminaries

2.1. Some Basic Concepts of q-Calculus

The quantum calculus (q-calculus) is the q extension of ordinary calculus. The theory of q-calculus operators has recently been applied in modern sciences (e.g., see [21,22,23]).

Now, we recall some of the basic concepts and related details of q-calculus which will be used in the current study.

Definition 1

(See [24]). Given a value of  $0<\left|q\right|<1$, the q-number  ${\left[\kappa \right]}_q$ is defined as

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\left[\kappa \right]}_q:=\left\{\begin{array}{cc}\frac{q^{\kappa }-1}{q-1},\hfill & q\ne 1,\hfill \\ \\ \kappa ,\hfill & q=1,\kappa \in \mathbb{R}.\hfill \end{array}\right.\end{array}\end{array}$$

(1)

For  $\kappa \in \mathbb{N}$, the q-factorial  ${\left[\kappa \right]}_q!$ is given by

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\left[\kappa \right]}_q!=\left\{\begin{array}{cc}{\left[\kappa \right]}_q{[\kappa -1]}_q\dots {\left[1\right]}_q,\hfill & \kappa \in \mathbb{N},\hfill \\ \\ 1,\hfill & \kappa =0.\hfill \end{array}\right.\end{array}\end{array}$$

(2)

Definition 2

(See [24]). For  $0<\left|q\right|<1$, the q-analogue of the common Pochhammer symbols (q-shifted factorials) are given as follows:

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\left(\theta ;q\right)}_{\kappa }=\left\{\begin{array}{cc}{\prod }_{ı=0}^{\kappa -1}\left(1-\theta q^{ı}\right),\hfill & \kappa \in \mathbb{N},\hfill \\ \\ 1,\hfill & \kappa =0,\hfill \end{array}\right.\end{array}\end{array}$$

(3)

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\left(\theta ;q\right)}_{\infty }=\prod _{ı=0}^{\infty }\left(1-\theta q^{ı}\right),\theta \in \mathbb{C}\end{array}\end{array}$$

(4)

and

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\left(\theta ;q\right)}_{\eta }=\frac{{\left(\theta ;q\right)}_{\infty }}{{\left(\theta q^{\eta };q\right)}_{\infty }},\theta ,\eta \in \mathbb{C}.\end{array}\end{array}$$

(5)

Definition 3

(See [24]). For  $0<\left|q\right|<1,$ the q-analogue of  ${(u-v)}_q^{\kappa }$ is given as follows:

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\left(u-v\right)}_q^{\kappa }=u^{\kappa }{\left(v/u;q\right)}_{\kappa }=u^{\kappa }\frac{{\left(v/u;q\right)}_{\infty }}{{\left(v{q}^{\kappa }/u;q\right)}_{\infty }},u\ne 0.\end{array}\end{array}$$

(6)

Definition 4

(See [24,25]). For  $0<\left|q\right|<1,$ the q-gamma function is defined by

$$\begin{array}{cc}\hfill {\Gamma }_q\left(\theta \right)=& \frac{{(q;q)}_{\infty }}{\left(q^{\theta };q\right)\infty }{(1-q)}^{1-\theta }=\frac{{(q;q)}_{\theta -1}}{{(1-q)}^{\theta -1}}Re\left(\theta \right)>0,\hfill \end{array}$$

(7)

where  ${lim}_{q\to 1}{\Gamma }_q\left(\theta \right)=\Gamma \left(\theta \right)$ and  ${\Gamma }^{}\left(\theta \right)$ is the classical gamma function

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\Gamma }^{}\left(\theta \right)={\int }_0^{\infty }{\tau }^{\theta -1}{e}^{-\tau }d\tau ,\theta \in \mathbb{C}\setminus {\mathbb{Z}}^{-}.\end{array}\end{array}$$

(8)

Note that

$$\begin{array}{c}\hfill {\Gamma }_q(\theta +1)={\left[\theta \right]}_q{\Gamma }_q\left(\theta \right)={\left[\theta \right]}_q!Re\left(\theta \right)>0\end{array}$$

Moreover, the q-analogue of Stirling’s asymptotic formula for the q-gamma function (Equation (7)) is defined in [25] by

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\Gamma }_q\left(\theta \right)\sim {(1+q)}^{\frac{1}{2}}{\Gamma }_{q^2}\left(\frac{1}{2}\right){(1-q)}^{\frac{1}{2}-\theta }{e}^{{{\rm Y}}_q\left(\theta \right)},\end{array}\end{array}$$

(9)

where  ${{\rm Y}}_q\left(\theta \right)=\frac{\phi q^{\theta }}{1-q-q^{\theta }},0<\phi <1.$

Definition 5

(See [24,25]). For  $0<q<1$, the q-beta function is defined by

$$\begin{array}{c}\hfill B_q(\theta ,\phi )=\frac{{\Gamma }_q\left(\theta \right){\Gamma }_q\left(\phi \right)}{{\Gamma }_q(\theta +\phi )},Re\left(\theta \right)>0,Re\left(\phi \right)>0.\end{array}$$

(10)

Definition 6

(See [24]). The q-derivative  ${\mathbf{D}}_q\Psi$ of a function  $\Psi \left(\eta \right)$ is defined as

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\mathbf{D}}_q\Psi \left(\eta \right)=\frac{\Psi \left(\eta \right)-\Psi \left(q\eta \right)}{\eta (1-q)},\eta \ne 0,q\ne 1,\end{array}\end{array}$$

(11)

and  $\left({\mathbf{D}}_q\Psi \right)\left(0\right)={\Psi }^{\prime }\left(0\right)$, provided  ${\Psi }^{\prime }\left(0\right)$ exists.

We observe that

$$\begin{array}{c}\hfill {\mathbf{D}}_q^s\left({\eta }^r\right)=\frac{{\Gamma }_q(r+1)}{{\Gamma }_q(r-s+1)}{\eta }^{r-s},r,s\in \mathbb{N},\end{array}$$

(12)

And if  $\Psi \left(\eta \right)$ is differentiable, then one has

$$\begin{array}{c}\hfill \underset{q\to 1}{lim}{\mathbf{D}}_q\Psi \left(\eta \right)=\underset{q\to 1}{lim}\frac{\Psi \left(q\eta \right)-\Psi \left(\eta \right)}{\eta (q-1)}=\frac{d\Psi \left(\eta \right)}{d\eta }.\end{array}$$

Definition 7

(See [24]). For  $0<\left|q\right|<1,$ the q-Jackson integral is defined as

$$\begin{array}{c}\hfill {\int }_0^{\eta }f\left(\tau \right)d_q\tau =(1-q)\eta \sum _{\iota =0}^{\infty }{q}^{\iota }f\left(q^{\iota }\eta \right),\end{array}$$

(13)

provided the sum converges absolutely.

Definition 8

(See [21,22]). The Riemann–Liouville fractional q-integral operator is defined by

$$\begin{array}{c}\hfill I_q^{\nu }f\left(\eta \right):=\frac{1}{{\Gamma }_q\left(\nu \right)}{\int }_0^{\eta }{\left(\eta -qu\right)}_q^{\nu -1}f\left(u\right)d_{q}u,Re\left(\nu \right)>0.\end{array}$$

(14)

It is clear that

$$\begin{array}{c}\hfill I_q^{\nu }\left[{\eta }^{\mu }\right]=\frac{{\Gamma }_q(1+\mu )}{{\Gamma }_q(1+\mu +\nu )}{\eta }^{\mu +\nu },\mu \in (-1,\infty ).\end{array}$$

(15)

Definition 9

(See [21,22]). The Kober fractional q-integral operator is defined by

$$\begin{array}{c}\hfill I_q^{\mu ,\nu }f\left(\eta \right):=\frac{{\eta }^{-\mu -\nu }}{{\Gamma }_q\left(\nu \right)}{\int }_0^{\eta }{\left(\eta -qu\right)}_q^{\nu -1}{u}^{\mu }f\left(u\right)d_{q}u,Re\left(\mu \right)>0.\end{array}$$

(16)

Definition 10

(See [21,22]). The Caputo fractional q-derivative is defined by

$$\begin{array}{c}\hfill {}^c{D}_q^{\nu }\left[f\right]\left(\eta \right)=I_q^{1-\nu }{D}_q\left[f\right]\left(\eta \right)=\frac{1}{{\Gamma }_q(1-\nu )}{\int }_0^{\eta }{(\eta -q\tau )}_q^{-\nu }{D}_{q}f\left(\tau \right)d_q\tau ,Re\left(\nu \right)>0.\end{array}$$

(17)

Definition 11

(See [21,22]). The Hilfer fractional q-derivative of an order ν and of type μ is defined by

$$\begin{array}{c}\hfill D_q^{\mu ,\nu }\left[f\right]\left(\eta \right)=I_q^{\mu (1-\nu )}{D}_q{I}_q^{(1-mu)(1-\nu )}\left[f\right]\left(\eta \right),0<\nu <1,0⩽\mu ⩽1.\end{array}$$

(18)

Definition 12.

For  $0<\left|q\right|<1,$ the two q-analogues of the exponential function are given in [22,24] as follows:

$$\begin{array}{c}\hfill {\mathfrak{E}}_q\left(\eta \right)=\sum _{\iota =0}^{\infty }\frac{{\eta }^{\iota }}{{(q;q)}_{\iota }}=\frac{1}{{(\eta ;q)}_{\infty }},\left|\eta \right|<1,\end{array}$$

(19)

and

$$\begin{array}{c}\hfill {\mathbb{E}}_q\left(\eta \right)=\sum _{\iota =0}^{\infty }{(-1)}^{\iota }{q}^{\iota (\iota -1)/2}\frac{{\eta }^{\iota }}{{(q;q)}_{\iota }}={(\eta ;q)}_{\infty },\left|\eta \right|<\infty .\end{array}$$

(20)

2.2. Some Versions of the Mittag-Leffler Function

Following Swedish mathematician Gosta Mittag-Leffler [26], the Mittag-Leffler (M-L) function is defined by

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\mathbf{E}}_{\alpha }\left(\xi \right)=\sum _{k=0}^{\infty }\frac{{\xi }^k}{\Gamma (k\alpha +1)},\xi ,\alpha \in \mathbb{C},Re\left(\alpha \right)>0,\end{array}\end{array}$$

(21)

where  ${\Gamma }^{}(.)$ is the gamma function (Equation (8)).

This function is the direct generalization of the exponential function

$$\begin{array}{c}\hfill {\mathbf{E}}_1\left(\xi \right)=e^{\xi }=\sum _{k=0}^{\infty }\frac{{\xi }^k}{k!}\xi \in \mathbb{C}.\end{array}$$

(22)

Various generalizations of the M-L function are provided below.

Generalizations of  ${\mathbf{E}}_{\alpha }\left(\xi \right)$ first appeared in a 1905 work by Wiman [27] as

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\mathbf{E}}_{\alpha ,\beta }\left(\xi \right)=\sum _{k=0}^{\infty }\frac{{\xi }^k}{\Gamma (k\alpha +\beta )}\xi ,\alpha ,\beta \in \mathbb{C},Re\left(\alpha \right)>0,Re\left(\beta \right)>0.\end{array}\end{array}$$

(23)

The function  ${\mathbf{E}}_{\alpha ,\beta }\left(\xi \right)$ is the natural extension of certain functions like hyperbolic trigonometric functions, among others (refer to [26,27]).

Prabhakar [28] introduced the following extension of Equation (23):

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\mathbf{E}}_{\alpha ,\beta }^{\delta }\left(\xi \right)=\sum _{k=0}^{\infty }\frac{{\left(\delta \right)}_k}{\Gamma (k\alpha +\beta )}\frac{{\xi }^k}{k!},(\xi ,\alpha ,\beta ,\delta \in \mathbb{C},Re\left(\alpha \right)>0,Re\left(\beta \right)>0),\end{array}\end{array}$$

(24)

where  ${\left(\delta \right)}_m$ denotes the usual Pochhammer symbol defined by

$$\begin{array}{c}\hfill {\left(\delta \right)}_k=\frac{\Gamma (\delta +m)}{\Gamma \left(\delta \right)}=\left\{\begin{array}{cc}\delta (\delta +1)...(\delta +m-1),\hfill & m\in \mathbb{N},\delta \in \mathbb{C}\hfill \\ \\ 1,\hfill & m=0;\delta \in \mathbb{C}\setminus \left\{0\right\},\hfill \end{array}\right.\end{array}$$

(25)

and  $\Gamma (.)$ is the gamma function defined by Equation (8).

In [29], Shukla and Prajapati presented the following function:

$$\begin{array}{c}\hfill \begin{array}{cc}& {\mathbf{E}}_{\alpha ,\beta }^{\delta ,\nu }\left(\xi \right)=\sum _{k=0}^{\infty }\frac{{\left(\delta \right)}_{\nu k}}{\Gamma (k\alpha +\beta )}\frac{{\xi }^k}{k!},\hfill \\ & (\xi ,\alpha ,\beta ,\delta \in \mathbb{C},Re(\alpha )>0,Re(\beta )>0,\nu \in (0,1)\cup \mathbb{N}),\hfill \end{array}\end{array}$$

(26)

where  ${\left(\delta \right)}_{\nu k}$ is the generalized Pochhammer symbol given as follows:

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\left(\delta \right)}_{\nu k}=\frac{\Gamma (\delta +k\nu )}{\Gamma \left(\delta \right)}=k^{\nu k}\prod _{j=1}^k\left(\frac{\delta +j-1}{k}\right),k\in \mathbb{N}.\end{array}\end{array}$$

(27)

Another generalization of Equation (23) according to Salim [30,31] is given by

$$\begin{array}{c}\hfill \begin{array}{cc}& {\mathbf{E}}_{\alpha ,\beta }^{\delta ,\nu }\left(\xi \right)=\sum _{k=0}^{\infty }\frac{{\left(\delta \right)}_k}{\Gamma (k\alpha +\beta ){\left(\nu \right)}_k}\frac{{\xi }^k}{k!},\hfill \\ & (\xi ,\alpha ,\beta ,\delta ,\nu \in \mathbb{C},Re(\alpha )>0,Re(\beta )>0,Re(\delta )>0,Re(\nu )>0).\hfill \end{array}\end{array}$$

(28)

On the other hand, we will mention some q-analogues of the Mittag-Leffler-type functions functions.

Jain [32] introduced a q-analogue of the generalized Mittag-Leffler function:

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\mathbb{E}}_{\alpha }(w;q)=\sum _{h=0}^{\infty }\frac{w^h}{{\Gamma }_q(\alpha h+1)},(Re\left(\alpha \right)>0).\end{array}\end{array}$$

(29)

Mansour [33] introduced a new form of a q-analogue of the Mittag-Leffler function, given by

$$\begin{array}{c}\hfill \begin{array}{c}\hfill {\mathbb{E}}_{\alpha ,\delta }(w;q)=\sum _{h=0}^{\infty }\frac{w^h}{{\Gamma }_q(h\alpha +\delta )},\left(\alpha >0,\delta \in \mathbb{C},\left(\right|w|<{(1-q)}^{-\alpha })\right).\end{array}\end{array}$$

(30)

Meanwhile, Sharma and Jain [34] presented and established the q-analogue of the Mittag-Leffler function below:

$$\begin{array}{c}\hfill {\mathfrak{E}}_{\alpha ,\delta }^{\mathsf{\gamma }}(w;q)=\sum _{h=0}^{\infty }\frac{{(q^{\mathsf{\gamma }};q)}_h}{{\Gamma }_q(h\alpha +\delta )}·\frac{w^h}{{(q;q)}_h},\left({\left|w\right|<(1-q)}^{-\alpha }\right),\end{array}$$

(31)

where  $\alpha >0,\mathsf{\gamma }>0,$ $\delta ,w\in \mathbb{C}$, and  $\left|q\right|<1.$

Furthermore, Purohit and Kalla [35] defined and investigated another generalization of Equation (30) in the form

$$\begin{array}{c}\hfill {\mathbb{E}}_{\alpha ,\delta }^{\mathsf{\gamma }}(w;q)=\sum _{h=0}^{\infty }\frac{q^{h(h-1)/2}{(q^{\mathsf{\gamma }};q)}_h}{{\Gamma }_q(h\alpha +\delta )}·\frac{w^h}{{(q;q)}_h},\left({\left|w\right|<(1-q)}^{-\alpha }\right),\end{array}$$

(32)

where  $\alpha >0,\mathsf{\gamma }>0,$ $\delta ,w\in \mathbb{C}$, and  $\left|q\right|<1.$

Further recent generalizations and applications of the q-Mittag-Leffler function can be found in [36,37,38,39].

3. The Generalized Hyper  $\mathit{q}$-Mittag-Leffler Functions

Motivated by the studies quoted in Section 2.2, we introduce here the two q-analogues of the generalized hyper Mittag-Leffler function in the following definition:

Definition 13.

Let  $\alpha , \beta , {\lambda }_s, {\delta }_r\in \mathbb{C}$ such that  $Re\left(\alpha \right)>0,$   $Re\left({\lambda }_s\right)>0$, and  $Re\left({\delta }_r\right)>0$ for the finite sets  $s\in \left\{1,2,\cdots ,d\right\}$ and  $r\in \left\{1,2,\cdots ,h\right\}.$ Let  $\mathit{\lambda }=({\lambda }_1,{\lambda }_2,\cdots ,{\lambda }_d)$,  $\mathit{\delta }=({\delta }_1,{\delta }_2,\cdots ,{\delta }_h)$, and  $0<\left|q\right|<1$. Then, the function

$$\begin{array}{c}\hfill {\mathfrak{E}}_{\alpha ,\beta ,\mathit{\delta }}^{\mathit{\lambda }}(w;q)=\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\ell }}{{(q;q)}_{\ell }},\end{array}$$

(33)

is called the small generalized hyper q-Mittag-Leffler function. Similarly, the large generalized hyper q-Mittag-Leffler function is defined as follows:

$$\begin{array}{c}\hfill {\mathbb{E}}_{\alpha ,\beta ,\mathit{\delta }}^{\mathit{\lambda }}(w;q)=\sum _{\ell =0}^{\infty }\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\ell }}{{(q;q)}_{\ell }}.\end{array}$$

(34)

Alternatively, in light of the definition of the q-gamma function in Equation (7), we can rewrite Equations (33) and (34) in the respective forms

$$\begin{array}{c}\hfill {\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)=\sum _{\ell =0}^{\infty }\frac{{\left(q^{\ell \alpha +\beta };q\right)}_{\infty }{\Pi }_{r=1}^h{\left(q^{{\delta }_r+\ell };q\right)}_{\infty }}{{\Pi }_{s=1}^k{\left(q^{{\lambda }_s+\ell };q\right)}_{\infty }{\Pi }_{s=1}^k{\Pi }_{r=1}^h{\left(q^{{\delta }_r-{\lambda }_s+\ell \alpha +\beta };q\right)}_{\infty }}·\frac{w^{\ell }}{{(q;q)}_{\ell }}\end{array}$$

(35)

and

$$\begin{array}{c}\hfill {\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)=\sum _{\ell =0}^{\infty }\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}{\left(q^{\ell \alpha +\beta };q\right)}_{\infty }{\Pi }_{r=1}^h{\left(q^{{\delta }_r+\ell };q\right)}_{\infty }}{{\Pi }_{s=1}^k{\left(q^{{\lambda }_s+\ell };q\right)}_{\infty }{\Pi }_{s=1}^k{\Pi }_{r=1}^h{\left(q^{{\delta }_r-{\lambda }_s+\ell \alpha +\beta };q\right)}_{\infty }}·\frac{w^{\ell }}{{(q;q)}_{\ell }}.\end{array}$$

(36)

Remark 1.

We note the following special cases for Equations (33) and (34):

• By putting  $\lambda ={\lambda }_1$ and  $\delta ={\delta }_1$ in Equations (33) and (34), we obtain the following generalized q-Mittag-Leffler functions:

  $$\begin{array}{c}\hfill {\mathit{e}}_{\alpha ,\beta ,{\delta }_1}^{{\lambda }_1}(w;q)=\sum _{\ell =0}^{\infty }\frac{{\left[{\Gamma }_q({\lambda }_1+\ell )\right]}^d}{{\Gamma }_q(\ell \alpha +\beta ){\left[{\Gamma }_q({\delta }_1+\ell )\right]}^h}·\frac{w^{\ell }}{{(q;q)}_{\ell }},\end{array}$$

  (37)

  and

  $$\begin{array}{c}\hfill {\mathit{E}}_{\alpha ,\beta ,{\delta }_1}^{{\lambda }_1}(w;q)=\sum _{\ell =0}^{\infty }\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}{\left[{\Gamma }_q({\lambda }_1+\ell )\right]}^d}{{\Gamma }_q(\ell \alpha +\beta ){\left[{\Gamma }_q({\delta }_1+\ell )\right]}^h}·\frac{w^{\ell }}{{(q;q)}_{\ell }}.\end{array}$$

  (38)
• By setting  $q\to 1$ in Equations (33) and (34), we reduce them to the following generalized hyper Mittag-Leffler function:

  $$\begin{array}{c}\hfill {\mathcal{E}}_{\alpha ,\beta ,\delta }^{\lambda }\left(w\right)=\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d\Gamma ({\lambda }_s+\ell )}{\Gamma (\ell \alpha +\beta ){\Pi }_{r=1}^h\Gamma ({\delta }_r+\ell )}·\frac{w^{\ell }}{\ell !},\end{array}$$

  (39)

  which is a generalization of the Mittag-Leffler functions defined by Equations (21)–(24) and (26).
• By taking  $d=h=1$ in Equations (37) and (38), we give the following generalized q-Mittag-Leffler functions with four parameters:

  $$\begin{array}{c}\hfill {\mathrm{e}}_{\alpha ,\beta ,{\delta }_1}^{{\lambda }_1}(w;q)=\sum _{\ell =0}^{\infty }\frac{{\Gamma }_q({\lambda }_1+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Gamma }_q({\delta }_1+\ell )}·\frac{w^{\ell }}{{(q;q)}_{\ell }}\end{array}$$

  (40)

  and

  $$\begin{array}{c}\hfill {\mathrm{E}}_{\alpha ,\beta ,{\delta }_1}^{{\lambda }_1}(w;q)=\sum _{\ell =0}^{\infty }\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}{\Gamma }_q({\lambda }_1+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Gamma }_q({\delta }_1+\ell )}·\frac{w^{\ell }}{{(q;q)}_{\ell }}.\end{array}$$

  (41)

  which generalize the Mittag-Leffler functions defined by Equations (20), (29), and (32).
• When  $d=h=1$ and  ${\delta }_1={\lambda }_1$ in Equations (40) and (41), we obtain the following generalized q-Mittag-Leffler functions:

  $$\begin{array}{c}\hfill {\mathfrak{E}}_{\alpha ,\beta }^{}(w;q)=\sum _{\ell =0}^{\infty }\frac{1}{{\Gamma }_q(\ell \alpha +\beta )}·\frac{w^{\ell }}{{(q;q)}_{\ell }}\end{array}$$

  (42)

  and

  $$\begin{array}{c}\hfill {\mathbb{E}}_{\alpha ,\beta }^{}(w;q)=\sum _{\ell =0}^{\infty }\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}}{{\Gamma }_q(\ell \alpha +\beta )}·\frac{w^{\ell }}{{(q;q)}_{\ell }}.\end{array}$$

  (43)

We now discuss the convergence properties of the generalized hyper q-Mittag-Leffler functions defined in Equations (33) and (34). First, the following theorem ensures the convergence and order of Equation (34):

Theorem 1.

Under the conditions and hypothesis in Definition 13,  ${\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)$ is an entire function for  $0<\left|q\right|<1$, and its order is equal to zero.

Proof. 

For this proof, we can rewrite Equation (34) in the form

$$\begin{array}{c}\hfill {\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)=\sum _{\ell =0}^{\infty }\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\ell }}{{(q;q)}_{\ell }}=\sum _{\ell =0}^{\infty }{\mathbb{U}}_{\ell }^q{w}^{\ell },\end{array}$$

(44)

where

$$\begin{array}{c}\hfill {\mathbb{U}}_{\ell }^q=\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{1}{{(q;q)}_{\ell }}.\end{array}$$

(45)

Using the relation in Equation (9), and after simplifications, we arrive at the radius of regularity R of the above series as follows:

$$\begin{array}{cc}\hfill \frac{1}{\mathrm{R}}& =\underset{\ell \to \infty }{lim}\sqrt[\ell ]{\left|{\mathbb{U}}_{\ell }^q\right|}\hfill \\ \hfill & \simeq \underset{\ell \to \infty }{lim}|\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}}{{(1+q)}^{\frac{1}{2}}{\Gamma }_{q^2}\left(\frac{1}{2}\right){(1-q)}^{\frac{1}{2}-(\ell \alpha +\beta )}{e}^{{{\rm Y}}_q(\ell \alpha +\beta )}}{|}^{\frac{1}{\ell }}\hfill \\ \hfill & \times |\frac{{\Pi }_{s=1}^d{(1+q)}^{\frac{1}{2}}{\Gamma }_{q^2}\left(\frac{1}{2}\right){(1-q)}^{\frac{1}{2}-({\lambda }_s+\ell )}{e}^{{{\rm Y}}_q({\lambda }_s+\ell )}}{{\Pi }_{r=1}^h{(1+q)}^{\frac{1}{2}}{\Gamma }_{q^2}\left(\frac{1}{2}\right){(1-q)}^{\frac{1}{2}-({\delta }_r+\ell )}{e}^{{{\rm Y}}_q({\delta }_r+\ell )}}{|}^{\frac{1}{\ell }}\hfill \\ \hfill & \times |\frac{1}{{(1+q)}^{\frac{1}{2}}{\Gamma }_{q^2}\left(\frac{1}{2}\right){(1-q)}^{\frac{-1}{2}}{e}^{{{\rm Y}}_q(\ell +1)}}{|}^{\frac{1}{\ell }}\hfill \\ \hfill & \simeq \underset{\ell \to \infty }{lim}\left|q^{(\ell -1)/2}\right|=0,0<\left|q\right|<1,\hfill \end{array}$$

where  ${{\rm Y}}_q$ is defined in Equation (9). Therefore, we have  $R=\infty$, and hence the function  ${\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)$ is an entire function. Now, we calculate the order of the entire function  ${\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)$ as follows (see [21,22]):

$$\begin{array}{c}\hfill \rho \left({\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)\right)=\underset{\ell \to \infty }{lim\; sup}\frac{\ell ln\ell }{ln\left(\frac{1}{\left|{\mathbb{U}}_{\ell }^q\right|}\right)},\end{array}$$

which given from Equation (45) as

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill ln\left(\frac{1}{\left|{\mathbb{U}}_{\ell }^q\right|}\right)& =ln\left[\frac{{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell ){(1-q)}^{\ell }{\Gamma }_q(\ell +1)}{q^{\ell (\ell -1)/2}{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}\right]\hfill \\ & =ln\left|{\Gamma }_q(\ell \alpha +\beta )\right|+ln\left|{\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )\right|+ln\left|{\Gamma }_q(\ell +1)\right|\hfill \\ & +\ell ln|1-q|-\ell (\ell -1)/2ln\left|q\right|-ln\left|{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )\right|.\hfill \end{array}\end{array}$$

(46)

Then, in light of Equation (7), we see that

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill ln\left|{\Gamma }_q(\ell \alpha +\beta )\right|=& ln\left|\frac{{(q;q)}_{\infty }}{{\left(q^{(\ell \alpha +\beta )};q\right)}_{\infty }}{(1-q)}^{1-(\ell \alpha +\beta )}\right|\hfill \\ \hfill =& ln\left|{(q;q)}_{\infty }\right|-ln\left|{\left(q^{(\ell \alpha +\beta )};q\right)}_{\infty }\right|\hfill \\ \hfill +& (1-(\ell Re\left(\alpha \right)+Re\left(\beta \right)\left)\right)ln\left|\right(1-q\left)\right|.\hfill \end{array}\end{array}$$

(47)

According to [21], we observe that

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill ln\left|{\left(q^{(\ell \alpha +\beta )};q\right)}_{\infty }\right|& =ln\left[\prod _{m=0}^{\infty }|1-q^{(\ell \alpha +\beta +m)}|\right]\hfill \\ \hfill =& ln\left[\underset{\omega \to \infty }{lim}\prod _{m=0}^{\omega }|1-q^{(\ell \alpha +\beta +m)}|\right]\hfill \\ \hfill =& \underset{\omega \to \infty }{lim}\sum _{m=0}^{\omega }ln|1-q^{(\ell \alpha +\beta +m)}|\hfill \\ \hfill =& \sum _{m=0}^{\infty }ln|1-q^{(\ell \alpha +\beta +m)}|\hfill \\ & ⩽\sum _{m=0}^{\infty }{q}^{(\ell Re(\alpha +\beta )+m)}=\frac{q^{(\ell Re(\alpha +\beta \left)\right)}}{1-q}.\hfill \end{array}\end{array}$$

(48)

Consequently, from Equation (47), we have

$$\begin{array}{c}\hfill \underset{\ell \to \infty }{lim}\frac{ln\left|{\Gamma }_q(\ell \alpha +\beta )\right|}{\ell ln\ell }=0.\end{array}$$

(49)

Thus, using the above technique in Equation (46), one would obtain

$$\begin{array}{c}\hfill \underset{\ell \to \infty }{lim}\frac{ln\left(\frac{1}{\left|{\mathbb{U}}_{\ell }^q\right|}\right)}{\ell ln\ell }=\infty .\end{array}$$

(50)

Hence,  $\rho \left({\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)\right)=0.$ This completes the proof of Theorem 1. □

Second, the following result includes the absolute convergence of Equation (33):

Theorem 2.

Under the conditions and hypothesis in Definition 13, the function  ${\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)$ (Equation (33)) converges absolutely for  ${\left|w\right|<(1-q)}^{-\alpha }$ with  $0<\left|q\right|<1.$

Proof. 

To prove that we can rewrite Equation (33) as

$$\begin{array}{c}\hfill {\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)=\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\ell }}{{(q;q)}_{\ell }}=\sum _{\ell =0}^{\infty }{\mathbb{U}}_{\ell }^q{w}^{\ell },\end{array}$$

(51)

where

$$\begin{array}{c}\hfill {\mathbb{U}}_{\ell }^q=\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{1}{{(1-q)}^{\ell }{\Gamma }_q(\ell +1)}.\end{array}$$

(52)

using the relation in Equation (9), and after calculations, it follows that

$$\begin{array}{cc}\hfill \frac{1}{\mathrm{R}}& =\underset{\ell \to \infty }{lim}\sqrt[\ell ]{\left|{\mathbb{U}}_{\ell }^q\right|}\hfill \\ \hfill & \simeq \underset{\ell \to \infty }{lim}|\frac{1}{{(1+q)}^{\frac{1}{2}}{\Gamma }_{q^2}\left(\frac{1}{2}\right){(1-q)}^{\frac{1}{2}-(\ell \alpha +\beta )}{e}^{{{\rm Y}}_q(\ell \alpha +\beta )}}{|}^{\frac{1}{\ell }}\hfill \\ \hfill & \times |\frac{{\Pi }_{s=1}^d{(1+q)}^{\frac{1}{2}}{\Gamma }_{q^2}\left(\frac{1}{2}\right){(1-q)}^{\frac{1}{2}-({\lambda }_s+\ell )}{e}^{{{\rm Y}}_q({\lambda }_s+\ell )}}{{\Pi }_{r=1}^h{(1+q)}^{\frac{1}{2}}{\Gamma }_{q^2}\left(\frac{1}{2}\right){(1-q)}^{\frac{1}{2}-({\delta }_r+\ell )}{e}^{{{\rm Y}}_q({\delta }_r+\ell )}}{|}^{\frac{1}{\ell }}\hfill \\ \hfill & \times |\frac{1}{{(1+q)}^{\frac{1}{2}}{\Gamma }_{q^2}\left(\frac{1}{2}\right){(1-q)}^{\frac{-1}{2}}{e}^{{{\rm Y}}_q(\ell +1)}}{|}^{\frac{1}{\ell }}\hfill \\ \hfill & \simeq \underset{\ell \to \infty }{lim}|{\Pi }_{r=1}^h{\Pi }_{s=1}^d{(1-q)}^{-{\lambda }_s+{\delta }_r+\ell \alpha -1}{|}^{\frac{1}{\ell }}\simeq \underset{\ell \to \infty }{lim}|{(1-q)}^{\alpha }|.\hfill \end{array}$$

Thus, the function  ${\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)$ is convergent for  ${\left|w\right|<(1-q)}^{-\alpha }$ if  $0<\left|q\right|<1.$ We thus arrive at the required result. □

Remark 2.

According to Remark 1, many particular cases can be inserted by Theorem 1 and Theorem 2.

4. Fractional  $\mathit{q}$-Calculus Approach

In this section, we present some of the fractional q-operators of the generalized hyper q-Mittag-Leffler functions in Equations (33) and (34) in the following theorems:

Theorem 3

(Caputo fractional q-derivatives). The Caputo fractional q-derivatives of the generalized hyper q-Mittag-Leffler functions in Equations (33) and (34) are the following respective equations:

$$\begin{array}{c}\hfill {}^c{D}_q^{\nu }\left\{{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w^{\nu };q)\right\}=w^{-\nu }\sum _{\ell =0}^{\infty }\frac{{\Gamma }_q(1+\nu \ell ){\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Gamma }_q(1-\nu +\nu \ell ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\nu \ell }}{{(q;q)}_{\ell }},\end{array}$$

(53)

and

$$\begin{array}{c}\hfill {}^c{D}_q^{\nu }\left\{{\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w^{\nu };q)\right\}=w^{-\nu }\sum _{\ell =0}^{\infty }\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}{\Gamma }_q(1+\nu \ell ){\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Gamma }_q(1-\nu +\nu \ell ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\nu \ell }}{{(q;q)}_{\ell }}.\end{array}$$

(54)

Proof. 

By applying Equation (17) to Equation (33), we have

$$\begin{array}{cc}\hfill {}^c{D}_q^{\nu }\left\{{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w^{\nu };q)\right\}& =\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}{·}^c{D}_q^{\nu }\left[w^{\nu \ell }\right]\hfill \\ \hfill & =\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·I_q^{\nu -1}{D}_q\left[w^{\nu \ell }\right].\hfill \end{array}$$

Using Equations (12) and (15), one finds

$$\begin{array}{cc}\hfill {}^c{D}_q^{\nu }\left\{{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w^{\nu };q)\right\}& =\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·{[\nu \ell ]}_q{I}_q^{\nu -1}\left[w^{\nu \ell -1}\right]\hfill \\ \hfill & =\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{{[\nu \ell ]}_q{\Gamma }_q(\nu \ell )}{{\Gamma }_q(1-\nu +\nu \ell )}{w}^{\nu \ell -\nu }\hfill \\ \hfill & =\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{{\Gamma }_q(1+\nu \ell )}{{\Gamma }_q(1-\nu +\nu \ell )}{w}^{\nu \ell -\nu }.\hfill \end{array}$$

This is the required result (Equation (53)). Similarly, Equation (54) can be proven. □

Theorem 4

(Hilfer fractional q-derivatives). The Hilfer fractional q-derivatives of the generalized hyper q-Mittag-Leffler functions in Equations (33) and (34) are the following respective equations:

$$\begin{array}{c}\hfill D_q^{\nu ,\mu }\left\{{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w^{\nu };q)\right\}=w^{-\nu }\sum _{\ell =0}^{\infty }\frac{{\Gamma }_q(1+\nu \ell ){\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Gamma }_q(1-\nu +\nu \ell ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\nu \ell }}{{(q;q)}_{\ell }},\end{array}$$

(55)

and

$$\begin{array}{c}\hfill D_q^{\nu ,\mu }\left\{{\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w^{\nu };q)\right\}=w^{-\nu }\sum _{\ell =0}^{\infty }\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}{\Gamma }_q(1+\nu \ell ){\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Gamma }_q(1-\nu +\nu \ell ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\nu \ell }}{{(q;q)}_{\ell }}.\end{array}$$

(56)

where  $D_q^{\nu ,\mu }$ is the Hilfer fractional q-differential operator (Equation (18)).

Proof. 

By applying Equation (18) to the function  $f\left(w\right)=w^{\nu \ell }$, we obtain

$$\begin{array}{cc}\hfill D_q^{\nu ,\mu }\left[w^{\nu \ell }\right]& =I_q^{\mu (1-\nu )}{D}_q{I}_q^{(1-\mu )(1-\nu )}\left[w\right]\hfill \\ \hfill & =I_q^{\mu (1-\nu )}{D}_q\frac{{\Gamma }_q(1+\nu \ell )}{{\Gamma }_q\left(1+\nu \ell +(1-\mu )(1-\nu )\right)}{w}^{(1-\mu )(1-\nu )+\nu \ell }\hfill \\ \hfill & =I_q^{\mu (1-\nu )}\frac{{\Gamma }_q(1+\nu \ell )}{{\Gamma }_q\left(1+\nu \ell +(1-\mu )(1-\nu )\right)}\hfill \\ \hfill & \times {\left[(1-\mu )(1-\nu )+\nu \ell \right]}_q{w}^{(1-\mu )(1-\nu )+\nu \ell -1}\hfill \\ \hfill & =\frac{{\Gamma }_q(1+\nu \ell )}{{\Gamma }_q\left(1+\nu \ell +(1-\mu )(1-\nu )\right)}{\left[(1-\mu )(1-\nu )+\nu \ell \right]}_q\hfill \\ \hfill & \times \frac{{\Gamma }_q((1-\mu )(1-\nu )+\nu \ell )}{{\Gamma }_q\left((1-\mu )(1-\nu )+\nu \ell +\mu (1-\nu )\right)}{w}^{(1-\mu )(1-\nu )+\nu \ell +\mu (1-\nu )-1}\hfill \\ \hfill & =\frac{{\Gamma }_q(1+\nu \ell )}{{\Gamma }_q\left(1-\nu +\nu \ell \right)}{w}^{\nu (\ell -1)}\hfill \end{array}$$

This leads directly to Equations (55) and (56). □

Theorem 5

(Kober q-integrals). If  ${\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)$ and  ${\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)$ are the generalized hyper q-Mittag-Leffler functions of Equations (33) and (34), respectively, then their Kober fractional q-integrals are

$$\begin{array}{c}\hfill I_q^{\nu ,\mu }\left\{{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)\right\}=\sum _{\ell =0}^{\infty }\frac{{\Gamma }_q(1+\mu +\ell ){\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Gamma }_q(1+\mu +\nu +\ell ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\ell }}{{(q;q)}_{\ell }},\end{array}$$

(57)

and

$$\begin{array}{c}\hfill I_q^{\nu ,\mu }\left\{{\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)\right\}=\sum _{\ell =0}^{\infty }\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}{\Gamma }_q(1+\mu +\ell ){\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Gamma }_q(1+\mu +\nu +\ell ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\ell }}{{(q;q)}_{\ell }}.\end{array}$$

(58)

where  $I_q^{\nu ,\mu }$ is the Kober fractional q-differential operator (Equation (16)).

Proof. 

By applying Equation (16) to Equation (33), it follows that

$$\begin{array}{cc}\hfill I_q^{\nu ,\mu }\left[{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)\right]& =\frac{w^{-\mu -\nu }}{{\Gamma }_q\left(\nu \right)}{\int }_0^w{(w-qt)}_q^{\nu -1}{t}^{\mu }{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)d_{q}t\hfill \\ \hfill & =\frac{w^{-\mu -\nu }}{{\Gamma }_q\left(\nu \right)}\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}{\int }_0^w{(w-qt)}_q^{\nu -1}{t}^{\mu +\ell }{d}_{q}t\hfill \\ \hfill & =\frac{w^{-\mu -\nu }}{{\Gamma }_q\left(\nu \right)}\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}\hfill \\ \hfill & \times (1-q)w\sum _{k=0}^{\infty }{q}^k{(w-q^{k+1}w)}_q{\left(q^{k}w\right)}^{\nu +\ell }\hfill \\ \hfill & =\frac{1}{{\Gamma }_q\left(\nu \right)}\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}\hfill \\ \hfill & \times (1-q)w^{\ell }\sum _{k=0}^{\infty }{q}^k{(1-q^{k+1})}_q{q}^{k(\nu +\ell )}\hfill \\ \hfill & =\frac{1}{{\Gamma }_q\left(\nu \right)}\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}{w}^{\ell }{\int }_0^1{(1-qt)}_q^{\nu -1}{t}^{\mu +\ell }{d}_{q}t\hfill \end{array}$$

By using Equation (10), one obtains

$$\begin{array}{cc}\hfill I_q^{\nu ,\mu }\left[{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)\right]& =\frac{1}{{\Gamma }_q\left(\nu \right)}\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}{B}_q(1+\ell +\mu ,\nu )w^{\ell }\hfill \\ \hfill & =\frac{1}{{\Gamma }_q\left(\nu \right)}\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}\frac{{\Gamma }_q(1+\ell +\mu )}{{\Gamma }_q(1+\mu +\nu +\ell )}{w}^{\ell }.\hfill \end{array}$$

This is the required result (Equation (57)). Similarly, Equation (58) can be proven. □

Theorem 6

(Riemann–Liouville q-integrals). If  ${\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)$ and  ${\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)$ are the generalized hyper q-Mittag-Leffler functions in Equations (33) and (34), respectively, then their Riemann–Liouville fractional q-integrals are

$$\begin{array}{c}\hfill I_q^{\nu }\left\{{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)\right\}=w^{\nu }\sum _{\ell =0}^{\infty }\frac{{\Gamma }_q(1+\ell ){\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Gamma }_q(1+\nu +\ell ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\ell }}{{(q;q)}_{\ell }},\end{array}$$

(59)

and

$$\begin{array}{c}\hfill I_q^{\nu }\left\{{\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q)\right\}=w^{\nu }\sum _{\ell =0}^{\infty }\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}{\Gamma }_q(1+\ell ){\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Gamma }_q(1+\nu +\ell ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{w^{\ell }}{{(q;q)}_{\ell }}.\end{array}$$

(60)

where  $I_q^{\nu }$ is the Riemann–Liouville fractional q-integral (Equation (14)).

Proof. 

The required results can be directly obtained by applying the fact that

$$\begin{array}{cc}\hfill I_q^{\nu }\left[w^{\ell }\right]& =\frac{1}{{\Gamma }_q\left(\nu \right)}{\int }_0^w{(w-qt)}_q^{\nu -1}{t}^{\ell }{d}_{q}t\hfill \\ \hfill & =\frac{w^{\nu +\ell }}{{\Gamma }_q\left(\nu \right)}{\int }_0^1{(1-qt)}_q^{\nu -1}{t}^{\ell }{d}_{q}t\hfill \\ \hfill & =\frac{{\Gamma }_q(1+\ell )}{{\Gamma }_q(1+\nu +\ell )}{w}^{\nu +\ell }.\hfill \end{array}$$

Therefore, the details are omitted. □

Remark 3.

Several special addenda of the outcomes in [21,22] can be obtained through Theorems 3–6. Furthermore, we can deduce many corollaries from Theorems 3–6 by using Remark 1.

5. Solutions to the Generalized Fractional  $\mathit{q}$-Kinetic Equations Pertaining to the Generalized Hyper  $\mathit{q}$-Mittag-Leffler Functions

In this section, we investigate some of the solutions to generalized fractional q-kinetic equations that encompass the generalized hyper q-Mittag-Leffler functions by applying the q-Sumudu transform.

First, we find the q-Sumudu transform of the generalized hyper q-Mittag-Leffler functions Equations (33) and (34) in Theorem 7 below.

Definition 14

(See [18]). A q-analogue of the Sumudu transform is defined by

$$\begin{array}{c}\hfill {\mathcal{S}}_q[f\left(t\right);s]=\frac{1}{(1-q)s}{\int }_0^s{E}_q\left(\frac{q}{s}t\right)f\left(x\right)d_{q}t,s>0.\end{array}$$

(61)

where  $E_q\left(w\right)$ is defined by Equation (20).

By taking the limit to be  $q\to 1$, Equation (61) leads to the well-known Sumudu transform [21,22] defined by

$$\begin{array}{c}\hfill \mathcal{S}[f\left(t\right);s]=\frac{1}{s}{\int }_0^s{e}^{\frac{-t}{s}}f\left(x\right)dt,s>0.\end{array}$$

(62)

and the q-Sumudu transforms for some functions that we need later are given as follows:

$$\begin{array}{c}\hfill {\mathcal{S}}_q[t^{\alpha -1};s]=s^{\alpha -1}{(1-q)}^{\alpha -1}{\Gamma }_q\left(\alpha \right),Re\alpha >0,\end{array}$$

(63)

and

$$\begin{array}{c}\hfill {\mathcal{S}}_q\left[\frac{1-q}{t}{E}_{\alpha ,0}\left(-c{t}^{\alpha };q\right);s\right]=\frac{1}{s\left(1+c{s}^{\alpha }{(1-q)}^{\alpha }\right)},\end{array}$$

(64)

where  $E_{\alpha ,0}$ is the q-Mittag-Leffler function defined by Equation (30) for  $\delta =0$.

Furthermore, the q-Sumudu transform for the Riemann–Liouville fractional q-integral operator in Equation (14) is

$$\begin{array}{c}\hfill {\mathcal{S}}_q\left[I_q^{\alpha }f\left(x\right);s\right]=s^{\alpha }{(1-q)}^{\alpha }{\mathcal{S}}_q[f\left(x\right);s],\alpha >0.\end{array}$$

(65)

Moreover, using Equation (63) leads to the following results:

Theorem 7.

The q-Sumudu transforms of Equations (33) and (34) are the following respective equations:

$$\begin{array}{c}\hfill S_q\left[{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q),w\to s\right]=\sum _{\ell =0}^{\infty }\frac{{\Gamma }_q(1+\ell ){\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{{(1-q)}^{\ell }{s}^{\ell }}{{(q;q)}_{\ell }},\end{array}$$

(66)

and

$$\begin{array}{c}\hfill S_q\left[{\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(w;q),w\to s\right]=\sum _{\ell =0}^{\infty }\frac{{(-1)}^{\ell }{q}^{\ell (\ell -1)/2}{\Gamma }_q(1+\ell ){\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{{(1-q)}^{\ell }{s}^{\ell }}{{(q;q)}_{\ell }}.\end{array}$$

(67)

where  $S_q$ is the q-Sumudu transform (Equation (61)).

Furthermore, the following lemma gives an application of the q-Sumudu transform pertaining to the fractional q-kinetic equations:

Lemma 1

(See [18]). The solution to the q-kinetic equation

$$\begin{array}{c}\hfill N_q\left(t\right)-N_0{f}_q\left(t\right)=-c{I}_q^{\alpha }{N}_q\left(t\right)\end{array}$$

(68)

via the q-Sumudu transform is

$$\begin{array}{c}\hfill N_q\left(t\right)=N_0{\int }_0^t\frac{1}{y}{f}_q(t-qy)E_{\alpha ,0}(-c{y}^{\alpha };q)d_{q}y,\end{array}$$

(69)

where  $E_{\alpha ,0}$ is defined by Equation (30).

Next, the main results of this section are as follows:

Theorem 8.

The solution to the q-kinetic equation

$$\begin{array}{c}\hfill N_q\left(t\right)-N_0{t}^{\beta -1}{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }\left(-c{t}^{\alpha };q\right)=-c{I}_q^{\alpha }{N}_q\left(t\right)\end{array}$$

(70)

is given by

$$\begin{array}{c}\hfill N_q\left(t\right)=N_0{t}^{\beta -1}\sum _{k=0}^{\infty }{(-c)}^k{t}^{\alpha k}{\mathfrak{E}}_{\alpha ,\beta +\alpha k,\delta }^{\lambda }\left(-c{t}^{\alpha };q\right)\end{array}$$

(71)

where  $c>0$ and  ${\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }(t;q)$ is the generalized hyper q-Mittag-Leffler function defined by Equation (33).

Proof. 

Setting

$$\begin{array}{cc}\hfill f_q\left(t\right)& =t^{\beta -1}{\mathfrak{E}}_{\alpha ,\beta ,\delta }^{\lambda }\left(-c{t}^{\alpha };q\right)\hfill \\ \hfill & =\sum _{\ell =0}^{\infty }\frac{{(-c)}^{\ell }{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}·\frac{t^{\alpha \ell +\beta -1}}{{(q;q)}_{\ell }},\hfill \end{array}$$

and

$$E_{\alpha ,0}(t;q)=\sum _{k=0}^{\infty }\frac{t^k}{{\Gamma }_q\left(k\alpha \right)},$$

in Equation (69) gives

$$\begin{array}{cc}\hfill N_q\left(t\right)=N_0& \sum _{\ell =0}^{\infty }\frac{{(-c)}^{\ell }{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}\sum _{k=0}^{\infty }\frac{{(-c)}^k}{{\Gamma }_q\left(k\alpha \right)}\times \hfill \\ \hfill & \times {\int }_0^t{(t-qy)}^{\alpha \ell +\beta -1}{y}^{\alpha k-1}{d}_{q}y,\hfill \end{array}$$

which is the solution to the q-kinetic equation (Equation (70)). Simplifying the calculation above using the convergence of the power series gives

$$\begin{array}{cc}\hfill N_q\left(t\right)& =N_0\sum _{\ell =0}^{\infty }\sum _{k=0}^{\infty }\frac{{(-c)}^{\ell +k}{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{(q;q)}_{\ell }{\Gamma }_q\left(k\alpha \right){\Gamma }_q(\ell \alpha +\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}\times \hfill \\ \hfill & \times \frac{{\Gamma }_q(\alpha \ell +\beta ){\Gamma }_q\left(\alpha k\right)}{{\Gamma }_q(\alpha \ell +\alpha k+\beta )}{t}^{\alpha \ell +\alpha k+\beta -1}\hfill \\ \hfill & =N_0{t}^{\beta -1}\sum _{k=0}^{\infty }{(-c)}^k{t}^{\alpha k}\sum _{\ell =0}^{\infty }\frac{{\Pi }_{s=1}^d{\Gamma }_q({\lambda }_s+\ell )}{{\Gamma }_q(\alpha \ell +\alpha k+\beta ){\Pi }_{r=1}^h{\Gamma }_q({\delta }_r+\ell )}\frac{{(-c{t}^{\alpha })}^{\ell }}{{(q;q)}_{\ell }}\hfill \\ \hfill & =N_0{t}^{\beta -1}\sum _{k=0}^{\infty }{(-c)}^k{t}^{\alpha k}{\mathfrak{E}}_{\alpha ,\beta +\alpha k,\delta }^{\lambda }\left(-c{t}^{\alpha };q\right),\hfill \end{array}$$

which is the required result (Equation (71)). □

In the same way as in the previous proof, the validity of the following results can be proven:

Theorem 9.

The solution to the q-kinetic equation

$$\begin{array}{c}\hfill N_q\left(t\right)-N_0{t}^{\beta -1}{\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }\left(-c{t}^{\alpha };q\right)=-c{I}_q^{\alpha }{N}_q\left(t\right)\end{array}$$

(72)

is given by

$$\begin{array}{c}\hfill N_q\left(t\right)=N_0{t}^{\beta -1}\sum _{k=0}^{\infty }{(-c)}^k{t}^{\alpha k}{\mathbb{E}}_{\alpha ,\beta +\alpha k,\delta }^{\lambda }\left(-c{t}^{\alpha };q\right)\end{array}$$

(73)

where  $c>0$ and  ${\mathbb{E}}_{\alpha ,\beta ,\delta }^{\lambda }(t;q)$ is the generalized hyper q-Mittag-Leffler function defined by Equation (34).

Corollary 1.

The solution to the q-kinetic equation

$$\begin{array}{c}\hfill N_q\left(t\right)-N_0{t}^{\beta -1}{\mathit{e}}_{\alpha ,\beta ,\delta }^{\lambda }\left(-c{t}^{\alpha };q\right)=-c{I}_q^{\alpha }{N}_q\left(t\right)\end{array}$$

(74)

is given by

$$\begin{array}{c}\hfill N_q\left(t\right)=N_0{t}^{\beta -1}\sum _{k=0}^{\infty }{(-c{t}^{\alpha })}^k{\mathit{e}}_{\alpha ,\beta +\alpha k,\delta }^{\lambda }\left(-c{t}^{\alpha };q\right)\end{array}$$

(75)

where  $c>0$ and  ${\mathit{e}}_{\alpha ,\beta ,\delta }^{\lambda }(t;q)$ is the generalized q-Mittag-Leffler function defined by (37).

Corollary 2.

The solution to the q-kinetic equation

$$\begin{array}{c}\hfill N_q\left(t\right)-N_0{t}^{\beta -1}{\mathit{E}}_{\alpha ,\beta ,\delta }^{\lambda }\left(-c{t}^{\alpha };q\right)=-c{I}_q^{\alpha }{N}_q\left(t\right)\end{array}$$

(76)

is given by

$$\begin{array}{c}\hfill N_q\left(t\right)=N_0{t}^{\beta -1}\sum _{k=0}^{\infty }{(-c)}^k{t}^{\alpha k}{\mathit{E}}_{\alpha ,\beta +\alpha k,\delta }^{\lambda }\left(-c{t}^{\alpha };q\right)\end{array}$$

(77)

where  $c>0$ and  ${\mathit{E}}_{\alpha ,\beta ,\delta }^{\lambda }(t;q)$ is the generalized q-Mittag-Leffler function defined by Equation (38).

By setting  $q\to 1$ in Theorems 8 and 9, we obtain a generalized fractional kinetic equation involving the generalized hyper Mittag-Leffler function in Equation (39). The solution can be obtained by using the Sumudu transform (Equation (62)). This result is given in the following corollary:

Corollary 3.

The solution to the kinetic equation

$$\begin{array}{c}\hfill N\left(t\right)-N_0{t}^{\beta -1}{\mathcal{E}}_{\alpha ,\beta ,\delta }^{\lambda }\left(-c{t}^{\alpha };q\right)=-c{I}^{\alpha }N\left(t\right)\end{array}$$

(78)

is given by

$$\begin{array}{c}\hfill N\left(t\right)=N_0{t}^{\beta -1}\sum _{k=0}^{\infty }{(-c)}^k{t}^{\alpha k}{\mathcal{E}}_{\alpha ,\beta +\alpha k,\delta }^{\lambda }\left(-c{t}^{\alpha }\right)\end{array}$$

(79)

where  $c>0$ and  ${\mathcal{E}}_{\alpha ,\beta ,\delta }^{\lambda }\left(t\right)$ is the generalized Mittag-Leffler function defined by Equation (39).

6. Graphical Representations

In this section, we offer a graphical elucidation of the determined results with different values for the parameters in Equations (71) and (73). Figure 1 and Figure 2 are the visual representations of Equations (71) and (73), respectively, using the parameters  $N_0=1,c=1,t\in [0,5]$, and  $\lambda =\delta .$ In Figure 1a and Figure 2a, we fix the values of q and  $\beta$ and drew graphs for different values of  $\alpha$. In Figure 1b and Figure 2b, we set the values of  $\alpha$ and q and drew plots for different values of  $\beta .$ In Figure 1c and Figure 2c, we fixed the values of  $\alpha$ and  $\beta$ and drew graphs for different q values. All the different values for the parameters were chosen according to the conditions of convergence of the functions in Equations (33) and (34). Furthermore, we could obtain several results from the graphical solutions to Equations (71) and (73) by selecting different values for  $\lambda ,\delta ,\alpha ,\beta ,$ and q that satisfied convergence of the functions in Equations (33) and (34).

7. Conclusions

Motivated by the established applications and future applications of fractional analysis in many branches of science, as well as the advent of a significantly larger spectrum of generalized Mittag-Leffler functions in statistics [30], physics [7], engineering [8], and other fields [31], this manuscript presents new forms of generalized hyper q-Mittag-Leffler functions in Equations (33) and (34). Then, we provided some of the convergence properties, fractional q-derivative, and q-integral operators of these functions. Furthermore, we investigated the solutions to the q-fractional kinetic equations in terms of the decided generalized hyper q-Mittag-Leffler functions by applying the defined q-Sumudu transform. Also, some exceptional cases for our principal outcomes were archived. In addition, we introduced graphical representations of the solutions, demonstrated in Section 5, using the program MATLAB.