Refined compartmental models, asymptomatic carriers and COVID-19

1.1. Context

Compartmental models, based on ordinary differential equations, have a long history of use in epidemiology. Indeed, it is now almost a century since William Ogilvy Kermack and Anderson Gray McKendrick introduced the well known SIR model, [MK], on which most compartmental modeling elaborates on. Their relevance stems for their extreme simplicity and ability to capture important qualitative behaviors, rather than their limited capacity to make quantitative predictions. As a matter of fact, the current outbreak of COVID-19, caused by SARS-CoV-2, is teaching us that most models, not only the compartmental ones, are of very limited quantitative predictive capacity. Nevertheless, there are important qualitative features of outbreaks that one may attempt to capture, for instance: under which conditions can one guarantee that there will be only one peak? This stands as a very important question which ought to be answered if one is to effectively prevent a second, or third wave of infections. Specially, having in mind the concerns to be taken once a first wave have passed. The goal of this article is to propose and analyze some compartmental models which have been constructed to capture some important features of diseases such as COVID-19. Its two major characteristics which are not captured by the standard compartmental models are:

1. The long incubation period, during which the individuals that have been exposed do not yet develop symptoms but are already contagious (at least in the second half of that period). Even though the standard compartmental models do not take into account the contagiousness of the exposed individuals, they are easy to modify in order to account for it.

2. There is a non-negligible fraction of the contagious population which either has mild symptoms or never develops them, thus passing undetected. These are the so-called asymptomatic carriers. After some period of time they are no longer contagious and, as the infected individuals, developed anti-bodies. This suggests they are, at least for now, immune. The standard compartmental models which have a carrier mode are not suitable for modeling such a behavior.

In this article, namely in Sections 3 and 4, we analyze some compartmental models that have been specially designed to handle these two features. The resulting compartmental models are more complicated than the standard ones and using a mix of precise and numerical results, we shall answer some fundamentally important questions which we now pose:

Question 1. Regarding each of the models:

1. Are there disease-free equilibrium states?

2. Having answered positively to (a), one may ask whether there are necessary or sufficient conditions for these disease free equilibrium states to be stable in some sense.

Question 2. Again, for each of the models one can ask:

1. Are there criteria guaranteeing that any outbreak of the epidemic will have a unique peak?

2. If these criteria are violated, can the outbreak have other separate peaks, i.e. second, third waves and so?

3. Can these possible later peaks be higher than the first?

Depending on which model one considers, these two questions will be answered separately in Sections 3 and 4. For instance, using the more refined model we give in Corollary 1 and Remark 7 a quantitative condition ensuring only one peak will form. The broad conclusion is that if all transmission rates are kept sufficiently small, then only one peak form, which answers item (a) of question 2. If this condition is not met, more peaks can form as we show in example 5 and these can be higher, thus answering items (b) and (c) of question 2 above.

As for the next question, it will not be answered in any precise way but we shall give some evidence on how to answer it based on numerical results. This will be done in the examples 3, 4 and 5 presented in subsection 4.2.

Question 3. What is the role of asymptomatic carriers? Can they help the spread of the disease and thus make the outbreak worse by anticipating and increasing its peak and/or can they shield the rest of the population by creating herd immunity.

One other major question which one poses now is related to the quantifying by how much specific policies reduce the rate of transmission.

Question 4. Is there a strategy to effectively reduce the rate of transmission in a manner that can be quantified?

Broadly speaking, for the general strategy or policy it is hard, if not impossible, to make such a quantitative analysis. For instance, how can we predict by how much the use of masks by the general population reduces the rate of transmission? or washing hands? Indeed, it is impossible to say it a-priori. We only know that these measures do reduce the rate of transmission but not by how much. Having this in mind, we propose one further measure whose contribution to the reduction of the transmission rate can be quantified using the models. The “idea” is to split the population into n different groups which are supposed to never interact, this reduces the contact rate of any individual by a factor of n thus reducing the overall transmission rates by a factor of n. This intuitive argument can be made more precise using the models we use. The last section of this article, Section 5, is dedicated to making such an argument and running a simulation. We shall also explain how this relates to our answer to the questions 1 and 2 above.

1.2. A note on the timing of this article

