P.Veeresha,Esin Ilhan,D.G.Prakasha,Haci Mehmet Baskonus and Wei Gao

1Department of Mathematics,CHRIST(Deemed to be University),Bengaluru,560029,India

2Faculty of Engineering and Architecture,Kirsehir Ahi Evran University,Kirsehir,40500,Turkey

3Department of Mathematics,Faculty of Science,Davangere University,Shivagangothri,Davangere,577007,India

4Department of Mathematics and Science Education,Faculty of Education,Harran University,Sanliurfa,63050,Turkey

5School of Information Science and Technology,Yunnan Normal University,Kunming,650500,China

ABSTRACT In this paper, we analyze the behaviour of solution for the system exemplifying model of tumour invasion and metastasis by the help of q-homotopy analysis transform method (q-HATM) with the fractional operator.The analyzed model consists of a system of three nonlinear differential equations elucidating the activation and the migratory response of the degradation of the matrix, tumour cells and production of degradative enzymes by the tumour cells.The considered method is graceful amalgamations of q-homotopy analysis technique with Laplace transform(LT),and Caputo-Fabrizio(CF)fractional operator is hired in the present study.By using the fixed point theory, existence and uniqueness are demonstrated.To validate and present the effectiveness of the considered algorithm,we analyzed the considered system in terms of fractional order with time and space.The error analysis of the considered scheme is illustrated.The variations with small change time with respect to achieved results are effectively captured in plots.The obtained results confirm that the considered method is very efficient and highly methodical to analyze the behaviors of the system of fractional order differential equations.

KEYWORDS Tumour cell; invasion and metastasis; q-homotopy analysis transform method; Caputo-Fabrizio derivative

1 Introduction

The existence of cancer among the people around the globe has been developed swiftly, and globally it becomes the second foremost cause of demise after cardiovascular diseases [1].The process of spreading and formation of secondary tumours is known as Metastasis and this nature of cancer cells is the main reason for the death in cancer patients.In addition, the prediction of size, stage, and evolution of a tumour is very critical for the treatment of cancer.Moreover,mathematics plays as an essential tool and aid us to analyze the behaviour of the tumour.The tumour growth has been mathematically modelled by the number of researchers and which are appeared in the literature [2-4].The growth of the tumour invasion and metastasis are described by PDEs.Particularly, deterministic diffusion-reaction equations and these equations are employed to model the spatial spread of tumours at the initial development and later invasive moments [5].

The fractional-order derivatives are introduced by Leibnitz soon after the classical order derivatives.As compared to classical calculus, it was soon discovered that fractional calculus (FC)is more appropriate for describing real-world problems [6-10].The calculus of arbitrary order turned out one of the most essential tools to describe biological phenomena.The human diseases which are modelled through derivative having fractional-order help us to incorporate the information about its present and past states [11-16].Moreover, these models demonstrate the non-local distributed effects, hereditary properties and system memory.These properties are necessary to describe the biological models.In connection with this, recently many authors established the arbitrary-order model to analyze the diffusion equation and to forecast the effect of the tumour and they applied many powerful methods to find the solution for these models [17-32].The pivotal aim of generalizing the integer to fractional order is to capture consequences related non-locality,long-range memory and time-based properties and also anomalous diffusion aspects.Many realworld problems exemplified with complex and nonlinear models are effectively, systematically and accurately illustrated and investigated by the aid of theory and fundamental concept of FC.Many pioneers nurtured and developed novel and distinct notions of fractional order for both differential and integral operators.Most familiarly hired operators to analyze many models are Riemann,Liouville, Caputo, Fabrizio and others.However, researchers pointed out some limitations while generalizing the system with these notions.The Riemann-Liouville derivative fails to elucidate the essence of initial conditions; the Caputo derivative has overcome this limitation and latter it has been widely applied to the numerous classes of mathematical models exemplifying real-world problems.But this fractional operator is unable to describe the singular kernel of the systems or problems.However, Caputo et al.[33] in 2015 overcome the foregoing obliges and then the number of authors employed CF derivative to investigate and study wide classes of complex and nonlinear problems.It has been proved by many researchers; CF fractional operator as great results compared to other fractional operators.

The study of mathematical models effectively exemplified diverse phenomena.However, as much as important of nurturing the phenomena with the system of equations finding the corresponding solution is also very vital and difficult.In this path, many authors examined diverse phenomena, for instance, the structured predator-prey model with prey refuge [34],COVID-19 [35-39], Zika virus transmission [40], planar system-masses in an equilateral triangle [41], a harmonic oscillator with position-dependent mass [42], time fractional Burgers equation [43], fractional optimal control problems [44], Emden-Flower type equations [45], and many others [46-53].These investigations help researchers to understand the importance of generalizing the classical concept into fractional operators, and efficiency and difference between diverse schemes.

