Diverse Exact Soliton Solutions of the Time Fractional Clannish Random Walker ’ s Parabolic Equation via Dual Novel Techniques

In this article, we acquire a variety of new exact traveling wave solutions in the form of trigonometric, hyperbolic, and rational functions for the nonlinear time-fractional Clannish Random Walker ’ s Parabolic (CRWP) equation in the sense of beta-derivative by employing the two modi ﬁ ed methods, namely, modi ﬁ ed ð G ′ / G 2 Þ − expansion method and modi ﬁ ed F − expansion method. The obtained solutions are veri ﬁ ed for aforesaid equations through symbolic soft computations. To promote the essential propagated features, some investigated solutions are exhibited in the form of 2D and 3D graphics by passing on the precise values to the parameters under the constrain conditions. The obtained solutions show that the presented methods are e ﬀ ective, straight forward, and reliable as compared to other methods. These methods can also be used to extract the novel exact traveling wave solutions for solving any types of integer and fractional di ﬀ erential equations arising in mathematical physics.

The time-fractional Clannish Random Walker's Parabolic (CRWP) equation [31,32] is a model that can determine the behavior of two species A and B of random walker who execute a concurrent one-dimensional random walk characterized by an intensification of the clannishness of the members of one species A at point x at a time t, uðx , tÞ, can be expressed by the time-fractional CRWP equation as where α 1 is a parameter describing the order of the fractional time derivative and 0 < α 1 ≤ 1:. The major concern of this existing study is to utilize the novel meanings of fractional-order derivative, named beta fractional derivative [4], for time-fractional CRWP equation, and to find the novel comprehensive exact traveling wave solutions in the form of hyperbolic, trigonometric, and rational functions by employ two modified methods, modified ð G ′ /G 2 Þ − expansion method [33] and modified F − expansion method [34]. Beta-derivative has some interesting consequences in diverse areas including fluid mechanics, optical physics, chaos theory, biological models, disease analysis, and circuit analysis. To the best of our knowledge, the obtained solutions are more general and in different form which have never been reported in previously published studies [31,32]. Our results also enrich the variety of the dynamics of higher-dimensional nonlinear wave field. It is hoped that these results will provide some valuable information in the higher-dimensional nonlinear field.
By using the modified ðG′/G 2 Þ − expansion method [33], traveling wave solutions have been found for the nonlinear Schrödinger equation along third-order dispersion. Different types of traveling wave solutions of the Fokas-Lenells equations have been determined in [35] by this method. Alyahdaly found the general exact traveling wave solutions to the nonlinear evolution equations in [36]. Gepreel and Nofal [37] obtained the analytical solutions for nonlinear evolution equations in mathematical physics. Siddique and Mehdi found the exact traveling wave solutions for two prolific conformable M-fractional differential equations in [23]. Exact solutions for nonlinear integral member of Kadomstev-Petviashvili hierarchy differential equations have been determined by Gepreel [38].
A modified F-expansion method is proposed by taking full advantages of F expansion method and Riccati equation in seeking exact solutions of nonlinear PDEs. Darvishi and Najafi [39] used a modified F-expansion method to handle the foam Drainage equation. Aasaraai [40] used this method to construct new solutions of the nonlinear ð1 + 2Þ-dimensional Maccari's system. Aasaraai and Mehrlatifan [41] applied this method to coupled system of equation. Ali et al. [42] derived dispersive analytical soliton solutions of some nonlinear wave's dynamical model with the help of modified F-expansion method. Darvishi et al. [43] found traveling wave solutions for the ð3 + 1Þ-dimensional breaking soliton equation.
This article organized it as follows: in Section 2, we present beta-derivative and its properties. The descriptions of strategies are given in Section 3. In Section 4, we present a mathematical analysis of the models and its solutions via proposed methods. In Section 5, some graphical representations for some analytical solutions are presented. Some conclusions are drawn in the last section.

Beta-Derivative and Its Properties
Definition: suppose a function hðxÞ that is defined ∀ nonnegative x. Therefore, the beta-derivative of the function hð xÞ is given as [4]: Properties: assuming that a and b are real numbers, gðxÞ and hðxÞ are two functions β − differentiable and β ∈ ð0, 1, then, the following relations can be satisfied ii: iii: 3. Description of Strategies 3.1. The Modified ðG ′ /G 2 Þ-Expansion Method. Let us consider the nonlinear PDE is in the form where u = uðx, tÞ is an unknown function, and Q is a polynomial depending on uðx, tÞ and its various partial derivatives.
Step 1. By wave transformation Here, ν is the speed of traveling wave. The wave variable permits us to reduce Eq. (8) into a nonlinear ordinary differential equation (ODE) for U = Uð ηÞ: where R is a polynomial of UðηÞ and its total derivative with respect to η.