Once we are through the worst part of the first wave of the COVID-19 outbreak with relative success, meaning that only a small fraction of the population got infected and so we cannot count with group immunity, the question we must pose is: How can we proceed and go back to a relative normality so that a putative second wave is as mild as possible?

I hope that the answers to question 2 above, stated as Proposition 1 and as Corollary 1 depending on the model, may serve as good indication that it is possible, but difficult, to have only one peak. They quantify and show what one could already suspect, that the way to do so is to reduce all the transmission rates. It is enough that one transmission rate be high for the outbreak to be out of control, see Conclusion 2. How to keep the transmission rates small remains a formidable challenge which I do not attempt at answering.

As said before, there is no way to quantify how most policies reduce the transmission rates. In that direction Section 5 explains a strategy, which may be difficult to implement, but whose effect on the transmission rate can be pined down.

Disclaimer 1. I understand that the strategy of splitting the population into non-interacting groups is of difficult implementation. I do believe it cannot be literally implemented in a society which respects our civil liberties and so would remain at the level of a request to the population, which one would hope to be understood. Of course, in that scenario the different groups would interact, even if mildly. Having such considerations into account leads to modifications of the proposed solutions which will consequently modify the models one should use into mathematically more intractable ones. We may hope that the increased mathematical sophistication of those models will not lead to substantially different qualitative outcomes but instead quantitative ones. In any case, one must understand that the models are simply that: attempts to model a reality which is much more complex. They are, by no means, a truthful reflection of the real world.

Acknowledgment

I want to thank Álvaro Kruger Ramos for his comments on an earlier version of this manuscript.

Gonçalo Oliveira is supported by Fundação Serrapilheira 1812-27395, by CNPq grants 428959/2018-0 and 307475/2018-2, and by FAPERJ through the grant Jovem Cientista do Nosso Estado E-26/202.793/2019.

2.1. The modified SEIR model

We shall start by constructing a version of the SEIR model where S, E, I, R represent the fraction of the population which is susceptible, exposed, infected and recovered respectively. In order to keep the system as simple as possible one may attempt to regard the asymptomatic carriers as being part of the exposed population. For this one must regard each exposed individuals as representing an average of someone that is in the incubation period and an asymptomatic carrier. We are then led to the following system for some positive functions of time β_e, β_i, a, γ. The parameters β_e, β_i represent the transmission rates for the exposed and infected respectively, a denotes the rate of transition from the exposed to infected state and γ the rate of recovery. This is not the standard SEIR model as we are forcefully making the exposed interact with the susceptible as the first of these are assumed to be contagious. The analysis of this model is the scope of Section 3.

2.2. The SEIAR model

There are at least two easy ways to further refine this model. First, we may split the exposed into two groups, those which are not yet infectious and those which already are. Secondly, we may have some of these exposed passing directly to the recovered state to account for the asymptomatic carriers which never develop symptoms and have therefore passed undetected. Thirdly, we may distinguish the asymptomatic carriers in a more effective way as they may remain so for a longer period of time.

2.2.2. Having the asymptomatic carriers separate (SEAIR)

The goal of this model is to better capture the effects of asymptomatic carriers. In the case of COVID-19 we do know this to be a non-trivial part of all those which are contagious. However, given the small number of anti-bodies tests that have been made so far, not much quantitative knowledge regarding these is known.

In order to account for the asymptomatic carriers we consider one more compartment of the population, denoted by A, which represents the fraction of the population which is an asymptomatics carrier. Then, we propose the following model for some positive functions of time β_i, β_e, β_a, a_i, a_a, γ_i, γ_a, µ which need not be constants. This system will be carefully analyzed in Section 4. For now, we shall simply comment on the biological meaning of all the parameters on the model. Similar remarks hold for the parameters of the simpler SEIR type model. The β_i, β_e, β_a respectively represent the transmission rates of the infected, exposed (in th incubation period), and asymptomatic individuals. The parameters a_i, a_a represent the rate at which the exposed individuals respectively pass to the infected and asymptotic compartments. In particular, and respectively represent the fraction of the population which if exposed will eventually become infected or asymptomatic. The quantity represents the average time an infected individual takes to recover, while is the average time an asymptomatic individual takes to stop being contagious. Finally, µ stands for the death rate associated with the disease. For convenience, we show in figure 1 a diagramatic presentation of this model. A somewhat similar model have also recently been used in [Re] and in [Vo]. The first of these attempted at predicting the fraction of undocumented infections for COVID-19 in China. One other similar version of it using Markov Chains have also recently been used to predict the spatiotemporal spread of COVID-19 in Spain, see [Ae] for this work.

