Consequences of delays and imperfect implementation of isolation in epidemic control.

For centuries isolation has been the main control strategy of unforeseen epidemic outbreaks. When implemented in full and without delay, isolation is very effective. However, flawless implementation is seldom feasible in practice. We present an epidemic model called SIQ with an isolation protocol, focusing on the consequences of delays and incomplete identification of infected hosts. The continuum limit of this model is a system of Delay Differential Equations, the analysis of which reveals clearly the dependence of epidemic evolution on model parameters including disease reproductive number, isolation probability, speed of identification of infected hosts and recovery rates. Our model offers estimates on minimum response capabilities needed to curb outbreaks, and predictions of endemic states when containment fails. Critical response capability is expressed explicitly in terms of parameters that are easy to obtain, to assist in the evaluation of funding priorities involving preparedness and epidemics management.

Introduction

Roughly 400 infectious diseases have been identified since 1940. New pathogens are emerging at higher rates, despite the increase in awareness and vigilance. A grave public health concern is when and how the next outbreak will occur[1]; threats of imminent global outbreaks are real[2]. The 1917 Spanish influenza, which killed 50 million people, was the worst-ever pandemic on record – and that was back at a time when travel by ship was the fastest means of transportation around the globe. In today’s tightly connected world, an epidemic can potentially travel at jet-speed. Indeed, Swine flu was first detected in April of 2009 in Mexico and within a week it turned up in London.

In spite of advances in medicine, there is still no known vaccine for many infectious diseases, and even when a vaccine is known, it is generally impractical to maintain an adequate stockpile. Preventive measures may fall out of fashion when the disease is absent for long periods, only to hit harder when it returns due to the lower level of immunity. A standard recommendation of health organizations in the face of an outbreak is the isolation of infected individuals[3-5]. This control strategy dates back centuries and its usefulness has not diminished with time, as evidenced in the recent outbreaks of SARS in Asia and Ebola in West Africa. Research has progressed in the meantime; mathematical models to understand the impact and effectiveness of isolation and quarantine protocols have been proposed and analyzed[6-9]; it has been shown that the effectiveness of isolation protocols can be increased by contact tracing and isolation of individuals prior to symptoms[10,11].

Perfect implementation of isolation, however, can never be achieved. Among the obstacles are the following: Physical facilities to house the isolated population, and medical personnel to care for them, have to be available[11,12]. Even then, the identification of infected hosts, a crucial first step in any isolation strategy, can be a formidable task. At early stages of an outbreak, the general population – even local health authorities – may fail to recognize the symptoms of an infectious disease (which are often ambiguous), or they may not be cognizant of the potential dangers of runaway exponential growth of infectious diseases[9]. For reasons of their own, such as distrust of authorities, infected individuals may choose not to seek medical attention. The cost of preparedness and surveillance, and the maintenance of adequate infrastructure and personnel that can be deployed at a moment’s notice, can be prohibitively high[13].

A case in point was the initial less-than-adequate emergency response in the devastating Ebola epidemics in West Africa. It was not until Oct 2014 that billions of dollars were allocated for training of personnel, construction and maintenance of holding centers and treatment units. Once the response capability improved, and delays in identifying infected individuals were shortened, the spreading of the disease was more effectively curtailed[12,14]. Detrimental effects of delay were also reported during the SARS outbreak, which became more effectively controlled only when the time between onset of symptoms and isolation was shortened[12,15]. Though recent studies have used linear analysis[6,16] and simulations[17] to investigate the role of isolation, a general comprehensive theory of isolation as a means of epidemic control that focuses on consequences of imperfect implementation is lacking, as are studies of the nonlinear effects of isolation protocols.

This paper addresses fundamental issues pertaining to imperfections in the implementation of isolation strategies as a means of disease control. First and foremost is the minimum response capability required to contain an outbreak. We determine the isolation probability, the critical fraction of the infected population that must be isolated, and the identification time, how quickly local authorities must act to prevent a full-blown outbreak. Our model predicts the nature of the endemic state that follows should the response be inadequate, and the extent to which a stronger response can ameliorate the severity of an outbreak. It offers crucial information to health authorities such as estimates on the fraction of the population that can be expected to fall ill, and answers theoretical questions such as the potential benefits of longer durations of isolation.

