Fractional q -Calculus Operators Pertaining to the q -Analogue of M -Function and Its Application to Fractional q -Kinetic Equations

: In this article, the authors introduce the q -analogue of the M -function, and establish four theorems related to the Riemann–Liouville fractional q -calculus operators pertaining to the newly defined q -analogue of M -functions. In addition, to establish the solution of fractional q -kinetic equations involving the q -analogue of the M -function, we apply the technique of the q -Laplace transform and the q -Sumudu transform and its inverse to obtain the solution in closed form. Due to the general nature of the q -calculus operators and defined functions, a variety of results involving special functions can only be obtained by setting the parameters appropriately.


Introduction
Fractional differential equations have a remarkably significant role in applied sciences, especially in dynamic systems, mathematical physics, control systems, and engineering, and they are used to build mathematical models of various physical phenomena.Solving fractional differential equations can be challenging because fractional derivatives are nonlocal operators, and analytical solutions are often not readily available.Numerical methods and special functions, like the Mittag-Leffler function, are commonly used to approximate solutions to fractional differential equations.Fractional calculus, which deals with fractional integrals and derivatives, provides a mathematical framework for studying and solving fractional differential equations.Sharma and Jain [1] recently established the generalized M-series as a function and stated it, using power series for z, ς, ϑ ∈ C, ℜ(ς) > 0, as: where (c i ) n , d j n ; i = 1, 2, . . ., r; j = 1, 2, . . ., s are known as Pochhammer symbols.The M-functions have some special cases, for example, if we set ϑ = 1, the M-function reduces in the M-series defined by Sharma [2], for ς = ϑ = 1, and the generalized M-function is the well-known generalized hypergeometric function.Also, when ς = ϑ = 1, and r = s = 0, then we obtain the Mittag-Leffler function [3] with one parameter, for r = s = 0, and after that we have a Mittag-Leffler function of two indexes established by Wiman [4].Several characteristics and applications of the M-series (1) and its particular cases have been studied by Miller [5], Saxena et al. [6], Purohit et al. [7], Mishra et al. [8], and Saigo and Kilbas [9].The generalized M-function is crucial in the solution of fractional order integral and differential equations.
Researchers continue to develop new techniques and tools to better understand and solve fractional differential equations due to their importance in complex physical and biological systems.More specifically, kinetic equations are the fundamental equations of natural sciences and mathematical physics because they define the continuity of the motion of the matter.They can be applied to qualitative and quantitative descriptions of biological, chemical, physical, social, and other processes.They are called "master equations" since the kinetic equation-based mathematical simulation of evolutionary processes is productive and powerful.The system of Vlasov's equation, Liouville equation, Fokker-Planck radiation transfer equations, Vegner equations, the Boltzmann equation, and the chain of Bogoljubov's equations for the matrix of density and diffraction are significant kinetic equations.In 2000, Haubold and Mathai [10] developed a fractional differential equation involving the reaction's rate of change, the destruction rates, and the production rate.We discovered several papers in the literature that investigate extensions and generalizations of these equations utilizing different fractional calculus operators, for example, the work of , Kumar et al. [14,15], Suthar et al. [16,17], Habenom et al. [18], Baleanu et al. [19], and Gupta and Parihar [20].
As a result of the applications of q-calculus in physics, statistics, and mathematics, there has recently been a significant increase in activity in the q-calculus area of research.Particularly in partition theory, exactly solvable models in statistical mechanics, combinatorics, optimal control problems, computer algebra, geometric function theory, q-transform analysis, q-integral equations, q-analysis, and q-difference have found many useful applications.Discussions of q-analysis can be found in the works of Annaby and Mansour [21], Kac and Cheung [22], and Gasper and Rahman [23].For more information, we refer to Abdi [24] and Al-Omari et al. [25] for applications of the q-Laplace transform to solve the q-difference equations.Albayrak et al. [26] established the q-analogue of Sumudu transforms and derived some fundamental characteristics of q-Sumudu transforms such as differentiation, shifting theorems, linearity, integration, and some fascinating connection theorems that involve the q-Laplace transform.
Garg and Chanchlani [27] recently computed the following fractional q-kinetic equation in closed form using q-Laplace transforms and its inverse: Moreover, Purohit and Ucar [28] have proposed a new solution to the q-kinetic equation with the Riemann-Liouville fractional (RLF) q-integral operator.Using the inverse q-Sumudu transform, the solution is found using the q-Mittag-Leffler function.
Motivated by the above series representation which is defined in Equation ( 1) and its importance and useful applications in the field of applied mathematics, we established the q-analogue of the generalized M-functions and determined four theorems on the RLF q-calculus operators pertaining to the above mentioned function.Furthermore, we calculate the solution of the fractional q-kinetic Equation (2) involving the RLF q-integral operator by using the q-Laplace and Sumudu transform.