Many physicists, engineers and mathematicians recently proposed and modified diverse solution procedure with a different approach with respect to increasing in accuracy and methodology,to reduces the complexity, many additional assumptions and consideration, huge time for evaluation and to save computer memory.Moreover, each method is suitable for some specific family of problems and they have their own limitations, including conversion of nonlinear to linear,partial to ordinary differential equations, splitting complex and nonlinear term to simple parts terms.In this connection, with the help of topological concept called homotopy,Liao Shijunwho is a Chinese Mathematician proposed algorithm called homotopy analysis method (HAM) and illustrated to confirm it overcomes almost all the limitations raised while we solving nonlinear systems exists in sciences and other disciplines associated to mathematics [54].The most familiar thing of employing this method by many authors is including it solves nonlinear problems without linearization and perturbation.

As science and technology-enhanced, mankind always expecting new tools or modifications in existing tools to improve the accuracy and reduces the time taken for finding needful.In this regard, some scholars pointed out similar things in HAM and suggested to union with existing and familiar transformation.Authors in [55] modifiedq-HAM with the help of Laplace transform(LT) and manifest new modified scheme is calledq-HATM.This method is perceptible includes all merits which are achieved by HAM and also it attracted many researchers to analyze the diverse class of models and systems.For instance, the model exemplifying three Lakes pollution with the newly proposed fractional operator is investigated by authors in [56], fractional vibration equation is analyzed by authors in [57] with some interesting results, authors in [58] presented the efficiency of the projected scheme while analyzing Swift-Hohenberg equation having arbitrary order, the accuracy of the hired scheme in comparison with existing results is illustrated by authors in [59] with respect to the physical model, the convergence analysis of the considered method for Lienard’s equation is demonstrated in [60], many others analyzed various biological and physical phenomena by the assist of the projected scheme [61,62].

In the present study, we find the series solution for a system of nonlinear differential equations describing the model of tumour invasion and metastasis usingq-HATM with the help of a novel fractional operator.By using the important results of fixed-point theory for the projected system the existence and uniqueness are demonstrated.The novelty of the projected scheme gives more freedom to choose the initial conditions and the novelty is it offers a simple solution procedure and associated with parameters to provide the swift convergence.Further, it contains the results achieved by other classical methods including ADM, HPM,q-HAM and some other methods [63-72].In the present study, we analyzed the system describing the tumour invasion and metastasis with different time and space for different fractional-order usingq-HATM within the frame of the novel fractional operator which can describe the singular kernel.This study can help us to analyze more complex and nonlinear mathematical models describing the deadly virus or diseases.

The rest of the manuscript is organized as follows: The basic and fundamentals are presented in the next section, the hired model is exemplified in Section 3, the basic procedure of theq-HATM is presented in Section 4, and its algorithm is illustrated for the considered model in Section 5.The existence and uniqueness for the archived results and error analysis are respectively presented in Sections 6 and 7.Moreover, with the aid of behaviour captured for the obtained result, the corresponding comments and conclusion are respectively exemplified in Sections 8 and 9.

2 Preliminaries

The basic definitions are presented in this segment for the FC and Laplace transform.Specifically, we recall the notions related to Caputo-Fabrizio fractional operator [33,73].

Definition 1.The CF fractional derivative forf∈H1(a,b)is presented as [33]

whereM(α)is a normalization function and admitsM(0)=M(1)=1.Further, if1(a,b)then we have

Definition 2.The CF fractional derivative forf∈H1(a,b)is presented as [73]

Note:According to [73], the following must hold

Definition 3.TheLTfor a CF derivativeis presented as [33] below

3 Mathematical Model

On the basis of generic solid tumour growth and assuming it is in avascular stage, the mathematical model has been proposed [74].In this stage, most of the tumours are asymptomatic and further there is a possibility of cells to migrate and escape to the lymph nodes.The considered system exemplifies the interfaces of the surrounding tissue with the tumour and it can be extended to incorporate tumour and the vasculature.Here, the projected system of equations illustrates the interactions of the matrix-degrading enzymes (MDE, signifies byE), extra cellular matrix(ECM, symbolised byC) and tumour cells (denoted byT).With respect to ECM, most of the macromolecules are essential for cell motility, spreading and adhesion.Further, the ECM associated with many macromolecules, for instances collagen, laminin and fibronectin.During the various stages of metastasis, invasion and turn our growth, MDEs play a vital role.The ECM locally degrades by MDEs which are produced by tumour cells.Further, the method wherein they interact with tumour cells, growth factors and inhibitors are highly intricate.The tumour cells in the considered system as haptotaxis and in order to integrate this concept in the model, the hypotactic flux is considered as [75,76]