• Download figure
• Open in new tab

Figure 1.

Diagram of the SEIAR model.

2.3. Refining the models

There are a few natural ways to further refine the models proposed above.

3.1. The linearized system

From the last equation above, i.e. 3.4, a critical point will certainly have I_c = 0 which then also gives E_c = 0 from the semi-last equation 3.3. Then, the reaming equations 3.1 and 3.2 we conclude that R_c and S_c can be any two constants so that R_c + S_c ≤ 1. These are disease free equilibrium points and are obviously not isolated. Thus, the Grobman-Hartmann theorem cannot be applied but I believe the linearization to still carry important information. Thus, we shall linearize the system around such an equilibrium and use it to define the following notion.

Definition 1. Given a disease free equilibrium point (S_c, 0, 0, R_c) is called infectiously-stable if the associated linear system is such that any of its solutions (s, e, i, r) satisfies

Intuitively, we may think of this condition as saying that perturbing away from the equilibrium always makes the disease to get extinct. In order to study their stability, we shall now linearize the equations around these equilibria. The linearized system is given by

Clearly, the dynamics of (e, i) controls the whole dynamics which evolve independently. The stability is therefore dependent on the sign of the eigenvalues. As we want both of these two be negative we need to have det(A) > 0 and tr(A) < 0, where A is the matrix which controls the system. Thus, the condition for stability is which we can rewrite as

The second of this equations says that the rate at which the exposed infect other people must be smaller than the rate of recovery (plus the death rate). The first says that not only the exposed must be infecting people slower than they pass to the infected state (β_eS_c < a) but also, the already infected must infect others at a very very slow rate when compared to the recovery rate. More precisely, we have obtained the following result.

Conclusion 1. It is possible to reduce the population of those infected keeping many people S_c < S(0) still susceptible (without ever having been infected). One way to do that is to have the transmission rate of those exposed small when compared to the rate of recovery, i.e.

Furthermore, one must have the rate of transmission of those infected even smaller when compared with rate of recovery, namely that

which in particular also assumes that β_e < a, i.e. that the rate at which the exposed infect the susceptible is small when compared with the rate of transition from the exposed to infected state. We may also rewrite this last condition as

which suggests regard the left-hand-side as an analogue of the basic reproduction rate for this model. Achieving these conditions may, however, be a difficult task to accomplish. If they are not, then the final population of susceptible may be very low which is quite bad.

3.2. Examples and numeric simulations

Example 1. Consider this model with S(0) = 0.99, I(0) = 0, E(0) = 0.01 and R(0) = 0 (for small nonzero I(0), E(0) the outcome seems to not be very dependent on the precise value) whose simulation is shown in figure 2. In this precise case, only a very small fraction of the population got away without being infected. Indeed, form equation 3.5 for the stability of the disease free equilibrium we find that it must satisfy

which is in agreement with what we see in figure 2.

• Download figure
• Open in new tab

Figure 2.

Example when β_i = 0.9, β_e = 2.5, a = 1, γ = 0.9, µ = 0.1 which has and initial conditions S(0) = 0.99, I(0) = 0, E(0) = 0.01 and R(0) = 0. The red, blue, green and purple curves respectively denote the fraction of the population composed of susceptible, exposed, infected and recovered.

On the other hand, we may consider the same model but supposing that social isolation and other measures have been put in place to reduce the transmission rate. In figure 3 we run one such example. In this case, equation 3.5 gives that the stable disease free equilibrium to which the model converges satisfies which is compatible with what can be inferred from figure 3.

• Download figure
• Open in new tab

Figure 3.

Example when β_i = 0.45, β_e = 1.25, a = 1, γ = 0.9, µ = 0.01 which has and initial conditions S(0) = 0.99, I(0) = 0, E(0) = 0.01 and R(0) = 0. The color code is that of figure 2.

