arxiv: v1 [physics.soc-ph] 9 Mar 2017

Transcription

1 Memory effes on epidemi evoluion: The susepible-infeed-reovered epidemi model arxiv:73.39v [physis.so-ph] 9 Mar 27 M. Saeedian, M. Khalighi, N. Azimi-Tafreshi, 2 G. R. Jafari,,3,4 and M. Ausloos 5,6,7 Deparmen of Physis, Shahid Beheshi Universiy, G.C., Evin, Tehran 9839, Iran 2 Physis Deparmen, Insiue for Advaned Sudies in Basi Sienes, Zanjan, Iran 3 Shool of Biologial Sienes, Insiue for Researh in Fundamenal Sienes (IPM), Tehran, Iran 4 Cener for Nework Siene, Cenral European Universiy, H-5, Budapes, Hungary 5 GRAPES, rue de la Belle Jardinière 483, B-43, Angleur, Belgium 6 Shool of Managemen, Universiy of Leieser, Universiy Road, Leieser, LE 7RH, Unied Kingdom 7 ehumaniies group, Royal Neherlands Aademy of Ars and Sienes, Joan Muyskenweg 25, 96 CJ, Amserdam, The Neherlands (Daed: Marh, 27) Memory has a grea impa on he evoluion of every proess relaed o human soieies. Among hem, he evoluion of an epidemi is direly relaed o he individuals experienes. Indeed, any real epidemi proess is learly susained by a non-markovian dynamis: memory effes play an essenial role in he spreading of diseases. Inluding memory effes in he susepible-infeed-reovered (SIR) epidemi model seems very appropriae for suh an invesigaion. Thus, he memory prone SIR model dynamis is invesigaed using fraional derivaives. The deay of long-range memory, aken as a power-law funion, is direly onrolled by he order of he fraional derivaives in he orresponding nonlinear fraional differenial evoluion equaions. Here we assume fully mixed approximaion and show ha he epidemi hreshold is shifed o higher values han hose for he memoryless sysem, depending on his memory lengh deay exponen. We also onsider he SIR model on sruured neworks and sudy he effe of opology on hreshold poins in a non- Markovian dynamis. Furhermore, he lak of aess o he preise informaion abou he iniial ondiions or he pas evens plays a very relevan role in he orre esimaion or prediion of he epidemi evoluion. Suh a onsrain is analyzed and disussed. I. INTRODUCTION The sudy of epidemiology, onerning he dynamial evoluion of diseases wihin a populaion, has araed muh ineres during he reen years []. Mahemaial models of infeious diseases have been developed in order o inegrae realisi aspes of disease spreading [2 5]. A simple and ommonly sudied model, inrodued by Kermak and MKendrik, is he susepible-infeedreovered(sir) model[6]. In his model, populaions an be in eah of hree saes: susepible, infeed and reovered (removed), denoed by S, I and R, respeively. Originally, i is assumed ha susepible individuals beome infeed wih a rae proporional o he fraion of infeed individuals in he overall populaion(fully mixed approximaion) and infeed individuals reover a a onsan rae. The epidemi proess presens a (perolaion) ransiion beween a phase, in whih he disease oubreak reahes a finie fraion of he populaion, and a phase wih only a limied number of infeed individuals [7, 8]. The model has also been invesigaed for populaion on laies (e.g. [9 ]) or on neworks (e.g. [2 4]). For simpliiy, we will keep he medial epidemi voabulary hereafer. However, he model has also been ineresing for desribing nonmedial epidemis, suh as for finanial bubbles [5, 6], migraion [7], opinion formaion [8, 9] or inerne worm propagaion [2 22]. SIR models wih disribued delay and wih disree delay have also been sudied [23, 24]. In he usual SIR model, i is assumed ha all onas ransmi he disease wih he same probabiliy. Moreover, he ransmission and reovery oeffiiens are onsan. Hene he sae of sysem a eah ime does no depend on he previous hisory of he sysem: i is a memoryless, so-alled Markovian, proess. However, real surveys show evidene of a non-markovian spreading proess[25, 26] in agreemen wih ommon expeaion. The epidemi proesses evoluion and onrol, in human soieies, anno be onsidered wihou any memory effe. When a disease spreads wihin a human populaion, he experiene or knowledge of individuals abou ha disease should affe heir response [5]. If people know abou he hisory of a erain disease in he area where hey live, hey use differen preauions, suh as vainaion, when possible. Thus, some endogenous onrolled suppression of he spreading is expeed, alhough oher faors an help [27 29]. However, knowledge abou he hisory of a disease does no have he same influene a all imes. Experiene abou he prevalene of a disease and preauions relaed o he old imes are no always appliable or reommended, hene people end o follow new sraegies agains he diseases. In oher words, memory of he earlier imes ould have less effe on he presen siuaion, as ompared o more reen imes. I an be expeed ha long-range memory effes deay in ime more slowly han an exponenial deay, bu an ypially behave like a power-law damping funion. While muh effor has been made so far o deermine

