
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39300092
70788
10.1038/s41598-024-70788-7
Article
Predictive modeling of hepatitis B viral dynamics: a caputo derivative-based approach using artificial neural networks
Turab Ali 1
Shafqat Ramsha ramshawarriach@gmail.com

2
Muhammad Shah 3
Shuaib Mohammad 4
Khan Mohammad Faisal 5
Kamal Mustafa 5
1 https://ror.org/01y0j0j86 grid.440588.5 0000 0001 0307 1240 School of Software, Northwestern Polytechnical University, 127 West Youyi Road, Beilin District, Xi’an 710072 China
2 https://ror.org/051jrjw38 grid.440564.7 0000 0001 0415 4232 Department of Mathematics and Statistics, The University of Lahore, Sargodha, 40100 Pakistan
3 https://ror.org/02f81g417 grid.56302.32 0000 0004 1773 5396 Department of Mathematics, College of Science, King Saud University, P.O. Box 2455, Riyadh, 11451 Saudi Arabia
4 https://ror.org/00et6q107 grid.449005.c 0000 0004 1756 737X Department of Mathematics, Lovely Professional University, Phagwara, Punjab India
5 https://ror.org/05ndh7v49 grid.449598.d 0000 0004 4659 9645 Department of Basic, Sciences College of science and theoretical studies, Saudi Electronic University, Riyadh, 11673 Saudi Arabia
19 9 2024
19 9 2024
2024
14 2185314 5 2024
21 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
A fractional model for the kinetics of hepatitis B transmission was developed. The hepatitis B virus significantly affects the world’s economic and health systems. Acute and chronic carrier phases play a crucial part in the spread of the HBV infection. The Hepatitis B infection can be spread by chronic carriers even though they show no symptoms. In this article, we looked into the Hepatitis B virus’s various stages of infection-related transmission and built a nonlinear epidemic. Then, a fractional hepatitis B virus model using a Caputo derivative and vaccine effects is created. First, we determined the proposed model’s essential reproductive value and equilibria. With the aid of Fixed Point Theory, a qualitative analysis of the problem’s approximative root has been produced. The Adams-Bashforth predictor-corrector scheme is used to aid in the iterative approximate technique’s evaluation of the fractional system under consideration that has the Caputo derivative. In the final section, a graphical representation compares various noninteger orders and displays the discovered scheme findings. In this study, we’ve utilized Artificial Neural Network (ANN) techniques to partition the dataset into three categories: training, testing, and validation. Our analysis delves deep into each category, comprehensively examining the dataset’s characteristics and behaviors within these divisions. The study comprehensively analyzes the fractional HBV transmission model, incorporating both mathematical and computational approaches. The findings contribute to a better understanding of the dynamics of HBV infection and can inform the development of effective public health interventions.

Keywords

Mittag-Leffler
Vaccination
Fractional derivative
Adams-Bashforth
Mild solution
Numerical simulations
Artificial neural networking
Subject terms

Engineering
Mathematics and computing
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The hepatitis B virus (HBV) is a serious viral disease affecting millions worldwide. The virus also causes other illnesses related to human health, such as liver damage. Hepatitis B epidemiology proposes two human transmission routes, such as horizontal and vertical. Although vertical transmission occurs when a woman who has the virus passes it on to her kid, horizontal transmission involves sharing objects like razors, needles, and sexual intercourse. 2015 saw a significant increase in the number of HBV infections associated with liver infections and fatalities, and several nations throughout the world, including China, are still dealing with this problem.

More than 350 million people worldwide have chronic HBV bacterial infections1; these disorders have a 26–41% risk of developing fibrosis, hepatocellular carcinoma, and severe liver diseases2. For this, HBV infection is a major global health concern. According to this information, a person infected with the HBV infection has a plasma weight of 2 × 1011 to 3 × 10123. Thus, a liver weighs 1.5 kg, nearly the same as the human body’s total number of liver cells. Because of its size, we can only employ the routine rate of occurrence in turn of the incidences for adequate pressure. HBV is directly produced after two or three days after therapy4. It is included in the article for no more than 7 to 13 months5.

It is now widely acknowledged that epidemiological formulations create conflicting effects, appropriately in the safety and management of viral disorders6–8. Complex mathematical models with various premises and characteristics are created by examining the literature to study various illnesses. Generally, the intensity of viruses in people’s viral, bacterial infections is affected by probabilistically disturbed components (such as moisture, rainfall, and heat dissipation on someone else). Moreover, we can probabilistically examine the idealized models using the same methodologies in the formal constructs. It is acknowledged that such actual models can also include the impact of changing environmental conditions. The oscillation of factors provided by the outside atmosphere or the sounds created by systems in the real world to make up factors affections9. The stochastic models with a probabilistic system will be close to precise results in the comparative and provide an increased layer of options. For detailed testing and a more extensive discussion of the literature on probabilistic representations of modeling interrupted by Zig-Zag phenomena, see10. Compared to predictable models for different processes, the outputs of probabilistic solutions are closer to those of big issues.

Because of the various options of orders of derivatives, fractional differential and integral equations that depict real-world occurrences are essential in studying nonlinear dynamics. The current calculus of fractional differentiation and integral has been shown to have numerous applications for various issues in a more exact and accurate style11. The conclusions drawn through noninteger order products are viewed to be highly accurate and have an extra bit of dynamics, especially with traditional models. This allows for exploring epidemiology and different ecological, social, physiological, and financial realities using worldwide fractional-order derivatives. Academics and scientists are more interested in studying the new notion of differentiation and integrals than they are in studying traditional derivatives. Specialists and mathematicians are attempting to focus on numerous sectors of noninteger order derivatives for outcomes due to the more actual results of present calculus. Various real-world issues, including modeling of viral infections and the fundamental Lotka-Volterra system, have been addressed using variable orders integral equations and derivative equations12,13. The modeling of noninteger order differential equations is covered in the research using a variety of approximation and exact techniques, including Predictor-Corrector, Adams-Bashforth, and other well-known approaches14,15. Numerous academics have developed mathematical models of COVID-19’s viscous flow; for instance, see16. Additionally, a wealth of data has been published on using current calculus for tuberculosis models, and numerous researchers have explored biology and physics issues due to FDEs17.

The study of infectious disease models has gained significant attention in recent years, with researchers exploring various mathematical frameworks to understand disease dynamics and inform public health interventions. In particular, the application of fractional-order derivatives has emerged as a promising approach to model complex epidemiological phenomena. Several studies have investigated fractional-order models for diseases such as toxoplasmosis18,19, HIV/AIDS20, tuberculosis21, pine wilt disease22,23, typhoid fever24, and malaria25. These models have demonstrated the ability to capture the non-integer order dynamics inherent in many biological systems. Additionally, researchers have explored the use of stochastic and discrete-time models to study disease propagation, accounting for the inherent randomness and heterogeneity in real-world scenarios26,27. Furthermore, the incorporation of advanced machine learning techniques, such as ensemble classifiers28, has shown potential in improving disease classification and prediction. The growing body of literature in this field highlights the multifaceted and interdisciplinary nature of epidemiological modeling, providing a strong foundation for the current study29. Fractional-order models have emerged as a promising approach for studying the complex dynamics of infectious diseases. Recent research has explored the application of these models to analyze diseases such as human papillomavirus30 and disease dynamics in prey-predator systems31. The incorporation of fractional calculus and fractal theory in these models has provided valuable insights into the underlying mechanisms and transmission patterns31. Continued advancement in this field is crucial for developing innovative techniques to better understand and manage infectious disease challenges.

