A Fractional-Order Epidemic Model with Quarantine Class and Nonmonotonic Incidence: Modeling and Simulations.

In any outbreak of infectious disease, the timely quarantine of infected individuals along with preventive measures strategy are the crucial methods to control new infections in the population. Therefore, this study aims to provide a novel fractional Caputo derivative-based susceptible-infected-quarantined-recovered-susceptible epidemic mathematical model along with a nonmonotonic incidence rate of infection. A new quarantined individual compartment is incorporated into the susceptible-infected-recovered-susceptible compartmental model by dividing the total population into four subpopulations. The nonmonotonic incidence rate of infection is considered as Monod-Haldane functional type to understand the psychological effects in the population. Qualitative analysis of the study shows that the model solutions are well-posed i.e., they are nonnegative and bounded in an attractive region. It is revealed that the model has two equilibria, namely, disease-free (DFE) and endemic (EE). The stability analysis of equilibria is investigated for local as well as global behaviors. Mathematical analysis of the model reveals that DFE is locally asymptotically stable when the basic reproduction number ( R 0 ) is lower than one. The basic reproduction number R 0 is computed using the next-generation matrix method. The existence of EE is shown and it is investigated that EE is locally asymptotically stable when R 0 > 1 under some appropriate conditions. Moreover, the global stability behaviors of DFE and EE are analyzed under some conditions using R 0 . Finally, some numerical simulations are performed to interpret the theoretical findings.

Introduction

Mathematical models are playing a crucial role in understanding and controlling the spread of infectious diseases since a very long time. Epidemiological researchers have developed numerous models to understand the dynamics of infectious diseases, such as and many others (where stands for susceptible, for alert, for Exposed, for vaccinated, for Infected, and for recovered individuals) (Gumel et al. 2006; Goel et al. 2020a, 2020b; Kumar and Nilam 2019; Kumar et al. 2019; Haung et al. 2010; Goel and Nilam 2019; Alexander et al. 2004; Xu and Ma 2009; Wang 2002; Michael et al. 1999; Zhou and Wang 2022). In mathematical epidemiological literature, many studies are based on classical and integer-order derivatives, which are taken as a vital tool for model formulation; for example, see Gumel et al. (2006), Goel et al. (2020a), Kumar and Nilam (2019), Kumar et al. (2019), Alexander et al. (2004) and Wang (2002). In the last few years, researchers are focusing on developing more realistic models using Fractional-order derivatives instead of integer order derivatives which provide more rich dynamics of epidemic models; for example, see Matignon (1996), Rostamy and Mottaghi (2016), Ahmed et al. (2006), Ye and Xu (2019), Khan et al. (2020), Swati (2022), Cui et al. (2022), Chatterjee et al. (2022) and Zhou and Wang (2022).

Fractional calculus is a special branch of mathematics that deals with derivatives of fractional order, unlike the ordinary derivative which is a local operator in nature. The fractional-order derivative has a vital property called the memory effect. As we all know that every biological phenomenon has long-range temporal memory, therefore, the use of the fractional-order derivative in epidemic models instead of the integer-order derivative provides the systems with long-term memory, which has been proven more realistic. There are several definitions of fractional order derivatives (Camargo and Oliveira 2015; Saad et al. 2018). The majorly used definition of fractional order derivative is proposed by Caputo in 1967. Caputo fractional-order derivative-based differential equations system is one of the most used in the field of mathematical modeling. In Caputo derivative-based equation systems, the initial conditions can be expressed in the same manner as in the integer-order based differential equations system. Also, the Caputo fractional derivative has the same property as an integer-order derivative, i.e., the derivative of a constant is zero. This makes the Caputo fractional derivative more popular and useful over other fractional derivatives such as Riemann-Liouville, Grunwald-Letnikov fractional derivatives, etc. Below is the definition and expression of the Caputo fractional order derivative (Lu and Zhu 2018; Ye and Xu 2019):