Remark 1. It would be interesting to actually get a lower bound on S_c rather than an upper bound.

3.2.1. Example with two peaks

We could now try and make β_i and β_e periodic in a way that, on average, the inequalities in Conclusion 1 hold true. In that case, I expect that, after the peak, the number of infected people will also decrease, in average, with time. Of course, it may be that, with no extra measures, the actual values of β_e and β_i are so high that by opening up (even if only for a short period of time) will make the task of making the inequalities above hold in average very difficult. It is exactly this scenario that we explore in the next example where after a strict control of the transmission rate, it is once again allowed to become large. This is an alert to the dangers of a non-careful reopening, after social distancing measures had been put in place, for an insufficient time so that herd immunity develops. In this example we can actually see that the second peak can be made larger than the first which provides an answer to the items (b) and (c) in question 2.

Example 2. Figure 4 shows an example where after seeing positive signs in the decrease in the number of infected, the rate of transmission is allowed to increase. This leads to a formation of a second peak of infections which ends up infecting almost everyone.¹ Indeed, we can see in the plot in figure 5 that the conditions in Conclusion 1 are violated right before time t = 5. This leads to a change of tendency which will makes the curve of infected restart increasing and eventually leading to the formation of the second peak. This is an example of something we want to avoid.

• Download figure
• Open in new tab

Figure 4.

Example of an un-careful reopening which leads to a second peak.

• Download figure
• Open in new tab

Figure 5.

Plots of and which one can see to surpass the value 1 before t = 5.

Remark 2. This example leads to a fundamental question. Is there a criteria which guarantees that no second peak will form? This is precisely the problem raised in th problem (a) by question 2. Preferably a criteria depending on β_e and β_i which are the quantities that can be somewhat controlled using social distancing and other measures. We shall answer this question next in Proposition 1.

3.3. Rigorous qualitative results

First, notice that from the first equation 3.1 we find that so S(t) is decreasing. Consider the asymptotic fraction of susceptible population and recall that it is to our interest that S_c > 0 and as big as possible (of course S_c ≤ S(0) ≤ 1). Similarly, we shall define I_c and E_c as the asymptotic values of I and E if they exist.

Lemma 2. Suppose that S_c > 0, then I_c = 0 = E_c.

Proof. Again, from the first equation 3.1 we have which will converge to zero if I and E do not themselves uniformly converge to zero.

Political measures can more easily affect the transmission rate than the recovery and death rate which are more dependent on medical conditions. Thus, in the next result we find a, potentially useful, criteria for having only one peak which does not assume the transmission rates β_i, β_e to be constant. This provides a rigorous answer to the item (c) of question 2.

Proposition 1. Suppose that a, γ, µ are all constant and (S, E, I, R) is a nonconstant solution of 3.1–3.4. If there is a time t_∗ such that İ(t_∗) < 0 and

then, for t ≥ t_∗, the fraction of infected population I(t) has at most one critical point. Moreover, if any such critical point exists, then it is a maximum. In particular, if

then I(t) at most a unique maximum.

Proof. First notice that if I and E both vanish at some point, then the solution is constant. Then, we may suppose this is not the case and differentiate the third equation 3.3, using the second 3.2 to substitute for Ė. This gives which at a critical point of I yields

As between any two maxima there must be a minimum, in order for I to only have one maximum it is enough that it has no minimum. At a minimum we would have Ï ≥ 0 which from the previous computation is ruled out if

However, the condition that we are at a critical point of I(t), i.e. İ(t) = 0, yields that aE = (γ + µ)I which upon inserting above gives aβ_iS + (β_eS − a)(γ + µ) < 0, which in terms of S reads as

Thus, if S always satisfies this inequality past some time, then there is at most one critical point of I and this must a maximum, if it exists at all.

Remark 3. The previous result gives a criteria to guarantee that a second peak will not form. For that, one must try and control the β_e and β_i so that equation 3.7 holds.

4.1. Linear analysis