Preliminaries
For a complex or real m, the q-shifted factorial is given by in [23] as: The q-beta and q-gamma functions [23] are defined as and An analogue of q of the power function (x − τ) m is given by: Moreover, Rajkovic et al.
[29] stated the q-derivative for function f (x) as and A q-extension of an exponential function is stated as 3. q-Analogue of M-Function Motivated by Equation ( 1) and their related literature, we introduce the q-analogue of the M-function where Γ q (.) is the q-gamma function and (q c i , q) υ , (q converges for all z if r ≤ s, and for |z| < We describe the relationships as specific instances of ς,ϑ r M s ( q, z) with other special functions as given below: For r = s = 0, the M-function becomes the q-Mittag-Leffler function defined by Mansoor [30] as: When ϑ = 1, r = s = 0, the M-function becomes the q-Mittag-Leffler function defined by Jain [31] as: (14), we obtain the generalized small q-Mittag-Leffler introduced by Purohit and Kalla [32] as follows: On putting r = s = 1, c 1 = γ , d 1 = δ in ( 14), we obtain the generalized q-Mittag-Leffler introduced by Sharma and Jain [33] as follows: Putting r = 1, s = 0, c 1 = δ in ( 14), we obtain where the function 1 ϕ 0 (q δ ; −; q, z) = (1 − z) −δ can be called the q-binomial function.
Finally, in view of a relations lim and lim Furthermore, note that lim

Fractional q-Integration and q-Differentiation Approach
To begin, we would like to define these q-fractional operators in this section.
Agarwal [34] introduced a left-sided Riemann-Liouville (R-L) q-fractional operator for a = 0 as follows: Furthermore, Mansour [35] stated a right-sided R-L q-fractional operator by: From [35], for ε > 0 and ⌈ε⌉ = m, the right-and left-sided R-L fractional q-derivatives of order ε are given as and The left-and right-sided Caputo fractional q-derivatives of order ς are given as This part would determine the following interesting results in the form of theorems.Here, we present q-analogues of M-functions in view of the fractional q-integration and q-differentiation representations.

Generalized Fractional Kinetic Equation
The fractional order of differential equations have recently proven to be useful modeling instruments in a wide range of physical phenomena.This is due to the fact that the feasible modeling of a physical phenomenon depends on the instant time, which is additionally achieved effectively by using fractional calculus.Their treatment from the point of view of the q-calculus by using fractional q-calculus operators of the R-L type can also open up new avenues.
If p is the production rate on N, let N(t) be an arbitrary time-dependent reaction, d be the destruction rate, and then a mathematical description of the three quantities is as follows: A differential fractional equation for the quantities N(t), p, and d, as introduced by Haubold and Mathai [10], is as follows: where N t (t * ) = N(t − t * ) for t * > 0. In addition, Haubold and Mathai [10] discovered that when the spatial shifts or inhomogeneities in quantities N(t) are ignored, Equation (48) becomes: For the starting situation, let N i (t = 0) = N 0 be a number density of species i at a time t = 0; the constant c i > 0, is called the standard kinetic equation.So, the solution of Equation ( 49) is provided by: Or else, when the index i is removed and the standard kinetic equation is integrated (49), we obtain: where 0 D −1 t is the standard operator of the integral.Haubold and Mathai [10] stated the fractional generalization of Equation (51) by: The authors' q-Laplace transform of the R-L fractional q-integral operator in [27] is as follows: where the R-L fractional q-integral operator [37] is provided as Using Equation (62), the power function of q-Laplace transform and q L κ N q (t) = N * q (κ) in Equation ( 61), above, we obtain After rearrangement of the terms, we use the geometric series expansion formula, and we obtain Using the inverse q-Laplace transform on both sides of (64), we have We obtain the required result by interpreting the result in (65) in light of (15), that is, ) υ e ς,ςυ+ϑ (−ωt ς ; q).

Alternative Solution of Generalized Fractional Kinetic Equations by Using the Sumudu Transform
In this section, we look at how a q-Sumudu transform can be used to solve a generalized fractional q-kinetic equation involving the q-analogue of the M-functions (14).
Albayrak et al. [26] established and studied a q-Sumudu transform to aid in the process of solving integral and differential equations in the domain time, as well as for use in different applications of systems in applied physics and engineering.The q-Sumudu transform has very unique and helpful characteristics, and it can be used to solve problems governing kinetic equations in science and engineering.The q-Sumudu transforms are the theoretical counterpart of the q-Laplace transforms.
The primary type of q-analogue of the Sumudu transform defined and represented by Albayrak et al. [26] (see also [28]) is as follows: When applied cross the set of functions, we have The following are the q-Sumudu transforms of some unique functions used in the sequel: S q (1 − q) x E ς,0 (−ωx ς ; q); σ = 1 Albayrak et al. [26] provided a q-Sumudu convolution theorem as follows: S q f * q g (x); σ = σ S q { f (x); σ}S q {g(x); σ}, (70) when f * q g are the q-convolutions of two analytical functions f (t) and g(t), provided by where Discussion 1.Let 0 < |q| < 1, ω > 0, ς > 0; the solution to the generalized fractional q-kinetic Equation (58) is provided by (59).
Using the q-Sumudu applied to both sides of Equation (58), it becomes The q-Sumudu transform of the R-L fractional q-integral operator is defined by Purohit and Ucar [28] as: where the R-L fractional q-integral operator was established by Al-Salam [37] and was provided independently by Agarwal [34] as: Now, using Equations ( 14) and (76) in Equation (72), we have the following: Substituting S q N q (t); σ = N(σ), we have Using Equation (68), we obtain After rearrangement of the terms, we use a geometric series expansion formula and obtain Using the inverse q-Sumudu transform with both sides of (73), we have We easily arrive at the desired result (59) by interpreting the result (77) in the context of (15).
Furthermore, by focusing the variables involved in the main theorem, we can solve some other fractional q-kinetic equations and also the ordinary q-kinetic equations.
Corollary 13.Let ω > 0, ς > 0, 0 < |q| < 1; then the following fractional q-kinetic equation N q (t) − N 0 t ϑ−1 E δ ς,ϑ (q, −ωt ς ) = −ω I ς q N q (t) , (90) has the following solution N q (t) = N 0 t ϑ−1 ∞ ∑ υ=0 q δ ; q υ (q; q) υ (−ωt ς ) υ e ς,ςυ+ϑ (−ωt ς ; q).(91) If we take q → 1 in our Theorem 5 and using Equations ( 15)-( 17), we also determine that the findings presented here are general enough to produce solutions to a lot of known or novel fractional kinetic equations involving such other unique functions as (for example) those examined by Saxena et al. [12], as well as solutions to the problems examined by Saxena et al. [11,13] and Haubold and Mathai [10].Also, we conclude by noting that the q-Sumudu transform method is an extremely successful and practical alternative method for solving fractional q-differential equations.The extended kinetic equations are expected to have potential applications in nuclear power, nuclear physics, astrophysics, and other related fields.