These and other questions of a highly nonlinear nature are best tackled using dynamical systems techniques[18,19]. In the present work, we analyzed the continuum limit of an isolation strategy. Under some assumptions on identification and isolation times, we obtained a system of Delay Differential Equations (DDE) describing the time course of an infection following an outbreak. These equations are amenable to mathematical analysis[20,21], enabling us to understand rigorously many of the issues above. Our model is complex enough to capture key ideas, and simple enough to permit the representation of relevant quantities explicitly in terms of disease and control parameters. We view these as among our model’s greatest strengths.

Results

SIQ: an epidemic model with (partial) isolation

Consider, to begin with, a random network of N nodes; each node represents a host, and nodes that are linked by edges are neighbors. The mean degree at a node is given by 〈k〉, and an important parameter for us is the density of connections m = 〈k〉/N at each node. In the absence of any control mechanism, the situation is as follows. Each host is in one of four discrete states: healthy and susceptible (S), incubating (E), infectious (I), and recovered (R). Infectious hosts infect their neighbors at rate β; infected nodes incubate the disease for a period of σ units of time during which they are assumed to be neither symptomatic nor infectious. At the end of the incubation period they become infectious. Infectious nodes remain in that state until they recover and gain immunity for a period that on average lasts 1/α units of time at the end of which they rejoin the susceptible group.

To the setting above, we introduce the following isolation scheme: If a host remains infectious for τ units of time without having recovered, it enters a new state, Q (for quarantine), with probability p. The hosts that do not enter state Q at time τ remain infectious until they recover on their own. We use the terms isolation and quarantine interchangeably, and use the letter Q to better distinguish it from the infectious state I. A host in state Q does not infect others; it remains in this state for κ units of time, at the end of which it is discharged and joins the recovered class.

Figure 1 depicts the model states and main parameters: The five states of this model are S, E, I, Q and R, representing respectively the healthy and susceptible, incubating, infectious, isolated and recovered groups within the population. The numbers σ, τ, κ > 0 and are the new parameters of the model, with σ being the incubation period, τ representing the length of time between a host’s becoming infectious and entering isolation if it is identified, and κ the isolation time. The number p is the probability that an infectious host is identified and isolated. We refer to p as the isolation probability, and τ as the identification time.

Figure 1

The SIQ model. Healthy individuals in contact with an infected node become infected at rate β. Once infected, a node incubates the disease for σ units of time, during which it is not symptomatic and cannot transmit the disease. At the end of incubation, these nodes become infectious. Infectious nodes recover on their own at rate γ gaining immunity against the disease. Immunity is lost at a rate α. Those infectious nodes that are not recovered after τ units of time are, with probability p, identified and isolated. They remain in isolation for a duration of κ at the end of which they are released and rejoin the healthy group.

Table 1 summarizes the main parameters of the SIQ model and their meaning.

Table 1

Main parameters of the SIQ model.

Parameter	Meaning
β	transmission rate
γ	recovery rate
α	rate of immunity loss
m	density of contacts
p	probability to identify and isolate an infected individual
τ	time elapsed between infection and identification
σ	incubation period
κ	time spent in isolation (quarantine)

See Table 2 for specific values of parameters for given diseases.

Table 2

Critical response capability p and critical identification time τ (in days) for various diseases with basic reproductive number r, and σ as well as 1/γ (in days).

	r	σ	1/γ	p c	τ c
H1N1 2016 [Brazil]	1.7	4	7.0	0.41	4.7
Ebola 2014 [Guin./Lib.][28]	1.5	20	12.0	0.33	10.5
Ebola 2014 [Sierra Leone][29]	2.5	20	12.0	0.6	3.5
Spanish Flu 1917[29]	2	4	7.0	0.5	3.3
Influenza A[29]	1.54	0.23	3.0	0.35	1.0
Hepatitis A[11]	2.25	29.1	13.4	0.56	4.89
SARS[11]	2.90	11.8	21.6	0.66	4.31
Pertussis[11]	4.75	7.00	68.5	0.79	0.91
Smallpox[11]	4.75	11.8	17.0	0.79	0.26