2 2 exa epidemi hresholds in Markovian epidemi models [3 34], few works have been devoed o sudy he non- Markovian aspes of epidemi proesses [35, 36]. Furhermore, in his work we fous on long-range memory effes, whih means arbirarily long hisory an be inluded. Tha is in onras o shor-erm memory effes whih have been exensively sudied. For insane, Dodds and Was [37] inrodued a general model of onagion onsidering memory of pas exposures o a onagious influene. The auhors have argued ha heir model an fall ino one of hree universal lasses, due o he behavior of fixed poin urves. Also, in [38 4], he auhors onsider implii memory by applying asynhronous adaping in disease propagaion. They show ha his ype of memory an lead o a firs-order phase ransiion in oubreaks, hus hyseresis an arise in suh models [4]. I is here briefly realled ha fraional alulus is a valuable ool o observe he influene of memory effes on he dynamis of sysems [4 44], and has been reenly used in epidemiologial models [45 48]. Typially, he evoluion of epidemiologial models is desribed wih differenial equaions, he derivaives being of ineger order. By replaing he ordinary ime derivaive by a fraional derivaive, a ime orrelaion funion or memory kernel appears, hereby making he sae of he sysem dependen on all pas saes. Thus, i seems ha suh a mehod based on derivaives wih non-ineger order, as inrodued by Capuo for geophysis problems [49], is a very proper formalism for suh non-markovian problems. Moreover, Capuo s formalism provides he advanage ha i is no neessary o define he fraional order iniial ondiions, when solving suh differenial equaions [49 52]. Furhermore, he ime orrelaion funion, in he definiion of Capuo fraional derivaive, is a powerlaw funion, whih is flexible enough o refle he fa ha he onribuion of more early saes is noieably less relevan han he onribuion of more reen ones on he presen sae of he dynamial sysem. Mos of he previous works have sudied he epidemiologial models wih fraional order differenial equaions, from a mahemaial poin of view. They mainly foused on presening effeive a mahemaial mehods in order o solve he orresponding differenial equaions [47, 53 55]. For insane in [55] a mahemaial ool (he muli sep generalized differenial ransform mehod) is inrodued o approximae he numerial soluion of he SIR model wih fraional differenial equaions. Also in [46] he auhors use fraional order differenial equaions for epidemi models and onenrae on he equilibrium poins of he models and heir asympoi sabiliy of differenial equaions of fraional order. Oher variaions of he SIR model wih fraional derivaives have also been sudied. For insane, Seo, e al. inrodued he SIR epidemi model wih square roo ineraion of he susepible and infeed individuals and disussed he loal sabiliy analysis of he model [54]. Also in [47], numerial soluion of he SIR epidemi model of fraional order wih wo levels of infeion for he ransmission of viruses in a ompuer nework has been presened. In all previous works, he auhors rarely disuss he effe of fraional order differenial equaions and memory on he epidemi hresholds and he marosopi behavior of epidemi oubreaks. Hene, one quesion remains; we address i in his paper: How does he sysem robusness hange if memory is inluded in he SIR model? We also use he fraional differenial equaions, desribing he SIR model on sruured neworks, o see he effe of opology on he evoluion of he SIR model inluding memory effes. Furhermore he lak of aess o aurae informaion on iniial ondiions someimes leads o doub abou epidemi evoluion prediions [56]. The same ype of diffiuly ours in relaed problems, suh as in opinion formaion [57, 58]. Moreover, i may also happen in erain ases ha individuals do no believe in old sraegies in order o avoid he disease. This means ha he iniial ime for aking ino aoun he disease onrol memory is shifed oward more reen imes: hereafer, he dynamis is evolving wih a new fraion of susepible and infeed individuals, differen from ha predied by he soluion of he differenial equaions. In onras, he fraional alulus mehod allows us o hoose any arbirary iniial ime a whih he effe of iniial ondiions an be inrodued on he spreading dynamis wih a memory onen. The ineres of fraional alulus will appear hrough suh aspes in he ore of he paper. Thus, he paper is organized as follows. In Se. II, following Capuo s approah, we onver he differenial equaions of he sandard SIR model o he fraional derivaives, hereby allowing us o onsider memory effes. Using numerial analysis resuls (Se. III), we disuss he influene of memory on he epidemi hresholds in Se. IIIA. We also disuss he dynamis of a non- Markovian epidemi proess, when hoosing differen iniial ondiions or modifying he proporions of agens a a given ime in Se. IIIB. To omplee ourdisussion, we sudy he dynamis of he model on sruured neworks in Se. IV. We also poin ou ha we have observed qualiaively similar resuls for he SIS (susepible-infeedsusepible) epidemi model. The onlusions are found in Se. V. II. MEMORIAL PROCESS TO FRACTIONAL EQUATION The evoluion of he sandard SIR model is desribed by a se of oupled ordinary differenial equaions for susepible (S), infeed (I), and reovered(r) individuals,

