New Numerical Results for the Time-Fractional Phi-Four Equation Using a Novel Analytical Approach

: This manuscript investigates the fractional Phi-four equation by using q -homotopy analysis transform method ( q -HATM) numerically. The Phi-four equation is obtained from one of the special cases of the Klein-Gordon model. Moreover, it is used to model the kink and anti-kink solitary wave interactions arising in nuclear particle physics and biological structures for the last several decades. The proposed technique is composed of Laplace transform and q -homotopy analysis techniques, and fractional derivative deﬁned in the sense of Caputo. For the governing fractional-order model, the Banach’s ﬁxed point hypothesis is studied to establish the existence and uniqueness of the achieved solution. To illustrate and validate the e ﬀ ectiveness of the projected algorithm, we analyze the considered model in terms of arbitrary order with two distinct cases and also introduce corresponding numerical simulation. Moreover, the physical behaviors of the obtained solutions with respect to fractional-order are presented via various simulations.


Introduction
Recently, many new fractional-order operator types have been introduced to the literature by some scholars, and established its fundamental properties [1][2][3][4][5][6]. Fractional theory is widely used to explain many properties of phenomena such as nanotechnology [7], optics [8], human diseases [9], chaos theory [10], and others . As an example, the tumor-immune surveillance model has been investigated in [33]. Some important properties of dengue fever have been effectively investigated in [34,35]. The analytical as well as numerical solutions for the equations illustrating the above cited phenomena, have an important role in describing the behavior of nonlinear models ascends [36][37][38][39][40]. For the differential equation, symmetry is a transformation that keeps its family of solutions invariant, and further, its analysis can be applied to examine and illustrate various classes of differential equations. The study of fractional calculus associated to symmetry recently attracted many researchers from different disciplines in order to present their viewpoints while analyzing real-word problems.
Here, λ 1 and λ 2 are real-valued constants. Here, v(x, t) is the normalized propagation of distance and retracted time t.
Many physicists and mathematicians recently proposed very accurate and more effective algorithms to find and examine the solution for equations describing nonlinear mechanisms arising in science and engineering. As an approximate method, the homotopy analysis method (HAM) has been introduced by L. Shijun [43,44], which is formed to describe deformation from zero to 1. Recently, it has been efficiently and beneficially employed to examine the nature of nonlinear problems without linearization and perturbation. However, HAM necessitates huge time and computer memory for computational work. Therefore, there is an essence of the consolidation of this scheme with familiar transform algorithms. The projected solution procedure is the mixture of q-HAM with LT [45]. Since the future scheme is a modified method of HAM, it does not kernel discretization, perturbation or linearization [46][47][48][49][50][51][52][53][54][55][56].
The considered nonlinear problem recently fascinated authors from various areas of science. Due to the numerous applications of the proposed model and the important role it plays in describing various nonlinear phenomena, several researchers found and investigated the solution numerically and also analytically; for instance, researchers in [57] find the new exact solution for the fractional case of considered nonlinear problem, and the spectral collocation method is considered in [58] to find the solution for PF equation. The considered nonlinear PF equation has attracted the attention of researchers, and, in order to present their viewpoint, many techniques are considered [59][60][61][62][63][64][65], presenting some interesting results. Recently many researchers from all over the world have investigated Fractional Calculus (FC) as an effective tool by comparing with other operators some important physical problems [66][67][68][69][70][71][72][73][74][75][76][77][78][79]. For instance, authors in [80] presented the reflection symmetry in fractional calculus, and proposed some interesting properties and corresponding consequences.
The main aim of this paper is to investigate the numerical solutions for the considered fractional Phi-Four equation using the novel technique. The proposed algorithm provides more liberty to consider the diverse class complex as well as nonlinear problems and the initial guess.
In this paper, we try to capture in the physical behaviors obtained numerical solutions with distinct arbitrary order and the parameters accessible by the projected algorithm. Moreover, in order to present accuracy and effectiveness, we present the numerical simulations. By the help of these results, we try to explain the diverse model exemplifying numerous phenomena arising in daily life. The rest of the paper is organized as follows: the basic properties of FC are presented in Section 2, the basic solution algorithm of the projected method for the considered problem is defined in Section 3, the convergence analysis of the considered scheme is illustrated in Section 4, the solution for FPF equation and its corresponding consequences in terms plots and also numerical simulation are respectively cited in Sections 5 and 6. Finally, some important findings of the paper are presented as a conclusion section.

Preliminaries
In this subsection of the paper, it is presented the basic definitions of fractional calculus. (2)

Definition 2.
Fractional derivative of f ∈ C n −1 is given in the sense of Caputo as the following: where V(x, s) is LT of v(x, t).