“The Caputo fractional derivative of order of a function is defined as where is the Gamma function and is the positive integer such that . When , , we have.

The spread of any infectious disease mainly depends on the incidence rate. In epidemiology, the number of individuals who become infected per unit of time is known as the incidence rate (Kermack and Mckendrik 1927; Kumar and Nilam 2021; Goel and Nilam 2019). Kermack and McKendrick (1927) introduced the first infectious disease model using a bilinear incidence rate of the form . This bilinear incidence rate is based on the law of mass action, which is inconsiderable for a large population. Therefore, in lateral years, many researchers have modified this bilinear incidence rate to the nonlinear incidence rates such as the Holling type II, Crowley–Martin functional type, Beddington–DeAngelish functional type, etc., to study the dynamics of infectious diseases for a large population (Kumar et al. 2019; Goel and Nilam 2019; Dubey et al. 2016; Shi et al. 2011; Liu et al. 1987). In 1986 Liu et al. (Liu et al. 1986) suggested a general incidence rate of the form:

In the above expression, is denoting the nonmonotonic force of infection at the rate , that is, is increasing when is small but it is decreasing when is large. The incidence rate can be used to interpret “psychological effects” by susceptibles, i.e., the infection force may get reduced as the number of infected increases for a large number of infected population because in such a situation the number of contacts may tend to reduce for the case . The psychological effects can be observed on the general public during the epidemic outbreak of SARS; aggressive measures and policies, such as border screening, mask-wearing, maintaining hygiene, etc., have been proven very effective (Xiao and Ruan 2007) in reducing the infective rate at the later stage of SARS outbreak, even when the number of infective individuals was getting increasing. In Eq. (2) when we consider then Eq. (2) reduces to the following form:where measures the infection force of the disease and describes the psychological effect from the behavioral change of susceptibles when the number of infectives is very high. The parameter measures the psychological effect of disease on the population when the infective individuals become sufficiently larger in society. The expression of Eq. (3) is a nonmonotonic functional which is also known as the simplified Monod–Haldane incidence rate (Kumar and Nilam 2021).

Whenever an epidemic outbreak emerges, the quarantine of infected people is an effective tool to reduce the spread of the epidemic. Hence, the role of quarantine people cannot be ignored in the spread control of the disease. In this study, to understand the impact of quarantine people, a quarantine class of individuals is introduced into the standard epidemic model for disease transmission dynamics. The word quarantine is simply meant to the situation of forced isolation or some restrictions on the movement of infected individuals to control the disease from spreading further. To reduce the new infection cases among susceptibles, one can quarantine those infectives who have tested positive. In the case of covid-19, quarantined people are those positively tested infected individuals who are staying in a separate room in the household or a repurposed hotel, depending on the severity of infection and risk factors for developing the severity of the disease as per the WHO report (World Health Organization 2021). For the present study, in case of an outbreak of infectious disease, infected individuals who are in hospitals as per the severity of symptoms are also considered as quarantined individuals. In 2002, Hethcote et al. studied the effect of quarantine class with a bilinear incidence rate for infectious disease in six epidemic models. Later, in 2017 Erdem et al., proposed the disease dynamics with the imperfect quarantine of individuals for altering the influenza infection in the population. In our study, it is supposed that quarantined individuals do not mix up with other individuals so they do not infect them, hence, they will not cause further infection in the population, and do help in reducing the incidence rate for new infection cases.