3 3 respeively given by, ds() di() dr() = βs()i(), = βs()i() γi(), = γi(), () in whih, β and γ are infeion and reovery oeffiiens, respeively. The infeed individual makes β onas per uni ime produing new infeions wihin a mean infeious ime of order /γ. The evoluion of he model is onrolled by quaniy β/γ, suh ha above he epidemi hreshold, (β/γ), he disease spreads among a finie fraion of individuals. These (ordinary) differenial equaions desribe a Markov epidemi proess, in whih he sae of individuals a eah ime sep does no depend on previous seps. The se of Eqs. () an be solved ieraively unil ime. In pariular, he fraion of susepible individuals a ime, denoed as S, an be deermined. In fa, S is he size of oubreaks, i.e. he populaion ha has or has had he disease unil ime. In order o observe he influene of memory effes, firs we rewrie he differenial equaions () in erms of ime dependen inegrals as follows, ds() di() dr() = β κ( )S( )I( ), ( ) = κ( ) βs( )I( ) γi( ), = γ κ( )I( ), (2) in whih, κ( ) plays he role of a ime-dependen kernel and is equal o a dela funion δ( ) in a lassial Markov proess. In fa, any arbirary funion an be replaed by a sum of dela funions, hereby leading o a given ype of ime orrelaions. A proper hoie, in order o inlude long-erm memory effes, an be a power-law funion whih exhibis a slow deay suh ha he sae of he sysem a quie early imes also onribues o he evoluion of he sysem. This ype of kernel guaranees he exisene of saling feaures as i is ofen inrinsi in mos naural phenomena. Thus, le us onsider he following power-law orrelaion funion for κ( ): κ( ) = Γ(α ) ( ) α 2, (3) in whih < α and Γ(x) denoes he Gamma funion. The hoie of he oeffiien /Γ(α ) and exponen (α 2) allows us o rewrie Eqs. (2) o he form of fraional differenial equaions wih he Capuo-ype derivaive. If his kernel is subsiued ino Eqs. (2), he righ hand side of he equaions, by definiion are fraional inegrals of order (α ) on he inerval [,], denoedby D (α ). Applying a fraional Capuo derivaive oforderα onboh sides ofeaheq.(2), and using he fa he Capuo fraional derivaive and fraional inegral are inverse operaors, he following fraional differenial equaions an be obained for he SIR model: D α S() = βs()i(), D α I() = βs()i() γi(), D α R() = γi(), (4) where, D α denoes he Capuo derivaive of order α, defined for an arbirary funion y() as follows [49], D α y() = Γ(α ) y (τ)dτ ( ) α (5) Hene, he fraional derivaives, when inroduing a onvoluion inegral wih a power-law memory kernel, are useful o desribe memory effes in dynamial sysems. The deaying rae of he memory kernel (a imeorrelaion funion) depends on α. A lower value of α orresponds o more slowly-deaying ime-orrelaion funions (long memory). Hene, in some sense, he srengh (hrough he lengh ) of he memory is onrolled by α. As α, he influene of memory dereases: he sysem ends oward a memoryless sysem. Noe ha for simpliiy, we assume he same memory onribuions (same value of α) for differen saes of S, I and R. Obviously, more ompliaed funions han Eq. (3) and aking ino aoun differen α i (i =,2,3) ould be invesigaed in furher work o ake ino aoun differen ime sales. Alhough analyial soluions of Eqs. (4) are hard o obain for he general ase, hey an be obained a he early sage of he epidemi under a linearizaion approximaion. In his ase, i urns ou ha he number of infeed individuals behaves as a Miag-Leffler funion [5]: I() = E α,ζ () k ((β γ) α ) k Γ(αk +ζ) (6) in whih, ζ is a onsan, - whih depends on he iniial ondiions [5]. In pariular, for α = ζ =, he Miag- Leffler funion is he exponenial funion. Thus, in he early sage of epidemi dynamis, he growh rae of he infeed populaion in Eq. (6) is posiive, if β γ >. Therefore, he number of infeed individuals grows exponenially in suh a ase, for β > γ, as of ourse i is expeed for he sandard memoryless SIR model. The same reasoning applies in order o deermine he epidemi hreshold for α <. III. NUMERICAL RESULTS LeibereemphasizehaEqs.(4)onsisinasysem of oupled non-linear differenial equaions of fraional