Let us now find the critical points of this system. From equations 4.3, 4.4 and 4.5 we find that for the generic values of a_a, a_i, γ_i, γ_a, µ all E, A, I must vanish at a critical point. From the remaining equations we find that S and R can be any arbitrary constants as long as S + R ≤ 1. We shall write such a critical point as c = (S, E, A, I, R) = (S_c, 0, 0, 0, R_c). Notice in particular that any such critical point is disease free.

Then, the linearised system at such a critical point is given by

In this case, we have that the equations for where x:= (e, a, i) do control the whole system. This subsystem is then given by where

While it is possible to compute the eigenvalues of this system in general, the formulas are extremely unwieldy and so it is hard to find a general statement which decides on the stability of the system in general. Nevertheless, as necessary condition for stability is that det(A) < 0 and tr(A) < 0 which respectively give

If β_e < a_i + a_a + γ_a + γ_i + µ, then the first equation above suggests that the quantity may behave as the basic reproduction rate for this model. Hence, and in order to have large (closed to 1) values of the asymptotic fraction of susceptible, we would like to keep 4.7 below 1. This is then compatible with having S_c close to 1, at least from the point of view of the first equation in 4.6.

Conclusion 2. The above discussion suggests that for an outbreak to be under control we must have β_e < a_i + a_a + γ_a + γ_i + µ and the quantity in equation 4.7 must be kept below 1. Furthermore, from formula 4.7 one can immediately see that it is enough that one of the transmission rates β_e, β_i, or β_a to be out of control, i.e. to be large, for the outbreak to be out of control!

4.2. Controlling an outbreak

The main problem when faced with an outbreak is to keep it under control. But what does that mean in the context of the model we are considering. One possible definition is that the asymptotic value os susceptible is large. Alternatively, we may attempt at keeping the total number of infected people, i.e. as small as possible. The proof of Lemma 2 show that these two definitions are equivalent for the model 3.1–3.4. However, in model at hand there is the possibility that as susceptible passes to the recovered state taking the asymptomatic route which makes this model fundamentally different. Indeed, from integrating the first equation 4.1 we find that

This shows that maximizing S_c is equivalent to minimizing

From the point of view of the second definition, one regards the asymptomatic to be innocuous, which seems reasonable from a purely medical point of view as they never become sick. Furthermore, it is conceivable that they can create herd immunity which shields the rest of the population. However, intuitively speaking, having a relatively large number of asymptomatic may also be bad, at least initially while S is large, as they may move undetected and “infect” even more people. The numerical conclusion to which we arrive is that depending on the exact parameters of the model the asymptomatic can carry both roles.

4.2.1. Examples

Below we shall run some simulations which seem to suggest the previous intuitive reasoning to be correct in different settings. In order to this we shall fix the parameters β_i, β_e, β_a, γ_i, γ_a, µ and vary a_i and a_a. The examples in this section have been constructed in order to provide answers to the question 3 raised in the Introduction.

Remark 4. A few words must be said about the choice of the parameters β_i, β_e, β_a, a_i, a_a, γ_i, γ_a, µ ruling the model. It is at the date medically impossible to determine the values of all these and so it is meaningless to assign them specific values and pretend we are modeling the true outbreak. Nevertheless, I must justify the values I will be assigning in my model. I shall suppose that

and recall that these encode the number of infections that the infected, exposed and asymptomatic cause by unit time. This choice may seem strange but there is a reason for this choice. I assume that the infected, while possibly being extremely infectious, are dealt with carefully and so do not infect as many other people by unit time as the exposed and asymptomatic. This justifies β_i < β_a and also β_i < β_e. The exposed, even though not having been detected yet, are not as infectious as they will eventually become later thus β_e < β_a. Notice that we do not have β_e < β_i as the infected individuals will be detected and dealt with in a way that avoids further infections. Without this extra care that inequality would hold. From all this we then have: β_i < β_e < β_a, as claimed.

Next, we shall assume in our simulations that

which means that the average time an asymptomatic individual takes to reduce its viral charge to zero is smaller than the time an infected indivual needs . This is highly debatable and, even though seemingly reasonable, I have no way to better justify this choice.

The simulation in the next example 3 illustrates a case where having a larger fraction of asymptomatics carriers did not make the disease develop faster and actually lead to a smaller peaks of infection. This need not always be the case as we shall later see in example 4.