The critical τ is calculated assuming that 80% of infected individuals are identified and isolated.

Continuum limit

We assume the total population is constant, and let S(t), E(t), I(t), Q(t) and R(t) denote the fractions of susceptible, incubating, infectious, isolated and recovered nodes at time t, so that S(t) + E(t) + I(t) + Q(t) + R(t) = 1.

Assuming that links between the infectious and susceptible nodes are uncorrelated, we obtain, by moment closure, the following system of DDEs in the continuum limit as :

Each term in the rate equation model (1)–(5) can be attributed to a specific transition between the compartments, see Fig. 1. For instance, the terms βmS(t)I(t) in Eqs (1) and (2) correspond to the transition from S (susceptible) to E (incubating). It is proportional to the product S(t)I(t) that is, the contact probability of an infected and a susceptible node, the transmission rate β, and the density of contacts m. After time σ, these incubating nodes undergo transition from E (incubating) to I (infectious), which is reflected in the terms βmS(t − σ)I(t − σ). The infectious nodes recover with the rate γ, therefore, out of the βmS(t − σ − τ)I(t − σ − τ) nodes entering I at time t − τ, only the fraction e−βmS(t − σ − τ)I(t − σ − τ) remains infectious at time t. These nodes are then undergo transition from I (infectious) to Q (isolated) with probability p corresponding to the terms βmpe−S(t − σ − τ)I(t − σ − τ) in Eqs (3) and (4). The derivation of system (1)–(5) is analogous to that in[22] where we provided a detailed derivation for a simpler model. The results below pertain to the continuum limit system (1)–(5).

Disease-free states and minimum response needed to squash a small outbreak

Without isolation and assuming everyone is susceptible, whether or not an outbreak will spread is determined by the disease reproductive numberthat is, an infection spreads if and only if r > 1[23]. We are interested in the case r > 1, so that if no measures are taken the infection will spread and reach an endemic state.

Our first result is that when the isolation probability is p and identification time is τ, the effective disease reproductive number r becomes

The number ε can be thought of as the response capability of the affected community. Upon the sudden appearance of a small group of infected individuals, the incipient epidemic is squashed if and only if r(1 − ε) < 1, i.e., the isolation protocol as described above reduces the disease reproductive number by rε. This is under the assumption that the outbreak starts near the disease-free equilibrium . The case where a fraction of individuals is immune to begin with will be discussed later.

We now consider separately the two components in the notion of response capability ε.

Minimum isolation probability

An outbreak can be prevented only if p > p where

The number p is a theoretical minimum: to have a chance to quell the outbreak, health authorities must have the ability to identify – and the facilities to isolate – a fraction p of the infectious population immediately after an individual becomes infectious.

Critical identification time

Possessing the ability to isolate individuals with probability p > p alone is not enough; one must be able to identify infectious hosts with sufficient speed. We prove that for each p > p, there is a critical identification timeso that the infection dies out if τ < τ.

Observe that in terms of what it takes to squash an outbreak, neither σ, the incubation period, nor α, the rate of recovery, play a role. The only requirement on κ is that it should be sufficiently long for the individual in isolation to recover. It is strictly the capability to quickly identify a large enough fraction of the infectious and to isolate it from the general population.

Figure 2 shows the results of simulations based on the set of delay differential equations above. For illustration, we consider r = 2.5 and 1/γ = 12 as in the 2014 Ebola outbreak in Sierra Leone (see Table 2). We also consider full identification probability p = 1. Then, starting from a small number of infected nodes, one observes two distinct asymptotic regimes: the system tends to a disease-free equilibrium for as we have shown; otherwise it tends to an endemic state. In the left panel, we considered an incubation time σ = 2.4 days and an isolation time κ = 60 days, in the middle panel an isolation time κ = 60 days and α = 1/1200. Finally, in the right panel, α = 1/1200 and σ = 2.4. The demarkation between the two asymptotic regimes is consistent with our analytical result that critical identification time is independent of incubation time σ, isolation time κ and the rate of immunity acquisition α. Observe though that for τ > τ the asymptotic states do depend on these values.