4 4 S α=.4, = α=.5, = α=.6, = α=.7, = α=.8, = α=.9, = α=., = β/γ (β/γ).3.. =2 =2 =2 α FIG. : Oubreak size S for a SIR sysem having evolved unil ime =, vs. he parameer defining he hreshold: β/γ, when inluding memory effes. Eah urve orresponds o a differen value of α, as indiaed in he inse. As α dereases, he epidemi hreshold (β/γ) shifs o higher values. FIG. 2: Variaion of hreshold poin vs. α for differen finie imes = 2,2,2. For eah ime, he epidemi hreshold is shifed o higher values wih dereasing α. The axes are logarihmi and he numbers are presened as -base exponenial noaion. order, in he following general form D α y (i) () = f (i) (,y () (),y (2) (),y (3) ()) y (i) ( ) = y (i), (7) where, i =,2,3 and y (),y (2),y (3) denoe S,I,R ases respeively. Also, y (i) are onsans whih indiae he iniial ondiions. To solve he equaions, we use he he predior orreor algorihm, whih is well known for obaining a numerial soluion of firs order problems [59 6]. I is assumed ha here exis a unique soluion for eah ofy (i) on he inerval [,T] for a given se of iniial ondiions. Considering a uniform grid { n = nh : n =,,2,...,N}, in whih N is an ineger and h T/N, eah Eq. (7) an be rewrien in a disree form, n y n (i) = y(i) +h α b n k f (i) k, (8) k= where he oeffiiens b n k refer o he onribuion of eah of he n pas saes on he presen sae of n. The oeffiiens are given by b n k = (n k)α (n k) α Γ(α+) (9) Thereby afer solving Eq. (7), numerially, he influene of memory on he evoluion of he SIR epidemi model an be analyzed. As menioned in he Inroduion, le us onsider wo perinen aspes suessively: he finie ime behavior and he role of hanging iniial ondiions. Eq. (7) wih iniial ondiions y () = S = ǫ, y (2) = I = ǫ. Fig. shows he size of he oubreak for differen values of α, measured uil = and for ǫ = 4. The size of oubreak, S, is zero (wih auray 4 ) for small values of β/γ. The speifi value of β/γ, in whih he epidemi size sars o ge a non-zero value, is idenified as he epidemi hreshold poin. Thesaionaryimeforamemorylesssysem(α = )is =. Wih dereasing he value of α (inluding memory) he sysem needs muh ime o reah he saionary sae. Hene a =, he hreshold poin is shifed o he higher value of β/γ. Figure. 2 shows ha he hreshold poin is inreased wih dereasing of α for a finie ime. Furhermore, as i an be seen in Fig., he size of oubreaks dereases for dereasing α. Le he inerval [,] be he ime inerval in whih memory effes are aken ino aoun. In Fig. 3, we ompare he evoluion ofhemodelwih memoryfordifferenvaluesofhe finie ime. The memory effes are onsidered for a weigh α =. I is seen ha as ime evolves he influene of memory dereases, sine memory effes deay in ime like a power-law funion. Hene, he epidemi hreshold shifs o lower values of effeive infeion rae β/γ and approahes he hreshold of he memoryless model (α = ). The urves for α = a = 2 and = 2 are hardly disinguishable from he urve a = 2 and are no drawn for beer readabiliy. The variaion of hreshold poin, wih inreasing finie ime, is shown in Fig. 4. Furhermore, for a given β/γ value, i appears ha here is more ime available for disease spreading, whene more individuals beome infeed. A. Epidemi hreshold a finie imes Le us ompare he evoluion of a sysem inluding memory effes wih he memoryless ase. We solve B. Iniial ondiions Reall ha he dynamis of a non-markovian proess is direly influened by all evens from he beginning of he

5 5 S α=., =2 α=, =2 α=, =2 α=, =2 Susepible Populaion (a) β/γ FIG. 3: Order parameer S for a SIR sysem having evolved unil ime, when inluding muh memory (α = ). Eah urve orresponds o a differen finie ime, as indiaed in he inse. The hreshold values an be ompared wih ha of he orresponding epidemi hreshold for a memoryless sysem, i.e. when α = (and = 2). The urves for α = a = 2 and = 2 are no drawn for beer readabiliy. Infeed Populaion (b) (β/γ).3.. α= α=.5 α= Removed Populaion () FIG. 4: Variaion of hreshold poin vs. for differen values of α =,.5,.8. For eah α, he epidemi hreshold is shifed o lower values wih inreasing finie ime. The axes are logarihmi and he numbers are presened as base exponenial noaion. proess. However, some loss of informaion abou some period of ime in he pas may lead one o onsider ha he influene of memory migh no need o be onsidered as oninuous. I may happen, in many soial neworks, ha individuals do no have enough informaion abou he hisory of a disease, as reen ases and sudies indiae; e.g. see [62 65]. Only afer several individuals have already been infeed, do people sar o inrease heir knowledge abou he disease and ake differen preauions. The quesion arises on how he iniial ime a whih a non-markovian proess is sared, affes he subsequen dynamis of he proess. If wo Markovian proesses sar a wo differen imes, he evoluion of boh proesses is idenial. However, he senario is quie differen for a non-markovian pro FIG. 5: Effe of differen iniial imes on he dynamis of a non-markovian proess. The urves denoe he fraion of (a) susepible, (b) infeed and () removed individuals. Dashed and solid lines orrespond o Markovian and non-markovian proesses, respeively, sared from =. The urves wih symbols, orrespond o he dynamis of non-markovian proesses, sared from non-zero iniial imes wih differen iniial ondiions. ess, i.e. in whih he memory plays a role. This is illusraed hrough Fig. 5 where he fraions of susepible, infeed, and removed individuals are ompared in he ase of wo Markov and non-markov epidemi proesses. Coninuous and dashed (blak) lines orrespond o a sysem wih and wihou inluding memory effes, respeively, evolving from he same iniial ime =. As i an be seen he fraion of susepible individuals is greaer in a sysem wih inlusion of memory effes wih respe o ha ignoring he memory (Fig. 5(a)). In