Example 3. In all the simulations we shall run in this example we shall have β_i = 0.9, β_e = 1.5 and β_a = 2.8 while γ_i = 0.8, γ_a = 1.2 and µ = 0.01.

In Figure 6 we shall plot the number of infected people, these are those which actually get sick, as a function of time where we compare different simulations obtained with the same value of a_i + a_a = 1 but different values of a_i and a_a individually. Indeed, from this simulation we can infer that the smaller a_a is, or equivalently the larger a_i is, the further high and late the peak is. Intuitively we can justify this from the fact that a large a_a gives a lot of asymptomatic carriers, which when recovered create a group immunity shielding the rest of the population.

• Download figure
• Open in new tab

Figure 6.

In black, green, blue, red and purple the values of I(t) are plotted for the simulations obtained from a_a being 0.1, 0.3, 0.5, 0.7 and 0.9 respectively.

I believe the values of the transmission rates β_e, β_i, β_a in the previous example 3 are reasonable and somewhat realistic. Nevertheless, in order to better understand the role that asymptomatic carriers can carry it is actually more convenient to slightly exaggerate the values by having β_a ≫ β_e. We shall do so in the next example which illustrates one other role that asymptomatic carries can have. Namely, that in an early stage, i.e. when S is still large, having more asymptomatic carriers can cause a faster increase of the infected population anticipating and increasing the peak.

Example 4. In this example we use the same constants of the previous example except for the values of β_i and β_a which we shall set to be β_i = 0.5 and β_a = 7. As explained before, the fact that β_a is so big when compared to β_i means not only that the disease is extremely contagious but also that those individuals which are showing symptoms are isolated and handled extremely carefully while the asymptotic ones are not. In figure 7 we have plotted in black and green the simulation obtained when a_a = 0.1 (and a_i = 0.9) and a_a = 0.2 (and a_i = 0.8) respectively. In this very particular example, we see that the higher probability of an exposed person to become asymptomatic leads to higher degree of contagion which itself leads to a higher earlier peak of the infected population. It is natural to inquire about the actual population of asymptomatic carriers in both of these simulations, these are plotted in figure 8.

• Download figure
• Open in new tab

Figure 7.

In black and green are the plots of I(t) for the simulations obtained from a_a being 0.1 and 0.2 respectively.

• Download figure
• Open in new tab

Figure 8.

In black and green are the plots of A(t) for a_a being 0.1, 0.2 respectively.

4.2.2. Second waves

In the simpler model given by the system 3.1–3.4 we saw in Example 2 that a second wave can form. In the same manner as in that case, we expect this will happen when the quantity derived in subsection 4.1, rises above 1. We shall now see an example of a solution to the model 4.1–4.5 for which a second wave forms. In fact, we shall see that the second peak can be larger than the first and thus this example will serve as an answer to the items (b) and (c) of question 2 using this more refined model. Also, we shall use this opportunity to further elaborate on the role of asymptomatic carriers giving further input towards question 3 in the Introduction.

Example 5. Let β_i(t) = 0.9f (t), β_e(t) = 1.5f (t), β_a(t) = 2.8f (t), a_a = 0.2 and a_i = 0.8, γ_i = 0.8, γ_a = 1.2 and µ = 0.01, for some function f (t) which was carefully designed² but whiose specific form is unimportant for now. As a matter of showing the importance of the quantity 4.7 we shall plot this in figure 9 which we can see has two large plateaus, well above 1, connected by a region where it is substantially smaller. The simulation run for these values is presented in figure 10 and is made with the initial conditions S(0) = 0.99, E(0) = 0.01, I(0) = 0 = R(0).

• Download figure
• Open in new tab

Figure 9.

Plot of .

• Download figure
• Open in new tab

Figure 10.

Plot of I(t) and A(t) in green and black. These represent the infected and the asymptomatic carriers respectively.

One sees that the first peak is too small for a sufficient fraction of the population to acquire immunity and the high increase of the rate of transmissions, here enconded in the quantity 4.7 leads to a second peak of the outbreak. Indeed, this second wave of the outbreak creates an even larger peak than the first. Of course, the actual values of the transmission rates we used in this example may not be well adapted for a realistic modeling.³ Nevertheless, this serves as an example to stress the possibility of having a second outbreak that may be worse than the first.