Figure 2

Critical identification time τ is independent incubation σ, isolation time κ and rate of immunity loss α. Colors represent the asymptotic equilibrium value of infectious individuals Ieq. We consider r = 2.5 and 1/γ = 12 (as for the 2014 Ebola outbreak in Sierra Leone) with full identification probability p = 1. Our predicted critical time is days independent of κ, σ and α, in agreement with simulations. For τ > τ, the disease-free stable state gives place to an endemic state. Left panel: α versus τ corresponds to the fixed values σ = 2.4 days and κ = 60. Center panel: σ versus τ corresponds to the fixed values κ = 60 days and α = 1/1200. Right panel: κ versus τ corresponds to the fixed values σ = 2.4 days and α = 1/1200.

Predicting critical isolation probabilities and identification times from data

The critical isolation probability p and identification time τ = τ(p, r) are readily computed once the disease reproductive number r is known. Table 2 shows p and τ for p = 0.8, that is, τ = τ(0.8, r), for several diseases.

As is evident from the data shown, even with the ability to identify 80% of infectious individuals, which requires considerable resources, the critical identification time τ can be as short as 3 days for severe outbreaks such as the Spanish Flu and the Ebola in Sierra Leone. That is to say, to prevent such an epidemic from spiraling out of control, infected individuals must be identified and isolated within 3 days of the end of their incubation period, which is typically when symptoms begin to appear and not always easy to recognize.

While the critical isolation probability p depends solely on the basic reproductive number r, the critical delay depends on both r and γ, so diseases having the same r such as Pertussis and Smallpox can have quite distinct τ. Recall that r = βm/γ. Among diseases with the same r, for the same value of p > p, our formula for τ implies that ; that is to say, the delay one can afford is inversely proportional to the transmission rate β of the disease for m fixed, and inversely proportional to the density of links m for fixed β.

Immunity at the time of an outbreak. The discussion above assumed that at time 0, when the infection first appears, one is near the disease-free equilibrium ; in particular, everyone is susceptible. We now consider the case where a fraction of the population has immunity at time 0, i.e., R(0) > 0, possibly from a previous exposure or from vaccination. Assuming , and that the dynamics that determine whether or not the outbreak can be contained is relatively fast, the situation can be approximated by one in which for the time period under consideration. Under these conditions, our results for the case R(0) = 0 together with a rescaling of the total population that participate in the dynamics yields the result that the outbreak is contained if and only if

That is, the effective disease reproductive number is scaled down by a factor of 1 − R(0) as expected. Having a fraction of the population with immunity allows greater leeway in terms of minimum isolation probability and critical delay can be computed explicitly as before.

Of major concern is what happens if the response capability is inadequate, that is, if the isolation probability falls below the required minimum, or if the required identification time is not met. Our next result addresses this scenario.

Prediction of endemic states

When the response capability as defined above is inadequate, the infection will persist. That is, for parameters p, τ, κ and σ such that , starting from the sudden appearance of a small group of infectious individuals, i.e., from an initial condition near , S = 1 − I and E = Q = R = 0, the I-component of the solution will grow; see Methods for more detail. For such initial conditions, we determine that the solution tends to an endemic equilibrium

A precise statement is as follows: Given p, τ, κ, σ and an initial condition as above, we prove analytically that there is a unique endemic equilibrium to which the solution can converge, and that this endemic equilibrium tends to (Seq, Eeq, Ieq, Qeq, Req) given by the formulas above as the initial condition tends to (1, 0, 0, 0, 0). The fact that orbits starting from near (1, 0, 0, 0, 0) will approach an endemic equilibrium is difficult to prove; this part of our prediction is supported by numerical simulations. Figure 3 shows the behavior of such solutions tending to their respective endemic equilibrium, with I > 0.