6 6 oher words, he experiene and knowledge whih individuals have abou he disease are obviously helping hem o proe hemselves agains he disease. Equivalenly, in a sysem inluding he memory effes, he infeion grows more slowly as seen in Fig. 5(b). Thereafer, onsider ha a non-markovian proess, inluding memory effes, has evolved unil a speifi ime. Le he proess be oninuing is evoluion, bu le he memory of he sysem be removed a ha ime. This orresponds o having a new iniial ime and new iniial ondiions for he epidemis spreading. The proess an be oninued wihou or wih memory. The Markovian ase is rivial hereafer and hus no disussed. Insead, onsider ha memory effes are only aken ino aoun a his saring new iniial ime. In oher words, le he populaion ignore he disease onrol hisory (memory) unil ; le he sysem oninue is dynamis bu aking ino aoun memory effes hereafer from. The iniial ondiions for he evoluion of he sysem are now a fraion of susepible and infeed individuals a ime. The urves wih square symbols in Fig. 5, orrespond o wha happens for differen new iniial imes = 3,7, for he dynamis of suh a non-markovian epidemi proess. As i an be seen, a he beginning of he dynamis, he fraion of susepible individuals is redued, sine people do no know abou he disease. However, as soon as i is influened by memory, he sysem beomes more resilien o he spreading. Hene, he fraion of S individuals remains greaer as ompared o ha wih a memoryless sysem, having sared a =. In a similar manner, he fraion of infeed and removed individuals deviae from he original one and end oward he populaed saes of a memoryless sysem when he memory from furher pas imes is inluded. In his ase, he urves beome loser o he dashed urve orresponding o a memoryless sysem. Tha means ha he sysem loses he informaion relaed o pas imes and ends o presen a behavior similar o a memoryless sysem. Finally, one an onsider o remove he memory of an epidemi proess a various imes. A eah ime sep, he sysem is supposed o lose (or praially negae) he informaion abou he disease before some re-awareness ime (see also [66]) and o oninue is dynamis regardless of he pas. For illusraion, onsider he ase of suh a sudden awareness and is impa on epidemi oubreaks ase hrough Fig. 5; he sysem loses is memory a imes = 3 and 2 = 7, i.e. he dynamis is sopped a = 3, hen is oninued unil 2 = 7, removing all he hisory of he sysem before ha ime, nex reinroduing he memory dynamis again a 2 = 7: see he (red urves wih) riangular symbols in Fig. 5 orresponding o his ase of a double loss of memory. Noie ha in his pariular illusraive ase, he behavior of he sysem is seen o be lose o he dynamis of a memoryless sysem, sine onribuions of he memory of he sysem are someimes removed. Suh an illusraion poins o he ineres of he model in order o ompare i wih he ase of epidemis spreading waves [63]; for ompleeness, le i be poined ou ha he onneion of periodi epidemis o SIR models has been already menioned [67]: flu is yearly reurren. Noie also ha he value of α ould be modified a eah new awareness ime, bu his invesigaion goes ouside he presen paper. IV. THE MODEL ON STRUCTURED NETWORKS So far we have onsidered he fully mixed approximaion, suh ha an infeed individual is equally likely o spread he disease o any oher individual. However, in he real world an individual onnes o a small fraion of people. Hene, as is well known, more realisi modeling an be sudied hrough neworks, where heir opology has a signifian effe on he epidemi proess [68 7]. For homogeneous neworks, eah individual has he same number of onneions k k and disease propagaes wih spreading rae β k. In his ase, i is obvious ha he epidemi hreshold ( β γ ) is simply replaed by ( β γ ) k. I is also rue for he ase of fraional differenial Eqs. 4. In oher words, hreshold poin in Fig. for eah value of α is shifed o ( β γ ) k. Now, le us onsider heerogeneous sale free neworks wih degree disribuion P(k) k λ. In heerogeneous mean field approximaion, i is assumed ha all nodes are saisially equivalen and hus one an onsider groups of nodes wih he same degree k. Wih his assumpion he ordinary differenial equaions desribing he SIR model are given by where ds k () di k () dr k () Θ k () = = βks k ()Θ k (), = βks k ()Θ k () γi k (), = γi k (), () k kp(k)i k() k kp(k). () and i k,s k and r k denoe he densiy of infeed, susepible and removed nodes in eah group, respeively. I was urned ou ha in sale-free neworks haraerized by a degree exponen 2 < λ 3, here is no epidemi hreshold [7]. Following he same proedure, presened in Se.II, we an rewrie Eqs. o he fraional derivaives, as follows D α s k() = βs k ()Θ k (), D α i k () = βs k ()Θ k () γi k (), D α r k() = γi k (), (2)