Using the findings of the prior studies, we create a viral illness method for HBV utilizing the idea of the fractional operator by using a mild solution with the Vaccination class and Asymptomatic class in this article. To investigate the epidemic more thoroughly, we build the characteristics of HBV using Vaccinated and Asymptomatic class carriers, which have multiple elements for examining HBV dynamical behavior in the research. The authors of17 have investigated the HBV virus using a brand-new mathematical framework with an Asymptomatic class. Individuals with HBV infection may remain symptom-free for 31 to 181 days and spread the illness to others, which can result in death in human civilizations32,33. Our assessment attempts to develop a brand-new Hepatitis (HBV) model that includes vaccination and asymptomatic carriers based on the research on HBV modeling34. Initially, using the numerals of the Atangana-Baleanu fractional operator, a unique constructive derivative, the behavior of the dynamic system is produced in the context of integer or natural order derivatives and spreads to noninteger order derivatives. The main goal of this research project is to easily compare human health issues related to vaccination that cause a delay in the effective management of the HBV epidemic through mathematical models. We will finish developing effective vaccines to prevent illnesses from emerging in the actual world.

The document’s arrangement continues as follows: A brief overview of the spread of HBV and the contributions of well-known experts in this field are provided in the section under “Introduction.” Section “HBV mathematical model” presents the construction of the model mathematically as well as its fundamental features.The domestic and international stability of the stochastic model are presented in Section “Stability Analysis.” We reconstructed the model with fractional derivative in Section “Fractional order HBV model.” In the part under “Preliminaries,” we also include a few fundamental lemmas and definitions. The examination of the recommended fractional model’s solution is presented in the section under “Existence of Solution.” The suggested fractional model’s Hyers-Ulam stability is covered in Section “Hyers-Ulam stability.” We discussed the proposed model’s mathematical setup and visual outcomes in Section “Numerical Solutions.” The result of the work done in the paper is briefly covered in the section under “Conclusion.”

Model for the HBV virus

Vaccinations and time-delay inspections are price strategies for reducing the HBV viral sickness from the perspective of actual medical examinations. The World Health Organization correctly advised that effective HBV immunization must be used according to step-by-step criteria in areas with high infection rates. Several scientists have created stochastic HBV epidemic simulations using vaccination strategies as a result5. The model for the nonlinear occurrence of HBV transmission is given. The cases are split into five groups based on the features of HBV: the susceptible population S(t), the vaccination class V(t), the acute infectious class A(t), the chronic infection C(t), and the immunized R(t). N(t)=S(t)+V(t+A(t)+R(t)+C(t) is the population density over all time t; see Fig. 1.Fig. 1 Illustration of the transmission dynamics of Hepatitis.

Furthermore, we presumed A1 : The variables and stochastic processes of the simulation are good.

A2 : Because the duration of acute infection is relatively brief if treatment is unsuccessful, the infected individual will revert to the chronic compartment.

A3 : We have excluded the mortality rate associated with HBV if liver histological has reduced due to therapy.

A4 : An individual is considered susceptible forever once they have received a vaccination or been cured after treatment.

A5 : The exterminated species possesses protection.

The situation under consideration is created and enhanced by Din et al.35. Numerous formulas of epidemic diseases with a prevalence rate of the Harmonic Mean type (HMTIR) have been researched in the literature36,37. Therefore, the paper examines a formulation used in epidemiology to describe the incidence of HBV in a society. We create the HMTIR between good health and acutely affected participants and between acutely infected and vaccinated people to determine which contagious patients go rapidly. The normality test of a data collection is measured as the average of the two numbers. The geometric mean is utilized for average ratios or frequencies by varying the provided data. Despite being considered, a harmonic mean performs less well for some massive quantities than geometric or arithmetic means. Occasionally, it is used when the variance of the parameters is substantial. The probability that the density will achieve habitat loss in a specific amount of time but far faster than other event ratios is known as the harmonic occurrences. Whenever a number determines a person’s value, the harmonic mean is a decent incidence frequency compared to similar units. The Harmonic incidence rate also has a genuine mean of the rates and fractions in various studies. It is susceptible to a single number that is below average. Depending on those above (A1-A5) hypotheses and other associated constraints, the behaviors of HBV are described by the simulation tool:

Several practical issues, including epidemiological dynamical scenarios, demonstrate how dependent systems with global features are crucial to the validity of natural problems. Owing to fractional-order derivatives’ time-dependent kernels, systems with noninteger-order derivatives represent these problems well. The Caputo noninteger order derivative is the more generic of the various fractional operators described in the literature. The significance and benefit of using Caputo fractional-order derivatives is that they use begin assumptions comparable to modern differential calculus, i.e., they do not necessitate noninteger order initial assumptions14.

Motivated by38, our study’s goal is to examine a newly developed System 1 that uses noninteger order differentiation operators with a non-local and singular kernel since it produces more accurate outcomes than naturalistic orders11. We discover the existence of our model under consideration’s root using the fixed point theory notion. An arbitrary-orders operator system in the meaning of the suggested approach is examined:1 cDtϑS(t)=κΛ-2βS(t)A(t)s(t)+A(t)+V(t)ρ-S(t)(μ+η),cDtϑV(t)=Λ(1-κ)+μS(t)-βV(t)A(t)(1-τ1)2V(t)+A(t)-V(t)(ρ+η),cDtϑA(t)=2βS(t)A(t)S(t)+A(t)+2(1-τ1)βV(t)A(t)V(t)+A(t)-(α1+α2+η)A(t),cDtϑC(t)=α1A(t)-(γ1+γ2+η)C(t),cDtϑR(t)=γ1C(t)+α2A(t)-ηR(t).

We will investigate the model underlying initial conditions for ϑ∈(1,2),S(0)=S0,V(0)=V0,A(0)=A0,C(0)=C0,R(0)=R0,S′(0)=S0+n(S0),V′(0)=V0+n(V0),A′(0)=A0+n(A0),C′(0)=C0+n(C0),R′(0)=R0+n(R0).

The governing equations are described in Table 1. The findings show that the medicines are not pure because they only partially cure illnesses. By doing so, those who have received vaccinations can contract illnesses after coming into contact with seriously ill people. It should be remembered that 0<τ1<1(τ1=1) indicates a high vaccine, whereas τ1=0 indicates a vaccine that provides no safety at any time.

In this part, we determine the figure’s viability, stability point, and fundamental biological value for 1 with the help of Table 1:Table 1 Parameters used in the numerical simulations of model 1.

Parameters	Description	
Λ	Death rate among infants	
β	Infectiousness level	
μ	The percentage of nearby people who are ill	
γ1	Percentage of Chronically Affected Patients Recovering	
γ1	There is a significant danger owing to the frequency of people who have HBV infection that is recurrent.	
α1	Rapid group transfers to Excessively afflicted at a higher percentage.	
α2	Individuals affected by a serious illness have an abysmal recovery rate.	
κ	Numbers of young babies born despite having received any form of effective immunization	
ρ	Hepatitis vaccination frequency	

Region of invariant

We view all of the model’s variables as positive. Since there are real Humans in the environment, the state variables are believed to be favorable at time υ=0. In the form of differential equations, the behaviors of the entire human population are as regards:2 DtϑN(t)=κΛ-ηN(t)γ2C(t)≤κΛ-ηN(t).

The system’s orientation equilibrium isN=κΛη

Consider the following biologically viable area Δ:Δ={(S,V,A,C,R)∈R+5,N≤κΛη}.

When equation 2 is resolved, we get3 N=N(0)e-ηt+κΛη(1-e-ηt).

SoN→κΛηυ→∞.

Because of this, the positively invariant domain is Δ, and the model is theoretically and epidemiologically bounded.

Reproduction number R0

To compute the basic reproduction number for our system, we follow as38:

Let x=(A,C) is the infectious compartment, then it give from model 1 that:4 cDtϑA(t)=2βS(t)A(t)S(t)+A(t)+2(1-τ1)βV(t)A(t)V(t)+A(t)-(α1+α2+η)A(t),cDtϑC(t)=α1A(t)-(γ1+γ2+η)C(t).

Now, using the next generation matrix approach, we have the following:5 dxdt=ξ-v,

whereξ=2βS(t)A(t)S(t)+A(t)+2(1-τ1)βV(t)A(t)V(t)+A(t)0,v=(α1+α2+η)A(t)-α1A(t)+(γ1+γ2+η)C(t).

The spectral radius R0 of the matrix is6 R0=2β+2β(1-τ1)(α1+α2+η)

and it is the required basic reproduction number R0 for the system 1.

Preliminaries

Below is a list of the results we can still recall.

Lemma 2.1

(Gronwall lemma)39 Assume μ , y∈H([0,1],R+) and let μ be increasing. If u∈H([0,1],R+) satisfiesu(t)⩽μ(t)+∫0ty(s)u(s)ds,t∈[0,1],

thenu(t)⩽μ(t)exp∫0ty(s)u(s)ds,t∈[0,1].

Definition 2.1

[40] The fractional RL derivative is defined asaDωpχ(ω)=1Γ(n-p+1)(ddω)n+1∫aω(ω-τ)n-pχ(τ)dτ,n⩽p⩽n+1.

Definition 2.2

40] The Caputo fractional derivatives aCDωαχ(ω) of order α∈R+ are defined byaCDωαχ(ω)=aDωα(χ(ω)-∑ȷ=0k-1χ(ȷ)(a)ȷ!(ω-a)ȷ),