Solution Procedure for Fractional Phi-Four Equation
Here, we consider the nonlinear FPF equation to illustrate the basic solution algorithm of the projected method with associated initial conditions as follows: and v(x, 0) = g 1 (x) and v t (x, 0) = g 2 (x).
Here, D µ t v(x, t) signifies the fractional Caputo-derivative of v(x, t). Here, v(x, t) is a bounded function. On using the LT on Equation (5), and then by the help of Equation (6), it is given as following By the assist of Equation (7), the nonlinear operator is given as follows: Now, the homotopy structure is defined by where q ∈ 0, 1 n (n ≥ 1) and 0 are respectively the embedding and auxiliary parameter. Further, these parameters help us to control and adjust the convergence region of the obtained solution. For the proper choice of n and , the obtained solution quickly converges to the exact solution in an admissible domain. The following results are respectively true for q = 0 and q = 1 Near to q, we define ϕ(x, t; q) in series form by help of the Taylor theorem and then one can obtain where For the proper chaise of v 0 (x, t), n and ; the series (11) converges at q = 1 n . Then Dividing by m! after m-times differentiating Equation (9) with q and then putting q = 0, it simplifies to where k m and v m are defined by and On applying inverse LT on Equation (14), produces With the help of Equations (17) and (18) Using Equation (19), we achieve the series of v m (x, t). Lastly, the series q-HATM solution is defined as

Convergence Analysis of the Proposed Method
In this part of the paper, for the considered problem, we introduce the convergence analysis as follows. then there is a fixed point for F with the help of Banach's fixed point theory [81][82][83]. Then, for u 0 , v 0 ∈ E the sequence obtained by q -HATM solution converges to a specified point of F, and Proof. Let us consider a Banach space (C[J], . ) with the norm defined as g(t) = max t∈J g(t) for all For LT, with the help of the convolution theorem, we get Then, the above inequality reduces to . It completes our required proof.  (5), the q -HATM solution is unique wherever 0 < α < 1, where α = (k m + ) + δ 2 + λ 1 + λ 2 P 2 + Q 2 + PQ T . (5), the solution is illustrated as following

Proof. For fractional PF equation cited in Equation
Suppose v and v be the two different solutions for Equation (5) such that |v| ≤ P and v ≤ Q, then using the above relation, it is obtained as following For Laplace transform, we have, by the help of the convolution theorem, Since 0 < α < 1, therefore v − v = 0, which gives v = v . This completes the required proof.

Solution for Fractional PF Equation
In order to validate the applicability and efficiency of the future technique, here we consider the two distinct examples.

Example 2. Consider the FPF equation at
subjected to the initial conditions v(x, 0) = tan h By using the proposed algorithm and Equation (28), the Equation (27) becomes Symmetry 2020, 12, 478 10 of 16 We evaluate the iterative terms of v m (x, t) with the assistance of Equation (29), and as follows Now, the series form of q-HATM solution is as follows:  Table 2).

Numerical Results and Discussion
Here, we have demonstrated the numerical simulations of the Phi-four equation having fractional-order by using q-HATM. The fourth order solution (i.e., up to ( , )) we consider to illustrate the behavior of the obtained solution in terms of plots and tables. In Table 1 and Table 2, the numerical studies have been conducted to ensure the exactness of the projected algorithm, and further from the tables, it is confirmed that the considered model conspicuously depends on the time. The surfaces of the q-HATM solution, exact solution and absolute error are captured in Figure 1 for

Numerical Results and Discussion
Here, we have demonstrated the numerical simulations of the Phi-four equation having fractional-order by using q-HATM. The fourth order solution (i.e., up to v 4 (x, t)) we consider to illustrate the behavior of the obtained solution in terms of plots and tables. In Tables 1 and 2, the numerical studies have been conducted to ensure the exactness of the projected algorithm, and further from the tables, it is confirmed that the considered model conspicuously depends on the time. The surfaces of the q-HATM solution, exact solution and absolute error are captured in Figure 1 for the FPF equation considered in Example 1. To show the projected problem remarkably depends on arbitrary order and the response of q-HATM solution for different µ, the plots have been drowned and are cited in Figures 2 and 5 for Examples 1 and 2., respectively. Similarly, nature obtained and exact solutions in association with absolute error for Example 2, are presented in Figure 4. The -curves have been drowned to analyze the behavior of the achieved solution related homotopy parameter with different µ for both examples and have been respectively captured in Figures 3 and 6. These curves can help us to adjust and regulate the region of convergence for the q-HATM solution. For unsuitable , the q-HATM solution swiftly converges to exact solution. The demonstrated plots help us to better understand the nature of the fractional Phi-four equation when temporal-spatial variables vary in comparison with arbitrary order.

Conclusions
In the present framework, q-HATM has been employed to find the numerical solutions for the fractional Phi-four equation. With the help of Banach's fixed point theory, the convergence analysis projecting the nonlinear problem has been presented. To present the efficiency as well as the applicability of the projected algorithm, we have considered two distinct cases. The present study confirms that the projected nonlinear problem is remarkable with the time instant as well as the time history and these can be effectively exemplified by employing the concept of fractional calculus. The governing model plays a vibrant role while analyzing many physical phenomena, and thus, for future work, this can be examined by using recently proposed and nurtured effective and accurate methods [87,88]. With the help of the obtained results, we can capture more interesting consequences. Finally, it has been observed that the considered scheme is highly effective, more accurate and extremely methodical, and this can be employed to exemplify the various classes of nonlinear models that exist in science and technology. The results obtained demonstrate that the considered method is more effective and easy to employ to scrutinize the behaviors of fractional differential equations with multi-dimensions arising in associated areas of science and technology.