key: cord-0723954-nd7r03cy
authors: Kumar, Sunil; Chauhan, R. P.; Aly, Ayman A.; Momani, Shaher; Hadid, Samir
title: A study on fractional HBV model through singular and non-singular derivatives
date: 2022-03-02
journal: Eur Phys J Spec Top
DOI: 10.1140/epjs/s11734-022-00460-6
sha: f68ca86febad1dd51ec3454b41a9fd56aef64387
doc_id: 723954
cord_uid: nd7r03cy

The current study’s aim is to evaluate the dynamics of a Hepatitis B virus (HBV) model with the class of asymptomatic carriers using two different numerical algorithms and various values of the fractional-order parameter. We considered the model with two different fractional-order derivatives, namely the Caputo derivative and Atangana–Baleanu derivative in the Caputo sense (ABC). The considered derivatives are the most widely used fractional operators in modeling. We present some mathematical analysis of the fractional ABC model. The fixed-point theory is used to determine the existence and uniqueness of the solutions to the considered fractional model. For numerical results, we show a generalized Adams–Bashforth–Moulton (ABM) method for Caputo derivative and an Adams type predictor-corrector (PC) algorithm for Atangana–Baleanu derivatives. Finally, the models are numerically solved using computational techniques and obtained results graphically illustrated with a wide range of fractional-order values. We compare the numerical results for Caputo and ABC derivatives graphically. In addition, a new variable-order fractional network of the HBV model is proposed. Considering the fact that most communities interact with each other, and the rate of disease spread is affected by this factor, the proposed network can provide more accurate insight for the modeling of the disease.

Hepatitis B virus is an infectious disease that infects a large number of people globally. The virus causes liver damage as well as other diseases that are harmful to human health. Hepatitis B epidemiology indicates two transmission routes in humans: horizontal and vertical. Horizontal transmission occurs when people share blades, needles, unprotected intercourse, and so on, while vertical transmission occurs when an infected mother passes the virus on to her kid. In 2015, many HBV cases of liver infection and death were registered and still, several countries around the world, including China, are dealing with this problem with many Guest editors: Oscar Castillo, Hadi Jahanshahi, Amin Yousefpour. a e-mail: skumar.rs.apm@itbhu.ac.in (corresponding author) b e-mail: aymanaly@tu.edu.sa cases and deaths. Acute HBV clinical information is age-specific. HBV infection is normally asymptomatic in newborns, young children (under the age of 10) and immunocompromised adults. Adults and older children with 30-50% infection are symptomatic. Sometimes the hepatitis B-infected person who does not show symptoms for more than 6 months may be prone to liver damage or even death, although some research shows that an HBV-infected individual without having symptoms has the potential to infect another individual and could lose their life [1] [2] [3] . Many mathematical models have been developed to study HBV and its characteristics by researchers from all over the world. For instance, HBV with control strategies and analysis with different age structures with population of New Zealand presented in Refs. [4, 5] , control study of hepatitis B via vaccination [6] , vaccination effects on HBV transmission presented in Ref. [7] , for HBV disease a multi-group model presented in Ref. [8] , HBV model with cost-effective analysis presented in Ref. [9] , significance of vaccination on HBV transmission analyzed via a mathematical model in Ref. [10] , HBV model with diffusion and delay studied in Ref. [11] and model with time delay studied in Ref. [12] , the HBV transmission dynamics and its control via a mathematical model studied in Ref. [13] , migration effect on HBV transmission presented in Ref. [14] , and an HBV model with three controls discussed in Ref. [15] . The above-described HBV models are in the frame of integer-order derivative. As the models with ordinary order do not provide the memory and the genetic properties involved in the fractional-order models. Therefore, in recent years, researchers have given more attention to fractional-order derivatives for mathematical modeling in mathematical biology because the fractional derivatives reflect memory effects and genetic properties in models.