in which k=[α]+1.

Definition 2.3

41 The Wright function, ψα is defined byψα(κ)=∑ȷ=0∞(-κ)ȷȷ!Γ(-αȷ+1-α)=1π∑ȷ=1∞(-κ)ȷ(ȷ-1)!Γ(ȷα)sin(ȷπα),υ∈(0,1),κ∈C.

Lemma 2.2

42 Consider {H(ω)}ω∈R as a strongly continuous cosine category in X fulfilling ‖H(ω)‖Lb(X)≤MeΩ|ω|,ω∈R, and A as an infinitesimal generator of {H(ω)}ω∈R. Then for Re(λ)>ω and λ2∈ρ(A),λR(λ2;A)x=∫0∞e-λωH(ω)xdω,R(λ2;A)x=∫0∞e-λωS(ω)xdω.

Theorem 2.3

If X(t) fulfills (1) of equation (4) then Y(t) is given byY(t)=Y0+(Y1+m(Y1))t+1Γ(β)∫0t(t-s)q-1Pq(t-s)Ω(x,Y(x))dx,t∈[0,∞),

holds, thenY(t)=Cq(t)(Y0+m(Y0))+Kq(t)(Y1+m(Y1))+∫0t(t-s)q-1Pq(t-s)Ω(x,Y(x))dx,

∀t∈[0,b], such thatCq(t)=∫0∞Mq(θ)C(tqθ)dθ,Kq(t)=∫0tCq(s)ds,Pq(t)=∫0∞qθMq(θ)C(tqθ)dθ,

where Cq(t) and Kq(t) are continuous with K(0)=I and C(0)=I, |Cq(t)|⩽c,c>1 and |Kq(t)|⩽c,c>1,∀t∈[0,T].

Existence of solution

To consider the possibility that the presented problem has a solution, we first construct the mapping as regards:7 Θ1(S,V,A,C,R,t)=κΛ-2βS(t)A(t)s(t)+A(t)+V(t)ρ-S(t)(μ+η),Θ2(S,V,A,C,R,t)=Λ(1-κ)+μS(t)-βV(t)A(t)(1-τ1)2V(t)+A(t)-V(t)(ρ+η),Θ3(S,V,A,C,R,t)=2βS(t)A(t)S(t)+A(t)+2(1-τ1)βV(t)A(t)V(t)+A(t)-(α1+α2+η)A(t),Θ4(S,V,A,C,R,t)=α1A(t)-(γ1+γ2+η)C(t),Θ5(S,V,A,C,R,t)=γ1C(t)+α2A(t)-ηR(t).

The created problem can be expressed as presented using 7:8 cDtϑZ(t)=Θ(t,Z(t)),t∈[0,ℑ],1≤ϑ≤2,Z(0)=Z0Z′(0)=Z0+nZ0.

Using Theorem 2.3, equation 8 becomes9 Z(t)=Z0(t)+(Z0+mZ0)+[Θ(t,Z(t))-Θ0(t)]1-ϑM(ϑ)+ϑM(ϑ)Γ(ϑ)×∫0t(t-y)ϑ-1Θ(y-Z(y))dy,1≤ϑ≤2,