Motivated by the work of Kumar (2020); Michael et al. 1999; Erdem et al. 2017), the proposed work aims to study the effect of nonmonotonic incidence and quarantine class in disease transmission dynamics to achieve adequate understanding in implementing measures to prevent and control infectious diseases among people. For this purpose, a susceptible-infected-recovered-susceptible epidemic compartmental model is extended and proposed by introducing a quarantine class into it. In the model, it is assumed that recovered individuals are having temporary immunity, and with time they will lose the temporal immunity and will become susceptible again. Some infected individuals who are asymptomatic or do not seek medical attention or have not tested positive remain in the infected class for their entire infectious period. Further, it is assumed that some infected individuals move to the recovered class due to the auto-recovery or by adopting alternate therapies other than medical intervention while other infected individuals who have tested positive and seeking sufficient medical attention move into the quarantined class and will remain there until they are no longer infectious and move to recovered class. This paper also aims to provide a new Caputo derivative-based fractional-order susceptible-infected-quarantined-recovered-susceptible epidemic model along with a nonmonotonic incidence rate and its mathematically rich dynamics.

The rest of the paper is structured as follows: in Sect. 2, the fractional-order mathematical epidemic model is proposed, and some of its basic properties are presented. In Sect. 3, the model equilibria and their local stability behaviors are analyzed with the help of threshold value . The Global stability behaviors of the model equilibria are investigated using and fractional order Lyapunov function method (Boukhouima et al. 2020) in Sect. 4. In Sect. 5, numerical simulations are performed using Matlab 2012b to validate the analytical studies. Finally, Sect. 6 is dedicated to the discussion, and conclusion of the present study.

Fractional-Order Epidemic Model and Its Properties

This section is devoted to the formulation of a non-integer Caputo derivative-based mathematical epidemic model. For the development of the model, it is assumed that the total population size is . Further, the total population is divided into four subpopulations or classes, namely; susceptible S , infected , quarantine and recovered , respectively. The susceptible class consists of those individuals who can catch the infection by coming in close contact with infected people. The infected class consists of those individuals who have been caught by the infection and can transfer the infection to susceptibles under appropriate conditions. Quarantine class consists of those infected people who are tested positive and staying at their homes or hospitalized according to the level of severity of infection. More explanation about quarantine class has mentioned in Sect. 1. Recovered class consists of those individuals who have moved from and compartments due to auto-recovery and proper medical treatment, respectively. It is assumed that recovered individuals cannot have immunity for a long period, therefore, they can recatch the infection, and hence, they will become susceptible after a certain period. Let denote the constant recruitment rate of susceptible. Let be the transmission rate of disease among the susceptible population, be the psychological effect, be the rate at which recovered individuals move in compartment. Let be the natural death rate, and parameters denote the disease-induced mortality rates in and classes, respectively. The parameters , and denote the quarantine rate, the auto-recovery rate of infectives, and the recovery rate of quarantined people after medical attention, respectively. Under the above assumptions, the following fractional-order epidemic model is proposed:

Subject to the conditions:

Here, is taken in the form of the Caputo derivative of order such that .

In the proposed model (4), the term denotes the non-monotonic force of infection at the rate .

For biological reasons, all model parameters are supposed to be positive. The list of all variables and parameters of the model (4) is given below in Table 1.

Now, the following lemmas are being provided of the fractional calculus to derive the basic properties of the model (4).

Lemma 1

(Odibat and Shawagfeh 2007), Consider that and the functions and its fractional derivative . Then, we get that for all ,

, where .

It is observed that for above Lemma 1reduces to the classical mean value theorem.

Lemma 2

(Podlubny 1999), Consider and . Define , where denotes the two-parameter Mittag–Leffler function along with parameters and .Then the Laplace transformation of is given by

Lemma 3

(Podlubny 1999), Assume that is an arbitrary real number. If , then there exists a constant such that, for all in the complex plane,

Now, we state and prove the following theorems for the basic properties of the model (4) below:

Theorem 1

All the solutions of the proposed mathematical epidemic model (4) with initial conditions given by Eq. (5) are non-negative.

Proof

Consider that at . Firstly, we have to prove that . On contrary, assume that is not true. Then, there exists a such that

where is assumed sufficiently small. From the first equation of the model (4), we get, . From Lemma 1, for any , we have.consequently, we get , which is a contradiction of Eq. (6) that for . Therefore, we get that .