7 7 Infeed Populaion.5..5 α=. α=.8 α=.5 α= FIG. 6: Fraion of infeed individuals versus ime for he SIRmodel onasale-free neworkwih degree exponenλ = 3 and for differen values of α. S x β/γ α=, = α=.4, = α=.6, = α=., =...3 β/γ FIG. 7: Oubreak size, S, for he SIR model on a salefree nework wih degree exponen λ = 3 in erms of β/γ. The dynamis is evolved unil ime =, when inluding memory effes. Eah urve orresponds o a differen value of α, as indiaed in he inser. V. CONCLUSION Memory plays a signifian role in he evoluion of many real dynamial proesses, inluding he ases of epidemi spreading. Here we have repored a sudy on he evoluion of he SIR epidemi model, onsidering memory effes. Using he fraional alulus ehnique, we show ha he dynamis of suh a sysem depends on he srengh of memory effes, onrolled by he order of fraional derivaives α. A finie imes, inluding memory effes, he epidemi hreshold (β/γ) is shifed o higher values han hose for memoryless sysems, a values depending on he memory deay rae α. In he ase ha he model evolves on heerogeneous sale-free neworks wih 2 < λ 3, he hreshold poin is always zero. However, he fraion of individuals who are infeed or reovered, is redued if he memory lengh inreases. Hene, memory renders he sysem more robus agains he disease spreading. If he epidemi proess evolves furher in ime, for a fixed memory srengh, (i) he disease an infe more individuals and (ii) he epidemi hreshold is shifed o smaller values and ends o he memoryless ase values. Furhermore, we have shown he following resul: he evoluion of an epidemi proess, inluding memory effes, muh depends on he fraion of infeed individuals a he beginning of he memory effe inserion in he evoluion. During a non-markovian epidemi proess, if he sysem abruply loses is memory a a definie ime andiffromha ime on, oneles he non-markovianproess oninue again, saring wih he number of infeed individuals a ha ime, he dynamis of he sysem deviaes from he basi ase, in whih he sysem oninuously inludes memory effes from he beginning of he proess. For a nework wih degree exponen λ = 3, we solve Eqs.(2) numerially. Figure 6 shows he evoluion of he fraion of oal infeed individuals i() = k P(k)i k(), wih onsidering memory effes wih differen values of α. While for a memoryless SIR model (α = ), he sysem reahesasaionarysae afer ashorime ( 2), he saionary ime is inreased wih dereasing he value of α. Furhermore, we obain he size of oubreaks a a finie ime. Figure 7 shows S, measured wih auray 5 unil = for differen values of α. As we an see he epidemi hreshold is always zero, as i is for Markov epidemi spreading on sale-free neworks wih λ = 3. However he size of epidemi dereases wih dereasing α. The same resuls are obained for neworks wih 2 < λ < 3. However for λ > 3, he epidemi hreshold is shifed wih inluding he memory, similar o wha is observed for he homogenous neworks. Our observaions are obained from a simple epidemiologial model: he SIR model. Obviously many parameers are here assumed o be onsan. We are aware ha some, e.g. poliy, feedbak migh influene he parameer values. They may depend on spae, groups, and ime. Exernal field ondiions may also surely influene real aspes. However, we guess ha many qualiaive behaviors as hose presened here are likely o be quie generally found in realiy. More advaned epidemi models, based on various ypes of omplex neworks are surely ineresing subjes for furher invesigaions, in line wih invesigaions suh as, e.g., in [3 34]. We also wish o poin ou ha we have observed qualiaively similar resuls for he SIS epidemi model. Finally, we may laim ha our resuls are no limied o he epidemiologial ( medial ) models bu also an be exended for analogous epidemi spreading of rumors, gossip, opinions, religions, and oher opis perinen o epidemis on many soial neworks.