As a way of comparison, and to further investigate on the role of asymptomatic carriers, we present in figure 11 a simulation of the same model except that we now have a much higher fraction of asymptomatic carriers, in this case a_a = 0.8 and a_i = 0.2. This means that each individual has a 80 per cent probability of becoming asymptomatic when exposed to the disease. In this case, we see that the large quantity of asymptomatic carriers keeps the peaks of infection much smaller. In fact, we also find that in this case the second peak is smaller. This can be interpreted in terms of the herd immunity acquired by the asymptomatic carriers which shields the rest of population while keeping the number of infections smaller.

• Download figure
• Open in new tab

Figure 11.

Plot of I(t) and A(t) in green and black respectively.

Remark 5. Example 5 raises the question of whether it is possibly to find conditions under which an outbreak will develop a unique peak. This is precisely the constant of item (a) in question 2. As we shall see in Corollary 1 of the next section such a criteria can be found and cast in terms of the quantity 4.7.

4.3. Precise qualitative results

As in Lemma 2 we find that in order for to be positive we must have all A(t), I(t) and E(t) converging to zero as t + → ∞ The other natural easy question to investigate is that of finding criteria for which there will be only one peak of the outbreak, very much as done in Proposition 1 for the simpler model 3.1–3.4. Here, we shall find an analogue of that which holds under more restrictive hypothesis.

Proposition 2. Let (S, E, A, I, R) be a solution to 4.1–4.5 with a_i, a_a, γ_i, µ constant and S_c > 0. Then, there is a constant c > 0 so that I for all t ≥ 0. Furthermore, if

the function I(t) has at most one critical point and this is a maximum if it exists.

Proof. First notice that as S_c > 0 we have that both A and I converge to zero and so there is a constant c as in the statement. Thus, in order to prove the result we shall follow the strategy of Proposition 1. This consists in finding conditions which guarantee that Ï < 0 at any critical point of I. If that is the case, then I can have no minimum and thus will have at most one maximum, as in-between any two maxima there must be a minima. Thus, we start by computing which at a critical point of I further satisfies İ = 0, i.e. a_iE = (γ_i + µ)I, and so we can rewrite

Thus, we want to find conditions under which the following inequality holds

If we further have , then for inequality 4.8 to hold, it is enough that which we can rewrite as

As S ≤ 1 we find that this condition is immediate from that in the statement.

Remark 6. From the proof of Proposition 2 we find that, in fact, for I(t) to only have at most one maximum it is enough to show that

or even

Of course, one would like to keep S(t) as close to one as possible and so the first of this seems to be a more reasonable test. Furthermore, it does not require to know the value of S(t) but that of the ratio A(t)/I(t) which may be obtained by sampling with anti-body tests for example.

Of course, the previous proposition would be much more useful if for any given solution (S, E, A, I, R) one can determine the value of c.

Lemma 3. Let (S, E, A, I, R) be a solution to 4.1–4.5. Then, for all t ≥ 0 we have

Furthermore, let A(0), I(0) respectively denote the initial value of asymptomatic and infected individuals. If a_iA(0) ≤ a_aI(0) and γ_a ≥ γ_i + µ, then

Proof. Equations 4.3, 4.4 can be equivalently rewritten as from which we infer that and using the hypothesis that γ_a ≥ γ_i + µ

Integrating this yields which upon rearranging gives the result in the statement.

Putting this Lemma 3 together with Proposition 2 we find the following Corollary. This serves as a rigorous answer to the inquire raised in item (a) of question 2.

Corollary 1. Let (S, E, A, I, R) be a solution to 4.1–4.5 with a_i, a_a, γ_i, µ constant and S_c > 0. Suppose that a_iA(0) ≤ a_aI(0), γ_a ≥ γ_i + µ and

Then, the function I(t) has at most one critical point and this is a maximum if it exists.

Remark 7. In the previous Corollary 1, the hypothesis that γ_a > γ_i + µ holds immediately if the average time for an asymptomatic carrier to stop being contageous is smaller than that required by the average infected individual. Still under this hypothesis, if the condition in equation 4.9 holds, then we also have

which should be compared with Conclusion 2.