Fractional Calculus (FC) is developing as a robust tool in applied mathematics to solve real-world issues. The fractional-order (FO) models are very suitable for understanding the memory effect and inherited properties than integer-order models. In addition, the integer-order derivative does not examine the dynamics between two different points in the real-world description. Several non-local FO differentiation concepts have been proposed in the existing article to address the limitations of classical local differentiation. The fractional derivatives such as Caputo [16] , Caputo-Fabrizio [17] and Atangana-Baleanu (AB) [18] are receiving much attention from the researchers. Because of its Mittag-Leffler (ML) memory and crossover properties for mean-square displacement, the ABC derivative offers an excellent description. In recent years, the Atangana-Baleanu derivative has been applied to the epidemiological model extensively [19] [20] [21] [22] [23] [24] [25] . Several other recent research on application of fractional derivatives to the non-linear models from different engineering and science fields have been studied in Refs. [26] [27] [28] [29] [30] [31] [32] . A novel coronavirus model is investigated with ABC fractional derivative in Ref. [26] . The transmission of Hookworm infection is studied via the Caputo operator in Ref. [28] . In Ref. [33] , COVID-19 model is investigated through Caputo fractional derivative. A pine wilt disease model is examined under the Caputo fractional operator in Ref. [34] . In Ref. [35] , the authors studied a combined co-infection model of HIV-COVID-19 through ABC fractional derivative. In Ref. [36] , wastewater treatment model is examined with Caputo fractional derivative. Some fractional models regarding HBV can be found in the existing literature. A FO HBV model with optimal control is investigated in Ref. [37] , and optimal control on HBV model with vaccination, treatment is discussed in Ref. [38] . Using the ABC operator, the HBV model with hospitalized class is studied in Ref. [39] and HBV model with treatment is studied in Ref. [40] . HBV model with hospitalization using Caputo derivative is studied in Ref. [41] . A FO HBV model with antibody immune response is discussed in Ref. [42] .

Recently, an HBV model with asymptomatic carriers was studied using Caputo derivative in Ref. [43] .

Thus, inspired by recent studies, we consider a FO HBV model [43] with Caputo and ABC derivatives.

We study the dynamics of the proposed model using two different numerical techniques. The proposed work is divided as follows: some basic preliminaries of FC are presented in Sect. 2. The FO HBV model in the frame of Caputo and ABC derivatives along with some mathematical results are presented in Sect. 3. The existence and uniqueness of solutions of FO model are presented in Sect. 4. The numerical schemes and obtained numerical results are shown in Sect. 5. A variable-order fractional network of the HBV model is proposed in Sect. 6, also the numerical simulation is presented for the model. Lastly, the outcomes are discussed in Sect. 7.

Here, we present a quick glance of some basic concept of FC [16, 18] .

C 0 D p t χ(t) = 1 Γ(1 − θ) t 0 (t − υ) −p χ (υ)dυ, 0 < p < 1. (2.1)

where J(p) is the normalization function satisfying J(0) = J(1) = 1, and E p is ML function

The corresponding fractional integral of ABC derivative is given by 

(2.5)

Here, we have considered an HBV model with class of asymptomatic carriers [43] . In the model, total popula-tion represented by N (t) divided into six classes, susceptible individuals (S(t)), exposed population (E(t)), acute infected population (A(t)), asymptomatic carrier (A c (t)), chronic infected individuals (C(t)) and the recovered population (R(t)), so

Initially, the HBV model in frame of Caputo derivative is given as follows:

Further, the HBV model in frame of ABC derivative is given as follows:

where the initial conditions S(0) ≥ 0, E(0) ≥ 0, A(0) ≥ 0, A c (0) ≥ 0, C(0) ≥ 0 and R(0) ≥ 0, and the parameter Λ 1 is birth rate of the susceptible individuals. The effective contact rate and natural fatality rate are described by the parameters β 1 and μ 1 , respectively. The exposed individuals are infected at the speed ψ 1 (1 − θ 1 ), where the number of individual move to class A at the rate of ψ 1 θ 1 while a portion of ψ 1 (1 − θ 1 ) enters A c and become asymptomatically infected. Individuals in acute and asymptomatic classes become carriers at rates of η 1 and τ 1 , respectively. The recovery rate of acute, asymptomatic and the carriers individuals are denoted by κ 1 , ν 1 and σ 1 , respectively. The death rate because of disease at acute and chronic classes are represented by γ 1 and ε 1 respectively. The coefficient of asymptomatic and carriers individuals are shown by φ 1 (the infectiousness of asymptomatic relative to acute infections) and 1 (the infectiousness of carriers relative to acute infections), respectively.

is positively invariant for the ABC fractional model (3.2) .

Proof The dynamics of the total population is obtained by summing all equations of the model (3.2)

Hence,

employing the Laplace transform on Eq. (3.5), we have

after simplification we get

and E p,q is two-parameter ML function and describe as

Due to asymptotic nature of ML function, we obtain

Thus, the entire solutions of the proposed model for the initial conditions belonging to Θ stay in Θ for every t > 0. Hence, Θ is a positively invariant region and will attract all the solutions in R + 6 .

For the model (3.2), there exist two equilibrium points, disease free equilibrium point (DFE)ξ 0 and the endemic equilibrium (EE)ξ 1 . The DFE point is given as

The basic reproduction number R 0 of the proposed model (3.2) is achieved using next-generation matrix technique [44] . For the reproduction number R 0 , we use [44] and obtain the matrices F and V as follows:

Therefore, using the ρ(FV −1 ), we obtained the basic reproduction number for the model (3.2)

The EE point of the model (3.2) is given byξ

using these in first equation of model (3.2), we get 

Here, we will examine the existence of unique solution of fractional-order model (3.2) . Assume that H (J) is a Banach space containing real-valued continuous functions defined on the interval J = [0, T ] with sup norm

Now, using the definition (2.4) on above equation, we get

A c (t)), Ξ 5 (t, C(t)) and Ξ 6 (t, R(t)) are said to satisfy the Lipschitz condition if and only if S(t), E(t), A(t), A c (t), C(t) and R(t) have an upper bound. Let S(t) and S 1 (t) be two functions, then we get

(4.5)

In a similar manner, one can obtain

(4.6)

Thus, the Lipschitz condition is satisfied for all the six functions Ξ 1 , Ξ 2 , Ξ 3 , Ξ 4 , Ξ 5 and Ξ 6 where 1 , 2 , 3 , 4 , 5 and 6 are the corresponding Lipschitz constants. Using the recursive principle, the Eq. (4.2) becomes

where the initial conditions are defined as follows:

Taking the difference of successive terms, we obtained the following expression:

In addition, using Eqs. (4.5)-(4.6) and taking into consideration that Ψ S,n−1 (t) = S n−1 (t)−S n−2 (t), Ψ E,n−1 (t)

we have the following: 

Proof It is evident that S(t), E(t), A(t), A c (t), C(t) and R(t) are the bounded functions. Furthermore, we can identify from Eqs. (4.5)-(4.6), the expressions Ξ 1 , Ξ 2 , Ξ 3 , Ξ 4 , Ξ 5 and Ξ 6 satisfy the Lipschitz condition. Thus, utilizing Eq. (4.9) along with a recursive formula, we arrive at

(4.11) Therefore, the above specified sequences exist and satisfy Ψ S,n (t) → 0, Ψ E,n (t) → 0, Ψ A,n (t) → 0, Ψ Ac,n (t) → 0, Ψ C,n (t) → 0, Ψ R,n (t) → 0, as n → ∞. Moreover, from Eq. (4.11) and using the triangle inequality for any m, we get

where

J(p)Γ(p) i < 1 by hypothesis. Therefore S n , E n , A n , A cn , C n , R n are described as Cauchy sequences in the Banach space H (J). Therefore, they are uniformly convergent as mentioned in [45] . Applying the limit theorem on Eq. (4.7) when n → ∞ indicates that the limit of these sequences is the unique solution of the model (3.2) . Finally, the existence of a unique solution for model (3.2) has been achieved under the condition (4.10).

are constantly looking for a suitable numerical approach to solve them. Now, we will present the solution of the mentioned fractional-order HBV model.

Here, we present ABM method (see for more details on the ABM scheme and its convergence) to the approximate solution of the proposed Caputo HBV model. Following from Refs. [46] [47] [48] , setting h = T N , t n = nh, n = 0, 1, 2, . . . , N, the model solution can be written as follows:

j,n+1G1(tj , Sj, Ej, Aj, Ac j , Cj, Rj), 

j,n+1G2(tj , Sj, Ej, Aj, Ac j , Cj, Rj),

j,n+1G3(tj , Sj, Ej, Aj, Ac j , Cj, Rj),

j,n+1G4(tj , Sj, Ej, Aj, Ac j , Cj, Rj),

j,n+1G5(tj , Sj, Ej, Aj, Ac j , Cj, Rj),

j,n+1G6(tj , Sj, Ej, Aj, Ac j , Cj, Rj).

and

Here, we will show a Adams type PC numerical technique [49] with Atangana-Baleanu fractional integral operator for numerical results of the proposed HBV model (3.2) . For the sake of simplicity, the model (3.2) is given as Now, we consider the system (5.3) with Atangana-Balenau (AB) integral and get the following:

S, E, A, Ac, C, R)

S, E, A, Ac, C, R) 

S, E, A, Ac, C, R)

S, E, A, Ac, C, R) 

S, E, A, Ac, C, R) We need now a numerical approximation technique for AB fractional integral. For this, we utilize the Adams type PC numerical approach for AB fractional integral. Therefore, we have mula of integral form of derivative is provided as follows: The predictor value χ P h (t n+1 ) is defined by following:

where,

Operating above Adams type PC technique on the HBV model (3.2), we get the following iterative formula:

(e) (f) Fig. 6 Numerical simulation for the HBV model (3.2) for p = 0.8, when R0 < 1