8 8 Homogeneous Time Axes τ where Time Slowing Down T T (τ) = Γ(α+) {α (τ ) α } (A3) FIG. 8: Shemai omparison beween homogeneous and fraional ime axes. Appendix A I ould be insruive o sudy fraional order operaors wihin a geomeri inerpreaion (see ineresing referenes in [5]). Here we ompare he ime sales in fraional- and ineger-order dynamis. To image a geomeri inerpreaion, le us onsider he righ-sided fraional inegral of order α, x α = Γ(α) and wrie i in he form x α = v(τ)( τ) α dτ v(τ)dt (τ) (A) (A2) If we ompare Eq. A2 wih is ounerpar x α = v(τ)dτ, in he homogeneous ime sheme, he main differene is relaed o he differen ime variables, T and τ. Noie ha ime variable T (τ) has a saling propery. If we ake = k and τ = kτ, hen T (τ ) = k α T (τ). Hene, in he fraional order dynamis, he ime is aeleraing in he early ime and afer ha i is slowing down, as skehed in Fig. 8. In his ase, he passing ime in he wo axes of ime is no he same. For his reason, in epidemi fraional dynamis, he epidemi hreshold is shifed o he higher values. A lower α indiaes a sronger (long lasing) memory and a more pronouned shif of hreshold poin. However, if one wais long enough, he same behavior will be observed in boh fraional and usual homogeneous ime. In fraional dynamis, afer a very long ime, he hreshold poin oinides wih he one appearing in ineger-order dynamis. [] R. Pasor-Saorras, C. Casellano, P. Van Mieghem, and A. Vespignani, Rev. Mod. Phys. 87, 925 (25). [2] F. Brauer and C. Casillo-Chávez, Mahemaial Models in Populaion Biology and Epidemiology, Tex in Applied Mahemais 4 (Springer, New York, NY, USA, 22). [3] N. C. Grassly and Ch. Fraser, Na. Rev. Mirobiol. 6, 477 (28). [4] C. Ash, Siene 347, 23 (25). [5] N. K. Vianov and M. R. Ausloos, Knowledge epidemis and populaion dynamis models for desribing idea diffusion, in Models of Siene Dynamis: Enouners Beween Complexiy Theory and Informaion Sienes, A. Sharnhors, K. Boerner, and P. van den Besselaar, Eds. (Springer Verlag, Berlin Heidelberg, 22), Ch. 3, pp [6] W. O. Kermak and A. G. MKendrik, Pro. R. So. Lond. A 5, 7 (927). [7] P. Grassberger, Mah. Biosi. 63, 57 (983). [8] D. Sauffer and A. Aharony, Inroduion o Perolaion Theory 2nd ed. (Taylor and Franis, London, 994). [9] C.J. Rhodes and R. M. Anderson, J. Theor. Biol. 8, 25 (996). [] L.M. Sander, C.P. Warren, and I.M. Sokolov, Physia A 325, (23). [] T. Tomé and R.M. Ziff, Phys. Rev. E 82, 592 (2). [2] R.M. May and A.L. Lloyd, Phys. Rev. E 64, 662 (2). [3] E. Volz and L.A. Meyers, Pro. R. Soi. (London) B 274, 2925 (27). [4] R. Parshani, S. Carmi, and S. Havlin, Phys. Rev. Le. 4, 2587 (2). [5] G. Roundo, Logisi funion in large finanial rashes, in The Logisi Map and he Roue o Chaos: From he Beginning o Modern Appliaions, M. Ausloos and M. Dirikx, Eds., (Springer-Verlag, Berlin, Heidelberg, 25) pp [6] G. Roundo, Physia A 344, 77 (24). [7] N. K. Vianov, M. Ausloos, G. Roundo, Advanes in Complex Sysems, 5, Suppl. No., 2549 (22). [8] L. Zhao, H. Cui, X. Qiu, X. Wang, and J. Wang, Physia A 392, 995 (23). [9] S. Nizamani, N. Memon, and S. Galam, Physia A 46, 62 (24). [2] M. Liljensam, Y. Yuan, B. J. Premore, and D. Niol, A mixed absraion level simulaion model of large-sale Inerne worm infesaions, in Proeedings. h IEEE Inernaional Symposium on Modeling, Analysis and Simulaion of Compuer and Teleommuniaions Sysems MASCOTS 22, A. Boukerhe, S. K. Das, and S. Majumdar Eds., (IEEE, 22) pp [2] J. Kim, S. Radhakrishnan, and S. K. Dhall, Measuremen and analysis of worm propagaion on Inerne nework opology, in Proeedings. 3h Inernaional Conferene on Compuer Communiaions and Neworks, (ICCCN, 24) pp [22] B. K. Mishra and S. K. Pandey, Nonlinear Anal. Real World Appl., 4335 (2). [23] E. Berea and Y. Takeuhi, J. Mah. Biol. 33, 25 (995). [24] C. C. MCluskey, Nonlinear Anal. Real World Appl., 55 (2). [25] R. M. Yulmeyev, N. A. Emelyanova, S. A. Demin, F. M. Gafarov, P. Hänggi, and D. G. Yulmeyeva, Physia A 33, 3 (24).