Jhapto=χT∇C,

whereχ ＞0 denotes haptotactic coefficient and which is constant.The random motion is another contribution to tumour cell motility and it helps to study ECM in isolation.Moreover, the flux is defined for the tumour cells with exemplified random motility is

Jrandom=-D(C,E)∇T,

where ∇Tis the chemokinetic response,D(C,E)is the function of either the ECM or MDE concentration, or constant.

For the tumour cell density(T), the conservation equation is presented as

and for the cell proliferation absence, the equation describing tumour cell motion is defined as

For the notation, the random motility coefficient of tumour cell is considered asD(C,M)=DTand which is constant.Therefore, the degradation process is exemplified by the subsequent equation with positive constantδ

Active MDEs are formed byT, experience some form of decay and diffuse throughout the tissue.The equation modelling the evolution of MDE concentration is presented with MDE diffusion coefficientDEas

whereg=μTandh=λE,handgare the functions respectively describing the MDE decay and the production of active MDEs.Moreover, in the surrounding tissues there is a linear relationship between the level of active MDEs and the density of tumour cells.

From the above description, the system is presented as [18,54,57]:

Here, with appropriate initial conditions, Eq.(10) is assumed to satisfy on a region of tissue or domainΩ.Moreover, the model is nurtured so that the MDEs and tumour cells remain inside the domain of tissue within deliberation and hence no-flux boundary conditions are executed on∂Ω.The terms contained in the above system with ECM density(C0), tumour cell density(T0)and MDE concentration(E0)by setting

wherelsignifies scale length andτis the time.Then we have a scaled system of equations by dropping the tildes for notational convenience [18,74,77]

wheredT=DT/D,γ=χC0/D,η=τE0δ,dE=DE/D,α=τμT0/E0andβ=τλ.

The projected model can be protracted to integrate interactions between blood vessels and the tumour cells [74].

Now, we modify the time derivative by the CF derivative in Eq.(12) and given by

whereαis fractional order.The associated initial conditions are

4 Fundamental Idea of the Considered Scheme

In this section, we hired the differential equation to present the basic procedure of the projected scheme with initial conditions

and

We obtained by applying LT on Eq.(15)

Forφ(x,t;q),Nis contracted as follows:

whereLis signifyingLT.Forq=0 andthe following conditions satisfies

By using Taylor theorem we get

where

For the proper chaise ofv0(x,t),nand ħ the series (13) converges atq=Then

After differentiating Eq.(19)m-times withqand multiplying byand substitutingq=0,one can get

where the vectors are defined as

Eq.(24) reduces after employing inverseLTto

where

and

Here,Hmis homotopy polynomial and presented as

By the help of Eqs.(26) and (27), we found

By the help ofq-HATM, the series solution is

5 Implementation of the q-Homotopy Analysis Transform Method

Consider the system of equation cited in Eq.(13) describing the tumour invasion and metastasis in CF fractional derivative

Applying Laplace transform on Eq.(32) and then with the help of Eq.(13), we get

The non-linear operatorNdefined as

Them-thorder deformation equation by the projected scheme atH(x,t)=1 is given by

where

On employing inverse LT on Eq.(35), it simplifies to

By usingT0(x,t),C0(x,t)andE0(x,t)and then solving the forgoing equations, we can obtain the terms of

6 Existence and Uniqueness of Solutions

In this section, the existence and uniqueness are illustrated for the considered system with the assist of fixed-point theory.We consider the Eq.(32) as follows:

Now, using Eq.(32) and results derived in [53], we obtained

Then we have from [73] as follows:

Theorem 1.The kernelG1admits the Lipschitz condition and contraction if 0 ≤＜1 satisfies.

Proof.Let us consider the two functionsuandu1to prove the theorem, then

Eq.(43) provides the Lipschitz condition forG1.Similarly, we can see that if 0 ≤dTδ2-γ(λ1δ+λ2)＜1, then it implies the contraction.Similarly, we can prove

By the assist of the above equations, Eq.(41) simplifies to

Then we get the recursive form as follows:

The associated initial conditions are

Now, between the terms, the successive difference is defined as

Notice that

Then we have

Application of the triangular inequality, Eq.(50) reduces to

The Lipschitz condition satisfied by the kernelT1, so, we have

Similarly, we have

By the help of above result, we state the following theorem:

Theorem 2.If we have specifict0, then the solution for Eq.(32) will exist and unique.Further,we have

Proof.LetT(x,t),C(x,t)andE(x,t)be the bounded functions admitting the Lipschitz condition.Then, we get by Eqs.(52) and (53)

Therefore, for the obtained solutions, continuity and existence are verified.Now, to prove the Eq.(54) is a solution for Eq.(32), we consider