Journal of Function Spaces
Step 2. Extend the solutions of Eq. (9) in the following form where α i ði = 0, 1, 2, 3, ⋯, mÞ are constants and find to be later. It is important that α i ≠ 0: Integer m can be determined by considering the homogenous balance between the governing nonlinear terms and the highest order derivatives in Eq. (9). The function G = GðηÞ satisfies the following Riccati equation, where λ 0 and λ 1 are constants. We gain the below solutions to Eq. (11) due to different conditions of λ 0 : When λ 0 λ 1 < 0, When λ 0 λ 1 > 0, When λ 0 = 0 and λ 1 ≠ 0, where C 1 and C 2 are arbitrary constant.
Step 4. By substituting Eq. (10) of which α i , ν and other parameters that are found in step 3 into Eq. (8), we get the solutions of Eq. (7).

The Modified
Here, we will describe the basic steps of F − expansion method [34].
where a 0 ,a i , and b i are constants to be determined. FðηÞ sat-isfies the Riccati equation: where A,B, and C are constants to be determined. The prime denotes d/dη. Integer m can be determined by considering the homogenous balance between the governing nonlinear terms and the highest order derivatives of UðηÞ in Eq. (9).
Given different values of A, B, and C, the different Riccati function solution FðηÞ can be obtained from Eq. (16) (see Table 1).
Step 2. Extend the solution of Eq. (9) in the following form Step 3. Substituting Eq. (15) along with Eq. (16) into Eq. (9) and collect coefficients of F i ðηÞ to zero yields a system of algebraic equations for a i and b i : Step 4. Solve the system of algebraic equations, probably with the aid of Mathematica. a i and b i can be expressed by A, B, and C (or the coefficients of Eq. (9)). Substituting these results into (16), we can obtain the general form of traveling wave solutions to Eq. (9).
Step 5. Selecting A, B, C, and FðηÞ from Table 1 and substituting them along with a i and b i into Eq. (15), a series of soliton-like solutions, trigonmetric function solutions, and rational solutions to Eq. (7) can be obtained. The modified F-expansion method is more effective in obtaining the soliton-like solution, trigonometric function solutions, exponential solutions, and rational solutions of the nonlinear partial deferential equations. This method will yield more rich types solutions of the nonlinear partial deferential equations. It shows that the modified F-expansion method is more powerful in constructing exact solutions of nonlinear PDEs.

Application
Time-fractional Clannish Random Walker's Parabolic equation: Let us assume the transformations: where c is constant. By using Eq. (17) into Eq. (1), we get the following ordinary differential equation.
In the following subsections, the proposed methods are applied to extract the required solutions: where a 0 and a 1 are unknown parameters. By using Eq. (19) along with Eq. (11) into Eq. (18) and summing up all coefficients of same order of ðG ′ /G 2 Þ, we get the set of algebraic equations involving a 0 , a 1 and other parameters. Solving the obtained set of algebraic equations with Mathematica, we reach the following results: Now we use the Eqs. (19) and (12)- (14) into Eq. (19) and set the below cases.
if λ 0 λ 1 < 0, then If λ 0 λ 1 > 0,then If λ 0 = 0, λ 1 ≠ 0,then where a 0 and a 1 are unknown parameters. By using Eq. (24) along with Eq. (16) into Eq. (18) and summing up all the coefficients of same order of F, we get the set of algebraic equations involving a 0 , a 1 and other parameters. Solving the obtained set of algebraic equations with Mathematica, we reach the following results.  [34].

Journal of Function Spaces
Family-III For Family-II  Family-III For Family-II Family-III For Family-II Family-III For A = 0, B = 0.

Results and Discussion
In this section, we discuss the graphical interpretation of obtained results. Two powerful analytical methods, namely, modified ðG ′ /G 2 Þ − expansion method and modified F − expansion method, are used to extract the trigonometric, hyperbolic, and rational wave solutions of the governing model. The physical significance of these solutions is shown by assigning particular values of free parameters. The solutions to Eqs.
, (32), and (34). Finally, we explain the dynamic performance of the trigonometric function answers of Eqs. (22) and (34) in Figures 7 and 8, which depict the 3D shape and 2D graph for different values of the fractional parameter β for the trigonometric function answers of Eqs. (22) and (34). The implemented mathematical simulations acknowledge that the answers are of periodic wave shapes and of rational, hyperbolic, and trigonometric categorizations. Furthermore, through observing the construction of the acquired solutions, it could be understood that the parameter β of fractional derivatives has important role in the formulation of all the solutions.

Conclusions
In this work, we applied the modified ðG′/G 2 Þ-expansion method and modified F − expansion method in a satisfactory way to find the novel exact traveling wave solutions of the time-fractional CRWP equation in the sense of betaderivative. Various obtained solutions are in the form of hyperbolic, trigonometric, and rational forms. To describe the physical phenomena of the time-fractional CRWP model, some solutions are plotted in the form of 2D and 3D by assigning the specific value to the parameters under the constrain conditions. All algebraic computations and graphical representations in this work have been provided for the obtained solutions at various parameters values with the help of Mathematica. It is essential to note that these new solutions of the time-fractional CRWP equation have not been exposed in literature by employed our two analytical modified mathematical methods. Lastly, the studied methods can be potentially applied to solve various nonlinear PDEs that are apparent in many important nonlinear scientific phenomena in physics and engineering.

Data Availability
No data were used to support this study.

Conflicts of Interest
The authors declare that they have no conflicts of interest.