Next, we have to show that . Again, we shall complete this by the method of contradiction. Suppose that is not true. Then, there exists a such thatwhere is assumed sufficiently small. From the third equation of the proposed model (4), we get, > 0. From Lemma 1, for , we have.

As a result, we obtain , which is a contradiction of Eq. (7) that is for . Therefore, we get that .

Similarly, we can show that and for all . Therefore, as a consequence, all the solutions of the model (4) with initial conditions (5) are non-negative.

Theorem 2

The region is an attractive region of the model (4), where is a constant described in Lemma 3.

Adding all the equations of the model (4), we get the following fractional order differential equation: where . Now, as , we get

Now, take the initial value problem

From the comparison principle (Lu and Zhu 2018), the solution of Eqs. (9) and (10) satisfy the following inequality:

Take the Laplace Transform of Eq. (10), we get

Also, from Lemma 2, we get

Taking the inverse Laplace transform of the above, we get

Using Eq. (9), we have

Using Lemma 3, we obtain

where, is a constant, described in Lemma 3. So, as t , we obtain where .

Therefore, region is attracting all the solutions of the model (4).

Equilibria and Their Stability Analysis

This section is devoted to the existence of the model equilibria and their stability analysis. The model equilibria are obtained by setting all the equations of the model (4) to zero, which are as given below:

Disease-free equilibrium (DFE),.

Endemic equilibrium (EE), .

A detailed discussion on the existence of endemic equilibrium is given in section 3.3.

To investigate the stability behavior of both equilibria, we need to compute the basic reproduction number (Driessche and Watmough 2002). The basic reproduction number is obtained in Sect. 3.1 using the next-generation matrix method (Driessche and Watmough 2002; Ye and Xu 2019).

Computation of the Basic Reproduction Number :

For the computation of the basic reproduction number , first, we consider thatwhere , is the matrix of new infection term and is the matrix of transfer term into the compartments and out of the compartments. Jacobian matrices of and at are given by:

For the next-generation matrix, we need to calculate the inverse of matrix . The inverse of the matrix is given in "Appendix A".

Now, the next-generation matrix is:

On simplification, we get

From Ye and Xu (2019), the basic reproduction number is the spectral radius of the next generation matrix that is denoted by . Hence, the basic reproduction number of the model (4) is:

In the above Eq. (12), the term denotes the average life expectancy of infectious individuals. The right-hand side term of provides the number of secondary infections of susceptible individuals that one infected individual can produce in a disease-free population.

Further, the local stability behaviors of the equilibrium point of the model (4) are discussed with the help of the basic reproduction number .

Local Stability Analysis of Disease-Free Equilibrium ()