Let us consider

This process gives

Similarly, att0we can obtain

Asn→∞and from Eq.(57), ‖K1n(x,t)‖→0.Similarly, we can verify for ‖K2n(x,t)‖ and‖K3n(x,t)‖.

Next, for the solution of the projected model, we prove the uniqueness.SupposeT*(x,t),C*(x,t)andE*(x,t)be the set of other solutions, then

Now, employing the norm on the above equation we get

On simplification

From the above condition, it is clear thatT(x,t)=T*(x,t), if

Hence, Eq.(61) proves our required result.

7 Error Analysis of the q-Homotopy Analysis Transform Method

Theorem 3.Let(B[0,T],‖·‖)be a Banach space and supposevn(x,t)andv(x,t)define in the that, then the solution defined in Eq.(31) converges to the solution of Eq.(15), if 0＜λ1＜1.

Proof:Let {Sn} be a sequence of partial sum of Eq.(31).Then, we need to prove {Sn} is Cauchy sequence in(B[0,T],‖·‖).Now, consider

For everyn,m∈N(m≤n), now we have

But 0＜λ1＜1, therefore limn,m→∞‖Sn-Sm‖=0.Hence, {Sn}is the Cauchy sequence.

Theorem 4.The maximum absolute error for the series solution of the Eq.(15) defined in Eq.(31) is determined as

Proof:By using Eq.(62), we get

But 0＜λ1＜0 ⇒1-λ1n-m ＜1.Hence, we have

This ends the proof.

Table 1: Description of parameters presented in the projected system [74]

8 Results and Discussion

Here, we demonstrate the future scheme is efficient and reliable and evaluate the approximate results for the system of partial differential equations representing a model of tumour invasion and metastasis.In the present study, we find the fourth-order solution to present the nature of the system.In Tab.1, we present the specific values of the parameters cited in Fig.1 captures the behaviour ofq-HATM solution for tumour cells(T), extra cellular matrix(C), and matrix degrading enzymes(E)in 3D plots by using the Tab.1 and the combined surface for the three components at the initial stage (i.e.,t=0) is cited in Fig.2.By generalizing the system with a newly nurtured fractional operator, it aids us to capture more interesting consequences associated with singular kernel.In the present work, we demonstrated the nature ofq-HATM results for districtαboth in the change ofxandt, and which are presented in Fig.3.From these curves,we can observe that, as varying in both time and space with fractional order, the obtained results show noticeable vicissitudes in the behaviour.Specifically, extra cellular matrix and matrix degrading enzymes show stimulating behaviour for the changeα.

Figure 1: Surfaces of q-HATM solution for (a) tumour cells (T), (b) extra cellular matrix(C),(c) matrix degrading enzymes (E) at n=1,α=1 and ħ=-1 and using Tab.1

Figure 2: Surface of q-HATM solution for Eq.(32) at n=1,α=1,ħ=-1 and using Tab.1

The behaviors have been captured for different fractional Brownian motions and standard motion(α=1)with the change in ħ.In Fig.4, we drowned the ħ-curves for the obtained solutions forT(x,t),C(x,t)andE(x,t)with the appropriate value of ħ.The ħ-curves aid to adjust and control the convergence province of the achieved results.Fig.5 presents the 2D plots of an analytic-approximate solution for Eq.(32) at a distinct time.By the plots we can see that, the tumour cells and matrix degrading enzymes are also increases while time increases, but the extra cellular matrix decreases.Moreover, these types of investigation can open the door for analyses the stimulating models exemplifying deadly disease by incorporating diffusion co-efficient.

Figure 4: ħ-curves drown for q-HATM solution of (a) T(x,t),(b) C(x,t),(c) E(x,t) for distinct α at t=0.01,x=0.1,n=1 and using Tab.1

Figure 5: Response of obtained solution for the considered model with varying in x at (a) initial time (t=0), (b) t=0.1, (c) t=0.5 and (d) t=1 with n=1,α=1,ħ=-1 and Tab.1

9 Conclusion

In the present study, we analyzed and capture the behaviour of the nonlinear fractional model of tumour invasion and metastasis by using the fractional operator and efficient analytical technique.The existences and uniqueness are demonstrated with the assist of a fixed point hypothesis.The plots captured in the present investigation display the stimulating behaviour and these can help scholars for some essential and interesting consequence of the hired system.The present study shows, the phenomena conspicuously be contingent on the time history and the time instant and, these can be proficiently studied using fundamental perceptions of FC and newly proposed fractional operator.The investigations of these types of models can provide new notions to analyze more real-world problems and it opens the door for employing an efficient method to study complex phenomena associated with science and technology.

Funding Statement:The authors received no specific funding for this study.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.