Figure 3

Transient dynamics. (a) Time evolution of epidemic from outbreak to endemic state when isolation protocol fails. (b) Zoom-in of (a) to show events in the first two cycles. These two panels show the solution starting from initial condition (S(t), E(t), I(t), Q(t), R(t)) = (1, 0, 0, 0, 0) for t < 0 and (S(0), E(0), I(0), Q(0), R(0)) = (0.99, 0, 0.01, 0, 0) corresponding to a small initial infection, with parameters r = 2.5, α = 1/1200, γ = 1/12, σ = 2.4, κ = 60, p = 1 and τ = 12. Panels (c–e) show the dependence of Ipeak, Qpeak (solid lines) and endemic equilibrium values Ieq, Qeq (dashed lines) on α, σ, and ε respectively. Unspecified parameters are as in (a and b).

Observe that ε, the response capability, not only determines whether or not a small outbreak can be contained, it influences strongly the eventual endemic state when the outbreak happens: Seq, the fraction of the population that remains healthy in the eventual endemic steady state, is inversely proportional to (1 − ε), and the fraction Ieq can be rewritten as

Observe in the numerators of the quantities above that r = 1, equivalently , is the tipping point on one side of which the initial outbreak is contained and on the other the system tends to an endemic equilibrium with Ieq > 0. While the parameters σ, κ and α appear in the formula for Ieq, hence in Eeq, Qeq and Req as can be seen in Fig. 2, we observe that the dependence on σ and κ is weak when , that is, when acquired immunity lasts for a significantly longer time period than incubation, recovery or isolation, which is the case for most diseases. Assuming that, we have

Here the effect of ε is particularly clear: Since is the number ε must beat in order to contain the initial outbreak, it can be seen as the minimal response capability required to prevent a full-blown epidemic. We therefore have the interpretation that the fraction of infectious in the endemic state is proportional to how far short the actual response capability ε is from the minimum response needed.

Infection cycles between outbreak and endemic equilibrium

Here we consider the situation where an outbreak occurs at a time no one in the population has immunity, such as after the disease has vanished for a couple of generations and suddenly returns, and the response is inadequate for curbing the initial outbreak. Figure 3(a,b) depict the time evolution of the epidemic from outbreak to near-endemic equilibrium assuming p, τ, and κ remain fixed throughout. It is evident that the infection occurs in waves, or cycles, at least initially, with the size of the infection decreasing with each cycle leveling off as it approaches endemic state. It is also quite evident that even as Ieq is small in the endemic state, the first wave of infection is orders of magnitude worse.

The oscillatory behavior seen is known to occur in the SIRS model[23] here it occurs for the same reason: as a large enough fraction of the population acquires immunity, and immunity is lost at a slow rate, the size of the susceptible fraction of the population decreases causing the rate of infection, which is proportional to S(t)I(t), to slow down – only to flare up again when the R-part of the population rejoins the S-part. What is different here is that we have assumed a constant response ε, so that for small I(t), the situation is similar to that discussed in the paragraph on R(0) > 0 in the first part of Results. In particular, if S(t) is such that r(1 − ε)S(t) < 1, the infection will temporarily abate, to flare up again when S(t) is large enough to bring this quantity above 1. Notice that the x-axis in Fig. 3(a) is in units t/(τ + σ + κ), so that the timescales for the cycles of infection are significantly longer than any of these quantities. Rather it is indexed to α and is affected by the fraction of individuals infected during the course of one cycle. A smaller group is infected in the second cycle than the first because a fraction of the population is immune at the time of the second outbreak.