whereZ(t)=S(t)V(t)A(t)C(t)R(t),Z0(t)=S0V0A0C0R0,Z0(t)+nZ0(t)=S0+nS0(t)V0+nV0(t)A0+nA0(t)C0+nC0(t)R0+nR0(t),

   10 Θ(t,Z(t)Θ1(t,S,V,A,C,R)Θ2(t,S,V,A,C,R)Θ3(t,S,V,A,C,R)Θ4(t,S,V,A,C,R)Θ5(t,S,V,A,C,R),Θ0(t,Z(t)Θ1(0,S0,V0,A0,C0,R0))Θ2(0,S0,V0,A0,C0,R0)Θ3(0,S0,V0,A0,C0,R0)Θ4(0,S0,V0,A0,C0,R0)Θ5(0,S0,V0,A0,C0,R0).

Using 9 and , defines two different operators F and G, by 911 FZ=Z0(t)+(Z0+mZ0)+[Θ(t,Z(t))-Θ0(t)]1-ϑM(ϑ),GZ=ϑM(ϑ)Γ(ϑ)×∫0ϵ(ϵ-y)ϑ-1Θ(y-Z(y))dy,for1≤ϑ≤2.

Additionally, we assume that the parameters are satisfied:

(D1) If constant β1 and β2, fulfilled as given|Θ(ϵ,Z(ϵ))|≤β1|Z(ϵ)|+β2.

(D2) If constant κ>0∀Z,Z1∈X, fulfilling as12 |Θ(t,Z(t))-Θ(t,Z1(t))|≤κ‖Z-Z1‖.

Theorem 3.1

If (D1) and (D2) hold, therefore 9 must include their root, which implies that the root of the specified system 1 exists. Similarly, if(1-ϑ)κM(ϑ)<1.

Proof

Consider Z1∈B as B={Z∈Y:‖Z‖≤r,r>0} is a near convex set to demonstrate that F is a contraction. Using the research of F as shown in 11, we obtain13 ‖FZ-FZ1‖=(1-ϑ)M(ϑ)maxt∈[0,T]‖Θ(t,Z(t))-Θ(t,Z1(t))‖≤(1-ϑ)pM(ϑ)‖Z-Z1‖.

F is a contraction.

We demonstrate that G has upper and lower bounds and must be continuous to deduce that G demonstrates relative compactness. To do this, we write as follows:

As long as Θ is continuous, Mapping G acts constantly. As an afterthought for u∈B, we put as14 |G(Z)|=supϵ∈[0,T]MM(t)Γ(t)‖∫0t(t-y)ϑ-1Θ(y,Z(y))dy‖≤MM(t)Γ(t)∫0ℑ(ℑ-y)t-1|Θ(y,Z(y))|dy≤MℑtM(t)Γ(t)[β1r+β2].

Therefore, 12 indicates that G has upper and lower bounds. Take the expression t1>t2∈[0,ℑ] as15 |GZ(t1)-GZ(t2)|=tM(t)Γ(t)|∫01(t1-y)ϑ-1Θ(y,Z(y))dy-∫02(t2-y)ϑ-1Θ(y,Z(y))dy|≤[β1r+β2]M(ϑ)Γ(ϑ)[t1ϑ-t2ϑ].

As t1→t2, result 15 finishes, also G is continuous so|GZ(t1)-GZ(t2)|→0,ast1→t2.

Because G has upper and lower bounds and is continuous, it is also bounded, continuous, and uniformly continuous. G exhibits relative compactness and is entirely continuous, according to Arzela-Ascoli results. According to Theorem 3.1, problem 9, which is mentioned above, has a root. □

For one root, we proceed as follows.

Theorem 3.2

Equation 9 has one root due to the requirement (D2), which suggests that the current proposal 1 has just one result.[(1-ϑ)κM(ϑ)+ϑℑϑκM(ϱ)Γ(ϑ)]<1.

Proof

Let define ℑ:Y→Y by16 ℑZ(t)=Z0(t)+[Θ(t,Z(t))-Θ0(t)]1-ϑM(ϑ)+ϑM(ϑ)Γ(ϑ)×∫0t(t-y)ϑ-1Θ(y,Z(y))dy,t∈[0,ℑ].

Let Z,Z1∈Y, then17 ‖ℑZ-ℑZ1‖≤(1-ϑ)M(Γ(ϑ))maxϑ∈[0,ℑ]|Θ(t,Z(t))-Θ(ϑ,Z1(t))|+ϑM(ϑ)Γ(ϑ)maxt∈[0,ℑ]|∫0t(t-y)t-1Θ(y,Z(y))dy-∫0t(t-y)ϑ-1Θ(y,Z1(y))dy|≤[(1-ϑ)κM(ϑ)+ϑℑϑκM(t)Γ(t)]‖Z-Z1‖≤τ‖Z-Z1‖,

where18 τ=[(1-ϑ)κM(ϑ)+ϑℑϑκM(t)Γ(t)].

ℑ is contracted from 17. The solution 9, therefore, only has one root. Simulation 1 has one zero as a result. □

Hyers-Ulam stability

Definition 4.1

When γ≥0∀ϵ>0, as well as any root Y∈C1(G,R) of inequality are present, system 1 is Hyers-Ulam stable43:19 |cDλY(t)-ϑ(t,Y(t))|≤ϵ,t∈G,

there will be a distinct root for problem 1Ξ∈C1(G,R) as20 |Y(t)-Ξ(t)|≤γt,t∈G.

Definition 4.2

A fixed Xi>0,i∈N5 trying to satisfy the following conditions is considered Hypers-Ulam stable for the Caputo fractional integration model described by equation 9. For every γi>0,i∈N5, with|S(t)-1-ϑBϑK1(ϑ,t,S(t))+ϑB(ϑ)Γ(ϑ)×∫0t(t-X)ϑ-1K1(ϑ,X,S(X))dX|≤γ1,|V(t)-1-ϑBϑK2(ϑ,t,V(t))+ϑB(ϑ)Γ(ϑ)×∫0t(t-X)ϑ-1K2(ϑ,X,V(X))dX|≤γ2,|A(t)-1-ϑBϑK3(ϑ,t,V(t))+ϑB(ϑ)Γ(ϑ)×∫0t(t-X)ϑ-1K3(ϑ,X,A(X))dX|≤γ3,|C(t)-1-ϑBϑK4(ϑ,t,C(t))+ϑB(ϑ)Γ(ϑ)×∫0t(t-X)ϑ-1K4(ϑ,X,C(X))dX|≤γ4,|R(t)-1-ϑBϑK5(ϑ,t,R(t))+ϑB(ϑ)Γ(ϑ)×∫0t(t-X)ϑ-1K5(ϑ,X,R(X))dX|≤γ5,

there exist (S^(t),V^(t),A^(t),C^(t),J^(t),R^(t)) which are fulfillingS^(t)=1-ϑBϑK1(ϑ,t,S(t))+ϑB(ϑ)Γ(ϑ)×∫0t(t-X)ϑ-1K1(ϑ,X,S^(X))dX|≤γ1,V^(t)=1-ϑBϑK2(ϑ,t,V(t))+ϑB(ϑ)Γ(ϑ)×∫0t(t-X)ϑ-1K2(ϑ,X,V^(X))dX|≤γ2,A^(t)1-ϑBϑK3(ϑ,t,V(t))+ϑB(ϑ)Γ(ϑ)×∫0t(t-X)ϑ-1K3(ϑ,X,A^(X))dX|≤γ3,C^(t)=1-ϑBϑK4(ϑ,t,C(t))+ϑB(ϑ)Γ(ϑ)×∫0t(t-X)ϑ-1K4(ϑ,X,C^(X))dX|≤γ4,R^(t)=1-ϑBϑK5(ϑ,t,R(t))+ϑB(ϑ)Γ(ϑ)×∫0t(t-X)ϑ-1K5(ϑ,X,R^(X))dX|≤γ5,

21 |S(t)-(^S)(t)|≤ς1γ1,|V(t)-(^V)(t)|≤ς2γ2,|A(t)-(^A)(t)|≤ς3γ3,|C(t)-(^C)(t)|≤ς4γ4,|R(t)-(^R)(t)|≤ς5γ5.

Theorem 4.1

The proposed solution 1 is Hyers-Ulam stable under Y=C[0,ℑ].

Proof

Theorem 5 demonstrates that the given issue 1 has a single root (S(t), V(t), A(t), C(t), R(t)) that satisfies the system’s equation 9. Now comes22 ‖S(t)-S^(t)‖≤1-ϑB(t)‖K1(ϑ,ℑ,S(ϑ))-K1(ϑ,t,S^(t))‖+ϑB(ϑ)Γ(ϑ)∫0t(t-X)ϑ-1‖dX≤[1-ϑB(ϑ)+ϑB(ϑ)Γ(σ)]Ω1‖S-S^‖,‖V(t)-V^(t)‖≤1-ϑB(t)‖K2(ϑ,ℑ,V(ϑ))-K2(ϑ,t,V^(t))‖+ϑB(ϑ)Γ(ϑ)∫0t(t-X)ϑ-1‖dX≤[1-ϑB(ϑ)+ϑB(ϑ)Γ(σ)]Ω1‖V-V^‖,‖A(t)-A^(t)‖≤1-ϑB(t)‖K3(ϑ,ℑ,A(ϑ))-K3(ϑ,t,A^(t))‖+ϑB(ϑ)Γ(ϑ)∫0t(t-X)ϑ-1‖dX≤[1-ϑB(ϑ)+ϑB(ϑ)Γ(σ)]Ω1‖CA^‖,‖C(t)-C^(t)‖≤1-ϑB(t)‖K4(ϑ,ℑ,C(ϑ))-K4(ϑ,t,C^(t))‖+ϑB(ϑ)Γ(ϑ)∫0t(t-X)ϑ-1‖dX≤[1-ϑB(ϑ)+ϑB(ϑ)Γ(σ)]Ω1‖CC^‖,‖R(t)-R^(t)‖≤1-ϑB(t)‖K5(ϑ,ℑ,R(ϑ))-K5(ϑ,t,R^(t))‖+ϑB(ϑ)Γ(ϑ)∫0t(t-X)ϑ-1‖dX≤[1-ϑB(ϑ)+ϑB(ϑ)Γ(σ)]Ω1‖RR^‖.

Taking, γi=Ωi,Δi=1-ϑB(ϑ)+ϑB(ϑ)Γ(ϑ), this implies23 ‖S(t)-S^(t)‖≤γ1Δ1.

Thus, we get24 ‖V(t)-V^(t)‖≤γ2Δ2‖A(t)-A^(t)‖≤γ2Δ2‖C(t)-C^(t)‖≤γ4Δ4‖R(t)-R^(t)‖≤γ5Δ5.

The fractional model 1 is Hyers-Ulam stable because the issue 9 is Hyers-Ulam by equations 23 and 24. □

Numerical solutions

In this part of the study, we will find the approximations to the model’s solutions 1. The popular Fractional Adam-Bashforth approach has been used to enumerate the noninteger order integration18. Using beginning assumptions and the operator I0ϱ, we modify the noninteger order model 9 into fractional Volterra type integration equations:S(t)=S0+(S0+mS0)+1-ϑM(ϑ)K1(t,S(t))+ϑΓ(ϑ)M(ϑ)∫0tK1(σ,S(σ))(t-σ)ϑ-1dσ,V(t)=V0+(V0+mV0)+1-ϑM(ϑ)K2(t,V(t))+ϑΓ(ϑ)M(ϑ)∫0tK2(σ,V(σ))(t-σ)ϑ-1dσ,A(t)=A0+(A0+mA0)+1-ϑM(ϑ)K3(t,A(t))+ϑΓ(ϑ)M(ϑ)∫0tK3(σ,A(σ))(t-σ)ϑ-1dσ,C(t)=C0+(C0+mC0)+1-ϑM(ϑ)K4(t,S(t))+ϑΓ(ϑ)M(ϑ)∫0tK4(σ,C(σ))(t-σ)ϑ-1dσ,R(t)=R0+(R0+mR0)+1-ϑM(ϑ)K5(t,S(t))+ϑΓ(ϑ)M(ϑ)∫0tK5(σ,R(σ))(t-σ)ϑ-1dσ,

We obtain for t=tm+1:25 S(tm+1)=S0+(S0+mS0)+1-ϑM(ϑ)K1(tm,S(tm))+∑k=0mϑΓ(ϑ)M(ϑ)∫tkk+1K1(σ,S(σ))(tm+1-σ)ϑ-1dσ,V(tm+1)=V0+(V0+mV0)+1-ϑM(ϑ)K2(tm,V(tm))+∑k=0mϑΓ(ϑ)M(ϑ)∫tkk+1K2(σ,V(σ))(tm+1-σ)ϑ-1dσ,A(tm+1)=A0+(A0+mA0)+1-ϑM(ϑ)K3(tm,A(tm))+∑k=0mϑΓ(ϑ)M(ϑ)∫tkk+1K3(σ,A(σ))(tm+1-σ)ϑ-1dσ,C(tm+1)=C0+(C0+mC0)+1-ϑM(ϑ)K4(tm,C(tm))+∑k=0mϑΓ(ϑ)M(ϑ)∫tkk+1K4(σ,C(σ))(tm+1-σ)ϑ-1dσ,R(tm+1)=R0+(R0+mR0)+1-ϑM(ϑ)K5(tm,R(tm))+∑k=0mϑΓ(ϑ)M(ϑ)∫tkk+1K5(σ,R(σ))(tm+1-σ)ϑ-1dσ,

By using the bi-step lagrange’s polynomial interpolate to approximate K1 to K5 in [tk,tk+1] and re-inserting it into 25, we obtainS(tm+1)=S0+(S0+mS0)+1-ϑM(ϑ)K1(tm,S(tm))+ϑΓ(ϑ)M(ϑ)×∑k=0m(K1(tk,Sk)h∫tkk+1(σ-tk-1)(tm+1-σ)ϑ-1dσ),V(tm+1)=V0+(V0+mV0)+1-ϑM(ϑ)K2(tm,V(tm))+ϑΓ(ϑ)M(ϑ)×∑k=0m(K1(tk,Vk)h∫tkk+1(σ-tk-1)(tm+1-σ)ϑ-1dσ),A(tm+1)=A0+(A0+mA0)+1-ϑM(ϑ)K3(tm,A(tm))+ϑΓ(ϑ)M(ϑ)×∑k=0m(K1(tk,Ak)h∫tkk+1(σ-tk-1)(tm+1-σ)ϑ-1dσ),

26 C(tm+1)=C0+(C0+mC0)+1-ϑM(ϑ)K4(tm,C(tm))+ϑΓ(ϑ)M(ϑ)×∑k=0m(K1(tk,Ck)h∫tkk+1(σ-tk-1)(tm+1-σ)ϑ-1dσ),R(tm+1)=R0+(R0+mR0)+1-ϑM(ϑ)K5(tm,R(tm))+ϑΓ(ϑ)M(ϑ)×∑k=0m(K1(tk,Rk)h∫tkk+1(σ-tk-1)(tm+1-σ)ϑ-1dσ),

Upon combining it with the term in 26 and inserting it once more, we obtainS(tm+1)=S0+(S0+mS0)+1-ϑM(ϑ)K1(tm,S(tm))+ϑΓ(ϑ)M(ϑ)×∑k=0m(hϑK1(tk,Sk)Γ(ϑ+2)((m+1-k)ϑ(m-k-2+ϑ)-(-k+m)ϑ(2ϑ+m+2-k))-hϑK1(tk-1,Sk-1)Γ(ϑ+2)×((1-k+m)ϑ+1-(m-k)ϑ(ϑ+m+1-k))),V(tm+1)=V0+(V0+mV0)+1-ϑM(ϑ)K1(tm,V(tm))+ϑΓ(ϑ)M(ϑ)×∑k=0m(hϑK1(tk,Vk)Γ(ϑ+2)((m+1-k)ϑ(m-k-2+ϑ)-(-k+m)ϑ(2ϑ+m+2-k))-hϑK1(tk-1,Vk-1)Γ(ϑ+2)×((1-k+m)ϱ+1-(m-k)ϑ(ϑ+m+1-k))),A(tm+1)=A0+(A0+mS0)+1-ϑM(ϑ)K1(tm,A(tm))+ϑΓ(ϑ)M(ϑ)×∑k=0m(hϑK1(tk,Ak)Γ(ϑ+2)((m+1-k)ϑ

(m-k-2+ϑ)-(-k+m)ϑ(2ϑ+n+2-k))-hϑK1(tk-1,Ak-1)Γ(ϑ+2)×((1-k+m)ϑ+1-(m-k)ϑ(ϑ+m+1-k))),C(tm+1)=C0+(C0+mC0)+1-ϑM(ϑ)K1(tm,C(tm))+ϑΓ(ϑ)M(ϑ)×∑k=0m(hϑK1(tk,Ck)Γ(ϑ+2)((m+1-k)ϑ(m-k-2+ϑ)-(-k+m)ϑ(2ϑ+m+2-k))-hϑK1(tk-1,Ck-1)Γ(ϑ+2)×((1-k+m)ϑ+1-(m-k)ϑ(ϑ+m+1-k))),R(tm+1)=R0+(R0+mR0)+1-ϑM(ϑ)K1(tm,S(tm))+ϑΓ(ϑ)M(ϑ)×∑k=0m(hϑK1(tk,Rk)Γ(ϑ+2)((m+1-k)ϑ(m-k-2+ϑ)-(-k+m)ϑ(2ϑ+n+2-k))-hϑK1(tk-1,Rk-1)Γ(ϑ+2)×((1-k+m)ϑ+1-(m-k)ϑ(ϑ+m+1-k))),

where27 K1=κΛ-2βS(t)A(t)s(t)+A(t)+V(t)ρ-S(t)(μ+η),K2=Λ(1-κ)+μS(t)-βV(t)A(t)(1-τ1)2V(t)+A(t)-V(t)(ρ+η),K3=2βS(t)A(t)S(t)+A(t)+2(1-τ1)βV(t)A(t)V(t)+A(t)-(α1+α2+η)A(t),K4=α1A(t)-(γ1+γ2+η)C(t),K5=γ1C(t)+α2A(t)-ηR(t).

Discussions of numerical simulation

In this part, we use the approximate solution technique to demonstrate the evolving behaviors of the Caputo derivative generalized HBV illness issue 1. The initial approximations were S(0) = 1.45; V(0) = 1.154; A(0) = 1.125; C(0) = 1.105; and R(0) = 1.105. The parameters from 17 that we selected are shown in Table 2 and used in the numerical results. The timing scale is suggested to be 25 in the graphical numerical iterations. With the use of Matlab’s Adams-Bashforth Moulton iteration, the general fractional predictor-corrector method is used to create the graphical depiction. The agents are initially generated for a range of noninteger order values to understand the dynamics of the populations of the above agents based on their infection circumstances. Figs. 2-6 show the expatriation. In the numerical simulation procedures, the populations of Susceptible and Vaccinated people grow exponentially. A significant increase is observed in the population of those who have received vaccinations, whereas acute cases of chronic infection and recovery cases decrease. According to biological ideas, this effect results from lengthy immunization periods and a wealth of public knowledge. The graphical representation demonstrates that for the taken-in HBV system, the results of globalized operators appear to be more accurate. Figures 2-6 illustrate how the System 1 changes its behaviors when given different values of = 1.95, 1.85, 1.75, 1.65, and 1.55.Table 2 Parameters used in the numerical simulations of model 1.

Parameters description	Numerical value	
Λ	0.8	
μ	0.073	
ρ	0.005	
β	0.0576	
γ1	0.0734	
γ2	0.0734	
κ	0.0009	
λ	0.007	
α1	0.054	
α2	0.0003	
τ	0.967	
S(0)	1.45	
V(0)	1.154	
A(0)	1.125	
C(0)	1.105	
R(0)	1.105	

Fig. 2 Existence of susceptible human for ϑ∈(1,2).

The model’s data set is divided into three categories: 70% for training, 15% for testing, and 15% for validation. Every category has a unique fractional order. Figure 3a represent the approach of ANN with Adam Bashforth method. Figure 3b illustrates the performance of the considered model at epoch 29, with mean square errors of 1.4698e-11. The error histogram can be shown in Fig. 3c, where -4.1e-07 is the best value we obtain in this example. Figure 3d displays the errors in the exact figure and the best fit of the training and testing data. Figure 4 shows the regression of the model under consideration for all data, testing data, and training data. The data is clearly on the regression line, demonstrating that the generated solution is accurately trained. The value of R is about 1.Fig. 3 Statistically, dynamics of the considered model with ANN (a) comparison, (b) mean square error, (c) error histogram, and (d) training fit.

Fig. 4 Dynamically representation of the regression with ANN for the considered system.

as shown in Fig. 5a, b illustrates the performance of the considered model at epoch 1000, with mean square errors of 3.1295e-05. The error histogram can be shown in Fig. 5c, where -0.00197 is the best value we obtain in this example. Figure 5d displays the best fit of the training and testing data and errors, as indicated in the exact figure. Figure 6 shows the regression of the model under consideration for all data, testing data, and training data. The data is clearly on the regression line, demonstrating that the generated solution is accurately trained. The value of R is about 1.Fig. 5 Statistically, dynamics of the considered model with ANN (a) comparison, (b) mean square error, (c) error histogram, and (d) training fit.

Fig. 6 Dynamically representation of the regression with ANN for the considered system.

as shown in Fig. 7a represent the approach of ANN with Adam Bashforth method. Figure 7b illustrates the performance of the considered model at epoch 1000, with mean square errors of 2.3443e-08. Figure 7c displays the error histogram, and 0.001564 is the best value we obtain in this example. Figure 7d displays the errors in the same figure and the best fit of the training and testing data. Figure 8 displays the regression of the model under consideration for all data, testing data, and training data. The data is clearly on the regression line, demonstrating that the generated solution is accurately trained. The value of R is about 1.Fig. 7 Statistically, dynamics of the considered model with ANN (a) comparison, (b) mean square error, (c) error histogram, and (d) training fit.

Fig. 8 Dynamically representation of the regression with ANN for the considered system.

as shown in Fig. 9a represent the approach of ANN with Adam Bashforth method. Figure 9b illustrates the performance of the considered model at epoch 193, with mean square errors of 1.4342e-12. Figure 9c displays the error histogram, and -61e-08 is the best value we obtain. Figure 9d displays the errors stated in the same figure and the best fit of the training and testing data. Figure 10 displays the regression of the model under consideration for all data, testing data, and training data. The data is clearly on the regression line, demonstrating that the generated solution is accurately trained. The value of R is about 1.Fig. 9 Statistically, dynamics of the considered model with ANN (a) comparison, (b) mean square error, (c) error histogram, and (d) training fit.

Fig. 10 Dynamically representation of the regression with ANN for the considered system.

as shown in Fig. 11a represent the approach of ANN with Adam Bashforth method. With mean square errors of 9.6556e-14 at epoch 8, the performance of the studied model is displayed in Fig. 11b, c displays the error histogram, and -8.3e-08 is the best value we obtain in this situation. Figure 11d displays the best fit of the training and testing data and errors in the same figure. The model under consideration’s regression is displayed in Fig. 12 for all data, testing, and training. These graphs show that the data is on the regression line, indicating that the derived solution was accurately trained. The value of R is also approximately 1.Fig. 11 Statistically, dynamics of the considered model with ANN (a) comparison, (b) mean square error, (c) error histogram, and (d) training fit.

Fig. 12 Dynamically representation of the regression with ANN for the considered system.

Application: informing hepatitis B control strategies

The comprehensive analysis of the fractional-order HBV transmission model offer important insights that can inform the development of effective public health strategies for controlling hepatitis B virus infection.

One key insight from the model is the identification of critical thresholds for the basic reproduction number (R0) and the fractional derivative order (α). The analysis revealed that when R0 is below a critical value and the fractional order α is within a certain range, the disease-free equilibrium is stable, indicating the potential for successful disease control. This information can guide public health authorities in setting targeted goals for reducing HBV transmission rates and designing vaccination programs to achieve the necessary reduction in R0 to drive the system towards disease elimination.

Furthermore, the model’s predictions on the long-term dynamics of HBV infection, including the existence of endemic equilibria and the potential for backward bifurcations, can inform the design of treatment protocols and public awareness campaigns. By understanding the conditions under which the disease may persist or resurge, interventions can be tailored to address the specific challenges posed by HBV transmission in different epidemiological settings.

The model’s ability to capture the non-linear and memory-dependent nature of HBV infection dynamics through the use of fractional derivatives also presents opportunities for more accurate forecasting of disease trends. By calibrating the model parameters to match real-world HBV incidence data, public health authorities can leverage the model’s predictive capabilities to anticipate future disease patterns and allocate resources more effectively for prevention, early detection, and treatment programs.

Overall, the insights gained from the fractional HBV transmission model can serve as a valuable tool for informing the development and implementation of comprehensive hepatitis B control strategies, ultimately contributing to the reduction of the global burden of this infectious disease.

Conclusion

The key findings and contributions of this study are:Developed a fractional model to describe the kinetics of hepatitis B (HBV) transmission, incorporating the Caputo derivative and vaccine effects.

Determined the essential reproductive value and equilibria of the proposed fractional HBV model.

Conducted a qualitative analysis of the problem’s approximative root using Fixed Point Theory.

Employed the Adams-Bashforth predictor-corrector scheme to evaluate the iterative approximate technique for the fractional system.

Provided graphical representations to compare the results for different noninteger orders.

Utilized Artificial Neural Network (ANN) techniques to partition the dataset into training, testing, and validation categories, with a comprehensive examination of the dataset’s characteristics and behaviors within these divisions.

Future research directions include:Incorporating additional factors specific to hepatitis.

Conducting sensitivity analysis to identify influential parameters.

Exploring spatio-temporal and multi-scale modeling approaches.

Investigating optimal control strategies and interventions.

Validating the model and quantifying uncertainties.

These avenues of research can further enhance our understanding of hepatitis transmission dynamics and inform effective decision-making to combat the disease.

Acknowledgements

This research is funded by “Researchers Supporting Project number (RSPD2024R733), King Saud University, Riyadh, Saudi Arabia

Author contributions

All author’s contributed equally to the writing of this paper. All author’s read and approved the final manuscript.

Data availability

The datasets used and/or analyzed during the current study available from the corresponding author on reason-able request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Hu M Chen W Assessment of total economic burden of chronic hepatitis B (CHB)-related diseases in Beijing and Guangzhou, China Value Health 2009 12 S89 S92 10.1111/j.1524-4733.2009.00636.x 20586991
Hu, M. & Chen, W. Assessment of total economic burden of chronic hepatitis B (CHB)-related diseases in Beijing and Guangzhou, China. Value Health 12, S89–S92 (2009).20586991
2. Nowak M May RM Virus Dynamics: Mathematical Principles of Immunology and Virology: Mathematical Principles of Immunology and Virology 2000 Oxford Oxford University Press
Nowak, M. & May, R. M. Virus Dynamics: Mathematical Principles of Immunology and Virology: Mathematical Principles of Immunology and Virology (Oxford University Press, Oxford, 2000).
3. Wodarz D May RM Nowak MA The role of antigen-independent persistence of memory cytotoxic T lymphocytes Int. Immunol. 2000 12 4 467 477 10.1093/intimm/12.4.467 10744648
Wodarz, D., May, R. M. & Nowak, M. A. The role of antigen-independent persistence of memory cytotoxic T lymphocytes. Int. Immunol. 12(4), 467–477 (2000).10744648
4. Nowak MA Bonhoeffer S Hill AM Boehme R Thomas HC McDade H Viral dynamics in hepatitis B virus infection Proc. Natl. Acad. Sci. 1996 93 9 4398 4402 10.1073/pnas.93.9.4398 8633078
Nowak, M. A. et al. Viral dynamics in hepatitis B virus infection. Proc. Natl. Acad. Sci. 93(9), 4398–4402 (1996).8633078
5. Din A Li Y Yusuf A Delayed hepatitis B epidemic model with stochastic analysis Chaos Solitons Fractals 2021 146 110839 10.1016/j.chaos.2021.110839
Din, A., Li, Y. & Yusuf, A. Delayed hepatitis B epidemic model with stochastic analysis. Chaos Solitons Fractals 146, 110839 (2021).
6. Ullah S Altaf Khan M Farooq M A new fractional model for the dynamics of the hepatitis B virus using the Caputo-Fabrizio derivative Eur. Phys. J. Plus 2018 133 1 14 10.1140/epjp/i2018-12072-4
Ullah, S., Altaf Khan, M. & Farooq, M. A new fractional model for the dynamics of the hepatitis B virus using the Caputo-Fabrizio derivative. Eur. Phys. J. Plus 133, 1–14 (2018).
7. Khan A Zarin R Hussain G Usman AH Humphries UW Gomez-Aguilar JF Modeling and sensitivity analysis of HBV epidemic model with convex incidence rate Res. Phys. 2021 22 103836
Khan, A. et al. Modeling and sensitivity analysis of HBV epidemic model with convex incidence rate. Res. Phys. 22, 103836 (2021).
8. Arif MS Raza A Rafiq M Bibi M A reliable numerical analysis for stochastic hepatitis B virus epidemic model with the migration effect Iran. J. Sci. Technol. Trans. A Sci. 2019 43 2477 2492 10.1007/s40995-019-00726-0
Arif, M. S., Raza, A., Rafiq, M. & Bibi, M. A reliable numerical analysis for stochastic hepatitis B virus epidemic model with the migration effect. Iran. J. Sci. Technol. Trans. A Sci. 43, 2477–2492 (2019).
9. Lu Q Stability of SIRS system with random perturbations Physica A 2009 388 18 3677 3686 10.1016/j.physa.2009.05.036
Lu, Q. Stability of SIRS system with random perturbations. Physica A 388(18), 3677–3686 (2009).
10. Wang F Khan MN Ahmad I Ahmad H Abu-Zinadah H Chu YM Numerical solution of traveling waves in chemical kinetics: Time-fractional fishers equations Fractals 2022 30 02 2240051 10.1142/S0218348X22400515
Wang, F. et al. Numerical solution of traveling waves in chemical kinetics: Time-fractional fishers equations. Fractals 30(02), 2240051 (2022).
11. Chu, Y. M., Shankaralingappa, B. M., Gireesha, B. J., Alzahrani, F., Khan, M. I., & Khan, S. U. (2022). RETRACTED: Combined impact of Cattaneo-Christov double diffusion and radiative heat flux on bio-convective flow of Maxwell liquid configured by a stretched nano-material surface.
12. He ZY Abbes A Jahanshahi H Alotaibi ND Wang Y Fractional-order discrete-time SIR epidemic model with vaccination: Chaos and complexity Mathematics 2022 10 2 165 10.3390/math10020165
He, Z. Y., Abbes, A., Jahanshahi, H., Alotaibi, N. D. & Wang, Y. Fractional-order discrete-time SIR epidemic model with vaccination: Chaos and complexity. Mathematics 10(2), 165 (2022).
13. Zhao TH Castillo O Jahanshahi H Yusuf A Alassafi MO Alsaadi FE Chu YM A fuzzy-based strategy to suppress the novel coronavirus (2019-NCOV) massive outbreak Appl. Comput. Math. 2021 20 1 160 176
Zhao, T. H. et al. A fuzzy-based strategy to suppress the novel coronavirus (2019-NCOV) massive outbreak. Appl. Comput. Math. 20(1), 160–176 (2021).
14. Atangana, A., & Baleanu, D. (2016). New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model. arXiv preprint arXiv:1602.03408.
15. Atangana A Koca I Chaos in a simple nonlinear system with Atangana-Baleanu derivatives with fractional order Chaos Solitons Fractals 2016 89 447 454 10.1016/j.chaos.2016.02.012
Atangana, A. & Koca, I. Chaos in a simple nonlinear system with Atangana-Baleanu derivatives with fractional order. Chaos Solitons Fractals 89, 447–454 (2016).
16. Din A Shah K Seadawy A Alrabaiah H Baleanu D On a new conceptual mathematical model dealing the current novel coronavirus-19 infectious disease Res. Phys. 2020 19 103510
Din, A., Shah, K., Seadawy, A., Alrabaiah, H. & Baleanu, D. On a new conceptual mathematical model dealing the current novel coronavirus-19 infectious disease. Res. Phys. 19, 103510 (2020).
17. Li, X. P. et al. A new Hepatitis B model in light of asymptomatic carriers and vaccination study through Atangana-Baleanu derivative. Res. Phys. 29, 104603 (2021).
18. Zafar ZUA Tunç C Ali N Zaman G Thounthong P Dynamics of an arbitrary order model of toxoplasmosis ailment in human and cat inhabitants J. Taibah Univ. Sci. 2021 15 1 882 896 10.1080/16583655.2021.1990603
Zafar, Z. U. A., Tunç, C., Ali, N., Zaman, G. & Thounthong, P. Dynamics of an arbitrary order model of toxoplasmosis ailment in human and cat inhabitants. J. Taibah Univ. Sci. 15(1), 882–896 (2021).
19. Zafar ZUA Ali N Baleanu D Dynamics and numerical investigations of a fractional-order model of toxoplasmosis in the population of human and cats Chaos Solitons Fractals 2021 151 111261 10.1016/j.chaos.2021.111261
Zafar, Z. U. A., Ali, N. & Baleanu, D. Dynamics and numerical investigations of a fractional-order model of toxoplasmosis in the population of human and cats. Chaos Solitons Fractals 151, 111261 (2021).
20. Zafar ZUA DarAssi MH Ahmad I Assiri TA Meetei MZ Khan MA Hassan AM Numerical simulation and analysis of the stochastic hiv/aids model in fractional order Res. Phys. 2023 53 106995
Zafar, Z. U. A. et al. Numerical simulation and analysis of the stochastic hiv/aids model in fractional order. Res. Phys. 53, 106995 (2023).
21. Zafar ZUA Zaib S Hussain MT Tunç C Javeed S Analysis and numerical simulation of tuberculosis model using different fractional derivatives Chaos Solitons Fractals 2022 160 112202 10.1016/j.chaos.2022.112202
Zafar, Z. U. A., Zaib, S., Hussain, M. T., Tunç, C. & Javeed, S. Analysis and numerical simulation of tuberculosis model using different fractional derivatives. Chaos Solitons Fractals 160, 112202 (2022).
22. Ahmad Z Bonanomi G di Serafino D Giannino F Transmission dynamics and sensitivity analysis of pine wilt disease with asymptomatic carriers via fractal-fractional differential operator of Mittag-Leffler kernel Appl. Numer. Math. 2023 185 446 465 10.1016/j.apnum.2022.12.004
Ahmad, Z., Bonanomi, G., di Serafino, D. & Giannino, F. Transmission dynamics and sensitivity analysis of pine wilt disease with asymptomatic carriers via fractal-fractional differential operator of Mittag-Leffler kernel. Appl. Numer. Math. 185, 446–465 (2023).
23. Shah K Wenqi L Raezah AA Khan N Khan SU Ozair M Ahmad Z Unraveling pine wilt disease: Comparative study of stochastic and deterministic model using spectral method Expert Syst. Appl. 2024 240 122407 10.1016/j.eswa.2023.122407
Shah, K. et al. Unraveling pine wilt disease: Comparative study of stochastic and deterministic model using spectral method. Expert Syst. Appl. 240, 122407 (2024).
24. Sinan M Shah K Kumam P Mahariq I Ansari KJ Ahmad Z Shah Z Fractional order mathematical modeling of typhoid fever disease Res. Phys. 2022 32 105044
Sinan, M. et al. Fractional order mathematical modeling of typhoid fever disease. Res. Phys. 32, 105044 (2022).
25. Sinan M Ahmad H Ahmad Z Baili J Murtaza S Aiyashi MA Botmart T Fractional mathematical modeling of malaria disease with treatment & insecticides Res. Phys. 2022 34 105220
Sinan, M. et al. Fractional mathematical modeling of malaria disease with treatment & insecticides. Res. Phys. 34, 105220 (2022).
26. Sami A Ali A Shafqat R Pakkaranang N Ur Rahmamn M Analysis of food chain mathematical model under fractal fractional Caputo derivative Math. Biosci. Eng. 2023 20 2 2094 2109 10.3934/mbe.2023097 36899524
Sami, A., Ali, A., Shafqat, R., Pakkaranang, N. & Ur Rahmamn, M. Analysis of food chain mathematical model under fractal fractional Caputo derivative. Math. Biosci. Eng. 20(2), 2094–2109 (2023).36899524
27. Anjam YN Shafqat R Sarris IE Ur Rahman M Touseef S Arshad M A fractional order investigation of smoking model using Caputo-Fabrizio differential operator Fractal Fract. 2022 6 11 623 10.3390/fractalfract6110623
Anjam, Y. N. et al. A fractional order investigation of smoking model using Caputo-Fabrizio differential operator. Fractal Fract. 6(11), 623 (2022).
28. Li B Eskandari Z Dynamical analysis of a discrete-time SIR epidemic model J. Franklin Inst. 2023 360 12 7989 8007 10.1016/j.jfranklin.2023.06.006
Li, B. & Eskandari, Z. Dynamical analysis of a discrete-time SIR epidemic model. J. Franklin Inst. 360(12), 7989–8007 (2023).
29. Zhu X Xia P He Q Ni Z Ni L Ensemble classifier design based on perturbation binary salp swarm algorithm for classification Comput. Model. Eng. Sci. 2023 135 1 653 671
Zhu, X., Xia, P., He, Q., Ni, Z. & Ni, L. Ensemble classifier design based on perturbation binary salp swarm algorithm for classification. Comput. Model. Eng. Sci. 135(1), 653–671 (2023).
30. Zafar ZUA Hussain MT Inc M Baleanu D Almohsen B Oke AS Javeed S Fractional-order dynamics of human papillomavirus Res. Phys. 2022 34 105281
Zafar, Z. U. A. et al. Fractional-order dynamics of human papillomavirus. Res. Phys. 34, 105281 (2022).
31. Ahmad Z Bonanomi G Cardone A Iuorio A Toraldo G Giannino F Fractal-fractional sirs model for the disease dynamics in both prey and predator with singular and nonsingular kernels J. Biol. Syst. 2024 6 1 34 10.1142/S0218339024400035
Ahmad, Z. et al. Fractal-fractional sirs model for the disease dynamics in both prey and predator with singular and nonsingular kernels. J. Biol. Syst. 6, 1–34 (2024).
32. Shepard CW Simard EP Finelli L Fiore AE Bell BP Hepatitis B virus infection: Epidemiology and vaccination Epidemiol. Rev. 2006 28 1 112 125 10.1093/epirev/mxj009 16754644
Shepard, C. W., Simard, E. P., Finelli, L., Fiore, A. E. & Bell, B. P. Hepatitis B virus infection: Epidemiology and vaccination. Epidemiol. Rev. 28(1), 112–125 (2006).16754644
33. Zhang J Zou S Giulivi A Epidemiology of hepatitis B in Canada Can. J. Infect. Dis. Med. Microbiol. 2001 12 6 345 350 10.1155/2001/790915
Zhang, J., Zou, S. & Giulivi, A. Epidemiology of hepatitis B in Canada. Can. J. Infect. Dis. Med. Microbiol. 12(6), 345–350 (2001).
34. Din A Li Y Yusuf A Ali AI Caputo type fractional operator applied to Hepatitis B system Fractals 2022 30 01 2240023 10.1142/S0218348X22400230
Din, A., Li, Y., Yusuf, A. & Ali, A. I. Caputo type fractional operator applied to Hepatitis B system. Fractals 30(01), 2240023 (2022).
35. Din A Li Y Stationary distribution extinction and optimal control for the stochastic hepatitis B epidemic model with partial immunity Phys. Scr. 2021 96 7 074005 10.1088/1402-4896/abfacc
Din, A. & Li, Y. Stationary distribution extinction and optimal control for the stochastic hepatitis B epidemic model with partial immunity. Phys. Scr. 96(7), 074005 (2021).
36. Alnahdi AS Shafqat R Niazi AUK Jeelani MB Pattern formation induced by fuzzy fractional-order model of COVID-19 Axioms 2022 11 7 313 10.3390/axioms11070313
Alnahdi, A. S., Shafqat, R., Niazi, A. U. K. & Jeelani, M. B. Pattern formation induced by fuzzy fractional-order model of COVID-19. Axioms 11(7), 313 (2022).
37. Abuasbeh K Shafqat R Alsinai A Awadalla M Analysis of the mathematical modelling of COVID-19 by using mild solution with delay Caputo operator Symmetry 2023 15 2 286 10.3390/sym15020286
Abuasbeh, K., Shafqat, R., Alsinai, A. & Awadalla, M. Analysis of the mathematical modelling of COVID-19 by using mild solution with delay Caputo operator. Symmetry 15(2), 286 (2023).
38. Liu P Din A Zarin R Numerical dynamics and fractional modeling of hepatitis B virus model with non-singular and non-local kernels Res. Phys. 2022 39 105757
Liu, P., Din, A. & Zarin, R. Numerical dynamics and fractional modeling of hepatitis B virus model with non-singular and non-local kernels. Res. Phys. 39, 105757 (2022).
39. Caputo M Fabrizio M A new definition of fractional derivative without singular kernel Prog. Fract. Diff. Appl. 2015 1 2 73 85
Caputo, M. & Fabrizio, M. A new definition of fractional derivative without singular kernel. Prog. Fract. Diff. Appl. 1(2), 73–85 (2015).
40. Caputo M Fabrizio M Applications of new time and spatial fractional derivatives with exponential kernels Prog. Fract. Diff. Appl. 2016 2 1 1 11 10.18576/pfda/020101
Caputo, M. & Fabrizio, M. Applications of new time and spatial fractional derivatives with exponential kernels. Prog. Fract. Diff. Appl. 2(1), 1–11 (2016).
41. El-Saka HAA The fractional-order SIS epidemic model with variable population size J. Egyptian Math. Soc. 2014 22 1 50 54 10.1016/j.joems.2013.06.006
El-Saka, H. A. A. The fractional-order SIS epidemic model with variable population size. J. Egyptian Math. Soc. 22(1), 50–54 (2014).
42. Din A Li Y Yusuf A Ali AI Caputo type fractional operator applied to Hepatitis B system Fractals 2022 30 01 2240023 10.1142/S0218348X22400230
Din, A., Li, Y., Yusuf, A. & Ali, A. I. Caputo type fractional operator applied to Hepatitis B system. Fractals 30(01), 2240023 (2022).
43. Ali Z Zada A Shah K Ulam stability to a toppled systems of nonlinear implicit fractional order boundary value problem Bound. Value Problems 2018 2018 1 16
Ali, Z., Zada, A. & Shah, K. Ulam stability to a toppled systems of nonlinear implicit fractional order boundary value problem. Bound. Value Problems 2018, 1–16 (2018).