The local stability behavior of the disease-free equilibrium of the model (4) can be analyzed by the linearization of the model (4) at the disease-free equilibrium (. Therefore, a linearized matrix of the model (4) is obtained as given below:

From the above, the linearized matrix of the system (4) at DFE ()is given as:

The characteristic equation of the matrix is as follows:

All characteristic roots of Eq. (13) are given by ,, , .

We see that the eigenvalues and of the matrix are of with negative sign and the eigenvalue will be negative when . Hence, for all .

Now, when, the eigenvalue, and in this situation Therefore, the disease-free equilibrium ( is unstable when. Hence, by Kumar (2020) and Matignon 1996) and the above discussion, we present the following theorem for the local stability behavior of the disease-free equilibrium:

Theorem 3

The disease-free equilibrium of the proposed model (4) is:

locally asymptotically stable when the basic reproduction number .

unstable when .

Figure 2 is the graphical representation of Theorem 3 using data as given in Table 2. It confirms that disease will die out when and disease will continue its propagation among the population when .

Fig. 2

Basic reproduction number versus Infected population

Table 2

Parameters and their numerical value

Parameters	A	μ	β	σ	γ	d_1	ω	θ	d_2	δ
Value	9	0.03	0.003	0.002	0.0002	0.04	0.08	0.005	0.04	0.007

Transfer diagram of the model (4)

Existence of the Endemic Equilibrium ( and Its Local Stability Analysis

In this section, we present the existence of a unique endemic equilibrium and discuss its local stability behavior. The endemic equilibrium of the model is obtained by taking all equations of the model (4) to zero as given below:

On solving Eq. (14), we get

and is the solution of the following quadratic equation:

where the coefficients and of Eq. (15) are as given below:

Now, we state the result in form of the theorem for the existence of a unique positive equilibrium.

Theorem 4

There is at most one positive equilibrium of the model (4) if .

Proof

Suppose that

The coefficients and of are negative, and coefficient is positive when as all parameters are assumed to be positive. From Eq. (16), it can be seen that is a continuous function in . Therefore, by Descartes’ Rule of the signs (Wang 2004), there exists at most one real positive root of Eq. (16). Thus, after obtaining the unique positive value of , we can calculate the value of and , respectively. Hence, there exists a unique positive endemic equilibrium of the model (4).

Now, the local stability behavior of the endemic equilibrium is investigated.

To investigate the locally asymptotical stability of the endemic equilibrium , we need the following lemma:

Lemma 4

(Ahmed et al. 2006; Naik et al. 2020) Define the following characteristic equation:

The following conditions make all the roots of the characteristics Eq. (17) satisfy the equation given below:

For , the condition for Eq. (17) is .

For , the conditions for Eq. (17) are either Routh-Hurwitz conditions or ,

For , if the discriminant of , namely is positive, the following conditions are the necessary and sufficient conditions to satisfy Eq. (18):

If then condition Eq. (18) is satisfied for all .

For general is necessary for Eq. (18).

Now, to analyze the local stability behavior of , we linearize the model (4) at , and obtained the following Jacobian matrix:

The characteristic equation of the matrix is:where the coefficients and are given by:

All the coefficients , , , of the characteristic Eq. (19) are of positive sign. Thus, from condition (v) of Lemma 4, the unique positive endemic equilibrium is locally asymptotically stable (Ahmed et al. 2006; Naik et al. 2020). Hence, this result is stated in form of the theorem below:

Theorem 5

The endemic equilibrium of the fractional-order epidemic model (4) is locally asymptotically stable when and unstable otherwise.

Global Stability Analysis

In this section, we discuss the global stability behavior of the disease-free and endemic equilibria .

For global stability analysis, we assume that

Now, we establish the following theorems corresponding to the global stability behavior of and .

Theorem 6

The disease-free equilibrium of the model (4) is globally asymptotically stable when the basic reproduction number is less than or equal to unity, that is, .

From the system of Eq. (4), we derive the following conditions at

Now, consider the Lyapunov function as:

The differentiation of Eq. (21) along the solution of the model (4) is

Now from above, it can be observed that if + and then . In addition, the maximum invariant set in is the singleton set . Therefore, according to Lassalle’s invariance principle (Kumar 2020; Salle 1976), the disease-free equilibria of the model (4) is globally asymptotically stable when .

Theorem 7

When the basic reproduction number is strictly greater than one, that is, then the endemic equilibrium of the model (4) is globally asymptotically stable.

It is derived the following conditions from the system of Eqs. (4) at

Now, to prove the global stability of , we define the following Lyapunov function:

The differentiation of along the solution of system (4) is:

Using the property of arithmetic mean, we get: and if

Hence, when the inequalities (23)-(28) are satisfied simultaneously then . Besides, the maximum invariant set in is the singleton set . Therefore, according to Lassalle’s invariance principle (Kumar 2020; Salle 1976), the endemic equilibria of the model (4) is globally asymptotically stable when .

Numerical Results

In this section, we present the numerical simulation of the model (4) in support of our analytical findings. The numerical simulations are performed using the predictor-corrector method (Diethelm et al. 2002; MathsWorks 2012). For the simulation, we use the experimental data as given in Table 2.

From the data as given in Table 2, the value of the coefficients of Eq. (16) is These coefficients value satisfy Theorem 4 for the existence of a unique positive equilibrium. So, the unique positive endemic equilibrium is calculated as for which the basic reproduction number is calculated as The numerical values of the coefficients of Eq. (19) are evaluated using the data given in Table 2 and are given as . The numerical value of ,, and satisfy the condition (v) of Lemma 4. Thus, the unique positive endemic equilibrium is asymptotically stable whenever the basic reproduction number . This confirms Theorem 5 numerically.

Table 1

Notations and descriptions of the model’s variables and parameters

Notations	Description
St	Susceptible population
Yt	Infected population
Q(t)	Quarantine population
Z(t)	Recovered population
A	Recruitment rate
μ	Natural death rate
β	Transmission rate of infection
σ	Rate of psychological effects
γ	Rate at which recovered move to S class
d_1	Death due to disease in Y class
ω	Quarantine rate
θ	Auto recovery rate
d_2	Death due to disease in Q class
δ	Recovery rate due to treatment

The initial values of the subpopulations for numerical simulations are considered as given below:

Figures 3, 4, 5, and 6 depict the impact of the fractional derivative of order on the susceptible, infected, quarantined, and recovered subpopulations, respectively. On the ground of Theorem 5, the endemic equilibrium is stable for every value of fractional order which is shown by these figures. Figures 3, 4, 5 and 6 also confirm that the disease is endemic as all the solution curves tend to the endemic equilibrium . Further, from Figures 3, 4, 5 and 6, we observed that for , the system approaches the stationary state in a short duration. The stationary time is increased when we decrease the value of .

Fig. 3

Impact of fraction order on susceptible population

Fig. 4

Impact of fraction order on infected population

Fig. 5

Impact of fraction order on quarantine population

Fig. 6

Impact of fraction order on recovered population

Figure 3 depicts that for different values of as time passes, the number of susceptible individuals decreases, and finally the population reaches endemic equilibrium with time. It is noticed that as approaches one, the system reaches the steady-state in a very short period.

Figure 4 shows the effect of fractional order derivative on the infected individuals and it is observed that for the lower value of steady time is increasing. As approaches one the steady-state is achieved in a very short duration for infected populations which means infected populations will need more time to vanish when decreases due to the presence of the memory effect.

Figure 5 depicts the effect of fractional order derivative on the quarantined individuals and it is noticed that as the quarantine takes into account, for the lower value of steady time is increased. As approaches to the value one the steady-state is achieved in a very short time for the quarantined populations. It can be concluded from this figure that timely and effective quarantine can be proved a vital factor in the recovery of infectives and hence control the spread of infections.

Figure 6 shows that, if we decrease the fractional-order, the total number of recovered cases is decreased and it takes more time to approach to a steady-state. It is observed that as increases then the number of recovered individuals also increases and as reaches one steady state is achieved in a very short time.

Figure 7 depicts the infected individuals at various initial values of infected individuals (. It shows that when the number of infectives is low initially, infection occurs at a low level in comparison to high initial values of infected populations. Moreover, with any values of , the infected population approaches the same steady-state as time passes.

Fig. 7

Impact of different initial values on the infected populations at

Figure 8 shows the impact of timely quarantine of individuals on the infected population. It is depicted that as the quarantine of individuals takes into account in the system then the infected population increases initially but later due to the timely quarantine of individuals, the infection reaches its peaks and then significantly decreases. Hence, timely quarantine of individuals is crucial to control the peak of infection during the epidemic outbreak.

Fig. 8

Quarantine population versus infected population at

Figure 9 shows the impact of timely quarantine of individuals on the recovered population. It is obvious that as the quarantine takes place timely the recovery increases and reaches its steady state. In Fig. 10, it is observed that, initially, as the number of infected individuals increases recovered population is at a low level, and as time passes recovery increases with the increase of infected individuals and reaches its steady state.

Fig. 9

Quarantine population versus recovered population

Fig. 10

Infected population versus recovered population

Figures 11 and 12 are portrayed to observe the impact of psychological effects on both susceptible and infected populations. As the value of the rate of psychological effects increases, the number of susceptibles increases, and the peak of infected individuals start decreasing as shown in Fig. 11 and Fig. 12, respectively. Hence, the psychological effects can be used as an effective method for reducing the peak of infection in society during an outbreak.

Fig. 11

Impact of psychological effects on the susceptible population at

Fig. 12

Impact of psychological effects on infected individuals at

Discussion and Conclusion

This article aimed to propose and mathematically analyze a fractional-order epidemic model to understand the role of quarantine individuals along with a deep understanding of psychological effects on susceptible individuals in disease transmission dynamics during an outbreak. Therefore, a class of quarantined individuals is incorporated into the standard susceptible-infected- recovered-susceptible compartmental epidemic model with Monod-Haldane type incidence rate of infection. Monod-Haldane type incidence rate is a nonmonotonic incidence rate that considers the psychological effects on susceptibles during an epidemic. The Caputo derivative is considered for the population rate of each subpopulation. The Caputo derivative is a fractional-order derivative that considers the memory effects which can be seen for any emerging disease. Whenever a disease emerges in society, at the beginning of the disease, we always use the knowledge which we have collected from previously emerged diseases. This collected information is always helpful for controlling the transmission of newly emerging diseases and this is the memory effect in the fractional-order epidemic model that we cannot measure by the integer-order derivative-based epidemic model. The model analysis has shown that the proposed model is well-posed i.e. the solutions of the model are non-negative and bounded in a compact region; and the model has two equilibria, namely disease-free and endemic. Further, we computed the basic reproduction number by the Next-generation matrix method and investigated that the disease-free equilibrium is locally asymptotically stable when the basic reproduction number is less than the unity and unstable when it is greater than unity as stated in Theorem 3. The existence of a unique positive equilibrium is established and the local stability behavior of endemic equilibrium is analyzed. We investigated that EE is locally asymptotically stable when is greater than unity and follows any one condition of Lemma 4 as stated in Theorem 5. Furthermore, the global stability behaviors of both equilibria DFE and EE have been analyzed by the Lyapunov function method and it is shown that EE is globally asymptotically stable when and DFE is globally asymptotically stable when as stated in Theorems 6 and 7, respectively.

Numerical results are proposed in support of analytical results. Figure 2 is plotted in support of Theorem 3 which has shown that when then the disease-free equilibrium is locally asymptotically stable and when slightly increases from one then endemic equilibrium arrises. Figures 3, 4, 5, and 6 are plotted to show the impact of fractional order on susceptible, infected, quarantine, and recovers subpopulations, respectively. These figures depict that as the value of fractional order increases to one, each subpopulation takes less time to reach its steady state. We noticed that various values of have no effect on the stability nature of the equilibrium points but affect only the time to reach the equilibrium state. These graphs support Theorem 5 which states that the EE is locally asymptotically stable when greater than one. It is also observed from the simulation that any initial population of infectives does not affect the steady-state of the infected population as shown in Fig. 7. Further, numerical simulations show that the timely quarantine of infected individuals helps us in reducing the new infection cases significantly, and also helps us in increment of recovered population. From graphs, it is observed that if we timely monitor the psychological effects and take the appropriate action to increase the value of the rate of psychological effects that helps to reduce the new infection significantly in society.

In conclusion, this study suggests that if health policymakers and professionals consider the memory effect, do the timely quarantine of infected individuals, and focus on the appropriate actions to increase the psychological effects on susceptibles then we can predict the disease more correctly, may reduce the impact of the epidemic, and this can also be helpful to reduce the period of the outbreak as well.