(e) (f) Fig. 7 Comparison of numerical result for the HBV model for p = 1, when R0 > 1

(e) (f) Fig. 8 Comparison of numerical result for the HBV model for p = 0.9, when R0 > 1

are given as 

δj,n+1G4(tj, Sj, Ej, Aj, Ac j , Cj, Rj),

G5(tn, Sn, En, An, Ac n , Cn, Rn)

δj,n+1G5(tj, Sj, Ej, Aj, Ac j , Cj, Rj),

G6(tn, Sn, En, An, Ac n , Cn, Rn)

δj,n+1G6(tj, Sj, Ej, Aj, Ac j , Cj, Rj).

(5.9)

In the present subsection, we have discussed the obtained numerical results for the Caputo model (3.1) and ABC model (3.2) with unit of time per year through the above suggested numerical techniques. The values of model parameter used in the numerical simulations are Λ 1 = 2, μ 1 = 1/67.7, β 1 = 0.042, φ 1 = 0.002, 1 = 0.002, ψ 1 = 0.004, θ 1 = 0.6, γ 1 = 0.001, η 1 = 0.02, κ 1 = 0.02, τ 1 = 0.02, ν 1 = 0.1, ε 1 = 0.003, σ 1 = 0.2. Using the above given parameter values except β 1 = 0.042, we show the numerical simulation for the model (3.2) in Figs. 1, 2, and 3, when R 0 > 1. From the presented graphics, we observe that the susceptible population decrease while the exposed individuals increase, the infected classes, acute, asymptomatic and chronic population increase, the increase is slower for the smaller values of FO parameter. Making the parameter β 1 = 0.00042 too small, we get R 0 < 1 and the graphical results are shown for this in Figs. 4, 5, and 6. The graphical solutions display that by decreasing the value of the non-integer-order parameter p, the number of individuals in the infected compartment decreases. In Figs. 7, 8, and 9 , the comparison of numerical result for Caputo and ABC operators are presented graphically. We observe that ABC provide better results for the smaller fractional order. In Fig. 10 , taking distinct values of parameter η 1 with p = 0.95, we show the behavior of the infected compartments, acute, asymptomatic and chronic individuals. In addition, we provided the CPU output time in Tables 1 and 2. Eur. Phys. J. Spec. Top. In this section, we propose a new mathematical model that considers the interaction of different communities in modeling the disease. Since the relationships and travels between two countries will definitely affect the rate of diseases in both of them, by embedding the interaction between different communities in the modeling of the disease, more realistic results will be obtained.

The mathematical model of the disease is considered as neurons affected by each other in a network. To construct a network of HBV model, using the nearest neighbor method some models are joined. We suppose that the interactions between neurons are through the susceptible as well as asymptomatic carrier components, i.e., the coupling is made on the variable S and A c . Consequently, the governing equation of the proposed network is given by

where p(t) stands for the time-varying fractional-order derivative; also, q is the number of nearest neighbors on each side, and l denotes the coupling strength. Figure 11 demonstrates the numerical results of the third neuron of the proposed network. In this case, the number of communities that are in the network is considered to be N S = 4. In addition, q = 2 and l = 1. The fractional-order derivative is considered to be a function of time as p(t) = 0.95 + 0.01 cos( t 20 ). By comparing the results of the HBV model in the network and the case that there is a sole community, it can be found out that the neighborhood communities play important roles in the spread of the disease.

In the present work, we explored the dynamics of the HBV model with the asymptomatic carrier-class using two different fractional operators. For the arbitrary order fractional HBV model, we considered the Caputo and ABC fractional derivatives and presented some mathematical analysis for the FO ABC model, such as invariant region and basic reproduction number. Furthermore, with fixed-point results, the existence and uniqueness of system solutions to the arbitrary order HBV model were investigated. The numerical solution was obtained via the presented numerical schemes and the numerical outcomes were explored graphically for different choices of FO parameter. The importance of arbitrarily chosen fractional order can be observed from our presented results. Thus, from the analysis given in this work, we can conclude that the models formulated in FO demonstrate some interesting results and several characteristics that cannot be recognized in integer order. Finally, a network of the HBV model was introduced to consider the effects of different societies on each other. In the proposed model, it is supposed that the value of the fractional derivative is changing with respect to time, and the populations are jointed through the nearest neighbor method. In a future study, we will investigate the control of the proposed variableorder network through an artificial intelligence-based approach. This model can also be studied with other fractional derivatives, such as generalized Caputo and conformable derivatives.

Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications

Introduction to Functional Analysis

Acknowledgements The research was supported by the Taif University Researchers Supporting Project number (TURSP-2020/77), Taif University, Taif, Saudi Arabia.