We introduce the quantity Ipeak = maxI(t), the maximal fraction of infected individuals that is reached during the transient to Ieq; Qpeak is defined analogously. These quantities are of great relevance from the practical standpoint for health authorities involved in the planning of facilities and care for patients during isolation. To get our hands on these quantities, we consider first the limiting case where α = σ = ε = 0. In this case, our model reduces to the classical SIR model, for which we obtainso that can never be larger than , although a direct connection between and Ipeak cannot be made. Figure 3(c–e) describe how Ipeak and Qpeak vary as we permit α, σ and ε to become positive. From Fig. 3(c), we see that while α has a strong effect on Ieq and Qeq, its effect on Ipeak and Qpeak are relatively small, because immunity loss occurs on a much longer timescale than the time it takes to reach Ipeak or Qpeak. Figure 3(d) shows that both Ipeak and Qpeak are smaller for diseases with longer incubation periods. This presumably is due to the fact that a long incubation slows down the rate of infection, while recovery occurs at the same rate. Figure 3(e) shows that Ipeak decreases as the response is stronger as expected. The dependence of Qpeak on ε is more complicated: In the figure, τ is fixed and p is increased, i.e., a larger fraction of the infectious is isolated and Qpeak increases with ε for this interval of ε, but this trend will be reversed as ε continues to increase and Ipeak continues to drop.

Discussion

Isolation of infected hosts disconnects the infection pathway for diseases that are transmitted human-to-human and should, in principle, be 100% effective, yet it has not always stopped incipient outbreaks from burgeoning into full-blown epidemics. In this paper we called attention to the fact that a major reason may be “too little, too late”, meaning the identification of infected hosts may not be sufficiently thorough and are too often subject to delays. We sought a mathematical model that could offer guidance on the response needed to contain an outbreak, as well as predict the consequences when one falls short. Our aim was to develop a general theory that is relevant rather than to carry out detailed studies of specific diseases. Such a model is necessarily a little idealized to be amenable to (rigorous) mathematical analysis. With these goals in mind we came up with the system of delay differential equations (1)–(5), the analysis of which produced the results reported here.

We have framed response capability in terms of isolation probability, representing the fraction of the infectious population that will be identified and isolated, and identification time, the duration an infectious individual can transmit the disease before being isolated. An isolation strategy can prevent an outbreak only if the response capability exceeds a minimum threshold value. This was evident in the Ebola outbreak in West Africa, where the epidemic had started to escalate but reversed course after response capability was increased through international financial support[12].

We have sought to build a theory from first principles. To stress the novelty of this class of models, which focuses on isolation including the costs of failure in its implementation, we have taken the liberty to name the class of models considered here SIQ, “Q” for “quarantine”, to distinguish it from standard SIR type models. In our model, we do not distinguish between isolation and quarantine. We are aware of a different use of the term “quarantine” in the literature to refer to the isolation of individuals who may have been infected but are not yet symptomatic[4,5].

An interesting direction of research is to consider the isolation of individuals that have been exposed but are not yet infectious, using for example contact tracing. These strategies have the effect of increasing the isolation probability, thereby leading to possibility larger tolerable implementation delays. The influence of the underlying spatial topology of the contact network here was traded off for a clear understanding of the effects of delays in isolation. The network heterogeneity can strongly influence the dynamics[24] and it can either be tackled by means of heterogeneous mean field approximation[25,26], or by a direct study of networks with adaptive change of the structure due to network rewiring caused by the isolation protocol[27]. Such approaches should enable one to clarify the role of highly connected nodes in the isolation protocol.

Methods

Dynamical Systems Framework

Here we present a mathematical framework for the study of the delay system (1–5) [20]. Time delay implies that this system is infinite-dimensional, i.e. its state can be described by the history function, where the variables S, I, E, Q, R are prescribed on an interval of the length σ + τ + κ of maximal time delay. More specifically, let be the Banach space of -valued continuous functions with the sup-norm. Given an initial function , the solution , t ≥ 0, to the initial value problem exists and is unique[20]. Using the notationit is known that defines a C1 semi-flow on .

The conservation of mass property, namely S′(t) + E′(t) + I′(t) + Q′(t) + R′(t) = 0 for all t ≥ 0, implies that if ϕ = (ϕ, ϕ, ϕ, ϕ, ϕ) and x(t; ϕ) = (S(t), E(t), I(t), Q(t), R(t)), then S(t) + E(t) + I(t) + Q(t) + R(t) = ϕ(0) + ϕ(0) + ϕ(0) + ϕ(0) + ϕ(0) for all t ≥ 0. To obtain biologically relevant solutions we further restrict to the subset on which S(t), E(t), I(t), Q(t), R(t) ≥ 0.