9 9 [26] S. P. Blyhe and R. M. Anderson, Mah. Med. Biol. 5, 8 (988). [27] D. J. Eobihon, Biomed. Environ. Si. 3, 27 (99). [28] D.T. Halperin e al., Fuure HIV Therapy 2, 399 (28) [29] A. Legrève and E. Duveiller, Prevening poenial disease and pes epidemis under a hanging limae, h.4 in Climae hange and rop produion, M. P. Reynolds, Ed. (CABI Climae Change Series, Wallingford, 2) pp [3] Y. Moreno, R. Pasor-Saorras, and A. Vespignani, Eur. Phys. J. B 26, 52 (22). [3] M. A. Serrano and M. Boguñá, Phys. Rev. Le. 97, 887 (26). [32] M. Boguñá, R. Pasor-Saorras, and A. Vespignani, Phys. Rev. Le. 9, 287 (23). [33] M. E. J. Newman, Phys. Rev. E. 66, 628 (22). [34] W. Wang, Q. -H Liu, L. -F Zhong, M. Tang, H. Gao, and H. E. Sanley, Prediing he epidemi hreshold of he susepible-infeed-reovered model arxiv [35] P. V. Mieghem, and R. van de Bovenkamp, Phys. Rev. Le, 87 (23). [36] M. Boguñá, L. F. Lafuerza, R. Toral, and M. Ángeles Serrano, Phys. Rev. E 9, 428 (24). [37] P.S. Dodds, D. J. Was, Phys. Rev. Le. 92, 287 (24). [38] L. Chen, F. Ghanbarnejad, W. Cai, and P. Grassberger, Europhys. Le. 4, 5 (23). [39] W. Cai, L. Chen, F. Ghanbarnejad and P. Grassberger, NATURE Physis 25, 3457 (25). [4] L. Chen, F. Ghanbarnejad, D. Brokmann, arxiv: (26). [4] R. Herrmann, Fraional Calulus: An Inroduion for Physiiss, 2nd Ediion, World Sienifi, River Edge, NJ, 24. [42] P.L. Buzer, U. Wesphal, J. Douglas, W.R. Shneider, G. Zaslavsky, T. Nonnemaher, A. Blumen, and B. Wes, Appliaions of fraional alulus in physis, (World Sienifi, Singapore, 2). [43] R. Mezler and J. Klafer, Phys. Rep. 339, (2). [44] H. Safdari, M. Z. Kamali, A. H. Shirazi, M. Khaliqi, G. Jafari, and M. Ausloos, PLoS One. (5): e54983 (26). [45] E. F. D. Goufo, R. Mariz, and J. Munganga, Adv. Differenial Equaions. 278, (24). [46] H. A. A. El-Saka, Mah. Si. Le. 2, 95 (23). [47] A. A. M. Arafa, M. Khalil, and A. Hassan, J. Fra. Cal. Appl. 6, 28 (25). [48] A. A. M. Arafa, S. Z. Rida and M. Khalil, In. J. Biomah, 7, 4536 (24). [49] M. Capuo, Geophys. J. R. Asron. So. 3, 529 (967). [5] I. Podlubny, Fraional Differenial Equaions (Aademi Press, New York, 999). [5] I. Podlubny, Fraional Calulus and Applied Analysis, 5, 367 (22). [52] A. A. A. Kilbas, H. M. Srivasava, and J. J. Trujillo, Theory and appliaions of fraional differenial equaions 24 Norh-Holland Mahemais Sudies (Elsevier Siene B.V., Amserdam, 26). [53] F. Awawdeh, A. Adawi, and Z. Musafa, Chaos, Solions & Fraals 42, 347 (29). [54] Y.I. Seo, A. Zeb, G. Zaman, and I.H. Jung, Appl. Mah. 3, 882 (22). [55] A. A. Freiha and A. H. Handam, Appl. Appl. Mah (24). [56] A. H. Shirazi, A. Namaki, A. A. Roohi, and G. R. Jafari, Journal of Arifiial Soiey and Soial Simulaions 6, (23). [57] L. F. Caram, C. F. Caiafa, A. N. Proo, and M. Ausloos, Physia A 389, 2628 (2). [58] L. F. Caram, C. F. Caiafa, A. N. Proo, and M. Ausloos, Phys. Rev. E 92, 2285 (25). [59] K. Diehelm and A. D. Freed, The FraPECE subrouine for he numerial soluion of differenial equaions of fraional order, in Forshung und wissenshaflihes Rehnen 998 S. Heinzel and T. Plesser, Eds., 52 (GWDG- Berihe, Gesellshaf fr wissenshaflihe Daenverarbeiung, Gingen, 999), pp [6] K. Diehelm, N. J. Ford, and A. D. Freed, Numer. Algorihms 36, 3 (24). [6] R. Garrappa, In. J. Compu. Mah. 87, 228 (2). [62] L. Eihelberger, So. Si. Med (27). [63] B. K. Johns, Changing Waves: The Epidemis of 832 and 854, in Ch2olera: Hamilon s Forgoen Epidemis, D. A. Herring and H. T. Bales, Eds., (MMaser Univ., Hamilon, CND, 22) pp [64] D. M. Morens and J. K. Taubenberger, Lane Infe. Dis. 5, 852 (25). [65] O. Tomori, BMC Med. 3, 6 (25). [66] S. Funk, E. Gilad, C. Wakins, and V. A. A. Jansen, Pro. Nal. Aad. Si. USA 564 (29). [67] D. Greenhalgh, IMA J. Mah. Appl. Med. Biol. 5, 8 (988). [68] C. Moore and M. E. J. Newman, Phys. Rev. E 6, 5678 (2). [69] M. Kuperman and G. Abramson, Phys. Rev. Le. 86, 299 (2). [7] R. Pasor-Saorras and A. Vespignani, Phys. Rev. Le. 86, 32 (2).