key: cord-0765987-ec23tirk authors: Mandal, Manotosh; Jana, Soovoojeet; Nandi, Swapan Kumar; Kar, T. K. title: Modelling and control of a fractional-order epidemic model with fear effect date: 2020-09-26 journal: Energy Ecol Environ DOI: 10.1007/s40974-020-00192-0 sha: 08b6d8cc8d6d0b76fcb6935fbaba20707315c18f doc_id: 765987 cord_uid: ec23tirk In this paper, we formulate and study a new fractional-order SIS epidemic model with fear effect of an infectious disease and treatment control. The existence and uniqueness, nonnegativity and finiteness of the system solutions for the proposed model have been analysed. All equilibria of the model system are found, and their local and also global stability analyses are examined. Conditions for fractional backward and fractional Hopf bifurcation are also analysed. We study how the disease control parameter, level of fear and fractional order play a role in the stability of equilibria and Hopf bifurcation. Further, we have established our analytical results through several numerical simulations. Infectious diseases have become one of the most threatening issues in todays lifestyle. The infectious diseases including chickenpox, measles, cholera, tuberculosis, influenza, SARS, COVID-19, etc., have massive impact than the other types of noninfectious diseases as these types of diseases can be transmitted from one individual to another, some are spread out by bites from insects or animals, and some diseases are acquired by consuming flagitious water or food or being exposed to organisms in the environment. Hence, an infectious disease may spread in a huge region throughout the globe within a very short time period. Due to the improvement of lifestyle of common people and massive enhancement in transportation and globalization, an infectious diseases become pandemic in a less amount of time compared to earlier days (for example Spanish flu in twentieth century had taken a long period of time to become a pandemic compared to the pandemic due to COVID-19 in ongoing time period). Thus, not only to study the dynamics of an infectious disease, but also to control or determining the procedure to control an infectious disease, researchers from various fields have engaged themselves. Several mathematical tools induce a great affect in regulating many infectious diseases. Mathematical modelling is one of the popular and commonly used mathematical tools that can be utilized efficaciously to monitor several infectious diseases. The first effective research work on studying an infectious diseases using mathematical model was probably first done by Kermack and Mckendric (1927) . There are several research works where the authors have been studied mathematical models on controlling infectious diseases (see Zhou et al. 2014; Jana et al. 2016 Jana et al. , 2017a . In his book, Murray (2002) has analysed theoretical works on some simple SI, SIS, SIR epidemic models. Present days mathematicians formulate more complex epidemic model systems which are almost identical to the real-world problems. In mathematical studies of an epidemic model, proper control strategies are significant tools in monitoring and controlling the infectious diseases. In this regard, two most effective and commonly used control tools used in controlling diseases are vaccination and treatment (Kar and Jana 2013a, b) . However, isolation (see Jana et al. 2017b) , insecticide (to control vector-borne disease)(see Kar and Jana 2013b) , etc., are also used in recent days. In the present work, we consider fear factor due to the infectious disease which is an another important parameter for epidemic model. People generally got scared and make a significant distance to prevent the infectious disease. Hence, fear induced by the infectious disease compelled susceptible population to isolate which actually decreases the birth rate of population and also survival of adults is affected consequently. Due to the SARS outbreak in HongKong (SARS started on November 2002; peaked on March 2003; eliminated on June, 2003) , the birth rate dramatically had fallen from 8.742 (2002) to 8.436 (2003) and again it increased to 8.558 in 2004 (see Worldbank (2018) ). In our work, we formulate a new fractional-order SIS epidemic model with fear effect (Clinchy et al. 2013; Wang et al. 2016) and treatment control to eradicate the disease. Mathematical models based on ordinary differential equations depict the interactions between the population classes (e.g. susceptible-infected), and it is a classical approach in theoretical epidemiology. In recent days, researchers are concerned in developing mathematical model by the fractional-order differential equations because it is an important apparatus for the study of the memory and also some hereditary properties of several biological components and also it has a very close relation to the fractal theory. These are the main advantages of the fractional-order derivatives which are not included in the models based on the ordinary-order derivatives. Therefore, fractional-order derivatives are more naturalistic than the ordinary derivatives. The models constructed by the fractional-order derivatives have been widely applied in different field of research after the some famous books and research works on fractional-order differential equations (Diethelm and Braunschweig 2003; Hilfer 2000; Kilbas et al. 2006; Miller and Ross 1993; Petras 2011; Podlubny 1999; Sabatier et al. 2007; Sengupta et al. 2020; Yadav et al. 2020; Karthikeyan and Arul 2020) . However, a fractional-order derivative may be defined in several ways. The most popular and commonly useful definitions of fractional-order derivatives are in the sense of Riemann-Liouville, Grünwald-Letnikov and Caputo definitions (Petras 2011; Podlubny 1999) . Since, in the definition due to Caputo, the initial conditions can be expressed in a similar fashion as the integer-order differentiation, it is the most frequently used fractional-order derivative in mathematical modelling. There are very few theories to analyse the dynamical behaviour of the mathematical models with fractional-order derivatives (Delavari et al. 2012; Deshpande et al. 2017; Guo 2014; Li et al. 2016; Liang and Wu 2015; Li et al. 2010) . In recent days, the linearization, Lyapunov method and Lyapunov direct method (due to Lyapunov function Garrappa 2010) have established to study the local and global stability of equilibrium point(s) of a fractional-order system. Our research article studies the modelling and control of an epidemic problem using some advanced mathematical tools. In particular, this article focuses on approaches that integrate ecology and environment considerations and promote multidisciplinary solutions for tackling environmental problems and complex socio-economic problems. Present global scenarios show that human epidemiology has a significant impact on both the ecology and the environment. As this research article studies the modelling and control of an epidemic problem, the theoretical analysis of epidemiology is relevant to both ecology and environment (see Venturino et al. 2016; Maji et al. 2019; Khatua et al. 2020, etc.) . Moreover, the modelling approach using fractional-order differential equations gives an additional novelty and dimension of the article. The rest of the article is prepared as follows: in Sect. 2, we formulate an SIS-type fractional-order epidemic system and then discuss some preliminaries on fractional-order derivatives in Sect. 3. In Sect. 4, we describe the dynamics of the proposed model, including equilibria, their existence, bifurcation and both the local and the global behaviours. In Sect. 5, we present some computer simulation works to validate our theoretical observations. In Sect. 6, we provide some of the central findings of the paper and the last section is devoted to draw the conclusions of the paper. In this section, we construct a new fractional-order SIS (susceptible-infected-susceptible)-type epidemic model with treatment control and fear effect due to the infectious disease. We separate the total populations into two timedependable mutually exclusive classes, viz. susceptible S(t) (vulnerable for the infectious disease) and infected I(t) (already have pathogen of the disease and infectious for the population). Assume that birth rate of the susceptible population be r, natural death rate of both susceptible and infected population be a and d be the death rate due to the disease. Also, consider b as the death rate of susceptible population due to intraspecific competition (Ghirlanda et al. 2009 ). Note that b is smaller as compared to a, but it is not negligible for a large population. Further, let us take fear function as f ðb; IÞ which reduces the growth rate of susceptible population. Again, we consider that the disease transmission function is g(I) and the disease control function is h(u). Consider m as the natural recovery rate of the infected population from the disease. Finally, assume that the recovery from infectious disease is not permanent and hence all the recovered individuals will be reinfected, i.e. they are susceptible to the infectious disease. We formulate our proposed model in the following form by the above conditions: Here t 0 D a t ; a 2 ð0; 1 is the fractional-order derivative (aorder) in the sense of Caputo (Petras 2011; Podlubny 1999) . In our present model, we consider fear function as f ðb; IÞ ¼ 1 1þbI with the level of fear as b and the disease transmission function as gðIÞ ¼ cI 1þdI , where c is disease transmission rate and d is half-saturation incidence rate. Further, we assume that the disease control function is hðuÞ ¼ u. Here u is the treatment control. So our model becomes System (2.2) for different fractional differentiations a 1 ; a 2 2 ð0; 1 is known as incommensurate fractionalorder system. It is presented by following way System (2.2) is known as commensurate fractional-order system. Note that all the parameters considered in the formulations of the model are nonnegative. In this section, we recall some essential definitions, effective lemmas for both types of commensurate and incommensurate fractional-order systems. where C n ½t 0 ; 1Þ is a space of n times continuously differentiable functions on ½t 0 ; 1Þ, t [ t 0 and CðÁÞ be the Gamma function with n 2 Z þ (the set of all positive integers) such that n À 1 . Particularly, if 0\1, the definition becomes: Then we have the inequality: Here E a is the Mittag-Leffler function of one parameter. Lemma 3.4 (Petras 2011) Assume the following fractional-order system of the order a 2 ð0; 1 where x 2 R n and M 2 R nÂn (set of all n 9 n real matrices). The system is said to be stable asymptotically iff j argðkÞj [ ap 2 satisfy for each of the eigenvalues k of the matrix M and the system is said to be stable only iff j argðkÞj ! ap 2 for each of the eigenvalues of matrix M with the eigenvalues satisfying the critical condition j argðkÞj ¼ ap 2 must have geometric multiplicity one. Lemma 3.5 (Petras 2011) Let us assume the fractionalorder system of the order a 2 ð0; 1 A stationary point of the system is called to be locally asymptotically stable iff j argðk k Þj [ ap 2 for all eigenvalues k k ðk ¼ 1; 2. . .nÞ of Jacobian matrix J ¼ of ox calculated at the corresponding stationary point. Modelling and control of a fractional-order epidemic model with fear effect 4 The fractional-order system Here, first we assume commensurate fractional-order system (2.2) by considering a 1 ¼ a 2 ¼ a 2 ð0; 1 in the incommensurate system (2.3). Now we analyse the existence criteria of nonnegative equilibria of constructed system (2.2). The system possesses the equilibrium points which are trivial equilibrium (TE) E 0 ð0; 0Þ, disease-free equilibrium (DFE) E 1 ððr À aÞ=b; 0Þ and three endemic equilibria Þ depending upon the positive roots I à , I à 1 and I à 2 of the following cubic equation. where R 0 ¼ cr bðaþdþmþuÞþca is known as the basic reproduction number which can be defined as the number of new infections made by an infected individual to the susceptible population during its infection period. Also, We rewrite the cubic equation (4.1) in following standard form of strum's function : The strum's functions of equation (4.2) may be represented as We use the famous Strum's theorem to verify the existence of number of distinct positive roots of equation (4.2) and hence of equation (4.1). It can be easily verified that Here the number of changes of signs of the successive strum's functions is zero. Here the number of changes of signs of the strum's functions depends on the sign of B. If B, then the number of changes sign of the successive strum's functions is one, and hence equation (4.2) possesses exactly one positive root. If B [ 0, then the number of changes sign of the successive strum's functions is two, and hence equation (4.2) possesses exactly two positive roots. Since a 0 ; a 1 [ 0 for any choice of parameters, a 3 ; a 2 [ 0 and 2a 2 1 \9a 0 a 2 together imply that With this above discussion, we state the next theorem: From above Theorem 4.1, we can detect that system (2.2) possesses two endemic equilibrium points if R 0 \1, the basic reproduction number is less than unity, and this confirms that the system passes through a backward bifurcation at R 0 ¼ 1. Hence R 0 \1 is not a sufficient condition to eradicate the disease from our system (Fig. 1) . The backward bifurcation diagram shows that when R 0 \1, there exist two endemic equilibrium points and one of which is asymptotically stable and other is unstable. Theorem 4.2 (El-Saka 2015) Fractional-order system (2.2) undergoes backward bifurcation at R 0 ¼ 1 iff a 2 and a 3 ðbÞ ¼ 0. Here we choose b as backward bifurcation parameter. If the condition (ii) of Theorem 4.1 is true (the conditions for the existence of two EE points), then system (2.2) possesses backward bifurcation at R 0 ¼ 1. This implies that system (2.2) goes through two EE points in the interval ½R à 0 ; R 0 ¼ 1 and notice that one EE is asymptotically stable and other is unstable. Here (2015) implies that system (2.2) possesses a backward bifurcation at R 0 ¼ 1 and also the system has two EE for the parameter b in an interval. With the help of Note: A similar condition of Theorem 4.2 can be found if we consider a 3 as the function of disease control parameter, i.e. a 3 ðuÞ. In this stage, we state and prove the following lemma with the help of Li et al. (2010) . Lemma 4.3 Let us assume the following fractional-order differential equation of the order of differentiation a 2 ð0; 1: where f : ½t 0 ; 1Þ Â R ! R n ; R & R n be a function, if the Lipschitz condition is satisfied by f(t, x) with respect to x, then system (4.5) possesses unique solution on the interval ½t 0 ; 1Þ Â R. Proof We consider the region ½t 0 ; U  R to prove the existence and uniqueness criterion of the solutions of system (2.2), where R ¼ fðx; yÞ 2 R 2 : maxfjSj; jIjg Vg, and U and V are two finite positive real numbers. Let X ¼ ðS; IÞ and Y ¼ ðS 1 ; I 1 Þ be two points in R and set the mapping E : R ! R 2 by EðXÞ ¼ ðE 1 ðXÞ; E 2 ðXÞÞ where f g Therefore, the Lipschitz's condition is satisfied by the function E(X) with respect to the variable X ¼ ðS; IÞ 2 R. Lemma 4.3, implies that system (2.2) exhibits a unique solution X 2 R in respect of the initial conditions Accordingly we can state the next theorem: Theorem 4.4 For any initial point X t 0 ¼ ðS t 0 ; I t 0 Þ 2 R system (2.2) has a unique solution X t ¼ ðSðtÞ; IðtÞÞ 2 R for any time t [ t 0 . The nonnegative and boundedness solutions are the only useful solutions in mathematical epidemiology. Let us consider R þ ¼ fðS; IÞ 2 R : S; I 2 R þ g, where R þ is the set of all positive real numbers. Theorem 4.5 Each solution of system (2.2) starting in R þ 2 is nonnegative and uniformly bounded. Proof Nonnegativity: First consider that X t 0 ¼ ðS t 0 ; I t 0 Þ 2 R þ be an initial solution of system (2.2). Then we establish that any solution XðtÞ 2 R þ is nonnegative. Assume T be a real number such that t 0 t\T and The first equation of (2.2) implies that t 0 D a t SðtÞj SðTÞ ¼ 0. Lemma 3.2 implies that SðT þ Þ ¼ 0;, and this contradicts to the assumption SðT þ Þ. Hence for any t 2 ½t 0 ; 1Þ, we have SðtÞ ! 0. Similarly, we can also prove that IðtÞ ! 0 for all t 2 ½t 0 ; 1Þ. Uniform boundedness: To prove the uniform boundedness of the solutions, assume the succeeding function HðtÞ ¼ SðtÞ þ IðtÞ. Then, ; where Hence, the eigenvalues of system (2.2) at TE, E 0 ð0; 0Þ are k 1 ¼ Àða À rÞ and k 2 ¼ Àða þ d þ m þ uÞ. See that j argðk 1 Þj ¼ p [ ap 2 if a [ r, otherwise j argðk 1 Þj ¼ 0\ ap 2 , a 2 ð0; 1 and j argðk 2 Þj ¼ p [ ap 2 , where a 2 ð0; 1. Hence by Lemma 3.4, system (2.2) is locally stable around the trivial equilibrium asymptotically iff a [ r, i.e. if natural death rate is higher than the growth rate. Again the eigenvalues of system (2.2) at DFE, E 1 ððr À aÞ=b; 0Þ are k 1 ¼ Àðr À aÞ and k 2 ¼ crÀfcaþbðaþdþmþuÞg b . See that j argðk 1 Þj ¼ p [ ap 2 as r [ a(existence criterion for DFE), otherwise j argðk 1 Þj ¼ 0\ ap 2 , a 2 ð0; 1 and j argðk 2 Þj ¼ p [ ap 2 if crÀfcaþbðaþdþmþuÞg b , i.e. R 0 \1 , where a 2 ð0; 1. Hence by Lemma 3.4, the system is locally stable around the DFE asymptotically iff R 0 \1. Therefore we can state the following theorem relating the local stability of the TE and DFE points. (i) The TE, E 0 ð0; 0Þ of proposed model system (2.2) is locally asymptotically stable iff a [ r. (ii) The DFE, E 1 ððr À aÞ=b; 0Þ of model system (2.2) is locally asymptotically stable iff R 0 \1. Next our objective is to analyse the local stability behaviour of the unique EE point E à ðS à ; I à Þ. Characteristic polynomial of Jacobian matrix J at the EE, E à ðS à ; I Ã Þ point is given by Let k 1 and k 2 be the eigenvalues. Now, a 1 will be negative if Again, we consider Equation (4.6) as function of disease control parameter, u. Then we have Let k 1 and k 2 be the eigenvalues. Now, a 1 will be negative if The eigenvalues are The different values of k 1 and k 2 are depending on the coefficients a 1 and a 2 . (i) If a 2 1 ! a 2 and a 1 , then both the eigenvalues k 1 , k 2 are negative. Therefore j argðk 1;2 Þj ¼ p [ ap 2 , a 2 ð0; 1Þ and hence by Lemmas 3.4 and 3.5, the EE is asymptotically stable. (ii) If a 2 1 ! a 2 and a 1 ! 0, then one of the eigenvalues k 1 or k 2 will be nonnegative. Therefore j argðk i Þj ¼ 0\ ap 2 for i ¼ 1 or 2, where a 2 ð0; 1Þ and hence by Lemmas 3.4 and 3.5, the EE is unstable. (iii) If a 2 1  2 and a 1 [ 0, then both the eigenvalues k 1 and k 2 will be complex conjugate. is the imaginary unit: Therefore, j argðk 1;2 Þj ¼ tan À1 j ffiffiffiffiffiffiffiffiffi Thus we can state the next theorem regarding asymptotic local stability of the EE point. Theorem 4.7 The asymptotic local stability conditions for EE point E à ðS à ; I Ã Þ of system (2.2) are followed by (i) If a 2 1 ! a 2 and a 1 , then the EE is asymptotically stable. (ii) If a 2 1 ! a 2 and a 1 ! 0, then the EE is unstable. (iii) If a 2 1  2 and a 1 [ 0, then the EE is asymptotically stable. (iv) If a 2 1  2 and a 1 , then the EE is asymptotically stable. (v) If a 2 1  2 and a 1 ¼ 0, then the EE is asymptotically stable. Existence criteria of Hopf Bifurcation We can rewrite our proposed model system (2.2) as: where a 2 ð0; 1 and f ðb;xÞ ¼ f ðS; IÞ ¼ rS 1þbI À aS À bS 2 À cSI 1þdI þ ðm þ uÞI; cSI 1þdI À ðm þ u þ a þ dÞIÞ, is a function from ½t 0 ; 1Þ Â R to R 2 with R & R 2 with parameter b 2 R. For EE, E à point of system (4.10), suppose k 1 ðbÞ and k 2 ðbÞ are complex conjugate eigenvalues of the Jacobian matrix J of model system (4.10) at E à (which exists if a 2 1  2 and a 1 [ 0). System (4.10) goes through a Hopf bifurcation at critical value b ¼ b 0 [by equation (4.7)] such that following singularity condition[(i)] and transversality condition [(ii)] are satisfied (Deshpande et al. 2017) . Note: a similar condition of Hopf bifurcation can also be found if we consider system (2.2) as function of the disease control parameter, u. In this portion, first we only state some lemmas in order to study the global asymptotic stability criterion of DFE and the EE points. Let R be a closed and bounded subset in R and t : R ! R be a continuous function which satisfies the condition: t 0 D a t t 0. Then for any solution of the model system (4.5) starting from a point x 2 R & R n always remains in R. Also assume that P is the largest invariant subset of the set S ¼ fx 2 R : t 0 D a t t ¼ 0g. Then any solution x(t) of the model (4.5) starting in R converges to P as t ! 1. In peculiar case if P ¼ f0g, then xðtÞ ! 0 as t ! 1. Now we introduce and establish the theorems corresponding to the global asymptotic stability criterion for DFF and EE points. Theorem 4.10 If the basic reproduction number of system (2.2) is strictly less than unity, i.e. R 0 \1, then DFE point E 1 ððr À aÞ=b; 0Þ is asymptotically stable globally. Proof Using Lemma 4.9, we prove this theorem. We formulate the following positive definite Lyapunov function: Then with the help of a; the order, the constructed positive definite function V and the solution of system (2.2), we have So by Lemma 4.9, any solution starting in R converges to the largest invariant set S ¼ fðS; IÞ 2 R : t 0 D a t V ¼ 0g. Hence lim t!1 IðtÞ ¼ 0. Moreover, if IðtÞ ¼ 0, then system (2.2) reduces to single eqation t 0 D a t SðtÞ ¼ ðr À a À bSÞS ð4:11Þ The solution of (4.11) is ! SðtÞE a Àðr À aÞt a ½ As t ! 1, we have SðtÞ ¼ rÀa b . Therefore the equilibrium point of limit set (4.11) is asymptotically stable globally, and hence DFE point of system (2.2), E 1 is asymptotically stable globally. Hence the proof. h Theorem 4.11 The EE point, E à ðS à ; I Ã Þ of system (2.2) is asymptotically stable globally if it is asymptotically stable locally. Proof We construct the following Lyapunov function in order to prove this theorem: LðS; IÞ ¼U I À I à 1 þ ln It can be verified easily that the Lyapunov function L(S, I) is positive definite for all S [ 0; I [ 0. Since EE, E à ðS à ; I Ã Þ is a fixed point of constructed system (2.2), we can calculate the following: Thus by Lemma 4.8, it follows that M. Mandal et al. þ UcðS À S à ÞðI À I à Þ; by the Theorem4:5 Thus t 0 D a t LðS; IÞ 0 if 4aðaÀrÞ r 2 þ 2 [ 0. Now see that 4aðaÀrÞ r 2 þ 2 ¼ 2 r 2 ½2aða À rÞ þ r 2 ¼ 2 r 2 ½a 2 þ a 2 À 2arþ r 2 ¼ 2 r 2 ½a 2 þ ða À rÞ 2 [ 0 for any choice of parameters. Hence the proof. h Here, we assume incommensurate fractional-order system (2.3), a i 2 ð0; 1 for i ¼ 1; 2 and n ¼ 1 where M ¼ diagðk Aa 1 ; k Aa 2 Þ and J is the Jacobian matrix calculated at that fixed point E(S, I). To perform the numerical simulations of our model, we use modified predictor-corrector method (Diethelm and Braunschweig 2003) , and we mainly emphasize the effects of the parameters: level of fear b, treatment control u, the fractional-order a of system (2.2) and the fractional orders a 1 and a 2 of model system (2.3). First, we choose the parameters of our proposed model system (2.2) as r ¼ 0:85; a ¼ 0:1; b ¼ 0:1; b ¼ 0:12; c ¼ 0:1; m ¼ 0:01; d ¼ 0:045; d ¼ 0:6; u ¼ 0:06 and we see that backward bifurcation occurs at R 0 ¼ 1. Figure 1 depicts that if R 0 2 ðR à 0 ¼ 0:542; R 0 ¼ 1Þ, then there exist two EE points and one of them is stable and another is unstable which means that only the condition R 0 \1 is not sufficient to eliminate the disease from system. Again with the help of the parameters set r ¼ 0:66; a ¼ 0:67; b ¼ 0:6; b ¼ 0:1; c ¼ 0:51; m ¼ 0:1; d ¼ 0:1; d ¼ 0:23; u ¼ 0:01, we solve the system (2.2). Here as a [ r, Theorem 4.6 confirms that the TE, E 0 ð0; 0Þ is asymptotically stable. Fig. 5 shows that the EE(0.508, 0.034) is asymptotically stable (green curves). But if the a increases to the value a 0 ¼ 0:95, then we will see limit cycles are formed around EE for this particular order of differentiation and hence it proves the existence of Hopf bifurcation. If the a further increases, then limit cycles will be bigger. Figure 5 depicts these behaviours (red and blue curves). Choose the next parameter set: r ¼ 0:6; a ¼ 0:55; b ¼ 0:8; b ¼ 0:001; c ¼ 1; m ¼ 0:154; d ¼ 0:1; d ¼ 0:13; u ¼ 0:01; with fractional-order a ¼ 0:99. Figure 6 shows that the EE(0.819, 0.039) is asymptotically stable (blue curves). But if the b decreases to the value b 0 ¼ 0:3, then we will see limit cycles are formed around EE for this particular b ¼ b 0 and hence it proves the existence of Hopf bifurcation. Figure 6 shows this behaviour (red curves). Again, assume the parameter set: r ¼ 0:6; a ¼ 0:55; b ¼ 0:4; b ¼ 0:001; c ¼ 1; m ¼ 0:154; d ¼ 0:1; d ¼ 0:13; u ¼ 0:2; with fractional-order a ¼ 0:99, Fig. 7 shows that the EE(0.703, 0.043) is asymptotically stable (blue curves). But if the u decreases to the value u 0 ¼ 0:01, then we will see limit cycles are formed around EE for this particular u ¼ u 0 and hence it confirms that the system undergoes a Hopf bifurcation. Figure 7 shows these phenomena (red curves). Next, we choose the parameter set: r ¼ 0:5; a ¼ 0:5; b ¼ 0:5; b ¼ 0:01; c ¼ 0:9; m ¼ 0:01; d ¼ 0:1; d ¼ 0:13; u ¼ 0:01; with a 1 ¼ 0:99; a 2 ¼ 0:98, and using this, we solve system (2.3). Since for these parameters R 0 ¼ 0:986\1 the DFE(0.17,0) is asymptotically stable. In this present paper, we have proposed and analysed a new fractional-order SIS-type epidemic model with fear effect of the infectious disease and also the effect of treatment control. The fractional-order system is more naturalistic, and it brings on more affects as compared to the ordinary system or commonly known as integer-order system. The features of most of the real-world problems (or systems) assure that their model possibly more exactly formed with the fractional-order system. A higher-order system can also be modelled to a lower-order system by using fractional calculus. In case of control theory, the systems which are modelled by fractional calculus produce much better affect than the ordinary or integer-order system of differential equations. In our proposed model systems, i.e. commensurate system (2.2) and incommensurate system (2.3), we deduced the criterion for the existence and uniqueness of various equilibrium points of the systems in several fractional-order cases. Then we examined the conditions for nonnegativity of solutions and also uniform boundedness of that solutions. Furthermore, we established sufficient conditions to check that DFE, E 0 and the EE, E à of system (2.2) are asymptotically stable globally by developing suitable Lyapunov functions. Further, the theoretical results are asserting by some numerical simulations. We restrict our numerical works at the EE point of proposed systems (2.2) and (2.3). We have studied that R 0 , the basic reproduction number of our model acts a significant purpose in local and also global stability of DFE point (Fig. 3) . Also, treatment control u and the parameter of the level control b have key role to local and also the global asymptotic stability of EE point (Fig. 4) . We have demonstrated the criterion for the existence of Hopf bifurcation for commensurate system (2.2), against the parameter of the level of fear b and also noticed that similar types of conditions can be found for the existence of Hopf bifurcation against the disease control parameter u. We have also presented the Hopf bifurcation against the parameters a; b and u in Figs. 5, 6 and 7, respectively. In case of incommensurate system (2.3), Fig. 8 depicts the dynamical behaviour for the different fractional orders a 1 and a 2 . In this numerical simulation, we have applied the predictor-corrector PðECÞ m E (Predict, multi-term (Evaluate, Correct), Evaluate) method Diethelm et al. (2002) , Diethelm and Braunschweig (2003) and Garrappa (2010) with the help of MATLAB software. We use the implicit fractional linear multistep methods (FLMMs) as solver function for our proposed systems (2.2) and (2.3) with the fractional differential equations, and applying this method, we have described the phase portrait. Notice that PðECÞ m E is modified iteration formula of the method, PECE (Predict, Evaluate, Correct, Evaluate) and Adams-Moulton algorithm. We conclude that our proposed model yields efficient results and the model may serve as a tool for a wide range of utilizations in the mathematical epidemiology. In our future works, we develop and study several fractional-order epidemic models with the help of various numerical methods like Taylor series approximations method, homotopy perturbation method, Diethelma's method, Adomian decomposition method, etc. Also in future, we Modelling and control of a fractional-order epidemic model with fear effect stretch SIS-type epidemic model to SIR, SEIR, SEIRS model with treatment and vaccination controls in the presence of fear of the disease and apply fractional optimal control analysis of the disease control parameters which helps us to optimize the cost used in the controlling of the disease. Moreover, this mathematical model may be applied to examine the nature and also dynamics of almost every infectious disease. Predator-induced stress and the ecology of fear Stability analysis of Caputo fractional-order non linear system revisited On Hopf bifurcation in fractional dynamical systems Efficient solution of multi-term fractional differential equations using P(EC)mE methods A predictor-corrector approach for the numerical solution of fractional differential equations Backward bifurcations in fractional-order vaccination models Dynamic analysis of fractionalorder predator-prey biological economic system with Holling type II functional response On linear stability of predictor-corrector algorithms for fractional differential equations Intraspecific competition and coordination in the evolution of lateralization The Stability of Solutions for a Fractional Predator-Prey System Complex dynamics of an SIR epidemic model with saturated incidence rate and treatment Optimal control and stability analysis of an epidemic model with population dispersal Mathematical analysis of an epidemic model with isolation and optimal controls The effect of vaccines on backward bifurcation in a fractional-order HIV model A theoretical study on mathematical modelling of an infectious disease with application of optimal control Application of three controls optimally in a vector-borne disease-a mathematical study Uniqueness and stability results for non-local impulsive implicit hadamard fractional differential equations Contribution to the mathematical theory of epidemics Impact of human mobility on the transmission dynamics of infectious diseases Mittag-Leffler stability of fractional-order nonlinear dynamic systems Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag-Leffler stability Dynamical analysis of a fractional-order predator-prey model incorporating a prey refuge Laplace transform of fractional-order differential equations Bifurcation and global stability in an eco-epidemic model with refuge An introduction to the fractional calculus and fractional differential equations Fractional-order nonlinear systems: modeling analysis and simulation Advances in fractional calculus: theoretical developments and applications in physics and engineering Prediction of ventricular hypertrophy of heart using fractional calculus A model for the control of the mosaic virus disease in Jatropha curcas plantations Modelling the fear effect in predator-prey interactions Fertility rate, total (births per woman Stability analysis, control of simple chaotic system and its hybrid projective synchronization with fractional Lu system Optimal Vaccination Policies for an SIR Model with limited resources Acknowledgements The work of Soovoojeet Jana is financially supported by Dept of Science and Technology & Biotechnology, Govt. of West Bengal (vide memo no. 201 (Sanc.)/ST/P/S&T/16G-12/2018 dt 19-02-2019). Also the authors are very much grateful to the anonymous reviewers and Prof. Bin Chen, Editor in Chief of the journal, for their constructive comments and useful suggestions to improve both the quality and presentation of the manuscript significantly. Conflict of Interest On behalf of all authors, the corresponding author states that there is no conflict of interest.