In Results, we considered initial conditions ϕ representing the sudden appearance of a small group of infectious individuals. These initial conditions should be near (1, 0, 0, 0, 0), the disease-free equilibrium with constant coordinate functions. Two examples are (i) ϕ = (ϕ, ϕ, ϕ, ϕ, ϕ) where ϕ = 1 − c, ϕ = c and ϕ = ϕ = ϕ = 0, all constant functions with , and (ii) ϕ(t) = (1, 0, 0, 0, 0) for t < 0 and ϕ(0) = 1 − c, ϕ(0) = c, ϕ(0) = ϕ(0) = ϕ(0) = 0. While such a ϕ is not in , for t ≥ σ + τ + κ.

Numerically, DDEs are solved using the method of steps: Given an initial function ϕ we obtain the solution x(t, ϕ) to the DDE in the interval [0, τ + σ + κ] as the solution of the nonautonomous Ordinary Differential Equation using a 4th order Runge-Kutta with integration step 10−3 where the initial function is a known time-dependent input.

Disease-free states

For fixed α, σ, p, τ and κ, the set of initial conditions ϕ for which all component functions are constant and ϕ is identically equal to zero defines the “disease-free” equilibria. These comprise the two-parameter family of constant functions (S(t), E(t), I(t), Q(t), R(t)) = (S(0), E(0), 0, Q(0), 0) for all t, where S(0) + E(0) + Q(0) = 1. To investigate the stability of a given equilibrium solution, we linearize the system to obtain the following quasi-polynomial characteristic equation

Two zero eigenvalues correspond to directions along the manifold of equilibria. An analysis shows that all disease-free equilibria are linearly stable if S(0) ≤ 1/r(1 − ε), and unstable otherwise. At this bifurcation value, a third eigenvalue crosses the imaginary axis with nonzero speed and gives rise to a different class of endemic (nonzero infected) equilibria. If S(0) ≤ 1/r(1 − ε) no further bifurcations occur, so that these disease-free equilibria are linearly stable independently of the incubation time and isolation protocol.

The analysis above is carried out for fixed α, σ, p, τ and κ. We now turn the question around and ask for which values of these parameters is a stable equilibrium. This is how the minimum isolation probability p and critical identification times τ(p) announced in Results were deduced.

Endemic states

In addition to the disease-free states, system (1)–(5) possesses a two-parameter family of endemic equilibria consisting of constant functions (S(t), E(t), I(t), Q(t), R(t)) = (Seq, E(0), I(0), Q(0), γα−1I(0)/(1 − ε)) with I(0) > 0, Seq = 1/r andwhere I(0) and Q(0) can be considered as the free parameters. One sees immediately that Seq is not affected by the recovery rate α, the incubation period σ and isolation time κ. In particular, it is constant for the whole family of endemic equilibria. The linear stability analysis of these equilibria is very challenging from an analytical point of view. They are stable for reasonable values of κ, destabilizing as κ increases. A detailed analysis will be published elsewhere. The formal limiting case α → ∞ (immediate waning of immunity) has been studied in detail[22].

The main concern of this paper is the effect of isolation on disease outcome starting from a small outbreak. To this end, the following conserved quantities will be of major help:

The conservation of these quantities along solutions can be proved by direct computation. Restricting our attention to initial functions of the type (1, 0, 0, 0, 0) for and (1 − I0(0), 0, I0(0), 0, 0) at t = 0 for some I0(0), we have H(x(ϕ)) = H(ϕ) = 0. This implies

Equations (12) and (13) can be used to find the corresponding endemic equilibrium (Seq, Eeq, Ieq, Qeq, Req) reached from the considered initial function. More specifically, we determine Eeq = βmσSeqIeq and Qeq = βmεκSeqIeq. Inserting these values in Eq. (11), we obtain Ieq, and, hence, the unique endemic equilibrium (6)–(10).